This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Strand Aware Intersect With Bx-Python

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?

python bed intersect

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.

0 answers

No answers yet.

Log in to answer this question.