getseqfromNCBIbygi

本脚本通过Perl实现,能根据输入的GI号列表从NCBI批量下载对应的基因序列,并将序列信息以GenBank格式保存。同时,脚本记录了成功下载的GI号。

#! /usr/bin/perl -w
#get sequences from NCBI according to gi
use strict;
use Bio::Perl;
use Bio::SeqIO;
use Bio::Seq;

open(GIIN,"$ARGV[0]") or die "Cannot open file";
open(GIED,"> $ARGV[1]") or die "cannot create file";

my @allgi=<GIIN>; #read all of gies
my $seqout = new Bio::SeqIO(-file =>'>> nameofafile', -format => "genbank");

foreach my $allgi(@allgi){
    my $seq = get_sequence('genbank', $allgi);
    if($seq){
        $seqout->write_seq($seq);    
        my $gied = $allgi;
        print GIED "$gied"; #put gies of downloaded sequences into a file
    }
}

close GIIN;
close GIED;

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值