Entering edit mode
2.6 years ago
rob234king
▴
610
Is there a tool that will give you from a bam, binned distribution of how many alignments reads in a bam have, as below for a bam file?
A. Reads with one alignment 100
B. Reads with two alignments 50
C. Reads with three alignments 10
etc
What aligner did you use? Almost no aligner will report all potential alignments - even the multimapping aligners have limits to prevent the data size explosion of reporting tens to hundreds of thousands alignment locations per read that maps to highly repetitive sequence.
Also, what underlying question are you trying to answer? There will be a better way of answering whatever (non-homework) question you have than a read multimapping histogram.