Tophat multiple alignment and mapping rates
0
3
Entering edit mode
7.1 years ago
anlin00007 ▴ 90

Hello,

I am using tophat to map the ribosome profiling data. First I limited the multiple hits report up to 2, which is:

-g/--max-multihits <int> = 2

Then, I changed this parameter to 5, with no more other changes. However, from the align_summary file, I can see the mapping rates increased, which means more reads mapped to the reference. Based on the tophat manual, this parameter should just change the number of multiple hits in the report which has no effect on the alignment?

I am confusing why the mapping rate changed.

RNA-Seq next-gen alignment • 3.4k views
ADD COMMENT
1
Entering edit mode

Tophat is one of those programs where a lot is going on under the hood that isn't documented. I wonder if this parameter is getting passed to bowtie2 at some point and that's limiting it to return reads that have at most this number of alignments. If you grep through the BAM file and look at the NH:i:3, NH:i:4 and NH:i:5 alignments and add them up (only counting once per read), is that perhaps the difference in number of aligned reads?

ADD REPLY
0
Entering edit mode

The 'NH:i:3' flags refer to the multiple alignments, not the reads. I am wondering why there are more reads get aligned with no change in the parameters relate to alignment strength.

ADD REPLY

Login before adding your answer.

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