[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