Bowtie2 alignment: Reads that aligned concordantly >1 times
1
0
Entering edit mode
6.8 years ago
elfeusilva • 0

Hello Guys!

I would like to extract reads that "concordantly >1 times" from SAM output file and make a new file with only this reads [reads that aligned concordantly > 1 times]. I did a alignment on Bowtie2 software.

Thanks very much.

35133445 reads; of these:
  35133445 (100.00%) were paired; of these:
    33688527 (95.89%) aligned concordantly 0 times
    1373944 (3.91%) aligned concordantly exactly 1 time
    70974 (0.20%) aligned concordantly >1 times [make a new file with this reads.]

 

bowtie bowtie2 alignment score map • 4.4k views
ADD COMMENT
3
Entering edit mode
6.8 years ago

Something like samtools view -h foo.bam | grep -v "NH:i:1" | samtools view -Sbo foo.multimappers.bam - will probably work. You could also just pipe to a SAM file rather than samtools if you want.

ADD COMMENT

Login before adding your answer.

Traffic: 2344 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6