[Bioperl-l] SeqFeature::Annotated broken
Scott Cain
cain at cshl.org
Fri Nov 19 16:05:03 EST 2004
Allen,
I think Bio::SeqFeature::Annotated is broken, though I am not really
sure. When I run this test script:
#!/usr/bin/perl -w
use strict;
use lib '/home/scott/bioperl-live';
use Bio::FeatureIO;
my $in = Bio::FeatureIO->new(-file=> "test.gff", -format=>"gff");
while (my $feature = $in->next_feature()) {
print join ("\t",($feature->display_name,
$feature->start,
$feature->end)),"\n";
my $ac = $feature->annotation;
foreach my $key ($ac->get_all_annotation_keys() ) {
my @values = $ac->get_Annotations($key);
foreach my $value (@values) {
print print "Annotation ",$key," stringified value ",$value->as_text,"\n"; }
}
}
and a small test GFF file, I get the following output:
-------------------- WARNING ---------------------
MSG: '##sequence-region' directive handling not yet implemented
---------------------------------------------------
-------------------- WARNING ---------------------
MSG: seq_id() is deprecated, use id()
---------------------------------------------------
Can't locate object method "remove_Annotations" via package "Bio::SeqFeature::Annotated" at /home/scott/bioperl-live/Bio/SeqFeature/Annotated.pm line 174, <GEN0> line 4.
Which is weird, because the remove_Annotations method is there. Do you
have any insights as to what the problem might be?
Thanks,
Scott
--
------------------------------------------------------------------------
Scott Cain, Ph. D. cain at cshl.org
GMOD Coordinator (http://www.gmod.org/) 216-392-3087
Cold Spring Harbor Laboratory
More information about the Bioperl-l
mailing list