Question: Bowtie & Tophat
gravatar for k.nirmalraman
5.8 years ago by
k.nirmalraman980 wrote:

My understanding to alignment tools are quite limited. Any suggestions are quite welcome. There are quite a lot of articles and questions already on this, but I am not able to find an answer in this context.

I am trying to understand what difference does it make to align the reads by using bowtie and tophat and cufflinks.

My data: Illumina short reads, single end reads . I am not interested in isoforms.

Currently I am using bowtie and aligning with 2/3 mismatches. (I am trimming away the barcodes, adapter sequences, also I am trimming the far end of sequences... At the end I have reads of length either 35 to 40). On an average I have about 78 - 80% of reads aligning to genome.

Do you see any limitations in directly using Bowtie over the other tools.

tophat cufflinks bowtie • 9.8k views
ADD COMMENTlink modified 22 months ago by mrmkanni930 • written 5.8 years ago by k.nirmalraman980

./configure --prefix/opl/tophat_2.1.1 --with-boost=/opt/boost_1_62_0 running this command to install tophat2bu its showing an error bash: ./configure: No such file or directory can any one plzzz give me a solution.

ADD REPLYlink written 22 months ago by mrmkanni930

This is not an answer to this discussion and should be opened as a separate question. I moved this to a comment now.

You are in the wrong directory for that command. First cd into your downloaded directory.

Finally: tophat is superceded by hisat and therefore shouldn't be used anymore. Also STAR is a better alternative.

ADD REPLYlink modified 22 months ago • written 22 months ago by WouterDeCoster36k
gravatar for Federico Giorgi
5.8 years ago by
Columbia University
Federico Giorgi440 wrote:

Use Tophat if your reference is the genome. Use Tophat if your reference is the transcriptome (e.g. a list of cDNA sequences, one per gen or one per isoform).

ADD COMMENTlink written 5.8 years ago by Federico Giorgi440

So use tophat in both cases? Is that what you meant to write?

ADD REPLYlink written 5.8 years ago by Neilfws48k
gravatar for Ashutosh Pandey
5.8 years ago by
Ashutosh Pandey11k wrote:

Bowtie is used for genomic reads alignment against the reference genome. Tophat is used for transcriptomic reads alignment against the reference genome and it uses bowtie in the first phase of alignment. But in addition to bowtie, tophat can also align reads that span exon-exon junctions. so you will get more read aligned against the reference genome if you use tophat for RNAseq data. Cufflinks is not an alignment tool. It is primarily used to calculate transcript abundance and abundance of different isoforms of the same gene.

Check the pages for bowtie, tophat and cufflink for more information.

ADD COMMENTlink written 5.8 years ago by Ashutosh Pandey11k
gravatar for Geparada
5.8 years ago by
Geparada1.4k wrote:

If you use directly bowtie to align RNA-seq data to reference genome you won't map any read to splice junctions of the mRNAs. Tophat cat do the gaped alignments required to map reads to the splice junctions, but it's much more slower than bowtie.

As your reads are shorts, you probably get similar results with both programs due to the low probability that a short read has to map to a splice junction. So it's a good "quick n' dirty" approximation to start.

ADD COMMENTlink written 5.8 years ago by Geparada1.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1120 users visited in the last hour