multifastaから複数の配列を抜き出したい
Q&A
解決したいこと
multifastaファイル(IRGSP-1.0_2kb-upstream_2022-09-01.fasta.gz)から、ヘッダーのOsXXtXXXXXXXの文字列リストファイル(riceFREND_TOC1_coexpression_MRlist_foranalysis.txt)を用いて該当するヘッダーとシーケンスをriceFREND_TOC1_coexpression_MRl
ist_2000bpupstreem.faというファイルに抜き出したいです。
該当するソースコード
ソースコードを入力
grep -f riceFREND_TOC1_coexpression_MRlist_foranalysis.txt -A 9 IRGSP-1.0_2kb-upstream_2022-09-01.fasta.gz > riceFREND_TOC1_coexpression_MRl
ist_2000bpupstreem.fa
### 自分で試したこと
ヘッダー1行とシーケンスが9行だったので、grepでその後ろの9行まで抜き出そうとしました
0