INNER CODE UNIT · Python

af

aquaskyline/Clairvoyante · clairvoyante/callVar.py:124

                    af = sum(XBatch[j,param.flankingBaseNum+1,:,2]) / dp
                    # infer the deletion length
                    if varLength == maxVarLength:
                        for k in range(param.flankingBaseNum+1, 2*param.flankingBaseNum+1):
                            if k < (param.flankingBaseNum + maxVarLength) or sum(XBatch[j,k,:,2]) >= (inferIndelLengthMinimumAF * sum(XBatch[j,k,:,0])):
                                inferredIndelLength += 1
                            else:
                               break
                    # deletions longer than (param.flankingBaseNum-1) are marked SV
                    coordination = int(coordination)
                    if inferredIndelLength >= param.flankingBaseNum:
                        refBase = refSeq[param.flankingBaseNum]
                        altBase = "<DEL>"
                        info.append("SVTYPE=DEL")
                    elif varLength != maxVarLength:
                        refBase = refSeq[param.flankingBaseNum:param.flankingBaseNum+varLength+1]
                        altBase = refSeq[param.flankingBaseNum]
                    else:

View source record →