我們需要使用seqkit的grep功能來實現。首先官方的語句是這樣的：
$ zcat hairpin.fa.gz | seqkit grep -f list > new.fa（
這是針對linux系統環境下。如果實在windows環境下，則要使用語句：
type non_snare.fasta|seqkit grep -f non_snareind.txt > new.fa
即：將zcat命令換成type命令，同時需注意zcat後面是壓縮的fasta檔案而type命令後面是fasta檔案。
根據id提取fasta序列
bioperl讀入寫出fasta 要求根據序列id，從fasta檔案中提取目標序列並輸出 資料序列id fasta檔案 思路以序列id為鍵，構建雜湊 用bioperl讀入fasta，獲得序列id 如果id存在於雜湊中，輸出序列 die perl 0 unless ar 3 0程式名 use bio ...
FASTA Q序列處理神器 seqkit
該軟體對於處理fasta q十分方便，省去自己編寫指令碼 1 conda install seqkit 1 取方向序列 2 seqkit seq test.fa r test re.fa34 取互補序列 5 seqkit seq test.fa p test com.fa67 取方向互補序列 8 s...
Docker根據名稱查詢容器ID映象ID並停止刪除
根據容器名稱查詢容器id並刪除 第一種寫法 docker stop docker ps a grep test project awk docker rm docker ps a grep test project awk 第二種寫法 docker stop docker ps aq filter ...