bioerl 获取gi号

代码示例:

use Bio::DB::EUtilities;

my @ids = qw(CAB02640 EAS10332 YP_250808 NP_623143 P41007);

my $factory = Bio::DB::EUtilities->new(-eutil => 'efetch',
                                       -db      => 'protein',
                                       -id      => \@ids,
                                       -email   => 'mymail@foo.bar',
                                       -rettype => 'gi');

my @gis = split(m{\n},$factory->get_Response->content);

print join(',',@gis), "\n";

参考资料:

    http://bioperl.org/howtos/EUtilities_Cookbook_HOWTO.html

 

posted on 2017-07-05 19:00  庐州月光  阅读(283)  评论(0编辑  收藏  举报