Thankyou again for a great program, soon to be cited in my paper!
UGENE produces gb files of hits found, for instance when searching for repats or pattern matches. I've been trying to further treat these results by parsing them with bio-perl/python, but the parsers of these two languages are not compatible with the kind of genbank output I get from UGENE.
I notice now I'm using 1.4.1 (on windows, 1.4 gHz single processor on tc4200 laptop) not the most recent verison, but I can't see anything about this problem in the later bug fixes.
The perl code that produced the latest error looked like this:
Code:use Bio::SeqIO;
### This is a file from UGENEs find pattern tool:####
my $filenme="C:/Documents and Settings/Carl/My Documents/NewBioInf/myGB/hitsOniS3.gb";
$seqio_obj = Bio::SeqIO->new(-file => $filenme, -format => "genbank" );
$seq_obj = $seqio_obj->next_seq;
And here's the error:
Code:C:\Documents and Settings\Carl>perl readGB.pl
--------------------- WARNING ---------------------
MSG: Unknown alphabet:
---------------------------------------------------
--------------------- WARNING ---------------------
MSG: Unexpected error in feature table for Skipping feature, attempting to recover
---------------------------------------------------
------------- EXCEPTION -------------
MSG: Alphabet '1' is not a valid alphabet ('dna','protein','rna') lowercase
STACK Bio::PrimarySeq::alphabet C:/Perl/site/lib/Bio/PrimarySeq.pm:571
STACK Bio::PrimarySeq::new C:/Perl/site/lib/Bio/PrimarySeq.pm:208
STACK Bio::Seq::new C:/Perl/site/lib/Bio/Seq.pm:484
STACK Bio::Seq::RichSeq::new C:/Perl/site/lib/Bio/Seq/RichSeq.pm:110
STACK Bio::Seq::SeqFactory::create C:/Perl/site/lib/Bio/Seq/SeqFactory.pm:116
STACK Bio::Factory::ObjectFactoryI::create_object C:/Perl/site/lib/Bio/Factory/ObjectFactoryI.pm:102
STACK Bio::Seq::SeqBuilder::make_object C:/Perl/site/lib/Bio/Seq/SeqBuilder.pm:337
STACK Bio::SeqIO::genbank::next_seq C:/Perl/site/lib/Bio\SeqIO\genbank.pm:717
STACK toplevel readGB.pl:5
-------------------------------------
Best regards,
Theo Vigil