{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "cb38eda5-f260-4dd1-9f89-179e70fa379e",
   "metadata": {},
   "source": [
    "# An introductory TRGT/TRVZ tutorial"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6fd80e6f-1205-491e-a55d-e7068c055a16",
   "metadata": {},
   "source": [
    "In this tutorial, you will learn the basics of TRGT and TRVZ by analyzing a tiny example dataset included in this repository."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6b5e159b-d769-4361-9dca-65a2f39bae72",
   "metadata": {},
   "source": [
    "## Overview of the input files"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "46f2585a-49e8-411d-89fb-7a742a13bc41",
   "metadata": {},
   "source": [
    "Our example dataset consists of a reference genome (`reference.fasta`), a file with aligned reads (`sample.bam`), and the repeat definition file (`repeat.bed`) all stored in the `input/` directory."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "cf7618b1-0b2f-4b1c-9d36-791598d19c47",
   "metadata": {},
   "outputs": [],
   "source": [
    "genome  = \"input/reference.fasta\"\n",
    "repeats = \"input/repeat.bed\"\n",
    "reads   = \"input/sample.bam\""
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5d23d501-6dc3-4002-8d0c-754408c804f5",
   "metadata": {},
   "source": [
    "The repeat definitions are stored in a BED file that, among other information, contains repeat coordinates, repeat identifiers, and motifs:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8822ace1-6366-4794-9e54-308d46b82b50",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "chrA\t10001\t10061\tID=TR1;MOTIFS=CAG;STRUC=(CAG)n\n"
     ]
    }
   ],
   "source": [
    "! cat {repeats}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "46408d32-3420-49e1-a86e-b28f891e3f4a",
   "metadata": {},
   "source": [
    "## Genotype repeats with TRGT"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "30c5e89d-6be1-4af8-99e4-70f9a9f8285c",
   "metadata": {},
   "source": [
    "We can genotype the repeat by runing TRGT as shown below. Note that TRGT also has a `--karyotype` parameter that is used to specify sex of the input sample."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5f41642b-dc2f-4df3-ac36-0a113db6b287",
   "metadata": {},
   "outputs": [],
   "source": [
    "! trgt --genome {genome} --reads {reads} --repeats {repeats} --output-prefix output/sample "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5eceb27b-025f-48d1-a696-6e1f1b755aa5",
   "metadata": {},
   "source": [
    "The output consists of files `sample.vcf.gz` and `sample.spanning.bam`. The VCF file contains repeat genotypes while the BAM file contains pieces of HiFi reads that fully span the repeat sequences."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7c5203ca-c624-4ae0-b872-0c29b03e5265",
   "metadata": {},
   "source": [
    "For example, here is the first entry of the VCF file:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "863264a4-691b-4287-bb88-bddbe50defa4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "chrA\t10002\t.\tCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAG\tCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAG\t0\t.\tTRID=TR1;END=10061;MOTIFS=CAG;STRUC=(CAG)n\tGT:AL:ALLR:SD:MC:MS:AP:AM\t1/1:33,33:30-39,30-39:15,14:11,11:0(0-33),0(0-33):1,1:.,.\n"
     ]
    }
   ],
   "source": [
    "! bcftools view --no-header output/sample.vcf.gz | head -n 1"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1f72817f-c88e-4a5c-9401-3d413733f051",
   "metadata": {},
   "source": [
    "It says that:\n",
    "- There is a tandem repeat starting at position 10001 of chrA\n",
    "- The reference sequence of this repeat is CCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAG\n",
    "- This repeat has two alleles spanning 16 CAG motifs each"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bdfa1ec8-ea2f-4017-983a-7d6c07720819",
   "metadata": {},
   "source": [
    "## Sort and index the outputs"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9cf72cc8-dce9-4ad3-ac66-caf12c048207",
   "metadata": {},
   "source": [
    "TRGT outputs are not sorted. So we first sort and index the VCF:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "dab8ffc2-06f7-4e8a-8f35-d58ba2f2aea9",
   "metadata": {},
   "outputs": [],
   "source": [
    "! bcftools sort -Ob -o output/sample.sorted.vcf.gz output/sample.vcf.gz\n",
    "! bcftools index output/sample.sorted.vcf.gz"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "430d40c1-4664-492d-aee5-3709c8a5cea0",
   "metadata": {},
   "source": [
    "and then the BAM:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "52814e32-dc28-48be-9ad5-ba1b50b076f1",
   "metadata": {},
   "outputs": [],
   "source": [
    "! samtools sort output/sample.spanning.bam -o output/sample.bam\n",
    "! samtools index output/sample.bam"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "23bffe26-c3af-4cba-bb81-3c9a9a407ca1",
   "metadata": {},
   "source": [
    "And that's it! The output files are now ready for downstream analysis."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5dcbb185-2811-49ad-b660-9ded2b437118",
   "metadata": {},
   "source": [
    "## Visualize a repeat with TRVZ"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "80c7f269-cf49-4340-a810-c1866614bbbf",
   "metadata": {},
   "source": [
    "To visualize the repeat with the identifier \"TR1\", run:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f5e07518-8492-4fd1-9cf6-8b2ce42bdcec",
   "metadata": {},
   "outputs": [],
   "source": [
    "! trvz --genome {genome} --repeats {repeats} --vcf output/sample.sorted.vcf.gz --spanning-reads output/sample.bam --repeat-id TR1 --image output/TR1.png"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "57a668f1-879d-4454-a3db-f357578adad5",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwIAAAGNCAYAAACixcacAABkLklEQVR4Ae3gAZAkSZIkSRKLqpm7R0REZmZmVlVVVVV3d3d3d/fMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMzMdHd3d3dXV1VVVVVmZkZGRIS7m5kKz0xmV3d1d3dPz8zMzMxMomybq6666qqrrrrqqquuuur/E2TbXHXVVVddddVVV1111VX/nyDb5qqrrrrqqquuuuqqq676/wTZNlddddVV/412d3f5mI/5GAC+6qu+iuPHj3PVv92tt97K93zP93C/t3qrt+KlX/qluepf56//+q/5nd/5HS5evMhnf/Znc7/d3V0+53M+h1tvvZXP+qzP4qVf+qW56qqrrvpfCNk2V1111VX/jR7ykIfwUR/1Uezu7vI7v/M7/NZv/RZX/dv99m//Nu/93u/Ne7/3ewPw3u/93jz4wQ/mqn+d137t1+Z+v/3bv8393vu93xuAl37pl+ZrvuZrePrTn85VV1111f9CyLa56qqrrvpv8td//de893u/N3/9138NgCRsc9W/3W//9m/zPu/zPrz3e783D37wg3mv93ovrvq3+e3f/m0++7M/m9/+7d/mfpK4ePEix48f56Vf+qX56q/+al77tV+bq6666qr/ZZBtc9VVV1313+S3f/u3+ezP/mx++7d/G4DXfu3X5rM/+7N57dd+ba76t/nrv/5rPvqjP5rXfu3X5qd/+qd567d+az77sz+bq/71fvu3f5vP/uzP5rd/+7e5nyRsA/Dar/3afPZnfzav/dqvzVVXXXXV/zLItrnqqquu+m/y27/923zMx3wMf/VXfwXAiRMnePrTn87x48e56t/vt3/7t/nsz/5sfvu3f5ur/vV++7d/m8/+7M/mt3/7t7mfJGwD8Nqv/dp89md/Nq/92q/NVVddddX/Msi2ueqqq676b/TgBz+Yr/7qr+YZz3gG3/Vd38Vf//Vfc9W/3Xd/93dz4sQJXuqlXorv+Z7v4a/+6q/46Z/+aa761/nrv/5rfud3fofv/u7v5qu/+qt5qZd6KY4fP85rv/Zr8zqv8zq81mu9Fm/91m/NrbfeyvHjx7nqqquu+l8G2TZXXXXVVf+N/vqv/5r3fu/35vjx43z1V381L/3SL81V/3Y//dM/zVd/9Vezu7vLgx/8YL77u7+b48ePc9W/zlu/9Vuzu7vL/b76q7+al37pl+av//qv+eiP/mh2d3f57u/+bl76pV+aq6666qr/hZBtc9VVV1111VVXXXXVVVf9f4Jsm6uuuuqq/2Z//dd/DcBLv/RL82/1O7/zOzzoQQ/iwQ9+MP9Wv/M7v8ODHvQgHvzgB/Nv9Tu/8zs86EEP4sEPfjD/Fru7u/zN3/wNL/VSL8Xx48f5t9jd3eVv/uZveKmXeimOHz/Ov8Xu7i5/8zd/w0u91Etx/Phx/i12d3f5m7/5G17qpV6K48eP829x66238oxnPIPXeq3X4t/q1ltv5RnPeAav9Vqvxb/VrbfeyjOe8Qxe67Vei6uuuuqq/wOQbXPVVVdd9d/stV/7tQH47d/+bf6tJPFZn/VZfPZnfzb/VpL4rM/6LD77sz+bfytJfNZnfRaf/dmfzb/Fb//2b/M6r/M6/NZv/Rav/dqvzb/Fb//2b/M6r/M6/NZv/Rav/dqvzb/Fb//2b/M6r/M6/NZv/Rav/dqvzb/Fb//2b/M6r/M6/NZv/Rav/dqvzb/FZ3/2Z/M5n/M52Obf6rM/+7P5nM/5HGzzb/XZn/3ZfM7nfA62ueqqq676PwDZNi/A3t4eP/qjP8of/dEf8bSnPY2rrrrqqv8sf/3Xfw3AS7/0S/Nv9du//ds8+MEP5sEPfjD/Vr/927/Ngx/8YB784Afzb/Xbv/3bPPjBD+bBD34w/xa7u7v89V//NS/90i/N8ePH+bfY3d3lr//6r3npl35pjh8/zr/F7u4uf/3Xf81Lv/RLc/z4cf4tdnd3+eu//mte+qVfmuPHj/Nvceutt3Lrrbfy2q/92vxb3Xrrrdx666289mu/Nv9Wt956K7feeiuv/dqvzVVXXXXV/2QPfehDeZVXeRXe8R3fkZ2dHV4AZNs8H3t7e7z92789f/zXf00eP8FUO6666qqr/rMMT3oCIDYe8lL8Wx0+5S/oT95Ad/J6/q0On/IX9CdvoDt5Pf9Wh0/5C/qTN9CdvJ5/i1wesLzziSxufBSx2OLfIpcHLO98IosbH0Ustvi3yOUByzufyOLGRxGLLf4tcnnA8s4nsrjxUcRii3+L8cLdDBfuYvPhL8e/1XjhboYLd7H58Jfj32q8cDfDhbvYfPjL8W+1uvNJOJLulgdz1VVXXfWfpU4jsXuRV37pl+bHf/zH2dnZ4flAts3z8RVf8RV8zpd8Cf37fiCzl3ppblwsuOqqq676z/JX7/7OyOK1vub3+bf69Q94KR76Fh/MQ9/yQ/i3+vUPeCke+hYfzEPf8kP4t/r1D3gpHvoWH8xD3/JD+Le4+MQ/4y++/P15uY//dk486hX4t7j4xD/jL778/Xm5j/92TjzqFfi3uPjEP+Mvvvz9ebmP/3ZOPOoV+Le4+MQ/4y++/P15uY//dk486hX4t3jaz34TT/u5b+b1v+1v+Ld62s9+E0/7uW/m9b/tb/i3etrPfhNP+7lv5vW/7W/4t/qLL38/WrfmUZ/3RVx11VVX/We5c7lk/Td/zfCd38pnfdIn8XEf93E8H8i2eT4e9ehHc+8NN7H5AR8MwCuePMlVV1111X+W3/nIDyWmypt+wY/xb3Xv4/6UrTM3snnmRv6t7n3cn7J15kY2z9zIv9W9j/tTts7cyOaZG/m3GA73uPiMJ3DiQY+m39zh32I43OPiM57AiQc9mn5zh3+L4XCPi894Aice9Gj6zR3+LYbDPS4+4wmceNCj6Td3+Lc4PHsnB2fv5NrHviL/Vodn7+Tg7J1c+9hX5N/q8OydHJy9k2sf+4r8W/3q570XrQ68+pd+JVddddVV/1n+9MIFAA6/7Zu59q47eOITnsDzgWyb5+PMddcxvP4bsfVmbwHAK548yVVXXXXVf5bf/8SPpUw9b/gZ38NVV/1f9auf9160OvDqX/qVXHXVVVf9Z/nTCxcAOPiFn6P/9V/h7D338Hwg2+b5OHPddQyv/0ZsvdlbAHDjYsFVV1111X+WJ37Gp1DGGS/38d/BVVf9X/UXX/5+tG7Noz7vi7jqqquu+s9y53IJwMEv/Bz9r/8KZ++5h+cD2TbPx5nrrmN4/Tdi683egquuuuqq/2znP++zqeOcB3/Yt3DVVf9X3foNH8TUrTj1GZ/NVVddddV/toNf+Dn6X/8Vzt5zD88Hsm2ej0c+6lHcd+PNbH7AB3PVVVdd9Z/t/Od9NnWc8+AP+xauuur/qlu/4YOYuhWnPuOzueqqq676z3b4bd/MNXfezpOe+ESeD2TbPB9f8RVfwed8yZfQv+8HMnupl+aqq6666j/T+c/7bOo458Ef9i1cddX/Vbd+wwcxdStOfcZnc9VVV131n2n9N3/N+J3fymd+0ifxcR/3cTwfyLZ5Pvb29ni7t3s7/uRv/oY8cZKpVq666qqr/rOMz7gVZTC/8ZFcddX/Vas7n4Qj6R70YK666qqr/rPUaSIuXuCVXuql+Imf+Al2dnZ4PpBt8wLs7e3xIz/yI/zRH/0RT3va07jqqquuuuqqq6666qqr/md76EMfyqu8yqvwTu/0Tuzs7PACINvmqquuuuqqq6666qqrrvr/BNk2V1111VVXXXXVVVddddX/J8i2eSE+66/+hl++/S7uODzkqquuuuo/053v8vb8e/zZHfu8z48/iauu+p/s4rVP56qrrrrqP9NNm5u88c038Dkv81K8EMi2eQE+8y//hu944pMBuHGx4KqrrrrqP9OfvvWb8e/x53fs8ym//HSuuup/skunb+eqq6666j/TncslAO/3qEfwuS/7UrwAyLZ5PvbGiVf6mV9kbxwAeMWTJ7nqqquu+s/0U2/0uvx7/OWd+3z1793BVVf9T7Z/4m6uuuqqq/4z/emFCwDsdD1/8lZvyk5XeT6QbfN8/OF9Z3mH3/gd7veKJ09y1VVXXfWf6eMf/XL8e9yzP/Izf3+Wq676n2z/xN1cddVVV/1n+tMLF7jfj73ea/Gq15zh+UC2zfPxh/ed5R1+43e4342LBVddddVV/5me9vcn+fd4+Ru3aM1cddX/ZJdO385VV1111X+mO5dL7vdjr/davOo1Z3g+kG3zfPzhfWd5h9/4Ha666qqr/qv86Ku+Af8et18a+PzfvI2rrvqf7OK1T+eqq6666r/Kj73ea/Gq15zh+UC2zfOxN0684s/8AvvjyFVXXXXVf4U73+Xt+ff4szv2eZ8ffxJXXfU/2cVrn85VV1111X+F7a7jT9/qzdjpKs8Hsm1egM/8y7/hO574ZK666qqr/ivc+S5vz7/Hn92xz/v8+JO46qr/yS5e+3Suuuqqq/4rfOyLP5aPe4nH8gIg2+aF+My//Bt++Y47ufPwiKuuuuqq/0x3vsvb8+/xZ3fs8z4//iSuuup/sovXPp2rrrrqqv9MN25u8MY33cjnvuxL8UIg2+aqq6666qqrrrrqqquu+v8E2TZXXXXVVVddddVVV1111f8nyLZ5IT7rr/6GX779Lu44POSqq676/+OVrznDT7zea/HcvuLvHsdX/v3j+M9w4t6HcNVV/9ddvPbp/Gf42Bd/LB/3Eo/lub3db/wOf3zfWa666qr/H+5613cA4DP/8q/5nJd5KV4IZNu8AJ/5l3/DdzzxyQDcuFhw1VVX/f/x4idP8J2v+ao8t6/7hyfwfU9+Kv8Zjp27mauu+r/u0unb+c/wHo94GB/xYo/mub3v7/4hf3/hIlddddX/D3/2Nm8OwA0/+GO836Mewee+7EvxAiDb5vnYGyde6Wd+kb1xAOAVT57kqquu+v/jwdtbfNWrviLP7Vsf9yR+6fY7+M+wffF6rrrq/7r9E3fzn+FNbr6JD3zsI3luH/OHf8qt+wdcddVV/z/89Bu/HgA3/OCPsdP1/MlbvSk7XeX5QLbN8/GH953lHX7jd7jfK548yVVXXfX/x/UbC97p4Q/luf3JvffxR/ee5T/D9sXrueqq/+v2T9zNf4ZXvOYMr3rdNTy3H3/qrdxxeMhVV131/8NPv/HrAXDDD/4YAD/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4X43LhZcddVV/3+cns/5pTvv4rm9zYNu5o6DQ/4zHDt3M1dd9X/dpdO385/h9GLOL91xF8/trW66kbuWS6666qr/H/7sbd4cgBt+8McA+LHXey1e9ZozPB/Itnk+/vC+s7zDb/wOV1111f9PL3f6FB/74o/huf3anXfz3U9+Kv8ZTtz7EK666v+6i9c+nf8M7/2Ih/EGN17Pc/vKv3s8f3H+PFddddX/D3e96zsAcMMP/hgAP/Z6r8WrXnOG5wPZNs/H3jjxij/zC+yPI1ddddX/P698zRl+4vVei+f2FX/3OL7y7x/Hf4YT9z6Eq676v+7itU/nP8PHvvhj+biXeCzP7e1+43f44/vOctVVV/3/cNe7vgMAN/zgj7HddfzpW70ZO13l+UC2zQvwmX/5N3zHE5/MVVdd9f/PK19zhp94vdfiuX3F3z2Or/z7x/Gf4cS9D+Gqq/6vu3jt0/nP8LEv/lg+7iUey3N7u9/4Hf74vrNcddVV/z/c9a7vAMANP/hjfOyLP5aPe4nH8gIg2+aF+My//Bt++Y47ufPwiKuuuur/j1e+5gw/8XqvxXP7ir97HF/594/jP8OJex/CVVf9X3fx2qfzn+FjX/yxfNxLPJbn9na/8Tv88X1nueqqq/5/uOtd3wGAz/iLv+ZzX/aleCGQbXPVVVddddVVV1111VVX/X+CbJurrrrqqquuuuqqq6666v8TZNu8EJ/1V3/DL99+F3ccHnLVVVdd9Z/pxL0P4aqr/q+7eO3Tueqqq676z3TT5iZvfPMNfM7LvBQvBLJtXoDP/Mu/4Tue+GQAblwsuOqqq676z3Ts3M1cddX/dZdO385VV1111X+mO5dLAN7vUY/gc1/2pXgBkG3zfOyNE6/0M7/I3jgA8IonT3LVVVdd9Z9p++L1XHXV/3X7J+7mqquuuuo/059euADATtfzJ2/1pux0lecD2TbPxx/ed5Z3+I3f4X6vePIkV1111VX/mbYvXs9VV/1ft3/ibq666qqr/jP96YUL3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO9zvxsWCq6666qr/TMfO3cxVV/1fd+n07Vx11VVX/We6c7nkfj/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4aqrrrrqv8qJex/CVVf9X3fx2qdz1VVXXfVf5cde77V41WvO8Hwg2+b52BsnXvFnfoH9ceSqq6666r/CiXsfwlVX/V938dqnc9VVV131X2G76/jTt3ozdrrK84FsmxfgM//yb/iOJz6Zq6666qr/CifufQhXXfV/3cVrn85VV1111X+Fj33xx/JxL/FYXgBk27wQn/mXf8Mv33Endx4ecdVVV131n+nEvQ/hqqv+r7t47dO56qqrrvrPdOPmBm9804187su+FC8Esm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXXXVVVddddVVV/1/gmybF+Kz/upv+OXb7+KOw0Ouuuqqq/4znbj3IVx11f91F699Oldd9T/Z9z/sFspqxYviVyl81513c9UL9tbXnOaduuBF9S533seL6vsfdgtlteK5fcLuIW988w18zsu8FC8Esm1egM/8y7/hO574ZABuXCy46qqrrvrPdOzczVx11f91l07fzlVX/U/2lTdfTx0GXhS/bvjZ+85z1Qv2+idP8JZd8KL6yHvP86L6ypuvpw4Dz+3tbr0DgPd71CP43Jd9KV4AZNs8H3vjxCv9zC+yNw4AvOLJk1x11VVX/Wfavng9V131f93+ibu56qr/yT75utP008SL4jfG5A8u7XHVC/bKO1u8YV95UX3uuV1eVJ983Wn6aeK5vfmTngbATtfzJ2/1pux0lecD2TbPxx/ed5Z3+I3f4X6vePIkV1111VX/mbYvXs9VV/1ft3/ibq666n+yT77uNP008aL4jTH5g0t7XPWCvfLOFm/YV15Un3tulxfVJ193mn6aeG5v/qSncb8fe73X4lWvOcPzgWyb5+MP7zvLO/zG73C/GxcLrrrqqqv+Mx07dzNXXfV/3aXTt3PVVf+TfeXN11PHgRfFryf87H3nueoFe/2TJ3jLPnhRfeQ953lRfeXN11PHgef2dk+/g/v92Ou9Fq96zRmeD2TbPB9/eN9Z3uE3foerrrrqqv8qJ+59CFdd9X/dxWufzlVX/U/2/Q+7hbJa8aL4VQrfdefdXPWCvfU1p3mnLnhRvcud9/Gi+v6H3UJZrXhu73Lnfdzvx17vtXjVa87wfCDb5vnYGyde8Wd+gf1x5Kqrrrrqv8KJex/CVVf9X3fx2qdz1VX/k33/w26hrFa8KH6VwnfdeTdXvWBvfc1p3qkLXlTvcud9vKi+/2G3UFYrntu73HkfANtdx5++1Zux01WeD2TbvACf+Zd/w3c88clcddVVV/1XOHHvQ7jqqv/rLl77dK666n+y73/YLZTVihfFr1L4rjvv5qoX7K2vOc07dcGL6l3uvI8X1fc/7BbKasVze5c77wPgY1/8sXzcSzyWFwDZNi/EZ/7l3/DLd9zJnYdHXHXVVVf9Zzpx70O46qr/6y5e+3Suuup/su9/2C2U1YoXxa9S+K477+aqF+ytrznNO3XBi+pd7ryPF9X3P+wWymrFc/v43QPe+KYb+dyXfSleCGTbXHXVVVddddVVV1111VX/nyDb5qqrrrrqqquuuuqqq676/wTZNi/EZ/3V3/DLt9/FHYeHXHXVVVf9Zzpx70O46qr/6y5e+3Suuuqqq/4z3bS5yRvffAOf8zIvxQuBbJsX4DP/8m/4jic+GYAbFwuuuuqqq/4zHTt3M1dd9X/dpdO3c9VVV131n+nO5RKA93vUI/jcl30pXgBk2zwfe+PEK/3ML7I3DgC84smTXHXVVVf9Z9q+eD1XXfV/3f6Ju7nqqquu+s/0pxcuALDT9fzJW70pO13l+UC2zfPxh/ed5R1+43e43yuePMlVV1111X+m7YvXc9VV/9ftn7ibq6666qr/TH964QL3+7HXey1e9ZozPB/Itnk+/vC+s7zDb/wO97txseCqq6666j/TsXM3c9VV/9ddOn07V1111VX/me5cLrnfj73ea/Gq15zh+UC2zfPxh/ed5R1+43e46qqrrvqvcuLeh3DVVf/XXbz26Vx11VVX/Vf5sdd7LV71mjM8H8i2eT72xolX/JlfYH8cueqqq676r3Di3odw1VX/11289ulcddVVV/1X2O46/vSt3oydrvJ8INvmBfjMv/wbvuOJT+aqq6666r/CiXsfwlVX/V938dqnc9VVV131X+FjX/yxfNxLPJYXANk2L8Rn/uXf8Mt33Mmdh0dcddVVV/1nOnHvQ7jqqv/rLl77dK666qqr/jPduLnBG990I5/7si/FC4Fsm6uuuuqqq6666qqrrrrq/xNk21x11VVXXXXVVVddddVV/58g2+aF+Ky/+ht++fa7uOPwkKuuuuqq/0wn7n0IV131f93Fa5/OVVdd9W/z1tec5p264EX1Lnfex/9HN21u8sY338DnvMxL8UIg2+YF+My//Bu+44lPBuDGxYKrrrrqqv9Mx87dzFVX/V936fTtXHXVVf82r3/yBG/ZBS+qj7z3PP8f3blcAvB+j3oEn/uyL8ULgGyb52NvnHiln/lF9sYBgFc8eZKrrrrqqv9M2xev56qr/q/bP3E3V1111b/NK+9s8YZ95UX1ued2+f/oTy9cAGCn6/mTt3pTdrrK84Fsm+fjD+87yzv8xu9wv1c8eZKrrrrqqv9M2xev56qr/q/bP3E3V1111b/NK+9s8YZ95UX1ued2+f/oTy9c4H4/9nqvxatec4bnA9k2z8cf3neWd/iN3+F+Ny4WXHXVVVf9Zzp27mauuur/ukunb+eqq676t3n9kyd4yz54UX3kPef5/+jO5ZL7/djrvRaves0Zng9k2zwff3jfWd7hN36Hq6666qr/KifufQhXXfV/3cVrn85VV131b/PW15zmnbrgRfUud97H/3c/9nqvxatec4bnA9k2z8feOPGKP/ML7I8jV1111VX/FU7c+xCuuur/uovXPp2rrrrq3+atrznNO3XBi+pd7ryP/8+2u44/fas3Y6erPB/ItnkBPvMv/4bveOKTueqqq676r3Di3odw1VX/11289ulcddVV/zZvfc1p3qkLXlTvcud9/H/2sS/+WD7uJR7LC4BsmxfiM//yb/jlO+7kzsMjrrrqqqv+M5249yFcddX/dRevfTpXXXXVv81bX3Oad+qCF9W73Hkf/x/duLnBG990I5/7si/FC4Fsm6uuuuqqq6666qqrrrrq/xNk21x11VVXXXXVVVddddVV/58g2+aF+Ky/+ht++fa7uOPwkKuuuuqq/0wn7n0IV131f93Fa5/OVVddddV/pps2N3njm2/gc17mpXghkG3zAnzmX/4N3/HEJwNw42LBVVddddV/pmPnbuaqq/6vu3T6dq666qqr/jPduVwC8H6PegSf+7IvxQuAbJvnY2+ceKWf+UX2xgGAVzx5kquuuuqq/0zbF6/nqqv+r9s/cTdXXXXVVf+Z/vTCBQB2up4/eas3ZaerPB/Itnk+/vC+s7zDb/wO93vFkye56qqrrvrPtH3xeq666v+6/RN3c9VVV131n+lPL1zgfj/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4X43LhZcddVVV/1nOnbuZq666v+6S6dv56qrrrrqP9OdyyX3+7HXey1e9ZozPB/Itnk+/vC+s7zDb/wOV1111VX/VU7c+xCuuur/uovXPp2rrrrqqv8qP/Z6r8WrXnOG5wPZNs/H3jjxij/zC+yPI1ddddVV/xVO3PsQrrrq/7qL1z6dq6666qr/Cttdx5++1Zux01WeD2TbvACf+Zd/w3c88clcddVVV/1XOHHvQ7jqqv/rLl77dK666qqr/it87Is/lo97icfyAiDb5oX4zL/8G375jju58/CIq6666qr/TCfufQhXXfV/3cVrn85VV1111X+mGzc3eOObbuRzX/aleCGQbXPVVVddddVVV1111VVX/X+CbJurrrrqqquuuuqqq6666v8TZNu8EJ/1V3/DL99+F3ccHnLVVVdd9Z/pxL0P4aqr/q+7eO3Tueqqq676z3TT5iZvfPMNfM7LvBQvBLJtXoDP/Mu/4Tue+GQAblwsuOqqq676z3Ts3M1cddX/dZdO385VV1111X+mO5dLAN7vUY/gc1/2pXgBkG3zfOyNE6/0M7/I3jgA8IonT3LVVVdd9Z9p++L1XHXV/3X7J+7mqquuuuo/059euADATtfzJ2/1pux0lecD2TbPxx/ed5Z3+I3f4X6vePIkV1111VX/mbYvXs9VV/1ft3/ibq666qqr/jP96YUL3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO9zvxsWCq6666qr/TMfO3cxVV/1fd+n07Vx11VVX/We6c7nkfj/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4aqrrrrqv8qJex/CVVf9X3fx2qdz1VVXXfVf5cde77V41WvO8Hwg2+b52BsnXvFnfoH9ceSqq6666r/CiXsfwlVX/V938dqnc9VVV131X2G76/jTt3ozdrrK84FsmxfgM//yb/iOJz6Zq6666qr/CifufQhXXfV/3cVrn85VV1111X+Fj33xx/JxL/FYXgBk27wQn/mXf8Mv33Endx4ecdVVV131n+nEvQ/hqqv+r7t47dO56qqrrvrPdOPmBm9804187su+FC8Esm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXXXVVVddddVVV/1/gmybF+Kz/upv+OXb7+KOw0Ouuuqqq/4znbj3IVx11f91F699OlddddVV/5lu2tzkjW++gc95mZfihUC2zQvwmX/5N3zHE58MwI2LBVddddVV/5mOnbuZq676v+7S6du56qqrrvrPdOdyCcD7PeoRfO7LvhQvALJtno+9ceKVfuYX2RsHAF7x5Emuuuqqq/4zbV+8nquu+r9u/8TdXHXVVVf9Z/rTCxcA2Ol6/uSt3pSdrvJ8INvm+fjD+87yDr/xO9zvFU+e5KqrrrrqP9P2xeu56qr/6/ZP3M1VV1111X+mP71wgfv92Ou9Fq96zRmeD2TbPB9/eN9Z3uE3fof73bhYcNVVV131n+nYuZu56qr/6y6dvp2rrrrqqv9Mdy6X3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO1x11VVX/Vc5ce9DuOqq/+suXvt0rrrqqqv+q/zY670Wr3rNGZ4PZNs8H3vjxCv+zC+wP45cddVVV/1XOHHvQ7jqqv/rLl77dK666qqr/itsdx1/+lZvxk5XeT6QbfMCfOZf/g3f8cQnc9VVV131X+HEvQ/hqqv+r7t47dO56qqrrvqv8LEv/lg+7iUeywuAbJsX4jP/8m/45Tvu5M7DI6666qqr/jOduPchXHXV/3UXr306V1111VX/mW7c3OCNb7qRz33Zl+KFQLbNVVddddVVV1111VVXXfX/CbJtrrrqqquuuuqqq6666qr/T5Bt80J81l/9Db98+13ccXjIVVddddV/phP3PoSrrvq/7uK1T+eqq6666j/TTZubvPHNN/A5L/NSvBDItnkBPvMv/4bveOKTAbhxseCqq6666j/TsXM3c9VV/9ddOn07V1111VX/me5cLgF4v0c9gs992ZfiBUC2zfOxN0680s/8InvjAMArnjzJVVddddV/pu2L13PVVf/X7Z+4m6uuuuqq/0x/euECADtdz5+81Zuy01WeD2TbPB9/eN9Z3uE3fof7veLJk1x11VVX/Wfavng9V131f93+ibu56qqrrvrP9KcXLnC/H3u91+JVrznD84Fsm+fjD+87yzv8xu9wvxsXC6666qqr/jMdO3czV131f92l07dz1VVXXfWf6c7lkvv92Ou9Fq96zRmeD2TbPB9/eN9Z3uE3foerrrrqqv8qJ+59CFdd9X/dxWufzlVXXXXVf5Ufe73X4lWvOcPzgWyb52NvnHjFn/kF9seRq6666qr/CifufQhXXfV/3cVrn85VV1111X+F7a7jT9/qzdjpKs8Hsm1egM/8y7/hO574ZK666qqr/iucuPchXHXV/3UXr306V1111VX/FT72xR/Lx73EY3kBkG3zQnzmX/4Nv3zHndx5eMRVV1111X+mE/c+hKuu+r/u4rVP56qrrrrqP9ONmxu88U038rkv+1K8EMi2ueqqq6666qqrrrrqqqv+P0G2zVVXXXXVVVddddVVV131/wmybV6Iz/qrv+GXb7+LOw4Pueqqq676z3Ti3odw1VX/11289ulcddVVV/1numlzkze++QY+52VeihcC2TYvwGf+5d/wHU98MgA3LhZcddVVV/1nOnbuZq666v+6S6dv56qrrrrqP9OdyyUA7/eoR/C5L/tSvADItnk+9saJV/qZX2RvHAB4xZMnueqqq676z7R98Xquuur/uv0Td3PVVVdd9Z/pTy9cAGCn6/mTt3pTdrrK84Fsm+fjD+87yzv8xu9wv1c8eZKrrrrqqv9M2xev56qr/q/bP3E3V1111VX/mf70wgXu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gd7nfjYsFVV1111X+mY+du5qqr/q+7dPp2rrrqqqv+M925XHK/H3u91+JVrznD84Fsm+fjD+87yzv8xu9w1VVXXfVf5cS9D+Gqq/6vu3jt07nqqquu+q/yY6/3WrzqNWd4PpBt83zsjROv+DO/wP44ctVVV131X+HEvQ/hqqv+r7t47dO56qqrrvqvsN11/OlbvRk7XeX5QLbNC/CZf/k3fMcTn8xVV1111X+FE/c+hKuu+r/u4rVP56qrrrrqv8LHvvhj+biXeCwvALJtXojP/Mu/4ZfvuJM7D4+46qqrrvrPdOLeh3DVVf/XXbz26Vx11VVX/We6cXODN77pRj73ZV+KFwLZNlddddVVV1111VVXXXXV/yfItrnqqquuuuqqq6666qqr/j9Bts0L8Vl/9Tf88u13ccfhIVddddVV/5lO3PsQrrrq/7qL1z6dq6666qr/TDdtbvLGN9/A57zMS/FCINvmBfjMv/wbvuOJTwbgxsWCq6666qr/TMfO3cxVV/1fd+n07Vx11VVX/We6c7kE4P0e9Qg+92VfihcA2TbPx9448Uo/84vsjQMAr3jyJFddddVV/5m2L17PVVf9X7d/4m6uuuqqq/4z/emFCwDsdD1/8lZvyk5XeT6QbfN8/OF9Z3mH3/gd7veKJ09y1VVXXfWfafvi9Vx11f91+yfu5qqrrrrqP9OfXrjA/X7s9V6LV73mDM8Hsm2ejz+87yzv8Bu/w/1uXCy46qqrrvrPdOzczVx11f91l07fzlVXXXXVf6Y7l0vu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gdrrrqqqv+q5y49yFcddX/dRevfTpXXXXVVf9Vfuz1XotXveYMzweybZ6PvXHiFX/mF9gfR6666qqr/iucuPchXHXV/3UXr306V1111VX/Fba7jj99qzdjp6s8H8i2eQE+8y//hu944pO56qqrrvqvcOLeh3DVVf/XXbz26Vx11VVX/Vf42Bd/LB/3Eo/lBUC2zQvxmX/5N/zyHXdy5+ERV1111VX/mU7c+xCuuur/uovXPp2rrrrqqv9MN25u8MY33cjnvuxL8UIg2+aqq6666qqrrrrqqquu+v8E2TZXXXXVVVddddVVV1111f8nyLZ5IT7rr/6GX779Lv7kLd+E5/aH953lHX7jd7jqf7bhcf/Auc//bGaPeTFOfcZnc9VV/1OduPchXHXV/3UXr306V1111VX/mW7a3OSNb76Bz3mZl+KFQLbNC/CZf/k3fMcTnwzAne/y9jy3P77vLB/5h3/KVf+z7f/93/GEz/gUdl78JXjU530RV131P9Wxczdz1VX/1106fTtXXXXVVf+Z7lwuAXi/Rz2Cz33Zl+IFQLbN87E3TrzSz/wie+MAwJ3v8vY8tz+57xxf/Fd/y1X/s53727/h9z/xYzn9ki/Fq3/pV3LVVf9TbV+8nquu+r9u/8TdXHXVVVf9Z/rTCxcA2Ol6/uSt3pSdrvJ8INvm+fjD+87yDr/xO9zvB1/7NXhudx4e8mNPvZWr/mc797d/w+9/4sdy+iVfilf/0q/kqqv+p9q+eD1XXfV/3f6Ju7nqqquu+s/0pxcucL8fe73X4lWvOcPzgWyb5+MP7zvLO/zG73C/u46WPLdXOn2adHLV/2z7f/93POEzPoWdF38JHvV5X8RVV/1PdezczVx11f91l07fzlVXXXXVf6Y7l0vu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gd7vcDr/3qPLdnHBzyqX/+V1z1P9vwuH/g3Od/NrPHvBinPuOzueqq/6lO3PsQrrrq/7qL1z6dq6666qr/Kj/2eq/Fq15zhucD2TbPx9448Yo/8wvsjyMAd77L2/Pc/vC+s7zDb/wOV/3PNjzuHzj3+Z/N7DEvxqnP+Gyuuup/qhP3PoSrrvq/7uK1T+eqq6666r/Cdtfxp2/1Zux0lecD2TYvwGf+5d/wHU98MgB3vsvb89z+8L6zvMNv/A5X/c82PO4fOPf5n83sMS/Gqc/4bK666n+qE/c+hKuu+r/u4rVP56qrrrrqv8LHvvhj+biXeCwvALJtXojP/Mu/4ZfvuJM/fcs35bn94X1neYff+B2u+p9teNw/cO7zP5vZY16MU5/x2Vx11f9UJ+59CFdd9X/dxWufzlVXXXXVf6YbNzd445tu5HNf9qV4IZBtc9VVV1111VVXXXXVVVf9f4Jsm6uuuuqqq6666qqrrrrq/xNk27wQn/VXf8Mv334XdxwectVVV131n+nEvQ/hqqv+r7t47dO56qqrrvrPdNPmJm988w18zsu8FC8Esm1egM/8y7/hO574ZABuXCy46qqrrvrPdOzczVx11f91l07fzlVXXXXVf6Y7l0sA3u9Rj+BzX/aleAGQbfN87I0Tr/Qzv8jeOADwiidPctVVV131n2n74vVcddX/dfsn7uaqq6666j/Tn164AMBO1/Mnb/Wm7HSV5wPZNs/HH953lnf4jd/hfq948iRXXXXVVf+Zti9ez1VX/V+3f+Jurrrqqqv+M/3phQvc78de77V41WvO8Hwg2+b5+MP7zvIOv/E73O/GxYKrrrrqqv9Mx87dzFVX/V936fTtXHXVVVf9Z7pzueR+P/Z6r8WrXnOG5wPZNs/HH953lnf4jd/hqquuuuq/yol7H8JVV/1fd/Hap3PVVVdd9V/lx17vtXjVa87wfCDb5vnYGyde8Wd+gf1x5Kqrrrrqv8KJex/CVVf9X3fx2qdz1VVXXfVfYbvr+NO3ejN2usrzgWybF+Az//Jv+I4nPpmrrrrqqv8KJ+59CFdd9X/dxWufzlVXXXXVf4WPffHH8nEv8VheAGTbvBCf+Zd/wy/fcSd3Hh5x1VVXXfWf6cS9D+Gqq/6vu3jt07nqqquu+s904+YGb3zTjXzuy74ULwSyba666qqrrrrqqquuuuqq/0+QbXPVVVddddVVV1111VVX/X+CbJsX4rP+6m/45dvv4o7DQ6666qqr/rPc9a7vAMCLfdWf83/Rrd/wQRw+5S94yId9CxsPfzn+v/uHj3l5AF7sq/6c/28uXvt0rrrqqqv+M920uckb33wDn/MyL8ULgWybF+Az//Jv+I4nPhmAGxcLrrrqqqv+s/zZ27w5AK//bX/D/0V/8eXvx8Un/jkv9/HfzolHvQL/3/36B7wUAK//bX/D/zeXTt/OVVddddV/pjuXSwDe71GP4HNf9qV4AZBt83zsjROv9DO/yN44APCKJ09y1VVXXfWf5aff+PUAeI8fehz/F/3q570X9z7uz3jDz/hurn3sK/L/3fe9y2MBeI8fehz/3+yfuJurrrrqqv9Mf3rhAgA7Xc+fvNWbstNVng9k2zwff3jfWd7hN36H+73iyZNcddVVV/1n+ek3fj0A3uOHHsf/Rb/6ee/FvY/7M97wM76bax/7ivx/933v8lgA3uOHHsf/N/sn7uaqq6666j/Tn164wP1+7PVei1e95gzPB7Jtno8/vO8s7/Abv8P9blwsuOqqq676z/Jnb/PmALz+t/0N/xf9xZe/Hxef+Oe83Md/Oyce9Qr8f/frH/BSALz+t/0N/99cOn07V1111VX/me5cLrnfj73ea/Gq15zh+UC2zfPxh/ed5R1+43e46qqrrvqvcNe7vgMAL/ZVf87/Rbd+wwdx+JS/4CEf9i1sPPzl+P/uHz7m5QF4sa/6c/6/uXjt07nqqquu+q/yY6/3WrzqNWd4PpBt83zsjROv+DO/wP44ctVVV131n+2ud30HAF7sq/6c/4tu/YYP4vApf8FDPuxb2Hj4y/H/3T98zMsD8GJf9ef8f3Px2qdz1VVXXfVfYbvr+NO3ejN2usrzgWybF+Az//Jv+I4nPpmrrrrqqv9sd73rOwDwYl/15/xfdOs3fBCHT/kLHvJh38LGw1+O/+/+4WNeHoAX+6o/5/+bi9c+nauuuuqq/wof++KP5eNe4rG8AMi2eSE+8y//hl++407uPDziqquuuuo/y13v+g4AvNhX/Tn/F936DR/E4VP+god82Lew8fCX4/+7f/iYlwfgxb7qz/n/5uK1T+eqq6666j/TjZsbvPFNN/K5L/tSvBDItrnqqquuuuqqq6666qqr/j9Bts1VV1111VVXXXXVVVdd9f8Jsm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXXXVVVddddVVV/1/gmybq6666qqrrrrqqquuuur/E2TbXHXVVVddddVVV1111VX/nyDb5qqrrrrqqquuuuqqq676/wTZNlddddVVV1111VVXXXXV/yfItrnqqquuuuqqq6666qqr/j9Bts1VV1111VVXXXXVVVdd9f8Jsm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXfXfaHd3l4/5mI8B4Ku+6qs4fvw4V/3b3XrrrXzP93wP93urt3orXvqlX5qr/nX++q//mt/5nd/h4sWLfPZnfzb3293d5XM+53O49dZb+azP+ixe+qVfmquuuuqq/4WQbXPVVVdd9d/oIQ95CB/1UR/F7u4uv/M7v8Nv/dZvcdW/3W//9m/z3u/93rz3e783AO/93u/Ngx/8YK7613nt135t7vfbv/3b3O+93/u9AXjpl35pvuZrvoanP/3pXHXVVVf9L4Rsm6uuuuqq/yZ//dd/zXu/93vz13/91wBIwjZX/dv99m//Nu/zPu/De7/3e/PgBz+Y93qv9+Kqf5vf/u3f5rM/+7P57d/+be4niYsXL3L8+HFe+qVfmq/+6q/mtV/7tbnqqquu+l8G2TZXXXXVVf9Nfvu3f5vP/uzP5rd/+7cBeO3Xfm0++7M/m9d+7dfmqn+bv/7rv+ajP/qjee3Xfm1++qd/mrd+67fmsz/7s7nqX++3f/u3+ezP/mx++7d/m/tJwjYAr/3ar81nf/Zn89qv/dpcddVVV/0vg2ybq6666qr/Jr/927/Nx3zMx/BXf/VXAJw4cYKnP/3pHD9+nKv+/X77t3+bz/7sz+a3f/u3uepf77d/+7f57M/+bH77t3+b+0nCNgCv/dqvzWd/9mfz2q/92lx11VVX/S+DbJurrrrqqv9GD37wg/nqr/5qnvGMZ/Bd3/Vd/PVf/zVX/dt993d/NydOnOClXuql+J7v+R7+6q/+ip/+6Z/mqn+dv/7rv+Z3fud3+O7v/m6++qu/mpd6qZfi+PHjvPZrvzav8zqvw2u91mvx1m/91tx6660cP36cq6666qr/ZZBtc9VVV1313+iv//qvee/3fm+OHz/OV3/1V/PSL/3SXPVv99M//dN89Vd/Nbu7uzz4wQ/mu7/7uzl+/DhX/eu89Vu/Nbu7u9zvq7/6q3npl35p/vqv/5qP/uiPZnd3l+/+7u/mpV/6pbnqqquu+l8I2TZXXXXVVVddddVVV1111f8nyLa56qqrrvpv9td//dcAvPRLvzT/Vr/zO7/Dgx70IB784Afzb/U7v/M7POhBD+LBD34w/1a/8zu/w4Me9CAe/OAH82+xu7vL3/zN3/BSL/VSHD9+nH+L3d1d/uZv/oaXeqmX4vjx4/xb7O7u8jd/8ze81Eu9FMePH+ffYnd3l7/5m7/hpV7qpTh+/Dj/FrfeeivPeMYzeK3Xei3+rW699Vae8Yxn8Fqv9Vr8W91666084xnP4LVe67W46qqrrvo/ANk2V1111VX/zV77tV8bgN/+7d/m30oSn/VZn8Vnf/Zn828lic/6rM/isz/7s/m3ksRnfdZn8dmf/dn8W/z2b/82r/M6r8Nv/dZv8dqv/dr8W/z2b/82r/M6r8Nv/dZv8dqv/dr8W/z2b/82r/M6r8Nv/dZv8dqv/dr8W/z2b/82r/M6r8Nv/dZv8dqv/dr8W3z2Z382n/M5n4Nt/q0++7M/m8/5nM/BNv9Wn/3Zn83nfM7nYJurrrrqqv8DkG3zAuzt7fGjP/qj/NEf/RFPe9rTuOqqq676z/LXf/3XALz0S780/1a//du/zYMf/GAe/OAH82/127/92zz4wQ/mwQ9+MP9Wv/3bv82DH/xgHvzgB/Nvsbu7y1//9V/z0i/90hw/fpx/i93dXf76r/+al37pl+b48eP8W+zu7vLXf/3XvPRLvzTHjx/n32J3d5e//uu/5qVf+qU5fvw4/xa33nort956K6/92q/Nv9Wtt97Krbfeymu/9mvzb3Xrrbdy66238tqv/dpcddVVV/1P9tCHPpRXeZVX4R3f8R3Z2dnhBUC2zfOxt7fH27/92/PHf/3X5PETTLXjqquuuuo/y/CkJwBi4yEvxb/V4VP+gv7kDXQnr+ff6vApf0F/8ga6k9fzb3X4lL+gP3kD3cnr+bfI5QHLO5/I4sZHEYst/i1yecDyzieyuPFRxGKLf4tcHrC884ksbnwUsdji3yKXByzvfCKLGx9FLLb4txgv3M1w4S42H/5y/FuNF+5muHAXmw9/Of6txgt3M1y4i82Hvxz/Vqs7n4Qj6W55MFddddVV/1nqNBK7F3nll35pfvzHf5ydnR2eD2TbPB9f8RVfwed8yZfQv+8HMnupl+bGxYKrrrrqqv8sf/Xu74wsXutrfp9/q1//gJfioW/xwTz0LT+Ef6tf/4CX4qFv8cE89C0/hH+rX/+Al+Khb/HBPPQtP4R/i4tP/DP+4svfn5f7+G/nxKNegX+Li0/8M/7iy9+fl/v4b+fEo16Bf4uLT/wz/uLL35+X+/hv58SjXoF/i4tP/DP+4svfn5f7+G/nxKNegX+Lp/3sN/G0n/tmXv/b/oZ/q6f97DfxtJ/7Zl7/2/6Gf6un/ew38bSf+2Ze/9v+hn+rv/jy96N1ax71eV/EVVddddV/ljuXS9Z/89cM3/mtfNYnfRIf93Efx/OBbJvn41GPfjT33nATmx/wwQC84smTXHXVVVf9Z/mdj/xQYqq86Rf8GP9W9z7uT9k6cyObZ27k3+rex/0pW2duZPPMjfxb3fu4P2XrzI1snrmRf4vhcI+Lz3gCJx70aPrNHf4thsM9Lj7jCZx40KPpN3f4txgO97j4jCdw4kGPpt/c4d9iONzj4jOewIkHPZp+c4d/i8Ozd3Jw9k6ufewr8m91ePZODs7eybWPfUX+rQ7P3snB2Tu59rGvyL/Vr37ee9HqwKt/6Vdy1VVXXfWf5U8vXADg8Nu+mWvvuoMnPuEJPB/Itnk+zlx3HcPrvxFbb/YWALziyZNcddVVV/1n+f1P/FjK1POGn/E9XHXV/1W/+nnvRasDr/6lX8lVV1111X+WP71wAYCDX/g5+l//Fc7ecw/PB7Jtno8z113H8PpvxNabvQUANy4WXHXVVVf9Z3niZ3wKZZzxch//HVx11f9Vf/Hl70fr1jzq876Iq6666qr/LHculwAc/MLP0f/6r3D2nnt4PpBt83ycue46htd/I7be7C246qqrrvrPdv7zPps6znnwh30LV131f9Wt3/BBTN2KU5/x2Vx11VVX/Wc7+IWfo//1X+HsPffwfCDb5vl45KMexX033szmB3wwV1111VX/2c5/3mdTxzkP/rBv4aqr/q+69Rs+iKlbceozPpurrrrqqv9sh9/2zVxz5+086YlP5PlAts3z8RVf8RV8zpd8Cf37fiCzl3pprrrqqqv+M53/vM+mjnMe/GHfwlVX/V916zd8EFO34tRnfDZXXXXVVf+Z1n/z14zf+a185id9Eh/3cR/H84Fsm+djb2+Pt3u7t+NP/uZvyBMnmWrlqquuuuo/y/iMW1EG8xsfyVVX/V+1uvNJOJLuQQ/mqquuuuo/S50m4uIFXumlXoqf+ImfYGdnh+cD2TYvwN7eHj/yIz/CH/3RH/G0pz2Nq6666qqrrrrqqquuuup/toc+9KG8yqu8Cu/0Tu/Ezs4OLwCyba666qqrrrrqqquuuuqq/0+QbXPVVVddddVVV1111VVX/X+CbJsX4rP+6m/45dvv4o7DQ6666qqr/jOduPchPLcPfeXr+dBXvoEXxZ/dsc/7/PiTuOqq/6n+4WNeHoAbfvDHuOqqq676z3LT5iZvfPMNfM7LvBQvBLJtXoDP/Mu/4Tue+GQAblwsuOqqq676z3Ts3M08t3d6qTO8/ytcz4viz+/Y51N++elcddX/VL/+AS8FwCv81M9z1VVXXfWf5c7lEoD3e9Qj+NyXfSleAGTbPB9748Qr/cwvsjcOALziyZNcddVVV/1n2r54Pc/t9R9xgvd8uet4Ufzlnft89e/dwVVX/U/1fe/yWADe+pd/g6uuuuqq/yx/euECADtdz5+81Zuy01WeD2TbPB9/eN9Z3uE3fof7veLJk1x11VVX/Wfavng9z+1lb9riFW7e4UVxz/7Iz/z9Wa666n+q73uXxwLw1r/8G1x11VVX/Wf50wsXuN+Pvd5r8arXnOH5QLbN8/GH953lHX7jd7jfjYsFV1111VX/mY6du5nndmqz49efusuL4uVv3KI1c9VV/1P9+ge8FACv8FM/z1VXXXXVf5Y7l0vu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gdrrrqqqv+q5y49yE8t3d9qTO81kOP8aK4/dLA5//mbVx11f9U//AxLw/ADT/4Y1x11VVX/Vf4sdd7LV71mjM8H8i2eT72xolX/JlfYH8cueqqq676r3Di3ofw3D70la/nQ1/5Bl4Uf3bHPu/z40/iqqv+p/qHj3l5AG74wR/jqquuuuo/23bX8adv9WbsdJXnA9k2L8Bn/uXf8B1PfDJXXXXVVf8VTtz7EJ7bh77y9XzoK9/Ai+LP7tjnfX78SVx11f9U//AxLw/ADT/4Y1x11VVX/Wf72Bd/LB/3Eo/lBUC2zQvxmX/5N/zyHXdy5+ERV1111VX/mU7c+xCe24e+8vV86CvfwIviz+7Y531+/ElcddX/VP/wMS8PwA0/+GNcddVVV/1nuXFzgze+6UY+92VfihcC2TZXXXXVVVddddVVV1111f8nyLa56qqrrrrqqquuuuqqq/4/QbbNC/FZf/U3/PLtd3HH4SFXXXXV/x+vfM0ZfuL1Xovn9hV/9zi+8u8fx3+GE/c+hKuu+r/u4rVP56r/OfZ/4kfZ/4kfY/vt3oHtt3tHrrrqv9ud7/L2PLc/vO8s7/Abv8OL6qbNTd745hv4nJd5KV4IZNu8AJ/5l3/DdzzxyQDcuFhw1VVX/f/x4idP8J2v+ao8t6/7hyfwfU9+Kv8Zjp27mauu+r/u0unbuep/jrt++Ae580d+kBvf6V254Z3flauu+u/2p2/9Zjy3P77vLB/5h3/Ki+rO5RKA93vUI/jcl30pXgBk2zwfe+PEK/3ML7I3DgC84smTXHXVVf9/PHh7i6961VfkuX3r457EL91+B/8Zti9ez1VX/V+3f+Jurvqf4wnf/z084fu/l0e/+3vy6Hd/L6666r/bT73R6/Lc/uS+c3zxX/0tL6o/vXABgJ2u50/e6k3Z6SrPB7Jtno8/vO8s7/Abv8P9XvHkSa666qr/P67fWPBOD38oz+1P7r2PP7r3LP8Zti9ez1VX/V+3f+Jurvqf4wnf/z084fu/l0e/+3vy6Hd/L6666r/bR7/ki/Hc7jw85Meeeisvqj+9cIH7/djrvRaves0Zng9k2zwff3jfWd7hN36H+924WHDVVVf9/3F6PueX7ryL5/Y2D7qZOw4O+c9w7NzNXHXV/3WXTt/OVf9z3PXDP8idP/KD3PhO78oN7/yuXHXVf7c/O3+B5/ZKp0+TTl5Udy6X3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO1x11VX/P73c6VN87Is/huf2a3fezXc/+an8Zzhx70O46qr/6y5e+3Su+p9j/yd+lP2f+DG23+4d2H67d+Sqq/67/cBrvzrP7RkHh3zqn/8V/xY/9nqvxatec4bnA9k2z8feOPGKP/ML7I8jV1111f8/r3zNGX7i9V6L5/YVf/c4vvLvH8d/hhP3PoSrrvq/7uK1T+eq/zn2f+JH2f+JH2P77d6B7bd7R6666r/bne/y9jy3P7zvLO/wG7/Dv9Z21/Gnb/Vm7HSV5wPZNi/AZ/7l3/AdT3wyV1111f8/r3zNGX7i9V6L5/YVf/c4vvLvH8d/hhP3PoSrrvq/7uK1T+eq/zn2f+JH2f+JH2P77d6B7bd7R6666r/bne/y9jy3P7zvLO/wG7/Dv9bHvvhj+biXeCwvALJtXojP/Mu/4ZfvuJM7D4+46qqr/v945WvO8BOv91o8t6/4u8fxlX//OP4znLj3IVx11f91F699Olf9z7H/Ez/K/k/8GNtv9w5sv907ctVV/93ufJe357n94X1neYff+B1eVDdubvDGN93I577sS/FCINvmqquuuuqqq6666qqrrvr/BNk2V1111VVXXXXVVVddddX/J8i2eSE+66/+hl++/S7uODzkqquuuuo/04l7H8JVV/1fd/Hap3PVVVdd9Z/pps1N3vjmG/icl3kpXghk27wAn/mXf8N3PPHJANy4WHDVVVdd9Z/p2Lmbueqq/+sunb6dq6666qr/THculwC836Mewee+7EvxAiDb5vnYGyde6Wd+kb1xAOAVT57kqquuuuo/0/bF67nqqv/r9k/czVVXXXXVf6Y/vXABgJ2u50/e6k3Z6SrPB7Jtno8/vO8s7/Abv8P9XvHkSa666qqr/jNtX7yeq676v27/xN1cddVVV/1n+tMLF7jfj73ea/Gq15zh+UC2zfPxh/ed5R1+43e4342LBVddddVV/5mOnbuZq676v+7S6du56qqrrvrPdOdyyf1+7PVei1e95gzPB7Jtno8/vO8s7/Abv8NVV1111X+VE/c+hKuu+r/u4rVP56qrrrrqv8qPvd5r8arXnOH5QLbN87E3Trziz/wC++PIVVddddV/hRP3PoSrrvq/7uK1T+eqq6666r/Cdtfxp2/1Zux0lecD2TYvwGf+5d/wHU98MlddddVV/xVO3PsQrrrq/7qL1z6dq6666qr/Ch/74o/l417isbwAyLZ5IT7zL/+GX77jTu48POKqq6666j/TiXsfwlVX/V938dqnc9VVV131n+nGzQ3e+KYb+dyXfSleCGTbXHXVVVddddVVV1111VX/nyDb5qqrrrrqqquuuuqqq676/wTZNi/EZ/3V3/DLt9/FHYeHXHXVVVf9Zzpx70O46j/Op75sx03s8ty+9c6T/PW9a67673Hx2qdz1VVXXfWf6abNTd745hv4nJd5KV4IZNu8AJ/5l3/DdzzxyQDcuFhw1VVXXfWf6di5m7nqP86HPrZwvfZ5bt979zGefGHkqv8el07fzlVXXXXVf6Y7l0sA3u9Rj+BzX/aleAGQbfN87I0Tr/Qzv8jeOADwiidPctVVV131n2n74vVc9R/nnR8K15clz+1H797izoPGVf899k/czVVXXXXVf6Y/vXABgJ2u50/e6k3Z6SrPB7Jtno8/vO8s7/Abv8P9XvHkSa666qqr/jNtX7yeq/7jvPND4fqy5Ln96N1b3HnQuOq/x/6Ju7nqqquu+s/0pxcucL8fe73X4lWvOcPzgWyb5+MP7zvLO/zG73C/GxcLrrrqqqv+Mx07dzNX/cf50McWro99ntv33nWMJ18Yueq/x6XTt3PVVVdd9Z/pzuWS+/3Y670Wr3rNGZ4PZNs8H39431ne4Td+h6uuuuqq/yon7n0IV/3H+dSX7biJXZ7bt955kr++d81V/z0uXvt0rrrqqqv+q/zY670Wr3rNGZ4PZNs8H3vjxCv+zC+wP45cddVVV/1XOHHvQ7jqP86nvmzHTezy3L71zpP89b1rrvrvcfHap3PVVVdd9V9hu+v407d6M3a6yvOBbJsX4DP/8m/4jic+mauuuuqq/won7n0IV/3H+dSX7biJXZ7bt955kr++d81V/z0uXvt0rrrqqqv+K3zsiz+Wj3uJx/ICINvmhfjMv/wbfvmOO7nz8Iirrrrqqv9MJ+59CFf9x/nUl+24iV2e27feeZK/vnfNVf89Ll77dK666qqr/jPduLnBG990I5/7si/FC4Fsm6uuuuqqq6666qqrrrrq/xNk21x11VVXXXXVVVddddVV/58g2+aF+Ky/+ht++fa7uOPwkKuuuuqq/0wn7n0IV131f93Fa5/OVVddddV/pps2N3njm2/gc17mpXghkG3zAnzmX/4N3/HEJwNw42LBVVddddV/pmPnbuaqq/6vu3T6dq666qqr/jPduVwC8H6PegSf+7IvxQuAbJvnY2+ceKWf+UX2xgGAVzx5kquuuuqq/0zbF6/nqqv+r9s/cTdXXXXVVf+Z/vTCBQB2up4/eas3ZaerPB/Itnk+/vC+s7zDb/wO93vFkye56qqrrvrPtH3xeq666v+6/RN3c9VVV131n+lPL1zgfj/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4X43LhZcddVVV/1nOnbuZq666v+6S6dv56qrrrrqP9OdyyX3+7HXey1e9ZozPB/Itnk+/vC+s7zDb/wOV1111VX/VU7c+xCuuur/uovXPp2rrrrqqv8qP/Z6r8WrXnOG5wPZNs/H3jjxij/zC+yPI1ddddVV/xVO3PsQrrrq/7qL1z6dq6666qr/Cttdx5++1Zux01WeD2TbvACf+Zd/w3c88clcddVVV/1XOHHvQ7jqqv/rLl77dK666qqr/it87Is/lo97icfyAiDb5oX4zL/8G375jju58/CIq6666qr/TCfufQhXXfV/3cVrn85VV1111X+mGzc3eOObbuRzX/aleCGQbXPVVVddddVVV1111VVX/X+CbJurrrrqqquuuuqqq6666v8TZNu8EJ/1V3/DL99+F3ccHnLVVVdd9Z/pxL0P4aqr/q+7eO3Tueqqq676z3TT5iZvfPMNfM7LvBQvBLJtXoDP/Mu/4Tue+GQAblwsuOqqq676z3Ts3M1cddX/dZdO385VV1111X+mO5dLAN7vUY/gc1/2pXgBkG3zfOyNE6/0M7/I3jgA8IonT3LVVVdd9Z9p++L1XHXV/3X7J+7mqquuuuo/059euADATtfzJ2/1pux0lecD2TbPxx/ed5Z3+I3f4X6vePIkV1111VX/mbYvXs9VV/1ft3/ibq666qqr/jP96YUL3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO9zvxsWCq6666qr/TMfO3cxVV/1fd+n07Vx11VVX/We6c7nkfj/2eq/Fq15zhucD2TbPxx/ed5Z3+I3f4aqrrrrqv8qJex/CVVf9X3fx2qdz1VVXXfVf5cde77V41WvO8Hwg2+b52BsnXvFnfoH9ceSqq6666r/CiXsfwlVX/V938dqnc9VVV131X2G76/jTt3ozdrrK84FsmxfgM//yb/iOJz6Zq6666qr/CifufQhXXfV/3cVrn85VV1111X+Fj33xx/JxL/FYXgBk27wQn/mXf8Mv33Endx4ecdVVV131n+nEvQ/hqqv+r7t47dO56qqrrvrPdOPmBm9804187su+FC8Esm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXXXVVVddddVVV/1/gmybF+Kz/upv+OXb7+KOw0Ouuuqqq/4znbj3IVx11f91F699OlddddVV/5lu2tzkjW++gc95mZfihUC2zQvwmX/5N3zHE58MwI2LBVddddVV/5mOnbuZq676v+7S6du56qqrrvrPdOdyCcD7PeoRfO7LvhQvALJtno+9ceKVfuYX2RsHAF7x5Emuuuqqq/4zbV+8nquu+r9u/8TdXHXVVVf9Z/rTCxcA2Ol6/uSt3pSdrvJ8INvm+fjD+87yDr/xO9zvFU+e5KqrrrrqP9P2xeu56qr/6/ZP3M1VV1111X+mP71wgfv92Ou9Fq96zRmeD2TbPB9/eN9Z3uE3fof73bhYcNVVV131n+nYuZu56qr/6y6dvp2rrrrqqv9Mdy6X3O/HXu+1eNVrzvB8INvm+fjD+87yDr/xO1x11VVX/Vc5ce9DuOqq/+suXvt0rrrqqqv+q/zY670Wr3rNGZ4PZNs8H3vjxCv+zC+wP45cddVVV/1XOHHvQ7jqqv/rLl77dK666qqr/itsdx1/+lZvxk5XeT6QbfMCfOZf/g3f8cQnc9VVV131X+HEvQ/hqqv+r7t47dO56qqrrvqv8LEv/lg+7iUeywuAbJsX4jP/8m/45Tvu5M7DI6666qqr/jOduPchXHXV/3UXr306V1111VX/mW7c3OCNb7qRz33Zl+KFQLbNVVddddVVV1111VVXXfX/CbJtrrrqqquuuuqqq6666qr/T5Bt80J81l/9Db98+13ccXjIVVddddV/phP3PoSrrvq/7uK1T+eqq6666j/TTZubvPHNN/A5L/NSvBDItnkBPvMv/4bveOKTAbhxseCqq6666j/TsXM3c9VV/9ddOn07V1111VX/me5cLgF4v0c9gs992ZfiBUC2zfOxN0680s/8InvjAMArnjzJVVddddV/pu2L13PVVf/X7Z+4m6uuuuqq/0x/euECADtdz5+81Zuy01WeD2TbPB9/eN9Z3uE3fof7veLJk1x11VVX/Wfavng9V131f93+ibu56qqrrvrP9KcXLnC/H3u91+JVrznD84Fsm+fjD+87yzv8xu9wvxsXC6666qqr/jMdO3czV131f92l07dz1VVXXfWf6c7lkvv92Ou9Fq96zRmeD2TbPB9/eN9Z3uE3foerrrrqqv8qJ+59CFdd9X/dxWufzlVXXXXVf5Ufe73X4lWvOcPzgWyb52NvnHjFn/kF9seRq6666qr/CifufQhXXfV/3cVrn85VV1111X+F7a7jT9/qzdjpKs8Hsm1egM/8y7/hO574ZK666qqr/iucuPchXHXV/3UXr306V1111VX/FT72xR/Lx73EY3kBkG3zQnzmX/4Nv3zHndx5eMRVV1111X+mE/c+hKuu+r/u4rVP56qrrrrqP9ONmxu88U038rkv+1K8EMi2ueqqq6666qqrrrrqqqv+P0G2zVVXXXXVVVddddVVV131/wmybV6Iz/qrv+GXb7+LOw4Pueqqq676z3Ti3odw1VX/11289ulcddVVV/1numlzkze++QY+52VeihcC2TYvwGf+5d/wHU98MgA3LhZcddVVV/1nOnbuZq666v+6S6dv56qrrrrqP9OdyyUA7/eoR/C5L/tSvADItnk+9saJV/qZX2RvHAB4xZMnueqqq676z7R98Xquuur/uv0Td3PVVVdd9Z/pTy9cAGCn6/mTt3pTdrrK84Fsm+fjD+87yzv8xu9wv1c8eZKrrrrqqv9M2xev56qr/q/bP3E3V1111VX/mf70wgXu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gd7nfjYsFVV1111X+mY+du5qqr/q+7dPp2rrrqqqv+M925XHK/H3u91+JVrznD84Fsm+fjD+87yzv8xu9w1VVXXfVf5cS9D+Gqq/6vu3jt07nqqquu+q/yY6/3WrzqNWd4PpBt83zsjROv+DO/wP44ctVVV131X+HEvQ/hqqv+r7t47dO56qqrrvqvsN11/OlbvRk7XeX5QLbNC/CZf/k3fMcTn8xVV1111X+FE/c+hKuu+r/u4rVP56qrrrrqv8LHvvhj+biXeCwvALJtXojP/Mu/4ZfvuJM7D4+46qqrrvrPdOLeh3DVVf/XXbz26Vx11VVX/We6cXODN77pRj73ZV+KFwLZNlddddVVV1111VVXXXXV/yfItrnqqquuuuqqq6666qqr/j9Bts0L8Vl/9Tf88u13ccfhIVddddVV/5lO3PsQrrrq/7qL1z6dq6666qr/TDdtbvLGN9/A57zMS/FCINvmBfjMv/wbvuOJTwbgxsWCq6666qr/TMfO3cxVV/1fd+n07Vx11VVX/We6c7kE4P0e9Qg+92VfihcA2TbPx9448Uo/84vsjQMAr3jyJFddddVV/5m2L17PVVf9X7d/4m6uuuqqq/4z/emFCwDsdD1/8lZvyk5XeT6QbfN8/OF9Z3mH3/gd7veKJ09y1VVXXfWfafvi9Vx11f91+yfu5qqrrrrqP9OfXrjA/X7s9V6LV73mDM8Hsm2ejz+87yzv8Bu/w/1uXCy46qqrrvrPdOzczVx11f91l07fzlVXXXXVf6Y7l0vu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gdrrrqqqv+q5y49yFcddX/dRevfTpXXXXVVf9Vfuz1XotXveYMzweybZ6PvXHiFX/mF9gfR6666qqr/iucuPchXHXV/3UXr306V1111VX/Fba7jj99qzdjp6s8H8i2eQE+8y//hu944pO56qqrrvqvcOLeh3DVVf/XXbz26Vx11VVX/Vf42Bd/LB/3Eo/lBUC2zQvxmX/5N/zyHXdy5+ERV1111VX/mU7c+xCuuur/uovXPp2rrrrqqv9MN25u8MY33cjnvuxL8UIg2+aqq6666qqrrrrqqquu+v8E2TZXXXXVVVddddVVV1111f8nyLZ5IT7rr/6GX779Lu44POSqq6666j/TiXsfwlVX/V938dqnc9VVV131n+mmzU3e+OYb+JyXeSleCGTbvACf+Zd/w3c88ckA3LhYcNVVV131n+nYuZu56qr/6y6dvp2rrrrqqv9Mdy6XALzfox7B577sS/ECINvm+dgbJ17pZ36RvXEA4BVPnuSqq6666j/T9sXrueqq/+v2T9zNVVddddV/pj+9cAGAna7nT97qTdnpKs8Hsm2ejz+87yzv8Bu/w/1e8eRJrrrqqqv+M21fvJ6rrvq/bv/E3Vx11VVX/Wf60wsXuN+Pvd5r8arXnOH5QLbN8/GH953lHX7jd7jfjYsFV1111VX/mY6du5mrrvq/7tLp27nqqquu+s9053LJ/X7s9V6LV73mDM8Hsm2ejz+87yzv8Bu/w1VXXXXVf5UT9z6Eq676v+7itU/nqquuuuq/yo+93mvxqtec4flAts3zsTdOvOLP/AL748hVV1111X+FE/c+hKuu+r/u4rVP56qrrrrqv8J21/Gnb/Vm7HSV5wPZNi/AZ/7l3/AdT3wyV1111VX/FU7c+xCuuur/uovXPp2rrrrqqv8KH/vij+XjXuKxvADItnkhPvMv/4ZfvuNO7jw84qqrrrrqP9OJex/CVVf9X3fx2qdz1VVXXfWf6cbNDd74phv53Jd9KV4IZNtcddVVV1111VVXXXXVVf+fINvmqquuuuqqq6666qqrrvr/BNk2L8Rn/dXf8Mu338Udh4dcddVVV/1nOnHvQ7jqqv/rLl77dK666qqr/jPdtLnJG998A5/zMi/FC4FsmxfgM//yb/iOJz4ZgBsXC6666qqr/jMdO3czV131f92l07dz1VVXXfWf6c7lEoD3e9Qj+NyXfSleAGTbPB9748Qr/cwvsjcOALziyZNcddVVV/1n2r54PVdd9X/d/om7ueqqq676z/SnFy4AsNP1/MlbvSk7XeX5QLbN8/GH953lHX7jd7jfK548yVVXXXXVf6bti9dz1VX/1+2fuJurrrrqqv9Mf3rhAvf7sdd7LV71mjM8H8i2eT7+8L6zvMNv/A73u3Gx4KqrrrrqP9Oxczdz1VX/1106fTtXXXXVVf+Z7lwuud+Pvd5r8arXnOH5QLbN8/GH953lHX7jd7jqqquu+q9y4t6HcNVV/9ddvPbpXHXVVVf9V/mx13stXvWaMzwfyLZ5PvbGiVf8mV9gfxy56qqrrvqvcOLeh3DVVf/XXbz26Vx11VVX/VfY7jr+9K3ejJ2u8nwg2+YF+My//Bu+44lP5qqrrrrqv8KJex/CVVf9X3fx2qdz1VVXXfVf4WNf/LF83Es8lhcA2TYvxGf+5d/wy3fcyZ2HR1x11VVX/Wc6ce9DuOqq/+suXvt0rrrqqqv+M924ucEb33Qjn/uyL8ULgWybq6666qqrrrrqqquuuur/E2TbXHXVVVddddVVV1111VX/nyDb5oX4rL/6G3759ru44/CQq6666qr/TCfufQhXXfV/3cVrn85VV1111X+mmzY3eeObb+BzXualeCGQbfMCfOZf/g3f8cQnA3DjYsFVV1111X+mY+du5qqr/q+7dPp2rrrqqqv+M925XALwfo96BJ/7si/FC4Bsm+djb5x4pZ/5RfbGAYBXPHmSq6666qr/TNsXr+eqq/6v2z9xN1ddddVV/5n+9MIFAHa6nj95qzdlp6s8H8i2eT7+8L6zvMNv/A73e8WTJ7nqqquu+s+0ffF6rrrq/7r9E3dz1VVXXfWf6U8vXOB+P/Z6r8WrXnOG5wPZNs/HH953lnf4jd/hfjcuFlx11VVX/Wc6du5mrrrq/7pLp2/nqquuuuo/053LJff7sdd7LV71mjM8H8i2eT7+8L6zvMNv/A5XXXXVVf9VTtz7EK666v+6i9c+nauuuuqq/yo/9nqvxatec4bnA9k2z8feOPGKP/ML7I8jV1111VX/FU7c+xCuuur/uovXPp2rrrrqqv8K213Hn77Vm7HTVZ4PZNu8AJ/5l3/DdzzxyVx11VVX/Vc4ce9DuOqq/+suXvt0rrrqqqv+K3zsiz+Wj3uJx/ICINvmhfjMv/wbfvmOO7nz8Iirrrrqqv9MJ+59CFdd9X/dxWufzlVXXXXVf6YbNzd445tu5HNf9qV4IZBtc9VVV1111VVXXXXVVVf9f4Jsm6uuuuqqq6666qqrrrrq/xNk27wQn/VXf8Mv334XdxwectVVV/3/MDzuHzj3+Z/N7DEvxqnP+Gz+q5y49yFcddX/dRevfTpXXXXVVf+Zbtrc5I1vvoHPeZmX4oVAts0L8Jl/+Td8xxOfDMCNiwVXXXXV/w/7f/93POEzPoWdF38JHvV5X8R/lWPnbuaqq/6vu3T6dq666qqr/jPduVwC8H6PegSf+7IvxQuAbJvnY2+ceKWf+UX2xgGAVzx5kquuuur/h3N/+zf8/id+LKdf8qV49S/9Sv6rbF+8nquu+r9u/8TdXHXVVVf9Z/rTCxcA2Ol6/uSt3pSdrvJ8INvm+fjD+87yDr/xO9zvFU+e5Kqrrvr/4dzf/g2//4kfy+mXfCle/Uu/kv8q2xev56qr/q/bP3E3V1111VX/mf70wgXu92Ov91q86jVneD6QbfN8/OF9Z3mH3/gd7nfjYsFVV131/8P+3/8dT/iMT2HnxV+CR33eF/Ff5di5m7nqqv/rLp2+nauuuuqq/0x3Lpfc78de77V41WvO8Hwg2+b5+MP7zvIOv/E7XHXVVf//DI/7B859/mcze8yLceozPpv/KifufQhXXfV/3cVrn85VV1111X+VH3u91+JVrznD84Fsm+djb5x4xZ/5BfbHkauuuur/l+Fx/8C5z/9sZo95MU59xmfzX+XEvQ/hqqv+r7t47dO56qqrrvqvsN11/OlbvRk7XeX5QLbNC/CZf/k3fMcTn8xVV131/8vwuH/g3Od/NrPHvBinPuOz+a9y4t6HcNVV/9ddvPbpXHXVVVf9V/jYF38sH/cSj+UFQLbNC/GZf/k3/PIdd3Ln4RFXXXXV/w/D4/6Bc5//2cwe82Kc+ozP5r/KiXsfwlVX/V938dqnc9VVV131n+nGzQ3e+KYb+dyXfSleCGTbXHXVVVddddVVV1111VX/nyDb5qqrrrrqqquuuuqqq676/wTZNlddddVVV1111VVXXXXV/yfItrnqqquuuuqqq6666qqr/j9Bts1VV1111VVXXXXVVVdd9f8Jsm2uuuqqq6666qqrrrrqqv9PkG3zInjxr/4L/jX+/qNfjquuuuqqq6666qqrrrrqfyRk27wIXvyr/4J/jb//6Jfjqquuuuqqq6666qqrrvofCdk2L4IX/+q/4F/j7z/65Xig3d1dvud7vofd3V1e+qVfmpd6qZcC4MEPfjC7u7v8zd/8Dfd70IMexIMf/GCe26233srP/MzPsLu7y0u/9EvzUi/1Uhw/fpzjx49z1VVXXXXVVVddddVVV73IkG3zInjxr/4L/jX+/qNfjvv99m//Nm/zNm/D7u4uD/RWb/VW/PRP/zSf/dmfzed8zudwv7d6q7fip3/6p3mg7/7u7+Z93ud9eG6f9VmfxWd/9mdz1VVXXXXVVVddddVVV73IkG3zInjxr/4L/jX+/qNfDoDd3V0e8pCH8FIv9VL89E//NMePH+ev//qveeu3fmse/OAH89u//dvc7yEPeQi2ecYznsHTn/50HvzgBwPw13/917zMy7wM7/Ve78VXf/VXc/z4cW699VZe+qVfmo/+6I/msz/7s7nqqquuuuqqq6666qqrXmTItnkRvPhX/wX/Gn//0S8HwFd/9VfzMR/zMVy8eJHjx49zv6/+6q8G4KM/+qMB+Omf/mne5m3eht/6rd/ird/6rfnoj/5oPvuzPxuAj/7oj+anf/qnufXWW3mgn/7pn+bBD34wL/3SL81VV1111VVXXXXVVVdd9SJDts2L4MW/+i/41/j7j345AD77sz+bn/7pn+av//qveWHe+q3fmltvvZW//uu/5r3f+735nd/5HZ7+9KcD8Nqv/doA/PZv/zZXXXXVVVddddVVV1111b8bsm1eBC/+1X/Bv8bff/TLAfDZn/3Z/MzP/Ax/9Vd/xQty66238pCHPIQHP/jBPPjBD2Z3d5e//uu/5qd+6qd467d+a177tV8bgN/+7d/mr//6r3md13kddnd3ATh+/DgXL17kqquuuuqqq6666qqrrnqRIdvmRfDiX/0X/Gv8/Ue/HADf/d3fzfu8z/vwV3/1V7z0S78099vd3WV3d5cHP/jBfPZnfzZf/dVfzUu/9Etzv7/+67/mtV/7tfnpn/5p3vu935vf+Z3f4elPfzoAX/3VX83u7i7f/d3fzUu/9Evz0z/901x11VVXXXXVVVddddVVLzJk27wIXvyr/4J/jb//6Jfjfi/90i/NpUuX+OzP/mwe/OAHs7u7y2d/9mdz8eJFvvu7v5v3eZ/3AeC7vuu7eO3Xfm1uvfVW3ud93off/u3f5ru/+7t50IMexOu8zuvw1m/91nz0R380ALfeeiuf/dmfzXu913vx2Z/92Vx11VVXXXXVVVddddVVLzJk27wIXvyr/4J/jb//6Jfjfrfeeisf/dEfzc/8zM9wvwc96EH89E//NO/93u/N3/zN3wDwVm/1Vvz0T/80X/3VX83HfMzHAHDs2DF2d3f56Z/+aT76oz+aZzzjGdzvpV7qpfjqr/5qXvu1X5urrrrqqquuuuqqq6666kWGbJsXwYt/9V/wr/H3H/1yPD+//du/zYMf/GAe/OAH829x6623cuutt/LgBz+YBz/4wVx11VVXXXXVVVddddVV/2rItrnqqquuuuqqq6666qqr/j9Bts1VV1111VVXXXXVVVdd9f8Jsm2uuuqqq6666qqrrrrqqv9PkG1z1VVXXXXVVVddddVVV/1/wj8CIB1K/qqj+MwAAAAASUVORK5CYII=",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import Image\n",
    "Image(filename=\"output/TR1.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4fdef400-7dcc-43a0-ac38-9d3bd923302a",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}