Question: tophat running does not complete properly.
0
gravatar for yff
22 months ago by
yff0
yff0 wrote:

Hello everyone,

I am trying to use tophat to align my fastq reads to the reference genome.

I first used bowtie to create index files. I have trimmed my read1 and read2 fastq files. All output files of bowtie index and read files are in the same directory. Here you can see the screen shot of my directory.

enter image description here

Then I am running the following command to run tophat.

tophat -p 8 -G genome.gtf -o sample_thout genome SAMPLE_R1.fq SAMPLE_R2.fq

I do not know what I am doing wrong since my command runs but after almost 45 minutes it stops.

enter image description here

The sample_thout directory only contains the following files.

enter image description here

Inside the tmp folder I can see the following:

enter image description here

inside the logs folder I can see the following files:

enter image description here

So I can not run tophat properly to get the desired output. I am pretty new to RNA seq and using tophat. Does anybody have any idea what is my mistake?

bowtie rna-seq alignment tophat • 641 views
ADD COMMENTlink written 22 months ago by yff0

Please run STAR or HISAT2 or so. Tophat is not recommended these days.

ADD REPLYlink written 22 months ago by Jeffin Rockey1.1k

Thanks for your reply. Actually I have to use tophat because it is part of my research pipeline. Do you have any idea what is the issue?

ADD REPLYlink written 22 months ago by yff0
1

OK. What is the genome given in the command line?

tophat -p 8 -G genome.gtf -o sample_thout **genome** SAMPLE_R1.fq SAMPLE_R2.fq
ADD REPLYlink modified 22 months ago • written 22 months ago by Jeffin Rockey1.1k

The genome is my index which is called ToxoDB35. I mean that is the name for out out of bowtie index.

The exact command that I am running is as follows:

tophat -p 8 -G ToxoDB35.gtf -o sample_thout ToxoDB35 SAMPLE_R1.fq SAMPLE_R2.fq
ADD REPLYlink modified 22 months ago • written 22 months ago by yff0

DON'T USE TOPHAT - LIOR PACHTER TWEET

Hello yasaman.rezvani66,

You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.

ADD REPLYlink written 22 months ago by lakhujanivijay5.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1518 users visited in the last hour