diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/varianteval/VariantEvalWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/varianteval/VariantEvalWalker.java index 78da993b4..b4f1fe89c 100755 --- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/varianteval/VariantEvalWalker.java +++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/varianteval/VariantEvalWalker.java @@ -155,8 +155,16 @@ public class VariantEvalWalker extends RefWalker { AllelicVariant eval = (AllelicVariant)tracker.lookup("eval", null); //if ( eval!=null ) System.out.printf("Eval: %f %d %b%n", eval.getVariationConfidence(), minDiscoveryQ, eval.getVariationConfidence() < minDiscoveryQ); - if ( eval != null && eval.getVariationConfidence() < minDiscoveryQ ) - eval = null; + if ( eval != null ) { + out.println(eval); + if ( evalContainsGenotypes ) { + // Genotyping - use best vs. next best lod + if ( eval.getConsensusConfidence() < minDiscoveryQ ) eval = null; + } else { + // Variant discovery - use best vs. reference lod + if ( Math.abs(eval.getVariationConfidence()) < minDiscoveryQ ) eval = null; + } + } // update stats about all of the SNPs updateAnalysisSet(ALL_SNPS, eval, tracker, ref, context);