From eb2f0ebd62a5b8764a4fd8b746e5345d4af1783a Mon Sep 17 00:00:00 2001
From: kiran <kiran@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Wed, 8 Apr 2009 19:58:18 +0000
Subject: [PATCH] If the first base of a read is 'N', and the alignment cigar
 says every base matches, samtools calls shennanigans.  Now I just output an
 A, but the real way to do this is to modify the cigar string accordingly.

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@337 348d0f76-0448-11de-a6fe-93d51630548a
---
 .../playground/fourbasecaller/RepairBadlyCombinedSamFile.java   | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/java/src/org/broadinstitute/sting/playground/fourbasecaller/RepairBadlyCombinedSamFile.java b/java/src/org/broadinstitute/sting/playground/fourbasecaller/RepairBadlyCombinedSamFile.java
index 474d2e2f3..2b09b718a 100644
--- a/java/src/org/broadinstitute/sting/playground/fourbasecaller/RepairBadlyCombinedSamFile.java
+++ b/java/src/org/broadinstitute/sting/playground/fourbasecaller/RepairBadlyCombinedSamFile.java
@@ -50,7 +50,7 @@ public class RepairBadlyCombinedSamFile extends CommandLineProgram {
                         case 'C': newbase = 'G'; break;
                         case 'G': newbase = 'C'; break;
                         case 'T': newbase = 'A'; break;
-                        default: newbase = '.';
+                        default: newbase = 'A'; // hack.  This is so downstream applications don't barf on an ambiguous base 'matching' the reference
                     }
                     newbases[i] = newbase;
                     newquals[i] = quals[bases.length - i - 1];