view stack-histogram.py @ 13:c4bde687c846 draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 0cfd88d16b0b20a66d4b21b037ddee6a8c440d88
author iuc
date Wed, 12 Apr 2023 15:09:05 +0000
parents 4b519282a05b
children
line wrap: on
line source

#!/usr/bin/env python
import logging
import sys

import pyBigWig

logging.basicConfig(level=logging.INFO)
log = logging.getLogger()

if __name__ == "__main__":
    files = sys.argv[1:]
    bws = [pyBigWig.open(x) for x in files]

    # obtain some chroms. Hope all sets are identical!
    k = list(bws[0].chroms().keys())

    # do magic?
    # nah.
    # just assert that intervals are identical.
    # and crash otherwise.
    # sorry not sorry.

    for chrom in k:
        for interval_set in zip(*[bw.intervals(chrom) for bw in bws]):
            (start, end) = interval_set[0][0:2]
            values = ",".join(map(str, [x[2] for x in interval_set]))
            sys.stdout.write("%s\t%s\t%s\t%s\n" % (chrom, start, end, values))