[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