Home > スポンサー広告 > BioPerlでFASTAファイルから指定したIDの配列を抜き出す

スポンサーサイト

スポンサード リンク
上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

Comments:-

コメント欄

Trackback:-

この記事へのTrackBack URL
http://bioinformatics.blog7.fc2.com/tb.php/485-01da0af2
スポンサーサイト from バイオ やばいぉ・・・

Home > スポンサー広告 > BioPerlでFASTAファイルから指定したIDの配列を抜き出す

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

上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。