updated reverse complement to handle ambiguous bases

shendurelab · Feb 23, 2014 · ec0acf8 · ec0acf8
1 parent 71d8129
commit ec0acf8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/MinusSVMipv4.cpp b/MinusSVMipv4.cpp
@@ -21,6 +21,8 @@ string reverse_comp (string seq) {
                         case 'T':
                                 rev_comp.append("A");
                                 break;
+                        default:
+                                rev_comp.append(string(1, seq[i]));
           	 	 }
         	}
        	return rev_comp;

diff --git a/mipgen b/mipgen