The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Last chance to register for the GATK workshop next week in Basel, Switzerland! http://www.sib.swiss/training/upcoming-training-events/training/gatk-workshop-lecture

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: 543Member, Dev ✭✭✭✭

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