def local2global(local_coord, start, end, strand):
    """Return global coordinate within a region from a local coordinate"""

    # swap if strands disagree
    if strand == 1:
        return local_coord + start
    else:
        return end - local_coord