Hi,
I am trying to ryn geneprimp, gene prediction - validation tool. with below command line.
I am not sure if it would run for genemark.hmm.
./geneprimp -if ~/Documents/file.embl.txt -it EMBL -idp 'gene ' -idlen 5 -euk N -nr ~/ -jobid NM -tool genemark.hmm -self 'neisseria meningitidis' -jobdir ~/Documents/ -blastbin ~/bin/blastOld -bl2seqbin ~/bin/blastOld
Even I replace genemark.hmm with prodigal and the output file from tools
I get the error:
MSG: Can't see new qualifier in: rbs_spacer=None;score=1.11;cscore=-1.56;sscore=2.68;
from:
/ID=1_1;partial=10;start_type=Edge;rbs_motif=None;
rbs_spacer=None;score=1.11;cscore=-1.56;sscore=2.68;
rscore=0.00;uscore=3.18;tscore=0.00
/source="Prodigal_v2.50"
/score=1.1
STACK: Error::throw
STACK: Bio::Root::Root::throw Bio/Root/Root.pm:368
STACK: Bio::SeqIO::embl::_read_FTHelper_EMBL Bio/SeqIO/embl.pm:1294
STACK: Bio::SeqIO::embl::next_seq Bio/SeqIO/embl.pm:392
STACK: Parser:
arseGBorEMBLfile Parser.pm:263
STACK: PAR::_run_member script/geneprimp.pl:145
STACK: PAR::_run_member script/main.pl:26
STACK: PAR::import /jgi/tools/lib/perl5/site_perl/5.10.1/PAR.pm:428
STACK: __par_pl::BEGIN -e:953
STACK: script/geneprimp.pl:0
I have bioperl installed in my system.
How do I resolve this error?
I am trying to ryn geneprimp, gene prediction - validation tool. with below command line.
I am not sure if it would run for genemark.hmm.
./geneprimp -if ~/Documents/file.embl.txt -it EMBL -idp 'gene ' -idlen 5 -euk N -nr ~/ -jobid NM -tool genemark.hmm -self 'neisseria meningitidis' -jobdir ~/Documents/ -blastbin ~/bin/blastOld -bl2seqbin ~/bin/blastOld
Even I replace genemark.hmm with prodigal and the output file from tools
I get the error:
MSG: Can't see new qualifier in: rbs_spacer=None;score=1.11;cscore=-1.56;sscore=2.68;
from:
/ID=1_1;partial=10;start_type=Edge;rbs_motif=None;
rbs_spacer=None;score=1.11;cscore=-1.56;sscore=2.68;
rscore=0.00;uscore=3.18;tscore=0.00
/source="Prodigal_v2.50"
/score=1.1
STACK: Error::throw
STACK: Bio::Root::Root::throw Bio/Root/Root.pm:368
STACK: Bio::SeqIO::embl::_read_FTHelper_EMBL Bio/SeqIO/embl.pm:1294
STACK: Bio::SeqIO::embl::next_seq Bio/SeqIO/embl.pm:392
STACK: Parser:
arseGBorEMBLfile Parser.pm:263STACK: PAR::_run_member script/geneprimp.pl:145
STACK: PAR::_run_member script/main.pl:26
STACK: PAR::import /jgi/tools/lib/perl5/site_perl/5.10.1/PAR.pm:428
STACK: __par_pl::BEGIN -e:953
STACK: script/geneprimp.pl:0
How do I resolve this error?
