Like randomintersection, but return the bp overlap instead of the number of intersecting intervals.
pybedtools
BedTool
BEDTools
pybedtools.contrib