def get_bases_at_site_fetch(alignment_obj, chromosome, position,
    reference_base, alternative_base):
    """

    So to convert INS from 0 to 1-based system:
    start = start + 1; end = end + 1
    and to convert INS from 1 to 0-based system:
    start = start - 1; end = end - 1

    :param alignment_obj:
    :param chromosome:
    :param position:
    :return:
    """
    position = int(position)
    base_list = []
    counts = alignment_obj.count_coverage(chromosome, position, position + 1)
    base_list = []
    print(counts[0][0])
    """
    # for read in alignment_obj.fetch(chromosome,  position, position + 1, multiple_iterators=True):
        print(read.__dict__)
        if not read.is_duplicate:
            for i in range(len(read.query_sequence)):
                tmp_global_pos = read.pos + i
    """
    return base_list
