code: def print_first_difference(arec, brec, ignore_case=False, ignore_N=False,
          rc=False, report_match=True):
      
      plus_match = _print_first_difference(arec, brec, ignore_case=ignore_case,
              ignore_N=ignore_N, report_match=report_match)
      if rc and not plus_match:
          logging.debug("trying reverse complement of %s" % brec.id)
          brec.seq = brec.seq.reverse_complement()
          minus_match = _print_first_difference(arec, brec,
                  ignore_case=ignore_case, ignore_N=ignore_N,
                  report_match=report_match)
          return minus_match
  
      else:
          return plus_match
docs: Returns the first different nucleotide in two sequence comparisons
      runs both Plus and Minus strand