ALLHiC相当于Lachesis的升级版能够分型,他俩的特点是对于contig基因组拼接的准确性要求非常高,出现拼接错误很容易出现染色体聚类错误的情况,最快速的方法是用3ddna,一定能获得想要的染色体数目,其次是用ALLHiC流程中加入自带的纠错脚本这一步结果看运气,再其次是用3ddna纠错后的基因组再跑ALLHiC。虽然 3ddna排序错误率较高,但从染色体长度与基因长度来看并不影响后续基因结构预测的准确性。
您好,我按照流程处理自己的数据,在ALLHiC_partition -b sample.clean.bam -r draft.asm.fasta -e DpnII -k 15这一步,因为测序的物种是15条染色体,但运行后只有1个group出现:sample.clean.counts_GATC.15g1.txt,没有出现预期的15个group,请问是什么问题?谢谢