/home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/motif.py:126: RuntimeWarning: divide by zero encountered in log2
  ic=2+np.sum(p*np.nan_to_num(np.log2(p)), axis = 1)
/home/jvierstra/.local/miniconda3/envs/py3.9_default/lib/python3.9/site-packages/descartes/patch.py:63: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.
  concatenate([asarray(t.exterior)[:, :2]] +