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 1a083dbd8..5ca89856f 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java @@ -241,25 +241,27 @@ public class AlleleFrequencyWalker extends LocusWalker= 0 && baseIndex < 4); break; - } + int called2num = QualityUtils.compressedQualityToBaseIndex(hex_qual); - double qual2 = QualityUtils.compressedQualityToProb(hex_qual); - //System.out.printf("2ND %x %d %f\n", hex_qual, called2num, qual2); - quals[i][called2num] = qual2; + /* + int crossTalkPartnerIndex = BaseUtils.crossTalkPartnerIndex(callednum); + + if (called2num == crossTalkPartnerIndex) { + double nonref_quals = (1.0 - quals[i][callednum]) / 3.0; + for (int b=0; b<4; b++) + if (b != callednum) + quals[i][b] = nonref_quals; + } else { + */ + 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]) / 2; - for (int b=0; b<4; b++) - if (b != callednum && b != called2num) - quals[i][b] = nonref_quals; + double nonref_quals = (1.0 - quals[i][callednum] - quals[i][called2num]) / 2.0; + for (int b=0; b<4; b++) + if (b != callednum && b != called2num) + quals[i][b] = nonref_quals; + //} } } @@ -569,18 +571,24 @@ public class AlleleFrequencyWalker extends LocusWalker