提取可靠基因家族成员的结构域序列

运行代码如下

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

不提示报错但输出文件没有东西。
输入文件如下attachments-2024-09-Fqb0Hhj766e4fcc575c8d.png
attachments-2024-09-i16NXTCU66e4fcffa3e58.png
请先 登录 后评论

1 个回答

每天学习一点点

不应该会没有结果,你的list文件里面的基因id应该来自于你的hmmer文件,你要确保两个文件中的id对得上

请先 登录 后评论
  • 1 关注
  • 0 收藏,495 浏览
  • 树心台 提出于 2024-09-14 11:03

相似问题