What libraries are there for handling logic on intervals? I really want all the set operations on the intervals at genome scale. Bonus points if it has some persistence mechanism and can track the ancestor of an operation (i.e. interval X is the result of the intersection of Y and Z).
I looked at bx python, but the documentation is not great.