SeqKit根据ID提取序列

我们需要使用SeqKit的grep功能来实现。首先官方的语句是这样的：

$ zcat hairpin.fa.gz | seqkit grep -f list > new.fa（https://bioinf.shenwei.me/seqkit/usage/#seqkit）

这是针对Linux系统环境下。如果实在Windows环境下，则要使用语句：

TYPE non_snare.fasta|seqkit grep -f non_snareind.txt > new.fa

即：将zcat命令换成TYPE命令，同时需注意zcat后面是压缩的fasta文件而TYPE命令后面是fasta文件。