Number of secondary alignments per read
0
0
Entering edit mode
4.5 years ago
am3 ▴ 80

Hello! If I have a SAM file produced by aligning using BWA-MEM, is there a way to tell how many secondary alignments each read produced, or how many reads produced how many secondary alignments? Let's say I do samtools flagstat and see that I have 100 secondary alignments. I want to know whether that's because one of my reads mapped to 100 different places, or because 100 reads each mapped to one secondary location, or anything in between. How can I extract that information from the SAM file? Thank you!

bwa bwa-mem samtools alignment • 1.1k views
ADD COMMENT
0
Entering edit mode

if you combine the read ID with the match ID (+position) you should be able to deduct this. Your first situation would look it one to many, the other many to many .

ADD REPLY

Login before adding your answer.

Traffic: 2613 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