48 个回答
T2T基因组调研图中jellyfish -h和genomescope.R -m 都是设置的10...
设置依据是前面那一步，jellyfish stats -o illu.stat illu，这一步的结果illu.stat里面有个Max_count，-h和-m设置的比这个数字大就可以了。不同样本的差异根据这个统计结果判断。
回答于 5天前
verkko组装后怎么生成.hic和.assembly文件
verkko组装结果里没有这个，需要你自己用HiC数据比对到基因组上，然后生成.hic和.assembly文件
回答于 2025-04-14 14:22
关于填充gap后的polish
polish用的是什么软件呢？只使用了三代数据进行抛光吗？抛光完成之后和参考基因组或者和之前的组装做比对的话，共线性差异大吗？基因组减小了多少呢？可以采用racon对7个gap的组装重新做一下抛光。
回答于 2025-04-14 14:05
老师，我在做dadi分析的时候，发现没有代码里显示的easySFS.py的...
easySFS.py是easySFS这个软件里的： https://github.com/isaacovercast/easySFS
回答于 2025-04-01 10:45
老师我在构建泛基因列表，合成last.txt文件，结果只有两列
是的，经过python处理之后的结果确实只有两列。你可以不做python3 $script/find_1vs1_v2.py last.txt 1vs1.txt，直接把现在的last.txt作为1vs1.txt继续往下走就可以了。 如果想要执行python3 $script/find_1vs1_v2.py last.txt 1vs1.txt需要多几步的linux操作，不是很有必要。
回答于 2025-03-07 09:08
老师我在做泛基因家族分析里面的泛基因列表构建，生成last.txt文...
不需要执行这么复杂的操作，“cat ../01.jcvi/pangene_filter.*.last | grep -v "^#" > last.txt”，这一步是为了合并所有的比对文件并且删除掉"#"开头得行，主要是为了下一步去执行python的脚本。那你不合并，直接从原来的last文件里面删除 "#"开头的行，然后执行python的脚本就可以了。 for i in `find ../01.jcvi -nam...
回答于 2025-03-06 09:54