More filtering!
Awesome! Now that you've had some practice with filtering reads, let's use the function polynFilter()
. This function selects reads that contain less than a given number of duplicate nucleotides. For example, polynFilter(threshold = 20, nuc = c("A"))
will select all reads that contain less than 20 A's. The parameter nuc
is a character vector containing IUPAC symbols for nucleotides or the value "other"
for all non-nucleotide symbols.
The fqsample
object is available in your workspace.
This exercise is part of the course
Introduction to Bioconductor in R
Hands-on interactive exercise
Have a go at this exercise by completing this sample code.
# Check reads of fqsample
___
# Create myFil using polynFilter
myFil <- ___
# Check myFil
myFil