使用sklearn

使用bedtools intersect -v 找B.bed未覆盖A.bed的位置。输出结果里面，理论上A的总长度应该等于B的总长度加上新生成的C的总长度，由于结果总是不对，去官方文档看了一下，发现原来-v参数不能做到我想要的操作。

新硬盘组raid是怎么回事呢？

想使用命令行来操作IGV，然后输出截图。但是好像只有新版的[igvtools](https://software.broadinstitute.org/software/igv/igvtools_commandline)能做到，由于igvtools依赖java 11，而服务器中java版本是java 8，为了避免影响环境变量因此放弃使用。

这里以稀释曲线的方式求香农指数。可以查看结果是否到达平台期来判断香农指数是否达到最大。