Entering edit mode
Learner • 0@a923bf09
Last seen 5 weeks ago
I'd like to use Rsamtools to remove duplicate reads from a bam file. I'm looking for a solution similar to samtools markdup -r -s in.bam out.bam. Can anyone tell me how to do this in R, preferably with Rsamtools.
For completeness, an initial answer on a GitHub issue on the same topic.