Problems with refseq annotation when running Islandpath-DIMOB
0
0
Entering edit mode
4.3 years ago
hugo.avila ▴ 490

Hi, I'm trying to run Islandpath-DIMOB (conda version=1.0.4) on refseq genomes and I get this error massage on some files (.gbff):

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: start [5340971] is greater than end [1335]. 
If you want to truncated and reverse complement, 
you must call trunc followed by revcom. Sorry.
STACK: Error::throw
STACK: Bio::Root::Root::throw /home/user/miniconda3/envs/IsL/lib/site_perl/5.26.2/Bio/Root/Root.pm:447
STACK: Bio::PrimarySeqI::trunc /home/user/miniconda3/envs/IsL/lib/site_perl/5.26.2/Bio/PrimarySeqI.pm:456
STACK: GenomeUtils::read_and_convert /home/user/miniconda3/envs/IsL/opt/islandpath/lib/GenomeUtils.pm:422
STACK: /home/user/miniconda3/envs/IsL/bin/islandpath:129
-----------------------------------------------------------

Apparently, the error is caused by the fact that some orfs have an start coordinate value greater than the end in the .gbff file. I work around the problem by re-annotating all genomes with Prokka, but the process is resource intensive with larger data sets. I need a way to directly repair .gbff files (start < end) or fix Islandpath-dimob to allow these refseq annotations.

software error bacteria Prokka Islandpath-dimob • 956 views
ADD COMMENT

Login before adding your answer.

Traffic: 1891 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