プロトコル

[bioinfo]multiple fasta形式のデータから必要なデータのみを取り出す

まず、必要なfastaデータのIDを改行で区切ってテキストファイルに出力しておき、seqkitのgrepをwhileループで動かす。

while read line; do seqkit grep ./multiple.fasta -np "$line" -o "$line".fasta; done < id_list

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です