运行代码如下
cat ../id.txt|while read ind;do
grep -f ./check/${ind}.WRKYgene.list ${ind}.WRKY_hmmerOut.cut_e_len.txt > ${ind}.WRKY_hmmerOut.final.txt
done
不应该会没有结果,你的list文件里面的基因id应该来自于你的hmmer文件,你要确保两个文件中的id对得上