From b83df5616afe020686cd9229dc0b6fb8daa17f9b Mon Sep 17 00:00:00 2001
From: chartl <chartl@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Tue, 20 Oct 2009 17:36:31 +0000
Subject: [PATCH] Change for lower-case references (always compare upper case
 bases)

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1888 348d0f76-0448-11de-a6fe-93d51630548a
---
 .../walkers/BaseTransitionTableCalculatorJavaWalker.java     | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/java/src/org/broadinstitute/sting/playground/gatk/walkers/BaseTransitionTableCalculatorJavaWalker.java b/java/src/org/broadinstitute/sting/playground/gatk/walkers/BaseTransitionTableCalculatorJavaWalker.java
index fb3080fee..cd11176ae 100644
--- a/java/src/org/broadinstitute/sting/playground/gatk/walkers/BaseTransitionTableCalculatorJavaWalker.java
+++ b/java/src/org/broadinstitute/sting/playground/gatk/walkers/BaseTransitionTableCalculatorJavaWalker.java
@@ -138,7 +138,7 @@ public class BaseTransitionTableCalculatorJavaWalker extends LocusWalker<Referen
 
     public boolean useRead( SAMRecord read, int offset, ReferenceContext ref ) {
 
-        if ( Character.toUpperCase(read.getReadBases()[offset]) == ref.getBase() ) {
+        if ( Character.toUpperCase(read.getReadBases()[offset]) == Character.toUpperCase(ref.getBase()) ) {
             return false;
         } else if ( read.getMappingQuality() <= minMappingQuality ) {
             return false;
@@ -300,6 +300,9 @@ class BaseTransitionTable {
     }
 
     public void update(char observedBase, char refBase ) {
+        //if ( observedBase == refBase ) {
+        //    throw new StingException("BaseTransitionTable received equal observed and reference bases, which should not happen.");
+        //}
         table[BaseUtils.simpleBaseToBaseIndex(observedBase)][BaseUtils.simpleBaseToBaseIndex(refBase)] ++;
     }