def get_pol_region(alignment):
    """
    Slices the alignment to get the pol region only. Uses HXB2 coordinatese for the selection.
    """
    HXB2_start = 2084
    HXB2_stop = 5095
    return alignment[:, HXB2_start:HXB2_stop + 1]