3334353637383940414243
c2r.init(reference); int i=0; for (SAMRecord read : reader) { int numMismatches = c2r.numMismatches(read); int numIndelBases = getNumIndelBases(read); int editDistance = numMismatches + numIndelBases; read.setAttribute("NM", editDistance);