Question: Looking For Reasons Of Why The Results Of A Tophat Fusion Post Is Empty
gravatar for KS
7.1 years ago by
KS370 wrote:

Hello Everyone,

I ran tophat-2.7 on the example data with --fusion-search , it ran successfully and produced the following output files: accepted hits.bam, deletions.bed, insertions.bed, junctions.bed, fusions.out

I am trying to run tophat-fusion-post, after I run it, the output is empty. Can anyone suggest the reason for it?

The command I used to run tophat (I am retuning this command under topdir directory):

tophat -o tophatMCF7 -p 8 --fusion-search /path/to/hsapiens/bowtieindex/hg19 SRR0642861.fastq SRR0642862.fastq

The structure I am using to run tophat is the following:


blasthuman (nt and humangenomic databases)
indexes (hg19)
tophatMCF7 (output directory): accepted_hits.bam, deletions.bed, insertions.bed, junctions.bed, fusions.out, logs,, unmapped.bam
tophat fusion • 4.4k views
ADD COMMENTlink modified 5.8 years ago by Paul1.4k • written 7.1 years ago by KS370

have you checked the log files? I did not put human blast databases inside folder blasthuman. I extracted all archives under blast directory.

ADD REPLYlink written 7.1 years ago by samsara600
gravatar for nicolas.robine
7.0 years ago by
United States
nicolas.robine100 wrote:

I think I found the problem (after many errors). On line 696 of the script tophat-fusion-post, I saw:

        if string.find(file, "tophat_") != 0:

And I rename the folders tophat_ (in the above case tophat_MCF7), it seems to work! I hope that can be useful to others!

ADD COMMENTlink written 7.0 years ago by nicolas.robine100

Thanks, finally a ray of hope :-)

ADD REPLYlink written 7.0 years ago by Mikael Huss4.7k
gravatar for Mikael Huss
7.1 years ago by
Mikael Huss4.7k
Mikael Huss4.7k wrote:

You're not the only one! See these SeqAnswers threads, maybe they can provide some help:

Repost: tophat-fusion outputs empty result

tophat-fusion-post result empty

I have been unable to solve the problem in my own case, and instead moved on to other fusion finding programs like FusionCatcher, ChimeraScan, FusionMap, and FusionHunter (which I have all managed to run, although not without some trial-and-error).

ADD COMMENTlink written 7.1 years ago by Mikael Huss4.7k
gravatar for Paul
5.8 years ago by
European Union
Paul1.4k wrote:

I have got it - mcl (known gene fusions) can be obtained from previous tophatfusion version (in annotation directory),

ADD COMMENTlink modified 6 months ago by RamRS26k • written 5.8 years ago by Paul1.4k
Please log in to add an answer.


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