From 11c000ab1166df0b0797548532de4c21c8ae75da Mon Sep 17 00:00:00 2001 From: takadonet Date: Tue, 27 Sep 2016 08:54:36 -0500 Subject: [PATCH] added 'N' => 'N' when doing reverse complement. Former-commit-id: 8ee756e3b47b1412e8030b437950c64827f3fd75 --- stringMLST.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/stringMLST.py b/stringMLST.py index f787329..0020c9c 100755 --- 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: