diff --git a/java/src/org/broadinstitute/sting/gatk/contexts/variantcontext/VariantContextUtils.java b/java/src/org/broadinstitute/sting/gatk/contexts/variantcontext/VariantContextUtils.java index 547d4d344..83a16f680 100755 --- a/java/src/org/broadinstitute/sting/gatk/contexts/variantcontext/VariantContextUtils.java +++ b/java/src/org/broadinstitute/sting/gatk/contexts/variantcontext/VariantContextUtils.java @@ -301,6 +301,8 @@ public static VariantContext simpleMerge(Collection unsortedVCs, setValue = "Intersection"; else if ( nFiltered == VCs.size() ) // everything was filtered out setValue = "FilteredInAll"; + else if ( nVariant == 0 ) // everyone was reference + setValue = "ReferenceInAll"; else { // we are filtered in some subset List s = new ArrayList(); for ( VariantContext vc : VCs ) diff --git a/java/test/org/broadinstitute/sting/gatk/walkers/variantutils/CombineVariantsIntegrationTest.java b/java/test/org/broadinstitute/sting/gatk/walkers/variantutils/CombineVariantsIntegrationTest.java index 31fdf389d..13a483eca 100755 --- a/java/test/org/broadinstitute/sting/gatk/walkers/variantutils/CombineVariantsIntegrationTest.java +++ b/java/test/org/broadinstitute/sting/gatk/walkers/variantutils/CombineVariantsIntegrationTest.java @@ -74,4 +74,20 @@ public class CombineVariantsIntegrationTest extends WalkerTest { @Test public void combineSNPsAndIndels() { combine2("CEU.trio.2010_03.genotypes.vcf.gz", "CEU.dindel.vcf4.low_coverage.2010_06.indel.genotypes.vcf", "a66a799b1ae9fd09b40f78af6ef538d8"); } @Test public void uniqueSNPs() { combine2("pilot2.snps.vcf4.genotypes.vcf", "yri.trio.gatk_glftrio.intersection.annotated.filtered.chr1.vcf", "a9126d1cbe1fdf741236763fb3e3461f"); } + + @Test public void threeWayWithRefs() { + WalkerTestSpec spec = new WalkerTestSpec( + baseTestString(" -B NA19240_BGI,VCF4,"+validationDataLocation+"NA19240.BGI.RG.vcf" + + " -B NA19240_ILLUMINA,VCF4,"+validationDataLocation+"NA19240.ILLUMINA.RG.vcf" + + " -B NA19240_WUGSC,VCF4,"+validationDataLocation+"NA19240.WUGSC.RG.vcf" + + " -B denovoInfo,VCF,"+validationDataLocation+"yri_merged_validation_data_240610.annotated.b36.vcf" + + " -setKey centerSet" + + " -variantMergeOptions UNION" + + " -priority NA19240_BGI,NA19240_ILLUMINA,NA19240_WUGSC,denovoInfo" + + " -genotypeMergeOptions UNIQUIFY -L 1"), + 1, + Arrays.asList("37c66ca16fb5f890bcaa8e226a5f2e61")); + executeTest("threeWayWithRefs", spec); + + } }