pydna 5.5.5__tar.gz → 5.5.7__tar.gz

{pydna-5.5.5 → pydna-5.5.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pydna
-Version: 5.5.5
+Version: 5.5.7
 Summary: Representing double stranded DNA and functions for simulating cloning and homologous recombination between DNA molecules.
 License: BSD
 License-File: LICENSE.txt
@@ -39,7 +39,7 @@ Requires-Dist: pydivsufsort (>=0.0.14)
 Requires-Dist: pyfiglet (==0.8.post1)
 Requires-Dist: pyparsing (>=2.4.7) ; extra == "download"
 Requires-Dist: pyperclip (>=1.8.2) ; extra == "clipboard"
-Requires-Dist: regex (>=2024.11.6,<2025.0.0)
+Requires-Dist: regex (>=2024.11.6,<2027.0.0)
 Requires-Dist: requests (>=2.26.0) ; extra == "download"
 Requires-Dist: scipy (>=1.11.3) ; (python_version >= "3.12") and (extra == "gel")
 Requires-Dist: scipy (>=1.9.3) ; (python_version < "3.12") and (extra == "gel")

{pydna-5.5.5 → pydna-5.5.7}/pyproject.toml

@@ -35,7 +35,7 @@ license = "BSD"
 name = "pydna"
 readme = "README.md"
 repository = "https://github.com/pydna-group/pydna/tree/master"
-version = "5.5.5"
+version = "5.5.7"
 [tool.poetry.dependencies]
 appdirs = ">=1.4.4"
 biopython = "1.85"
@@ -60,7 +60,7 @@ scipy = [
   { version = ">=1.9.3", python = "<3.12", optional = true },
 ]
 seguid = ">=0.0.5"
-regex = "^2024.11.6"
+regex = ">=2024.11.6,<2027.0.0"
 opencloning-linkml = "^0.4.9"
 [tool.poetry.extras]
 clipboard = ["pyperclip"]

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/__init__.py

@@ -49,23 +49,27 @@ functions with a lowercase letter:
          ├── amplify
          │         ├── Anneal
          │         └── pcr
+         │
          ├── assembly
          │          └── Assembly
+         │
          ├── design
          │        ├── assembly_fragments
          │        └── primer_design
-         ├── download
-         │          └── download_text
+         │
          ├── dseqrecord
          │            └── Dseqrecord
          ├── gel
          │     └── Gel
+         │
          ├── genbank
          │         ├── genbank
          │         └── Genbank
+         │
          ├── parsers
          │         ├── parse
          │         └── parse_primers
+         │
          └── readers
                    ├── read
                    └── read_primers
@@ -143,7 +147,7 @@ __license__ = "BSD"
 __maintainer__ = "Björn Johansson"
 __email__ = "bjorn_johansson@bio.uminho.pt"
 __status__ = "Development"  # "Production" #"Prototype"
-__version__ = "5.5.5"
+__version__ = "5.5.7"
 
 
 class _PydnaWarning(Warning):

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/all.py

@@ -18,7 +18,7 @@ ttt
 Dseqrecord(-3)
 >>> from pydna.all import __all__
 >>> __all__
-['Anneal', 'pcr', 'Assembly', 'genbank', 'Genbank', 'download_text', 'Dseqrecord',
+['Anneal', 'pcr', 'Assembly', 'genbank', 'Genbank', 'Dseqrecord',
 'Dseq', 'read', 'read_primer', 'parse', 'parse_primers', 'primer_design', 'assembly_fragments', 'eq', 'gbtext_clean']
 >>>
 """
@@ -30,20 +30,16 @@ __all__ = [
     "Assembly",
     "genbank",
     "Genbank",
-    "download_text",
     "Dseqrecord",
     "Dseq",
     "read",
     "read_primer",
     "parse",
     "parse_primers",
-    # "ape",
     "primer_design",
     "assembly_fragments",
-    # "circular_assembly_fragments",
     "eq",
     "gbtext_clean",
-    # "PrimerList",
 ]
 
 
@@ -52,20 +48,13 @@ from pydna.amplify import pcr
 from pydna.assembly import Assembly
 from pydna.genbank import genbank
 from pydna.genbank import Genbank
-from pydna.download import download_text
 from pydna.dseqrecord import Dseqrecord
 from pydna.dseq import Dseq
 from pydna.readers import read
 from pydna.readers import read_primer
 from pydna.parsers import parse
 from pydna.parsers import parse_primers
-
-# from pydna.editor import ape
 from pydna.design import primer_design
 from pydna.design import assembly_fragments
-
-# from pydna.design import circular_assembly_fragments
 from pydna.utils import eq
 from pydna.genbankfixer import gbtext_clean
-
-# from pydna.myprimers import PrimerList

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/assembly2.py

@@ -357,7 +357,18 @@ def common_sub_strings(
     return [r for r in results if r not in shifted_matches]
 
 
-def gibson_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
+def _get_trim_end_info(
+    end_info: tuple[str, str], trim_ends: str, is_five_prime: bool
+) -> int | None:
+    """Utility function to get the trim information for terminal_overlap."""
+    if end_info[0] == trim_ends:
+        return len(end_info[1]) if is_five_prime else len(end_info[1]) * -1
+    return 0 if is_five_prime else None
+
+
+def terminal_overlap(
+    seqx: Dseqrecord, seqy: Dseqrecord, limit=25, trim_ends: None | str = None
+):
     """
     Assembly algorithm to find terminal overlaps (e.g. for Gibson assembly).
     The order matters, we want alignments like:
@@ -382,6 +393,9 @@ def gibson_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
         The second sequence
     limit : int
         Minimum length of the overlap
+    trim_ends : str
+        The ends to trim, either '5' or '3'
+        If None, no trimming is done
 
     Returns
     -------
@@ -389,32 +403,64 @@ def gibson_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
         A list of overlaps between the two sequences
 
     >>> from pydna.dseqrecord import Dseqrecord
-    >>> from pydna.assembly2 import gibson_overlap
+    >>> from pydna.assembly2 import terminal_overlap
     >>> x = Dseqrecord("ttactaAAAAAA")
     >>> y = Dseqrecord("AAAAAAcgcacg")
-    >>> gibson_overlap(x, y, limit=5)
+    >>> terminal_overlap(x, y, limit=5)
     [(6, 0, 6), (7, 0, 5)]
-    >>> gibson_overlap(y, x, limit=5)
+    >>> terminal_overlap(y, x, limit=5)
+    []
+
+    Trimming the ends:
+    >>> from pydna.dseq import Dseq
+    >>> from pydna.dseqrecord import Dseqrecord
+    >>> from pydna.assembly2 import terminal_overlap
+    >>> x = Dseqrecord(Dseq.from_full_sequence_and_overhangs("aaaACGT", 0, 3))
+    >>> y = Dseqrecord(Dseq.from_full_sequence_and_overhangs("ACGTccc", 3, 0))
+    >>> terminal_overlap(x, y, limit=4)
+    [(3, 0, 4)]
+    >>> terminal_overlap(x, y, limit=4, trim_ends="5'")
+    [(3, 0, 4)]
+    >>> terminal_overlap(x, y, limit=4, trim_ends="3'")
     []
     """
 
-    # Because Gibson enzymes remove 5' overhangs, we remove them from the sequence
-    # when looking for homology, then we shift the location of the second fragment accordingly.
-    # This is only relevant for linear fragments, so we don't need to worry about
-    # shifting locations for circular fragments.
-    trim_x_left = -seqx.seq.ovhg if seqx.seq.ovhg < 0 else 0
-    trim_x_right = seqx.seq.watson_ovhg if seqx.seq.watson_ovhg < 0 else None
-    trim_y_left = -seqy.seq.ovhg if seqy.seq.ovhg < 0 else 0
-    trim_y_right = seqy.seq.watson_ovhg if seqy.seq.watson_ovhg < 0 else None
-
-    stringx = str(seqx.seq[trim_x_left:trim_x_right]).upper()
-    stringy = str(seqy.seq[trim_y_left:trim_y_right]).upper()
+    if trim_ends is not None and trim_ends not in ["5'", "3'"]:
+        raise ValueError("trim_ends must be '5' or '3'")
+
+    if trim_ends is None:
+        trim_x_left, trim_x_right, trim_y_left, trim_y_right = (0, None, 0, None)
+        stringx = str(seqx.seq).upper()
+        stringy = str(seqy.seq).upper()
+    else:
+        trim_x_right = _get_trim_end_info(
+            seqx.seq.three_prime_end(), trim_ends, is_five_prime=False
+        )
+        trim_y_left = _get_trim_end_info(
+            seqy.seq.five_prime_end(), trim_ends, is_five_prime=True
+        )
+
+        # I actually don't think these two are needed, since only the terminal
+        # join between x_right and y_left is tested, but maybe there is some edge-case
+        # that I am missing, so keeping them just in case.
+        trim_x_left = _get_trim_end_info(
+            seqx.seq.five_prime_end(), trim_ends, is_five_prime=True
+        )
+        trim_y_right = _get_trim_end_info(
+            seqy.seq.three_prime_end(), trim_ends, is_five_prime=False
+        )
+
+        stringx = str(seqx.seq[trim_x_left:trim_x_right]).upper()
+        stringy = str(seqy.seq[trim_y_left:trim_y_right]).upper()
+
     # We have to convert to list because we need to modify the matches
     matches = [
         list(m)
         for m in common_sub_strings_str(stringx, stringy, limit)
         if (m[1] == 0 and m[0] + m[2] == len(stringx))
     ]
+
+    # Shift the matches if the left end has been trimmed
     for match in matches:
         match[0] += trim_x_left
         match[1] += trim_y_left
@@ -423,6 +469,31 @@ def gibson_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
     return [tuple(m) for m in matches]
 
 
+def gibson_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
+    """
+    Assembly algorithm to find terminal overlaps for Gibson assembly.
+    It is a wrapper around terminal_overlap with trim_ends="5'".
+    """
+
+    return terminal_overlap(seqx, seqy, limit, trim_ends="5'")
+
+
+def in_fusion_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
+    """
+    Assembly algorithm to find terminal overlaps for in-fusion assembly.
+    It is a wrapper around terminal_overlap with trim_ends="3'".
+    """
+    return terminal_overlap(seqx, seqy, limit, trim_ends="3'")
+
+
+def pcr_fusion_overlap(seqx: Dseqrecord, seqy: Dseqrecord, limit=25):
+    """
+    Assembly algorithm to find terminal overlaps for PCR fusion assembly.
+    It is a wrapper around terminal_overlap with trim_ends=None.
+    """
+    return terminal_overlap(seqx, seqy, limit, trim_ends=None)
+
+
 def sticky_end_sub_strings(seqx: Dseqrecord, seqy: Dseqrecord, limit: bool = False):
     """
     Assembly algorithm for ligation of sticky ends.
@@ -796,7 +867,7 @@ def assemble(
         f_u = fragments[u - 1] if u > 0 else fragments[-u - 1].reverse_complement()
         f_v = fragments[v - 1] if v > 0 else fragments[-v - 1].reverse_complement()
         seq_u = str(loc_u.extract(f_u).seq)
-        seq_v = str(loc_v.extract(f_v).seq.rc())
+        seq_v = str(loc_v.extract(f_v).seq.reverse_complement())
         # Test if seq_u and seq_v anneal
         if not anneal_strands(seq_u, seq_v):
             raise ValueError("Mismatch in assembly")
@@ -1578,8 +1649,9 @@ class Assembly:
             fragment2[f2_1_start:f2_2_end]
         )
 
-        if overlap_diff == 0:
-            assert False, "Overlap is 0"
+        # Safeguard
+        if overlap_diff == 0:  # pragma: no cover
+            raise AssertionError("Overlap is 0")
 
         if overlap_diff > 0:
             new_loc_f1_1 = create_location(
@@ -1874,7 +1946,7 @@ class PCRAssembly(Assembly):
         results = super().assemble_linear(only_adjacent_edges, max_assemblies)
         for result in results:
             rp = self.fragments[2]
-            result.seq = result.seq[: -len(rp)] + Dseq(str(rp.seq.rc()))
+            result.seq = result.seq[: -len(rp)] + Dseq(str(rp.seq.reverse_complement()))
         return results
 
 
@@ -2077,7 +2149,9 @@ def in_fusion_assembly(
         List of assembled DNA molecules
     """
 
-    products = gibson_assembly(frags, limit)
+    products = common_function_assembly_products(
+        frags, limit, in_fusion_overlap, circular_only
+    )
     return _recast_sources(products, InFusionSource)
 
 
@@ -2101,7 +2175,9 @@ def fusion_pcr_assembly(
     list[Dseqrecord]
         List of assembled DNA molecules
     """
-    products = gibson_assembly(frags, limit)
+    products = common_function_assembly_products(
+        frags, limit, pcr_fusion_overlap, circular_only
+    )
     return _recast_sources(products, OverlapExtensionPCRLigationSource)
 
 

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/dseq.py

@@ -785,7 +785,7 @@ class Dseq(Seq):
         >>> ssobj
         Dseq(-7)
         GATTACA
-        <BLANKLINE>
+        |||||||
         >>> round(ssobj.mw(), 1)
         2184.4
         >>> ds_lin_obj2 = Dseq("GATZFCA")
@@ -908,7 +908,8 @@ class Dseq(Seq):
         w, c = representation_tuple(
             self._data.decode("ascii"), length_limit_for_repr=length_limit_for_repr
         )
-
+        w = w or "|" * len(c)
+        c = c or "|" * len(w)
         return pretty_str(header + "\n" + w + "\n" + c)
 
     def reverse_complement(self) -> "Dseq":
@@ -1574,7 +1575,7 @@ class Dseq(Seq):
         ctag
         >>> ds.nibble_five_prime_left(4)
         Dseq(-4)
-        <BLANKLINE>
+        ||||
         ctag
         >>> ds = Dseq.from_representation(
         ... '''
@@ -1601,9 +1602,8 @@ class Dseq(Seq):
             DESCRIPTION.
 
         """
-        recessed = copy.deepcopy(self)
         n += max(0, self.ovhg or 0)
-        recessed = Dseq(
+        return Dseq(
             self._data[:n]
             .translate(dscode_to_crick_table)
             .translate(complement_table_for_dscode)
@@ -1611,7 +1611,6 @@ class Dseq(Seq):
             .lstrip()
             + self._data[n:]
         )
-        return recessed
 
     def nibble_five_prime_right(self: DseqType, n: int = 1) -> DseqType:
         """
@@ -1657,7 +1656,7 @@ class Dseq(Seq):
         >>> ds.nibble_five_prime_right(4)
         Dseq(-4)
         gatc
-        <BLANKLINE>
+        ||||
         >>> ds = Dseq.from_representation(
         ... '''
         ... gatc
@@ -1668,18 +1667,16 @@ class Dseq(Seq):
         gatc
         ctag
         """
-        recessed = copy.deepcopy(self)
         n = len(self) - n
         ovhg = len(self) if self.right_ovhg is None else self.right_ovhg
         n -= max(0, ovhg)
-        recessed = Dseq(
+        return Dseq(
             self._data[:n]
             + self._data[n:]
             .translate(dscode_to_watson_table)
             .translate(dscode_to_watson_tail_table)
             .lstrip()
         )
-        return recessed
 
     exo1_front = nibble_five_prime_left  # TODO: consider using the new names
     exo1_end = nibble_five_prime_right  # TODO: consider using the new names
@@ -1728,7 +1725,7 @@ class Dseq(Seq):
         >>> ds.nibble_three_prime_left(4)
         Dseq(-4)
         gatc
-        <BLANKLINE>
+        ||||
         >>> ds = Dseq.from_representation(
         ... '''
         ...   gatc
@@ -1745,14 +1742,13 @@ class Dseq(Seq):
         """
         ovhg = len(self) if self.ovhg is None else self.ovhg
         n -= min(0, ovhg)
-        recessed = Dseq(
+        return Dseq(
             self._data[:n]
             .translate(dscode_to_watson_table)
             .translate(dscode_to_watson_tail_table)
             .lstrip()
             + self._data[n:]
         )
-        return recessed
 
     def nibble_three_prime_right(self: DseqType, n=1) -> DseqType:
         """
@@ -1797,7 +1793,7 @@ class Dseq(Seq):
         ctag
         >>> ds.nibble_three_prime_right(4)
         Dseq(-4)
-        <BLANKLINE>
+        ||||
         ctag
         >>> ds = Dseq.from_representation(
         ... '''
@@ -1812,7 +1808,7 @@ class Dseq(Seq):
         n = len(self) - n
         ovhg = len(self) if self.right_ovhg is None else self.right_ovhg
         n += min(0, ovhg)
-        recessed = Dseq(
+        return Dseq(
             self._data[:n]
             + self._data[n:]
             .translate(dscode_to_crick_table)
@@ -1820,7 +1816,6 @@ class Dseq(Seq):
             .translate(dscode_to_crick_tail_table)
             .lstrip()
         )
-        return recessed
 
     def no_cutters(
         self, batch: Union[RestrictionBatch, None] = None
@@ -1878,7 +1873,7 @@ class Dseq(Seq):
             """docstring."""
             w = f"{self.ovhg * '-'}{self.watson}{'-' * (-self.ovhg + len(self.crick) - len(self.watson))}".upper()
             c = f"{'-' * (self.ovhg + len(self.watson) - len(self.crick))}{self.crick}{-self.ovhg * '-'}".upper()
-            cs = ldseguid(w, c, alphabet="{DNA-extended}")
+            cs = ldseguid(w, c, alphabet="{DNA-extended},AU")
         return cs
 
     def isblunt(self) -> bool:
@@ -2250,7 +2245,7 @@ class Dseq(Seq):
             # argument is probably a RestrictionBatch
             enzymes = [e for e in enzymes[0]]
 
-        enzymes = flatten(enzymes)
+        enzymes = list(dict.fromkeys(flatten(enzymes)))  # remove duplicate enzymes
         out = list()
         for e in enzymes:
             # Positions of the cut on the watson strand. They are 1-based, so we subtract
@@ -2589,7 +2584,7 @@ class Dseq(Seq):
         >>> strands[0]
         Dseq(-2)
         ta
-        <BLANKLINE>
+        ||
         >>> ds = Dseq("tagaaptapgtatg")
         >>> ds
         Dseq(-14)
@@ -2602,7 +2597,7 @@ class Dseq(Seq):
         atctt    catac
         >>> strands[0]
         Dseq(-2)
-        <BLANKLINE>
+        ||
         at
         """
 
@@ -2636,7 +2631,7 @@ class Dseq(Seq):
         >>> strands[0]
         Dseq(-2)
         ta
-        <BLANKLINE>
+        ||
         >>> new, strands = ds.shed_ss_dna([],[(6, 8)])
         >>> new
         Dseq(-14)
@@ -2644,7 +2639,7 @@ class Dseq(Seq):
         atcttc  ccatac
         >>> strands[0]
         Dseq(-2)
-        <BLANKLINE>
+        ||
         at
         >>> ds = Dseq("tagaagtaggtatg")
         >>> new, (strand1, strand2) = ds.shed_ss_dna([(6, 8), (9, 11)],[])
@@ -2655,11 +2650,11 @@ class Dseq(Seq):
         >>> strand1
         Dseq(-2)
         ta
-        <BLANKLINE>
+        ||
         >>> strand2
         Dseq(-2)
         gt
-        <BLANKLINE>
+        ||
         """
 
         watson_cutpairs = watson_cutpairs or list()
@@ -2878,14 +2873,14 @@ class Dseq(Seq):
         >>> Dseq(parts.sticky_left5)
         Dseq(-3)
         GGG
-        <BLANKLINE>
+        |||
         >>> Dseq(parts.middle)
         Dseq(-3)
         ATC
         TAG
         >>> Dseq(parts.sticky_right5)
         Dseq(-3)
-        <BLANKLINE>
+        |||
         TCA
 
         Parameters

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/dseqrecord.py

@@ -26,7 +26,7 @@ from Bio.SeqFeature import CompoundLocation
 from Bio.SeqFeature import SimpleLocation
 from pydna.seqrecord import SeqRecord
 from Bio.Seq import translate
-from pydna.utils import identifier_from_string
+from Bio.Seq import Seq as BPSeq
 import copy
 import operator
 import os
@@ -453,9 +453,23 @@ class Dseqrecord(SeqRecord):
             feature.location += len(nucleotides)
         return newseq
 
-    def format(self, f="gb"):
+    def format(self, format: str = "gb"):
         """Returns the sequence as a string using a format supported by Biopython
         SeqIO [#]_. Default is "gb" which is short for Genbank.
+        Allowed Formats are for example:
+
+        * "fasta": The standard FASTA format.
+        * "fasta-2line": No line wrapping and exactly two lines per record.
+        * "genbank" (or "gb"): The GenBank flat file format.
+        * "embl": The EMBL flat file format.
+        * "imgt": The IMGT variant of the EMBL format.
+
+        The format string can be modified with the keyword "dscode" if
+        the underlying dscode string is desired in the output. for example:
+        ::
+
+            Dseqrecord("PEXIGATCQFZJ").format("fasta-2line dscode")
+
 
         Examples
         --------
@@ -477,6 +491,12 @@ class Dseqrecord(SeqRecord):
         ORIGIN
                 1 aaa
         //
+        >>> print(Dseqrecord("PEXIGATCQFZJ").format("fasta-2line"))
+        >id description
+        GATCGATCGATC
+        >>> print(Dseqrecord("PEXIGATCQFZJ").format("fasta-2line dscode"))
+        >id description
+        PEXIGATCQFZJ
 
 
         References
@@ -486,13 +506,19 @@ class Dseqrecord(SeqRecord):
 
 
         """
-
         record = copy.deepcopy(self)
-        if f in ("genbank", "gb") and self.circular:
+        if "dscode" in format:
+            format = format.replace("dscode", "")
+            obj = BPSeq("")
+            obj._data = record.seq._data
+            record.seq = obj
+        format = format.strip(" -")
+        if format in ("genbank", "gb") and self.circular:
             record.annotations["topology"] = "circular"
         else:
             record.annotations["topology"] = "linear"
-        return SeqRecord.format(record, f).strip()
+
+        return SeqRecord.format(record, format).strip()
 
     def write(self, filename=None, f="gb"):
         """Writes the Dseqrecord to a file using the format f, which must
@@ -851,15 +877,6 @@ class Dseqrecord(SeqRecord):
             return self.apply_cut(cut, cut)
         else:
             answer = Dseqrecord("")
-        identifier = "part_{id}".format(id=self.id)
-        if answer.features:
-            sf = max(answer.features, key=len)  # default
-            if "label" in sf.qualifiers:
-                identifier = " ".join(sf.qualifiers["label"])
-            elif "note" in sf.qualifiers:
-                identifier = " ".join(sf.qualifiers["note"])
-        answer.id = identifier_from_string(identifier)[:16]
-        answer.name = identifier_from_string("part_{name}".format(name=self.name))[:16]
         return answer
 
     def __eq__(self, other):

{pydna-5.5.5 → pydna-5.5.7}/src/pydna/seqrecord.py

@@ -21,7 +21,7 @@ from pydna.seq import ProteinSeq
 from pydna.common_sub_strings import common_sub_strings
 
 from Bio.Data.CodonTable import TranslationError
-from Bio.SeqRecord import SeqRecord
+from Bio.SeqRecord import SeqRecord as BioSeqRecordSeqRecord
 from Bio.SeqFeature import SimpleLocation
 from Bio.SeqFeature import CompoundLocation
 from pydna.seq import Seq
@@ -37,7 +37,7 @@ from warnings import warn
 import datetime
 
 
-class SeqRecord(SeqRecord):
+class SeqRecord(BioSeqRecordSeqRecord):
     """
     A subclass of the Biopython SeqRecord class.
 
@@ -86,7 +86,7 @@ class SeqRecord(SeqRecord):
         self.annotations = {ps(k): ps(v) for k, v in self.annotations.items()}
 
     @classmethod
-    def from_Bio_SeqRecord(clc, sr: SeqRecord):
+    def from_Bio_SeqRecord(clc, sr: BioSeqRecordSeqRecord):
         """Creates a pydnaSeqRecord from a Biopython SeqRecord."""
         # https://stackoverflow.com/questions/15404256/changing-the-\
         # class-of-a-python-object-casting

pydna-5.5.5/src/pydna/download.py

@@ -1,23 +0,0 @@
-#!/usr/bin/env python3
-# -*- coding: utf-8 -*-
-# Copyright 2013-2023 by Björn Johansson.  All rights reserved.
-# This code is part of the Python-dna distribution and governed by its
-# license.  Please see the LICENSE.txt file that should have been included
-# as part of this package.
-"""Provides a function for downloading online text files."""
-
-import textwrap
-
-from pydna._pretty import pretty_str as ps
-
-
-def download_text(url):
-    """docstring."""
-    import requests
-
-    req = requests.get(url)
-
-    result = textwrap.dedent(req.text).strip()
-    result = result.replace("\r\n", "\n").replace("\r", "\n")
-
-    return ps(result)