changed name of wc function to rc

Author: dave-doty

examples/16_helix_origami_barrel_from_algoSST_paper.py

@@ -234,7 +234,7 @@ def add_angle_inducing_insertions_deletions(design) -> None:
 # above is in order from right to left on helix 1, followed by left to right on helix 18
 
 def assign_dna_to_unzipper_toeholds(design: sc.Design) -> None:
-    uz_toes = [sc.wc(seq) for seq in uz_toes_wc]
+    uz_toes = [sc.rc(seq) for seq in uz_toes_wc]
 
     strands_h1 = design.strands_starting_on_helix(1)
     strands_h1.sort(key=lambda _strand: _strand.first_domain().offset_5p())

scadnano/scadnano.py

@@ -2821,18 +2821,18 @@ def _add_key_value_to_json_map_if_not_default(
 # Default Extension object to allow for access to default extension values. num_bases is a dummy.
 _default_extension: Extension = Extension(num_bases=5)
 
-_wctable = str.maketrans('ACGTacgt', 'TGCAtgca')
+_rctable = str.maketrans('ACGTacgt', 'TGCAtgca')
 
 
-def wc(seq: str) -> str:
+def rc(seq: str) -> str:
     """
-    Return reverse Watson-Crick complement of `seq`.
+    Return reverse complement of `seq`.
     For example, ``wc('AACCTG')`` returns ``'CAGGTT'``.
 
     :param seq: a DNA sequence
-    :return: reverse Watson-Crick complement of `seq`.
+    :return: reverse complement of `seq`.
     """
-    return seq.translate(_wctable)[::-1]
+    return seq.translate(_rctable)[::-1]
 
 
 @dataclass
@@ -4279,7 +4279,7 @@ def assign_dna_complement_from(self, other: Strand) -> None:
                     other_seq = domain_other.dna_sequence_in(overlap_left, overlap_right - 1)
                     if other_seq is None:
                         raise ValueError(f'no DNA sequence has been assigned to strand {other}')
-                    overlap_complement = wc(other_seq)
+                    overlap_complement = rc(other_seq)
                     domain_complement_builder.append(wildcards)
                     domain_complement_builder.append(overlap_complement)
                     start_idx = overlap_right
@@ -7959,7 +7959,7 @@ def to_oxview_json(self, warn_duplicate_strand_names: bool = True,
                                     (sc_strand2.dna_sequence is not None) and
                                     (sc_strand1.dna_sequence[d1] != DNA_base_wildcard) and
                                     (sc_strand2.dna_sequence[d2] != DNA_base_wildcard) and
-                                    (wc(sc_strand1.dna_sequence[d1]) != sc_strand2.dna_sequence[d2])):
+                                    (rc(sc_strand1.dna_sequence[d1]) != sc_strand2.dna_sequence[d2])):
                                 continue
 
                             oxv_strand1['monomers'][d1]['bp'] = s2_nuc_idx + d2