From 37b8d703212f4c1bbd24a5e99531449682e8060a Mon Sep 17 00:00:00 2001 From: Mark DePristo Date: Mon, 4 Jun 2012 07:04:14 -0400 Subject: [PATCH] Hidden option to SelectVariants to force the genotypes information to be decoded by computing AC --- .../gatk/walkers/variantutils/SelectVariants.java | 14 +++++++++++++- 1 file changed, 13 insertions(+), 1 deletion(-) diff --git a/public/java/src/org/broadinstitute/sting/gatk/walkers/variantutils/SelectVariants.java b/public/java/src/org/broadinstitute/sting/gatk/walkers/variantutils/SelectVariants.java index 3805e48df..690cc71ac 100755 --- a/public/java/src/org/broadinstitute/sting/gatk/walkers/variantutils/SelectVariants.java +++ b/public/java/src/org/broadinstitute/sting/gatk/walkers/variantutils/SelectVariants.java @@ -315,6 +315,11 @@ public class SelectVariants extends RodWalker implements TreeR @Argument(fullName="fullyDecode", doc="If true, the incoming VariantContext will be fully decoded", required=false) private boolean fullyDecode = false; + @Hidden + @Argument(fullName="forceGenotypesDecode", doc="If true, the incoming VariantContext will have its genotypes forcibly decoded by computing AC across all genotypes. For efficiency testing only", required=false) + private boolean forceGenotypesDecode = false; + + /* Private class used to store the intermediate variants in the integer random selection process */ private class RandomVariantStructure { private VariantContext vc; @@ -495,7 +500,14 @@ public class SelectVariants extends RodWalker implements TreeR } for (VariantContext vc : vcs) { - if ( fullyDecode ) vc = vc.fullyDecode(vcfRods.get(vc.getSource())); + // an option for performance testing only + if ( fullyDecode ) + vc = vc.fullyDecode(vcfRods.get(vc.getSource())); + + // an option for performance testing only + if ( forceGenotypesDecode ) + logger.info("forceGenotypesDecode with getCalledChrCount() = " + vc.getCalledChrCount()); + if ( IDsToKeep != null && ! IDsToKeep.contains(vc.getID()) ) continue;