Hamming distance match

hamming_match.py

'''
When you want to match files that only differ by one character: sample_1.fastq.gz sample_2.fastq.gz
'''
def is_hamming_match(seq1, seq2):
  if len(seq1) == len(seq2): # The sequences must be of equal length
    zipped = zip(seq1, seq2)
    is_mismatch = list([c1 != c2 for c1, c2 in zipped])
    if sum(is_mismatch) == 1: # The sequences must have only one mismatch
      mismatches = list([[c1, c2] for c1, c2 in zipped if c1 != c2])
      if mismatches[0] == ['1', '2'] or mismatches[0] == ['2', '1']:  # The mismatch must only be the numbers 1 and 2
        return True
  return False