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 10cffa09c..c0504b6a4 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/AlleleFrequencyWalker.java @@ -356,6 +356,7 @@ public class AlleleFrequencyWalker extends LocusWalker -5.0)) + String[] tokens; + tokens = alleleFreq.location.split(":"); + int current_offset = Integer.parseInt(tokens[1]); + + if (inside_confident_ref_interval && + ((alleleFreq.lodVsRef > -5.0) || (current_offset != last_position_considered + 1))) { // No longer hom-ref, so output a ref line. - String[] tokens; tokens = confident_ref_interval_start.split(":"); String contig = tokens[0]; int start = Integer.parseInt(tokens[1]); tokens = alleleFreq.location.split(":"); - int end = Integer.parseInt(tokens[1])-1; + int end = last_position_considered; double lod = confident_ref_interval_LOD_sum / confident_ref_interval_length; @@ -412,6 +418,8 @@ public class AlleleFrequencyWalker extends LocusWalker= 5) { this.output.print(alleleFreq.asGFFString()); } return "";