diff --git a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/ReadBackedPhasingWalker.java b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/ReadBackedPhasingWalker.java index 734ade376..dc5dfc907 100755 --- a/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/ReadBackedPhasingWalker.java +++ b/public/java/src/org/broadinstitute/sting/gatk/walkers/phasing/ReadBackedPhasingWalker.java @@ -129,6 +129,9 @@ public class ReadBackedPhasingWalker extends RodWalker unphasedSiteQueue = null; @@ -487,6 +490,13 @@ public class ReadBackedPhasingWalker extends RodWalker readsAtHetSites = null; + private void clearFields() { + hetGenotypes = null; + prevHetAndInteriorIt = null; + phasingSiteIndex = -1; + readsAtHetSites = null; + } + public boolean hasPreviousHets() { return phasingSiteIndex > 0; } @@ -513,12 +523,20 @@ public class ReadBackedPhasingWalker extends RodWalker= 10; cacheWindow = 20000; has inconsistent sites", spec); } + @Test + public void test7() { + WalkerTestSpec spec = new WalkerTestSpec( + baseTestString(hg18Reference, "phasing_test_chr20_332341_1332503.bam", "CEU.trio.2010_03.genotypes.hg18.vcf", 20000, 10, 10) + + " -L chr20:332341-802503", + 1, + Arrays.asList("c37548b333b65f58d0edfc5c2a62a28a")); + executeTest("Use trio-phased VCF, but ignore its phasing [TEST SEVEN]", spec); + } + + @Test + public void test8() { + WalkerTestSpec spec = new WalkerTestSpec( + baseTestString(hg18Reference, "phasing_test_chr20_332341_1332503.bam", "CEU.trio.2010_03.genotypes.hg18.vcf", 20000, 10, 10) + + " -L chr20:332341-802503" + " -respectPhaseInInput", + 1, + Arrays.asList("dfc7cdddd702e63d46d04f61a3ecd720")); + executeTest("Use trio-phased VCF, and respect its phasing [TEST EIGHT]", spec); + } + }