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文件。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章