[Bioperl-l] Write features in a gbk

Antony03 antony.vincent.1 at ulaval.ca
Thu Mar 6 01:36:59 UTC 2014


Hi,

I want to write some features to generate a new gbk file. I try like this:
#!/usr/bin/perl


 
use Bio::Seq;
use Bio::SeqIO;
use Bio::SeqFeature::Generic;


$feat = new Bio::SeqFeature::Generic(-start   => 40,
                                       -end     => 80,
                                       -strand  => 1,
                                       -primary => 'exon',
                                       -source  => 'internal' );

$seqobj = Bio::SeqIO->new(-file => '>test.gbk', -format => 'genbank');
$seqobj->add_SeqFeature($feat); # Add the SeqFeature to the Seq object

But I got this error:
Can't locate object method "add_SeqFeature" via package
"Bio::SeqIO::genbank" at ./annotation.pl line 17.

Thank!





--
View this message in context: http://bioperl.996286.n3.nabble.com/Write-features-in-a-gbk-tp17342.html
Sent from the Bioperl-L mailing list archive at Nabble.com.



More information about the Bioperl-l mailing list