请教Bio::SeqIO 写入文件

cyagen

请教一个BioPerl问题

目的： 有一条生物序列，包括序列、序列名称，以及feature信息，要把这些内容以genbank格式保存到文件。
代码如下：

1. use Bio::SeqIO;
   
2. use Bio::Seq;
   
3. use Bio::SeqFeature::Generic;
   
4. use strict;
   
5. 
   
6. my $seq = "AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA";
   
7. my $name = "test";
   
8. my $Obj_Seq = Bio::Seq->new( -display_id => $name, -seq => $seq);
   
9. my $feat = Bio::SeqFeature::Generic->new(
   
10.                 -start =>1, -end =>10, -strand =>1,
    
11.                 -primiary =>'CDS', -display_name =>'misc');
    
12. $Obj_Seq->add_SeqFeature($feat);
    
13. my $seq_io = Bio::SeqIO->new(-format => 'genbank', -file => "test.gb");
    
14. $seq_io->write_seq($seq_obj);
    

复制代码

但是test.gb内容为空，没有将序列写到文件中。
请问如何解决这个问题？
谢谢！

longbow0

1. 
   
2. my $seq_io = Bio::SeqIO->new(-format => 'genbank', -file => ">test.gb");
   

复制代码