[Bioperl-l] Bio::DB::SeqFeature spliced_seq()

Chris Fields cjfields at illinois.edu
Mon Mar 22 18:56:30 UTC 2010


I have just noticed that spliced_seq() is borked with
Bio::DB::SeqFeature and am thinking about implementing it.  Or is
similar functionality already implemented elsewhere?

Currently, it is calling entire_seq(), which I plan on avoiding simply
to prevent sucking in the entire sequence into memory.  This is
currently what happens:

<partial code>
---------------------------

my $it = $store->get_seq_stream(-type => 'mRNA');

my $ct = 0;
while (my $sf = $it->next_seq) {
    my $seq = $sf->spliced_seq;  # dies with exception
}

---------------------------

------------- EXCEPTION: Bio::Root::NotImplemented -------------
MSG: Abstract method "Bio::SeqFeatureI::entire_seq" is not implemented
by package Bio::DB::SeqFeature.
This is not your fault - author of Bio::DB::SeqFeature should be blamed!

STACK: Error::throw
STACK:
Bio::Root::Root::throw /home/cjfields/bioperl/live/Bio/Root/Root.pm:368
STACK:
Bio::Root::RootI::throw_not_implemented /home/cjfields/bioperl/live/Bio/Root/RootI.pm:739
STACK:
Bio::SeqFeatureI::entire_seq /home/cjfields/bioperl/live/Bio/SeqFeatureI.pm:325
STACK:
Bio::SeqFeatureI::spliced_seq /home/cjfields/bioperl/live/Bio/SeqFeatureI.pm:458
STACK: beestore.pl:17
----------------------------------------------------------------



chris




More information about the Bioperl-l mailing list