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

Hi,

One of my sequnce data is aligned with tophat, alignment rate is 6%. If the same sequence is aligned with bowtie2 means, alignment rate is 96%. why i am getting less alignment rate in tophat compared to bowtie2? Any suggestions please.

Thanks, Anitha

rna-seq alignment sequence
• 1,826 views • link
updated by Chirag Parsania Chirag Parsania • written by anitha.b • anitha.b

Is your data RNA-seq? I find it intriguing that you first tried Tophat (a splice aware aligner) and then Bowtie2 and then that the alignment rate with Bowtie2 is so high. Did you check that the Bowtie version, which TopHat is using, is the same as the one you used for 96% alignment ?

• 0 views • link
written by Santosh Anand Santosh Anand

1 answer

I dont know the reason why tophat is giving so less alignment. I would expact more aligned reads with tophat compare to bowtie. I would suggest you to use HISAT which new and alternative to tophat , developed by same group and they also recommend to use it over tophat.

• 0 views • link
written by Chirag Parsania Chirag Parsania

Log in to answer this question.

More posts like this
  • Low Alignment rate
    written by mavy • mavy
    Hello All I am working on homosapien chipset data I have done the indexing using the bowtie2 tool of the top-level assembly from the ensemble …
  • miRNA-seq alignment bowtie2
    written by imaparna27 • imaparna27
    I am trying to align miRNA-seq reads against miRBase reference using Bowtie2, but I am getting a very low overall alignment rate. using Bowtie2: > …
  • bowtie (1) 0% alignment on paired-end RNAseq data
    written by gulamaltab • gulamaltab
    Hi I am trying to align paired-end miRNA-seq data using bowtie, I have tried a number of times using different options, however, it seems like …
  • Low coverage for bases near the end of target reference sequence
    written by kspata kspata
    Hi All, I have viral samples sequenced using NextSeq PE 150 with over 90 million reads, I aligned these reads after trimming adapters and low-quality …
  • microbial RNA seq
    written by hothriananya hothriananya
    Hi, We recently extracted RNA seq data from human stool samples. We got reads that are 75bp and 150bp. These are paired-end reads. These are …
  • bowtie2 paired end insertion length
    written by jesselee516 jesselee516
    Hi, I am using bowtie2 to do alignment to reference genome. Actually, instead of getting real data, I simulated the paired-end reads(100bp) with insertion length …
  • Gene set enrichment Analysis
    written by anitha.b • anitha.b
    Hi, I am running GSEA for Gene Ontology(C5), Immunology(C7) and oncology(C6). But i am getting results as a pdfs only for Immunology and Oncology. why …
  • rna seq analysis
    written by anitha.b • anitha.b
    Hi, I have gencode long non coding(lnc) gtf and lnc reference fasta files.. if i run tophat using these two references, i am getting an …
  • alignment rate in bowtie2
    written by zizigolu zizigolu
    Hi, I used bowtie2 to map GSE69414 data sets but rate of alignment is 0 or max 0.02 How is that possible? Am I doing …
  • Bowtie2 settings in Tophat
    written by zioutisch • 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 …