Tophat alignments is low but STAR is high.
If I used Tophat to align reads to human genome and I got ~72% alignment rate. Then I realigned the reads using STAR and around 90% of the reads were mapped. I also used the same commands of tophat to align another dataset and got 90% alignments rate. So I don't think I am not using tophat incorrectly.
I just used trimgalore before alignment . This confuses me. Can we have explanation?
tophat -p 4 --no-mixed -G genes.gtf --library-type fr-unstranded -o out genome reads_val_1.fq.gz reads_val_2.fq.gz
2 answers
Tophat is an inferior and antiquated tool. Its inferiority has been born out in countless comparisons. You should never use tophat2 for new projects.
• 0 views
•
link
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
STAR Genome indexing (Homo_sapiens_assembly38.fasta vs. GRCh38.primary_assembly.genome.fa)
written by bsh •
-
featureCounts has low successfully assigned reads
written by tomas4482 •
-
Getting unmapped paired-end reads as fastq from unmapped.bam obtained with tophat2
written by cris.kgs •
-
Short Insert size and aggressive trimming
written by Sharon
-
RNA-seq read alignment
written by amy16 •
-
Tophat vs Tophat-fusion concordant alignment rate
written by IP
-
Cufflinks with STAR BAM and Tophat BAM produces different results
written by RT •
-
tophat2 aligment ratio reduced when increase insert size
written by 52Teth •
-
Tophat: Proper-Paired Reads
written by ai •
-
Strand Specific Tophat Output
written by GPR
