Home

バイオ やばいぉ・・・

スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。
  • Comments (Close): -
  • TrackBack (Close): -

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:
スポンサーサイト

Home

最新記事はRSSで
コメント
最新トラックバック
リンク
ブログで苗木 グリムス
blog parts

track feed

Powered by SEO Stats

ページランク

個人サイト検索 あったかNAVI

この日記のはてなブックマーク数

フィードメーター - バイオ やばいぉ・・・

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

Page Top

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