extract duplicated reads

ShouyongPengShouyongPeng BostonPosts: 1Member

Hi

It would be great if we can extract duplicated reads and check where they fall on the genome.

Shouyong

Tagged:

Answers

  • pdexheimerpdexheimer Posts: 419Member, GSA Collaborator ✭✭✭✭

    Since they're flagged in the BAM file, this isn't too hard. I use samtools for tasks like this: samtools view -f 0x400 -b -o Sample_dupes.bam Sample.bam

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 7,192Administrator, GATK Developer admin

    Ah, that's a good idea. Much simpler than writing a new read filter :)

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.