#根据序列ID特点筛选对应ID序列
$ zcat hairpin.fa.gz | seqkit grep -r -p ^hsa
>hsa-let-7a-1 MI0000060 Homo sapiens let-7a-1 stem-loop
UGGGAUGAGGUAGUAGGUUGUAUAGUUUUAGGGUCACACCCACCACUGGGAGAUAACUAU
ACAAUCUACUGUCUUUCCUA
>hsa-let-7a-2 MI0000061 Homo sapiens let-7a-2 stem-loop
AGGUUGAGGUAGUAGGUUGUAUAGUUUAGAAUUACAUCAAGGGAGAUAACUGUACAGCCU
CCUAGCUUUCCU
#多次使用
$ zcat hairpin.fa.gz | seqkit grep -r -p ^hsa -p ^mmu -v
#根据ID提取序列
$ seqkit grep -f id.txt seqs.fq.gz -o result.fq.gz
# ignore case
$ seqkit grep -i -f id.txt seqs.fq.gz -o result.fq.gz
更多使用方法见：https://bioinf.shenwei.me/seqkit/usage/
