extract duplicated reads

ShouyongPengShouyongPeng BostonMember Posts: 1

Hi

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

Shouyong

Tagged:

Answers

  • pdexheimerpdexheimer Member, Dev Posts: 534 ✭✭✭✭

    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 Administrator, Dev Posts: 10,287 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.