def rightmost_from_read(read):
    """Function that takes as input a read a returns its rightmost mapping position"""
    rightmost = 0
    for cigar in read.cigar:
        if cigar[0] == 0:
            rightmost += cigar[1]
        elif cigar[0] == 2:
            rightmost += cigar[1]
        elif cigar[0] == 3:
            rightmost += cigar[1]
    return read.reference_start + rightmost
