Different Results Of Tophat2 And Bowtie2

bbm · 01-09-2015, 07:30 AM

Hello all,

I ran the same fastq file with tophat2 and bowtie2, and got very different mapping rates. I wonder if any of you will help me on this puzzle. It is single-end, 60-100 bp after trimming. Thanks a lot!!
------------------------------

My code for tophat2 alignment:

tophat -G /gpfs_share/hlibyar/hlibyar/tophat/Am_genome.gff3 -o tophat_output_WO_48h --no-novel-juncs Am_genome /gpfs_share/hlibyar/hlibyar/trim/WO_48h_ACAGTG_L003_R1_001.trimmed60.fastq

The summary showed:
Reads:
Input : 19571887
Mapped : 1013220 ( 5.2% of input)
of these: 14267 ( 1.4%) have multiple alignments (14 have >20)
5.2% overall read mapping rate.

The code for bowtie2 alignment:
bowtie2 -x Am_genome -U /gpfs_share/hlibyar/hlibyar/trim/WO_48h_ACAGTG_L003_R1_001.trimmed60.fastq -S outputWO48trim60.sam

and the result is showing below:

19571887 reads; of these:

19571887 (100.00%) were unpaired; of these:

15594391 (79.68%) aligned 0 times

3948682 (20.18%) aligned exactly 1 time

28814 (0.15%) aligned >1 times

20.32% overall alignment rate

bbm · 01-09-2015, 07:42 AM

I also understand Tophat is splice-aware. My data is transcriptome (mRNA), I expect Tophat to map at a similar rate to Bowtie. One difference is tophat maps reads that span splice junctions, right?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Getting much higher coverage with bowtie2 than tophat2	Studentlost	RNA Sequencing	34	01-09-2015 07:39 AM
tophat2 : Bowtie2 not found...	bsmith030465	Bioinformatics	5	11-10-2014 11:28 AM
Tophat2 can't find Bowtie2	rnavon	Bioinformatics	5	02-27-2013 07:51 AM
Tophat2 very slow when running over Bowtie2	jdenvir	Bioinformatics	2	02-18-2013 06:28 AM
bowtie2 vs. Tophat2 RNA-Seq	tschauer	Bioinformatics	2	12-18-2012 03:45 AM

01-09-2015, 07:30 AM	#1
bbm Member Location: North Carolina Join Date: Sep 2011 Posts: 38	Different Results Of Tophat2 And Bowtie2 Hello all, I ran the same fastq file with tophat2 and bowtie2, and got very different mapping rates. I wonder if any of you will help me on this puzzle. It is single-end, 60-100 bp after trimming. Thanks a lot!! ------------------------------ My code for tophat2 alignment: tophat -G /gpfs_share/hlibyar/hlibyar/tophat/Am_genome.gff3 -o tophat_output_WO_48h --no-novel-juncs Am_genome /gpfs_share/hlibyar/hlibyar/trim/WO_48h_ACAGTG_L003_R1_001.trimmed60.fastq The summary showed: Reads: Input : 19571887 Mapped : 1013220 ( 5.2% of input) of these: 14267 ( 1.4%) have multiple alignments (14 have >20) 5.2% overall read mapping rate. The code for bowtie2 alignment: bowtie2 -x Am_genome -U /gpfs_share/hlibyar/hlibyar/trim/WO_48h_ACAGTG_L003_R1_001.trimmed60.fastq -S outputWO48trim60.sam and the result is showing below: 19571887 reads; of these: 19571887 (100.00%) were unpaired; of these: 15594391 (79.68%) aligned 0 times 3948682 (20.18%) aligned exactly 1 time 28814 (0.15%) aligned >1 times 20.32% overall alignment rate

01-09-2015, 07:42 AM	#2
bbm Member Location: North Carolina Join Date: Sep 2011 Posts: 38	I also understand Tophat is splice-aware. My data is transcriptome (mRNA), I expect Tophat to map at a similar rate to Bowtie. One difference is tophat maps reads that span splice junctions, right?