hmmsearch --domtblout WRKY_hmm_out.txt --cut_tc WRKY.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa 这是上课讲的hmm筛选出一个结构域，如果想筛选出两个结构域同时存在的应该怎么写脚本
可以换个hmm文件再搜索一次，得到的hmmer结果文件里面的基因取下交集就可以了；
如果觉得我的回答对您有用，请随意打赏。你的支持将鼓励我继续创作！
