fasta文件根据ID提取序列
2022-04-12 本文已影响0人
ayunga
把要提取的 序列 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
把要提取的 序列 ID 写入 id.txt , 一行一个ID
seqkit grep -f id.txt input.fa > output.fa
seqtk subseq input.fa id.txt > output.fa