将id.txt文件中不需要再跑的品系ID删除再循环跑就可以了。跑之前将前面生成的空文件删除。不过最好先看一下报错原因,为什么生成的空文件
回答于 2025-07-04 09:11
组装scaffold水平也能做,就是结果可能差点儿。 截图里面是字体样式报错,忽略就行
回答于 2025-06-27 09:45
构建泛基因列表,只需要把自己01.data_prepare时生成的longest_isoform.gff3和gene.cds.fasta文件链接过来作为数据去生成bed和cds文件就行。因为构建泛基因列表w.sh文件中的“# 提取最长的转录本”“# 提取cds”“# 改cds的名字为基因名字”步骤已经在01.data_prepare时完成了 并且不需要示例数据pangenlist里面的文件
回答于 2025-06-26 17:41
bed文件有问题吧,用“cat -A bed” 查看文件后缀是不是$。windows和linux的文字编码系统不一样,最好是将excel准备好的文件提交到Rstudio server 新建文件中进行粘贴保存。可以把结果截图发过来
回答于 2025-06-26 10:07
通过swiss-model得到模型,可以将模型和用的模板下载下来在pymol计算RMSD值。Python 得不到就用swiss-model方法,也推荐学习一下AlphaFold方法:https://zhuanlan.zhihu.com/p/699907182
回答于 2025-06-26 09:31
引物设计原则“引物应在核酸序列保守区内设计并有特异性;引物通常设计为跨内含子”,具体参考笔记:https://www.omicsclass.com/article/2227
回答于 2025-06-24 10:22
直接复制的代码,代码复制完全的话,应该没什么问题。 ${ind}.FAE1_hmmerOut.final.txt不是之前已经生成了?还有之前保留的备份吗?如果有备份就之间往下分析 这次的问题先检查这次问题涉及的文件 。用“cat -A .FAE1gene.list” 查看文件后缀是不是$。查看FAE1_hmmerOut.cut_e_len.txt 文件内容有没有问题。可以截图发过来...
回答于 2025-06-24 09:37