Opossum.py 文件源码

def ReadCigarLength(cigar, k=100) :

    if k < len(cigar) :
        cigarlist = cigar[:(k+1)]
    else :
        cigarlist = cigar

    l = sum([c[1] for c in cigarlist])

    return l


# Input cigar and pos denoting index within a cigar transformed into a string sequence.
# Return what is the index in the original cigar and how many extra positions from the beginning
# of this index is the current position.