diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java index 003b944cd..91ddf973c 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java @@ -150,11 +150,12 @@ public class AlleleFrequencyWalker extends LocusWalker 0); - Object SQ_field = reads.get(0).getAttribute("SQ"); + Object SQ_field = reads.get(i).getAttribute("SQ"); if (SQ_field == null || FORCE_1BASE_PROBS) { // Set all nonref qual scores to their share of the remaining probality not "used" by the reference base's qual double nonref_quals = (1.0 - quals[i][callednum]) / 3; @@ -168,13 +169,22 @@ public class AlleleFrequencyWalker extends LocusWalker= 0 && baseIndex < 4); break; + } + double qual2 = QualityUtils.compressedQualityToProb(hex_qual); //System.out.printf("2ND %x %d %f\n", hex_qual, called2num, qual2); quals[i][called2num] = qual2; - // - double nonref_quals = (1.0 - quals[i][callednum] - quals[i][called2num]) / 3; + double nonref_quals = (1.0 - quals[i][callednum] - quals[i][called2num]) / 2; for (int b=0; b<4; b++) if (b != callednum && b != called2num) quals[i][b] = nonref_quals;