def getBaseAt(read, pos):
    """
    Extract nucleotide within a read at a given set of positions

    read = read object [class pysam.AlignedSegment]
    pos = positions to convert [list]

    """
    nuc = []
    for p in pos:
        # print (p)
        nuc.append(read.seq[p])
    return nuc