你发的文件里面全是截图,我这边是完全看不清楚内容的,提问建议提供运行命令、输入文件截图和报错信息。
你第一次做的那个是划分成了chr的状态(蓝色框),第二个做的里面只有contig(绿框)而没有挂载到染色体上的信息,正常应该是这样:
在你没有chr(或者说是一条很大的染色体)的时候去绘图会报错,比如说你指定了需要画10条染色体但是只有1条,这是为什么共线性无法出图。
确认了酶切没有问题,是第一次没有问题还是第二次使用的没有问题?因为你第二次没有挂载成染色体。
左滑没有出现红色是指你在juice box里面滑动我下图绿色框框标出来的标尺,还是没有出现信号的意思?那这个时候你双击放大一下自己的图,再调整一下color的标尺。另外还需要注意:
1. 加载顺序:在Juicebox中,必须先加载 .hic文件,再加载 .assembly文件,才有信号;
2. 检查BAM文件。samtools flagstat就可以查看HiC数据的比对率。如果比对率低,说明数据质量可能有问题,或者基因组组装错误较多。
3. 检查有效互作对:out.links.txt如果很小说明有效互作对很少,需要重新检查HiC建库或测序数据质量
4. 检查挂载率:.agp文件里面可以看挂载情况。
