fasta文件根据ID提取序列

把要提取的 序列 ID 写入 id.txt , 一行一个ID

1.seqkit

seqkit grep -f id.txt input.fa > output.fa

2.seqtk

seqtk subseq input.fa id.txt > output.fa

你可能感兴趣的:(fasta文件根据ID提取序列)