パールじゃないけど
ＵＮＩＸのコマンドでＦＡＳＴＡフォーマットのファイル中にあるシークエンスの数を数える
grep -c というコマンドを使うと数を数える事ができる。
ＦＡＳＴＡフォーマットは各シークエンスに＞で名前がついているので、
'＞'の数を数えればそのファイルの中にいくつ'＞'があるかわかる。
grep -c '>' filename
’ ’の中身を変えれば別の言葉を数える事ができる。
たとえば’No hits found’とやるとblastの結果からヒットしなかった遺伝子の数を数える事ができる。
PR
