Question: Looking For Reasons Of Why The Results Of A Tophat Fusion Post Is Empty
0
gravatar for KS
6.3 years ago by
KS360
KS360 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:

top_dir:

blast
blasthuman (nt and humangenomic databases)
indexes (hg19)
mcl
refGene.txt
ensGene.txt
SRR0642861.fastq 
SRR0642862.fastq
tophatMCF7 (output directory): accepted_hits.bam, deletions.bed, insertions.bed, junctions.bed, fusions.out, logs, prep_reads.info, unmapped.bam
tophat fusion • 4.2k views
ADD COMMENTlink modified 5.0 years ago by Paul1.3k • written 6.3 years ago by KS360

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 6.3 years ago by samsara580
1
gravatar for nicolas.robine
6.2 years ago by
United States
nicolas.robine90 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 6.2 years ago by nicolas.robine90

Thanks, finally a ray of hope :-)

ADD REPLYlink written 6.2 years ago by Mikael Huss4.6k
0
gravatar for Mikael Huss
6.3 years ago by
Mikael Huss4.6k
Stockholm
Mikael Huss4.6k 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 6.3 years ago by Mikael Huss4.6k
0
gravatar for Paul
5.0 years ago by
Paul1.3k
European Union
Paul1.3k wrote:

I have got it - mcl (known gene fusions) can be obtained from previous tophatfusion version (in annotation directory), http://tophat-fusion.sourceforge.net...n-0.1.0.tar.gz

ADD COMMENTlink written 5.0 years ago by Paul1.3k
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: 1669 users visited in the last hour