Skip to content

Commit

Permalink
Merge pull request jordanlab#13 from Takadonet/fix_n_seq
Browse files Browse the repository at this point in the history
added 'N' => 'N' when doing reverse complement.

Former-commit-id: 2515ebf
  • Loading branch information
anujg30 authored Oct 3, 2016
2 parents cc0fce2 + 11c000a commit 1f2be4d
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion stringMLST.py
Original file line number Diff line number Diff line change
Expand Up @@ -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:
Expand Down

0 comments on commit 1f2be4d

Please sign in to comment.