Merge pull request #13 from Takadonet/fix_n_seq

added 'N' => 'N' when doing reverse complement.
jordanlab · Oct 3, 2016 · 2515ebf · 2515ebf
2 parents c6bfb7d + 8ee756e
commit 2515ebf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/stringMLST.py b/stringMLST.py
@@ -564,7 +564,7 @@ def printResults(results,output_filename,overwrite,timeDisp):
 """Returns the reverse complement of the sequence"""		
 def reverseComplement(seq):
 	seqU  = seq.upper()
-	seq_dict = {'A':'T','T':'A','G':'C','C':'G','Y':'R','R':'Y','S':'S','W':'W','K':'M','M':'K'}
+	seq_dict = {'A':'T','T':'A','G':'C','C':'G','Y':'R','R':'Y','S':'S','W':'W','K':'M','M':'K','N':'N'}
 	try:
 		return "".join([seq_dict[base] for base in reversed(seqU)])
 	except: