Removing duplicates
It is always a good practice to check that your sequence reads don't contain too many duplicates.
# Sample with duplicates of class: ShortReadQ
dfqsample
# Get the reads from dfqsample
mydReads <- sread(dfqsample)
# Counting duplicates
table(srduplicated(mydReads))
How would you go about removing duplicated reads in a file? Pay attention to what the condition should be in this filter.
This exercise is part of the course
Introduction to Bioconductor in R
Hands-on interactive exercise
Turn theory into action with one of our interactive exercises
