INNER CODE UNIT · Python
af
aquaskyline/Clairvoyante · clairvoyante/callVar.py:102
af = sum(XBatch[j,param.flankingBaseNum+1,:,1]) / dp
if varLength != maxVarLength:
for k in range(param.flankingBaseNum+1, param.flankingBaseNum+varLength+1):
altBase += num2base[np.argmax(XBatch[j,k,:,1])]
else:
for k in range(param.flankingBaseNum+1, 2*param.flankingBaseNum+1):
referenceTensor = XBatch[j,k,:,0]; insertionTensor = XBatch[j,k,:,1]
if k < (param.flankingBaseNum + maxVarLength) or sum(insertionTensor) >= (inferIndelLengthMinimumAF * sum(referenceTensor)):
inferredIndelLength += 1
altBase += num2base[np.argmax(insertionTensor)]
else:
break
coordination = int(coordination)
refBase = refSeq[param.flankingBaseNum]
# insertions longer than (param.flankingBaseNum-1) are marked SV
if inferredIndelLength >= param.flankingBaseNum:
altBase = "<INS>"
info.append("SVTYPE=INS")