makeblastdb -in ref.nbs.plant.fa -dbtype prot -out blastdb
blastp -num_threads 20 -db blastdb -query Arabidopsis_thaliana.TAIR10.pep.all.fa -outfmt 7 -seg yes > blastp.out &
cat blastp.out |awk '$3>75' |cut -f1 |sort -u > blastp_result_id.list
最后一步筛选ident>75的序列
comm -12 blastp_result_id.list final.NBS.list > common.list
这样可以取两个文件的交集，下载序列。
less Arabidopsis_thaliana.TAIR10.pep.all.fa.gz | /home/jun_qzt/miniconda3/envs/rna-seq/bin/seqkit grep -f common.list > final_searh_NBS-ARC_qua.fa
cut -f1 代表提取第一行
sort -u u代表唯一性
