[Bioperl-l] Error when using Bio::AlignIO on stockholm file from Pfam
Peng Yu
pengyu.ut at gmail.com
Thu Jan 3 22:41:37 UTC 2013
Hi,
I tried the following code on the Pfam stockholm file (downloaded at
ftp.sanger.ac.uk/pub/databases/Pfam/current_release/Pfam-A.full.gz).
But it gives me the following error. Is it a bug in bioperl or it is
an error in the Pfam-A.full.gz file?
~/linux/bin/xplat/src/stk/stk2acc/perl$ cat main.pl
#!/usr/bin/env perl
use strict;
use warnings;
#use FindBin;
#use lib "$FindBin::Bin/.";
use Bio::AlignIO;
use Data::Dumper;
my $in = Bio::AlignIO->new(
-format => 'stockholm',
-file => '/dev/stdin'
);
while(my $aln=$in->next_aln) {
print $aln->accession, "\n";
}
~/linux/bin/xplat/src/stk/stk2acc/perl$ cat main_full.sh
#!/usr/bin/env bash
zcat ~/db/Pfam-A.full.gz | ./main.pl
~/linux/bin/xplat/src/stk/stk2acc/perl$ ./main_full.sh
...
------------- EXCEPTION -------------
MSG: '1985D' is not an integer.
STACK Bio::Range::end /Users/py/Library/ActivePerl-5.14/lib/Bio/
Range.pm:227
STACK Bio::Annotation::Target::new /Users/py/Library/ActivePerl-5.14/
lib/Bio/Annotation/Target.pm:81
STACK Bio::AlignIO::Handler::GenericAlignHandler::_stockholm_target /
Users/py/Library/ActivePerl-5.14/lib/Bio/AlignIO/Handler/
GenericAlignHandler.pm:293
STACK Bio::AlignIO::Handler::GenericAlignHandler::data_handler /Users/
py/Library/ActivePerl-5.14/lib/Bio/AlignIO/Handler/
GenericAlignHandler.pm:73
STACK Bio::AlignIO::stockholm::next_aln /Users/py/Library/
ActivePerl-5.14/lib/Bio/AlignIO/stockholm.pm:470
STACK toplevel ./main.pl:17
-------------------------------------
More information about the Bioperl-l
mailing list