def get_ref(genome, chrom, start, end):
    """
    Fetch genome reference given chrom, start, end
    :param genome: pysam fasta file
    :param chrom:
    :param start:
    :param end:
    :return: DNA sequence
    """
    return genome.fetch(chrom, int(start) - 1, int(end) - 1)