This is a test version of Biostars. For the public version, visit https://www.biostars.org.
comparing mapping results in two tools

Friends,

I aligned my reads in both bowtie2 and tophat. Now how I can see and compare the mapping result in two mentioned tools?

Thank you

rna-seq tophat bowtie2

hello,

what version of bowtie you used in tophat? and what version of tophat you have used?

bowtie2 and tophat2

2 answers

Why do you want to do that ?. It is non sense..

TopHat will do a initial mapping using bowtie2 (so you get the same that using bowtie2 alone), and it will do a segmentation of the unmapped sequences to try a new mapping using bowtie2 again to detect spliced sequences

thank you Antonio,

sorry, how i can be aware of rate of the reads which aligned in multiple place on genome in tophat output?

Tophat2 will give them a MAPQ or 0, 1, 2, or 3. Non-multimappers are given a MAPQ of 50 (last time I looked at least).

thank youuuuuu Devon

Take a look to this page

Inside the SAM file, you can read a FLAG value which is interpreted by the page shown above. It will be useful for you

If you have a sam/bam file, you can compare which mapping tool could align more number of reads. properly and improperly paired reads and overall number of reads aligned would help you.

Log in to answer this question.