diff --git a/protected/gatk-protected/src/main/java/org/broadinstitute/sting/gatk/walkers/variantutils/GenotypeGVCFs.java b/protected/gatk-protected/src/main/java/org/broadinstitute/sting/gatk/walkers/variantutils/GenotypeGVCFs.java index 4620481af..00def65c0 100644 --- a/protected/gatk-protected/src/main/java/org/broadinstitute/sting/gatk/walkers/variantutils/GenotypeGVCFs.java +++ b/protected/gatk-protected/src/main/java/org/broadinstitute/sting/gatk/walkers/variantutils/GenotypeGVCFs.java @@ -207,11 +207,15 @@ public class GenotypeGVCFs extends RodWalker originalAttributes = combinedVC.getAttributes(); // only re-genotype polymorphic sites - if ( combinedVC.isVariant() ) - result = new VariantContextBuilder(genotypingEngine.calculateGenotypes(result)).attributes(originalAttributes).make(); + if ( combinedVC.isVariant() ) { + final VariantContext regenotypedVC = genotypingEngine.calculateGenotypes(result); + if ( regenotypedVC == null ) + return null; + result = new VariantContextBuilder(regenotypedVC).attributes(originalAttributes).make(); + } // if it turned monomorphic and we don't want such sites, quit - if ( result == null || (!INCLUDE_NON_VARIANTS && result.isMonomorphicInSamples()) ) + if ( !INCLUDE_NON_VARIANTS && result.isMonomorphicInSamples() ) return null; // re-annotate it