linux中awk命令从fasta中提取指定的scaffold

 

001、

[root@pc1 test4]# ls
test.fa
[root@pc1 test4]# cat test.fa    ## 测试fasta文件
>mmm
11111
>bbb
888888
7777
>kkk
6666666
99999
>qqq
0000
33333
[root@pc1 test4]# awk '/>/ {p=($0 ~ "kkk")} p' test.fa
>kkk
6666666
99999
[root@pc1 test4]# awk '/>/ {p=($0 ~ "qqq")} p' test.fa   
>qqq
0000
33333