[Bioperl-l] bp_genbank2gff.pl giving errors when using file....
Scott Cain
cain.cshl at gmail.com
Tue Jul 7 19:18:50 UTC 2009
Hi Gotham,
I was going to send you an email to complain to the author, until I
realized that it was me :-)
It has been quite a while since I looked at the code for this script,
as the one I typically use these days is bp_genbank2gff3.pl, but I
think I have a "fix". Try changing the name of the file to
NC_004329.gb or .gbk or .genbank. That is the (very weak) heuristic
that the script uses to determine if a file is genbank formated versus
embl (note that the error message says it's not an embl file--that's
why). If that doesn't do it, let me (and the mailing list) know.
Scott
PS: I wasn't really going to say to complain to the author directly--
that was just me trying to be funny.
PPS: As another side note, it is fairly funny to me that the code that
this script depends upon, Bio::DB::GFF::Adaptor::biofetch, says in the
documentation that it is proof-of-principle and should not be used in
production.
On Jul 7, 2009, at 1:59 PM, Gowthaman Ramasamy wrote:
>
> Hi All,
> I am trying to use bp_genbank2gff.pl script to convert a locally
> downloaded genbank file. It is throwing stack errors. But, the
> script works beautifully when I use --accession option to download
> and convert.
>
> Any suggestions? Thanks very much for checking this.
>
> the command i use:
> perl bp_genbank2gff.pl --stdout --file NC_004329.nb
>
> and i am getting the following exception message:
>
> ------------- EXCEPTION: Bio::Root::Exception -------------
> MSG: EMBL stream with no ID. Not embl in my book
> STACK: Error::throw
> STACK: Bio::Root::Root::throw /usr/lib/perl5/site_perl/5.8.5/Bio/
> Root/Root.pm:359
> STACK: Bio::SeqIO::embl::next_seq /usr/lib/perl5/site_perl/5.8.5/Bio/
> SeqIO/embl.pm:189
> STACK: Bio::DB::GFF::Adaptor::biofetch::load_from_file /usr/lib64/
> perl5/site_perl/5.8.5/x86_64-linux-thread-multi/Bio/DB/GFF/Adaptor/
> biofetch.pm:163
> STACK: bp_genbank2gff.pl:274
> -----------------------------------------------------------
>
>
> Many thanks in advance,
> Gowtham
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
-----------------------------------------------------------------------
Scott Cain, Ph. D. scott at scottcain dot net
GMOD Coordinator (http://gmod.org/) 216-392-3087
Ontario Institute for Cancer Research
More information about the Bioperl-l
mailing list