Off topic:Tophat-Fusion-Post Error
1
3
Entering edit mode
11.8 years ago
Richard ▴ 590

I ran tophat2 with the following command:

tophat-2.0.3.Linux_x86_64/tophat -p 6 --fusion-search -o outDir/ hg19 hg19.r1.fq hg19.r2.fq

I got all the usual tophat output files, including the fusions.out file. Next I want to run tophat-fusion-post.

I set up my folders as directed:

ls $
blast_human
ensGene.txt
mcl
refGene.txt
tophat_sim01
tophat_sim03

When I try and run my command I get a criptic error:

tophat-2.0.3.Linux_x86_64/tophat-fusion-post --output-dir tophatfusion_out ../tophat_fusion_refs/hg19
[Tue Jun 19 14:36:11 2012] Beginning TopHat-Fusion post-processing run (v2.0.3)
-----------------------------------------------
[Tue Jun 19 14:36:11 2012] Extracting 23-mer around fusions and mapping them using Bowtie
[Tue Jun 19 14:36:11 2012] Filtering fusions
        Processing: tophat_sim01/fusions.out
Traceback (most recent call last):
  File "/home/rcorbett/aligners/tophat2/tophat-2.0.3.Linux_x86_64/tophat-fusion-post", line 2083, in ?
    sys.exit(main())
  File "/home/rcorbett/aligners/tophat2/tophat-2.0.3.Linux_x86_64/tophat-fusion-post", line 2054, in main
    filter_fusion(bwt_idx_prefix, params)
  File "/home/rcorbett/aligners/tophat2/tophat-2.0.3.Linux_x86_64/tophat-fusion-post", line 695, in filter_fusion
    filter_fusion_impl(fusion_file, refGene_list, ensGene_list, seq_chr_dic, fusion_gene_list)
  File "/home/rcorbett/aligners/tophat2/tophat-2.0.3.Linux_x86_64/tophat-fusion-post", line 422, in filter_fusion_impl
    kmer_len = len(seq_chr_dic.keys()[0])
IndexError: list index out of range

In the output folder I have

ls $
fusion_seq.bwtout
fusion_seq.fa
fusion_seq.map
sample_list.txt

but only the sample_list.txt and fusion_seq.fa are non-empty files. Can anyone provide with with some help?

• 8.1k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2649 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6