INNER CODE UNIT · Python

find_intersect

cbg-ethz/V-pipe · workflow/scripts/alignmentIntervals.py:179

    def find_intersect(intervals_new, intervals_old):
        intervals_intersected = []

        for i_new in intervals_new:
            for i_old in intervals_old:
                lower_bound = max(i_old[0], i_new[0])
                upper_bound = min(i_old[1], i_new[1])

                if lower_bound < upper_bound:
                    intervals_intersected.append((lower_bound, upper_bound))

        return intervals_intersected

    idx = 0
    with open(input_file, "r") as infile:
        for line in infile:
            record = line.rstrip().split("\t")
            key = record[0]

View source record →