根据基因ID提取序列

seqkit grep -f bwa.txt pilon.fasta > filtered_sequences.fasta
seqkit grep -v -f bwa.txt pilon.fasta > filtered_sequences.fasta ##-v则是剔除文件中的序列

提取以字母开头的序列：

seqkit grep -r -p "^[A-Za-z]" sequences.fa > extracted_sequences.fa