老师我这次{ind}.FAE1_hmmerOut.final.txt都只有一个了,list文件正常,前一步生成的{ind}.FAE1gene.pep.fa也正常,这该怎么办

attachments-2025-06-HxdHqk6N685944d9c6e2d.png

attachments-2025-06-buRvlZPR685945f73420e.pngattachments-2025-06-9hRaz8Cy6859460066d9c.png我是把txt文件转成list文件然后上传到check里的,生成的{ind}.FAE1gene.pep.fa也正常,应该不是这的问题吧,这几步就这个还没检查

grep -f ./check/${ind}.FAE1gene.list ${ind}.FAE1_hmmerOut.cut_e_len.txt > ${ind}.FAE1_hmmerOut.final.txt这是用的代码,之前生成的文件都正常,今天想手动生成{ind}.domain_final.bed也不行了,这个问题还可以怎么解决么,卡在这里半个月了

请先 登录 后评论

1 个回答

rzx

直接复制的代码,代码复制完全的话,应该没什么问题。

${ind}.FAE1_hmmerOut.final.txt不是之前已经生成了?还有之前保留的备份吗?如果有备份就之间往下分析

这次的问题先检查这次问题涉及的文件 。用“cat -A .FAE1gene.list” 查看文件后缀是不是$。查看FAE1_hmmerOut.cut_e_len.txt 文件内容有没有问题。可以截图发过来。

这次的文件也可以手动生成,直接在本地打开FAE1_hmmerOut.cut_e_len.txt文件,将包含WRKYgene.list中基因的行保留下来就能得到FAE1_hmmerOut.final.txt文件

请先 登录 后评论

相似问题