qual_reads: Wrapper function to filter sequencing reads based on distance...

In arnederoeck/NanoSatelliteR: QC, plotting, and clustering of tandem repeat segmentation created by NanoSatellite.

Description

Returns a data.frame containing sequencing reads with a "flank" and/or "center" normalized dynamic time warping distance smaller or equal to a user specified cutoff value.

Usage

qual_reads(df, center_cutoff = 1, flank_cutoff = 1)

Arguments

df	A data.frame, generated by load_summary()
center_cutoff	A numeric value [0-1] corresponding to the maximum center normalized dynamic time warping distance
flank_cutoff	A numeric value [0-1] corresponding to the maximum flank normalized dynamic time warping distance

Examples

df <- load_summary("/storage/NanoSatellite_chunks/")
qc <- summary_qc(df)
df2 <- qual_reads(df,qc$center_cutoff)

arnederoeck/NanoSatelliteR documentation built on May 6, 2019, 1:35 a.m.