HMM方法

## hmmersearh进行搜索
hmmsearch --cut_tc --domtblout NBS-ABC.out NBS-ARC.hmm \\
             Arabidopsis_thaliana.TAIR10.pep.all.fa

# --cut_tc: use profile's TC trusted cutoffs to set all thresholding
# --domtblout <f>  : save parseable table of per-domain hits to file <f>

## 根据e-value 筛选高质量基因
grep -v "#" NBS-ABC.out|awk '($7 + 0) < 1E-20'|cut -f1 -d  " "|sort -u > NBS-ARC_qua_id.txt

blast

diamond makedb --in clean_as.fa --db green_cyto

makeblastdb -in gh.fa -dbtype prot -parse_seqids -title gh -out gh
blastp -query Lh_filter.fa -out lh_gh.txt -db gh -evalue 1e-5 -outfmt 6