3eea3791dc36fbfc521b45b9eebd854c73908b62,dataPrepScripts/CountNumInBed.py,,Calc,#Any#,15

Before Change

    logging.info("Loading BED file ...")
    tree = {}
    with open(args.bed_fn) as f:
        for row in f:
            row = row.strip().split()
            name = row[0]
            if name not in tree:
                tree[name] = intervaltree.IntervalTree()
            begin = int(row[1])
            end = int(row[2])
            tree[name].addi(begin, end)

    logging.info("Counting number of records in bed regions ...")

    a = 0
    o = 0

After Change

    logging.info("Loading BED file ...")
    tree = {}
    f = subprocess.Popen(shlex.split("gzip -fdc %s" % (args.bed_fn) ), stdout=subprocess.PIPE, bufsize=8388608)
    for row in f.stdout:
        row = row.strip().split()
        name = row[0]
        if name not in tree:
            tree[name] = intervaltree.IntervalTree()
        begin = int(row[1])
        end = int(row[2])
        tree[name].addi(begin, end)
    f.stdout.close()
    f.wait()

    logging.info("Counting number of records in bed regions ...")

Instances