Different Results Of Tophat And Bowtie2 · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Different Results Of Tophat And Bowtie2

I am trying to map RNASeq paired reads (60bp) (breast cancer) to human genome (without mt genome) using Tophat (v2.0.10) and I am returning an alignment rate of 64%. tophat -p 8 -r 200 --mate-std-dev 50 --b2-very-sensitive --phred64-quals -g 10 --microexon-search -G Homo_sapiens.GRCh37.74.gtf GRCh37.74 ff1.fastq ff2.fastq

I used the same reads in Bowtie version 2.1.0 and returned an alignment rate of 79.58% bowtie2 -p 8 --phred64 --very-sensitive -x GRCh37.74 -1 ff1.fastq -2 ff2.fastq -S out.sam

I checked my reads with fastqc and they are ok. but my problem is that why results of tophat and bowtie2 are different? I am really confused. if possible help me.

Thanks in advance

tophat tophat2 bowtie2 mapping aligner alignment

• 3,423 views • link

updated 12.5 years ago by Neilfws

• written 12.5 years ago by Mohammad Reza Bakhtiarizadeh

Shouldn't bowtie2 be ran with TopHat together?!

• 0 views • link

written 11.8 years ago by Parham

1 answer

It looks like you are providing a lot more parameters to TopHat. What happens if you run TopHat (or Bowtie) with just the default parameters?

• 137 views • link

written 12.5 years ago by Charles Warden

Thanks for your answer, but I try it and get a sane results. I really don't know what is the problem :(

• 0 views • link

written 12.5 years ago by Mohammad Reza Bakhtiarizadeh

Log in to answer this question.

More posts like this

Different version of Tophat gives me different number of mapped reads
written 8.0 years ago by BioDH •

I have paired-end sequenced RNA-seq files(Illumina, fastq). I trimmed the reads by trimmomatic java -jar <path>/trimmomatic-0.36.jar PE -threads 4 -phred33 1.fastq 2.fastq 1_trim1.fastq 1_unpaire d1.fastq …
Reducing error bars on TopHat RNA expression analysis
written 8.7 years ago by waterspring •

Using Tophat/Bowtie, I'm producing an RNA expression (FPKM) plot of a paralog in three different samples. The problem is the first specimen (blue bar in …
Low mapping rate for Tophat2
written 8.8 years ago by fadhil.abubaker •

Hi, I'm somewhat new to bioinformatics, so please bear with me. I'm running tophat2 on some fastq files using the HG38 as reference. This is …
sam2bed processing paired-end data and producing unexpectedly short intervals in bed file....
written 8.9 years ago by chrisclarkson100

I have downloaded the 2 paired end fastq files from ENA- documenting MNase Seq data for nucleosome occupancy across the mouse genome: http://www.ebi.ac.uk/ena/data/view/SRR2034494 I mapped …
RSEM ignores quality scores?
written 9.0 years ago by AW

Dear all, I want to quantify gene expression for my Trinity transcriptome assembly using RSEM. I am getting confusing results depending on how I specify …
Bowtie2/HISAT using single-end to generate FPKM???
written 9.2 years ago by jamieson.pierce •

Hi all, I am familiar with the fact that paired-end reads are used to generate FPKM, and single-end reads are used to make RPKM. I …
Very low concordant rate using TopHat2
written 10.8 years ago by AT •

I have paired end 76 bp reads for RNA seq data from FFPE samples and I am using tophat2 (version 2.1.0) using mostly default options: …
why tophat2 saying that --max--deletion-length not recognized???
written 10.9 years ago by zizigolu

Hi friends, I typed the below syntax: ``` [izadi@lbox161 bowtie2-2.2.5]$ $TOP/tophat2 --no-coverage-search --b2-sensitive --read-gap-length 12 --read-mismatches 2 --read-edit-dist 13 --max-insertion-length 20 --max--deletion-length 30 -p 8 …
Bowtie2 settings in Tophat
written 12.1 years ago by zioutisch •

<p>I tried to run Tophat with different bowtie2 options. I actually used the preset bowtie2 options (i.e. --b2-sensitive). Looking for differences between different runs, i …
Tophat Prep_Reads Error
written 14.1 years ago by Mohammad Reza Bakhtiarizadeh

<p>Hi I am running tophat on a test reads (<a href='http://tophat.cbcb.umd.edu/downloads/test_data.tar.gz'>Test data in tophat site</a> and the code was tophat -r 20 test_ref reads_1.fq reads_2.fq) …