docs/notebooks/Example_Gibson.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Example of a Gibson Assembly in pydna\n",
    "> Visit the full library documentation [here](https://pydna-group.github.io/pydna/)\n",
    "\n",
    "This example showcases a workflow of modelling Gibson assembly to clone gene fragments into plasmids for synthetic biology. The biological example is sourced [here](https://www.nature.com/articles/nmeth.1318#MOESM319), from the original Gibson assembly paper. This example constructs a synthetic pCC1BAC plasmid by joining sequence fragments from Ruminiclostridium (Clostridium) cellulolyticum. The R. cellulolyticum fragments joined are termed F1, F2, and F3, as in the paper.\n",
    "\n",
    "Source files can be found alongside this notebook, if you would like to follow along. Annotations are made alongside the code to describe key steps.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a target=\"_blank\" href=\"https://colab.research.google.com/github/pydna-group/pydna/blob/master/docs/notebooks/Example_Gibson.ipynb\">\n",
    "  <img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/>\n",
    "</a>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "%%capture\n",
    "# Install pydna (only when running on Colab)\n",
    "import sys\n",
    "if 'google.colab' in sys.modules:\n",
    "    %pip install pydna\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Importing all necessary classes and methods\n",
    "\n",
    "from pydna.parsers import parse\n",
    "from pydna.tm import tm_default\n",
    "from pydna.amplify import pcr\n",
    "from pydna.dseqrecord import Dseqrecord\n",
    "from pydna.assembly import Assembly\n",
    "from pydna.genbank import Genbank\n",
    "from pydna.gel import gel\n",
    "from pydna.ladders import GeneRuler_1kb_plus\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'molecule_type': 'DNA',\n",
       " 'topology': 'circular',\n",
       " 'data_file_division': 'BCT',\n",
       " 'date': '25-AUG-2017',\n",
       " 'accessions': ['CP001348', 'AAVC01000000', 'AAVC01000001-AAVC01000121'],\n",
       " 'sequence_version': 1,\n",
       " 'keywords': [''],\n",
       " 'source': 'Ruminiclostridium cellulolyticum H10',\n",
       " 'organism': 'Ruminiclostridium cellulolyticum H10',\n",
       " 'taxonomy': ['Bacteria',\n",
       "  'Bacillota',\n",
       "  'Clostridia',\n",
       "  'Eubacteriales',\n",
       "  'Oscillospiraceae',\n",
       "  'Ruminiclostridium'],\n",
       " 'references': [Reference(title='Complete sequence of Clostridium cellulolyticum H10', ...),\n",
       "  Reference(title='Direct Submission', ...)],\n",
       " 'comment': 'URL -- http://www.jgi.doe.gov\\nJGI Project ID: 4002584\\nSource DNA and bacteria available from Jizhong Zhou\\n(jzhou@rccc.ou.edu)\\nContacts: Jizhong Zhou (jzhou@rccc.ou.edu)\\n          David Bruce (microbe@cuba.jgi-psf.org)\\nAnnotation done by JGI-ORNL and JGI-PGF\\nFinishing done by JGI-LANL\\nFinished microbial genomes have been curated to close all gaps with\\ngreater than 98% coverage of at least two independent clones. Each\\nbase pair has a minimum q (quality) value of 30 and the total error\\nrate is less than one per 50000.\\nThe JGI and collaborators endorse the principles for the\\ndistribution and use of large scale sequencing data adopted by the\\nlarger genome sequencing community and urge users of this data to\\nfollow them. it is our intention to publish the work of this\\nproject in a timely fashion and we welcome collaborative\\ninteraction on the project and analysis.\\n(http://www.genome.gov/page.cfm?pageID=10506376).'}"
      ]
     },
     "execution_count": null,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Reading the R. cellulolyticum genome from GenBank\n",
    "gb = Genbank(\"example@example.com\")\n",
    "genome = gb.nucleotide(\"CP001348.1\")\n",
    "# Print the info of the genome\n",
    "genome.annotations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'molecule_type': 'DNA',\n",
       " 'topology': 'circular',\n",
       " 'data_file_division': 'SYN',\n",
       " 'date': '29-AUG-2024',\n",
       " 'accessions': ['.'],\n",
       " 'keywords': [''],\n",
       " 'source': 'synthetic DNA construct',\n",
       " 'organism': 'synthetic DNA construct',\n",
       " 'taxonomy': [],\n",
       " 'references': [Reference(title='Direct Submission', ...),\n",
       "  Reference(title='Direct Submission', ...),\n",
       "  Reference(title='Direct Submission', ...),\n",
       "  Reference(title='Direct Submission', ...)],\n",
       " 'comment': 'SGRef: number: 1; type: \"Journal Article\"; journalName: \"Submitted\\n(23-AUG-2007) 726 Post Road, Madison, WI 53713, USA\"\\nSGRef: number: 2; type: \"Journal Article\"\\nSGRef: number: 3; type: \"Journal Article\"'}"
      ]
     },
     "execution_count": null,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Reading the plasmid\n",
    "vector = parse(\"./pCC1BAC.gb\")[0]\n",
    "vector.annotations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Importing pre-designed primers for the PylRS insert fragment. \n",
    "\n",
    "F1_For = \"GCAGCTTCAAGTCCTGCAAACAAGGTGTACCAGGATCGTT\"\n",
    "F1_Rev = \"GATTTCAGTGTAGTTAGGGCCAGTTGAATTCAAACCTGCC\"\n",
    "F2_For = \"GGCAGGTTTGAATTCAACTGGCCCTAACTACACTGAAATC\"\n",
    "F2_Rev = \"CTTGGTGCCATCAGCATTGTTCTCTGTACCGCCCACTGTC\"\n",
    "F3_For = \"GACAGTGGGCGGTACAGAGAACAATGCTGATGGCACCAAG\"\n",
    "F3_Rev = \"CAGTTGAATAATCATGTGTTCCTGCGGCAAATGCAGTACC\"\n",
    "BACF1_For = \"AACGATCCTGGTACACCTTGTTTGCAGGACTTGAAGCTGCgcggccgcgatcctctagagtcgacctg\"\n",
    "BACF3_Rev = \"GGTACTGCATTTGCCGCAGGAACACATGATTATTCAACTGgcggccgccgggtaccgagctcgaattc\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "5210\n",
      "5384\n",
      "5172\n",
      "8221\n"
     ]
    }
   ],
   "source": [
    "# Getting the PCR products from the genome (might take a while since the genome is large)\n",
    "\n",
    "pcr_product_F1 = pcr(F1_For, F1_Rev, genome, limit=20)\n",
    "pcr_product_F2 = pcr(F2_For, F2_Rev, genome, limit=20)\n",
    "pcr_product_F3 = pcr(F3_For, F3_Rev, genome, limit=20)\n",
    "pcr_product_BAC = pcr(BACF1_For, BACF3_Rev, vector, limit=20)\n",
    "\n",
    "# Printing out the PCR fragment sizes\n",
    "print(len(pcr_product_F1))\n",
    "print(len(pcr_product_F2))\n",
    "print(len(pcr_product_F3))\n",
    "print(len(pcr_product_BAC))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCAFoAlgDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD5/ooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAK6/wCH3hWx8WeIrfT7+W4jhlcKWgZQ36g1yFel/Bf/AJHey/66D+dAGH8QfCtj4T8RXGn2EtxJDE5UNOylv0ArkK9L+NH/ACO97/10P8680oAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAK9L+C/wDyO9l/10H8680r0v4L/wDI72X/AF0H86AD40f8jve/9dD/ADrzSvS/jR/yO97/ANdD/OvNKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACvS/gv/wAjvZf9dB/OvNK9L+C//I72X/XQfzoAPjR/yO97/wBdD/OvNK9L+NH/ACO97/10P8680oAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAK7n4Y+INL8PeKbW91W6+z26OCz+Wz4H0UE1w1FAHc/E7xBpfiHxTdXulXX2i3dyVfy2TI+jAGuGoooAKKKKACinxQyzuEijeRz0VFJJq7faDrOmIj6hpN/aJINyNcWzxhh6jIGaAM+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDZ8Mf8AIat/98V7X8b/APkCaN/16JXinhj/AJDVv/viva/jf/yBNG/69EoA+fD1NJSnqaSgAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDZ8Mf8hq3/3xXtfxv/5Amjf9eiV4p4Y/5DVv/viva/jf/wAgTRv+vRKAPnw9TSUp6mkoAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooA2fDH/Iat/8AfFe1/G//AJAmjf8AXoleKeGP+Q1b/wC+K9r+N/8AyBNG/wCvRKAPnw9TSUp6mkoAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAoooHWgCa1u57OdZrd9kinIOAf51t65448R+JIIYdW1H7RHAgSMeRGm1R2+VRXY/Dn4baN4v07ULjULm/ie2gaRBbyIoJHrlTXnusWMWn6jLbxM7IjEAuQTQBn0UUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUDrRQOtAH0H8EP8AkB6z/wBej14r4n/5Ddx/vmvavgh/yA9Z/wCvR68V8T/8hu4/3zQBjUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUDrRQOtAH0H8EP+QHrP/Xo9eK+J/wDkN3H++a9q+CH/ACA9Z/69HrxXxP8A8hu4/wB80AY1FFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFA60UDrQB9B/BD/kB6z/16PXivif/AJDdx/vmvSvhZ448OeG9K1OHVtR+zyT27JGPIkfcx7fKpry/XbuC81Saa3ffGzEg4I/nQBmUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAPihlncJFG8jnoqKSTV2+0HWdMRH1DSb+0SQbka4tnjDD1GQM1Z8Mf8hq3/3xXtfxv/5Amjf9eiUAfPdFKeppKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKANnwx/yGrf/AHxXtfxv/wCQJo3/AF6JXinhj/kNW/8Aviva/jf/AMgTRv8Ar0SgD58PU0lKeppKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKANnwx/yGrf8A3xXtfxv/AOQJo3/XoleKeGP+Q1b/AO+K9r+N/wDyBNG/69EoA+fD1NJSnqaSgAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDZ8Mf8AIat/98V7X8b/APkCaN/16JXgNrdz2c6zW77JFOQcA/zrb1zxx4j8SQQw6tqP2iOBAkY8iNNqjt8qigDnz1NJRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAKBk4qaeFYgpUnkZ5qFfvCrV391PpQBUqWCMSyBWJx7VFViz/wBePrQAyeMRSFVJx71FVi8/15+tV6ACrVtbJMrFiwwM8VVq/Y/6t/8AdoApONrECm06T/WGm0AFdV4b8NWesW1xJcSzqY0LDy2A5/EGuVr0PwN/x43v/XE0AcJeQLb3LxoSQDjmq9XNU/4/pP8AeNU6AJ7aFZpQrEgH0ouYVhlKqSQPWpLD/Xr9aL//AF7fWgCpRRRQBLBGJJApJx7VLeWyW74QsfrTLT/XD61a1P8A1n4UAZ1FFFABRRRQAVLbxLLIFYkA+lRVZs/9ev1oAfe2kds+ELH6mqdaeq/638KzKACiiigB0ahnAPevR9C+Huk6n4XutTmuL1Z4QCqo6BT9crn9a85h/wBav1r3Lwj/AMk/1H/dFAHid9bpbXLxoWIBxzVar+rf8f0n1qhQAUUUUAFFFFABWppemw3r7ZWkA/2SP8Ky66Dw/wD638KAMm9tktpiiFiB61VrQ1T/AI+W+tZ9ABRRRQAUUUUAaWl6fFezKkjOAT/CR/hRqmnxWUzJGzkA/wARH+FWvD//AB9J9aPEH/H0/wBaAMOiiigAooooA6zwZ4XsvEWpxW13LcIjsATCyg/qDXT/ABJ+GujeDrjy9Pub+UbQf9IkRu3soqn8Lv8AkPW/+8K7746f8fo/3B/KgD57YYYihRlgKV/vmhPvCgDodK0C1voHeWSYFVyNrD/Cse9tUtpWRCxA9a63w7/x6Tf7hrmNV/4+W+tAGdSqMsAaSnR/fFAHS6R4cs7+3eSWScFVyNjAf0rF1C0jtZmRCxA/vGuy8Nf8eUv+4a5TWf8Aj6f60AZdFFFAGhplhFeTKkjOAT/CRW/rXhWx062jkhluGLLk72U/yFZWg/8AH0n1rsfFP/HjD/uCgDzaRQrkCm0+b/WGmUAKo3MAe9dt4V8Hafrk6JczXSBjz5TKP5qa4qP/AFi/WvVvh5/x9RfUUAc94q8Hafoc7pbTXThTx5rKf5KK4lhtYgdq9X+If/H1L9TXlMn+sb60ANooooAKKKKAFX7wq1d/dT6VVX7wq1d/dT6UAVKsWf8Arx9ar1Ys/wDXj60AF5/rz9ar1YvP9efrVegAq/Y/6t/92qFX7H/Vv/u0AUpP9YabTpP9YabQAV6H4G/48b3/AK4mvPK9D8Df8eN7/wBcTQBxGqf8f0n+8ap1c1T/AI/pP941ToAt2H+vX60X/wDr2+tFh/r1+tF//r2+tAFSiiigCe0/1w+tWtT/ANZ+FVbT/XD61a1P/WfhQBnUUUUAFFFFABVmz/16/Wq1WbP/AF6/WgC1qv8ArfwrMrT1X/W/hWZQAUUUUAPh/wBav1r3Lwj/AMk/1H/dFeGw/wCtX617l4R/5J/qP+6KAPGdW/4/pPrVCr+rf8f0n1qhQAUUUUAFFFFABXQeH/8AW/hXP10Hh/8A1v4UAZ+qf8fLfWs+tDVP+PlvrWfQAUUUUAFFFFAG54f/AOPpPrR4g/4+n+tHh/8A4+k+tHiD/j6f60AYdFFFABRRRQB6P8Lv+Q9b/wC8K7746f8AH6P9wfyrgfhd/wAh63/3hXffHT/j9H+4P5UAfPj/AHzQn3hQ/wB80J94UAdt4d/49Jv9w1zGq/8AHy31rp/Dv/HpN/uGuY1X/j5b60AZ1Oj++KbTo/vigDvfDX/HlL/uGuU1n/j6f611fhr/AI8pf9w1yms/8fT/AFoAy6KKKANrQf8Aj6T612Pin/jxh/3BXHaD/wAfSfWux8U/8eMP+4KAPN5v9YaZT5v9YaZQA6P/AFi/WvVvh5/x9RfUV5TH/rF+terfDz/j6i+ooAPiH/x9S/U15TJ/rG+terfEP/j6l+prymT/AFjfWgBtFFFABRRRQAq/eFWrv7qfSqq/eFWrv7qfSgCpViz/ANePrVerFn/rx9aAC8/15+tV6sXn+vP1qvQAVfsf9W/+7VCr9j/q3/3aAKUn+sNNp0n+sNNoAK9D8Df8eN7/ANcTXnleh+Bv+PG9/wCuJoA4jVP+P6T/AHjVOrmqf8f0n+8ap0AW7D/Xr9aL/wD17fWiw/16/Wi//wBe31oAqUUUUAT2n+uH1q1qf+s/Cqtp/rh9atan/rPwoAzqKKKACiiigAqzZ/69frVarNn/AK9frQBa1X/W/hWZWnqv+t/CsygAooooAfD/AK1frXuXhH/kn+o/7orw2H/Wr9a9y8I/8k/1H/dFAHjOrf8AH9J9aoVf1b/j+k+tUKACiiigAooooAK6Dw//AK38K5+ug8P/AOt/CgDP1T/j5b61n1oap/x8t9az6ACiiigAooooA3PD/wDx9J9aPEH/AB9P9aPD/wDx9J9aPEH/AB9P9aAMOiiigAooooA9H+F3/Iet/wDeFd98dP8Aj9H+4P5VwPwu/wCQ9b/7wrvvjp/x+j/cH8qAPnx/vmhPvCh/vmhPvCgDtvDv/HpN/uGuY1X/AI+W+tdP4d/49Jv9w1zGq/8AHy31oAzqdH98U2nR/fFAHe+Gv+PKX/cNcprP/H0/1rq/DX/HlL/uGuU1n/j6f60AZdFFFAG1oP8Ax9J9a7HxT/x4w/7grjtB/wCPpPrXY+Kf+PGH/cFAHm83+sNMp83+sNMoAdH/AKxfrXq3w8/4+ovqK8pj/wBYv1r1b4ef8fUX1FAB8Q/+PqX6mvKZP9Y31r1b4h/8fUv1NeUyf6xvrQA2iiigAooooAVThhU9xKjqu05wPSq9FABU1s6xyhmOBUNFAE1y6ySllORUNFFABVu1njjRg7YJHpVSigBzkFyR0ptFFABXZeFNd03TbS6S7ufLZ4yqjYxyfwFcbRQBZv5UmuneNsqTwcVWoooAsWkqRTBnOB9KLuVJZiyHI+lV6KACiiigCW3dUlBY4FWL6eOZ8xtkfQ1SooAKKKKACiiigAqa2kWOUFjgVDRQBf1C4inkzG+4fQiqFFFABRRRQA+NgrgnpmvVfDvjLQLDwfeWFzf7LmRQETyZDn8QuK8nooAt6hNHPdO8bblJ4OMVUoooAKKKKACiiigArY0e+trWTM0m0Y/uk/yrHooAuX88c87NG24E+mKp0UUAFFFFABRRRQBq6PeW9rcK00m1QeuCf5UaxeW91cM0Mm5SeuCP51lUUAFFFFABRRRQB2vgLXtM0XV4Z9QufJjVgS3ls38ga6/4r+OfDniW6D6RqP2ldoGfIkTt/tKK8booAViCxIoU4YE0lFAHVaLrFhaW8iTz7GK4A2Mf5CsLUJ4p52aNtwJ9CKpUUAFOQgMCabRQB1+h63p1nbSJPcbGKkAbGP8AIVz+p3MVxOzRPuBPoRVCigAooooA1NJu4LWdWmfaAeuCf5V0uv8AiDS720iS3ut7KgBHlsP5iuGooAfIwZyR0plFFADkIDgnpXoHg3xNpGk3EbXt35Sg8ny3b+QNee0UAeheMvE2katcSNZXfmqTwfLdf5gV5+5BckdKbRQAUUUUAFFFFABUzJbiyidZWNyZHEke3hUAXac9ySX4/wBketP057WPU7R76J5bRZkM8afeePcNwHI5IzXoGraZcWWi3V3qttYvpzwuLdbfw/JayLIVIjPmmFAAG2k5dsgEc5oA82oore8JJBNq7wPAstzLCVtDJatcokuVOWjUEsNocfdbkg4OKAMe6S3SVRbStJH5cZLMuCHKAuPoG3DPfFQ12njGyk07ToYNUgtv7TeYNFLa6S9ioiAbcGDRRbiSUx8pxg881xdABRRVmx06+1S4+z6fZXF3PjPl28TSNj1wATQBXZWRyjqVZTggjBBpKvaxb6pb6nM2sW1zb3szGWRbmIxuxYkk4IHU5qjQAUUVZsdOvtUuPs+n2Vxdz4z5dvE0jY9cAE0AV2Vkco6lWU4IIwQaSr2sW+qW+pzNrFtc297MxlkW5iMbsWJJOCB1Oao0AFKVZQpKkBhlSR1GccfkaWOOSaVYokZ5HIVVUZLE9gK09V0rXrC2tTq2mX9rBGnlQG5tmjGCzPgEgZ5Zj+NAGVRRRQAUu1ggcqdpJAOOCR1/mPzpOprXvNG8Q2ekRNe6VqMGno7SJJNaukYZwoJ3EY5Cr37UAZFFFFABRRRQAUUU6OOSaVYokZ5HIVVUZLE9gKAEKsoUlSAwypI6jOOPyNJWrqula9YW1qdW0y/tYI08qA3Ns0YwWZ8AkDPLMfxrKoAKKKKAClVWdwiKWZjgADJJpK0tFsNZu75J9Esby5ubZ1lU2sDSmNgcg4APcd6AM2iprq0uLK5e3u7eW3nQ4eKVCjL9QeRUNABRRRQAUUUUAFFFWbHTr7VLj7Pp9lcXc+M+XbxNI2PXABNAFdlZHKOpVlOCCMEGkq9rFvqlvqczaxbXNvezMZZFuYjG7FiSTggdTmqNABRRRQAUUUUAFFFTTpbrFbGGVnkaMmZSuAj72AA9RtCnPufSgCGiiigAooooAKKKmnS3WK2MMrPI0ZMylcBH3sAB6jaFOfc+lAENFFFABRU1mlvJewJdytFbNIolkVdxRM8kDuQM8VDQAUUUUAFFTWaW8l7Al3K0Vs0iiWRV3FEzyQO5AzxUNABRRRQAUUVMyW4sonWVjcmRxJHt4VAF2nPckl+P9ketAENFFFABRRRQAUUrKyHDKVOAcEY4IyKSgAooooAKKKKAL+iPdR6/pz2MSy3i3URgjf7rybhtB5HBOO9dlJHpMEGoz6JcQ3OqS2swniOpyygRlG8wqGt4w5C7iP3jdM84ritJTzdZsY/tf2PdcRr9pzjycsPnzkdOvXtXomqLfS6TfrNo+qaPGLeRjqU9vbRpcYUnYXSBGO/G3iRslu4oA8vroPCTMLy+RxGLKS0K3kjzmHy4t6chlDEHdsGArZzjBzXP10fgtJW1edoXnaRLdmFrAEL3fzKPKAcMp67sFW+50JoAueIEsbfw35GhyR3GmG7Rp5ftjzOkux9gw8MW0Eb+Qpzt68VyFdz4yS5bRIpLjT73Rdtwqrp11DDF5uVbMoWOKLO3GCSp++MHqK4agArqPCohk07U4NRaOLSXeHz5muWhZZAH2KCschbI3nGw/dzxiuXrrvBCT7dQkt4bq/dfLU6ZbRxSNcA7vnKyJICEx/cY/OOnNAEPioQx6dpkGnNHLpKPN5Ey3LTM0hCb1JaOMrgbDjYPvZ5zXL113jdJ9mnyXEN1YO3mKNMuY4o2twNvzhY0jAD5/uKfkPXiuRoAK6jwqIZNO1ODUWji0l3h8+ZrloWWQB9igrHIWyN5xsP3c8Yrl667wQk+3UJLeG6v3Xy1OmW0cUjXAO75ysiSAhMf3GPzjpzQBD4qEMenaZBpzRy6SjzeRMty0zNIQm9SWjjK4Gw42D72ec1y9dd43SfZp8lxDdWDt5ijTLmOKNrcDb84WNIwA+f7in5D14rkaANjwtJcR+IYGtkjd9kofzJDGojMbCRtw5XCbjkdMdD0rf1KLS7Xw9qa+Hp4rlJUT7bm+klaNBIhDKrQQj7+wZ+bGe2c1geFkaTxHaKl09s43srxlQzEIx2Ddxl8bOePm5z0rrPEq3knhy8a40jUdBRNjBLm3t4VvTvUbf3cERYjO/8AiHynpwaAPO6KKKAL+iPdR6/pz2MSy3i3URgjf7rybhtB5HBOO9dlJHpMEGoz6JcQ3OqS2swniOpyygRlG8wqGt4w5C7iP3jdM84ritJTzdZsY/tf2PdcRr9pzjycsPnzkdOvXtXomqLfS6TfrNo+qaPGLeRjqU9vbRpcYUnYXSBGO/G3iRslu4oA8vooooAKKKKACtjwtJcR+IYGtkjd9kofzJDGojMbCRtw5XCbjkdMdD0rHrY8LI0niO0VLp7ZxvZXjKhmIRjsG7jL42c8fNznpQBv6lFpdr4e1NfD08VykqJ9tzfSStGgkQhlVoIR9/YM/NjPbOa4ivRPEq3knhy8a40jUdBRNjBLm3t4VvTvUbf3cERYjO/+IfKenBrzugAooooAK6/w+ljceG/I1ySO30wXbtBL9seF3l2JvGEhl3ADZyVGN3XmuQruPB1rey6fAbeVJrZ7m4NzayQQThdkSMhVJVYKzklN2McD0oAx/FrMbyxRBGbKO0C2ciTmbzIt78lmVSTu3jBVcYxgYrn66Hxe5m1Czud8wFxaJKtvMULWylmAT5FVQCAHACrw447nnqACiiigAooooAK6jwqIZNO1ODUWji0l3h8+ZrloWWQB9igrHIWyN5xsP3c8Yrl667wQk+3UJLeG6v3Xy1OmW0cUjXAO75ysiSAhMf3GPzjpzQBD4qEMenaZBpzRy6SjzeRMty0zNIQm9SWjjK4Gw42D72ec1y9dd43SfZp8lxDdWDt5ijTLmOKNrcDb84WNIwA+f7in5D14rkaACiiigAooooAtadYTanepawFFdgzFpG2qiqpZmJ9AoJ/Cr91oCxWU91Z6vp2opbgNOtt5qtGpYKGxIiZG5lHy56ik8LI0niO0VLp7ZxvZXjKhmIRjsG7jL42c8fNznpXWeJVvJPDl41xpGo6CibGCXNvbwreneo2/u4IixGd/8Q+U9ODQB53RRRQAUUUUAWtOsJtTvUtYCiuwZi0jbVRVUszE+gUE/hV+60BYrKe6s9X07UUtwGnW281WjUsFDYkRMjcyj5c9RSeFkaTxHaKl09s43srxlQzEIx2Ddxl8bOePm5z0rrPEq3knhy8a40jUdBRNjBLm3t4VvTvUbf3cERYjO/8AiHynpwaAPO6KKKANXT9E+2WX2251Ky0+1MhiSS6Mh8x1ALALGjNwGXJxjkc1BqemSaXNErTQXEU0fmwzwMSkiZK5GQCOVYYIB4rq/BqXK6JLJb6fe61uuGVtOtYYZfKwq4lKyRS43ZwCFH3Dk9BWT40SVdXgaZ51ke3VjazhA9p8zDyiECqOm7AVfv8AQGgDnKKKKANXT9E+2WX2251Ky0+1MhiSS6Mh8x1ALALGjNwGXJxjkc1BqemSaXNErTQXEU0fmwzwMSkiZK5GQCOVYYIB4rq/BqXK6JLJb6fe61uuGVtOtYYZfKwq4lKyRS43ZwCFH3Dk9BWT40SVdXgaZ51ke3VjazhA9p8zDyiECqOm7AVfv9AaAOcooooAmtLWa+vYLS2TfPPIsUa5A3MxwBk+5rXfw2pimNprel3txDG0r20DSh9qgsxBeNUbABPDHgcZrO0lPN1mxj+1/Y91xGv2nOPJyw+fOR069e1eiaot9LpN+s2j6po8Yt5GOpT29tGlxhSdhdIEY78beJGyW7igDy+iiigAoorT0fSTqy36Rt/pENuJYY9yjzG8xFK89flZjgc/LQBRuLma6kEk7l3CJGCf7qKFUfgoA/Coq6nxL4dsdJtJzbi7jltbtbUtcSoy3WULGSNQAVAwvBLf6xeQeK5agAooooAKKKKALOnTW9vqdrNdwefbRzI80OceYgILL+IyK6P+2tDuri7N/ZxF0t51tLm2so4BIzRMqCSFTtGGIIZeQRzu4IKKAOTrT0O80yyup5NU01L+JoHWKN2cBJeCrHY6EjjB+boxPJAFFFABqOqWd7brHb6Dp1g4fcZbZ7hmIwflPmSuMc56Z4HPWsyiigAra0nUrCx0a+jubC1vJ5riDak6NkRhZd+11IKnJj6HnjIIFFFAFTUv7KbypdM+1x7s+ZBcbW8vpja4xuB56quMd6oUUUAFb3h5tDKNHrD+Xi8t5QfKZt8S7/MTK8ru3J27UUUAHiK+sb2C1MJtHu1kmEslpZi2jMWVEY2gAFuHOcZwwBJI4waKKALmlai2lanBfLbWt0YiT5N3CJYnBBBDKeD1/DqOa6PxF4stdU8P21lZW1tAXnlkni/s22jMQIi2hJI0U9UbJAUkEA5AySigDkKKKKALWmS2sGq2c19CZrSOdGniHV0DAsO3UZFdNLr+m3cV7Dc2ulJEttIga305Uku5vmEUi4UeUF+TKggEKeGJNFFAHH0UUUAFFFFABWhod1aWWrw3F9CJYFDggxCQKxUhW2MQG2sQ208HGD1oooA1JtW0a903U2nsY4tUeIRwS28CpFL+9Rt5QcRPsVh8vBBxgck83RRQAUUUUAFbem6ppFlo8kVzoNrf6h5+5ZrmSYL5RXBXEcqYIIyODncemBkooAz9RvIL24WS3021sECbTFbNKyk5PzHzHc55x1xwOOtVKKKACiiigAooooAK19G1DSrG1vhqGjw6hcOqG2MzyhEIPzAiORDgg9cnBUcck0UUAV9S1C2vvK+z6PZadszu+yvM3mZx18yR+mO2OpznjFCiigAooooAKKKKANPQLu1sdV+0XkUMsS284VJ4vMQyGFwmV/3yvPbrxin3k+i3lo80NrPYXoxiGJvNgk55xuO9MDJ5L59qKKAMmiiigAooooA0NGexTUf+JjxbNDMhbZv2u0TKjY9mKn8K3NUvdDGjTWVq9pMqW0AtjHZbJfPyPNdpSAxUgPwSR86gD5c0UUAcnRRRQB0Xh3xJDpVrLY3Wn2UsMr7xdNYQTzwnAHHmqQy8fd49iMnNbxVrA13xJfX0awrA80nk+VbRwZjLsVLBAMtg8k5PqTiiigDGooooA6LSNasNM0RkbT7O7u2mk3rc2yvuUoPLKucldjgkqMBg2DkdKviK6s7u7tntTbvKLcC6ktrcQRSS7mOVQBQBtKD7oyQTjnJKKAMeiiigCzp01vb6nazXcHn20cyPNDnHmICCy/iMiuj/ALa0O6uLs39nEXS3nW0ubayjgEjNEyoJIVO0YYghl5BHO7ggooA5OiiigArX0O/s9PTUpLq1trp3tQkEVzEXUv5sZPIwVOwP8wIPvzRRQBFqD6PPbrNYRXVrcFsPbORJHjByyvww5wNpB6/eNZtFFABRRRQAUUUUAf/Z",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAFoCAIAAAAElhK7AAAbC0lEQVR4Ae3dTaxdVdkAYG4vUOSvRQMUpI1JIRoSTUeYiDKhECKJRutPJCFGYiIJMSZqHDkyDMCBDgwTEhOVCQwalDAg0AjG3xAIN1ZhoBICJAoItJW/CqWujwN7n+909e5zetbad629nzvQdVbXfte7nnfv+3LvOZSTTvJFgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgACB4gVWNjDD1dXVU0455eyzz77ooot27dp19dVXf/nLXw753HXXXffff//a2tqzzz576NChN99888iRIxuYp60JECBAYMACmwZ8NkcjQIAAAQKdAhphJ5EFBAgQIDBkAY1wyNV1NgIECBDoFNAIO4ksIECAAIEhC2iEQ66usxEgQIBAp4BG2ElkAQECBAgMWUAjHHJ1nY0AAQIEOgU0wk4iCwgQIEBgyAIa4ZCr62wECBAg0CmgEXYSWUCAAAECQxbQCIdcXWcjQIAAgU4BjbCTyAICBAgQGLKARjjk6jobAQIECHQKaISdRBYQIECAwJAFNMIhV9fZCBAgQKBTQCPsJLKAAAECBIYsoBEOubrORoAAAQKdAhphJ5EFBAgQIDBkAY1wyNV1NgIECBDoFNAIO4ksIECAAIEhC2iEQ66usxEgQIBAp4BG2ElkAQECBAgMWUAjHHJ1nY0AAQIEOgU0wk4iCwgQIEBgyAIa4ZCr62wECBAg0CmgEXYSWUCAAAECQxbQCIdcXWcjQIAAgU4BjbCTyAICBAgQGLKARjjk6jobAQIECHQKaISdRBYQIECAwJAFNMIhV9fZCBAgQKBTQCPsJLKAAAECBIYsoBEOubrORoAAAQKdAhphJ5EFBAgQIDBkAY1wyNV1NgIECBDoFNAIO4ksIECAAIEhC2iEQ66usxEgQIBAp4BG2ElkAQECBAgMWUAjHHJ1nY0AAQIEOgU0wk4iCwgQIEBgyAIa4ZCr62wECBAg0CmgEXYSWUCAAAECQxbQCIdcXWcjQIAAgU4BjbCTyAICBAgQGLKARjjk6jobAQIECHQKaISdRBYQIECAwJAFNMIhV9fZCBAgQKBTQCPsJLKAAAECBIYsoBEOubrORoAAAQKdAhphJ5EFBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQqExgpbJ8E6W7srKyurq6efPmrVu3bt++fdeuXbt3796zZ08Iv3fv3n379q2trT3zzDMHDhw4fPjwkSNHjh49mmhnYQgQIECgLAF/xVpZ9ZANAQIECPQsoBH2DG47AgQIEChLQCMsqx6yIUCAAIGeBTTCnsFtR4AAAQJlCWiEZdVDNgQIECDQs4BG2DO47QgQIECgLAGNsKx6yIYAAQIEehbQCHsGtx0BAgQIlCWgEZZVD9kQIECAQM8CGmHP4LYjQIAAgbIENMKy6iEbAgQIEOhZQCPsGdx2BAgQIFCWgEZYVj1kQ4AAAQI9C2iEPYPbjgABAgTKEtAIy6qHbAgQIECgZwGNsGdw2xEgQIBAWQIaYVn1kA0BAgQI9CygEfYMbjsCBAgQKEtAIyyrHrIhQIAAgZ4FNMKewW1HgAABAmUJaIRl1UM2BAgQINCzgEbYM7jtCBAgQKAsAY2wrHrIhgABAgR6FtAIewa3HQECBAiUJaARllUP2RAgQIBAzwIaYc/gtiNAgACBsgQ0wrLqIRsCBAgQ6FlAI+wZ3HYECBAgUJaARlhWPWRDgAABAj0LaIQ9g9uOAAECBMoS0AjLqodsCBAgQKBnAY2wZ3DbESBAgEBZAhphWfWQDQECBAj0LKAR9gxuOwIECBAoS0AjLKsesiFAgACBngU0wp7BbUeAAAECZQlohGXVQzYECBAg0LOARtgzuO0IECBAoCwBjbCsesiGAAECBHoW0Ah7BrcdAQIECJQloBGWVQ/ZECBAgEDPAhphz+C2I0CAAIGyBDTCsuohGwIECBDoWUAj7BncdgQIECBQloBGWFY9ZEOAAAECPQtohD2D244AAQIEyhLQCMuqh2wIECBAoGcBjbBncNsRIECAQFkCGmFZ9ZANAQIECPQsoBH2DG47AgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIZBdYyb5DMRusrKxs2rTptNNO27Jly44dO3bt2rV79+49e/aEBPfu3btv3761tbWnn3764MGDb7zxxttvv3306NFicpcIAQIECOQS8Fes5ZIVlwABAgSqENAIqyiTJAkQIEAgl4BGmEtWXAIECBCoQkAjrKJMkiRAgACBXAIaYS5ZcQkQIECgCgGNsIoySZIAAQIEcglohLlkxSVAgACBKgQ0wirKJEkCBAgQyCWgEeaSFZcAAQIEqhDQCKsokyQJECBAIJeARphLVlwCBAgQqEJAI6yiTJIkQIAAgVwCGmEuWXEJECBAoAoBjbCKMkmSAAECBHIJaIS5ZMUlQIAAgSoENMIqyiRJAgQIEMgloBHmkhWXAAECBKoQ0AirKJMkCRAgQCCXgEaYS1ZcAgQIEKhCQCOsokySJECAAIFcAhphLllxCRAgQKAKAY2wijJJkgABAgRyCWiEuWTFJUCAAIEqBDTCKsokSQIECBDIJaAR5pIVlwABAgSqENAIqyiTJAkQIEAgl4BGmEtWXAIECBCoQkAjrKJMkiRAgACBXAIaYS5ZcQkQIECgCgGNsIoySZIAAQIEcglohLlkxSVAgACBKgQ0wirKJEkCBAgQyCWgEeaSFZcAAQIEqhDQCKsokyQJECBAIJeARphLVlwCBAgQqEJAI6yiTJIkQIAAgVwCGmEuWXEJECBAoAoBjbCKMkmSAAECBHIJaIS5ZMUlQIAAgSoENMIqyiRJAgQIEMgloBHmkhWXAAECBKoQ0AirKJMkCRAgQCCXgEaYS1ZcAgQIEKhCQCOsokySJECAAIFcAhphLllxCRAgQKAKAY2wijJJkgABAgRyCWiEuWTFJUCAAIEqBDTCKsokSQIECBDIJaAR5pIVlwABAgSqENAIqyiTJAkQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQILCAwMoCawe0dGVlZXV1dfPmzVu3bt2+ffuuXbt27969Z8+ecMS9e/fu27dvbW3tmWeeOXDgwOHDh48cOXL06NEBnd5RCBAgQKAV8FestRZGBAgQIDBCAY1whEV3ZAIECBBoBTTC1sKIAAECBEYooBGOsOiOTIAAAQKtgEbYWhgRIECAwAgFNMIRFt2RCRAgQKAV0AhbCyMCBAgQGKGARjjCojsyAQIECLQCGmFrYUSAAAECIxTQCEdYdEcmQIAAgVZAI2wtjAgQIEBghAIa4QiL7sgECBAg0ApohK2FEQECBAiMUEAjHGHRHZkAAQIEWgGNsLUwIkCAAIERCmiEIyy6IxMgQIBAK6ARthZGBAgQIDBCAY1whEV3ZAIECBBoBTTC1sKIAAECBEYooBGOsOiOTIAAAQKtgEbYWhgRIECAwAgFNMIRFt2RCRAgQKAV0AhbCyMCBAgQGKGARjjCojsyAQIECLQCGmFrYUSAAAECIxTQCEdYdEcmQIAAgVZAI2wtjAgQIEBghAIa4QiL7sgECBAg0ApohK2FEQECBAiMUEAjHGHRHZkAAQIEWgGNsLUwIkCAAIERCmiEIyy6IxMgQIBAK6ARthZGBAgQIDBCAY1whEV3ZAIECBBoBTTC1sKIAAECBEYooBGOsOiOTIAAAQKtgEbYWhgRIECAwAgFNMIRFt2RCRAgQKAV0AhbCyMCBAgQGKGARjjCojsyAQIECLQCGmFrYUSAAAECIxTQCEdYdEcmQIAAgVZAI2wtjAgQIEBghAIa4QiL7sgECBAg0ApohK2FEQECBAiMUEAjHGHRHZkAAQIEWgGNsLUwIkCAAIERCmiEIyy6IxMgQIBAK6ARthZGBAgQIDBCAY1whEV3ZAIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIEBi6wMnO+lZWVTZs2nXzyyZs3bz799NPDn7722muHDx9+66233n777aNHj86sP/ZliBAuP+WUU0477bQmwhtvvPHmm2+GIPNEODbmZGaS2Kmnnhoin3HGGWHy1VdfDZH/+9//TtI73oWd85OcQ9gtW7ZccMEFO3fuvPTSS8NVjz/++D/+8Y9//vOfBw8eDBvNn/+EMSAExve9731nnnlmiPbKK6+8/vrrATNQzIm5TuZBI3xNnGe2mGiHLcLXOhF6+KPgEL6a+yFU7ayzzgr7/uc//5nUrrkrlrkxjneQsHUgCruH+3BS2YsvvvhjH/vYZZddFi55+OGH//znP//973+f1Dfc5/Pf5MfbMdP85CCrq6vh5g+1DobhOO9///vDdi+99FK4OYNnuLXCg3DkyJFQ9ByYcx7tnbvy/27LyTeQs88+e+vWrR/4wAfC5S+++OKBAwcOHTo0+ZYyeQqy3qLBbYIWbryQxuTR/uhHPxpugMsvvzyk9Pvf/z7cBvv375885iG9cFtOGNc3nEQOt9Z0RcIxzz333BD2hRdeCIedrku4tUJp1o85LRwYQ+YTxpB8YDznnHNC/PPOOy8se/7550P8l19+OWCGhCffUialnw4yz7g5SPOdalKvsNG2bdtChH/961+T7QJO8x1s0eNMZ9J8Vwy38WSv888/P5Tmoosu2rFjR1j59NNPP/vss+GpfO655yb3zOT2nhxzOlTCsb9iLSGmUAQIECBQn4BGWF/NZEyAAAECCQU0woSYQhEgQIBAfQIaYX01kzEBAgQIJBTQCBNiCkWAAAEC9QlohPXVTMYECBAgkFBAI0yIKRQBAgQI1CegEdZXMxkTIECAQEIBjTAhplAECBAgUJ+ARlhfzWRMgAABAgkFNMKEmEIRIECAQH0CGmF9NZMxAQIECCQU0AgTYgpFgAABAvUJaIT11UzGBAgQIJBQQCNMiCkUAQIECNQnoBHWVzMZEyBAgEBCAY0wIaZQBAgQIFCfgEZYX81kTIAAAQIJBTTChJhCESBAgEB9AhphfTWTMQECBAgkFNAIE2IKRYAAAQL1CWiE9dVMxgQIECCQUEAjTIgpFAECBAjUJ6AR1lczGRMgQIBAQgGNMCGmUAQIECBQn4BGWF/NZEyAAAECCQU0woSYQhEgQIBAfQIaYX01kzEBAgQIJBTQCBNiCkWAAAEC9QlohPXVTMYECBAgkFBAI0yIKRQBAgQI1CegEdZXMxkTIECAQEIBjTAhplAECBAgUJ+ARlhfzWRMgAABAgkFNMKEmEIRIECAQH0CGmF9NZMxAQIECCQU0AgTYgpFgAABAvUJaIT11UzGBAgQIJBQQCNMiCkUAQIECNQnoBHWVzMZEyBAgEBCAY0wIaZQBAgQIFCfgEZYX81kTIAAAQIJBTTChJhCESBAgEB9AhphfTWTMQECBAgkFNAIE2IKRYAAAQL1CWiE9dVMxgQIECCQUEAjTIgpFAECBAjUJ6AR1lczGRMgQIBAQgGNMCGmUAQIECBQn4BGWF/NZEyAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQIECAAAECBAgQWETgc5/73CLLI2uXjxAJaooAAQIECPQj8NBDDy250fIRlkzA5QQIECBA4HgC3f89wt/97nfHu3jO+eUjzLmRZQQIECBAgAABAgQIECAwJoHcb0CmjZ82WrTOPWwR3XehyUKSLCSNhehmFg/gCDMn8pIAgYUFcr8BmTZ+2mhRrB62iO670GQhSRaSxkJ0M4sHcISZE3lJoH+B+HuEW7ZsueWWW/bv3//8O19hcOutt4bJ+fNbPsKce+V+AzJt/LTRokQ9bBHdd6HJQpIsJI2F6GYWD+AIMyfykkApAr/85S+/853vfPCDH9z0zlcYhJdhcv78lo8w/15WEiBAgACBxALRf8yMTh5v4+ji6OTxIpgnQGBgAqW9o5kpn0xhm5shd/xmozDoc6/pfXsex381+u9///vb3/72hRdeOPmJMAzCyzA5f3LLR5h/LysJEKhC4Fvf+lZReWbKJ1PYhi53/GajMOhzr+l9ex7HG+FXv/rVbdu23X///c+98xUG4WWYnD+55SPMv5eVBAhUIVDa74Qy5ZMpbFPi3PGbjcKgz72m9zUmQIAAAQIENlpg+c98Lh9hHYOswcO+aeOnjRZl6WGL6L4LTRaSZCFpLEQ3s3gAR5g5kZcEGoENeVcy/qvRn//85y+88MI111wTfiMavsIg/GsUYbLJtXOwfIR1tsgaPOybNn7aaFGWHraI7rvQZCFJFpLGQnQziwdwhJkTeUmgESjoXcno74Wjk032M4Po4ujkzIXzvIzGiU7OE+3YNdFQ0cljrz12JnphdPLYa+eciUaLTs4ZMMeyaD7RyRy7NzGjO0Ynm0tKG0SzjU6Wlrl8CHQK3HzzzZ1rki+I/0S4/Gc+l4+wzlGzBg/7po2fNlqUpYctovsuNFlIkoWksRDdzOIBHGHmRF4SaAS+//3vN+MNHoQ3IX74wx/+5S9/Cb8gDV9hEF6GyfnTWj7COntlDR72TRs/bbQoSw9bRPddaLKQJAtJYyG6mcUDOMLMibwkQIAAAQIECBAoTCD8I2fJf9fo8umt7502ftpo0cx72CK670KThSRZSBoL0c0srusIpWWbKZ9MYZvS547fbBQGfe41ve8GjuPvES7/sbTlI6yDkjV42Ddt/LTRoiw9bBHdd6HJQpIsJI2F6GYW13WE0rLNlE+msE3pc8dvNgqDPvea3re4cfQTaNHJ46UeXRydPF6EdeajcaKT6wRZ54+ioaKT6wRp/ih6YXSyuWTRQTRadHLRyAnXR/OJTibc9NhQ0R2jk8deW8hMNNvoZAkJRxOLTvaTbXTr6ORC+UQjRCcXCtssjoaKTjaXnPAgGjY6ecJblHZh/CfC5T+WtnyEdaSyBg/7po2fNlqUpYctovsuNFlIkoWksRDdzOK6jlBatpnyyRS2KX3u+M1GYdDnXtP7FjcOvyP2qdFlPjQ7XdHlMaejRcc9bBHdd6HJQpIsJI2F6GYW13WE0rLNlE+msE3pc8dvNgqDPvea3teYAAECBAgQKEkg/BOBT43u378//MVy4SsMbr311mByYiVaHrNz3x626Myhc0EhSTZ/b8XOnTt//etfh39NNvxvGHfmX86CQiTnBCkt20z5ZArbIOeO32wUBn3uNdm30Kdy+f++/PIRpgszM84aPOyVNn7aaDMUk5c9bBHdd6HJQpJ86KGHJmnfc889n//851dXVz/xiU/cd999C51lYxcXIjknQmnZZsonU9gGOXf8ZqMw6HOvyb6FPpXRDwhFJ6f5psfRxdHJ6avmHEfjRCfnDDizLBoqOjlzYfRl9MLoZPTyeSaj0aKT80TLtCaaT3QyUwKTsM0jN711M5l161TBpzNvYkYnmz/dwEE0sehkP0lGt45OLpRPNEJ0cqGwzeJoqOhkc8kJD6Jho5MnvMXMhc0DOL1LMzmzOMdLnxqNqKb90FTaaJF0K/mUVw8OUZyZyR07dtx4442f/exnH3/88Ztuuin8RBjGJ5988syykl8WIjknUWnZZsonU9gGOXf8ZqMw6HOvyb4b/lTGG+Hy/3355SNMF2ZmnDV42Ctt/LTRZigmL3vYIrrvQpOFJPmNb3zj1Vdf/fCHPxz+N7w1uH379q985Stf+9rXFjrLxi4uRHJOhNKyzZRPprANcu74zUZh0Odek30H8FROAxp3C/ziF7/oXmQFAQI1C3jMa67ee7lfffXVTzzxRPjA5E9/+tNzzz13Mn3vvfe+9+cb/P/hJ/cf/OAH5513XqY80n6E6fXXXz906NCB977eeuutyTBh8tdee+1TTz31t7/97corrwxP4IsvvviHP/zhIx/5SMItlg9V+E21/AFFiArkflqjm64zmfbpbjbK/Zj3+Yz3/6iWdpO8W9a1tbXwbfSMM84IPyOHf6/84x//ePiDhd66PHjw4Hvf+f/f/zf3zTKD8F3+61//+mOPPXb77beH33EtEyp6bXPSJB8s/PSnPx0S/tKXvjTZqwke3frEJh999NELLrhg9+7d4dd94R2v8L7Xpz71qQceeODEomW6avmbKkliWe/MJBl2BqnrCLmf1k6umQXNA5jk6W6C537M+3zG+39US7tJ3i3rb3/7202b3n37MLyP8sgjj4T3UZobqKn9OoOf/exn3/zmN9dZsMwfTTJZWVn5zGc+85vf/OZXv/rVFVdcsUzAmWubk6b6CNNZZ51122233X333RdeeGETfGbTZV42Me+4444mTpBpxiUMlr+pkpwi652ZJMPOIHUdIffT2sk1s6B5WFI93U38rI95k3YPz3j/j2ppN8m7Nb3hhhvCt9EvfvGLk9dnnnnm3r17X3nllabknYOtW7cGzW3btnWuPIEFzT0xuTb8wBrS++Mf/3gCoaKXPPnkk5MPFoafOJsPFk4/NtGrOic/+clP/ulPf/rrX//auXLRBXfdddd3v/vdU089dXJh+AzkF77whQcffHDROFnXL39TJUkv652ZJMPOIHUdIffT2sk1syDT093skukx7/MZ7/9RLe0maap5UnhrcPpNpvDjVzmfrLvzzjvbRN8bXXzxxe8Nl/3/q6666vrrr//e97734x//+Ec/+tGHPvShsOMll1yybNyTTtq8efN11123fJyZCKeffnrI9pxzzpnM/+QnPwnvFIZPJM8s2/CXJd9UG44z1ARyP62LuuV7uptMcjzmPT/jPT+qpd0kTSkNCBAgQIAAAQIECBAgQGDoAv8DNIwUyHDTngMAAAAASUVORK5CYII=",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=600x360>"
      ]
     },
     "execution_count": null,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Making a gel to show the PCR products\n",
    "im = gel(\n",
    "    [\n",
    "        GeneRuler_1kb_plus,\n",
    "        [pcr_product_F1],\n",
    "        [pcr_product_F2],\n",
    "        [pcr_product_F3],\n",
    "        [pcr_product_BAC],\n",
    "    ]\n",
    ")\n",
    "im.rotate(90, expand=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dseqrecord\n",
      "circular: True\n",
      "size: 23827\n",
      "ID: id\n",
      "Name: name\n",
      "Description: description\n",
      "Number of features: 28\n",
      "/molecule_type=DNA\n",
      "Dseq(o23827)\n",
      "GCAG..ccgc\n",
      "CGTC..ggcg\n"
     ]
    }
   ],
   "source": [
    "# Performing the Gibson Assembly. Note that the assembly class parameters should be given as a list.\n",
    "\n",
    "assembled = Assembly([Dseqrecord(pcr_product_F1), Dseqrecord(pcr_product_F2), Dseqrecord(pcr_product_F3), Dseqrecord(pcr_product_BAC)])\n",
    "assembled_circ = assembled.assemble_circular()\n",
    "\n",
    "# Printing out the Gibson Assembly product\n",
    "print(assembled_circ[0])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": ".venv",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}