Esempio n. 1
0
def updateAdjustedMismatches(fN, fF, guValue = .5, otherValue = 1.0):
    
    NX = Nexus(fN, fF)
    NX.load(['sigMask', 'adjustedNumMismatches'])

    while NX.nextID():

        mask = NX.sigMask
        numGU = mask.count('G')
        numGapAndMM = mask.count('X')
        
        NX.adjustedNumMismatches = (numGU * guValue) + (numGapAndMM * otherValue)

    NX.save()