I’m using this code for my trimmed paired-end RNA-Seq libraries
tophat2 -o outdir --mate-inner-dist 60 --mate-std-dev 50 -p 6 --GTF <annotation_file.gtf> -g 2 --library-type fr-secondstrand --no-discordant --no-mixed <bowtie2_index> file_trimmed_1.fastq file_trimmed_2.fastq
and prepare the count table with featuresCount from bam file for differential expression analysis.
Is it a correct way?
Thanks
1 answer
Unless you have good reasons to do so, you should run most tools with default parameters. In addition, you should also know that TopHat is in low maintenance and considered deprecated. It still works, but it's not considered the state of the art spliced aligner anymore. For that, have a look at HISAT2 or STAR.
Log in to answer this question.
-
Weird results in TopHat2 paired-end alignment
written by Aspire
-
tophat fusion error
written by dimitrischat
-
Getting unmapped paired-end reads as fastq from unmapped.bam obtained with tophat2
written by cris.kgs •
-
RNA-seq High Variation Between and Within Groups
written by jacobsen.jeremy •
-
Different version of Tophat gives me different number of mapped reads
written by BioDH •
-
error while running tophat ReadStream::getRead() called with out-of-order id#!
written by Kritika
-
Tophat2 and Bowtie2 : determining appropriate values for mate-inner-dist and mate-std-dev
written by irritable_phd_syndrome •
-
How to handle this error:tophat_reports,about tophat fusion
written by reany •
-
Pipeline for >100 RNA-Seq PE samples
written by umn_bist •
-
Error RNA-seq data, TOPHAT
written by David_emir
