[Bioperl-l] failure of add_seqfeature
Masa Masa
watashi at post.com
Thu Mar 6 12:06:42 UTC 2008
Dear experts,
Would anybody know why the following codes generate an error of:
------------- EXCEPTION -------------
MSG: Bio::SeqFeature::Generic=HASH(0x94583c0) is not contained within parent feature, and expansion is not valid
STACK Bio::SeqFeature::Generic::add_SeqFeature /usr/lib/perl5/site_perl/5.8.0/Bio/SeqFeature/Generic.pm:767
STACK toplevel test.pl:118
--------------------------------------
15616 15693
79568 83016
=================
use Bio::Graphics;
use Bio::SeqFeature::Generic;
use Bio::SeqIO;
my $bsg = 'Bio::SeqFeature::Generic';
my $unseqfea = $bsg->new( -start=>$from[$i], -end=>$to[$i], -display_name=>'U');
for (my $i=0; $i < @from; $i++) {
print "$from[$i] $to[$i]\n";
$unseqfea->add_SeqFeature($bsg->new(-start=>$from[$i],-end=>$to[$i]));
if ($i > 10) {
exit;
}
}
--
Want an e-mail address like mine?
Get a free e-mail account today at www.mail.com!
More information about the Bioperl-l
mailing list