[交流] 讨论下怎么通过gi号批量获得物种的definition 已有4人参与

#!/usr/bin/perl
use Bio::Seq;
use Bio::DB::GenBank;

$gb = new Bio::DB::GenBank;
$^I = ".bak";

while (<>) {
  $line = $_;
  if ( /gi\|(\d+)\|/ ) {
    $gi = $1;
    $seq_obj = $gb->get_Seq_by_gi ($1);
    $def = $seq_obj->desc;
  }
  $_ = $line;
  s#\t.*?$gi.*?\t#\t$def\t#;
  print;
}

2楼: Originally posted by libralibra at 2012-08-28 22:41:10
贴个例子看看

贴一个blast出来的结果(未处理的字符串)
贴一个你想要的结果(目标字符串)