{ "cells": [ { "cell_type": "markdown", "id": "b151e349", "metadata": {}, "source": [ "# JUMP CRISPR arm\n", "\n", "{func}`~mantispy.ds.jump_crispr` serves the CRISPR arm of `cpg0016-jump-assembled` {cite:p}`Chandrasekaran_2023`, 51,185 wells of knockouts in U2OS cells, corrected for well position and cell count and feature-selected by jump-profiling-recipe. Each gene is knocked out by several guides, and the guides are its replicates, so this is an arrayed, well-level CRISPR screen, the counterpart to the single-cell optical pooled screens {func}`~mantispy.ds.scallops_arv471` and {func}`~mantispy.ds.cp_posh`.\n", "\n", "**Use it if** your own data is an arrayed (well-level) genetic screen, knockout or otherwise, where several reagents share a target and act as its replicates. [CRISPR knockouts](../case_studies/crispr.ipynb) runs the knockout workflow on it, and [Published profiles and JUMP](../tutorials/data/profiles.ipynb) loads it." ] }, { "cell_type": "code", "execution_count": 1, "id": "34d3ea89", "metadata": { "execution": { "iopub.execute_input": "2026-10-21T00:58:14.478814Z", "iopub.status.busy": "2026-10-21T00:58:14.478443Z", "iopub.status.idle": "2026-10-21T00:58:19.567217Z", "shell.execute_reply": "2026-10-21T00:58:19.565583Z" } }, "outputs": [], "source": [ "import numpy as np\n", "import scanpy as sc\n", "\n", "import mantispy as mt" ] }, { "cell_type": "markdown", "id": "ad49b038", "metadata": {}, "source": [ "The well-level object is large, so this page uses the gene-level consensus, one `modz` profile per gene over its guides." ] }, { "cell_type": "code", "execution_count": 2, "id": "871e11cc", "metadata": { "execution": { "iopub.execute_input": "2026-10-21T00:58:19.571310Z", "iopub.status.busy": "2026-10-21T00:58:19.570822Z", "iopub.status.idle": "2026-10-21T00:58:19.948176Z", "shell.execute_reply": "2026-10-21T00:58:19.947078Z" } }, "outputs": [ { "data": { "text/plain": [ "AnnData object with n_obs × n_vars = 7974 × 595\n", " obs: 'Metadata_Gene', 'Metadata_ReplicateCount', 'Metadata_Source', 'Metadata_JCP2022', 'Metadata_Control_Type', 'Metadata_Perturbation', 'Metadata_Perturbation_Type', 'Metadata_Control', 'Metadata_ChromosomeArm'\n", " var: 'object', 'feature_group', 'feature', 'channel', 'scale', 'angle', 'gray_levels', 'radial_bin', 'params', 'is_feature', 'feature_kind'\n", " uns: 'mantispy'\n", " layers: None (.X)" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "adata = mt.ds.jump_crispr(aggregated=True)\n", "adata" ] }, { "cell_type": "markdown", "id": "0c8dcd98", "metadata": {}, "source": [ "## Controls against treatments\n", "\n", "`Metadata_Control_Type` labels each gene consensus a negative control (`negcon`), a positive control (`poscon`) or a treatment (`trt`). On the first two principal components the controls should sit apart from the bulk of the genes." ] }, { "cell_type": "code", "execution_count": 3, "id": "ce04f4dc", "metadata": { "execution": { "iopub.execute_input": "2026-10-21T00:58:19.951337Z", "iopub.status.busy": "2026-10-21T00:58:19.951185Z", "iopub.status.idle": "2026-10-21T00:58:20.573006Z", "shell.execute_reply": "2026-10-21T00:58:20.572250Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAApcAAAGuCAYAAADF8S1pAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjIsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvgI3uAAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAl/9JREFUeJzs3XV4VEcXwOHf7mbj7oFAIAQPXlyLQ4HitJQiNb56aakb9VKhRt0oFbRAcW1wL24hCSTE3XXlfn+E3LJNCAkkJMB5n6dPm3tn5s5dKHsYOaNRFEVBCCGEEEKIKqCt6Q4IIYQQQoibhwSXQgghhBCiykhwKYQQQgghqowEl0IIIYQQospIcCmEEEIIIaqMBJdCCCGEEKLKSHAphBBCCCGqjASXQgghhBCiykhwKYS4IcTHxzNv3jwSEhJquivX7NChQ/zxxx/MmzePwsJCtm3bxtKlSy3K7Ny5kyVLltRYH4UQ4mpp5IQekZqayqpVq+jQoQOtWrUCICwsjF27djFgwADq1q1bqs6ePXsIDQ1lwoQJ2NnZqW0A9OzZk0aNGpWqc/r0afbt2wfA2LFjcXR0tHh+Ca1Wi6urKx06dCjz2ZdjNps5evQo586dQ1EU6tevT/PmzXFycir1rv99Vrt27ahXr16FPpsSmZmZHDt2jLi4OHx9fQkODsbDw6PM+ld6t8r261a0efNmBgwYQEhICH369Knp7ly1+++/nxUrVjBw4EBsbGyYO3cuDzzwAHv37iUyMlItN3XqVDZv3kxMTEyN9lcIISpNEbe8AwcOKIDy3nvvqde+//57BVDWrVtXZp3p06crgBIfH2/RBqBMnTq1zDojRoxQy5w/f77U89u3b69MmTJFmTRpktKuXTtFq9UqDzzwgGIwGK74Dj/99JPi7++v+Pj4KHfeeacyfvx4pWnTpoqDg4Ny//33KwUFBZd9VseOHRWNRqNMnjy51LPK+mzMZrPy5ptvKk5OTkqLFi2UcePGKT169FCsra2VcePGKQkJCZV+t8r261Z08uRJZcqUKcrp06druitX7dixYwqgrFixwuL6V199pTzzzDMW16ZMmaLUrVv3OvdQCCGunQSXokqDy0aNGimOjo5KTk6ORfnExETFyspKCQoKumxweenzFUVRXnnlFQVQPvnkk3L7/8orrygajUZ59913SwVha9asUXx9fZXk5ORyn/Xuu+8qgPLxxx9f8bP57rvvFECZPXu2Rdn9+/crjRo1Uvbs2VPpd6tsv8SNaeHChQqgnDhx4oplJbgUQtyoZM2lqFITJ06koKCg1PqxX3/9FTs7O0aNGlXhtp555hkA1qxZc9kyBw4c4O233+Z///sfL774IlZWVhb3hw4dyrZt27Czsyv3WY8++igA69evv2K//vrrL6ysrHjqqacsrnfs2JFdu3ZVaCq/Iu9W2X6ViIqKYvHixaxbt47CwkIyMzOZN28e4eHhpcoWFhaydetW/vjjDzZu3EheXl6ptubNm0dqaiomk4ktW7awaNEiIiIiLvv8K7UJYDAY2LFjB4sWLWLbtm1kZWVd8b3KWnN5Nf273PsZDAY2btzI4sWLLaanL/eskJAQ5s+fb/F+oaGhLF68mKVLl3L+/HmL+suXL2fDhg0ArFy5knnz5ln889//Z8qTkJDAypUrWbhwIYcOHapwvcTERNavX8+SJUs4cuQI/10V9ffff7N8+XIA4uLiWLp0KWvXriU3N7dUW//884/a9/nz5/PXX3+VO4WflJTE6tWrWbJkCWfOnKnS9xJC1C5WFSgjRIV5e3szdOhQ5s2bx5QpU9Tr8+bNY8KECdjb21e4LVtbWzQaDQUFBZct8+233wIwY8aMy5Zp0qTJFZ9lbW2NRqOhqKjoimWdnZ0xmUxkZGTg7e1tcc/Hx+eK9angu1W2XwBvvvkmb775Jh06dMDX15dXX32VN954g2nTpvHzzz8TFBSkll25ciUPPfQQDg4OtGvXjtDQUJKSkvj999/p378/XAzep02bxurVq5k9ezYuLi6kpKSwf/9+Pv7441IBdkXaPHLkCMOGDcPOzo727duTmppKREQEb731FpMmTbrsu508eZJp06YREhKCr6/vVfXvv0rqL1u2jPfeew8fHx+1/ssvv8ybb75Zquzy5cuZPXs27u7unD59mr59+2Iymbj33nvZsGEDPXv2xGg0smvXLiZMmMD333+PjY0N+/fvJzQ0FID9+/fj4uKitr1t2zYURWHs2LHl9tdoNPLMM8/w9ddf07lzZzw9Pdm5cyetWrViyZIlpdb8Xmr27Nm88cYbdOzYET8/P8LDw9FoNCxatIjAwEAAPv/8c86cOUNBQQHvvPMOTZo04eDBgxQVFbFs2TK6deumtnfu3Dm2bt0KgMlk4vz58+zdu5fHH3+cTz75RC1nMBh45pln+Oabb2jTpg0BAQG8/vrrBAUFsXTpUqytra/pvYQQtVBND52KmleV0+JffPGFsmzZMkWj0ahT3yX3du3apbz++usVnhZfunSpApRai3apli1bKu7u7tf0roqiKMuXL1cA5bnnnrti+Y0bNyoajUZp166dsnjxYiUpKanSz/vvu1W2X2VZuXKlAigffvihei02Nlbp3bu3Aig///yzRb/0er3yv//9TzEajYpycS3po48+qjg5OSlxcXGKoijKkiVLFEDp2bOnEhkZqdZ/8MEHFXt7eyUxMbHSbfbt21fp2rWrYjKZ1LoZGRnKmjVryn2/TZs2KYASEhKiXqtM/8pSUr9Lly7KuXPn1OuffvqpAihLliwpVbZz585KRESEoiiKkp6ermRmZir33HOPYmtrq+zfv18tv2XLFkWn0ylPPvmkeq3k/6v/rhudMGGCEhAQYHGtrGnxmTNnKlZWVsrmzZvVa4mJiUpQUJAycuTIy75nbGysAiiff/65xfXDhw8rp06dUn++8847FQ8PD+W+++5Tl5jk5uYqvXr1Ury9vZXMzMxyP88NGzYogLJ69Wr12lNPPaXodDqLa8rF3695eXnX9F5CiNpJpsVFlRs2bBgeHh788ssvAPz88880btzYYtSjLCXTbD/++CMzZ85k8uTJNG/enOeff/6ydVJTU3F3d690Hy991rPPPsukSZPo0qVLuc8qMWDAANatW4ebmxv33HMP3t7e1K1bl2nTpnH8+PFrerdr6ddXX32Fv7+/xShunTp1GDduXKmys2fPxs7Ojk8++QSdTgeARqPh3XffJS8vj/nz51uUHzlyJAEBAerPDz74IHl5eezatavSbSYmJqojtyVcXFwYOnToFd/xcirSv/LceeedNGzYUP358ccfp1GjRsydO7dU2eHDh6sjfa6urhQVFbFgwQKmTJlCx44d1XJ9+/ZlzJgxfPvtt1ccoa6IrKwsvvjiCyZPnky/fv3U697e3jz77LOsWLHistPSSUlJAKWWh7Rt25bmzZtbXEtNTWXWrFnqEhN7e3vee+89kpKSWLx4sUXZoqIiduzYwYIFC5g3bx5xcXE4ODioI5oZGRl8+eWX3Hvvvdxxxx0WdYcPH46dnd01vZcQonaSaXFR5fR6PRMnTmT+/Pk8//zzLFiwgGefffaK9Uqm2bRaLS4uLnz//feMHj0aW1vby9ZxdHQscz1YRZ9VWFjIwYMH0Wq1vP766xUOVAcNGsSgQYMoKCjgzJkzhISE8Omnn/Lbb7+xfPlyhg0bdlXvdi39OnbsGG3atFEDuxJt27YtVXbfvn34+Pio6/wubu5DURScnZ05efJkuW34+/sDWHzpV7TNBx98kKeeeoqWLVsycuRIevfuTa9eva64LrY8Felfedq3b2/xs1arpV27dmzevPmKZU+cOIHZbLYILEt07NiRxYsXExoaSps2bSrUl8s5cuQIhYWFmM1mfvvtN3W9pKIo6vrOU6dOqe9+qeDgYLp27cr06dNZvnw5AwcO5Pbbb6dVq1YWQT6Al5dXqfRX7du3R6PRcPToUfXa2rVrmTp1KjY2NnTo0AEXFxc0Gg2KoqjrYo8ePYrBYCj3L5bX8l5CiNpJgktRJhsbG7g4MlGWwsJCuLh2sCxTp07l888/5+mnnyYzM5PJkydf8Znjxo3jhRdeqFQ/W7VqxfLly0lNTa3UuqxLn2U2m7nrrrsYOXIku3fvLhU8lMfW1pa2bdvStm1bJk6cSKNGjXj11VdLBZcVfbdr6ZfRaFR/3f7bx/8q+fUrK3gaMWJEqUDJ2dnZ4me9Xm/RTmXafPLJJ+nSpQuLFy9m8+bNfPDBBzg6OvLpp58yderUct/xcirSv/KU9bnZ2NhgMBhKXffy8rL4uaRMeZ99We1UVsm7nD9/vszPeMqUKXh6epZZ18rKim3btrFo0SLWrVvHnDlzeOqpp2jTpg1LliyhcePGatmy3kOv16PVatX3yM/P55577qFnz56sWLECrfbfSbAlS5aoAWJJ+fL+gngt7yWEqJ0kuBRlatCgAcBld92Gh4fj5uZmsSnhUu3ataNNmzZ8/fXXDBo0qFLJ0Ctj4sSJLF++nPnz5192U09ubi56vR5ra+sy72u1Wr755hu2bNnC9OnT2b9/f6nRnEulpKSU+WXn4+NDQEBAmTuNr0Zl+xUYGFjmjvCyrjVp0oTk5GTmzZtXJX2tbJudO3emc+fOACQnJzNhwgQefPBBi+T611NYWBi9e/cuda1k+rs8JWXOnj1b6l5oaCgajcZiyv1qlWxM69+/P6+88kql6+v1eiZNmqRumtq9ezcDBgzgpZdesjgJKD4+ntzcXBwcHNRr586dw2Qyqe8aGhpKRkYG48aNswgsw8LCLHbPN23aFC5uxqqu9xJC1D6y5lKUqUuXLgQEBPDLL7+Qn59vce/IkSPs2bOHCRMmlBvsvP7660yZMoXnnnuu2vo5evRoBg0axOuvv86ePXtK3Y+Pj2fIkCFXTHXj7u7Oc889x8GDB6+YEubhhx/mp59+KnU9NDSUsLAw2rVrdxVvcu39mjBhAidOnCAkJES9ZjKZ1LWvl3r44YcJDQ3ljz/+KHUvKyuL+Pj4Sve1om2W7Jgu4eXlpe6wzsnJqfRzq8L8+fMtRul37drF/v37mThx4hXrNmrUiM6dO/Pjjz+SkZGhXk9ISODXX39l8ODBVbLbOSAggCFDhvDll1+SmJhY6v5/P9dLJSQkkJmZaXGtW7du+Pn5WfSZixkKvvzyS4trc+bMQa/Xq+t3S/6yGBYWppZRFIV33nnHYpSyXr16DBo0iG+//Zbo6GiLNuPj4zEYDNf0XkKI2klGLoXq0kBRr9ezcOFCRowYQdu2bbn77rvx9vYmNDSUH374gS5dujB79uxy2xs1alSl8lpeDa1Wy59//snDDz9Mr169GDlyJD169MDa2ppjx47xxx9/4OPjU+60XInHH3+cTz/9lFdffZXRo0eXWrtYIigoiOnTp/Pll18yZswY3NzcOHv2LPPnz8fX17fUF/O1qmi/HnvsMdasWcPw4cN57LHH8PX1Zc2aNYwYMYJ169ZZ/PpOnDiR06dPM3nyZP766y+6dOmC2Wzm9OnTbNq0iQULFuDn51epfla0zXvvvRdXV1c1uAkPD+err77i4YcfVlMMXW+jR49m4MCBjBgxgpSUFObOnUvv3r15+umnK1T/559/pn///nTq1Ilp06ZhMpn47rvv8PLyUtNlVYVffvmF4cOH07JlS6ZNm0ZgYCCJiYns3buXpKSky+aGjIqKYty4cfTv35/g4GBsbGxYt24dcXFxfPPNNxZl69evT0JCAvfffz9t27Zl+/bt/Pnnn3z++efqCKyXlxcPPfQQH3zwAbm5udSvX5+1a9fSq1evUoH0vHnzGDx4MG3btuX+++8nICCAU6dOsXnzZg4fPoxer7/q9xJC1E4ycinUEY1Lp8G4OHoZHh7OCy+8QHp6OocPH8be3p6lS5eyfft2i3Vunp6eTJkypdTO0/9q27YtU6ZMsZj6LKn737O7K8rBwYH58+dz/PhxunbtSlhYGCdPnqRhw4asWLGC0NBQ9XnlPcve3p65c+fSpUsXTp06ddny7733HnFxcTz66KNkZWVx6NAh9Ho9X3/9NaGhoRafQUXfrbL9KouVlRXr16/nm2++ITMzk/j4eObMmcNtt90GFzc/Xeqtt97ixIkTtGvXjjNnzpCYmEiPHj04ceKEugGjQYMGZa55s7W1ZcqUKQQHB1e6zf379/Piiy+i0Wg4cuQIDg4O/P3333z11VflfkZ16tRhypQpFgFoZft3OT169OCrr74iKSmJ9PR0vvjiCzZv3myxlOJyzwJo3rw5p0+f5umnnyY6Opr4+HhmzZrFsWPHLDbHNGnShClTppRaTtK7d+9SOS579OjB+PHjLa55eXmxe/du5s+fj1ar5ciRI+h0OmbOnFluANa5c2dOnTpF7969iY2N5ezZswwcOJDIyEg1/+il5syZw+DBgzlz5gwNGjRgz549PPbYYxZlvv32W3777TdMJhNRUVE899xzvPzyy4wfP55evXqp5Xx9fTl48CDffPMNBoOB06dP07p1aw4dOqTmvb3a9xJC1E4a5b9HNIhbzgcffMDzzz/Ptm3bLL4UxM3hl19+YerUqRw5cuSadyzfbJYuXcq4ceM4cOCAGoTfykaOHMmZM2cue4KOEEJUhEyL38LWrVvHrl27+PTTT7n99tvp2bNnTXdJXKO4uDjq1Kmj/lxQUMDnn39O48aNad26dY32TQghxK1Bgstb2LFjx0hPT2fu3Lncfffd5W7OETeGRx99FDs7Ozp16kROTg6LFi0iKiqKtWvXyq+vEEKI60KCy1tYRU59ETeWJUuWsHz5cg4ePEhmZibTpk1j4sSJNbZRprYrbx3lrahv375XXDcthBBXImsuhRBCCCFElZHd4kIIIYQQospIcCmEEEIIIarMdV9zaTabiYuLw8nJSTYYCCGEEDcIRVHIzs6mTp06Fsd+CvFf1z24jIuLs0gqLIQQQogbR3R0NP7+/jXdDVGLXffg0snJCS7+5rz0hBchhBBC1F5ZWVnUq1dP/R4X4nKue3BZMhXu7OwswaUQQghxg5ElbeJKZNGEEEIIIYSoMhJcCiGEEEKIKiPBpRBCCCGEqDISXAohhBBCiCojwaUQQgghhKgyElwKIYQQQogqI8GlEEIIIYSoMhJcCiGEEEKIKiPBpRBCCCGEqDISXAohhBBCiCojwaUQQgghhKgyElwKIYQQQogqI8GlEOKmtztuNx8d+IikvKSa7kqVMZlNLDm7hB0xO2q6K0IIYcGqpjsghBDV7amQp8g35pNZlMlb3d+q6e5YyDXkYqW1wkZnU6l6y8OX8+aeN9FqtKwetZp6TvWqrY9CCFEZMnIphLjptfVqC0AbrzY13RULp1JPcfvi2xny5xBS81MrVdfPwQ+tRouLtQtOeqdq66MQQlSWjFwKIW563wz4hlxDLk7WtSsIC88IJ9+YT74xn4TcBDzsPCpct3vd7qwZtQZHvSOutq7llr3a0VEhhLgaElwKIW5IRrORl3a8xIXsC3zQ6wPqO9e/bFmtRlvrAkuAIQ2HkJSXhIuNCy09W1a6vr+T/xXLnEg5wbT103C2cWbp8KW42bpdZW+FEKJiJLgUQtyQLmRfYF3kOgA2RW3i/lb313SXKk2v1fNAqweq9Rln089SYCqgIK+ApLwkCS6FENVOgkshxA2pgXMDxjcZT1R2FEMbDq2xfpzPPI+dlR2+Dr411ofyDA8cTmp+Kh52HjR1b1rT3RFC3AI0iqIo1/OBWVlZuLi4kJmZibOz8/V8tBDiBnYk6QjzT81ndOPR9Kjbo6a7A8A/if9w34b7sNHZsPzO5dR1rFvTXRKi2sj3t6go2S0uhLghfPLPJ2yK2sTs/bNruiuqjIIMzIqZAmMBuYZcjicfZ/Cfg5m5bSZmxVzT3RNCiBohwaUQ4oYwsMFA7K3sGdRgUI0832Q2lbrWL6AfH/X+iG8GfEMTtyZsurCJ2JxYNkRuIKMwo0b6KYQQNU2CSyHEDeGe5vew7559PNbusSpvO7Mwk4TchMven7V7Fu1+bcfvp38vdW9Qg0F0q9MNgPFNxtPbvzePtH0Ed1t3tczGyI28vPNlorKi1GvrI9fTY2EP3t77dpW/jxBC1CQJLoUQt7TU/FSGLR/GkD+HcCDhQJlldsTuQEFhd9zuctvyd/Jnbr+5PNzmYYvrr+x6hZURK/nqyFcAHEw4yLt73yWzMJPV51ZX4dsIIUTNk+BSCHFLyyzKJKMwA6NiJCY7xuJeWHoYF7Iu8Fa3txgVNIqnOzx9Vc/oU68Peq2env49AfjxxI+kF6aj1+p5ufPLVfIeQghRW8hucSHELe/vC3+Tkp/CmMZj0Gl1cHEn+LT109Br9Sy7cxkBzgGVbnf+yflsiNzAjA4zuM33NvX6pqhNfHzwY+4MurPUKKcQtZV8f4uKkpFLIcQN54fjPzBo6SDWnV9XJe31rd+X8U3Hq4ElF49MVFAwmA0UGAuuqt2vjn7FsZRj/Hb6N4vrAwIGsH7Meh5u8zC7Y3czYOmAall7mW/M58+zfxKeHl7lbQshxOVIcCmEuOEsPLOQuNw4/gz7s9qe4e/kT996fXm1y6tXlXw815DLlBZTCHINYlyTcQAoisLZ9LMUmYrg4kai13e/TkJuAotDF1PVE0mfH/qcWXtmMW3DtCpvWwghLkeCSyHEDefJ9k/SybcTD7V6qEraKyvwemvPW/wd/TcLQhdUuq3HtjxGlz+6sCFyA0uGL6F73e4AfHjwQ8asHMMTIU8AsC1mGwl5xbvUO/t1RqPRYDAbWB+5nujs6Gt+r5Id6642rmg0GgAKjAX8ePxHQi6EXHP7QghRFjn+UQhxwxneaDjDGw2vkrbWnV/Hyztfpk+9PszpM0e93sqrFQcTD9LaszUp+SkUmYqo41jniu1FZ0WzLWYbAOcyz/H76d9ZcnYJYxqPIS4nDkD9d1e/rtR3qs+F7Avsjd/LydSTbIrcxI8nfsTd1p2t47eqQeHlxGTH4G3vjbXOutS9B1s/SNc6XanvXF+99seZP/j00KdoNVo2jd2Et713JT4tIYS4MgkuhRC3tJ2xOzGYDWyL3qZeyy7K5mzaWTr5dmJis4ncsewOikxF/DLkF1p7tSYpL4mMwgwauza2CP6e3vo0W6K20NKjJTE5MdjqbPno4EcAzPlnDlvGbqGddztcbFwYumwovf17M7vXbKasm4K93h4vOy9sdDYA2OhsrhhYzjsxj4//+Zi2Xm35deiv6vU8Qx7ZRdn4OPgQ7BlsUcdgMqDVaPFz8MNR71iqzfD0cOo41sFeb38Nn6oQ4lYmwaUQ4pYRnRVNgamAxm6N1Wv/a/M/NGjo4f/veeX7E/azK24XAN8c/YY8Yx4AyXnJHEk6wuR1k1FQGNFoBO/0eEettzN2J2bMuNu642ztzJ74Peq9jr4dKTIXMaXlFF7Y/gLR2dH8dvo3wtLDGBU0iic7PImTtRPT20ynvU97Grk2suh7riGXr498TQOXBowMGsmp1FOEpYfBxRHSEnmGPIYtG0ZKQQqf3f4Zt9e/Xb238MxC5h6Zi4OVA78N+a1UAPnLyV/46OBHBLkGsfzO5VXymQshbj0SXAohbgnRWdGM/GskBrOBHwf9SEffjgDUc6rH2z0sd2p39u1ML/9emMwmNl3YBEA3v270C+jHqohVKBSv0Vx7fi1Hko7wRd8vCHQN5P2e77M9ZjsHEw4Slf3vaTx/DvuTp7Y9xZBlQ3iv53u42bqp9/Yl7GNfwj4uZF+gR90etPBowdrza7mr6V142nnCxXWcD29+mMNJhwE4kHCAtefX0tKjJU+0e4LOfp0BOJlyknf3vUtyQTIAkVmRFu9Vct65RqPB1sq21GcUnxsPQGJuImbFjFYjy/KFEJUnwaUQ4paQb8rHYDagoJBTlFNuWUdrR77s92Xx5py/H+Nw4mH6BvTl0S2P0rdeXx5q/RCLQxeTUZjBhewLrD63mifaP0Hf+n3pW78vL+14iajsKDRoaOfdDl9HX/V4ycWhi9UgMcApgEJTIQl5CeyJ38Oe+D342vuSkJdAZGYkvwz5BYDEvES1jpuNG+czzwNwKvUUn97+Kb4OvsTmxDJ983QyCzMBcNI7MaLRCIv3urvZ3fg7+VPfqT6O1pZT4nmGPKa3nk49p3q09W4rgaUQ4qpJEnUhxE3JYDJgVIzYWdmp1w4kHCCnKMdiqriiZm6byYbIDdhZ2bH/nv38efZP3tn3Dm62bvww8AcaujRUy5oVMyn5KbjYuGCttUaj0bAjZgdLzi4huyibg4kHAXiu43MMDBjImJVjyCzKxFZnS2uv1uxP2M/j7R4n2DOYx7Y8ho3Ohh51exCZFck7Pd5h3bl1/HDiB7QaLRvHbCStII2719yNSTEB0MarDSdSTuBq48qKO1fgautKUl4S3x79llxDLnUc6/Bwm4fR6/RwccRz6vqpONs48+fwP3G1da2CXwFxs5Hvb1FRMnIphLjhJeclM2ntJIyKkV+H/IqdlR3jVo0jqyiLeYPn0cKjBVxc9/j3hb+ZsHoCE5pOYHTj0RV+xpCGQziYcJCBDQYCMKbJGMY0GVOq3MGEgyw+uxhrrTWh6aG82OlF2vu050jyEUKiQ7DWWjMicARWOiuWhy1na/RWVo5cSZGpiPcPvM+WC1uw09nRzrsdC84swGA2YDAbOJ5ynMS8RB7a+BDtfdrDxSB23fl1JOclq4Hlsx2fxd7KnqPJR0ktSCW9MB1XW1e+P/Y9i88uVvvZ0KWhuuP+bPpZCkwFFOQVsDFqI3qtnhGNRlgklRdCiIqS4FIIccMLTQ8lLrc4vc+p1FPUcaxDYl5i8b20UDW45OKmlVOpp3h337tsitzErG6z8HHwKbf9mdtmsit2F+/1fI8+9fqUun82/SzPbH2GQJdAYrJjOJtxVr239OxS2vu0x9/RHwAfBx9mdZ/F98e+Z1nGMgDe2/ceO2J3kGvMhYtT+PdtuI8+/sVnkhvMBmJzYgFILUglKS8Ja601ReYifjrxE+mF6erzWri3oK13WwqMBXjbe6sjqp39OrMwdCEAVhormrs3V+sMCxxGZFYkZ9PO8s7edzBjxqyYywyehRDiSiS4FELcEPKN+ayKWEUbrzalTszp6teVB1s9iNFspHe93ui1el7t8iqpBakMCxxmUfae5vdwIuUEhaZCdsbtZGnYUh5t+2iZz1wRvoILWRfYFLUJs2Lm6yNf88KOF2ju3pyE3ASe6vAU/er345cTvxCZFUlkViT2VvZo0dLRryOhqaGsObeGRi6NuL/1/bT3aY+nnSfR2dH8cPwH9Tnro9aX+fzjKcdxtXYluSAZDcVpiay0VhxNPgqAg5UDuYbigFSDhhkdZvDD8R/Yn7Afg9lAG6826khr/4D+9PLvxfaY7YwMGkmQWxBc3Cyk1+mJyIhgZ9xOtBotWrQWm46EEKIyZM2lEOKG8OGBD5l/aj5O1k7svGtnuRtOzqSdISw9jCENh2ClLf136BkhM9h8YTP2Vvb8NPgnWnq0tLifb8xnSegSPjz4IQBDGw4lpyiH8MxwNQE6QFuvthjNRk6knkCv0WNQDADcXu927KzsWHt+rUVZs2JmepvpoMCjf5cOaDVoqOdYjws5FwCw1dpSYC4+1zzIJYiYnBgKTJc/53x4w+GsOr/K4toDrR6gvlN9RjUeRZGpiHOZ52ji1oTXd7/OztidZBVm0cGnAw2cG7AgdAHN3ZsTmhaKn6Mfy0Ysk3yXQiXf36KiZORSCHFDKEnL42HrcdnAMqcoh1xDLpPXTSbfmE9iXiIPtHqgVLkPen9AZGYkjVwbldnW7P2zLc4t3x+/n5SCFOyt7LG3sifPmIdWo0Wn0XEk9QhcHAF0snai0FRIWHoYMTkxFm0eSS4u99iWx7DWFp+mo0WLGbNaRkFRA0tADSwBwjPDy/18dBod7X3as+nCJosAtGSENDQtlNZerdkbv5cpLaewInyFWmZP/B5GB41mXNNxbI/ezum008TnxpNdlM3y8OV8c/QbhgUOI8g1iDsC7ygzjZEQQpSQ4FIIcUOYFjyNLn5d8Hfyt7geciGEg4kHGR44nPs23IdJMWGjsyHfmI+ztTPR2dFsj9nOoAaD1ABVr9VbJFK/1IWsC+QZ8iyupRSkAKjJ1Bs4N6COYx0aOjfkn6R/ADBi5IWOL/DJP59YjG7+lwYNheZCAIvA8lqZFBNv7H3D4jkl+TgBfj/zO7+f+R2ApLwkWnq0pMBYwLnMcygohMSEMDhwcPFRkRpo4NIAHwcf1pxbQ0ZhBn+c+QOzYuZc5jme7fhslfVbCHHzkWlxIcQNKzo7mqHLhgLgbedNUn4SAO/2eJeGLg0J9gxm3KpxnEk7Q2//3sztN7fMdtIK0lh/fj0Gk4FPDn2CSTFxf/D9HE48zInUExSZiwCo61CXtMI0NIqGPFMeWrQ46h3JMmTR2LUxBaYCorOj1XZLNuPUBpdO219qSIMhpBakkpyXTFJ+Eu/1eI/b699OyIUQsg3ZeNt589vp3ziddpqkvCSebP9kmaPB4uYn39+iomTkUghxw3LSO6m7pvW64k08RrORYYHD1HO5fex9OJN2Bh97HwqMBTyz7RkyCzP5uPfH6i7xWbtnERIdYtG20WzEhImXu7zM67tfB8DV1pX1Y9dz/4b72Z+wHzNmdYd3WEYYzdyaWbRRz7Eennae7E/cf8V3KXmP6qLVaKGMoYT1kevRa/QUKcXPnrltJlNbTuW7498BcF/L+5h520zcbN24kHWBVl6tKvXc7THbcbVxpbVX66p5ESFErScjl0KIG1pMdgybIjcxuOFg/Bz9St0vMhURkRFBU/emHE0+yuR1kwGY1XWWmmrn7b1vsyh0kbqecmzjsfwV8RcGs4GGLg1xsXbhSPIRXunyCiODRtJ7UW9yDbl423kzIGCAOt0McJvPbaQXphOREXEdP4XL06LF085THdW9GlYaKxYPX4xZMfPmnjfp4NOBp297+or1NkVt4umtT6PT6Fg5cmXxlLu4Ycn3t6goGbkUQtQaR5KO8MjmRwh0DeTnQT+rJ8iUx9/Jn2mtpl32vrXOmuYexTkdgz2CGR44nIzCDPrW76uWebHTi4xqPIr6TvUpNBVyKPEQS8OWAtDKoxX5xnye6fAME5pOwGA2qCmAkvKT1NN2SuQb82tNYMnFdZ2XCyxberQkNT+VhLwEHPWOjAwaybKwZera0hJGxcjjWx7H0dqR0PRQjqUcY3qb6TjoHcp9to3OBi6mT9Jrr/xrKYS4OcjhsUKIWmNf/D6yDdkcTT5KSn6Kxb1VEav4/tj3GExXv4ZRr9Pzbs93+ar/VxZ5HHVaHS09WuJk7YSnnSctPVri71h8BndSXhKbLmzi438+JrMwsziHZtdX1bqX7py20dnwdIcrj+jVFidTT+JpX7zJyWA28Fi7xzAr/24yauDUgKktpgIQmxtLaHooAK08W10xsATo5d+LBXcs4M8Rf5Y5qiyEuDnJyKUQotaY0HQC8bnxNHJtZBGMnMs8x0s7XwLAxcaF8U3HV/mz8wx52FrZokHDush1aiqhC9kXsNJaYTQbeXbbs3w38Dvsrf7N/ZiQm6D+96Tmk+jk14m5t8/l8ZDH1d3aOnSYMFV5n6tCWGoYeo2eQlMh7+59l0JToXovMjuS38/8zrCGw9TAMiEvgec6PqeWSStIo9BYeNngMdgz+Dq8hRCiNpHgUghRa7jaujKr26xS1z1sPfCx9yG9IP2yKYSuxepzq3l558vUd6pPTE4MPer0wM/Bj8S8RLRo6eDdgX0J+4jLjcNoNlqcrlNSxoyZjZEbKTAV4GrjqgaWtjpbjGZjmZtpaoNC5d9gcuW5lepJQCUMZgOrz6+2uBaTHUNb77Yk5iYy6q9R5BnymBY8jeltpksOTCGEbOgRQtwYCk2FFJoKcbau+j83Sjb0lOSG1Gv1GM1GHmr9ECODRmKlteKzQ58xJGAImYZMdRS1hK3O1iJxuZPeiRxDDgpKqXyT/2WFFUaMVf5O1ampW1Mc9A6cSj1l8d56rZ4fBv5Ae5/2Ndo/UT3k+1tUlASXQoibUlRWFKvPrWZIwyEEugSWWzYlP4WfT/xMfaf6HEs5xq7YXaQWpNLVryvfDfyO/236H7vidsHFk3BMigktWvwc/IjNjaWxa2PCMsIs2tRr9ThbO5NakGpx3VzkQUH8GLS2sdj6rKmGN782l0uJ5Kx3xkHvQHxevMX12+vdzp64PWqQ2dStKUtHLL1u/RXXj3x/i4qSaXEhxE3p1V2vcjjpMLtid/HHHX+UW9bTzlM9dWYCE9gXv48159YwqcUkADIKM9SyJqV47WTnOp3xtPUk9lwsCbkJjG48mjMX7AhLyqbIaS0GDGquzUsZMjpgygvElBeItft2tPrsKn7za3O5XJtZhiyyDFnqz87WznT168q2mG0Wo5c9/XuiKAonUk7Q0KUhjtaO16XfQojaQ3aLCyGu2UcHPqL/kv6EXAipQOnro4lbE4CrWqPZ2a8zb3Z/U23jm/7f4GLjAiXJyIG0/DQCnAMA8HXw5a8z29l3pCVpcd0wpHcDIDU/tVTbVs6H0NrEoXc5gMaqdgWWlZFVlFWcmP6SwLKuQ1261unKV0e/YuLaiUxdP7VG+yiEqBkSXAohrtnis4tJzEtk9bnVFSh9fbzc+WU2jtnIrK6lNwj914mUE3x/7HsOJhws876rrSs6jQ5ATdUTmh5KVFYUy0YU54U0ajLR6NMAE1qb4qljH3ufUm3pbFNwCPwc2zp/UsbA5g2le93uzO45m44+HeFiuqKZW2ey5lzxdP/5rPO8v/99i/RGVSElP4WfTvxUq/KJCiH+JWsuhRDXbOnZpWyO2szj7R+npUfLCteLzIzk8b8fx8fBh7l951Zqp/G68+s4l3mO+4Lvw87K7ip7DukF6fRb0k89A/zZ254lozCDic0n4mHrwReHvyAlP4U+9frwdMjTFimF+vj3IbUgleMpx4svmK1QFD221mYMJgNmqjaoqm3sdHa42LiQkPdvOiYrjRVGxYiN1oZCc/FO9EXDFtHCowUFxoLL/hofSz5GWHoYIxqNuGLy/KdCnmLLhS0EOAewelTt+QvNzU6+v0VFyZpLIcQ1G9tkLGObjK10vZ2xO4nMilT/aeberAK1inNLPr/9eRQU7KzsuC/4PgAUReGxvx9jX/w+2nq15Y3ub1DXsW65bem1euys7DAUFQeXXx39ilxDLvG58UxqPonvj38PF6fZPe09ScxLVOsqKGpg2dqrNfMHz+do8lG+P/Y9O+N2Xv6ZGj0G5eqTwdcW+aZ88vPyLa4ZleKd73qtniJzEV52XhxLOsbLO18mIiOC17q+Vur3Sk5RDvdtuI9CUyFpBWnc3+p+souy1aUIXByt3BO3hz71+lDfqfgYST8HPw4nHaaNVxt1uYIQoubJ/41CiBpzR+AdDAwYyOQWk2nq1rTC9VxsXKjvXB8rrZVFQJpnzGN7zHYKTYXsS9jHj8d/LFXXZDaxKmIVe+L28MjmR5i5fSa/D/2d93q+x7zB82jl2QqAM2lnWB2xmmCPYHwdfAn2DFbXUDpZOwFwIP4AXAwWn2j7BLvjdlPXsS5Hk4+W2/+bIbC8nJJjHnOMxamYkvKTeGf/O4RnhKOgcCDhQOk6Or0aSHraefLI5kfosbAHf5z+dyPWk38/yUs7X+K1Xa8xo8MM/hr5F6n5qUxeN5nPD31+Hd9QCHElMnIphKgxbrZufNzn40rXs7OyY/mI5eSb8i3yXjroHXi1y6t8d+w70grS6OzbuVTdhaELeX//+2pKIYDwxuEMCxwGFzfvzNg6g5DoEMIzwvlh4A909itu58t+XxKWEcaCMwvILsomz1R8BrdBMfDApgeu+nO4WXT168qe+D1l3tNqtAyoP4An2z9JSn4KD216CCuNFd8O+BY3Wzf+HP4nSflJNHFrwocHPwT4d7kBYKcvXvpgr7dHo9EQ6BJItqF4Q1RWUVaZzxRC1AwJLoUQtd63R79lb/xeXuz8orqDW6/Tl7k2r+RoyLf2vsXcI3MZ0GCAxZSpu607XBx9bOnREq1GqwaPXDxnvJ13O0KiQ9Cixd/RX73XrW432vm045eTvwDQ178vITEhZSZJt9HZWByleCvYE7+Hlu4tOZl2stQ9s2Kmb/2+7Inbg4PegbD04rygx1OO08u/F662rrjaugLwaZ9P2Rm7U00FBfDZ7Z9xPOU47b3/TdD+48AfOZR0iIEBA6/L+wkhKkY29AgharUiUxEdfusAwLgm43it62tXrDPnnzn8fOJn7Kzs+G7Ad9R1rIuXvZd6PzQtFA87DzztPOHiueJv7HkDK60Vr3V9DRudDWfTz+Jm42ZRDyApL4l+S/rBxQDypc4v8dOJn4jKigK44ok8txI7nR35puI1mXqtXt00VdexLnUc6hDgEsBLnV664gYeUTvI97eoKBm5FELUatY6a+5udjf74vcxotGICtX5X+v/4WXnRVZhFveuuxdXG1fWjV6nJvRu6v7v+k5FUdgZu5O159cCMKjBIHr596KJWxOis6JJzE3Ex+HflEJajZbBDQazPnI9HnYejAoaxejGozmbdpZ1kevUc8d97H1IzEtEp9HRxa8LYelhJOUnVfGnU7uVBJZcHCkuMBaQZ8wjNieWIlMRX/X/SgJLIW5CElwKIWq9lzq/VIFS/7LX23Nvi3v5/fTvAOQb8zGaLc/vNpqNjPprFNHZ0Tx727M0cG5AfG48nxz8hECXQJLykrhvw31Y66xZcecK6jjW4Z/Ef3hg4wM4WzvzSZ9POJlykvf3v8/QhkNp4dGC6KxodBod/k7+LB2+lMjMSPyd/SkwFnD74tsBCHAKQEHhQvaFKvyEag8HnQP5pnyLNEx+9n7Y6+2JKIignmM9onOiSc5P5pWdr/BA6wcqnCVACHFjkOBSCHHTuqvpXbjZuBHgHKCu5yuRVZRFZFYkXNzkE5sTi8FsIDwznO0x2/Gy98KkmMg35pNdVLxxJCorCqPZSFpBGl8d+Uo9T/yPM3/QxK0JZ9PPAhDsGUxUVhS+Dr446h2x09nhZedFcn4yUdlR1/1zuJ5yTbmlrjVya0R6QTpcTLReYkPUBjZf2Mza0Wup41jnuvZTCFF9JLgUQlSbY8nH+Cv8L8Y1HVcjo1M6rY6hgUPLvBefG6+ujxweOJwlZ5eQkJdAc/fmDG44GHdbd2b3nI2TtZM6jT6i0QjyDHl42Xvx5eEvLdoLTw9X81euObeGNefW4Gbjxvox6zmbfpbk/OTr8s61iQ4dJkzsjN3JwICBdPLtxIIzCyyOjDQrZrQaLbE5sURlRdHVrysajYYiUxG743bT2qs1Oo2Ob499S3P35gxvNLxG30kIcWUSXAohqs0be97gbPpZwjLCmD9k/lW1kZSXxOR1kzErZuYPmY+vg2+V9C0yM1LdeNPdvzt3Nb+L1PxUGro0VMv8NzBVUDiZepK02OJzxc9nnaetZ1tsrGzYl7Cv1DGHGYUZrIxYyZeHv7wlN/qYMGGttabIXMTGqI1o0arT5SWfh4LCqohV/HLqFzILM3G3dcdaa42V1oqYnBiauTejt39vfj31Kxo0dK/bXd3xL4SonSS4FEJUm06+nTibfpZOvp2uuo0zaWeIzSmeSj2ZerLKgsvBDQaTmp+Ki42LemTlpTkz/+uP03+wJ24PW2O2AuCoL94clGXI4tMen/LE30/gbO3MsZRjap37gu9ja8xWMooyqqTPNyKzYlYDzEvXYSooaoD5x5k/KDQWp21KK0izqF9gLKCddzv0Wj1BrkFqAnshRO0lqYiEENWq0FSIjc7mqusbzUa+OPwFZsXME+2fUE+AqWwfPtj/ARmFGbjYuDCmyZhKnYGenJdM3yV9AWjg3ABXG1ciMiLINmTTyrMVPw/+WX3HRWcWcSL1BI+1fQxXW1dC00L59ui37Infg8FsUJO3a9Cg0+jU4xJvRe627qQVpKkjmi08WhCeHk6RuUgt082vG98O/JYiUxFn08/y2q7XKDIVUcexDu/2fFdNJyWqn3x/i4qS4x+FENXqWgJLACutFTM6zOCZ257hRMoJfj31K/nG/ArU/Nf2mO0sPruYjVEbWXJ2CW/vefuKdUxmE6/sfIX7NtxHgamA5u7N0Wl0DGowiF+H/sqMDjNw0DtwPOU4I5b9myJpQrMJ3B98P//b/D+6/tGVw4mHOZpyVM3xaFJM2OnsUFAIdAmkgVMDtLfoH8U+dsUpnqy0xZNoERkRFJmL6OjbkV7+vQBo6Fq8TMFaZ82ac2sIywgjKjuKPfF72By1uQZ7L4S4HJkWF0LcEAqMBTy48UEKTYWkF6TzRPsnrlhny4UtOFs708arDQ2cG5CQm0CBqYCOvh2vWDciM4K/Iv4CYFXEKrIKszApJuafms9j7R6jd73evLn3TQDi8uLYH7+fTn6deHnHy6w8t1Jt59NDn5Yancw35aNFS7Yhm/jc+Aq9f8lmoZuFvZU9Z9LPAFBkLuKtbm8xa88sAPwc/JjVdRaRWZEEuQZxNv0se+L2MDBgICdSTpBRmIGrjSt96vWp4bcQQpRFgkshxA3BSmuFr4MvUVlR+Dv5X7H85qjNzNg6A61Gy/IRy+lXvx8/nvgRnUbH/9r8r8w6BcYCPj74MUazkftb3c/QhkOJy4kjozBDTaFze73ifJXZRdnc7n87ITEhAGg0GsLTwy0CSy4mae/k24n9CfstrpsxVziw5OL55TeTPGOexc+z9szC3sqebEM2R5OOkpKfQmO3xgA8uuVREnITGBY4jF+H/npNz1UUhZ9P/oxZMXNf8H0WR4MKIaqGBJdCiBuCldaKxcMWk5yfTIBzwBXLO1o7okGDtdYaGysbmro3RYOGQNfAy07Vb4raxMLQhQCsPreadWPW4WnnyaHEQ6w5t4bm7s15p8c7LAtbxuu7X6euY13e7fEuXnZedPTtSFRWFFYaK4uRyjuD7sTT3pMDCQcY0mAI6yLX3XK7xivCpJiY0X4Gb+57k6jsKDZGbWRKyylw8bjIhNwEvOy8rtjOleyO280n/3wCQJBrkIx+ClENZEOPEILV51azKmIVD7d5mLbebWu6O1XmbPpZ7KzsqOdUD4DU/FScrJ2w1lkDEJ0Vzf0b78fZ2pmfBv9EZkEmE9ZMILsoG61Gy+qRq6nnXFw3LD2MAmMBcw/PpdBcyD+J/6iBZGvP1vx+R/FpQOcyzrE1ZiufH/ocs2JmasupbI/ZTkRmBI56R3IMORZ9HNJgCNlF2eyK21WhoPPSdD43mxbuLTBjxt7Knvd7vo+fox9c3FA1fvV4souy+W7Ad7TwaIGtle1VPSM+J5671tyFSTGxYOgC9ddXXJl8f4uKkpFLIQQfH/yYlPwUNBoN3/T/pqa7U2WauDWx+NnDzsPi53+S/iE+N5743HjC08Np79OeXXftYn3kejztPNXA47dTvzH7wGzsdHbqedkdfTtiZ2XH9pjt6vR2gbEAR2tH7gu+jzGNxzBu1Th+Pvmz+rwcQw5atDR2a8ydQXdS17Eufev3Zc4/c9gZt7NC73SzBpYAp9JOAXBH4B0WZ45nFmaSkp8CwIytM8gz5PHdwO9o592u0s/wc/Rj89jNKCjqXzKEEFVLgkshBKMbj2Z52HKGB958p5/MOzGPI8lHeLbjs9R1rGtxb2DAQI4mH8XZ2lkdsdVoNAxpOMSiXEJuAlzciFNCg4b3er7HyvCVBHsGYzAbGLdqHJFZkXTy7UShqbDUmkqdRse8wfNo5dkKnVbH3vi97Inbg6326kbhblZrzq0htyiXes71uL3e7XT07cib3d4kNC2U388UjxCfTDl5VcElYBG4CiGqnkyLCyFuWllFWXRf0B2Ae1vcy3Mdn6tU/dT8VNxs3SgwFrAifAVrz6/laPJRuvh14dUur1LfuT4fHviQ+afmM77JeJaFL8NottwZbm9lT5GpSF2HueCOBQR7BqsbjrjktJqyNHRuyPms81f5Cdy4bHW2FJgKcLd1Z9uEber1P8/+SWxOLA+2fhA7KzsA0gvSeWzLYwDM7TcXN1s3uJgjdVfsLpp7NMfb3ruG3uTmId/foqJk5FIIcdNy0jsxMGAgR5KP0K9+v0rV/frI13x19CsGNxjMh70/ZGLziUxoOoGMwgyL6fXDSYcBOJJ0hC/7fsmuuOJpdQcrBzr4dGBsk7EoKNy77l6sNFY46B0AiMuNU9sob61lVFbUVbz5jenSILvk/PGswiymrJvC9wO/x1pnzZgmY0rVO5x0WD0Z6VDSIfXX+vPDn/PziZ+p61iX9WPWX9d3EeJWJsGlEOKmpdFo+LjPx1dV90TqCbh45GQJnVZXat3m611fZ/qm6ZzNOMv2mO280PkFnu34rHp/wZkFGEwGnr/teX4++TPjV43ng94f8Ok/n1aoHzfzGstL6TV6utftTmp+KsdTj6vXjYqRQ0mHiMuJo4FLgzLrdqvTjeGBw0kvSKeB8yVlLsbs/z3zXQhRvWRaXAhxS1EUhTn/zCEmO4aXu7zMnrg97IzdyaNtH6W+c321XEx2DEvOLmFAwACCPYNLtZOcl0xGYQZ/hf/F/FPzUVDQarQcuOeAulFkX/w+Htj4QKm6Per0UDfwOFs7oygKhaZCi2MPS9hb2VvkhCxvCv1m9XDrh/Fx8KGxW2Nae7Uus8zGyI08s+0Z7K3sWTN6DZ52nhjNRnbE7CDILUjNGCCunnx/i4qSkUshxA0vNT8Vd1t3NBoNXEwb9Pmhz+ldrzdjm4y1KHsu8xzzTs4DoLVXaz4/9DlGxYheq+ftHv8eC+nv5M+MDjPKfF5aQRp3/nUn2UXZFtd9HXwtzj73d/LH2dqZAmOBReC4M24n9lb2OOgdyDHklDrOUoOGkUEjsdHZUGQqYlPUJrIN2QQ4BRCVffNNk2s12suOLuo0OnbG7eR4SvFo5rSW03j6tqfh4prKd/a9Q2ZhpnrqksFsUNe97o/fz+eHPyc8I5y7m93NS51fum7vJMStTIJLIcQNrWRDzZjGY5jVrfj4wJ9O/MTWmK3sjNtZKris71SfHnV7cDr1NAcTD9Kjbg/2JeyrVDLtfGM+eYbi0URve29S8lMwK2YCnAPUAHdb9DaWhS1jaMOhRGVFEZMTQ3R2tNpGnjGv1Ck1XAymdBoduUW5LL+w3OJeVHYUHrYepBakVvJTqt3Km7Y2KSY1sAT4+eTPZBRmcFezuyg0FbL07FK4ODX+6e2f4mvvi6+DLwDvH3if85nFm6H2xO2p9vcQQhSTc6+EEDe0Eyml10YOajAID1sPRgWNKlVer9PzSZ9PSCtIY3vMdnwdfNl/z376B/Sv8DPrOtbl+4Hf826Pd1kybAkt3Vtio7PhWPIxvjzyJQCzD8zm7+i/WRi6kD3xe9QdzPUd65fbtkkxUWQuYvOFzWXev1JgqdPocLdxr/C73CguPaZxefhyZoTMoKlbUzr7dqaxa2PmnZzHB/s/wNXWVS03MGAgNjobmrs358l2T/LqrlcZ9dcoTqeerqG3EOLWIMGlEOKGNqvbLKYFT+Pt7v9Oafep14etE7byWtfXyqxja2VL73q9cdA70Mu/V5llDCYDMdkxbI7azKIzizCZTRb3O/p2ZHij4cw7NY/jqccpNBWSa8hlRdgKAIY2HIqj3pG2Xm3xsfehnXc7dBoddZ3qYqX9d9JIQ/FIp7XWGke9o3rdxcYFvebfKfbOvp0r9Hn4OvgyotEItDfZH+8lo5sOVsW77Ru5NsJeb88Pg37ghU4vEJUVRVxuHEeSjqh1Hmv3GIuGLeJc5jme3/E8K8JXEJ4Rzrrz62rsPYS4Fci0uBDihtbQpSFPd3i60vW+6PsF7+17j2e3PUv3ut35qPdH6pQ2wP82/4/9CfvVn+30doxoNIJCUyFPhjxJQk4Cn/X9jDZebdBpdLhau5JWmIZWo+Xdfe8SnxvPqlGr8LTzJDQtlDkH56hTvL8P/Z27V9+NGTMKCnqNnhFBI3iy3ZPkG/NJLUjF39GfI8lHeGbrMxSZi8gqyrLof0uPlpxOO11qSjk2J5Z5p4rXlP73nPPayAorjFS8j74OvmQVZZFZlEm+MR87Kzs6+HRgWvA08gx5pUagU/JTKDQVAjCk4RBS81MZ3Xh0lb+HEOJfElwKIW4oYelh/BX+F8MbDaepe1OLe5mFmXx88GP8nfx5qPVD6vW4nDie3fYsPg4+zO45G71OT1pBGn+c+QOAjVEbySzMtJhSjciIAECv1aMoinq6T2RmJLtidwGwI2YHk1pM4r0e7/HcjuIE7XG5cSw4swAu7kzPKcrhcPJhFEWhl38vxjQeg6edJ8MDh/PXub8A+LTvp3Tx64K1zhpXXNUztXUaHc/e9iw5hhy+OPKF2rcg1yBe7Pwir+16jXOZ56jrUJfY3NhSn1VtDyyBCgeWeo0eg2IgIrP41yU5P5morCiauTdDp9Vd9i8Ynf0682GvD1FQSp28JISoHhJcCiFuKLP2zOJY8jEOJh5k4bCFFveWhy1neXjxJpj+Af0JdAkEICQ6pDjJdgo82OpBmns0x9XGlc6+nTmUdIghDYZYBJYAX/YrTog+MGAgdlZ2+Dj4ANDYrTGTmk8iITeBoYFDATAoBrg4xf1U+6c4lHSI6OxotsVss2jzvuD76ODTgZErRhKRGUFrz9YYzUYe3fIo/k7+rB65Gp1WB8D5zPM8uuXRUmmH6jrU5cXOLzJt/TQM5uLnphemV+lnXNs4WzvjbO1MTE4MXFz2UMehDn9H/Y0GTam/ZCTnJfP01qdp4taEV7u+yuCGg9V7S88uJT43ngdbPYitlRy7KUR1kOBSCHFDae3ZmmPJx8rMd9jZrzPutu74O/lTx6GOen1AwAC2RW9TcyVycYPID4N+uOxzWnq2pKVny1LXtRotz3d63uLa8MDh2FvZ42XvRRuvNtzHfQC8tOMl9ifsZ1rwNJq4NaGDTwe4OCIJ0NStKUvDinc7x2bHcijxEIn5iQxpMAQnayc1VREX12QWmYvoUqcLrTxbodfq1eDSYDZgo7NRp38v1cqzlcVu6xtRoamQuo51icmJwd3GnQdaPcBnhz7jQMIBvj3+LV/0/cJi7exLO1/iSPIRjiQf4c5Gd9Lau/j3yrmMc7yx5w0APO08ubvZ3TX2TkLczCSJuhDihpNRkFFqpLFEWkEaL+14CRcbF97u/jZ6nb7McjXphe0vsOb8GoY0HIKPvQ8rI1YyrOEwfjvzG2bFzAudXiCzIJP5p+aTa8wFoJFLIyIyI7CzsmNAwAAGNRjEC9tfwKSYMJqNZSZgB2jk3IiIrIjr/IaVp0OHCVOZ95q5NyOtIA0/ez9C00PVoyFLjAwayWNtH1NHl789+i1zj8xFr9Xz9/i/cbUp/r2SU5TDuFXjSC1I5fuB39PGq811eLObh3x/i4qSkUshxA3ncoElQMiFEHbFFa+JnNB0Au192lfZc3OKctgVt4ut0VuZ3GIyzT2aV6q+wWwg35hPUn4SXNx84+vgS1pBGiExIermHLPZzNfHvlbradDQ3qc9LtYuHEo+xMqIlaTkp5BtyL7iiT0+Dj5EZEWoI5+1VVmB5YsdX+Q2v9sYs7L4PHFvO2+LDUyetp608mpFcl4y/Zf256HWDzG4wWCauDVh85jNuNq5YqOzUcs7WjuyatQqDGYDdlZ26vWMggwe2vQQRaYivh3wrRqkCiGujgSXQoibSi//XrT3bo+zjTMtPFpUWbsp+SmM/mu0ur5xd9xutBotH/T6QD0dpjwGs4G7Vt9FREYEr3R+he51ujMwYCCv7HpFva9BQ3vv9nT0s2xPQWHJ2SXqz3qtnmCPYHbH7S7eba7V4+vga5GkHcDf0Z/k/GSAWh1YatGWeYZ6RHoEW2O2qj+fSTtD97rdic+Jp5NfJyY0nUADlwYMWz4MgOMpx/n11K/kG/N5pfMrTGg2oVSbVlori1RQXMyRejqtOPfl4aTDFms0hRCVJ8GlEOKm4mXvxS9DfrmmNgxmA1subKGle0vqORefSZ2an6oGlho0pBWkAfD3hb+p51SPd/a9QxO3Jjze7vEy28wz5BGREYFJMRGXG8cT7Z8A4NUur/Jn2J8MbzScINcgrHXWKIrC0x2eJiw9jFXnVpVqS6/V80jbRygyFzHv5Dz8HPwIcgsiKS/JYt1lyQaYirDR2lBoLr1mszr8d7S1rMASYHH4YoufjYpR3SQV7BVMA5cGAHzQ6wM2R21mcIPBTF0/FUA9370iOvl14q6md1FkLqrUSU1CiLLJmkshRK1iMBt4fvvzxObE8lGvj9Tgrjol5CYQmxOrbrj59J9P+fHEj7jbuhMyPkQ9HWbtubWkF6ZzV9O7+P307xxOOszMjjNZEb6Cb45+A8CmsZvU4wf/a0vUFk6mnmRa8DScrJ3KLLMifAVfHfmK+Nx4JjSdwOQWk9GgYfHZxRxMOIiTtRMzb5tJE/cmAMw9NJdvj38LQLB7MCfSTpT7rrU99+WV1ohq0OBh58Hb3d+me93upe7HZMeQkJvAbb63VXNPbz3y/S0qSkYuhRC1SlRmFJuiNgGw5cIWpgZPvab2Sk7WKUnx8195hjzGrRpHRmEGr3Z5lfFNx6vr9Kx11uoJOoCaeghgcsvJTG45GYB+9fux9txaGrs1xsvOq1T7JsWEk7UT/QL60S+gX7n9nb1/trpD/FTqKeo71+f3079jMBv4fuD3OFoXn+KTkJvA/vj9rIxYqdaNy4sDwFZny11N71KTqdvp7Hir+1tsvrCZ6Oxoi6Mya5uydrxzcX3loAaD2J+wn3pO9Wjt2RqD2YBea7lhy9/JH38nf/XnP07/wW+nf2N66+ncGXRntfdfCCHBpRCilgl0DWRck3HE5sRedu1baFooJsV0xTWVMdkx3LP2HrQaLQvuWFDmiKJZMVNkKl6PmG/MB2B6m+m082lHkGuQxak9l9PMvRlrRq8pdT05L5mxq8aSb8znl8G/VGgD0Migkaw+t5oOPh14vN3jRGVF8f7+9wHwsfdhXJNxOFo78uDGB4nMirSom1aQRgPnBnzY+0OauTfjt9O/YVSMWGmt2BK9hfWR66/4/JpkhdVlP++UghQWhS7CqBgJywjj9iW3Y2dlx+jGo3m4zcOlclaazMWnIf1x5g+is6NZHLpYgkshrhMJLoUQtcaCMwvYFbuLpzs8TaBrYJllzqSdYcLqCSiKwq9Dfy03nczZ9LPq2sjwjPAyg0tHa0d+G/obkVmR9KtfPKqo1WjRa/VkFmbiaeeplv3owEf8fvp3nurwFFNaTrni+yTkJqjPP5d5rkLB5fOdnrfIo5lvzKexW2Nis2NZEb6CTw99yshGI9XAMsg1yGI9aL/6/Wjm3ozjyccxKcWjtq90eYV8Y36ZZ2rbWdnRyKURJ1LLn06/Ht7p+Q6FpkJe2335M+FL8pQeTjpMoamQn078hIethzqKXOLdfe+y+OxiHPQOtPFqwwOtHrgu7yCEAG1Nd0AIIbg40vT+/vfZFrONeSfnlVtOURQUFIxmy7WDZsVMSn6K+nNv/9483OZhHm37KN3qdLtsm43dGjMgYIC6tnLLhS1MXT+VcavGEZcTp5bbGrPVYlPJlbTyasUb3d7gmQ7PMLhB+TuQjyUfY87BOZzPPM89a+6h/5L+nMs4h52VHctGLGPLuC1cyLqAWTFzKOkQXJz+Ds8IJ70wXU2tsyJ8RfHnpJjUTTMedh6ERIfAxWMUL5VvzKeNZ83ne2zv1Z4PDnzAJ4c+Ua/pNDp1iYJOoyPHkINZMdPRtyOtPFup5f57Qg/A2YyzAOQacnmk7SOlNuoYzUb2xO0hoyCjGt9KiFuTbOgRQtQar+9+nZ0xO3mz+5tlbtYocTjpMEazsVQKoJnbZrIhcgOPtH2Eh9s8fNX92BS1iae3Po2V1orVo1aTVZiFn4MfZ9LP8Ff4X0xqMYmWHqVP77kWvRb2Ir0wHQcrBzVx+kudX7I4RWbLhS0cTT7KkAZDWBG+gobODXn/wPuYFJO6A9vFxoWdd+0EYNGZRcw+MBsbnQ1Gs7FU8vESNZ0Ds4Fzg1JT/AD1neqjQUNzj+bkG/PZEbtD/YtFE7cmnE0/i5uNGyHjQ0qtqY3KjGL65ul42nnyTf9vcLR2ZFfsLqKzoxnTZAyf/vMp80/Np7FbY5aNWHYd3/bGJd/foqJkWlwIUWu80e2NCpVr592uzOsnU4o3qpxKOXVN/RgQMIDvB36Pm40bu2J38dbet/Cx92Hd6HW09mxNREYEZsWsjnSWJyE3gXf3vUuQa5Cafui/tkRtUTeylASWDnoHhgUO42DCQXQaHWfSzxDkGsTTHZ4mPiee02mnWRG+Qp361qJleKPhPNT6IbVda501BrNBPSbyci4XWF4pQbudzo58U/4VP4OyWOus1bWuUVlRZZa5kH0BACcbJxbcsYDcolz6Le1HriGXs+lnuT/4fsY3HY9Oq2PNuTX8feFvHmn7CI1cGxHgEsD6Mf+uMU3ITeCRLY9gVsyYFTMFxuJAu+TfQoiqI8GlEOKm8WHvD9kctZlxTcddVf0iUxELzyykgUsD9azqv6P/BiCrKAuD2cD0TdM5knyEqS2n8sxtz1yxzeVhywmJDiEkOoSxTcZSx7FOqTJv7n2TPGMeWo1WPYGmmXszTqae5MGND6pBnl6rJ2R8CBsiN3A46bBFG4MDB/NWj7csrg1vNJzjKcfZF7+PPEMeXep0Yd35dWpAqtPo+GnQT0xbPw0z5lLBpF6jp0i5/IhmZQLL7nW608mnE7vjd3Mm9QyZhkxsdbYUmAosnulq7UpGUfFUdUnapP71+6vPyzXkqmX3xO3hwdYPAvDGnjfIN+aj1Wj5sPeHpZ7voHfA3dadtII06jrWZXTj0dzmexttvdpW+B2EEBUjwaUQ4qYR7BlMsGfwVddfcGYBHx38CJ1Gx8axG/G29+aB4AfwtfelqXtT7PX26gad1PzUCrXZt35fVp9bTZBrEN723mWWGRAwgOVhyxneaDhFpiJ6+vekT70+pTbgeNh5YKOzYUCDAWyI3GCxCWdn7E7C0sPUDS9cPI1mfeR6souy4eIJNr4OvsTmxAJQ17Eu7X3aM6vbLFZGrMTD1oOdsTvV0dPyAsvKauvVlvta30cLrxY8uLE4ICxrmr7k2dZaa5bduYwA5wDWn1/P54c+Z2yTsRZnkJ9KO8XR5KN0q9ONIQ2HsCFyg7op67+crJ1YcecKcgw51HWsC8CQhkOq7P2EEP+S4FIIcV0UGAvIN+bjZutW0125rADnALQaLZ52njjqi/NJ6nV6RjUepZb5uv/X7I3fe8XApMBYgK2VLU3dm5aZpuhS9ZzqUWQuIiQ6hM1jN6PX6YlIj+CN3cXLBGa0n8Gi0EUUmApIykuivnN95g+dzz1r7iE0LRQzZjILM/nkn0/4qv9XFm239mytnrV+IesCeq0evVZPU7emnM88z6zds5h520wuZF+gqVtTBjcczIytM9T65Z3c42jliBkzrTxasS9xX6n7l46EFpgKGPznYLRocbF2IbMo0+J+ySilwWzgxU4v0tmvMwHOAaTkp/Dc9udQUNCgwcnGiczCTFp4tKC+U3018f0b3d644rIKFxsXXGxcyi0jhLh2ElwKcRMrMhVhNBux19vXaD/yDHmMXjmahNwE5vabS4+6PWq0P5fTp14fNozZgJO1k8VnpiiKmn+xvnN96jvXL7ed57Y9x7rIdbzQ6QWGNxrOo5sfpdBUyJf9vsTL3qtU+ZK1h0azUT0KcUPUBvW/MwoziM0tHm28c8WdaiqkhcMWklGQQe/FveFiCqX/+qr/VyTmJjLoz0EoKOr6Snsre3KNuaw5twYfex9+OP6DWufS6XlfB1+isv9dExnkGkR4Rjh2OjsKTAUYFSP7EvdZbAoqOSvc286bxPxEfOx92B23Wx0xBZjUbBLTgqeRb8znQvYFfj/9uxoEv7f/PTxsPdg0bhNO1k4EugQSlR1FoGsgOUU5KCgMbzSce5rfc8VfUyHE9SepiIS4SaXmpzLkzyH0WdyH06mna7Qv2UXZxOXEYVJMhKeHX3N7iqJQFYkuDCYDi0MXW6xf9HXwxUHvoP4clRVFn8V9GPznYHVKvDwnU06yO243AAcSDnA8+ThHko9wOu00BxMPllnnvuD7+LTPp/xxxx9q6p1Wnq3QoMFR78jU4Kk81PohXKxdMCpGNkdthotBoLudO4+0eYROvp2Y0WFGqba1Gi1+jn480OoBmrg2YXKLyczqOovnOj1H//r9GRo4lJ9O/IRWo8XLtjjwLQks3Wzc+KTPJ9RxKF4nek/ze1g4bCF3Nb2LfFO+xTGSDZwbEOAUgL2VPa42rgAk5ifS2bczMzrM4HSa5e/BxPxEvB28CXAJoKd/TwYEDLA4DSm1IJWNURux0dmwaPgivhvwHXE5cQS5BfFQq4eY0HSCRXvnM8/X+O9zIUQxGbkU4iaVkJtAUn4SAGEZYRVK4F1dfBx8mNOnOIfjXc3uuqa2kvKSmLhmIkazkd+G/mZx1J/BbGB52HLqO9eni1+XK7b1y6lf+OzQZ1hrrQmZEIKztTORmZG42Lio0/cnUk6oQeWTfz9JfG48c/rMobVX61Lt7YvfxwMbH0Cn0TG04VAea/sYPg4+jAoaRaGpkN7+vcvsh06rK3UsZE//nqwZvQYnvRP7EvZxIOEA9zS/hx2xO/B18CWrKAtn6+J0MA+3vXLapSfaP8GwwGGcTT+Lh50H9667lwbODcg15lJgKsDT1pN5Q+Zxx/I71DqPt3ucxu6NWTRsEVtjtuLv6M/ti27HTl+cU9PPwY8iYxGphalYW1nzeb/PeWDDA+pILBfzaL6w4wX1ZyuNFfe2uFc91jM+Jx69Ts+YJmPoVqcbnx36jM1Rm1FQ+CfhHwbWH8gL219g84XNahtN3Zpipf336+tC1gXGrByD0Wzk+4Hf09mvs3rPYDLw0KaHiMiIYG6/uWX+ugkhqpYEl0LcpFp6tuTVLq+SUZhRKzYu9A/oXyXthKeHk5iXCBdP67k0uCzJ62ilsWLj2I1lTkFfysfeBwA3WzdsdDaEXAjhyZAncbVxZdWoVbjYuDAgYABn0s4AqMndQ6JDygxSSoIqs2LG086TjMIM6jnX483ub5b5/EJTIQcSDtDGqw1O1k6l7tdzqgfAN0e/ITwjnOyibMIzwjmecpxA10CLXJ7Lw5Yz7+Q8PO08GddknMXRmXmGPGbtnsWmqE0YFSMdvDuQb8zndNpp5vadi4OVA3cE3kF95/o8e9uz/H3hb4YFDqOlZ0tG/TWKHEMOCbkJNHRuSLYhm2xDNl/2/ZL2Pu2JyIxgRfgKRgeN5vdTvxOTE6M+197Knt7+vTmWckz9vFeOXKkuOdgbt5fpm6ZjY2XDypEr8XP04/1e73Mw4SDTNkxjydkldPDpwLnMcxa/ZsMbDS/1uZckjf/v2eTJ+cnqiPHuuN0SXApxHUhwKcRNbHzT8TXdBdVf4X8REh3Co20ftdjRfCVHk4/ibeeNn6MfAF3qdOHhNg9jMBtKnbri5+CHBg0uNi7qiTXlGd5oOC09Wqq7sJPyklBQyCrKIt+Yj4uNC9Y6azXlkIuNC6dST5Waki3R078n3w34jkWhi5h/aj7Lwpax++7dlz0v+43db7Dq3Cpu87mNWd1moUFT5nrOu5vdzS8nf2Fis4ksPruYcxnnSuX6/Pnkz5zPPM+5zHMcSjzEgIABamLxkOgQ1kX+u/M8pSCFu5vdTSOXRvSu15ve9f4dUZ3ccrJ6lOKcg3MIz/h3GYOPvQ/d6nbD196XXvWKUzW18WpDG682vLDjBdacs9y45GzjjF7374lAt9e73WIt69dHv8aMmXxjPhmFGerxnA1dGuLn4EeeMY8WHi34ou8XvLHnDbr4dVFTD10qyC2IXwb/Qq4ht1Ty/TqOdZjRYQYRGRG16v8HIW5mElwKIa6Lt/e+TYGpAGutNR/0/qBCdVZFrOKlnS/hZO3EutHrcLFxQavR8kjbR8os3y+gH6tHrcbFxgVHa8cKPePSM8zHNhmLXqenrmPdMs8hr8j51F3rdOVI0hG2XNiCt733ZQNLLklenlmYycgVI0EDC+9YWOo4w/FNx6uB0dgmYzGYDVjrrC3KTGs5jW+OfUNCbgLd63a3OLGmo29Hmro1JTYnlhxDDk3cmvBS55fK7NO5zHPEZsfS078noxuP5mTqSbztvanjWIepLaeWOcLKxeltgCauTQhwDiC1IJVH2j6Cs7Uznnae1HOqx8yOM3l779tsjd7KrG6zaODcgENJh2jg3IBm7s3UtjzsPFg/Zj2Koqjv8eOgH8v93HMNuXx55EtismOY0Mwy+L8v+L5y6wohqpYEl0LcQkIuhBCeEc6kFpMqNLJXlYY1GsamqE0MbDCwwnXyDHlwybTnpYpMRWyI3EBLj5YWAeKVdnJvjNzItphtTG89vVRZnVbH6MajK9y/ywnLCIOLu7/L80a3N+hfvz96rZ6ntj4FCuQYcsqto9FoSgWWAKMaj2JU41Flnhzkbe/N0hFLyTPkcTT56GVPOMoszOSu1XeRb8zn1S6vMr7peD7o9QGD/xxMgakAf0d/RjUeRXh6OG62bnjYeah1Z/eazaaoTURmRrL47GIAbHQ2NPdoTsj44nPNFUXhz7N/YlSMrD23lje6vcGghoMI9gzGYDJgpbVSg3GtRguXicuLTEX8eOJHvOy8GNtkLADzT83neMpx4nPjSwWXQojrS4JLIW4RKfkpPLX1KcyKGY1GU6FRuKr0etfXeb3r65WqM77peDzsPKjrWBd3W3eLe18d+YofT/yIs7Uz2yZss9jgUZ5Xdr1CvjEfk2Li/Z7vV6o/FVWSpPvS9aBlcdA7qGsjv+z3JYCat/FqlXckpb3enq51ugKwNXor8bnxjGsyzuKzK6n/3r73iMmJYVrLaeo1K60Va8+t5fkdz+Ni7cKyO5epieF9HXzxsvPigwPFo9LO1s74OfhZPF+j0fBcp+fYFrONKS2noNfp6VanGwcSDvDw5ocJcA6w2DF/OSsjVvLVkeJ8nq08W9HUvSl3Nb2LpLwkRgaNvMpPTghRVSQVkRC3CEe9I/6O/ug0OoJcg2q6OxYyCzPZHbe71BnYGo2G/gH9y9zpXjI966h3vGxAtTt2N49ueZSdsTvVa33q9cFKa0XPuj0r1Levj3xN1z+68vvp3yv8Pk93eJrlI5bzRd8vLK7nG/M5n3m+zDq9/HupR05WREJuAmNWjmHKuinkFJU/2vlf0dnRPPH3E7y7712WhS1Tr7vYuLB42GKaujXFqBjZGLkRN1s3vuj7Bc90eIbhjYbzzdFvAMgsymTquqkWKaFKRsNtdDb8NvQ3fBx81HuKohCXE8ddTe/im/7fWEz9H046TKGpkLPpZ0kvSCfXkMsHBz5g0ZlFZfa/mXsz7Kzs8HPwU5cv3F7/dpbfuZwpLadU6rMQQlQ9GbkU4hZha2XLnyP+JNeQazGdWRtM3zSdk6knmdB0Aq90eaVCde5vdT8dfDqop+qU5dNDn3I67TQJuQlq4vYPen1gkRT9StaeX0uOIYcNkRsqnLRbo9EQ5FY6gJ+0dhJn08/yfMfnmdRiUoXaupx98fs4m34WgNNpp+no27HMciazyWL9JRdHFd1s3cgszCw1ulrfuT6vd32dX079wohGI8gz5DFj6wyyirKw1lmTWvDvsZepBakYFSN6TfGmnd71erPwjoU4WztTz7meRbuv736d5eHLCXINYtGwRRZT+3c3u5uMwgyCXIPwdfBl3ol5/HrqVwC61elWqq1gz2C2jt+KXqtn7fm1LD27lOltptfa5PxC3GokuBTiFmJrZYutlW1Nd6OUkvQx/00jcyVtvduWe394o+Ek5CYwotEIi+sVDSwBnuv4HH+G/cnkFpOvWHZX7C5MiqnMEUizYiYuJw6AiIwI7l17b7mn9lxJv/r92BW7CwdrB9p5t+NY8jFCokMIdAnkaPJRpracym+nf+OP03/wZPsnub/V/WpdFxsXlg5fyif/fML++P108etiEaC38mrFR70/govrXkuSqhvNRj67/TNCokNwsnaiR90e6LV6i3619GxZZn9LRo/DM8JZe24tIxv/O33tZO3Ecx2fU39u690WOys76jrWxdPes8z2Snadf330a2JzYvnh+A8SXApRS2iUqjjmohKysrJwcXEhMzMTZ2fn6/loIUQtlZibyKGkQ/Sp1+eyG40yCjJYHr6cLn5dajQh/OUcSTrCvevuBeDHgT/Sya9TqTLHk49zKOkQ/o7+xRt4gNk9ZzM0cOg1P7//kv4k5iWqZ3QPaTCE02mnicyKpKNvR34a9JNF+ZKd+AA/DfrpsiOfXDz9Jiorit7+vSsVmF+q02+dyDflAzC+yXhe7fpqueUNZgNWGqsrPm/+yfksCl3EY+0eqxX5XG9m8v0tKkrWXAohapyPgw9DGg4pdwf7hwc/ZM4/c3hkS9lpiMqz9OxSnt/+vDpyWNVWn1vNy7teRqfRodPocLB2KLNcC48WuNi4YGtly6igUQxpOKRUrs7ymBUzr+16jXvX3suFrAul2gYIcA5Ag4aOfh15vuPzDG0wlMENBnPPmntYcGaBWr6td1v8HPwIcg2iiVsTuBj8XromNDYnlnf2vkNsTix96vWpcGAZnxPPyBUjuWftPWQXZQMwpOEQNGjQarT09L/yele9Vl+h501uOZk1o9dIYClELSLT4kKIG0J9p+K0QSWn1lRUgbGAN/e8iYKCs7UzL3d5ucr7Nv/kfC5kXaC+U33m9ptLQ5eGZZZbFr6MN/e8iU6jY+3otTjoHXgy5EmsddZ82OtDiwTjZYnJjmF5+HIA1p1fx/Q209V7n/T5hPTCdDztPCk0FWI0Gxn11yiS8pK4kH2BE6kniM6O5u5md8PFz3Hj2I1q/e0x23l0y6NYa61ZNWoVdRzr8Nmhz1h3fh1rzq1h98TdFf489ifsJyIzAi6eotTRtyNvdH+DWd1mUWQuuuJucCHEjU2CSyHEDWF6m+kMbDCQOo51KlXP1sqW2+vdzoGEAxYn0VSlacHT+PnEz9zb4t7LBpYA3nbeaNDgoHdgc9RmkvOT2Ru/F4BDSYeuuGYwuygbvUaPTqsrNeKp0+rwtCten2ijsyEtP42E3AQUFJq6NyXPmFdq7emlSo6uNCkmNT9nR9+OrD+/vtwp87L0D+jPvvh9OFo7WuTU1Gg0ElgKcQuQNZdCCHEdRWZGsituF+/vfx8NGtp6t8XNxo33er53xZHLP07/wXv73wPgzxF/qtPZl5q9fzaHkg7xRrc3iM2OJSo7iknNJ5WZeP2/dsbuxNXGlWDPYIxmI1lFWdhb2dfKTWDi+pPvb1FRMnIphLhlFRgLWHp2KS08WtDep/11eWYDlwYk5SWh0+hwtHZkTp856ojjlYwMGkl0djQ+9j7kGfLYG7+XLn5d1PvZRdn8dvo3uHiW+/Odnq9U3y4dOZ26fipHk4/yRrc3quTUIiHErUM29Aghbhi743Zz34b7WHturcX17KJsTqeeprITMT8c/4HZB2bz0KaH1KMmr4dOfp1YO3otq0auqnBgycX0O893ep6udboyed1kHtz4ILtj/10L6WTtxJQWU2jt2bpSJ9WkF6SzPWY7BlNxEnuT2URoWigAp1JPlVt3VcQqvjn6TaXTSAkhbl4ycimEuGF8e/RbDiUdIi4nTk3fk5yXzKA/B2EwG7ij4R2836viRzqWHNPobe+NXqcvdX995Hpm75/NiEYjmNFhRhW+CZVeO3opvVaPTqPDpJhK9Xtmx5mVbu+hTQ9xJu2MmiJIp9Xxed/POZBwgHtb3HvZeheyLqjpjJysnSqcZF4IcXOT4FIIccMY22Qs8bnxTGg6Qb0WnxuvHhsZlhFWqfZGNR5Fe5/2eNp5lkoGDrA6YjUp+SksD1te5cHltQh0DWTJ8CUUmgtp6VF20nKAZWHLiMuJ44FWD5S7brJkA09sTiyT101mXJNxDG80XD2H/HJOpp7EQe9AkamIpm5Nyy17vWUUZPDD8R9o7dWagQ0G1nR3hLilSHAphChXZGYkL+54kUDXQN7q/tZlj1q8HoY3Gs7wRsMtrrX2as0jbR7haPJRXu1SfmLusgQ4B1z23kOtH8KoGLmj4R1X1d/qVNbxkpc6l3mO13e/DoC7rTsTm0+8bNnvBnzHkeQj/HrqVw4nHSYpL6nU51xi/sn5HE46zBPtn+ClnS9hNBsZ12Qct/nedo1vVLV+OvETv5z6Ba1GS7c63XC0dqzpLglxy5DgUghRrnWR6ziReoITqSf4X5v/VTrPJBeTf2cXZeNi41ItfXy47cPV0m5rr9Z80/+bamm7LBsiNzD38FzGNRnH5JZXPm6yPF52Xvg7+pOSn6ImWC9LvjGfd/e9i8Fs4M5Gd5Kan2oxMnypnKIcPjz4IVxMfN/SoyXHU45zm0/1BZZGs5FlYcvwc/CrUPL1Em2826A7paO5e/Nyk/MLIaqeBJdCiHINCxzGvvh9NHJphL+j/1W1MSNkBn9H/83M22YypeWUcssWGAtQUCwCgiNJR9gQuYHxTcazNGwpecY8nuv43A0bNKQXpDN1/VTyjfn09u+Ni40Lj7R9hEWhi4jMiuS3079ZBJeHEg/xyT+fMLjh4Aqta8wuymbp2aW81f0tWnq2LPdz+ifxHzZf2AwXf63XjF5z2bIOegcGBgzkcNJh+tXvx7O3PUueMQ8na6dKfwYVtSxsGW/tLR4xXzt6rbpO9kr61e/H7rt3Y2tlW6Oj7ULciiS4FEKUq55TPeYNnndNbRxLOQbA0eSj5ZaLyY7hrjV3YVbMLLxjIfWdi0/leWHHC8TmxHI0+SjHU44D0N67/WWnbmu70PRQzmWeA2Bh6EK4eBzjtJbTKDQVMrbxWIvy80/N50jyEcIywioUXH526DMWhS7C3sqe3XeXf7JOO+92dK/bHYPJYJHWqCwajYaP+3xsca06A0subnzSarS42rjiqK/c1PaV8oYKIaqHBJdCiGr3ce+PCYkOYWKzy6/74+Lu48zCTAAisyLV4LKddztic2Lp4teFIlMRuYbc65aXsjp09OnItOBppOansjt2NzqtjsaujfFx8Clz6ndM4zFEZERU+PzsHbE7ALC3sken1ZVb1kHvcF2n/iurR90erB29Fke9Y7UtqxBCVC05oUcIUWsoisKvp35FQWFyi8loNBr1erYhG2friv+ZkWvIZd7JeTRza0a/gH7V2OvK2x23m2Vhy5jUfBJtvNrAxVHByzmXcY7TaacZ2GBgmbva/6v3ot6kFaRxZ6M7ebvH21Xad3Hrku9vUVEycimEuO4OJByAi2dXX0qj0ZS5kUWj0VQqsOTiruZvjn6DBg0Ptn6QR9s+WmvW3r2//33OZ54nMTeRX4f+Wm5Zg9nApHWTyC7KJjIrkkfbPnrF9r8d8C374vdVKpG6EEJUldrxJ60Q4oZztZMeR5OPct+G+7hvw30cTjpc5f0q0dyjOVqNFgWF7459x/aY7dX2rMoaEDAAW50t/QP6X7GsFi32VsVrBx2sHK5YPqMgg7/C/8LLzkumkYUQNUKCSyFuAj8e/5E5/8y5bkfw/ZP4D51/78yI5SPIN+ZXqq5eq0er0aLVaLHWWldbH/vU68PiYYtxs3HD1caVRi6NqqztvfF7GbtyLN8f+/6q6j/e7nEOTDqg7pzPLMzknjX3MH7VeFLyUyzK6rQ6Fg5byC+Df2Fq8NQrtv3zyZ/57fRvvLjzRbKKsq6qf5VlNBvZHrOdpLyk6/I8IUTtJtPiQtzgjiUf49NDnwIQ6BJYrVOhCbkJvLP3HTILM8k35XM+6zwfHfiIV7tWPHl5C48WLBm+BEVRaOpefKpLniGPIlMRrrauVdrfpu5N2TJuCwoK1rqqC2QXhy4mND2UC9kXeLD1g9fc3tHko+qO+kOJh0qdKONp51nhM8jbe7dnvnY+zdya8dXhr0jMS+TVrq/ibut+zf28nC8Of8FPJ36irmNd1o9ZX23PEULcGCS4FOIGV9+pPv6O/mQVZRHsEVytz1oevpytMVsBsNXZUmAqKPdYwctp4tZE/e+MggzGrBxDemE6Pwz8gfY+7ckszORQ4iG61Olyzbksyzoz/Frd3exuYnNiGRhQNccKdvHrwsigkRjNxkolCi9L73q92TtxL5GZkYxdVZzSqL1P+3LPCL9WJUskCowFGMyGCm06EkLcvGRaXIgbnKutK2tHr2X7hO1XPBLwWvWt15cGzg3oV78ff935F5/d/hlPdXjqmtpMK0gjKT8Jg9lARGYEAE+GPMkTIU/w2q7XqqjnVaujb0cWDVvE/a3ur5L2rHXWvNX9Leo61mXosqGsilhV6TayirLYH78fk9mEjc6Ghi4N6VanGw1dGtKz7rUFrFfyePvHaeHegtSCVGbtnlWtzxJC1H4ycinETUCj0aDTlJ/PsCo0dW/KqlH/Bj51nOpcc5uBroHM7jmbxLxEdUq/ZCRM4bpmSqtxy8KWkZKfwupzqyudIP6BDQ9wOu00k5pP4vlOz2Ots+bbAd9WW18vpdfqsdIWf52U5CkVQty6JLgUQtS4oYFDLX7+vO/nHEw4SNc6XWusTzVh5m0zWXluJdNbT6903ZKNVbmG3Gro2ZXN6TOHHbE76Fe/duUUFUJcf5JEXQhRaYqisCJ8BS42LvSt37emuyOAuJw4/kn8h371+8mxh6JayPe3qCgZuRRCVNqGyA28trt4PeTS4UvVXd+3mqyiLBz1jrUiOXsdxzrUcbz2ZQpCCHGtav5PRCHEDeWtPW/x+eHP0aG7pc97/v3073Rf0J0n/n5CvZaQm8Cvp34lLieu3LrHko8x5+AcYrJjrkNPhRDi+pKRSyEEAGbFfMURuMzCTBafXQzA6KDRPNH+CTzsPK5TD2uXo8lHATiSfES99vLOl9mfsJ+NkRvLPdbxue3PEZsTy7nMc8ztN/e69LcmJOYm4mTtJNP0QtxiZORSCMGs3bNo92s7/jj9R7nlXGxcmBY8jfbe7ZnccnKVB5bRWdH8fOJnEnITqrTd6vB0h6eZ0mIKc3rPUa952XsB4OfgV27dYM/ifKStPFtVcy9rTsiFEAb+OZBRf40iz5BX090RQlxHMnIphGBHzA7MipldcbuY2HxiuWWf7vB0tfXjue3PcSL1BLvidvHDwB+q7TlVwdfBl5kdZ6o/r4xYydpza2nr1ZZ3er5Tbt0Pe33Iy51fxs3W7Tr0tGZEZkViVswk5iWSZ8yT0UshbiESXAoheLP7m6yPXM+UFlNqtB/+Tv6cSD2Bv6M/e+P38v2x7xnTeEypVEXVJaMgg5nbZqLX6fmo90c46B0qXPefxH9QUDiddhorTfl/tGo0mps6sASY2HwiWo2Whi4NK3x0pRDi5iCpiIQQtYbRbCQyM5KGLg15YOMDHEw8eF3Pq15/fj3Pbn8WgK/7f02Puj0qXDcxN5GfTvxEtzrd6F2vdzX2UoiaId/foqJk5FIIUWtYaa3UIyxHNx5NTE4M45qMu27P71a3G738e2Gjs6GDT4dK1fVx8OHFzi9WW9+EEOJGISOXQoirkpKfAiBTnkLcIuT7W1SU7BYXQlTa+czzDF02lCF/DiEsPaymuyOEEKIWkeBSCFFpyXnJ5BvzKTAVkJSXVNPdqbTkvGSGLx/OiBUj1BFYIYQQVUPWXApxA0nOS+bedfdiVszMHzIfXwffanuWwWRgefhymrg1oa13W4t7nfw68V7P9zArZrrX7V5tfagux1KOEZkVCcDJlJOyAUcIIaqQBJdC3EBOp50mNicWgJOpJ6s1uPz55M98cfgLbHQ2bB2/FUdrR4v7wwKHVduzq1uvur24u9ndaDVautXtVtPdEUKIm4oEl0LcQLrV6ca04GmYzWZ6+feq1md523sD4Gbrhl6nv6a28gx5fHH4C3wdfJnSsmZzaQLodXpe6vxSTXdDCCFuShJcCnEDsdJaVesJOZcaGTSSYI9gvOy9sNHZXFNbK8JX8Nvp3wDoWqcrTdyaVFEvhRBC1DayoUcIcVlBbkG42LhcczttvdvipHci0CWQuo51q6RvomJS8lPYFbsLk9lU010RQtwiZORSiBvInrg9uNq40tyjeU13pVJaeLRg59070Wqu399nDyYcZMuFLUxsPpF6TvWu23Nrm0lrJxGbE8sDrR7gyfZP1nR3RHUzGWHf1+DoA63H13RvxC1KgkshbhCbozYzY+sMrLRWrBy58oYLmK5nYAnw/I7nScpLIjYnls/7fn5dn12bGM1GuLj7X9wCji2Eja8U/7dPy+J/hLjOZFpciBtESXCmRYtOo7uuzw5PD+fu1Xfz1p63uM6Hel21Dt7FxzdW9hjHm838IfOZ02eOjFreKjybgM4GHLyLRy+FqAFy/KMQN5DDSYdxtnamkWuj6/rcT//5lB9P/AjA3+P+xsve67o+/1JZRVn8dPwngj2D6R/QX71uMpv4/PDnGM1Gnmr/FHqdnlxDLg56hxrra22xKmIVaQVp3NP8Hqy05U9YHUs+xt8X/mZC0wn4Ofpdtz6KKpSXBlY2YF21v/fl+1tUlEyLC3EDaefdrkaeOyJoBP8k/kMLjxY1GlgC/HziZ3488SM6jY4dd+3AydoJgD3xe/jpxE9wcQPRgIABElgCoWmhvLSzOO2Sq40rdwbdWW75mdtmEp8bz/nM83zW97Pr1EtRpezda7oH4hYnwaUQ4ooCXQL5deivNd0NAII9g9FpdDR2a4ydlZ16vZl7M+o51cNgNhDsEVyjfaxNPO08cbd1J6coh4YuDa9YvqVHS+Jz42npKWv1hBBXR6bFhRA3nOyibOys7K44xSuK5RnyMJgNFUorZVbMpBWk4WnneV36Jm4c1+v722RW2H8+jaTsArydbOnU0B2dVlNtz5s2bRr33nsvYWFh7N27F0dHRx566CFatWplUW7Hjh0sXLiQjIwMWrRowSOPPIKbm5t632w28+OPP/L333/j5+fHXXfdxYIFC+jbty/Dhw9Xy+3atYuFCxeSnp5O9+7deeihh9Dp/l1H/+eff7J69WoKCwvp3r07DzzwADY2/+Yarmh/a5Js6BFC3HCcrJ0ksKwEe719hfOVajVaCSxFjVl/Ip4es//m7u/38uTCI9z9/V56zP6b9Sfiq+2Z69atY/z48Rw+fJj+/fuTnZ1N165diYuLU8t89dVXjBs3jvr163PHHXdw8uRJ2rRpQ3p6ulpm5syZvPDCC3Tu3JnGjRtz11138f333xMREaGW+fzzz+nXrx+2trYMGjSI48ePM2PGDPX+rFmzuO+++2jevDm9e/fmk08+YeTIkZXub41TrrPMzEwFUDIzM6/3o4UQQghxlar7+3vd8TilwfOrlYD//NPg4j/rjsdVy3N9fHyUKVOmqD+bTCalXr16ytdff60oiqIkJycrtra2yq5duyzqdevWTXnzzTcVRVGU2NhYRa/XK6tXr1bvHzhwQAGUTz75RFEURUlISFCsra2VH374waKd1NRURVEUJT4+XrGxsVEWL16s3gsLC1N0Op2ybt26Cve3NpCRSyFEldsdt5tXd71KWHpYTXdFCHEDMJkV3lh1irLW6ZVce2PVKUzm6lnJ17NnT/W/tVotDRs2VEcCt23bRmFhIXPmzGHs2LGMGTOG0aNHExcXx/HjxwH4559/UBSFwYMHq+3cdttt+Pr6qj9v374dg8HAxIkTLZ7t7u6utlFUVMSdd/676S4oKIg2bdqwe/fuCve3NpB5JSFEldoavZXH/34cLh49+HX/r2u6S0KIWm7/+TTiMwsue18B4jML2H8+ja6NPKr8+ZeuaQTQaDSYzWYAMjIy0Ov1TJo0yaLM5MmT1eAxPT0dR0dHi7WTAC4u/y5HycrKws7ODjs7O8qSmpqKvb091tbWFtfd3NxITU2tcH9rAwkuhRBVakfMDvW/u/h2qdG+3Kpic2Kx0ljh4yBJtMWNISn78oHl1ZSrSgEBARQVFdGiRQuaNGly2TIZGRmkp6erm3wMBgOxsbFqmQYNGpCXl0d0dDT16pU+Ya1hw4bk5uaSmJiIj8+//+9GREQwaNCganm36iLT4kKIKnV/q/sZHjicN7q9wZTgKTXdnVvO8eTjDFs2jOErhhOdFV3T3RGiQrydbKu0XFXq3bs3QUFBzJgxg7y8PPX64cOH2bVrFwBdu3alXr16zJ49W73/5ZdfkpOTo/7cq1cvAgMDef755zEYio9jzc/PZ9WqVQB07tyZRo0a8fbbb6t15s+fT1xcHGPHjr0u71pVZORSCFGl6jjW4d2e79Z0N67ZucxzaNBUKDdkbZKSn4JRMWI0Gskqyqrp7ghRIZ0auuPnYktCZkGZ6y41gK9LcVqi602v17Nq1SrGjRtHw4YNad26NbGxsTg6OvLTT8UHN1hbW/Ptt98yfvx4NmzYgLW1NSaTibp166pT5Xq9nhUrVjB69GgCAwNp0qQJ4eHhvPfee2ob8+bNY8yYMYSEhODs7MzRo0eZO3cuDRveWH8OSZ5LIYT4j1Opp5i4ZiIaNCwctpCm7k1rukuVsvbcWuz19vSp16emuyJuItX9/b3+RDwP/3YILtnEw8XAEuDrSe0ZHFz1R5KuX7+e4OBg/P391Ws7d+7Ew8OD5s2bq9cUReHYsWPExcWpwaFGY5l/My0tjQMHDuDr60uLFi3w8fHh66+/ZsKECWoZk8nEP//8Q3Z2Nu3atVM39JTIz89n3759FBUV0aFDBzw8LNeYVrS/NUmCSyFuYRsiNxASHcL01tNvuBG66nQw4SDTNkwD4Nchv9LWu21Nd0mIGnc9vr/Xn4jnjVWnLDb3+LnY8vrwFtUSWFalXbt20aJFC3XN5dy5c3nmmWe4cOGCxRrKW4FMiwtxC3tt12vkGfMwm8180PuDmu5OheyP349Go6Gjb8dqe8ZtvrfxTf9v0KCRwFKI62hwsB8DWvhe1xN6qoqiKLRv354GDRqQlpZGZGQk8+bNu+UCSyS4FOLW1j+gPxsiN9CrXq+a7kqFHE46zP0b74frMKLYvW73amtbCHF5Oq2mWtINVbcePXpw4sQJDh8+jE6no2XLlrfsDK0El0Lcwt7p8Q7v9HgHgM8PfU5sTiwvdHoBN1u3K9atCdZaa7QaLRo02OhsKlBDCCGuHwcHB3r06FHT3ahxlQ4u8/PziY6OplGjRqWShUZHR1NUVESjRo2qso9C3PCyi7L5K/wvOvp2rJWbQyIyIvj++PcAtPBowZSWtTOFUEvPliwethitRktjt8Y13Z1a7UTKCTQaDS09WtZ0V4QQt5hK5bn85ZdfcHNzo2nTpnh7e/Prr79a3F+0aBFffy2ncQjxXx8f/JjZB2bz0KaHarorZfJ38qejb0fqONShW51uNd2dcjV1b3rZwPKH4z9w/4b7OZt+9rr3qzY5mnyUiWsmMnHNRE6mnqzp7gghbjEVDi5jY2N5+OGHmTVrFrt27WLatGlMmzaNd9+98fPZCVHd6jjWAcDPoXbudrTR2fDToJ/YMHZDtY8IKorCO3vfYer6qVWa5NtgMvDZoc/Yn7CfBWcWVFm7NyKzYkZBQVGK/xFCiOupwtPi27Zto3///rzwwgsAdOvWjb59+zJ27FjMZjOvvPJKdfZTiBvaQ60f4vZ6t+Pv5F+B0je3+Nx4FoYuBGDluZU82vbRKmlXr9MzoekE9sbvZVjgsKtqIyIjgoc3P4yPvQ/fDfwOO6uyzwCuLuvPr+d4ynEebPUgrrauV91OO+92zB8yH61GS7BncJX2UQghrqTCwWVOTg5+fpajLkOHDmXt2rXccccdKIpy2cPYhRDckmsEFUXhi8NfEJ8bzwudXsDFxgVfB1/ubHQnYRlhDGkwpEqf90qXa/tL7t74vcTnxhOfG09UVhTN3JtVWd+uJKsoi+d3PI9ZMWOltWJGhxnX1F4773ZV1jchhKiMCgeXHTp04Pvvvy91vU+fPmqAWa9ePe64446q7qMQooacSTvDc9ufo5lbM97v9T5aTaWWaXM2/azFRqF7W9yLVqPl7R5vX7FuTRgWOIxTqafwsfehqdv13Xhlb2VPC/cWnEk7I4GhEOKGVqng0tramj179tC1a1eLe71792bdunUMHTq0OvoohKghGyM3cj7zPOczz/PMbc/g41C5ZMABzgG0925PfG48Xf26VqBGzXKxcVFTM11vVlor/rjjDwpNhdha2dZIH4QQoipUKhXRypUrL7s4vGfPnpw4caLUOZtCiBvX6MajOZV2imZuzSodWALYWtnyy5BfqqVvNyONRiOBpRDihidniwtxHZxMPcni0MWMDBopU55CiBvSdfv+NpsgajfkJIKjDwR0A62uAhWvzvz58+natSuKonDo0CEcHR0ZNGgQer3eolxSUhIhISEUFhbStWtXGje2XEdvNpvZtm0b0dHRBAUF0bVrV4sBt7y8PEJCQkhPT6dLly4EBQVVqv2Sfmq1Wg4ePIijoyP9+vXD1rb2/YW0UsHlW2+9RdeuXenfv3+pe4cOHeLXX3/lk08+KbcNCS7FrejetfdyJPkIQa5BLL9zeU13RwghKu26fH+fWgnrn4esuH+vOdeBwbOhxYhqeaSvry+NGjUiPj6eLl26sHv3bry8vAgJCcHR0RGANWvWMH78eDp37oyLiwvr16/ntdde48UXX4SLB8z06dOHzMxMunTpwvnz57GysmLjxo3odDr27NnDqFGj8PHxITg4mEOHDvH444/zyCOPVKj9kn42btyY1NRUOnTowN69e7G3t2f//v3Y2NSyE8uUCoqPj1cCAwOVwsLCMu+bzWaldevWyqlTp8ptJzMzUwGUzMzMij5aiBveV0e+UoLnBSsf7v+wprsihBBXpdq/v0/+pSivuyjK687/+cel+J+Tf1XLY318fJTAwEAlIyNDURRFSU9PVwICApQ333xTURRFyc3NVerUqaO88sorap2//vpL0el0ysmTJxVFUZS1a9cqrq6uSn5+vlpm27ZtisFgUPLz8xU/Pz/lf//7n2I2mxVFURSj0ahs3769wu2X9LN79+5qHJaZmam4uLgof/zxR7V8Lteiwls/9+3bR9u2bbG2ti7zvkajoXv37uzYsaMqY18hbgoPt3mYw/ceZmbHmTXdFSGEqH3MpuIRS8qaTL14bf0LxeWqwb333ouLiwsArq6uTJkyhRUrVsDF+CcuLo4ZM/5NDzZixAgaNmzIypUr1Tp5eXns3LlT3ZvSq1cvrKys2Lp1K/Hx8bzxxhvqNLlOp6Nnz54Vbr/E3XffrcZhzs7OtGjRgvDw8Gr5TK5FhYPLCxcuULdu3XLL+Pv7c+HCharolxA3HSttpfbP1YgCYwFT1k1h0NJBhKfXvj+whBA3qajdllPhpSiQFVtcrhrUq1fP4uf69esTHV18glh0dDSOjo64u7tblAkICFBjnq5du/LOO+8wdepUfH19ufvuu9m6dStcPOHQwcEBb2/vMp9dkfZLlATAJaytrSksLLzq964uFQ4u3dzcrhg4RkVF4ebmVhX9EkLUgKisKA4lHSIuN4498XtqujtCiFtFTmLVlquk1NRUi59TUlLUYNDb25vc3Fzy8/MtyiQnJ+Pj828WjZkzZxITE8P27dtp2rQp/fv35+DBg7i7u5Obm0tOTk6Zz65o+zeSCgeXPXv2ZMOGDYSGhpZ5Pz4+nsWLF9OnT5+q7J8Q4jpq4taEB1s9yLDAYYxoVD2L54UQohTHCgZRFS1XSUuWLMFsNgNgMpks4pmOHTvi4ODAwoUL1fLHjx/n+PHjapkLFy6oI4hNmzZl1qxZuLu7c/z4cXr37o29vX2pg2hKBuwq0v6NpsLzdAEBAdxzzz307NmT119/ncGDB1O3bl0SEhLYunUrr732Gr1796ZDhw7V22MhRLXRaDQ80f6Jmu6GEOJWE9CteFd4Vvxl1l1qiu8HdKuWxyckJDBw4EAGDhzIhg0biImJUXdqe3h48O677/Lwww9z4sQJXFxc+Prrr5k4cSK9e/cG4OTJkwwYMIA77riD+vXrs337drRaLYMGDcLd3Z2vv/6aBx54gCNHjhAcHMzevXtp3Lgx77//foXav9FUKhVRQUEBjz76KD///HOpZOpjx47lp59+wsnJqdw2JBWREEIIceOp9u/vUyth8eSLP1waY1zMFTl+frWkI/L19eXDDz/E1tZWzR95//33U6dOHYtyu3btYtWqVRQVFdG9e3dGjx5tkccyOjqaJUuWEBMTQ4MGDZg4cSKenp7/vt6pUyxdupTs7Gx69OjBnXfeWan2n3/+ecaMGUOnTp3Ua3PmzCEoKIgRI2rXTNNVJVE/f/4827dvJyUlBVdXV3r06EHTphU7h1eCSyGur8zCTLQaLU7W5f/FTwghylNzeS7rwuD3qzXP5UcffcSkSZOqpf1bUaW2r6anp/Pxxx9z5swZGjRowNNPP10qshdC1B5h6WHcs/YerLRWLB2+lDqO8v+rEKIWazECmt1xXU/oEVWvwsGlyWSiR48exMbG0qJFC7Zu3cqCBQs4e/YsDg4O1dtLIcRVuZB9gXxj8Q7E+Nx4CS6FELWfVgcNe163x02ZMoUmTZpct+fdCiocXG7evJmcnBwiIiLw8PAgOzub7t27s3TpUqZMmVK9vRRCXJW+9fryYqcXsbWypYOPbLYTQoj/mj17dk134aZT4eAyIiKCoUOH4uHhAYCTkxNjxowhIiKiOvsnhLgGGo2Gic0n1nQ3hBBC3EIqnOcyLy+v1PS3g4MDeXl51dEvIYQQQghxA6rUhp7du3cza9Ysi59zcnIsrnXr1o2BAwdWbS+FEEIIIcQNocLBpYuLC3FxccybN6/UvUuvWVtbS3AphLhhKIpCkbkIG51NTXdFCCFuCv9v796DoyoPN44/m82NkGTJjZBIWC6ViyTUIF5+YokFMo6MVEV0qNFWUacytToKKnWKrZcCWtShWm3ryEVrjcVWGChSsSCCiJRAkASLqAkkEBJIJDcuuez7+0OTGhVIwps9m833M8Msu9mz77PMMOfJe855T6fWuTwbrHMJBJfPj36uBl+DhscPdzpKhzX5mnTz6pv13y/+q9//8Pf6QX//XaEKdDfsv9Fe7T7nEgC+6fPqz3Xdyut0w8oblFee53ScDqtvrFdhZaGafE3aUbHD6TgAEBQolwA6rbG5Uc2+ZhkZnWg64XScDvNEeDT3B3N104ib9NORLKkG9ESlpaWqq6vz23gNDQ0qLi6Wz+fz25hnYjsTh8UBnJW88jydaDqhseeMdToKgC4UrPvvgQMH6le/+pVuv/12v4yXn5+vzMxMHT58uM29x51kO1OHrhYHgG9icXYANjX7mrW9YrsOHzuspKgkje47Wu4uuv1jWVmZmpqaVFlZqeLiYrlcLnm9XpWUlCg+Pl69e/dWeXm5wsLCFB8f37pdVVWVwsLCFBMT0+bzfD6f9u/fL0kKCwtTSkqKQkL+d5C4ublZBw9+ed/0kpIS1dXVqXfv3oqNjVVZWZm8Xq9cLpcqKioUHx+v0NAva1pDQ4OOHj2qvn37nvK7nCrTyZMn23z2kSNHlJCQIJfLddpMSUlJnf+HNX5WXV1tJJnq6mp/Dw0AADrJH/vvtcVrzYS/TTDpS9Jb/0z42wSztnhtl4w3efJk43a7TUJCgvF6vWbIkCHGGGOSk5PNjBkzzODBg02/fv3Mww8/bIwxZvPmzWbUqFEmOjraxMbGmszMTJOXl9f6eVVVVcbr9Rqv12tSUlJMr169zM9//nPT0NBgjDGmrKzMpKSkGEkmLS3NeL1eM2PGDPPBBx8YSWbWrFkmMTHReDwek5iYaFavXm3mzJlj4uLiTGxsrBk0aJDZtWtXm+9wpkwtn/3YY4+ZpKQkExcXZ+Li4szy5ctPm+lsUC4BBJ0d5TvM7Pdmm7xDee14N4D26Or999ritSZjSUabYpm+JN1kLMkwGUsyuqxger1e8+KLL7Z5LTk52SQlJZmPPvqo9bV9+/YZj8djFi1aZIwxxufzmd/+9remf//+pra29js/u6ioyJx77rlm4cKFra/t2LHDSDKHDx9ufa2lAF577bWmvr7e+Hw+c9NNN5mwsDAzbdo0c/z4cdPU1GSmTJlirrjiig5lavnsq6++uvW13/zmNyYhIcGcPHnylJnOBhf0AAg687fO16rPV2ne1nlORwHQDs2+Zs3fOl9G374MpOW1J7Y+oWZfs98y/exnP1NGRkbr8xdffFHDhg1Tdna2SkpKVFJSopycHFVVVWnLli1tMxujiooKSdLkyZO1Zs2ado05d+5cRUVFyeVyadq0aWpsbNS8efMUGRkpt9utG264QXl5/1uZoyOZ5s2bp+joaEnSbbfd1noqQFfgnEsAQefytMu1p2qPLk+73OkoANphe8V2lR8rP+XPjYwOHTuk7RXbdWG/C/2SaciQIW2eFxQUaPfu3brsssvavJ6UlKSqqipJUlNTk+69914tWrRIERERio2NVU1NjdLS0to1Zv/+/Vv/3nLL7W++Vltb26FMLb6eoaVk1tTUtCtXR1EuAQSdO79/p+78/p1OxwDQToePHbb6PhtaLqZp4Xa7lZWVpVWrVp1ymz/+8Y9asWKF8vPzde6550qSHn/8ceXm5nZJxvZkcgKHxQEAgKOSotp3ZXJ739cRERERampqOuP7Lr30Um3cuFGHD7ctuMaY1vUhCwoKNHbs2NZiKUnvvvvut8bTV7OcZ6s9mdrDZiZRLgEAgNNG9x2t5KhkueT6zp+75FK/qH4a3Xe09bGHDh2qt99+W5988on27dt3yvfdcccdSk1N1aRJk/Svf/1LBQUFys3N1SWXXKLKykpJ0oUXXqg1a9Zo1apV2rlzp2bOnKn169e3+Zy0tDRFRkZq2bJlKioq+lYx7Ij2ZGoPm5lEuQQAAE5zh7g1+6LZ0ldF8utanj940YNdst7lE088ocbGRk2ePFkTJkyQvipbLecltoiJidHmzZs1fvx4PfTQQ/rxj3+slStX6vnnn29dE3L69OmaOXOm5syZo5ycHNXV1Wnu3Lk655xzWj8nOjpaixcv1muvvaYJEybo17/+tSIiIuT1etusiRkZGdm6NmWLqKgoeb3eDmX6rs8OCQmR1+ttnbH8rkxngzv0AAFk88HN2lmxUznn5Sg23Jn/Hz7j04aSDRrSZ4gGxA447XsP1h1UclRyly1wDCBw+GP//c6+dzR/6/w2F/f0i+qnBy96UBO9E7tkTNjHBT1AgDjZfFJ3/fsuNfoaVd9Yr1kXznIkx6KCRVq4faE8ER6tu36dwt3h3/m+Bf9ZoKW7l+rKgVfqyawn/Z7zZPNJRbgjVHGsQg+894ASIhM0/wfzFeYO83sWAHZM9E7UD9N+6Lc79KBrUC6BABEWEqYhfYZoT9UeDYsf5lgOt8vd+niq858k6dOjn0qS9h7d26lx1u5bq+fzn9fUoVOVMyKnQ9s+nfe0Fhcs1vT06Urtnaq88i/Xfbtl5C3KSMo44/YAApc7xO235YbQNSiXQIAIcYXor5P+quqGaiX2SvTr2IfqD+mRDx7RYM9gzRozS8PihmmQZ9BpZwEf/r+HtfzT5cr2ZndqzNf++5o+PfqpXi58ucPl8sOyD1sfnx3/rNbuW6v4XvEaHj+8U1kAAPZwziUA/Wnnn/Rc/nOSpNXXrlZabPsW/D0bG0o26IWdL2jq0KmaOnRqh7YtPFKov+/9u6479zqNTBzZZRnbo6G5QZXHK5USneJoDqCrsf9GezFzCfjBiaYTigyNdDrGKU30TtSa4jUa5Bnkt5KUlZalrLSsTm07MnGk46Wyxa1rbtVHRz7S7Itmd3gGFgCCEeUS6GIz3pmh9w+8r0fHPqprvneN03G+05A+Q/Tm1W86HaPbMcbos+rPJEmfHf3M6TgAEBBY5xLoYtsObZOR0bZD25yOAstcLpdemPiC7jr/Lt0z+h6n4wBAQGDmEuhiC7IW6P2D72t6+nSno6ALZPbNVGbfTKdjAEDAoFwCXexszi0EAKC74bA4AAAArKFcAgAAwBrKJQDryurK9NS2p7S9fLvTUQAAfka5BGDd77b9TksKl2jmhplORwEA+BnlEoB1oxJHtXkEAPQc3P4RgBW7K3ersLJQPxryI0W4I1R9slox4TEKcfE7LBAM2H+jvViKCMBZa2hu0K1rbtWxpmMqqyvT3aPvlifC43QsAIADmFIAcNbcLrfiIuMkSQm9EpyOAwBwEDOXAM6aO8St1696XYfqD2lY/DCn4wAAHES5BGCFJ8LDoXAAAIfFAQAAYA/lEgAAANZQLgFYdazxmHzG53QMAIBDKJcArFm/f73GvjZWN/7zRjX7mp2OAwBwAOUSgDX5h/PVZJr0cdXHOtZ0zOk4AAAHcLU4AGump09Xs69Z6UnpigmPcToOAMABlEsA1ngiPJp14SynYwAAHMRhcQAAAFhDuQQAAIA1lEsAAABYQ7kEAACANZRLAAAAWEO5BAAAgDWUSwBW3Lv+Xl306kV6t+Rdp6MAABxEuQS62I6KHbpp9U1aWrjU6ShdptHXqHUl63S86bg2lG5wOs5pHW86rpd2vaRNBzY5HQUAghKLqANd7C+7/6Kdh3dqT9Ue/XTkT52O0yXCQsL0yKWPaEvZFt2WfpvTcU5raeFS/SH/Dwp1hWr9DevVJ7KP05EAIKhQLoEudt3Q6/Tp0U+V7c12OkqXuuZ71+ia713jdIwzGuwZLJdcSo1OVVRYlNNxACDouIwxxp8D1tTUyOPxqLq6WrGxsf4cGsA3fF79uR7f8rjSE9N13wX3OR3Hb8rryxUbEateob2cjgJ0G+y/0V6ccwn0YG/ufVP/OfQfLS5YrKoTVU7H8Zvk3skUSwDoIhwWB3qwSYMm6f2D72tU4ijFRcQ5HQcAEAQol0APNiJhhP7xo384HQMAEEQ4LA4AAABrKJcAAACwhnIJAAAAayiXAAAAsIZyCQAAAGsolwAAALCGcgl0kM/49MWJL5yOAQBAQKJcAh10z7p7NO71cVpauNTpKAAABBzKJdBBu47savMIAAD+hzv0AB309OVP692Sd3XjiBudjgIAQMChXAIdNDp5tEYnj3Y6BgAAAYnD4gAAALCGcgkAAABrKJcAAACwhnIJAAAAayiXQDdXcaxCKz9bqbqGOqejAADA1eJAd3fXv+/Sx1Uf64qBV2hB1gKn4wAAejhmLoFuLios6svH0CinowAAwMwl0N09N/45FVYWanRf1t4EADiPcgl0c9Hh0bo45WKnYwAAIHFYHAAAADZRLgEAAGAN5RIAAADWUC4BAABgDeUSAAAA1lAuAQAAYA3lEujh6hrqVN9Y73QMAECQoFwCPVhRdZEmvjFR2W9kq7S21Ok4AIAgQLkEerD9NftV31iv2oZaldZRLgEAZ4879AA92Lj+4zT7otlyu9y6JOUSp+MAAIIA5RLowVwul3JG5DgdAwAQRDgsDgAAAGsolwAAALCGcgkAAABrKJcAAACwhnIJAAAAayiXAAAAsIZyCQAAAGsol+gWqk5UyWd8TscAAABnQLlEwFtUsEhZr2fp3vX3Oh0FAACcAeUSAa/gSEGbRwAAELi4/SMC3v1j7lf/mP4anzbe6SgAAOAMKJcIeCnRKbrvgvucjgEAANqBw+IATssYozXFa5Rfke90FABAN8DMJYDTWvHZCs15f45CXaH655R/KjU61elIAIAAxswlgNOKDY+VJEWERijcHe50HABAgGPmEsBpjR8wXssmL1OfiD5K7JXodBwAQICjXAI4o+Hxw52OAADoJjgsDgAAAGsolwAAALCGcgkAAABrKJc93NLCpbp73d0qqSlxOgoAAAgCXNDTg9U31mvBtgWSpJTeKfrlxb90OhIAAOjmmLnswXqH9daVA69U36i+muid6HQcAAAQBJi57OGezHrS6QgAACCIMHMJAAAAa4Jm5rLR16hntz8rl8ulX2T+QqEhQfPVAAAAuo2gaWAbSzdqceFiSdIFyRdoXP9xTkcCAADocYLmsPh5CecppXeKzok+RyPiRzgdBwgKPuPT3i/2qtHX6HQUAEA3ETQzl/1699PbU992OgYQVB7b8pje+OQNZXuz9fTlTzsdBwDQDQTNzCUA+w7UHvjyse6A01EAAN1E0MxcAvi2t4vf1p8/+rOmDZ+mqUOndnj7x8Y+preK3tKEARO6JB8AIPgwcwkEsVd2v6I9X+zR4oLFndo+uXeybkm/RWmxadazAQCCEzOXQBD7ycif6MRHJzRt2DSnowAAegjKJRDEsr3ZyvZmOx0DANCDcFgcAAAA1lAuAQAAYA3lEgAAANZQLgEAAGAN5RIAAADWUC4BAABgDeUSAAAA1lAuAQAAYA3lEgAAANZQLgEAAGAN5RIAAADWUC4BAABgDeUSAAAA1lAuAQAAYA3lEgAAANZQLgEAAGAN5RIAAADWUC4BAABgDeUSAAAA1lAuAQAAYA3lEgAAANZQLgEAAGAN5RIAAADWUC4BAABgDeUSAAAA1lAuAQAAYA3lEgAAANZQLgEEjNLaUj2w4QG98ckbTkcBAHQS5RJAwFhSuERvFb+lRz94VA3NDU7HAQB0AuUSQMDI6p+lmLAYTRgwQeHucKfjAAA6wWWMMf4csKamRh6PR9XV1YqNjfXn0AAAoJPYf6O9mLkEAACANZRLAAAAWEO5BAAAgDWUSwAAAFhDuQROwWd8WlSwSK/sfsXpKAAAdBuhTgcAAtX6/ev1TN4zkqShcUN1ccrFTkcCACDgUS6BUxjcZ7BiwmIUGhIqb6zX6TgAAHQLlEvgFAZ5BmndDevkcrkU4Y5wOg4AAN0C5RI4jcjQSKcjAADQrXBBDwAAAKyhXAIAAMAayiUAAACsoVwCAADAGsolAAAArKFcAgAAwBrKJQAAAKyhXAIAAMAayiUAAACsoVwCAADAGsolAAAArKFcAgAAwBrKJQAAAKwJ9feAxhhJUk1Njb+HBgAAndSy327ZjwOn4vdyWVtbK0lKS0vz99AAAOAs1dbWyuPxOB0DAcxl/PwriM/n08GDBxUTEyOXy+XPoQEAQCcZY1RbW6vU1FSFhHBWHU7N7+USAAAAwYtfPQAAAGAN5RIAAADW+P2CHgDdy6ZNm1RaWipJCg8P14ABA5SZmSm3293mfY2Njdq+fbvKysqUnJys88477ztP+t+1a5cKCws1btw4paam+u17AAD8g3IJ4LQWLFigbdu26bLLLlNDQ4O2bNmiPn36aPXq1Ro4cKAk6dVXX9XMmTMVGxurkSNH6siRIyoqKtK0adO0YMECSdLGjRv10EMPqby8XHv37tXKlSsplwAQhCiXAM5ozJgxys3NlSTV1dUpMzNT999/v5YtW6bc3FzdfPPNeumll3Trrbe2bnPy5Em9/PLLrc9ra2s1d+5cZWRkKC4uzpHvAQDoepRLAB0SHR2tiRMn6p133pExRg888ICuv/76NsVSkiIiInTHHXe0Pp80aZIk6ejRo37PDADwHy7oAdBhe/fuVUpKinbv3q2SkhJNmTLF6UgAgADBzCWAMzpw4IByc3PV0NCgdevWad26dVqxYoUOHTokSRowYIDTEQEAAYJyCeCMysrKtHz5coWHhystLU35+fkaNWqUtm7dKkmqrKx0OiIAIEBQLgGc0dcv6Pm6jIwMRUZGasuWLbrqqqscyQYACCyccwmg03r16qUZM2Zo4cKFKioq+tbP9+/f70guAIBzmLkEcFbmzZun4uJiXXDBBbr99tuVnp6uI0eOKC8vT4WFhcrPz5cklZaWatOmTTp27Jj01bqXdXV1Gj58uM4//3yHvwUAwBaXMcY4HQJA4HrqqackSTNnzjzt+9577z2tXr269Q49Y8aM0ZQpUxQa+uXvsB9++KGeeeaZb203efJk5eTkdFF6AIC/US4BAABgDedcAgAAwBrKJQAAAKyhXAIAAMAayiUAAACsoVwCAADAGsolAAAArKFcAgAAwBrKJQAAAKyhXAIAAMAayiUAAACsoVwCAADAGsolAAAArPl/sK644rqt0CAAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "scaled = adata.copy()\n", "sc.pp.scale(scaled)\n", "scaled.X = np.nan_to_num(scaled.X)\n", "sc.pp.pca(scaled, n_comps=10)\n", "sc.pl.pca(scaled, color=\"Metadata_Control_Type\", title=\"JUMP CRISPR genes in profile space\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.13.15" } }, "nbformat": 4, "nbformat_minor": 5 }