This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tophat2 Bowtie2 inchoerent results?

I am trying to align paired-end data from Illumina (300bp x 2) to a draft genome. First of all I use Trimmomatic to remove adapters and bad quality ends. later I tried to map the reads; first with Bowtie2 alone and later with TopHat2 (which uses Bowtie2).

Bowtie2 alone mapped about 68% of the reads while Tophat about 40%, with which parameters is TopHat calling Bowtie that the alignment rate decrease instead of increasing? Since the reads are quite long (up to 300nuc) should I modify the --read-edit-dist and -N parameters (in this way I get more percentage of mapped reads, but I don't know if I'm getting ad alignments)?

Thank you very much

tophat ngs bowtie2 rna-seq

AFAIK, tophat2 uses more strict default settings. For example tophat2 allows only 3 mismatches with default settings. There is a lot more stuff different but I can't tell from the top of my head. Therefore just try some different settings!

Best

0 answers

No answers yet.

Log in to answer this question.