[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。 FASTAシークエンスのいっぱい入ったファイルの中から、ある特定のFASTAシークエンスを抽出したい時につかえるプログラムを書いてみました。 ここまで書くには結構勉強した。 FASTAフォーマット >1 ATGATATGAGGATGCGTAGTA >2 AAAAATTTTGGGGCCCCC >3 TTTTTCCCGTGTAGTGATGTGTCGTGCTGATCGTACGTCG という遺伝子が入ったファイルがある。 その中から3のみを抽出したいときは、 >3 と書かれたファイルを用意して下のプログラムを実行すると。 FASTAフォーマットのファイルの中から >3 TTTTTCCCGTGTAGTGATGTGTCGTGCTGATCGTACGTCG だけを抽出してくれる。 今5万個とか扱