How to bin how many alignments a read maps
0
0
Entering edit mode
23 months 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

binning • 414 views
ADD COMMENT
0
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

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