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