From 0093482c628180a23b88a604b165fd5b0e79dd1b Mon Sep 17 00:00:00 2001
From: depristo <depristo@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Wed, 9 Sep 2009 21:19:36 +0000
Subject: [PATCH] N reference base fix for SSG

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1572 348d0f76-0448-11de-a6fe-93d51630548a
---
 .../genotyper/SingleSampleGenotyper.java      | 22 +++++++++++--------
 .../broadinstitute/sting/utils/BaseUtils.java |  4 ++++
 2 files changed, 17 insertions(+), 9 deletions(-)

diff --git a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SingleSampleGenotyper.java b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SingleSampleGenotyper.java
index be23c1b0e..c78c3ada5 100644
--- a/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SingleSampleGenotyper.java
+++ b/java/src/org/broadinstitute/sting/gatk/walkers/genotyper/SingleSampleGenotyper.java
@@ -92,18 +92,22 @@ public class SingleSampleGenotyper extends LocusWalker<SSGenotypeCall, SingleSam
      */
     public SSGenotypeCall map(RefMetaDataTracker tracker, ReferenceContext refContext, AlignmentContext context) {
         char ref = refContext.getBase();
-        DiploidGenotypePriors priors = new DiploidGenotypePriors(ref, heterozygosity, DiploidGenotypePriors.PROB_OF_TRISTATE_GENOTYPE);
+        if ( BaseUtils.isRegularBase(ref) ) {
+            DiploidGenotypePriors priors = new DiploidGenotypePriors(ref, heterozygosity, DiploidGenotypePriors.PROB_OF_TRISTATE_GENOTYPE);
 
-        // setup GenotypeLike object
-        GenotypeLikelihoods gl = GenotypeLikelihoodsFactory.makeGenotypeLikelihoods(baseModel, priors, defaultPlatform);
+            // setup GenotypeLike object
+            GenotypeLikelihoods gl = GenotypeLikelihoodsFactory.makeGenotypeLikelihoods(baseModel, priors, defaultPlatform);
 
-        gl.setVerbose(VERBOSE);
-        gl.setEnableCacheFlag(! disableCache);
-        ReadBackedPileup pileup = new ReadBackedPileup(ref, context);
-        gl.add(pileup, true);
-        gl.validate();
+            gl.setVerbose(VERBOSE);
+            gl.setEnableCacheFlag(! disableCache);
+            ReadBackedPileup pileup = new ReadBackedPileup(ref, context);
+            gl.add(pileup, true);
+            gl.validate();
 
-        return new SSGenotypeCall(!GENOTYPE, context.getLocation(), ref,gl, pileup);
+            return new SSGenotypeCall(!GENOTYPE, context.getLocation(), ref,gl, pileup);
+        } else {
+            return null;
+        }
     }
 
     /**
diff --git a/java/src/org/broadinstitute/sting/utils/BaseUtils.java b/java/src/org/broadinstitute/sting/utils/BaseUtils.java
index 68a65dd53..c974c95ce 100644
--- a/java/src/org/broadinstitute/sting/utils/BaseUtils.java
+++ b/java/src/org/broadinstitute/sting/utils/BaseUtils.java
@@ -133,6 +133,10 @@ public class BaseUtils {
         }
     }
 
+    static public boolean isRegularBase(char base) {
+        return simpleBaseToBaseIndex(base) != -1;
+    }
+
     /**
      * Converts a base index to a simple base
      *