Question: Kissplice2reftranscriptome error : IndexError
gravatar for jfo
9 months ago by
jfo10 wrote:


I am trying to run the last part of kissplice pipeline, which is kissplice2reftranscriptome and got this error:


kissplice2reftranscriptome -b ORP_corset_all_longest.fasta.transdeco der.bed -k kissplice_out_Fask41_coherents_type_0a.fa -t out_blat_SNP_ORP_ID_80.psl -o mainOutput.tsv -l lowQueryCOverageOutput.tsv


Opening all files and formatting tables...

Collecting data...

Traceback (most recent call last):

File "/home/user/src/anaconda2/bin/kissplice2reftranscriptome", line 4, in <module> __import__('pkg_resources').run_script('kissplice2reftranscriptomelib==1.3.2', 'kissplice2reftranscriptome')

File "/home/user/src/anaconda2/lib/python2.7/site-packages/pkg_resources/", line 748, in run_script self.require(requires)[0].run_script(script_name, ns)

File "/home/user/src/anaconda2/lib/python2.7/site-packages/pkg_resources/", line 1524, in run_script exec(script_code, namespace, namespace)

File "/home/user/src/anaconda2/lib/python2.7/site-packages/kissplice2reftranscriptomelib-1.3.2-py2.7.egg/EGG-INFO/scripts/kissplice2reftranscriptome", line 577, in <module> File "/home/usr/src/anaconda2/lib/python2.7/site-packages/kissplice2reftranscriptomelib-1.3.2-py2.7.egg/EGG-INFO/scripts/kissplice2reftranscriptome", line 223, in main

File "build/bdist.linux-x86_64/egg/kissplice2reftranscriptomelib/dataAnalysis/", line 101, in getGeneTable

IndexError: list index out of range

Additional Information:

The only deviation of my analysis from the guide/documentation is that I did not use Trinity for de novo assembly.

Thank you.

kissplice rna-seq snp k2rt • 274 views
ADD COMMENTlink modified 8 months ago by vincent.lacroix130 • written 9 months ago by jfo10
gravatar for vincent.lacroix
8 months ago by
vincent.lacroix130 wrote:

Dear user,

This bug should now be solved with version 1.3.3 of k2rt available here: When using >30 input files to KisSplice, the fasta header gets very long, and blat truncates it when creating .psl files, which causes downstream issues in k2rt. This is now fixed.

Regarding the use of another transcriptome assembler than Trinity, I think it should run fine. We only use Trinity header for filling columns 11 and 12 (SNP_in_multiple_assembled_genes, SNP_in_multiple_assembled_isoforms): c1_g1_i1 and c1_g1_i2 are interpreted as two isoforms of the same gene, while c1_g1_i1 and c2_g1_i1 are interpreted as two different genes. When using an alternative assembler, we will likely not be able to make this distinction, and most likely only one of these two columns will be interpretable.



ADD COMMENTlink modified 8 months ago • written 8 months ago by vincent.lacroix130

Thank you Vincent! Everything went smoothly now. :)

ADD REPLYlink written 8 months ago by jfo10
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: 1317 users visited in the last hour