同学,其实前面已经给你说过解决办法了,最好的办法就是你重新跑流程,不要乱改你的基因id,你改了之后我们也很难改流程去对你的数据进行调试,你可以把构建泛基因列表的代码这里改一下,然后回到你修改id之前重新跑一下,总之不要改动你原始的基因id,不然后面对应起来就是很麻烦,你按照流程跑: # 把第二列的前缀去掉 a...
回答于 2025-04-07 13:33
所以你在未做任何处理之前,你的gff文件里面基因ID的格式就是基因组名+基因id这样的是吗,我举个例子,比如说你的原始gff文件里面,你的基因编号就是C01_00G00000006对吗,要是你描述的这种情况,你需要把这条指令改成这样:awk 'FNR>1{print $1"\t"$3"_"$2"_"$3}' final.last.gene.pair.txt > yourfile.pan-genes.li...
回答于 2025-04-03 09:44
awk 'FNR>1{print $1"\t"$2"_"$3}' 26Pan-genes.list > 26mazi.pan-genes.list
回答于 2025-04-01 16:01
last报错应该是因为内存不足:https://www.omicsclass.com/article/1413
回答于 2025-03-31 14:10
你这应该是没有跑完,你看下结果有没有正常生成,建议你重新执行一下计算kaks的指令,另外你不用把之前的容器删除掉,你这个容器还是up的状态,直接attach重连就行: docker attach 你的容器id
回答于 2025-03-31 10:07
我们的代码里面有对文件进行正确拆分的部分: 你是从上一步完整的跑下来的吗,如果你的PAVlist和示例数据格式一样的话不应该会出现这样的问题
回答于 2025-03-27 15:47
你需要给我们提供一下明确的报错信息截图,你可以看一下发表问题的规范:https://www.omicsclass.com/article/82 因为软件报错的原因有很多种,需要有针对性解决
回答于 2025-03-27 09:21
如果你的vcf文件格式规范的话,你可以试一下vcftools:https://www.omicsclass.com/article/647
回答于 2025-03-27 09:19