Question: troubleshooting: error on run of kissplice2refgenome
0
gravatar for celine
23 months ago by
celine0
celine0 wrote:

I'm following a modified version of the pipeline outlined at http://kissplice.prabi.fr/tools/kiss2rt/ Except, that I have a reference genome.

Thus, I'm using kissplice2refgenome (and also I skip the trinity/transdecoder steps).

I've generated my .psl file through blat, and downloaded the annotations file in gtf format from http://genome.ucsc.edu/cgi-bin/hgTables.

With the following command,

kissplice2refgenome myPSLfile.psl myAnnotationsFile.gtf -o myOutputFile

I get the following error:

Run starts...
Reading annotations file...
Processing mapping information...
Traceback (most recent call last):
  File "/usr/local/bin/kissplice2refgenome", line 4, in <module>
    __import__('pkg_resources').run_script('kissplice2refgenome==1.0.0', 'kissplice2refgenome')
  File "/usr/lib/python2.7/dist-packages/pkg_resources/__init__.py", line 719, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/usr/lib/python2.7/dist-packages/pkg_resources/__init__.py", line 1511, in run_script
    exec(script_code, namespace, namespace)
  File "/usr/local/lib/python2.7/dist-packages/kissplice2refgenome-1.0.0-py2.7.egg/EGG-INFO/scripts/kissplice2refgenome", line 140, in <module>

  File "/usr/local/lib/python2.7/dist-packages/kissplice2refgenome-1.0.0-py2.7.egg/EGG-INFO/scripts/kissplice2refgenome", line 119, in main

TypeError: getInfoFile() takes exactly 18 arguments (16 given)

Thanks in advance for your help!

snp kissplice2refgenome • 819 views
ADD COMMENTlink modified 22 months ago by leandro.ishi.lima80 • written 23 months ago by celine0

Some more information:

My .psl file was generated with BLAT (I couldn't get the other options to work)

My annotations was downloaded from http://genome.ucsc.edu/cgi-bin/hgTables (http://bit.ly/2iDUiFB)

ADD REPLYlink modified 22 months ago • written 22 months ago by celine0
0
gravatar for leandro.ishi.lima
22 months ago by
leandro.ishi.lima80 wrote:

Hello,

Maybe this is a command-line error, I am not sure. Could you please try the following run:

kissplice2refgenome -a myAnnotationsFile.gtf -o myOutputFile myPSLfile.psl

If this still does not work, maybe it is a good idea to try with a .sam file. Transform your psl to sam using samtools/misc/psl2sam.pl (https://github.com/samtools/samtools) and try:

kissplice2refgenome -a myAnnotationsFile.gtf -o myOutputFile mySAMfile.sam

If it does not work even like this, please tell us.

Kind regards.

ADD COMMENTlink written 22 months ago by leandro.ishi.lima80
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: 1407 users visited in the last hour