所以你在未做任何处理之前,你的gff文件里面基因ID的格式就是基因组名+基因id这样的是吗,我举个例子,比如说你的原始gff文件里面,你的基因编号就是C01_00G00000006对吗,要是你描述的这种情况,你需要把这条指令改成这样:awk 'FNR>1{print $1"\t"$3"_"$2"_"$3}' final.last.gene.pair.txt > yourfile.pan-genes.li...
回答于 1天前
awk 'FNR>1{print $1"\t"$2"_"$3}' 26Pan-genes.list > 26mazi.pan-genes.list
回答于 3天前
你这应该是没有跑完,你看下结果有没有正常生成,建议你重新执行一下计算kaks的指令,另外你不用把之前的容器删除掉,你这个容器还是up的状态,直接attach重连就行: docker attach 你的容器id
回答于 4天前
我们的代码里面有对文件进行正确拆分的部分: 你是从上一步完整的跑下来的吗,如果你的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
“发现P-Value的原始文件有NaN的情况”,应该是你有些SNP位点缺失率过高之类的情况,大概是你前期的数据过滤不是很严格,你画图的时候直接删掉这样的数据就可以 关于你提出的曼哈顿图上面的点显示不全的问题,你可以对脚本进行如下的修改: scale_y_continuous(expand = c(0,0.05) 即把脚本中的y轴坐标轴上限增加5%
回答于 2025-03-19 16:32