Strand Aware Intersect With Bx-Python

Question

[Egl]

I've played a little with the bx-python library. This example script prints every region/line in fila_a overlapping at least one region in file_b. Is it possible to modify this script to also check if the strands of the overlapping regions are equal?

[Brad Chapman]

An easy way to do this is to subset the BED files into positive and negative strand files and do the intersections with those independently.