我有一个 FASTA 文件，它故意包含一些带有错误标头（即没有>）的序列和一些带有良好标头的序列。该文件格式良好，因为核苷酸序列在一行中。
例子 ：
2865958
AACTACTACAG
>hCoV-19/2832832
ACTCGGGGGG
28328332
ATTCCCCG
>hCoV-19/2789877
ACTCGGCCC
而且我只想保留带有正确标题（即以 开头的行>）的序列，如下所示：
>hCoV-19/2832832
ACTCGGGGGG
>hCoV-19/2789877
ACTCGGCCC
我为此尝试了各种方法（sed、grep、awk）但没有正确的结果：
awk '/^>/ { ok=index($0,"hCoV")!=0;} {if(ok) print;}' combined_v4.fa > combined_v5.fa
sed -n '/^>.*hCoV/,/^>/ {/^>.*hCoV/p ; /^>/! p}' combined_v4.fa > combined_v5.fa
grep -w ">" -A 1 combined_v4.fa > combined_v5.fa
你知道怎么做吗？