FC2ブログ

Home > BioPerl > BioPerlでFASTAファイルから指定したIDの配列を抜き出す

BioPerlでFASTAファイルから指定したIDの配列を抜き出す

スポンサード リンク
BioPerlでFASTAファイルから指定したIDの配列を抜き出すことをしたいと思います。
use Bio::DB::Fastaを使って、FASTAファイル中から指定したIDの配列をFASTA形式で出力させます。

use Bio::SeqIO;
use Bio::DB::Fasta;
$seq_out = Bio::SeqIO->new(-fh => \*STDOUT, '-format' => 'fasta');
$file = "sample.fasta";
$id = $ARGV[0];
$db = Bio::DB::Fasta->new($file);
$seqobj = $db->get_Seq_by_id($id);
$seq_out->write_seq($seqobj);


解説
01:use Bio::SeqIO;
02:use Bio::DB::Fasta;
03:$seq_out = Bio::SeqIO->new(-fh => \*STDOUT, '-format' => 'fasta');
04:$file = "sample.fasta";
05:$id = $ARGV[0];
06:$db = Bio::DB::Fasta->new($file);
07:$seqobj = $db->get_Seq_by_id($id);
08:$seq_out->write_seq($seqobj);

01:~02:おまじない
03:
04:
05:
06:
07:
08:

Comments:0

コメント欄

Trackback:0

この記事へのTrackBack URL
http://bioinformatics.blog7.fc2.com/tb.php/485-01da0af2
BioPerlでFASTAファイルから指定したIDの配列を抜き出す from バイオ やばいぉ・・・

Home > BioPerl > BioPerlでFASTAファイルから指定したIDの配列を抜き出す

あわせて読みたいブログパーツ

Page Top