Tophat2 Bowtie2 inchoerent results?
0
0
Entering edit mode
8.0 years ago
Bio_ysl ▴ 20

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 • 2.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1528 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6