diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantRecalibrator.java b/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantRecalibrator.java index 3fe3159f1..3a7dfc423 100755 --- a/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantRecalibrator.java +++ b/java/src/org/broadinstitute/sting/gatk/walkers/variantrecalibration/VariantRecalibrator.java @@ -245,7 +245,7 @@ public class VariantRecalibrator extends RodWalker vcsTruth = tracker.getVariantContexts(ref, truthNames, null, context.getLocation(), false, true); + final Collection vcsTruth = tracker.getVariantContexts(ref, truthNames, null, context.getLocation(), true, true); boolean isAtTruthSite = false; for( final VariantContext vcTruth : vcsTruth ) { if( vcTruth != null && vcTruth.isVariant() && !vcTruth.isFiltered() && (!vcTruth.hasGenotypes() || vcTruth.isPolymorphic()) ) { @@ -255,12 +255,14 @@ public class VariantRecalibrator extends RodWalker= QUAL_THRESHOLD ) { final VariantDatum variantDatum = new VariantDatum(); - variantDatum.isTransition = VariantContextUtils.getSNPSubstitutionType(vc).compareTo(BaseUtils.BaseSubstitutionType.TRANSITION) == 0; + + if (vc.isSNP()) + variantDatum.isTransition = VariantContextUtils.getSNPSubstitutionType(vc).compareTo(BaseUtils.BaseSubstitutionType.TRANSITION) == 0; final Collection vcsDbsnp = tracker.getVariantContexts(ref, "dbsnp", null, context.getLocation(), false, true); final Collection vcsHapMap = tracker.getVariantContexts(ref, "hapmap", null, context.getLocation(), false, true);