SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



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

Reply
 
Thread Tools
Old 01-09-2015, 06:30 AM   #1
bbm
Member
 
Location: North Carolina

Join Date: Sep 2011
Posts: 38
Default 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
bbm is offline   Reply With Quote
Old 01-09-2015, 06:42 AM   #2
bbm
Member
 
Location: North Carolina

Join Date: Sep 2011
Posts: 38
Default

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?
bbm is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:44 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO