"
],
"text/plain": [
"chain_status Extra pair Orphan VDJ Orphan VJ Single pair\n",
"has_contig \n",
"True 81 5 16 882"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pd.crosstab(adata.obs[\"has_contig\"], adata.obs[\"chain_status\"])"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now actually filter the AnnData object and run through a standard workflow starting by filtering genes and normalizing the data\n",
"\n",
"Because the 'filtered' `AnnData` object was returned as a filtered but otherwise unprocessed object, we still need to normalize and run through the usual process here. The following is just a standard scanpy workflow."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"adata = adata[\n",
" adata.obs[\"filter_rna\"] == \"False\"\n",
"] # from ddl.pp.recipe_scanpy_qc\n",
"# filter genes\n",
"sc.pp.filter_genes(adata, min_cells=3)\n",
"# Normalize the counts\n",
"sc.pp.normalize_total(adata, target_sum=1e4)\n",
"# Logarithmize the data\n",
"sc.pp.log1p(adata)\n",
"# Stash the normalised counts\n",
"adata.raw = adata"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Identify highly-variable genes"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABDMAAAGwCAYAAAC0Ky5qAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAA3U1JREFUeJzs3Xl4E9X6B/BvkpaWQBtIZSlSUqAqoiDFglorqCwVNxC1ULyuuFQvQgGFYFVArAYRqIpSRQW9ioIgXK9bWRSx1qWlbCIiWxcUBFtIgdIC6fn94W/GJE3STJq1/X6eJ49kMjN5s9hzcuY971EJIQSIiIiIiIiIiEKEOtABEBEREREREREpwcEMIiIiIiIiIgopHMwgIiIiIiIiopDCwQwiIiIiIiIiCikczCAiIiIiIiKikMLBDCIiIiIiIiIKKRzMICIiIiIiIqKQEhboAIJdXV0d/vjjD0RFRUGlUgU6HCIiogYJIXD8+HF06tQJajWvWzRX7MMQEVGoUdKH4WBGA/744w/ExcUFOgwiIiLFysvL0blz50CHQQHCPgwREYUqd/owHMxoQFRUFIC/38zo6OgAR0NERNSwqqoqxMXFyW0YNU/swxARUahR0ofhYEYDpLTM6OhodgSIiCikcGpB88Y+DBERhSp3+jCcSEtEREREREREIYWDGUREREREREQUUjiYQUREREREREQhhTUziChkWCwWnDlzJtBhEAWFFi1acNlVIiIiarY4mEFEQU8IgUOHDuHYsWOBDoUoaKjVanTt2hUtWrQIdChEREREfsfBDCIKetJARvv27aHVarlCAzV7dXV1+OOPP3Dw4EF06dKF/08QERFRs8PBDCIKahaLRR7IiImJCXQ4REGjXbt2+OOPP3D27FmEh4cHOhwiIiIiv+JkWyIKalKNDK1WG+BIiIKLNL3EYrEEOBIiIiIi/+NgBhGFBKbRE9ni/xNERETUnHEwg4iIiIiIiIhCCgcziIiIiIiIiCikcDCDiMhHrr76amRmZrrcR6VSYfXq1W6fc8OGDVCpVC6XqZ0xYwb69Onj9jn9acmSJWjTpo2iY+Lj45GTk+NyH6XvIxERERGFNq5mQkQUQAcPHkTbtm0DHYbfjBo1Ctdff32gwyAiIiKiEMfMjGYkNzcX8fHxyM3NDXQoRPT/OnbsiIiIiECH4RdnzpxBy5Yt0b59+0CHQkQhhn0YIiKyx8GMZsRkMqG0tBQmkynQoRA1G3V1dZgyZQr0ej06duyIGTNm2DxuPz2ioKAAffr0QWRkJJKSkrB69WqoVCps2bLF5rhNmzYhKSkJWq0WycnJ2LVrl8Pn37hxI8LDw3Ho0CGb7ZMnT8aAAQMcHpOeno7Ro0fbbDtz5gzOOeccLF68GADw5ZdfIiUlBW3atEFMTAxuvPFG7N27V96/pKQEKpUKy5cvx9VXX43IyEi899579aaZ7N27F8OHD0eHDh3QunVr9OvXD+vWrasX0/HjxzFmzBi0bt0anTp1wiuvvOIwdsnvv/+OUaNGoW3btoiJicHw4cNRUlLi8phPPvkE5513Hlq2bIlrrrkG77zzTr0pPQUFBRgwYABatmyJuLg4jB8/HidPnpQfj4+Px3PPPYf77rsPUVFR6NKlC9544w1FsW3YsAH9+/dHq1at0KZNG1x55ZUoLS11GTtRU8c+DBER2eNgRjNiNBphMBhgNBoDHQpRs/HOO++gVatW+PHHH/HCCy/gmWeewdq1ax3ue/z4cdx0003o1asXiouLMWvWLEydOtXhvllZWZg7dy6KiooQFhaG++67z+F+AwYMQLdu3fCf//xH3nb27Fm89957uPfeex0ec8cdd+CTTz7BiRMn5G15eXk4efIkbr31VgDAyZMnMWnSJBQWFmL9+vVQq9W45ZZbUFdXZ3OuqVOnYvz48di5cydSU1PrPdeJEydw/fXXY926ddi8eTNSU1Nx0003oayszGa/OXPmoHfv3iguLsa0adMwceJEp+9jdXU1rrnmGrRu3RobN25Efn4+Wrdujeuuuw6nT592eExJSQluu+02jBgxAlu2bMFDDz2ErKwsm322b9+O1NRUjBw5Etu2bcOyZcuQn5+PcePG2ew3d+5cJCUlYfPmzXjkkUfw8MMP49dff3UrtrNnz2LEiBEYOHAgtm3bhu+//x4PPvggl2GlZo99GCIiqkeQS2azWQAQZrM50KEQNUunTp0Sv/zyizh16pRXzrdw4UJhMBjEwoULvXI+VwYOHChSUlJstvXr109MnTpVvg9ArFq1So4tJibG5rUuWrRIABCbN28WQgjx9ddfCwBi3bp18j6fffaZACAfN336dHHJJZfIj8+ePVtceOGF8v3Vq1eL1q1bixMnTjiM+/Tp0+Kcc84R7777rrwtPT1d3H777U5f6+HDhwUAsX37diGEEPv37xcARE5Ojs1+ixcvFjqdzul5hBCiZ8+e4pVXXpHvGwwGcd1119nsM2rUKDFs2DD5vvX7+NZbb4kLLrhA1NXVyY/X1taKli1biry8PIfPOXXqVHHxxRfbbMvKyhIAxNGjR4UQQtx5553iwQcftNnn22+/FWq1Wn7vDQaD+Ne//iU/XldXJ9q3by9/3xqKraKiQgAQGzZscPkeCeH6/w22XSQEvwdERBR6lLRdzMwgombF36nKvXv3trkfGxuLw4cPO9x3165d6N27NyIjI+Vt/fv3b/C8sbGxAOD0vPfccw/27NmDH374AQDw9ttvIy0tDa1atXK4f3h4OG6//Xa8//77AP7Owvjvf/+LO+64Q95n7969GDNmDLp164bo6Gh07doVAOplVCQlJTl8DsnJkycxZcoU9OzZE23atEHr1q3x66+/1jvPFVdcUe/+zp07HZ5z06ZN2LNnD6KiotC6dWu0bt0aer0eNTU1NlNhrO3atQv9+vWz2Wb/3m/atAlLliyRz9m6dWukpqairq4O+/fvl/ez/mxUKhU6duwofzYNxabX63HPPffIGSovvfQSDh486PI9JCIiImqOuJoJETUrRqMRJpPJb6nK4eHhNvdVKlW9qRgSIUS96QRCiAbPKx3j7Lzt27fHTTfdhMWLF6Nbt274/PPPsWHDBpdx33HHHRg4cCAOHz6MtWvXIjIyEsOGDZMfv+mmmxAXF4dFixahU6dOqKurw8UXX1xvGoezARPJ448/jry8PLz44otISEhAy5YtcdtttzmdDmLN2dSLuro6XHrppfJgjLV27do5PMad976urg4PPfQQxo8fX+/4Ll26yP929Zm7E9vixYsxfvx4fPnll1i2bBmefPJJrF27FpdffrnD2ImIiIiaIw5mEFGzkpGRgYyMjECH4VCPHj3w/vvvo7a2Vl7hpKioyCvnvv/++zF69Gh07twZ3bt3x5VXXuly/+TkZMTFxWHZsmX44osvcPvtt6NFixYAgIqKCuzcuROvv/46rrrqKgBAfn6+R3F9++23uOeee3DLLbcA+LuGhqNCnVJWifX9Hj16ODxn3759sWzZMrRv3x7R0dFuxdGjRw98/vnnNtvs3/u+fftix44dSEhIcOucjYktMTERiYmJmDZtGq644gosXbqUgxlEREREVjjNhIgoSIwZMwZ1dXV48MEHsXPnTjljAXCeheCu1NRU6HQ6PPvss04Lf1pTqVQYM2YMcnNzsXbtWvzrX/+SH5NW4XjjjTewZ88efPXVV5g0aZJHcSUkJODjjz/Gli1bsHXrVvk9sPfdd9/hhRdewG+//YZXX30VH330ESZMmODwnHfccQfOOeccDB8+HN9++y3279+Pb775BhMmTMCBAwccHvPQQw/h119/xdSpU/Hbb79h+fLlWLJkifxeAH8XM/3+++/x73//G1u2bMHu3bvxySef4NFHH3X79TYU2/79+zFt2jR8//33KC0txZo1a/Dbb7/hwgsvdPs5iIiIiJoDDmYQEQWJ6Oho/O9//8OWLVvQp08fZGVl4emnnwYAmzoanlCr1bjnnntgsVhw1113uXXMHXfcgV9++QXnnnuuTSaHWq3Ghx9+iE2bNuHiiy/GxIkTMWfOHI/imj9/Ptq2bYvk5GTcdNNNSE1NRd++fevtN3nyZGzatAmJiYmYNWsW5s6d63B1FADQarXYuHEjunTpgpEjR+LCCy/Efffdh1OnTjnNhujatStWrFiBjz/+GL1798bChQvl1UykLJnevXvjm2++we7du3HVVVchMTERTz31lFyzxB0NxabVavHrr7/i1ltvxfnnn48HH3wQ48aNw0MPPeT2cxARERE1ByrhbEI2AQCqqqqg0+lgNpvdTlcmIu+pqanB/v370bVr10b/oA9F77//Pu69916YzWa0bNmyUed64IEH8Oeff+KTTz7xUnRNW3Z2NnJzc1FeXh7oUBxy9f8G2y4C+D0gIqLQo6TtCprMjI0bN+Kmm25Cp06doFKpsHr1apvH77nnHqhUKpubO/OHV65ciZ49eyIiIgI9e/bEqlWrfPQKiIga791330V+fj7279+P1atXY+rUqUhLS2vUQIbZbMa6devw/vvvK5oS0dy89tprKCwsxL59+/Cf//wHc+bMwd133x3osIiIiIjIgaApAHry5ElccskluPfee3Hrrbc63Oe6667D4sWL5ftSMTpnvv/+e4waNQqzZs3CLbfcglWrViEtLQ35+fm47LLLvBo/EZE3HDp0CE8//TQOHTqE2NhY3H777cjOzm7UOYcPH46ffvoJDz30EIYMGeKlSJue3bt349lnn0VlZSW6dOmCyZMnY9q0aYEOi4iIiIgcCMppJiqVCqtWrcKIESPkbffccw+OHTtWL2PDlVGjRqGqqgpffPGFvO26665D27Zt8cEHHzg8pra2FrW1tfL9qqoqxMXFMUWTKECa+zQTImc4zYQawu8BERGFmpCcZuKODRs2oH379jj//PPxwAMP4PDhwy73//777zF06FCbbampqSgoKHB6zPPPPw+dTiff4uLivBI7EREREREREXlHyAxmDBs2DO+//z6++uorzJ07F4WFhbj22mttsijsHTp0CB06dLDZ1qFDBxw6dMjpMdOmTYPZbJZvwVr4jYiIiIiIiKi5CpqaGQ0ZNWqU/O+LL74YSUlJMBgM+OyzzzBy5Einx6lUKpv7Qoh626xFRETIy/ARERERERERUfAJmcwMe7GxsTAYDNi9e7fTfTp27FgvC+Pw4cP1sjWIiIiIiIiIKHSE7GBGRUUFysvLERsb63SfK664AmvXrrXZtmbNGiQnJ/s6PCIiIiIiIiLykaCZZnLixAns2bNHvr9//35s2bIFer0eer0eM2bMwK233orY2FiUlJTgiSeewDnnnINbbrlFPuauu+7Cueeei+effx4AMGHCBAwYMACzZ8/G8OHD8d///hfr1q1Dfn6+318fEZG7rr76avTp0wc5OTmBDoWIiIiIKCgFTWZGUVEREhMTkZiYCACYNGkSEhMT8fTTT0Oj0WD79u0YPnw4zj//fNx99904//zz8f333yMqKko+R1lZGQ4ePCjfT05OxocffojFixejd+/eWLJkCZYtW4bLLrvM76+PiMjehg0boFKpcOzYsUCHQkREREQUUoImM+Pqq6+GEMLp43l5eQ2eY8OGDfW23XbbbbjtttsaExoRUcg7c+YMwsPDAx0GEREREZFXBE1mBhFRU1NbW4vx48ejffv2iIyMREpKCgoLCwEAJSUluOaaawAAbdu2hUqlwj333CMfW1dXhylTpkCv16Njx46YMWOGzbnNZjMefPBBtG/fHtHR0bj22muxdetW+fEZM2agT58+ePvtt9GtWzdEREQ4HTBetGgR4uLioNVqccstt2DevHlo06aNzT7/+9//cOmllyIyMhLdunXDzJkzcfbsWflxlUqFN998E7fccgu0Wi3OO+88fPLJJzbn+OWXX3D99dejdevW6NChA+6880789ddf8uMrVqxAr1690LJlS8TExGDw4ME4efKk2+83ERERETUfHMwgIvKRKVOmYOXKlXjnnXdQXFyMhIQEpKamorKyEnFxcVi5ciUAYNeuXTh48CBeeukl+dh33nkHrVq1wo8//ogXXngBzzzzjFzQWAiBG264AYcOHcLnn3+OTZs2oW/fvhg0aBAqKyvlc+zZswfLly/HypUrsWXLFocxfvfdd8jIyMCECROwZcsWDBkyBNnZ2Tb75OXl4V//+hfGjx+PX375Ba+//jqWLFlSb7+ZM2ciLS0N27Ztw/XXX4877rhDjufgwYMYOHAg+vTpg6KiInz55Zf4888/kZaWJj+enp6O++67Dzt37sSGDRswcuRIlxl7RERERNSMCXLJbDYLAMJsNgc6FKJm6dSpU+KXX34Rp06dCnQoipw4cUKEh4eL999/X952+vRp0alTJ/HCCy8IIYT4+uuvBQBx9OhRm2MHDhwoUlJSbLb169dPTJ06VQghxPr160V0dLSoqamx2ad79+7i9ddfF0IIMX36dBEeHi4OHz7sMs5Ro0aJG264wWbbHXfcIXQ6nXz/qquuEs8995zNPv/5z39EbGysfB+AePLJJ21ev0qlEl988YUQQoinnnpKDB061OYc5eXlAoDYtWuX2LRpkwAgSkpKXMZL/3D1/wbbLhKC3wMiIgo9StquoKmZQUTkD0VFRcjPz0dKSgqSkpJ89jx79+7FmTNncOWVV8rbwsPD0b9/f+zcubPB43v37m1zPzY2FocPHwYAbNq0CSdOnEBMTIzNPqdOncLevXvl+waDAe3atXP5PLt27bJZFQoA+vfvj08//VS+v2nTJhQWFtpkYlgsFtTU1KC6uhparbZezK1atUJUVJRNzF9//TVat25dL4a9e/di6NChGDRoEHr16oXU1FQMHToUt912G9q2besyfiIiIiJqnjiYQUTNSn5+PsxmM/Lz8306mCH+f3qESqWqt91+myP2xTpVKhXq6uoA/F1PIzY21mHRY+taF61atXIrTkcxWqurq8PMmTMxcuTIesdHRka6HfNNN92E2bNn1ztHbGwsNBoN1q5di4KCAqxZswavvPIKsrKy8OOPP6Jr164Nvg6ipiY+Ph6lpaX1tj/yyCN49dVXAxARERFRcOFgBhE1KykpKXJmhi8lJCSgRYsWyM/Px5gxYwD8vaJIUVERMjMzAQAtWrQA8HeWgxJ9+/bFoUOHEBYWhvj4+EbF2aNHD/z0008224qKiuo9365du5CQkODx8/Tt2xcrV65EfHw8wsIcNz0qlQpXXnklrrzySjz99NMwGAxYtWoVJk2a5PHzEoWqwsJCm78NP//8M4YMGYLbb789gFEREREFDw5mEFGzkpSU5NOMDEmrVq3w8MMP4/HHH4der0eXLl3wwgsvoLq6GmPHjgXw9zQQlUqFTz/9FNdffz1atmzpcBqGvcGDB+OKK67AiBEjMHv2bFxwwQX4448/8Pnnn2PEiBGKXt+jjz6KAQMGYN68ebjpppvw1Vdf4YsvvrDJ1nj66adx4403Ii4uDrfffjvUajW2bduG7du349lnn3Xref79739j0aJFSE9Px+OPP45zzjkHe/bswYcffohFixahqKgI69evx9ChQ9G+fXv8+OOPOHLkCC688EK3XwtRU2I/RcxkMqF79+4YOHBggCIiIiIKLlzNhIjIR0wmE2699Vbceeed6Nu3L/bs2YO8vDy5DsS5556LmTNnwmg0okOHDhg3bpxb51WpVPj8888xYMAA3HfffTj//PMxevRolJSUoEOHDopivPLKK5Gbm4t58+bhkksuwZdffomJEyfaTB9JTU3Fp59+irVr16Jfv364/PLLMW/ePBgMBrefp1OnTvjuu+9gsViQmpqKiy++GBMmTIBOp4NarUZ0dDQ2btyI66+/Hueffz6efPJJzJ07F8OGDVP0eoiaotOnT+O9997Dfffd53KaWm1tLaqqqmxuRERETZVK2E+OJhtVVVXQ6XQwm82Ijo4OdDhEzU5NTQ3279+Prl272vzAJt954IEH8Ouvv+Lbb78NdCjkgqv/N9h2NS3Lly/HmDFjUFZWhk6dOjndb8aMGZg5c2a97fweEBFRqFDSh2FmBhFRM/fiiy9i69at2LNnD1555RW88847uPvuuwMdFhH9v7feegvDhg1zOZABANOmTYPZbJZv5eXlfoqQiIjI/1gzg4iomfvpp5/wwgsv4Pjx4+jWrRtefvll3H///YEOi4gAlJaWYt26dfj4448b3DciIgIRERF+iIqIiCjwOJhBRNTMLV++PNAhEJETixcvRvv27XHDDTcEOhQiIqKgwmkmREREREGorq4Oixcvxt133+10SWMiIqLmioMZISQ3Nxfx8fHIzc0NdChEfsdaxUS2+P9E07du3TqUlZXhvvvuC3QoREREQYeDGSHEZDKhtLQUJpMp0KEQ+U14eDgAoLq6OsCREAWX06dPAwA0Gk2AIyFfGTp0KIQQOP/88wMdChERUdBhzmIIMRqNMJlMMBqNgQ6FyG80Gg3atGmDw4cPAwC0Wi1UKlWAoyIKrLq6Ohw5cgRarZbTD4iIiKhZYg8ohGRkZCAjIyPQYRD5XceOHQFAHtAgIkCtVqNLly4c3CMiIqJmiYMZRBT0VCoVYmNj0b59e5w5cybQ4RAFhRYtWkCt5mxRIiIiap44mEFEIUOj0bA+ABERERERsQAoEREREREREYUWZmYQERERecGuXbvwwQcf4Ntvv0VJSQmqq6vRrl07JCYmIjU1FbfeeisiIiICHSYREVGTwMwMIiIiokbYvHkzhgwZgksuuQQbN25Ev379kJmZiVmzZuFf//oXhBDIyspCp06dMHv2bNTW1gY6ZCIiopDHzAwiIiKiRhgxYgQef/xxLFu2DHq93ul+33//PebPn4+5c+fiiSee8GOERERETQ8HM4iIiIgaYffu3WjRokWD+11xxRW44oorcPr0aT9ERURE1LRxmgkRERFRI7gzkNGY/YmIiKg+ZmYQERERNcLLL7/s9r7jx4/3YSRERETNBwcziIiIiBph/vz5NvePHDmC6upqtGnTBgBw7NgxaLVatG/fnoMZREREXsJpJkRERESNsH//fvmWnZ2NPn36YOfOnaisrERlZSV27tyJvn37YtasWYEOlYiIqMngYAYRERGRlzz11FN45ZVXcMEFF8jbLrjgAsyfPx9PPvlkACMjIiJqWjiYEURyc3MRHx+P3NzcQIdCREREHjh48CDOnDlTb7vFYsGff/4ZgIiIiIiaJg5mBBGTyYTS0lKYTKZAh0JEREQeGDRoEB544AEUFRVBCAEAKCoqwkMPPYTBgwcHODoiIqKmg4MZQcRoNMJgMMBoNAY6FCIiIvLA22+/jXPPPRf9+/dHZGQkIiIicNlllyE2NhZvvvlmoMMjIiJqMlRCumxADlVVVUGn08FsNiM6OjrQ4RARETWIbVfg/fbbb/j1118hhMCFF16I888/3+8x8HtAREShRknbpWhpVrPZjFWrVuHbb79FSUkJqqur0a5dOyQmJiI1NRXJycmNCpyIiIioKYiPj4cQAt27d0dYmKLuVpOTm5sLk8kEo9GIjIyMQIdDRERNhFvTTA4ePIgHHngAsbGxeOaZZ3Dy5En06dMHgwYNQufOnfH1119jyJAh6NmzJ5YtW+brmImIiIiCUnV1NcaOHQutVouLLroIZWVlAIDx48c325pYrAlGRES+4NalgksuuQR33XUXfvrpJ1x88cUO9zl16hRWr16NefPmoby8HI899phXAyUiIiIKdtOmTcPWrVuxYcMGXHfddfL2wYMHY/r06c2yLpbRaJQzM4iIiLzFrZoZR44cQbt27dw+qdL9gxnnmxIRUahh2xU4BoMBy5Ytw+WXX46oqChs3boV3bp1w549e9C3b19UVVX5LRZ+D4iIKNQoabvcmmaidGDCk4GMjRs34qabbkKnTp2gUqmwevVq+bEzZ85g6tSp6NWrF1q1aoVOnTrhrrvuwh9//OHynEuWLIFKpap3q6mpURwfERERUUOOHDmC9u3b19t+8uRJqFSqAERERETUNLk1zeSTTz5x+4Q333yzR4GcPHkSl1xyCe69917ceuutNo9VV1ejuLgYTz31FC655BIcPXoUmZmZuPnmm1FUVOTyvNHR0di1a5fNtsjISI9iJCIiInKlX79++Oyzz/Doo48CgDyAsWjRIlxxxRWBDI2IiKhJcWswY8SIETb3VSoVrGenWF9psFgsHgUybNgwDBs2zOFjOp0Oa9eutdn2yiuvoH///igrK0OXLl2cnlelUqFjx44exURERESkxPPPP4/rrrsOv/zyC86ePYuXXnoJO3bswPfff49vvvkm0OERERE1GW5NM6mrq5Nva9asQZ8+ffDFF1/g2LFjMJvN+Pzzz9G3b198+eWXvo5XZjaboVKp0KZNG5f7nThxAgaDAZ07d8aNN96IzZs3u9y/trYWVVVVNjciIiIidyQnJ+O7775DdXU1unfvjjVr1qBDhw74/vvvcemllwY6PCIioiZD8cLnmZmZyM3NRUpKirwtNTUVWq0WDz74IHbu3OnVAB2pqamB0WjEmDFjXBYF6dGjB5YsWYJevXqhqqoKL730Eq688kps3boV5513nsNjnn/+ecycOdNXoRMREVET16tXL7zzzjuBDoOIiKhJcyszw9revXuh0+nqbdfpdCgpKfFGTC6dOXMGo0ePRl1dHV577TWX+15++eX417/+hUsuuQRXXXUVli9fjvPPPx+vvPKK02OmTZsGs9ks38rLy739EoiIiKiJuvbaax1eFDl69CiuvfbaAERERETUNCkezOjXrx8yMzNx8OBBeduhQ4cwefJk9O/f36vB2Ttz5gzS0tKwf/9+rF27VvEyY2q1Gv369cPu3bud7hMREYHo6GibGxEREZE7NmzYgAULFmDEiBE4efKkvP306dOsmUFERORFigcz3n77bRw+fBgGgwEJCQlISEhAly5dcPDgQbz11lu+iBHAPwMZu3fvxrp16xATE6P4HEIIbNmyBbGxsT6IkIiIiAhYt24dDh06hMsvv9wvWatERETNkeKaGQkJCdi2bRvWrl2LX3/9FUII9OzZE4MHD27U+uknTpzAnj175Pv79+/Hli1boNfr0alTJ9x2220oLi7Gp59+CovFgkOHDgEA9Ho9WrRoAQC46667cO655+L5558HAMycOROXX345zjvvPFRVVeHll1/Gli1b8Oqrr3ocJxEREZErsbGx+Oabb3DfffehX79++Oijj3DhhRcGOiwiIqImRfFgBvD3cqdDhw7FgAEDEBER0ahBDElRURGuueYa+f6kSZMAAHfffTdmzJiBTz75BADQp08fm+O+/vprXH311QCAsrIyqNX/JJscO3YMDz74IA4dOgSdTofExERs3LjR59NhiIiIqHmS+kQRERF4//338eyzz+K6667D1KlTAxwZERFR06ISQgglB9TV1SE7Oxu5ubn4888/8dtvv6Fbt2546qmnEB8fj7Fjx/oq1oCoqqqCTqeD2Wxm/QwiIgoJbLsCR61W49ChQ2jfvr28beXKlbj77rtx6tQpWCwWv8XC7wEREYUaJW2X4poZzz77LJYsWYIXXnhBnt4B/L0M2Ztvvqk8WiIiIqImYv/+/TjnnHNstt1666344Ycf8PbbbwcoKiIioqZH8WDGu+++izfeeAN33HEHNBqNvL1379749ddfvRocERERUSgxGAw2U14lF198Me6+++4AROR7ubm5iI+PR25ubqBDISKiZkRxzYzff/8dCQkJ9bbX1dXhzJkzXgmKiIiIKFSMHDkSS5YsQXR0NEaOHOly348//tjt8/7++++YOnUqvvjiC5w6dQrnn38+3nrrLVx66aWNDdmrTCYTSktLYTKZkJGREehwiIiomVA8mHHRRRfh22+/hcFgsNn+0UcfITEx0WuBEREREYUCnU4nF/7U6XReOefRo0dx5ZVX4pprrsEXX3yB9u3bY+/evWjTpo1Xzu9NRqMRJpMJRqMx0KEQEVEzongwY/r06bjzzjvx+++/o66uDh9//DF27dqFd999F59++qkvYiQiIiIKWosXL3b478aYPXs24uLibM4XHx/vlXN7W0ZGBjMyiIjI7xTXzLjpppuwbNkyfP7551CpVHj66aexc+dO/O9//8OQIUN8ESMRERFRs/LJJ58gKSkJt99+O9q3b4/ExEQsWrTI5TG1tbWoqqqyuRERkXcUFRUhJycHRUVFgQ6F/p/izAwASE1NRWpqqrdjISIiIgo5iYmJ8jSThhQXF7u13759+7Bw4UJMmjQJTzzxBH766SeMHz8eERERuOuuuxwe8/zzz2PmzJlux01ERO7Lz8+H2WxGfn4+kpKSAh0OwYPBjPvuuw8DBw6sV5G7qqoKmZmZXHaMiIiImpURI0Z4/Zx1dXVISkrCc889B+DvAZMdO3Zg4cKFTgczpk2bhkmTJsn3q6qqEBcX5/XYiIiao5SUFOTn5yMlJSXQodD/UwkhhJID1Go1WrZsibFjxyInJ0defuzPP/9Ep06dYLFYfBJooFRVVUGn08FsNiM6OjrQ4RARETWIbVfoMxgMGDJkCN58801528KFC/Hss8/i999/d+sc/B74X25urlwMlXVEiIiUU9J2Ka6ZAQCfffYZvvjiC6SmpuLo0aMeBUlEREREjl155ZXYtWuXzbbffvut3mpyFFysl6klIiLf8mgwo2fPnvjhhx9w5swZ9OvXDzt37vR2XEREREQhx2Kx4MUXX0T//v3RsWNH6PV6m5u7Jk6ciB9++AHPPfcc9uzZg6VLl+KNN97Av//9bx9GT41lNBphMBi4TC0RkR8oHsyQClzFxMRg3bp1uPrqq3H55Zfjk08+8XpwRERERKFk5syZmDdvHtLS0mA2mzFp0iSMHDkSarUaM2bMcPs8/fr1w6pVq/DBBx/g4osvxqxZs5CTk4M77rjDd8FTo2VkZKCkpIRTTIiI/MCjmhmHDh1C+/bt5W3z5s3D1KlTUVdXx5oZREREAca2K3C6d++Ol19+GTfccAOioqKwZcsWedsPP/yApUuX+i0Wfg+IiCjUKGm7FK9m8vXXX9dLk5w0aRJ69+6N7777TunpiIiIiJqMQ4cOoVevXgCA1q1bw2w2AwBuvPFGPPXUU4EMjYiIqElRPM1k4MCBCAurPwYyePBgTJ8+3StBEREREYWizp074+DBgwCAhIQErFmzBgBQWFiIiIiIQIZGRETUpLiVmTFp0iTMmjULrVq1slm/3JF58+Z5JTAiIiKiUHPLLbdg/fr1uOyyyzBhwgSkp6fjrbfeQllZGSZOnBjo8IiIiJoMtwYzNm/ejDNnzsj/dkYqDkpERETUHFkvyXnbbbehc+fOKCgoQEJCAm6++eYARkZERNS0KC4A2tyweBYREYUatl0E8HtAREShx6cFQImIiIjIud9//x3fffcdDh8+jLq6OpvHxo8fH6CoiIiImha3BjNGjhzp9gk//vhjj4MhIiIiCmWLFy9GRkYGWrRogZiYGJspuCqVioMZREREXuLWYIZOp/N1HEREREQh7+mnn8bTTz+NadOmQa1WvGgcERERucmtwYzFixf7Og4iIiKikFddXY3Ro0dzIIOIiMjH2NJSPbm5uYiPj0dubm6gQyEiIgopY8eOxUcffRToMIiIiJo8j1YzWbFiBZYvX46ysjKcPn3a5rHi4mKvBRcMmmMl8Pj4eJSWlsJgMKCkpCTQ4RARkULNse0KFhaLBTfeeCNOnTqFXr16ITw83ObxefPm+S0Wfg+IiCjUKGm7FGdmvPzyy7j33nvRvn17bN68Gf3790dMTAz27duHYcOGeRw0BQ+j0QiDwQCj0RjoUIiIiELKc889h7y8PPz555/Yvn07Nm/eLN+2bNkS6PCIiIiaDMWZGT169MD06dORnp6OqKgobN26Fd26dcPTTz+NyspKLFiwwFexBgSvahARUahh2xU4bdu2xfz583HPPfcEOhR+D4iIKOT4NDOjrKwMycnJAICWLVvi+PHjAIA777wTH3zwgQfhEhERETUNERERuPLKKwMdBoUQ1iojIvKM4sGMjh07oqKiAgBgMBjwww8/AAD2798PD8pvEBERETUZEyZMwCuvvBLoMCiEmEwmlJaWwmQyBToUIqKQ4tbSrNauvfZa/O9//0Pfvn0xduxYTJw4EStWrEBRURFGjhzpixiJiIiIQsJPP/2Er776Cp9++ikuuuiiegVAP/744wBFRsHKaDTCZDKxVhkRkUKKa2bU1dWhrq4OYWF/j4MsX74c+fn5SEhIQEZGBlq0aOGTQAOF802JiCjUsO0KnHvvvdfl44sXL/ZTJPweEBFR6FHSdinOzFCr1VCr/5mdkpaWhrS0NOVREhERETUhZ8+exdVXX43U1FR07Ngx0OEQERE1aYoHMwCgpqYG27Ztw+HDh1FXV2fz2M033+yVwIiIiIhCSVhYGB5++GHs3Lkz0KEQERE1eYoHM7788kvcdddd+Ouvv+o9plKpYLFYvBIYERERUai57LLLsHnzZhgMhkCHQkRE1KQpHswYN24cbr/9djz99NPo0KGDL2IiIiIiCkmPPPIIJk+ejAMHDuDSSy9Fq1atbB7v3bt3gCIjIiJqWhQXAI2OjsbmzZvRvXt3X8UUVFg8i4iIQg3brsCxrismUalUEEL4PYOV3wMiIgo1Pi0Aetttt2HDhg3NZjCDiIiIyF379+8PdAhERETNguLBjAULFuD222/Ht99+i169etVbP338+PFeC46IiIgolLBWBhERkX8oHsxYunQp8vLy0LJlS2zYsAEqlUp+TKVScTCjkXJzc2EymWA0GpGRkRHocIiIiEihvXv3IicnBzt37oRKpcKFF16ICRMmMKuViIjIi+pP7GzAk08+iWeeeQZmsxklJSXYv3+/fNu3b5/HgWzcuBE33XQTOnXqBJVKhdWrV9s8LoTAjBkz0KlTJ7Rs2RJXX301duzY0eB5V65ciZ49eyIiIgI9e/bEqlWrPI7RH0wmE0pLS2EymQIdChERESmUl5eHnj174qeffkLv3r1x8cUX48cff8RFF12EtWvXBjo8IiKiJkPxYMbp06cxatQohwWuGuPkyZO45JJLsGDBAoePv/DCC5g3bx4WLFiAwsJCdOzYEUOGDMHx48ednvP777/HqFGjcOedd2Lr1q248847kZaWhh9//NGrsXsiNzcX8fHxyM3NtdluNBphMBhgNBoDFBkRERF5ymg0YuLEifjxxx8xb948zJ8/Hz/++CMyMzMxderUQIdHRETUZChezWTixIlo164dnnjiCV/FBJVKhVWrVmHEiBEA/s7K6NSpk01HoLa2Fh06dMDs2bPx0EMPOTzPqFGjUFVVhS+++ELedt1116Ft27b44IMP3IrFV5XA4+PjUVpaCoPBgJKSEq+dl4iIiKtYBE5kZCS2b9+O8847z2b7b7/9ht69e6OmpsZvsfB7QEREocanq5lYLBa88MILyMvLQ+/evesVAJ03b57SUzZo//79OHToEIYOHSpvi4iIwMCBA1FQUOB0MOP777/HxIkTbbalpqYiJyfH6XPV1taitrZWvl9VVdW44J0wGo1ybQwiIiJqGtq1a4ctW7bUG8zYsmUL2rdvH6CoiIiImh7Fgxnbt29HYmIiAODnn3+2ecy6GKg3HTp0CADQoUMHm+0dOnRAaWmpy+McHSOdz5Hnn38eM2fObES07snIyGCBTyIioibmgQcewIMPPoh9+/YhOTkZKpUK+fn5mD17NiZPnhzo8IiIiJoMRYMZFosFM2bMQK9evaDX630Vk1P2gyVCiAYHUJQeM23aNEyaNEm+X1VVhbi4OA+idY4rlhARETVNTz31FKKiojB37lxMmzYNANCpUyfMmDGDK74RERF5kaIqnhqNBqmpqTCbzb6Kx6GOHTsCQL2MisOHD9fLvLA/TukxERERiI6Otrl5m7MVS5wVBSUiIqLQoFKpMHHiRBw4cABmsxlmsxkHDhzAhAkTfJbBSkRE1BwpXpKkV69ejVqC1RNdu3ZFx44dbZY0O336NL755hskJyc7Pe6KK66otwzamjVrXB7jD85WLOGyrERERE1HVFQUoqKiAh0GERFRk6R4MCM7OxuPPfYYPv30Uxw8eBBVVVU2N0+dOHECW7ZswZYtWwD8XfRzy5YtKCsrg0qlQmZmJp577jmsWrUKP//8M+655x5otVqMGTNGPsddd90lp3QCwIQJE7BmzRrMnj0bv/76K2bPno1169YhMzPT4zi9ISMjAyUlJfWmmLi7LCszOIiIiILTn3/+iTvvvBOdOnVCWFgYNBqNzY2IiIi8Q/HSrGr1P+Mf1umSUi0Ki8XiUSAbNmzANddcU2/73XffjSVLlkAIgZkzZ+L111/H0aNHcdlll+HVV1/FxRdfLO979dVXIz4+HkuWLJG3rVixAk8++ST27duH7t27Izs7GyNHjnQ7Ln8sa6a0hgaXdSUiIle4JGfgDBs2DGVlZRg3bhxiY2PrTS0ZPny4W+eZMWNGvYLkDRUxt8fvARERhRolbZfiwYxvvvnG5eMDBw5Ucrqg54+OgNLBCRYQJSIiV/gjNnCioqLw7bffok+fPo06z4wZM7BixQqsW7dO3qbRaNCuXTu3z8HvARERhRolbZfipVmb2mBFMDAajfLgRENyc3ORlZXlh6iIiIhIqbi4OCi8TuRUWFiYXASdiIiIbCmumQEAx44dw9y5c3H//ffjgQcewPz58/2+wklT4qyGhiMmkwmVlZWorKxUXCiUtTaIiIh8KycnB0aj0SvTQHfv3o1OnTqha9euGD16dIMF2Gtra71WyyzYsA9DRET2FA9mFBUVoXv37pg/fz4qKyvx119/Yd68eejevTuKi4t9EWOz5ajhNhqN0Ov10Ov1bmVyWONqKURERL41atQobNiwAd27d0dUVJTcZks3d1122WV49913kZeXh0WLFuHQoUNITk5GRUWF02Oef/556HQ6+RYXF+eNlxQU2Ich8r2ioiLk5OSgqKgo0KEQuUVxzYyrrroKCQkJWLRoEcLC/p6lcvbsWdx///3Yt28fNm7c6JNAAyWQ8029XeiTtTaIiJoH1koInHfeecfl43fffbdH5z158iS6d++OKVOmYNKkSQ73qa2tRW1trXy/qqoKcXFxfvse+LKfwT4Mke/l5OTAbDZDp9O5XP2xqKgI+fn5SElJQVJSkv8CpGbBpwVAW7Zsic2bN6NHjx4223/55RckJSWhurpaecRBzFcdQncaZTbcRETkCQ5mNE1DhgxBQkICFi5c6Nb+/v4ecLU1otDm7iCFu4MeRJ5Q0nYpnmYSHR2NsrKyetvLy8sRFRWl9HTNljvpkkpqaRAREVHTVVtbi507dyI2NjagcbiqXWE0GmEwGBRPg/Um1tYg8lxSUhIyMzMbzLZISUmBTqdDSkqKnyIjckzxYMaoUaMwduxYLFu2DOXl5Thw4AA+/PBD3H///UhPT/dFjE2Skgbflw0zG30iIqLg89hjj+Gbb77B/v378eOPP+K2225DVVWVx9NUvMXVxRh/XoRx1n9hbQ0i33N30IPI1xQPZrz44osYOXIk7rrrLsTHx8NgMOCee+7BbbfdhtmzZ/sixibJusFvaEDBmw2z/XOx0SciIgo+Bw4cQHp6Oi644AKMHDkSLVq0wA8//ACDwRDQuAKZfWHdh3HWfwmG7BAiIvIPxTUzJNXV1di7dy+EEEhISIBWq/V2bEHBHzUzsrKyUFlZCb1ej+zs7Hp1MrxZO8N+PivrchARNT2smUFA0/seWPdhjEYj+y9ERE2QT2tmSLRaLXr16oXevXs32YEMX3J2RcHRdimLA4DbU0KcZXvYX7FgXQ4iIiLvue+++3D8+PF620+ePIn77rsvABF5lz+mp7rTh2H/hYiIFGdmnDx5EiaTCevXr8fhw4dRV1dn8/i+ffu8GmCg+SMzA4DDf9s30EqqhId6RXFmjBARea6pXZEPJRqNBgcPHkT79u1ttv/111/o2LEjzp4967dYfPE9cNa/8GUWKQU39tmIyJt8ujRreno6vvnmG9x5552IjY2FSqWyeXzChAnKIw5iwdQhVNJYhHrDwo4MEZHngqntai6qqqoghEDbtm2xe/dutGvXTn7MYrHgf//7H4xGI/744w+/xuTt74Gz/oU32+1Q78M0N+yzEZE3KWq7hEI6nU7k5+crPSxkmc1mAUCYzWa/Pu/ChQuFwWAQCxcuVPRYU9EcXiMRka8Equ1qzlQqlVCr1U5vGo1GPPvss36NyZ/fA3fbbbbvTQ8/UyLyJiVtl+LMjK5du+Lzzz/HhRde6OlgS0jx9dUt66sPGzduxPLly5GWloa8vDy5KGhFRYXNMRwBJyIiV5iZ4X/ffPMNhBC49tprsXLlSuj1evmxFi1awGAwoFOnTn6NyR/fA6VZFOzDEBGRKz6dZvLee+/hv//9L955551mUfjT1x0B60b9wIEDsFgs0Gg00Ol0TgczmH5JRESucDAjcEpLSxEXFwe12uMa617jj++Bu4MTUt8lOTkZBQUF7MMQEZFDStquMKUnnzt3Lvbu3YsOHTogPj4e4eHhNo8XFxcrPWWzZr20mHVmxoABA2yKglrLyMhgB4CIiCgIGQwGHDt2DG+99RZ27twJlUqFnj174r777oNOpwt0eF5n3Y9xRVqt7fjx44iKivJTdERE1JQpzsyYOXOmy8enT5/eqICCTTCmaPr7fEREFFqYmRE4RUVFSE1NRcuWLdG/f38IIVBUVIRTp05hzZo16Nu3r99iCdT3wFE/RNp2/PhxVFZWOszkYP+FiIh8Os2kuQlUiqazBt2dhp7zUYmImjcOZgTOVVddhYSEBCxatAhhYX8nwJ49exb3338/9u3bh40bN/otlkB9D1z1Q1z1Y9h/ISIiJW2XTyZ0cnxEGaPRCIPBIKdo5ubmYty4cSgtLYXJZLLZV0rTtN/u6nxERETkH0VFRZg6dao8kAEAYWFhmDJlCoqKigIYmXfk5uYiPj4eubm5Tvdx1Q/JyMhASUmJwwsy7L8QEZESbg1mXHjhhVi6dClOnz7tcr/du3fj4YcfxuzZs70SXHMwZswYjBs3DsnJyXLDbjKZ5EKg9g26Ow29q44CERER+U50dDTKysrqbS8vL28StSLcuajiaT+E/RciIlLCrcGMV199FfPnz0eHDh0watQozJkzB++//z5WrlyJN998E5MmTUL//v2RmJgInU6HRx55xNdxNwm5ubn44IMPYLFYsGzZMnm7NGCxYMECNuhEREQhZNSoURg7diyWLVuG8vJyHDhwAB9++CHuv/9+pKenBzq8RmP2BBERBQtFNTMKCgqwbNkybNy4ESUlJTh16hTOOeccJCYmIjU1Ff/617/Qpk0bH4brf76cbyrNDQUAlUqF1157rcHBC84nJSKihrBmRuCcPn0ajz/+OHJzc3H27FkAQHh4OB5++GGYTCZERET4LZZgKGLOop5ERKQEC4B6ka86Arm5uZg8eTJOnToF4O86I9IAhauGn50CIiJqCAczAq+6uhp79+6FEAIJCQnQarV+j8HX3wOpxpfFYpGzNez7KLwIQ0RESgS8ACg1zGQyobq6GkIItG3b1iZlU5qPmpWVVa/Ilv18UleFuNwp0uWNY4iIiMiWVqtFr1690Lt374AMZPiS1FfIysqyqfHlqJ6G0WiEXq/H8ePHbfoW3u6/EBFR88PMjAb4MjMjKysLAJCdne1w+VVXa7FLXF3x8ORqCK+gEBGFPmZmBM7JkydhMpmwfv16HD58GHV1dTaP79u3z2+x+Op7IPUV9Ho9oqKi5EwMZ9mjjvoW3u6/EBFR08DMjBCQkZGBiooKVFRU1JsuImVfZGdnN1hky1UhLk+KdLGwFxERkefuv/9+vPXWW7jqqqswbtw4TJgwwebWFEh9hezsbJtsUWerkTjqW3i7/0JERM0PMzMa4O2rGkVFRcjPz0dKSgqSkpJsHlNaD4P1M4iIyBFmZgROmzZt8Nlnn+HKK68MdCh++x5Y90cAsO4XERF5jAVAvcjbHYGcnByYzWbodDpkZmbaPKY0rTImJgaVlZXQ6/WoqKhodGxERNQ0cDAjcLp27YrPP/8cF154YaBD8UsBUOtpsRqNBjqdTv63/RLz1v0cR8VCiYiIfDrNpLi4GNu3b5fv//e//8WIESPwxBNP4PTp08qjbWZSUlKg0+mQkpJS7zGmVRIREYW2WbNm4emnn0Z1dXWgQ/E5qeAnAGg0GlgsFpt/WxcCBWz7OY6KhRKR9xQVFSEnJwdFRUWBDsWrmurrIs8oHsx46KGH8NtvvwH4u4jV6NGjodVq8dFHH2HKlCleD7CpSUpKQmRkJG677bZ6VbqdzTW1J1X5Tk1NleesEqufExFR4M2dOxd5eXno0KEDevXqhb59+9rcmhLr2hkLFiyo92/7izPW/RxewCHyrfz8fJjNZuTn5wc6FK9qqq+LPKN4molOp0NxcTG6d++O2bNn46uvvkJeXh6+++47jB49GuXl5b6KNSB8kaLprAq40uNZ5dsW3xcior9xmkngzJw50+Xj06dP91Mk/vkesA4GUXByVacvlDXV10X/8GnNjOjoaGzatAnnnXcehgwZghtvvBETJkxAWVkZLrjgApw6dapRwQcbX3QElCy96up4dhxs8X0hIvobBzMI8M/3wNGFBGftMdtpIiJqiE8HM6699lrExcVh8ODBGDt2LH755RckJCTgm2++wd13393kroj7siPARp2IiHyBgxnBTQgBlUrl8+fxdR8mKysLNTU1iIyMlKe8urpYwwxKIiJqiE8LgObk5KC4uBjjxo1DVlYWEhISAAArVqxAcnKyZxE3M1Jth40bN+L48ePIyspinQciIqIQdeGFF2Lp0qUNFkLfvXs3Hn74YcyePdtPkflOVlYWKisrUV1djaioKGRkZNgUBHVUDyM5ORkajYb9RSIi8grFgxm9e/fG9u3bYTabbeZ9zpkzB++8845Xg2uqpMZ++fLlqKysRGVlpdNq3g0VtXSn6CULYxIREfnOq6++ivnz56NDhw4YNWoU5syZg/fffx8rV67Em2++iUmTJqF///5ITEyETqfDI488EuiQG62mpkb+tzQ4YV0Q1FFB84KCAlgsFhQUFABgH4aIiBpH8TQTADh27BhWrFiBvXv34vHHH4der0dxcTE6dOiAc8891xdxBowva2YkJycjLy8PAJCdne1wqklDKZnupGwyrZOIqHnhNJPAKCgowLJly7Bx40aUlJTg1KlTOOecc5CYmIjU1FT861//Qps2bfwWj6++B7m5uXjkkUcgdSGl/kVD02etHweAcePGwWKxsA9DREQyJW1XmNKTb9u2DYMGDUKbNm1QUlKCBx54AHq9HqtWrUJpaSneffddjwNvLjIyMtyukSGtxe5s6bLk5GQcOHDAZcpmQ+cgIiKixktOTm4WUyhMJpM8kKFWq+X+hZR5KmWb2g9sWPd/4uPjYbFYoNFoXPZP2IchIiJnFE8zmTRpEu69917s3r0bkZGR8vZhw4Zh48aNXg2uuXGUSmm9Jrsj9imbjjR0DiIiIiJ3GY1G6PV66PV6vPrqq3L/QppmIg1AWA9s2JPqZ6Slpbnsn7APQ0REzigezCgsLMRDDz1Ub/u5556LQ4cOeSUoR+Lj46FSqerd/v3vfzvcf8OGDQ73//XXX30Wo7uczf9sqOF3hMW0iIiIyJ8yMjJQUVGBiooKm0EG64EH64ENR9y5GENEROSK4sGMyMhIVFVV1du+a9cutGvXzitBOVJYWIiDBw/Kt7Vr1wIAbr/9dpfH7dq1y+a48847z2cxusvZoEVDDb8j7AwQERGRPxQVFSEnJwdFRUX1HrO/UOMqoyI3NxfHjx+HXq/n9BEiIvKY4sGM4cOH45lnnsGZM2cAACqVCmVlZTAajbj11lu9HqCkXbt26Nixo3z79NNP0b17dwwcONDlce3bt7c5TqPRuNy/trYWVVVVNjdvczZoIV3JMJlMblft9mQAhIiIiEip/Px8mM1m5Ofn13tMSXapyWRCZWWlvKQrBSeuJENEwU7xYMaLL76II0eOoH379jh16hQGDhyIhIQEREVFITs72xcx1nP69Gm89957uO+++6BSqVzum5iYiNjYWAwaNAhff/11g+d+/vnnodPp5FtcXJy3wpa5ulph3RlwpxHhXFIiIiLyh5SUFOh0OqSkpNR7TMnFFV6ICQ2eTH8moubHVdaerykezIiOjkZ+fj5WrlwJk8mEcePG4fPPP8c333yDVq1a+SLGelavXo1jx47hnnvucbpPbGws3njjDaxcuRIff/wxLrjgAgwaNKjBIqXTpk2D2WyWb+Xl5V6O3jUlxbOIiIgouBQXF2P79u3y/f/+978YMWIEnnjiCZw+fTqAkTWeqw6ruxdXGlq+lYIHB52IyB2usvZ8ToSgoUOHihtvvFHxcTfeeKO46aabFB1jNpsFAGE2mxU/nyOFhYXimWeeEUOHDhULFy50ue/ChQuFwWAQCxcutPm3q/38LZDPTUREjnm77SL3JSUliRUrVgghhNi7d6+IjIwU6enpIiEhQUyYMMGvsXj7e2AwGAQAodfrPW77pXMYDAYhBPsRREShrrCwUMyfP18UFhZ65XxK2i6VEP+/ULgC69evx/r163H48GHU1dXZPPb22297YYjFudLSUnTr1g0ff/wxhg8frujY7OxsvPfee9i5c6fbx1RVVUGn08FsNiM6OlppuPXk5OTAbDbj2LFjWLVqFUpKSty6ShEfH4/S0lIYDAaUlJS4/ZivBfK5iYjIMW+3XeQ+nU6H4uJidO/eHbNnz8ZXX32FvLw8fPfddxg9erTHGZ/PP/88nnjiCUyYMAE5OTluHePt74HUXzl+/DgqKys9avvt+zzsRxARkTUlbZfiaSYzZ87E0KFDsX79evz11184evSozc3XFi9ejPbt2+OGG25QfOzmzZsRGxvrg6jcl5KSArVajV9++QXJycmIj49HVlZWg9NJXKX6BTINkCmIRERE/xBCyBd61q1bh+uvvx4AEBcXh7/++sujcxYWFuKNN95A7969vRZnY6SmpkKv1+P48eOKi0PaT0dhP4KIiDymNO2jY8eO4t1331WeL+IFFotFdOnSRUydOrXeY0ajUdx5553y/fnz54tVq1aJ3377Tfz888/CaDQKAGLlypWKntOXqbp6vV4AEFqt1usplkzbJCJqvjjNJHCuueYacdddd4l3331XhIeHi927dwshhNiwYYM8tUKJ48ePi/POO0+sXbtWDBw4UNFUFV9NMzEYDPK/NRqN16fAsg9DRNR8KWm7FGdmnD59GsnJyd4cT3HbunXrUFZWhvvuu6/eYwcPHkRZWZl8//Tp03jsscfQu3dvXHXVVcjPz8dnn32GkSNH+jNkt0RGRnp9RRIWDyUiIvK/nJwcFBcXY9y4ccjKykJCQgIAYMWKFR71n/7973/jhhtuwODBgxvc19fLy1tnURiNRmg0GlgsFod9jcb0Q9iHISIidygezLj//vuxdOlSX8TSoKFDh0IIgfPPP7/eY0uWLMGGDRvk+1OmTMGePXtw6tQpVFZW4ttvv5VTPQNNWnI1NTUVBoPBZklb6bExY8Y0am1vpm0SERH5X+/evbF9+3aYzWZMnz5d3j5nzhy88847is714Ycfori4GM8//7xb+/t6eXnrKSIZGRlYsGCBT6bAenqsO0vaE/kDv4tE/qG4AOiECRPw7rvvonfv3ujduzfCw8NtHp83b55XAww0bxfPKioqwrJly7BhwwZs3rwZCxYssMnIkAphSVc7WBCLiIiUYgHQwDp27BhWrFiBvXv34vHHH4der0dxcTE6dOiAc889161zlJeXIykpCWvWrMEll1wCALj66qvRp08fpwVAa2trUVtbK9+vqqpCXFycX74HwbDkKouJUrAI9HexqKgI+fn5SElJQVJSkt+fn6gxfFoAdNu2bejTpw/UajV+/vlnbN68Wb5t2bLF05ibjfz8fLRu3RopKSkOUzOlqxFpaWnMrCAiIgox27Ztw3nnnYfZs2fjxRdfxLFjxwAAq1atwrRp09w+z6ZNm3D48GFceumlCAsLQ1hYGL755hu8/PLLCAsLg8ViqXdMREQEoqOjbW7e5OpqczBMDWFWKgWLQH8X8/PzYTabkZ+fH5DnJ/KXMKUHfP31176Io9lISUlBfn4+oqKiXP6RGzBgQMCm8xAREZFnJk2ahHvvvRcvvPACoqKi5O3Dhg3DmDFj3D7PoEGDsH37dptt9957L3r06IGpU6dCo9F4LWZ3WQ9Y2GdfGI1GOTMjUKTpL0SBFujvovR7IyUlJWAxEPmD4mkmkj179mDv3r0YMGAAWrZsCSEEVCqVt+MLOH+n6sbExKCyshJ6vR4VFRU+fz4iImp6OM0kcHQ6HYqLi9G9e3dERUVh69at6NatG0pLS3HBBRegpqbG43M3NM3Enre/Bw1NJQmGqSZERBTafDrNpKKiAoMGDcL555+P66+/HgcPHgTwd2HQyZMnexYxeYTFhYiIiIJLZGSkw1VEdu3ahXbt2gUgIu+xLgDqiHXmRkN9FPZhiIiosRQPZkycOBHh4eEoKyuDVquVt48aNQpffvmlV4Nrqpw14NJ9vV6P7OzsBhv6YJifSkRERP8YPnw4nnnmGZw5cwYAoFKpUFZWBqPRiFtvvbVR596wYYPbWRmBYF0noKE+CvswRMGrqKgIOTk5KCoqCnQoRC4pHsxYs2YNZs+ejc6dO9tsP++881BaWuq1wJoyZw24yWRCZWUloqKikJGR0WBDb91p4BUOIiKiwHvxxRdx5MgRtG/fHqdOncLAgQORkJCAqKgom6XYQ9G8efPw+OOPO125zjpzo6E+CvswRMGLBUQpVCiumREVFYXi4mKcd955NnNBCwsLcd111zW5Og++WJr1iy++QH5+Pm655RabVE37uaZK5p4GegkoIiIKHqyZEXhfffUViouLUVdXh759+2Lw4MF+j8Hb34PHH38crVu3xokTJzBnzhwA7tXJaKiPwj4MUXDh0q4USEraLsWDGTfccAP69u2LWbNmISoqCtu2bYPBYMDo0aNRV1eHFStWNCr4YOPtjkBOTg7MZjNOnDiB7t27e61AFotuERGRhIMZBHj/ezBv3jwcPHgQsbGx0Gq1MJlMOH78OCorK10ORLBwKBERucungxm//PILrr76alx66aX46quvcPPNN2PHjh2orKzEd999h+7duzcq+GDji8yMZcuWYcOGDThy5AivQBARkddxMCOw1q9fj/Xr1+Pw4cOoq6uzeeztt9/2Wxy+/B5Iq69ptVq0a9eOAxFEQYgZFhSKfLqaSc+ePbFt2zb0798fQ4YMwcmTJzFy5Ehs3ry5yQ1k+Ep0dDT0en2Da7FzDikREVFomTlzJoYOHYr169fjr7/+wtGjR21uTU1kZKRcJ4P9FqLgwtoX1NQpzsxobnw1zUSn0yEzM9PmMfs0S84hJSIiTzAzI3BiY2Pxwgsv4M477wx0KF7/Hkj9lOTkZOTl5QEAsrOz5YwMqd+i1+sRFRXFbI0gwqk8zRMzMygU+XSaybZt2xyfSKVCZGQkunTpgoiICCWnDGq+mGZi/UclNzcXWVlZ8uPW807Z8BARkSc4mBE4MTEx+Omnn4IiW9Xb34PU1FT07NkT3333HQoLC+tdbJH6Le7U0SD/4gUyIgoVPh3MUKvVUKlUAADpUOk+AISHh2PUqFF4/fXXERkZqTT2oOPrDqHUuADglQwiIvIKDmYEztSpU9G6dWs89dRTgQ7F69+DWbNmoa6uDmfOnMF7773Hgp4hhJ8JEYUKn9bMWLVqFc477zy88cYb2Lp1K7Zs2YI33ngDF1xwAZYuXYq33noLX331FZ588kmPX0BTl5ubi5iYGMTExCA5ORl6vR56vR7Z2dnyvFP7/TkHlYiIKPjV1NRg3rx5GDhwIB599FFMmjTJ5hbKhg0bBp1OhxEjRjjsr4wZMwZhYWHYuHGjw8cpcDIyMviZEFGTozgzo3///pg1axZSU1Nttufl5eGpp57CTz/9hNWrV2Py5MnYu3evV4MNBF9c3VKajcHUQCIiUoKZGYFzzTXXOH1MpVLhq6++8lss/v4ehIWFwWKxQKPR4OzZsz5/vsZgpoL/8T0nInf4NDNj+/btMBgM9bYbDAZs374dANCnTx8cPHhQ6ambjeTkZKjVami1WgBAaWkpTCaT/Lh9JobRaITBYGhw9RMiIiIKrK+//trpzZ8DGb7iLFs0NzcXERERUKvVSExMDPqMUpPJVK//Rb7F95yIvE3xYEaPHj1gMplw+vRpeduZM2dgMpnQo0cPAMDvv/+ODh06eC/KJqagoAB1dXWora1FampqvYEK+z/2TA0kIiIKLXv27EFeXh5OnToF4J86Y6HO2Q9Sk8mE6upqxMXF4ciRI0H/o5UXinzL0aBXIN5zTtUmatoUTzMpKCjAzTffDLVajd69e0OlUmHbtm2wWCz49NNPcfnll+M///kPDh06hMcff9xXcfuNL1I0c3NzMW7cOFgsFodTR4ItDS/Y4iEiItc4zSRwKioqkJaWhq+//hoqlQq7d+9Gt27dMHbsWLRp0wZz5871Wyze/h7MmzcPv//+O4qKipCenm7TJ7DuKwBgv6GZC5Yp0sESBxG5z6ermQDAiRMn8N577+G3336DEAI9evTAmDFjEBUV5XHQwcpXHcIxY8Zg+fLlSEtLw9KlS712Xl9gQ0BEFFo4mBE4d911Fw4fPow333wTF154IbZu3Ypu3bphzZo1mDhxInbs2OG3WLz9PXj88cfRunVrnDhxAnPmzPFChNRUBcuFsGCJg0JPUVER8vPzkZKSgqSkpECH06woabvCPHmC1q1b8w9CIxUUFMBisWD58uUYMGBAUL+fRqPR5moLERERObZmzRrk5eWhc+fONtvPO+88ufh3qIqJicGpU6cQExMT6FAoyGVkZARF3zZY4qDQk5+fD7PZjPz8fA5mBDHFNTPIO4xGIzQaDSwWi+I5pf6e/8eaHURERO45efKkXODb2l9//YWIiIgAROQ9kZGRUKvViIyMrPeYfd+EtQqImoaioiLk5OSgqKgo0KH4VUpKCnQ6HVJSUgIdCrnAwYwA0ul00Ov1ijMeWA2aiIgoOA0YMADvvvuufF+lUqGurg5z5sxxuWxrKLDv3Ofm5iImJgYxMTHIysqy6Zv4oq/CARIi/7POUGhOkpKSkJmZyayMIMfBjACQCoBWVlYiKirK7YwHqRFv164dNBoNkpOTPXpudgSIiIh8Y86cOXj99dcxbNgwnD59GlOmTMHFF1+MjRs3Yvbs2YEOr1HsO/cmkwmVlZWorKwEAJuVKnyxcgUv5hD5HzMUQlNzyahxazDj5ZdfRk1NDQCgrKysySwvFigmkwkWi0UekHB3cEFqxDdv3gyLxYKCggKPnpsdASIiIt/o2bMntm3bhv79+2PIkCE4efIkRo4cic2bN6N79+6BDq9RrDMxcnNzYTQaodfrodfrkZ2dbTMlNSMjQ6655a0LKFxOlcj/mKEQmppLRo1bq5mEhYXhjz/+QPv27aHRaHDw4EG0b9/eH/EFnK+WZpUKakqDC1qtFrW1tS5XN5GOS05ORkFBgUeVmVnVmYio6eNqJgR4/3sgrW4G/JOF4apPwdXQiIgCI5RXY/H60qxdunTBtGnTcP3116Nr164oKirCOeec43TfpsTXHUJpcKG8vBx1dXXQaDQ4e/as28d5Y1CCAxxERE0LBzMCZ9u2bQ63q1QqREZGokuXLn4rBOrt70Fubi4mT56MU6dOYfTo0SgoKKg3WGHdpwDg8/4F+zDNAz9noubD64MZb7zxBh599FGXP7KFEFCpVLBYLMojDmK+6hDa/1EeM2YMli9f7jIzw5o3r3bwygkRUdPCwYzAUavVUKlUACBPy5XuA0B4eDhGjRqF119/3eGqIN7ki++BdZ/BUWaGv/sU7MM0D/yciZoPJW2XWzUzHnzwQfz111/YunUrhBBYu3YtiouLbW6bN29GcXGxV15AcyBV/c7KyvLoeG/OG+UcVCIiIu9YtWoVzjvvPLzxxhvYunUrtmzZgjfeeAMXXHABli5dirfeegtfffUVnnzyyUCH6hGpTsbx48cBoN7S7f7uU7AP0zzwcyYiR9zKzLD2zjvvYPTo0SG/Vrq7fHV1KyYmBpWVldDr9aioqEBYWJhcFNRZBoy/UzeJiCg0MTMjcPr3749Zs2YhNTXVZnteXh6eeuop/PTTT1i9ejUmT56MvXv3+jQWX9X9GjduHCwWC6+SExGR13k9M8Pa3XffjYiICGzatAnvvfce3n//fWZkeCA7OxsGgwHZ2dkAgLS0NGg0GqSlpdXbV1pO1XoNd65KQkREFHy2b98Og8FQb7vBYMD27dsBAH369MHBgwf9HZpXSCuyAfBoiXgiIvpHc1lC1VcUD2YcPnwY1157Lfr164fx48dj3LhxSEpKwqBBg3DkyBFfxNgkZWRk2KRmDhgwAJ07d8aAAQPq7SsNXAD/VA9nuh0REVHw6dGjB0wmE06fPi1vO3PmDEwmE3r06AEA+P3339GhQ4dAhdgoRqMRGo0GABQvES9dnPHWUq1ERKGuuSyh6iuKBzMeffRRVFVVYceOHaisrMTRo0fx888/o6qqCuPHj/dFjE1KUVGRnH6am5vrMOvCnjRwYb2GuzQYAoAdAyIioiDx6quv4tNPP0Xnzp0xePBgDBkyBJ07d8ann36KhQsXAgD27duHRx55JMCReiYjI0POJlWamWGdVWo9sMFBDiJqatzNuEhJSYFOp0NKSoqfImtaFNfM0Ol0WLduHfr162ez/aeffsLQoUNx7Ngxb8YXcN6eb5qTkwOz2Yxjx45h1apVOH78OCorK6HVatGuXTvFNTCk6s56vR5RUVGsoUFERKyZEWAnTpzAe++9h99++w1CCPTo0QNjxoxBVFSUX+PwxdKsJpNJ7ru4WzNDOi45ORl5eXnydukcALhSBTUKl26lYCP95tPpdMjMzAx0OCHFpzUz6urqEB4eXm97eHg46urqlJ6u2UlJSYFarUZxcTGOHz+OmpoaAEBkZGS9iuDuXKmQsjYAKKqh4c65eaWEiIhIudatWyMjIwPz5s3D/Pnz8dBDD/l9IMMX7Ke9Jicnu9VPkI4rKChAVFQUKisr5XMonTrrj74J+z+hh7XkKNgw48JPhEI333yzGDBggPj999/lbQcOHBADBw4UI0aMUHq6oGc2mwUAYTabvXpeg8EgAAi9Xi8MBoNYuHCh032kx53tJ4Ro8HFX527MPkREFHx81XZRaPH298C+ryH1E1QqldDr9W71UdLT04VGoxHp6ekexaCkb6K0b+TJc1Bw8PSzJqLgo6TtUjyYUVZWJhITE0V4eLjo1q2b6N69uwgPDxd9+/YV5eXlHgUczHzVIXT2R9d6u/W/vd2wuvNHnw0DEVFo4mAGCeH7Pow0MAHA7T5KY/szSvomnj4X+z9E5InCwkIxf/58UVhYGOhQQpqStkvxNJO4uDgUFxfjs88+Q2ZmJsaPH4/PP/8cmzZtQufOnRubKOLUjBkzoFKpbG4dO3Z0ecw333yDSy+9FJGRkejWrVtIpAtKaXJZWVk2c/+8vXqJ/Woqnu5DREREzYtUtDwvLw8LFiyQpx+3a9euwWMb259R0jfx9LnceQ5ORamP7wk1d1yZxP8UD2ZIhgwZgkcffRTjx4/H4MGDvRmTUxdddBEOHjwo36T12h3Zv38/rr/+elx11VXYvHkznnjiCYwfPx4rV670S6wNkQYtHnnkEWg0GowZMwaA5zUwiIiIKDBefvlluQZWWVkZhLLa6iEtIyNDrpm2efNmt/Y3Go3yiia+js1XF2VYo6E+vifU3LFOhv95PJgRCGFhYejYsaN8c3UFIDc3F126dEFOTg4uvPBC3H///bjvvvvw4osv+jFi56TlzIQQqKurw/LlywH80/BmZ2fbXE1gA0FERBScJk2ahKqqKgBA165dceTIkQBH5FtSHyU7OxsA5KVa09LS6u0rXa0fM2aMfNW+KfRpvJ0x2xTwPaHmLikpCZmZmUhKSgp0KM1GWKADUGL37t3o1KkTIiIicNlll+G5555Dt27dHO77/fffY+jQoTbbUlNT8dZbb+HMmTMOV2QBgNraWtTW1sr3pc6JtxQVFSE/Px8VFRUAIE+Zse8AZGRk2FxJkK5isIEgIiIKLp06dcLKlStx/fXXQwiBAwcOyJka9rp06eLn6LzPuo+Sm5uLgoICLFiwwGEGhDRwceDAAVgsFrkvE+p9Gvt+GvE9ISL/C5nMjMsuuwzvvvsu8vLysGjRIhw6dAjJycnyoIC9Q4cOoUOHDjbbOnTogLNnz+Kvv/5y+jzPP/88dDqdfIuLi/Pq65DmUl188cXQ6/UYPXo04uLiMGDAAJfHsX4FERFRcHryySeRmZmJbt26QaVSoV+/fujatavNLT4+Hl27dnX7nAsXLkTv3r0RHR2N6OhoXHHFFfjiiy98+CoaVlRUhFmzZiE1NRW5ubnIzc3FuHHjXGZZSFfr09LS5Kv27NMQEZE3hMxgxrBhw3DrrbeiV69eGDx4MD777DMAwDvvvOP0GJVKZXNfmsNqv93atGnTYDab5Vt5ebkXov9HSkoKTpw4gY0bNwIAli9fLncCcnNzERMTg5iYGLmTEBMTg1atWsnbPOGsIBMLNRERETXegw8+iL/++gtbt26FEAJr165FcXGxzW3z5s0oLi52+5ydO3eGyWRCUVERioqKcO2112L48OHYsWOHD1+Ja/n5+airq0PPnj1hMplgMplgsVig0WicZllIAxdLly5VPIDhqp/izz5Mc+0vNdfXTUQhROlSKZs2bRLbtm2T769evVoMHz5cTJs2TdTW1io9XaMMHjxYZGRkOHzsqquuEuPHj7fZ9vHHH4uwsDBx+vRpt5/DF8uaSUuZabVaAUBoNBqb5VcBCL1eb7PcGayWFmtoyTBn68DbL03GddSJiJomLs0aOEuWLBE1NTU+OXfbtm3Fm2++6fTxmpoaYTab5Vt5eblXvweFhYXimWeeEQMGDBB6vV4kJSUJjUYj0tPT6+3ryfKm7vZfGnrM25prf0nJ6+ZythSMuFRqaFLSh1E8mJGUlCRWrFghhBBi7969IjIyUqSnp4uEhAQxYcIExcF6qqamRpx77rli5syZDh+fMmWKuPDCC222ZWRkiMsvv1zR8/iiQyg1Dnq93uYP/8KFC4Ver5dvAIRKpRItWrQQKpVKaLVam0EPR43LwoUL5UGQhgY/lA6KEBFRaOBgRuAVFRWJ//znP+K9994TmzZtatS5zp49Kz744APRokULsWPHDqf7TZ8+3eYiiHTz9vdA6ofY9zcc7SP1IxrqTyjpvzT0mLfZP1dz6R8peZ3NdcCHgtv8+fPFjBkzxPz58wMdCing08GM6OhosWfPHiGEECaTSQwdOlQIIUR+fr7o3Lmz0tO5bfLkyWLDhg1i37594ocffhA33nijiIqKEiUlJUIIIYxGo7jzzjvl/fft2ye0Wq2YOHGi+OWXX8Rbb70lwsPD5YEYd/miQyg1Dunp6fLAhX0DmZ6eLv/XOkOjoU6BdQejsY0sGyYiotDEwYzA+fPPP8U111wjVCqVaNu2rWjTpo1QqVTi2muvFYcPH1Z0rm3btolWrVoJjUYjdDqd+Oyzz1zu7+vMDIl9X8XRj3zrf7vTn/Bm/8VXpNciXXBi/+gfzWWAh0ILMzNCk08HM6KiosRvv/0mhPh7mkdOTo4QQojS0lIRGRmp9HRuGzVqlIiNjRXh4eGiU6dOYuTIkTZXJ+6++24xcOBAm2M2bNggEhMTRYsWLUR8fLxHf2C93SEsLCwUWVlZol+/fg6nkNhfyZAGMtRqtc2ghxCOGw5vNiZsmIiIQhMHMwInLS1NXHrppeKXX36Rt+3YsUMkJSWJ0aNHKzpXbW2t2L17tygsLBRGo1Gcc845LjMz7Pnqe+Csf+Aq89SdzIxg73NIr89Rn4yIiLzDp4MZ11xzjbjrrrvEu+++K8LDw8Xu3buFEH8PHDTFEWpvdwSkdKfMzEx5IEOtVsvzTR1dyXB2lYKZE/4VCh0tIiIhOJgRSNHR0eKnn36qt/3HH38UOp2uUeceNGiQePDBB93e31ffA6n/AUAkJSXJ2x1lLvi67fT3dBNXU2uIiBxhhogyStouxauZ5OTkoLi4GOPGjUNWVhYSEhIAACtWrEBycrLS0zU7KSkpMJvNyM/PBwBotVrU1dWhoKAAubm58rrrGRkZ8nJmztZulx5v7DrtrFbtHpPJ5HL5OSIiorq6OoSHh9fbHh4ejrq6ukadWwiB2traRp3DG6z7e0VFRfK/pZVLsrOz5f6Jr9tOf7bNGRkZWLBggVf6Xs0J+5nU3OXn59v8/iMv8tYIyqlTpxStEhIqfFUzQxrZt07FDFSmBTM83MPMDCIKFczMCJybb75ZDBgwQPz+++/ytgMHDoiBAweKESNGuH2eadOmiY0bN4r9+/eLbdu2iSeeeEKo1WqxZs0at8/h7e/B3LlzxWOPPSYGDBjgMDPDESVtpzdWQGnMucg32M+k5o6ZGcr4NDNDcvr0aRw4cABlZWUoKyvD4cOHcfDgwcaPrjRxK1euxOHDh/Hss89Cr9ejpqYGx48fB+D/TAtpv+TkZMXP2xxH2aUrTo6yZIiIiABgwYIFOH78OOLj49G9e3ckJCSga9euOH78OF555RW3z/Pnn3/izjvvxAUXXIBBgwbhxx9/xJdffokhQ4b4MHrXDh48iNatWyMpKQkGgwHp6ek4cuQIcnNz6/ULcnNzERMTg8mTJ8v9nIa4m2Vh/VxSJqvJZLLpkzCbMnh4q39LFKqSkpKQmZmJpKSkQIfS5KiEEELJAb/99hvGjh2LgoICm+1CCKhUKlgsFq8GGGhVVVXQ6XQwm82Ijo5u9PmeeeYZCCFQV1eHZ599Vk451Wg0TqeTWLOfiuJIfHw8SktLYTAYUFJS4vRc9vu5c26lz0FERP7n7baLlFu7di1+/fVXCCHQs2dPDB482O8xePt78PDDDyMyMhI1NTW45JJLMG7cOFgsFhgMBhw/fhyVlZXQ6/WoqKiQ+wkSd/oL7vZD7PsgjvokSvo05Bm+x0TkC4raLqVpH8nJyWLAgAHi888/F5s3bxZbtmyxuTU13k7RXLFihXj66afFyJEj5eKfKpVKnnLSEHdS9dxNrbTfT0kaINM3iYiCF6eZkBDe/x5Y9xPsi5RLRT+lVT70er3QarWiRYsWQqVSOSx07in7c/iiT8J+TsPsV+Dj+0VE3uDT1Uy0Wq3YuXOnR4GFIl8szfrMM8+I5ORkuXG37gA0xJeNBRsiIqKmgYMZJIT3vwfW/YT09HSh0WhcrsZmvbKJ1McJlfoJoRJnIDn7zImIGsOnNTN69uyJv/76S+lh9P/WrFmDuro6XHPNNRBC4IMPPkC3bt1gMBiQnZ3d4PG+rNvAmhBERETkTFJSEv71r3/h7bffxgcffACLxYK8vLx6+xmNRuj1ehw5cgRHjx6t95jBYEBycnJQ196yr/PQHGuFNcS63yh95sePH+d7RE3WypUr8cwzz2DlypWBDoX+n+LBjNmzZ2PKlCnYsGEDKioqUFVVZXMj186cOQMACAsLk7cVFxc3mUEENvZERERNU35+PsLDw3HllVfK22pqagDYFtzMyMhAVFQUqqurIYSARqORL9hIP4ALCgqCukCn/QWeUC8o6uv+mfSZV1ZWhux75GvsI3tfUVERcnJybJaI9qUdO3ZACIEdO3b45fmoYYoHMwYPHowffvgBgwYNQvv27dG2bVu0bdsWbdq0Qdu2bX0RY5OiVv/9llsXSo2MjAxUOF4X6o09EREROZaSkgKz2Yz8/Hx5m9SHSU5OhkajQXJyMoB/sjP0er3DAuehtsJFqMVrz75/5o0f1vbn8MZ71JR/8LOP7H35+fn1/ib50kUXXQSVSoWLLrrIL89HblA6h2XDhg0ub02Nt+ebPvPMM2LGjBniySeflAuAarVauWCWPVd1LIKxxkUwxkRE1NywZkbgbNq0SWzbtk2+v3r1ajF8+HAxbdo0UVtb69dYfPE9WLhwodBoNDbFP4UIrhoT7IvU15ii78744jMPpu+Rt/F76X2FhYVi/vz5orCwMNChkBcpabsUL83a3Hh7WbOZM2cCAOrq6vDMM89Ao9HIWRqOli1ztQQql0clIiJHuDRr4PTr1w9GoxG33nor9u3bh4suugi33HILCgsLccMNNyAnJ8dvsXjze1BUVIT169dDpVIhLCwMTz31FGpqajBq1CgsXbo0qJbpZP+oYd74vHzxmQfT94ial6KiIuTn5yMlJQVJSUluP0bep6TtUjzNBACOHTuGuXPn4v7778cDDzyA+fPnw2w2exRsc3PxxRcDAH799VeoVCqo1WpotVro9XqbtDwpzS45Odlpyl6opzwSERE1Nb/99hv69OkDAPjoo48wYMAALF26FEuWLAnponH5+fmoqanBqVOn8Pvvv6O6uhp1dXX44IMPEBMTAwAN1v9yNIXAF9MK2D9qmDeKvvuicDyL0VOguJqy4u/pLOQ+xYMZRUVF6N69O+bPn4/Kykr89ddfmDdvHrp3747i4mJfxNiklJeXAwBiY2MhhMCZM2dw6tSpevtJ8+oKCgps/qhbN/rSH3wATXZ+IRERUSgRQqCurg4AsG7dOlx//fUAgLi4uJBeDS4lJQVhYWEQQtTLdnC36KOjmgG+qCPAH8REpFRKSgp0Oh1SUlIUPUaBpXgwY+LEibj55ptRUlKCjz/+GKtWrcL+/ftx4403IjMz0wchNi1SBkubNm3kbUKIeh0B+6XLxowZg/j4eGRlZfmlI+ANTbmIExERkSNJSUl49tln8Z///AfffPMNbrjhBgDA/v370aFDhwBH1zgWiwUqlQpdunSBVqtFenq6XOTTnSwIRxkT7mRRsD9BRL6WlJSEzMxMh9NIXD3mTf5enaUp8CgzY+rUqTZLi4aFhWHKlCl8490gVf2urq6Wt6lUKmi1WpuG3H7psuXLl6O0tBQAQqYjEKyDLERERL6Sk5OD4uJijBs3DllZWUhISAAArFixQl7pIxTl5+dDCAEhBE6cOIHq6mrs3r0bFRUVqKiocCsLwlHGhDtZFN7uTwTz4Ih1bL6O093zB/P7pVRTei2k3MqVK/HMM88E7ZQ/TmdRTvFgRnR0NMrKyuptLy8vR1RUlFeCaspOnz4NAIiIiJC3CSEQGRnpsCGXBirS0tJgMBiQnZ3tUUdAyujIysry4qtxLVgHWYiIiHyld+/e2L59O8xmM6ZPny5vnzNnDt55550ARtY4Unq1SqVCp06dAMBvF7Hsl31trGC+2GIdm6/jdPf8nsYRjH28YP7syfd27NgBIQR27NgR6FAc4nQW5RQPZowaNQpjx47FsmXLUF5ejgMHDuDDDz/E/fffj/T0dF/E2KRI82g1Go3N9mPHjjn8Y5+RkQGj0YiCggJ5UCDYGgZnAnG1hYiIKBicPn0aBw4cQFlZGcrKynD48GEcPHgw0GF5LCkpCTqdDkIIecpseHi4R+dS+iO3oKAAFosFBQUFHj2fvWAuEGodm6/jdPf8nsYRjH28YP7syT2NmYpx0UUXQaVS4aKLLvJBZO5xFb+/prM0KUrXfa2trRXjx48XLVq0EGq1WqjVahERESEyMzNFTU2N0tMFPW+v0T5jxgwxY8YMMX36dKFSqYRarRYqlcrlmtrSmtt6vV5e213p+tvS2tbp6elBtcY119wmIvI+b7dd5L5du3aJlJQUuY8k3aQ235981Yd5+umnBQCRlJQkNBqNSE9Pr7evq/Zd6te425cJ1j6MLzWF/lFTeA0UfObPny9mzJgh5s+fH+hQPBLq8fuDkrZLJYQQngyCVFdXY+/evRBCICEhAVqt1kvDK8HFm2u0A8DMmTMB/D21ZM6cOaiuroZWq0W7du2crqmdm5uLrKwsHDt2DHV1ddBoNFiwYIFHVbq59joRUdPn7baL3HfllVciLCwMRqMRsbGxUKlUNo9fcsklfovF29+D5557DmfOnMHp06fx3HPPyds1Gg3Onj0r38/NzcW4ceNgsVgc9jdyc3NhMpmc9nucaU59mOb0Wn3F0+8ZBbeioiLk5+cjJSWl0RkMRUVF+OqrryCEwKBBg/ySEeHN+JsqJW2X4mkmEq1Wi169eqF3795NdiDDF3Q6HQDgxIkT8rbIyEin0zGkP8QAFA9kOErjDFR6XTDOmwwmfH+IiJqGLVu24PXXX8ewYcPQp08fXHLJJTa3UGaxWCCEwMmTJ222azQam/bLZDLBYrFAo9E47G+4Mw3V332YYGuHOR2i8aynuYwZMwZhYWEYM2ZMoMOiRmpoKoaSaSj5+fk4deoUampq/FZ0k1NJvMutzIyRI0diyZIliI6OxsiRI13u+/HHH3stuGDg66saWq0Wc+fOBYB6o8fWVzb0ej2ioqIUjS4H06h+MMUSjPj+ND8cmSdfYmZG4PTr1w/z588PigJu3v4ezJgxAyqVCkIIzJw5U/43AJv2S7oQk5ycjIKCAiQnJyMvLw8AkJ2d7VY/xt/tItth//N15oT1+aX+tH0WETU9OTk5MJvN0Ol0yMzMdLlvQ5kZ7KsFhtczM3Q6nZwmqdPpXN7ItTNnzgD4p2hWdXU1srKy5NVGrIskWV/ZcLSKSUOCaVQ/mGIJRnx/mh8uv0XUNM2ePRtTpkzBhg0bUFFRgaqqKptbKJP6glIx8/DwcLeXl6+srERlZaXbxSAdtYu+zJ5obDscbJkdoUBpgVCl77F1BlBaWho0Gg3S0tIaEzKFAFcrgthnbSQlJWHKlCmYOnWqPFhhvQ/7aiHAV4U7mgpfF8+SblqtVmg0GpGUlCQXS1q4cKHQ6/VCr9d7tXiSLwsyNXTu5ljEi8iRwsJCMX/+fFFYWBjoUKgJYgHQwFGpVHKxz6ZaAPSpp56y6cPo9XqH+1u3+d7ozygtHNpYzvo0jrb7OzYlgrUQp9I+YTC/xxQa3Cm+ab0P+2qBoaTtUjyYUV1dLU6ePCnfLykpEfPnzxd5eXlKTxUSfLmaSVJSkty46/V6AaDeaiX2f7i90SD5sjFo6NzS42q12mUHiIiIPMfBjMDZsGGDy5s/efN7UFhYKKZPny73YawHM9RqdYMXMaSLNI3pw/j7R7mzPo2j7cE6YCBE4/t9vv7c3I0vmN9jCg3uDE5wACPwfDqYMWTIEPmPyNGjR0X79u1F586dRWRkpHjttdeURxvkfDmYAUBezkzKwtBqtTZXLuz/cEt/8DUajeI/5v7IinA3M0Or1XIwg4jIRziYQUJ493sgXa207sNotVq/LC8vRGB+yCrJzPD2c3hTY5+jsYMhDR3PQQoisubTwYyYmBjx888/CyGEWLRokejdu7ewWCxi+fLlokePHsqjDXK+HsywTjl1p7FYuHCh3BnQ6/WK/vh7OyOjMY0PGy4iIt/hYEZgHT16VLz44oti7Nix4v777xfz5s0Tx44d83sc3s7MsO/DaDSaBi+QSBdrpIxMTy7GCNF0pxi4m4EbyH6T/XMrjSVU+nyhMLAUKpjdQI3h08GMli1bitLSUiGEELfffruYMWOGEEKIsrIy0bJlS6WnC3q+HszQarUNZkw4a0SkqSnuNuze/gPaVDsWREShjoMZgVNYWCj0er0499xzxS233CJGjBghOnfuLGJiYsSmTZv8Gou3vwdSdsbEiRMbzMiw7nM0JqvU0flClaPX4CwD152pLYESTLF4kz9eV1N97+y5U5uCyBklbZdbq5lYS0hIwOrVq1FeXo68vDwMHToUAHD48GEu/6aASqWCwWDA3Llz5WrOBQUFDlcskVY6mTx5MuLj4wEAJSUlyM7OVlR525113Z3xdL13VvcmIqLmZOLEibj55ptRUlKCjz/+GKtWrcL+/ftx4403NrhMYLCrrq4GALRs2RJCCKjVaiQnJ9dr66Wl5aWVKqT+woIFCzxegrMxfRhfc7ev42j1DvvX5axvFUyrnvk6lkD1Hf3xHgfT5+hLrlYUsWe/wojSx93hjXNQkFI6UvLRRx+J8PBwoVarxZAhQ+Ttzz33nLjuuuuUni7o+SozQ8poEaLhkXqpvoR0a2g019N0QFf7KR1J9jR7hIiIGo+ZGYETGRkpdu7cWW/7jh07/J7B6uvsUql9t+8jSPfVarXiFUw86cMEOmvD3T5SoOMMtjicaS7ZC6HMm9NIGsri8EaWh7cyRTh9xj98mplx2223oaysDEVFRfjyyy/l7YMGDcL8+fMbN7LSzDz88MMICwvDW2+9JW+TRqOlbAyTyYTIyEj5cY1GI4/mOhu5th/5l86VlZXlMh5X630rHUmWzgVAPo5ZGkRE1NRFR0ejrKys3vby8nJERUUFICLfUKvV0Ov1MBqNSE5OhlqtxpEjR5Cbmyv3Gdq0aYPKykqH/QpnPOnDuOq/+IO7faRAZ5c46mMGo+aSvRAqHGU15Ofnw2w2Iz8/v9HnbyiLQ0mWhyNFRUWora1Fy5YtPT6HxJuvm7xEySjJmTNnhEajEdu3b/d4pCXU+DIzIzMz0ybjQqvVyqPR1sU9pewMlUpls/qJs6rg9iPuUnZEQyuH+Lo6N0faiYj8g5kZgfPoo4+Kzp07iw8//FCUlZWJ8vJy8cEHH4jOnTuLCRMm+DUWb34PVqxYIS/N+vTTT8s1v4T4p323b+MdZVlIy9G7s5SrEO71YYI90yBYOOpjEjXEUVZDKGUoeLN+Ryi97lDm0wKg3bp1E1u2bPEosFDkyxTNpKQkm8EM/P9SrY4aGGcpnO4U0/JWI+9OJ8QfcRARkWsczAic2tpaMX78eNGiRQuhVquFWq0WERERIjMzU9TU1Pg1Fm9+D6wvxkjTTKQBhvT0dKFWq20GOByxHvRwtxioP/sO3uwv+TJmT8/v7+MCJdTiDXah/gM+1ONvjnw6mPH222+LYcOGiYqKCo+CCzW+HMzQaDRCr9eLFi1ayI279fxS6z/GjV0WyxucXXkh8hV2SHyDDXvTx8GMwDt58qTYtm2b2Lp1qzh58mRAYvDXYIaSmhHWy7QGW1+iMRmkjlZvcec8nrRz/s509efzeaPdZyYwUWjzac2Ml19+Gd9++y06deqECy64AH379rW5kfvUajUqKysRFhaG8PBwAEBdXZ08v9R6Dqj9PMuG5l36oj6F0WiEXq+X58g2R6z74V+BngfdVHHOJ5HvabVa9OrVC71794ZWqw10OF4lhIBer0dqairi4+ORnJzssMaBdZuZm5sLk8mE7OxsvPrqqy5rIoTiahbW7ZWS83jSznmzpoQ777U/a1h4o91v7jU3QmnljsbE6u3XGUrvG/1DJYQQSg6YOXOmy8enT5/eqICCTVVVFXQ6Hcxms1eWnpXePyGE/G+VSgXrj0GtVmPUqFEYMGCA3CgqLRYVHx+P0tJSGAwGlJSUNDpu+pvS91XqvHnyGRLfP18pKipCfn4+UlJSkJSUFOhwyAe83XaRayNHjsSSJUsQHR2NkSNHutz3448/9lNU3v0ePPfcczhz5ox8f/r06Q7bROu/21lZWaisrIRKpQLwd9/HnfbTH32YhtoXpe2Pp+1VoNu5YOsvBvr9aApycnJgNpuh0+kCvhx0Q/2NxsTq7dcZTO9bc6ek7VI8mNHc+GMwQ7piU1NTg8jISFRXV3vcqIwZMwbLly9HYmIijhw5ErDGoKk2RtavC0CDrzHYOglE1DxwMMO/7r33Xrz88suIiorCvffe63LfxYsX+ykq734PVq5ciZ9//lm+36FDB2RlZcl9l+zsbGRkZNi0e8ePH0dlZaV8jEajwYIFC5CRkVGvn5CbmyuvWJKamoqCggIkJyejoKDAJ32Jhtpnb7bfwdwnCubYyDPeuGDhrYseDQ0QNOZ5pGPj4uJQXl7e6Fh5oSd4KGq7PJnHcvToUbFo0SJhNBrl2hmbNm0SBw4c8OR0Qc0fa7SrVCpFlb5dkVY40Wg08vkCUXOgOcxXdOc1suYDEQUCa2aQEN79HsyePdumDxMeHl5vRTb7Pox0v0WLFg2uyOaoLpeSWhyu2lpHj7s6xvp1OCvMrkRz6BMFG/a/XGuodpa3VgDxR40ub8TKWmLBxac1M7Zt24bzzz8fs2fPxosvvohjx44BAFatWoVp06YpPR3h7ywNk8kkj44DkOtmKJ03mpaWBo1Gg7S0NADuzz1U+jwN7d8c5iu68xoDvaY8ERH516lTp1BdXS3fLy0tRU5ODtasWRPAqBpP2CXyWk85AYDq6mp5WklUVBQyMjKQkZGBiooKxMbGQgiBgoICAH/3TSwWCzQajdyGOqrLZd/OOut7NNTXcfS4q/bZZDLJr6OgoKDesUr7TI76C6zB5VuOPvNQes99Xb+hodpZKSkp0Ol0SElJadTzJCUlITMz0yuZDs7eE2/EylpiIUzpSMmgQYPE448/LoQQonXr1mLv3r1CCCG+++47n444P/fccyIpKUm0bt1atGvXTgwfPlz8+uuvLo/5+uuv6y19CkDs3LnT7ef1R2aGWq22qX5tvf53Y0fz3b0aofR5nO0fDCPhwbDyCxE5xysgvsfMjMAZMmSI3N4cPXpUtG/fXnTu3FlERkaK1157ze3zeNLvsefN78GKFSvE9OnTbVZkAyDatm1rk2mqVqtFenq6T1Zk87TvIT3ubpaFfex6vV5otVo5U8M+q8QTzNbwLW/0dQPJW5kRzjS2HQ5EO+7L94T9kuDi06VZo6OjxZ49e4QQtoMZJSUlIiIiQunp3JaamioWL14sfv75Z7FlyxZxww03iC5duogTJ044PUYazNi1a5c4ePCgfDt79qzbz+vrwQzrhtWTNEhvdQaUnsvZ/sHQUNjH4MuYAjFQwsEZCnW+7qQRBzMCKSYmRvz8889CCCEWLVokevfuLSwWi1i+fLno0aOH2+fxpN9jz5vfg+zsbIfTTOynm0jtrS/aXlcDJO4cp9frPYpJei34/2m80n8bMwVFaZ+PGi+U3l/rH9fB+EM7EO14ML4P5Bs+Hcxo3769KC4uFkLYDmbk5eWJzp07Kz2dxw4fPiwAiG+++cbpPtJgxtGjR90+b01NjTCbzfKtvLzcZ4MZ4eHhQqvV1qub4Yz9H2FPOwqedgaUnjtQ7K/AeGO+qzOBGCjx94BRMHym1LSwQ+J7HMwInJYtW4rS0lIhhBC33367mDFjhhBCiLKyMtGyZUuPz+tOv8eeN78HUv/FPrvU+uYsM8MXlLSF0r7Wma9KOKuh4e32OBguCFHwCcaBAyWDLWzzSSmf1swYPnw4nnnmGXmupEqlQllZGYxGI2699Valp/OY2WwGAOj1+gb3TUxMRGxsLAYNGoSvv/7a5b7PP/88dDqdfIuLi/NKvI6cOXMGp06dAvD3XNTJkyfXm8uXm5uLmJgYxMTEICsry2b+nzfqUjQ0z9TZ/EJn2xtTI8JbcxmlGKR5rgUFBT6rW+HL2iDOPht/1yPxxprvRNa8OYeWKNgkJCRg9erVKC8vR15eHoYOHQoAOHz4cKNWFHGn31NbW4uqqiqbm6+Fh4dDpVJBq9Wibdu2qKurQ0FBgU9qRln3E5S0hdK+2dnZHsUk1f6oqKjA0qVL5XN4uz1uDvXGSDlv1a9QoqEaEtbteEP7sh4F+ZQnIyVXXnmlaNOmjdBoNCIuLk6Eh4eLAQMGKEp9bIy6ujpx0003iZSUFJf7/frrr+KNN94QmzZtEgUFBeLhhx8WKpXK5VUNf2Zm4P+rf8PuyoZGo6mXfSHtq2SlE3eu7Dd05cTZlBR354s6O78/5jKGekZBsMQfLHEQkfuYmRE4H330kQgPDxdqtVoMGTJE3v7cc8+J6667zqNzutvvcZYx4Y3vwaxZsxxmZoSHh3ttRTZrgax5wHaPmjsl2RTMzCBv8+k0E8n69evFnDlzxOzZs8XatWs9PY1HHnnkEWEwGER5ebniY2+88UZx0003ub2/L2tmSIU/1Wq1w/mmQgibx6X0SOt/u2ponTX6nswzddSZcGe+qLMYnA2SsPNARNR4HMwIrIMHD4ri4mJhsVjkbT/++KOiAuTW3O33+PKCzMyZM51OM7G+OCLVpfB0SockkP0E636Ou88lFQaVlp9VSmmRUgp+7vyIl/aZO3dus/zcOdBBjvh0monk2muvxWOPPYYpU6Zg8ODBnp5GsUcffRSffPIJvv76a3Tu3Fnx8Zdffjl2797tg8g8J70OtVqN9PR0OcVQWqq1b9++0Gg0SE1NlVMQATSY+u8sXbGh1E/rNE5H+0rnXbBggcMly9yJwdF2LmPqW6G0JBkRUSg6e/YswsLC8NdffyExMRFq9T/drP79+6NHjx6Kz6mk3xMREYHo6Gibm7cIq6VZrf+t1WphNBrlKYkA3O6nuOLLfoKj9tB+CotGo4HFYnE7/uXLl8NisWD58uUexSS9f8uXL+fUzibCnekV0j4HDx5slp/7+vXrYTabsX79+kCHQiHKo8GM9evX48Ybb0T37t2RkJCAG2+8EevWrfN2bDaEEBg3bhw+/vhjfPXVV+jatatH59m8eTNiY2O9HJ1n6urqYDKZ8Pvvv8v3AcidAqlGxubNm2GxWJCXlweTyQSj0Yjs7OwG51V62ug3tDa39XmldeGPHz/u8Eeysxg4cOF/rH1BRORbYWFhMBgMsFgsjT6Xt/o9vqBSqaBWq5GUlISamhpMnjwZycnJNnUp7Psp9gMIDQ2w+7Kf4Kg9tN+m0+mg1+vdrl+RlpYGjUaDtLQ0j2KSBm/S0tKaVN2MoqIi5OTkoKioyOH9psxVrQvpfYiLi0PLli0RFRWFoUOHevy5O3tfg+39to9HpVLZ/JdIKcWDGQsWLMB1112HqKgoTJgwAePHj0d0dDSuv/56LFiwwBcxAgD+/e9/47333sPSpUsRFRWFQ4cO4dChQ3IBTQCYNm0a7rrrLvl+Tk4OVq9ejd27d2PHjh2YNm0aVq5ciXHjxvksTiVUKhWSk5NtGr5ly5bhkUceQWlpKWpqamwaNuCfqxxSIw/A61fbpQY1OTlZPrc0sJKVlWWzb0ZGBqKiolBZWYmsrCyPY2lM1oC7x/rjOYKZtwqLNYX3gpQJts4QUTB78sknMW3aNFRWVjbqPO70ewJJo9GguLgYdXV1qK6uRl5ensP9Nm7ciPj4+HpFzKXBg3HjxrnVnniz7XHUHlpvM5lMqKysRFRUlNuDKUuXLsXZs2exdOlSj2KS+nXWBUYB29ft6XvgjXN4yj47IRSLQXr6nrkqdi29D+Xl5WjRogWEEBg2bJjHg3fO3tdAv9/2/Qf7eK699lrodDpce+21AYmPmgClc1g6deokXnnllXrbFyxYIGJjY5Wezm1wUNAKgFi8eLG8z9133y0GDhwo3589e7bo3r27iIyMFG3bthUpKSnis88+U/S8vqyZIb2G9PR0kZ6eXq92hlRAS5pD19hiWErnmlqf23oOrLPzSvsomWOq9HW4+x54u3CYv4qOhQK+F81PIJaFo8ZhzYzA6dOnj2jdurWIiIgQ559/vkhMTLS5ucudfk9DfL00q3URc5VKZdMHsK47IfUfrNtlJcXEhfCs7XHUF3CntoUntcV8VevC+nV72v564xxCeFbfwP6YUKyR4It+j5KlTZWez53t3ji3O+z7DytWrBAzZ84UK1as8Dgeavp8WgC0devWYvfu3fW2//bbb6JVq1ZKTxf0/DGYAcCm4VepVDYDGa7+gCppcJUOfFhXJHenQrnSjoknr6MxAxeNKRzWXIqTuvM6m8t74alQ7Kg1pCm+pqaOgxmBY/2j39HNn3w9mLFw4cJ6q7FJba/UViQlJQmNRiOSkpLqXZxxNACgZBW0hjjqC0gxajSaRr8n1s/haf+nIQ1d1GrMOZSez9HAdnNoH4K93+PpZ+DucZ5e0CgsLBQmk0nMnj1bfg5eHCF3+HQwY8yYMeKFF16ot33OnDli9OjRSk8X9Hw5mCFdxbBv+N1p1F2RrjrYdxyULJcmZVmoVKp6y8TaN9TeaGjtNbYzE+wNT7DyxdWH5oYNNQUDDmaQEN79Hkh/26Q+jF6vF+np6XJfpm3btg77A1L2hpR9ap0d4GjFE6n/odVqPepbNNQnsc7M8EZfQUlmhtLn83VfRmmb7+jHb2PbvOYwGOJrnn4G7h7n6WfUXAe/qPF8Opgxa9YsodPpxPXXXy9mzZolZs2aJW644QbRpk0bMWvWLPHSSy/Jt6bA10uzqtVqodVqbdIupUbRPjPC1ai6dSMqXR2wnq5ifeVAyWCG9THOGlXrxrCxDa/9lBX+qPYvDgI1HhtqCgYczAiso0ePikWLFgmj0SgqKiqEEEJs2rRJHDhwwK9x+Dozw7q/IV38sG/HpcGO8PDweoMIjpZxtT+uoT6Mfbvl6Ad6Q/0XpdNjPW0rlQ4eNLR/YzM1vDE1prCwUMyePVuYTCa32z3rdpIXABrPUb9DybKwvuqvsD9EnvLpYEZ8fLxbt65du3oUfLDxxzQTrVYrXymQrmDYD0i0aNGiXqNufXXDOr1ROpf1MQ1NAbFvEBcuXGgzF9bVevHWx3raMZC4uloTDHw9N5aIyBs4mBE4W7duFe3atRMJCQkiLCxM7N27VwghxJNPPinuvPNOv8biy8EMKTPDut6XXq+XByGkzAqpvXR0kcLRBQz7jI6G+jD2P/gdtdPO6n5Zn1dJv0PpoIT96/VWZoancXh6nDNKBySs9+cPXt8I5kEifubUEJ8OZjQ3/hjMUKvVcmOqVqttMjPsC385Kwpq/+Pa0bQSVz/EHQ0iSNusi3w11PDZdzj8kVLpz2wCR4NHRETBhoMZgTNo0CDx+OOPCyH+rjMmDWZ89913fm8zfDmYIQ1idOnSRWi1WqFSqWwuglhPExHCdVvt6IKKs7pdjvo7rrIurAdLGlvEvDEXNHzRV3HnnN6awuyK0h+nvvox688fya6eKxh+rAdDDM4E80ALBQe/DmacPXtWbN68WVRWVjb2VEHJH4MZSUlJckcgKSnJpgF3VP/CmqMOgKMrII6yJ9y5OmJ9xcJZR8JeQ8/lbZ48RzCkZhIR+QoHMwInOjpa7NmzRwhhO5hRUlIiIiIi/BqLr6eZWNf9su8ruDNd1Lp+hfW/7dt16/sNtfmO2mlHF3gcHecoy9V61RPrCxqOpt+66lv4oz/kiKs+n6M4A/Ej2FvP6esaEO4+Vyj8WA/kKifBPNBCwUFJ26WGQpmZmXjrrbcAABaLBQMGDEDfvn0RFxeHDRs2KD1ds7dw4UIcOXIE1dXV6NKlC44cOYLKykpUVlYiKysLBQUFWLBgAQoLC1FSUgIANmuFjxs3Tl6nfcyYMfJ9ADAYDEhOTq63vruj9dUlqamp8mMZGRlYsGABDAYDsrOzERUVhcrKSixfvhylpaXIyspCfHw8xowZ43T9bem5pDg8Xdfc1Rrfrl6PM9L69tJ69+5ytg48ERERAERGRqKqqqre9l27dqFdu3YBiMg3VCoVAKBFixZyO5ydnY2KigpUVFQgOzu7wbZ5+fLlsFgsWL58uc2/7fsOycnJ8rkctfnWfQSpv/Phhx8iOTlZbuelPox9uy8dC0Du8xiNRpt4xowZg7CwMLRr1w4ajQYWiwUmk0l+rqysLAD/9C2k/pF1n8WTvor9a/OEo+d11QfKz8+H2WxGfn6+R8/nCW89Z0pKCnQ6HVJSUhr9fEVFRcjJyUFRUZHi53I3jkDyxnsunePzzz9HUVFRg++ZJCkpCZmZmUhKSvL4uYlkSkdKzj33XHkkbdWqVaJTp05i165dIisrSyQnJysfeglyvszMaNGiRb1MB6kyuEqlsql5IXF0dUK6WS83Zj9HVZo+4iyjwJ0rBvZXO+ynwTi6gmJ/fvx/JoqSlVXcjU8JFrokoqaMmRmB88ADD4gRI0aI06dPi9atW4t9+/aJ0tJSkZiYKCZMmODXWHydmSFllFr3DawLmFtvs86SsO7zOMrMcJQt6ihTwtF+1n0TaSqMdb+ooSLmer1eaLVaER4eLk+dsT6Pdfz201d8UcTcFxkdwZqZsWLFCr88tzuvMRSyKxrDW5kZM2fOlN+npv6ekf/4dJpJRESEKC8vF0L83WBLDfO+fftEVFSU0tMFPX9MM3G0/Km03X4wIykpSQB/L4EmrYTiaA1365RP6wbLUcMrNdxKBhekY6UGXq1Wu0zlXLjQdi166eZuwS1/zPkMBk3hNTWF1+BtTKkkf+NgRuCYzWZx5ZVXijZt2giNRiPi4uJEeHi4GDBggDhx4oTfY/HlYIZ1X8X+wor1YIKjlUmkH+eO2gxHfRjrCzbSMeHh4QL4e6UU64EU6xVU7C+42HM0Ndb+tWi12noDKc5id7XdE/5qU73VTjXmPMH0Y7ih18F2/W/W7wPfE/IWnw5mdOnSReTl5YmzZ8+KuLg48b///U8IIcTPP/8s2rRpozzaIOePwQzrRlaqDC6tcGI/0OCo0+Bo1N5Z42e/frurYl3usB8ccXUFwbr+h5I5tc44e65Q/jHtiysw/tYUXoO3G+dg6qBR88DBjMBbv369mDNnjpg9e7ZYu3ZtQGLwx2CG9ENfGrzQarVyv8W6jyH1FezbaKnNsC686agWhXWRUUcF0h1dyJH2U6vVbtf8sn4e6WKRo2yQUOxjuOKtdqox5/HHj2F/1+cItEANMHAwiBrDpzUz7r33XqSlpeHiiy+GSqXCkCFDAAA//vgjevToofR09P+kuYwFBQWoq6tDu3btsHv3blRWVqK6ulqe59myZUsA/8xVTUxMhNFohF6vx/Hjx+U5lVJtB6mmgzTnslu3btBoNACA0tJSnDp1So7h1KlTiutISHNis7Oz5ddhcDIndOnSpTh79iwKCwsVzal1xtlzeVoPIxg4e02NnTPrT57OCw4m1nNJvTGvNBTmzxKRd1177bV47LHHMGXKFAwePDjQ4XidSqWCXq/H8OHDodPpoNfrERkZierqakRGRiI7OxuRkZHQ6/VIT0+X+woZGRkwGo0wmUzIzc2F0WiU61BItSYA2PRhJk+ejOrqagghUF1dDQDQ6/Vo0aKFHI91TbDs7GyUlJTIfYw2bdrY1PySntu+7ldGRgbatWsHIQTatWuHI0eOwGKxoKCgwOZ5GtPHsG7P7dt2b7T1np7DW+1UY87jj1oK/q7P4Q3u1KJwts/69ethNpuxfv16X4dpo6H3ORC1WaiJ8mS05KOPPhLz5s2Tp5sIIcSSJUvE6tWrPTldUPNHZoaUHmk/TQMOrmZIVwSs118X4p8rG9JVEfslWQHH0zuk86jVapGUlCQf76hKtz1PlwRzxtF8WE8E01UTb8XSFLIdQgnTJinUMTMjsNatWyduuOEG0a1bN9G9e3dxww03BCQ7w5eZGVJbbd0+SX0KKdPTvt1yteKa9fQQ6+m19n0Y+xVP3Jmu6miFE/spJdJ2Z0vCurNiiTus3xf798gbbT37C66FYpvuThaIs31mz54tZsyYIWbPnu2z+By9p8zMoMbw69KsTZ2/pplIAwqOBjLCw8OdThexrn1hPU/VUeMo3aT0SfsBC+v9HDWy0vPY7++qwVTSqFrPh/VUMA1kCKHs9buKPdheFxEFNw5mBM4rr7wiwsLCxOjRo8VLL70kXnrpJZGeni7Cw8PFK6+84tdYfD3NRKPR2FxAsR7McHSBwrp/o1KpbB6znpKiVqsdFhuX9nFV0NNd9oMkjvo19su1eoN1e65kaVfpcXeWmG3M+xLIJTvJMVfvZ0PFU1esWCFmzpwpVqxY4bOYPJ1yw+8JOeP1wYyXXnpJnDp1Sv63q1tT48/BjIZuSUlJ8nkcNVbp6ek2+0udCanxlwYx7Ac7rM9pf0XC/jmsG3pvNLrWvJGZoeRKjT8o6VTwagoReQsHMwKnU6dODgctFixYIGJjY/0aiz9qZli3W/YrqNm3adYFOe0HDZz1fVq0aGGzyomzwQVPfsTb9zvcqecRSPYXnRriyXti/cNUyY9Nb/y4DRah9CO7offaV5+Fp98Tf8RGoc/rNTPmz5+PkydPyv92dsvJyXHndPT/XM0JVKvrfzQNrdtsPacTAKqrq/HBBx/I80vbtm2LuXPnQq/XQ6/Xy+u3W9fZkGpZZGRk1Ku7AQDJycnQaDRITk52+Djwz3zNrKwsVFZWIioqqt4+jvYHgM6dO2PAgAEuX6cr0nxZAEFRN8PZe+SIN2tNhFKNDV9xd71zIiJvqqqqwnXXXVdv+9ChQ1FVVRWAiHxLpVLJ7VZqaio0Gg1SU1MhhAAAlJWVISYmBrm5ubBYLDbHJScnA4BNWx0eHm5z/tOnT8NiseC///0vMjIysGDBAhgMBrkPI9W8mDx5MkpLS/HII4/Y1MFwJjc3FwUFBViwYAGWLl0KoH6bLbXLCxYsaLAdb2y76+z43NxcxMTEICYmBu3atYNKpYJWq3Wrr+BJfQ/rWhBK6hpY7xuMdaKU9AlCqZ6Do/fa+rX66rOwPq+ntU6C8XtCIcj3YyuhzVeZGTNmzBCZmZlOr0Q4mm5iXQ9DWu4sPDxcPvfChQvrXfWwvlkvASvEPyP8zpYcc8Q6Jmf7K82OsL7yAS9lJjhL3WxsWqq3+Wr6CLM8OOJPzRszMwJnzJgx4oUXXqi3fc6cOWL06NF+jcUfmRktWrSoN/3Uftqr1B5JmRBSlqjURknb27ZtK/dXHPVjrNt0++eQ+kXO+hPOsi681U66ez5n7b6z462zMRzVFPHkudzlaWaGLziaLqHkOZX0CUIpM8MR9n8o1LFmhhf5cppJUlKSyykl1g21s4ZdatylKR3WDbmjY6xrbEgDJlIHwFmtCuvG0P45GtrfFfsBBl8ONPhiwMQblEzbUYI1NkK/M0Lu42ddHwczAmfWrFlCp9OJ66+/XsyaNUvMmjVL3HDDDaJNmzZi1qxZfp2a64/BDKmvIV1wse97qFSqegU07aegOlpq1dFNarutf+Bb95uk8zrqT9hPGXHWTnpav8rdqbLOBi1cxSO9Lvvi774UbH9XZ86cKWbMmCFmzpwpb2tOAxRKNKfXSk2TkrZLJcT/5wG6MGnSpIZ2kc2bN8/tfUNBVVUVdDodzGYzoqOjG32+mTNnAgCEEPK/G8tgMOD48eOorKxscF+tVouTJ08iPj4epaWlUKvV0Gg0OHv2LFq2bIm5c+fWS6WU9pWmQEyePBk1NTUYNWqUnJrpCevzlpSUeHweV3Jzc2EymZCcnIyCggL5v0aj0a2pH/bnkdI6pX8rOUdD55XSQTUajVsprc2R9ful5P0pKiqSU1/dTYP05BgKjJycHJjNZuh0OmRmZgY6nKDg7baL3Ne1a1e39lOpVNi3b59PY/Hm98C6z6KkD6NWq1FXV4ekpCQUFhbW61OYTCb88ccfOHPmjM1xKpVKnqqiUqnQtm1beWnX3NxcPPLII7DuwjbUl8jNzcW4ceNgsVhctrNSfHq9HlFRUTbtvtTXcvRc7vZpHLVj9n2MrKwsAJBfr6tj7RUVFWH9+vVQqVS49tprPW6/gu3v6sqVK7Fjxw5cdNFFuPXWWwGwnW7ur5+aLkVtlzujI1dffbXNLSoqSmi1WpGYmCgSExNFq1atRHR0tLjmmmsaMQYTnIKpAKijm1TU07qwp3TTarVy1oV0pUS6QpKent5gSqbEG0u0eus4pcd4K43U+jzezKawr2pu/zkww8KWp5+nJymXTNMMHbwKVR8zM0gI/xUAdfcmhJD7H9K0WdhNEZFuer3eZuqsdXaCfSaHO32U9PR0p8u/Otrfeul6+8Km3l51zFEfw9O+i9R2uWq/3Pmb6cnf1ebwtziYpt6wn0JNldcLgH799dfy7aabbsLVV1+NAwcOoLi4GMXFxSgvL8c111yDG264wZ3TkRcJIVBdXY3IyEhoNBqbx2pra9G5c2cAQFxcHDQaDYQQqKysxPLly+X9VSoV0tLSbIpPWheiaqjQp1R0S7qSID0WExODVq1ayYW/rLfHxMQAgHxVxvpxVwW0lBaz8lZRTevzWP/bk+Ja1qyPty5sZn0lKJCFTB19HoEsLurp5+lJkanGFqZiAVL/8bT4GJE/WCwWbNmyBUePHg10KH5jXcRcpVLJ21q1aoUPP/wQdXV1qK6uRmVlJbRaLfr27SvvJ6mpqbEpGAoAlZWV8nHS/lqt1mWhbakdXb58uc1xzkh9nsjISABAXV0dgL8zP7Kzs22ey75P427Bb3v2fQypULsnfZeUlBS5iGpcXJzDfdwpcOnJ31VH521qbaGnRVF9gQU0ieBeZoa1Tp06iZ9//rne9u3bt/t9yTF/CPbMDOkmXXFw9rg0R9TRPtK8Ves5rO5cAZf2ka6oWB9r/zxSkVH7pWGtr3ZYn9PR1Q+lS726qzEZIt5Y594XWS/eYv09cFa4ranjkmMUipiZETgTJkwQb775phBCiLNnz4rk5GShUqlEq1atxNdff+3XWAKRmaHVakWXLl3qZZE621/KGHX0WHh4uNBoNCIpKUmuxyGdq0uXLjaZjQ1lTDirz+WqVoWU4WpdB8N6f+ssisYu3yq1NStWrLBpc9xtg6wLZDbU/vgqY8DReb3RFjoq/hko3sjM8Mb739A5mkOWDDVdXs/MsFZVVYU///yz3vbDhw/j+PHjSk9HXvLBBx84rJnRpUsXaDQaJCYmYvLkyTh27JjN43q9HtnZ2TCZTKisrERlZSWysrLkKwNHjhypl1kRHx+Pfv36oaysTL6iIi3FJl0BAf6+WgIALVq0QHV1NSwWC2pqalxecXC1tKoUY0NLvSrlSfaDdExBQYHbV2IcZTQ0tHRrQ48rzZJQur/11aLly5fLn6E3sl1ChadXVnjFhKh5WrFiBS655BIAwP/+9z+UlJTg119/RWZmpk0GY1NVXV2NsrIym23CRXk28f8Zo46Eh4ejc+fO+OWXX1BdXQ3xd+F6AEB5ebncPmZlZdlkiDrKLl26dKnNfwEgJiYGjzzyiNwHsM60qK6uxiOPPIK+ffuioKBAjsm6zyD1ldRqNRITE7F37163MxDsMxaktmbHjh02bY67bdCOHTsghMCOHTsabH+cZV00JovCWf0Gb7SF0mv7+eefA57lYf/euXrPnL3P3sjYaOgcobS8LFFjKB7MuOWWW3DvvfdixYoVOHDgAA4cOIAVK1Zg7NixGDlypC9ipEZQqVRYsGABNm/ejOrqajldUiINDBiNRpu00IyMDERFRcnpnNIPfakRLyoqghACtbW1OHLkCCwWi1xgUxrYaNeuHQDg7Nmz8nlHjRqFiooKVFRUICMjA9nZ2XLqpvS8JSUl8nbrH8yNmTLi6kd8Q+d1dKwnsfhiyoj1Od0ZqFAag/VgSlpaGjQaDUaNGuVxKq2v+WIKjKcdMU59aBr8mSLd1NKxm6u//voLHTt2BAB8/vnnuP3223H++edj7Nix2L59e4CjCw1ScXLg7wsb1dXVAGAzRaRly5b1jjt27JhcJLO0tBSPPPKIfEHGun0YM2YMHn74YVRWVsqDI8nJyTYXdg4ePIjWrVvj6quvhtFolP//zMjIwGOPPYbx48cjIyMDeXl5mDFjBlJTU9G6dWt89dVXbr1G+x+bUltz0UUX2bQ57rZBF110EVQqFS666CKnP7hXrlzp8m+Mq2kinhwLOG4Llf6tk15beHh40P1A/+qrr2A2m93+3AHvDPA0dA5eUKHmwq3VTKxVV1fjsccew9tvvy1Xng4LC8PYsWMxZ84ctGrVyieBBkoorGbiSlJSEjZt2lTvqohWq0W7du1sKmKPGTMGH374obyqCfBPRe3U1FTk5eWhpqYGkZGR6NatGzZv3oy0tDQMGDDA5mqTVOlbqilx5MgRVFdXQ6/Xo6KiwuMVKRqjMSunODtW6eto7OtuqAK61Hlz9RoD8d77kz9WyCHnmmJldX9W9Pfmc3E1k8AxGAxYtGgRBg0ahK5du+K1117DjTfeKF8x92ftjGBYzcRbVCoVWrZsierqaqjVarz66qs2baG0UonU/5DuA7DJ+DQYDDhw4ID8mLRqisFgQHJyMpYtW4bIyEjMmjULarVa/nsm/f8p7S/9fyptl0RGRmLq1KkNvh5f/r20P7ez2N2JqTHHOuPp37pgbGNmz54t943d+dy9KRjfDyJvUNJ2Kc7M0Gq1eO2111BRUYHNmzejuLgYlZWVeO2115rcQEZTIGVQSKSrGmfOnMGBAwewceNG+bGCggK5oGhWVpY8LSU1NRUFBQVy4a2oqCgUFhbi7Nmz8tKsZrNZThPV6/XylKOSkhIMHz5cztYAGp+h0NDVd29lUjR0bGOyHDzh6Pmsz+nOa2xsDMHOWwVfyTNNMa3Vn1e3eCWtabj33nuRlpaGiy++GCqVCkOGDAEA/Pjjj+jRo0eAowtd1n0ZjUaDcePGYcyYMQAgZw+qVCocOXIEALBgwQJ5WmtycjKOHDkCtVqN5ORkOdMwPT0dr732mtxvyfu/9u49Oqrq7B/4d2YyyUxCGAjXRCNgUKAJXkhaRYNIYkGgLBDaUrAUFfWHYhGpQoW2eKsIXoquFpS2C7RLBd9Ge0uA8CYFTNG3BENLIgIJEEYTBRMNEDK5zf79ke7Tc86cuSWTzAz5ftZiwUzOnLPnTJi9z3Oe/exdu+B2uzFo0CAsX75ck1HgL3MiIyMDDocDubm5hu3XZyPIjAUAIc/IKioqQkNDA4qKiny2Xc8oi8LXa9XvKZhsRG/fdf4yNsKR8eivTbm5uT4/92D2GWjGitxOZoWEu89lViGFU9CZGb1NtGRmyDXYvc07lbKysnD27Fk4nU643W6YTCZcccUV+OlPf4p9+/YpmRk2m03Zl1yPXb3uOQCPjAC5nTpDQH13RN4t72p2gNEa9ep99dTd+WDeRygyInr6eL0N73B0Hc9h5GBmRnj98Y9/hNPpxPe+9z1lVbHXX38d/fr1w8yZM3usHdGUmSGnlehXMFFLSkpCYmKi4Rhm9erVyrglPj5euVv+4osvaup5GWVZynFKUlKS8vwvf/lLpf/09d1WWlqK4uJiCCGQm5vrMZ1Cvk4Ge202G+Li4jyyJgLNUgjke3b9+vVoamoCAEyfPr3bvo9DnbnWk5lwgeqONnnbZ6DHktvpf5fCJRI/N4pu3ZqZQZHJZDKhvr5eWY4L6Ahc6JdrPXv2LE6dOoVx48YB6CiuJe/279q1C0IIuFwuTJkyRamhcdlll2Hx4sWaWhdyYLBkyRKcPXsWSUlJuP766/HQQw9h0KBBmkBDe3s7LBaLcrfcV3ZAIDUP/C2N6u/u/Pz58xETE6Pc0emsYLIc9O3sTG2HrhyvM8K5BGs4XIpZBT2NdUKIOnz3u9/FI488ogQyAGDhwoU9GsiINu3t7R6BDFlIXE3W5gI6amvol4Y3m81wuVzK8q8/+clPPJY7VfdvcpxiMpnw1VdfweVyKbXC5B3nwsJCNDQ0oLCwEOvWrcP69euVu9DFxcVoamqCy+VCUVGRYVFPGXyw2WxwuVyGtTJ8ZWSp73zrsy6M5OTkKJm4wfZpwdxl17e9q3foIzE7rTva5G2fgR5LbpebmxsRfW4kfm7UezAzw4/uzMzIz8/v1pQsk8mEmJgYpbYJAMybNw/79+9X6liYzWa43W5kZWXhxIkTmrsasuCWxWLRFPEEtHcyAHjMQZWrpAAIKkMg2KyKzmQgxMTEKAEW/fvq7D6DbWd3Z4+E4j2Euo2hOq/dlXXCrAK6lDAzo2e98soruP/++2Gz2fDKK6/43Hbp0qU91KroyswwIus0yAtz+W+jOmAzZ87E22+/DZPJhH79+im1Scxms0eQZMCAAaivr4fZbMbcuXOxf/9+w2wPGXgwYrfbldXa1OMsAModal+1K6ZNm+a1r8nLy0N5eTmsVitGjRqlrObhcDjQ0tKCpqYm2O12rFixwuu5C7RP89ZGeXFaVFQEk8mEnJyckNTCCGVfK/eVmpoKp9PJ/pvoEsHMjChgMpm6PYIphPDoYLdv366pDC5XN/noo4+UdgGAy+VSXvP973/fY9+LFy/WzEf96U9/qtlO3tXQZxP4utv/6quv4vz5816XbdVvO3z4cADwu3/1Umuvvvoqrr/+egBQ/tbrjlVH9OehO2s7hOpiP9RtDNV57Y7PB2BWARF13q9+9Ss0NjYq//b2Z8OGDeFtaJSRQQv1kqzqDFTp4sWL2L59u7Ktusjq3Llzve7f7XYrS6zPnTsXZrMZMTExSt0vfT0M9bGFEGhoaIDFYtGsrCLHd0YX7XJ/vgIZAFBeXg6go75ZeXm58t6zs7ORk5MDh8OBnJwcH2cu8D7N26oqcmqMy+VCU1NTQBkeRnfovS1BG4osSG/L2RJR78FgRpgIIbr1S1e9fJmafmlWtfr6euV1NpsNQEdnKIt8SjJgsG/fPiQmJioZGLt27YLVaoXJZFJWS5HTOeRr5FrwRheickk0uVysL94uaI2eVy+1tnr1apSVlQGAUiRMT30RH6qpFvr9dGchzu662O+qUAVHuisQ1FsLWAX7vnvreSLy5eTJkxgwYIDyb29/Tpw4EeaWRgeTyQSLxaJZMl5qaWkxfI16fGM2m2E2mzFv3jzDMYzL5VLGKrKvl1NtW1tblRsyWVlZyM7OVu76T548GTabDVarFe3t7bDZbMjNzUVsbKzS7mnTpgHoWJZXf4HtbclU/fdpTEyM5r0AHeOyzgbb8/Ly8NRTTyEvL8/jZzIAkZqaqgTbZBvl1Bi73R7QDTijAIqvYElnlJaWYv369Vi3bh1SU1MDKmwaSULdh7JPpt6OwYww6s4vHhnFNxoIGAU65CBA/i0zNz766CNNVgPw34vld955R5mr+uCDD6K+vh6tra0QQsBms6GsrAzt7e14++23lSAGAI8LUXmhf9NNNwV8kaq/oFXvIykpCWfPnlXarJ4rC8CjhoeeOtAQSGAgkIBHTwYYQnWxH+o2hyqA012BoN5aMyPY9y23Lygo4OCJiLqFEALt7e0e00m80Y9r3G433G433nnnHcyfP1/TRz/33HO4ePEihBAYNGgQAChjGHk8k8mk9KHq78isrCzExcWhtbUVra2tiIuLA9ARYImJiUFMTAyKiopQWFio7KulpcXrd6W3798pU6Yo7ykuLk6zWkZn+io5TaW8vBwbNmxAXl6exyokTqfTMPiycuVKrFixotOBFH3woqtZkCUlJUqNEqfTiWXLlmHOnDlRk1kZ6rFGbx27EEmsmeFHd9bMkH+Xl5fj3XffBQAlEn769GlcccUVKCkpQWlpKWbPno2MjAzU1NSgT58+mgj36dOnMXLkSAAdhahk5yR/tmfPHtTX1yvPqV8r95+VlaWkLBYXF+PMmTP49NNPNamdSUlJqKurU6Yx3HTTTdi/fz/Onz/vsYpKUlISpkyZgrfffhtAR1ClX79+hnU05NxVmc3RmekR6voOAHxWLQ9mCkYg2wdSWyIaVxeJxjZ3RW+tmRHs+y4tLUVBQYEyf5uVyyMTa2b0rOXLlwe87UsvvdSNLdGK9JoZWVlZmDx5MqxWK1pbW1FYWAgAynikpKQEn3zyCVwuF77xjW8gJycHsbGxMJvNqKmpQVJSEkwmE44fP46rrroKQgiUlZXh4MGD+OqrrzTBEPWKaurV2Z577jkMGjRIMyaS4x3A8ztSFuGUN0ZMJhOampoMa3lIFosF7e3tiImJweDBg1FTU6Nkj1itVsN6FN6+m9XPy3Pka5usrCzk5eWhoqJCqRUm26r+Dg+mLzA6Jz3Vf/paPaY7jxnqOh+hOle9dexCl7Zg+i4GM/zormCGmv4jUBe8Uv/M6Hl9cayWlhY8++yzWLZsGfr16+d1e7fbDYvFoqxeEhsbq1QHd7lcmDNnDrZv3449e/YodxSSkpKwa9cujw5sx44dqKmpweDBg1FeXo6dO3fixRdfRFZWFnbs2IG//vWvOHDggDKQGD58OAYNGoRbb70Vc+fOxZQpU1BfX49vfvObuPnmm/Hxxx9j165dQZ1X9YU3AM1AJdgL8WALSvW2i34iDp4iH4MZPWvSpEmaxwcPHkR7eztGjRoFADh27BgsFgsyMzNRXFzcY+2K9GDGqlWrlCkaQEdWaEtLC/r16+exrXqso3+s/rccw7z99tvYt28fSktLlaXjs7KylIKWaWlpqKqqQlNTE9xuN8xmM77++mts3rwZTz/9NMxms8/vOFnsEuiYApKUlISamhqPop3+dLZYpq/AsrptdrtdCZaEqmCmvtCnXArWX1HSSOftfOvfL/tAou7FYEYI9UQwIxRkR+52u1FQUIDbb79duWNgtJ1+UCBXNZHkNi6XCxaLBQcOHMAdd9yBhoYGpTjo9OnTlSXJ1K+bNm0aSkpK0NzcDJfLBbPZjN///vdYunSpcpza2lr06dMHDodDeRwfHw+z2QyTyYRf/OIXft+ztzsD+k462E5HPQgAjAcavb0j6+3vnyjSMZgRPi+99BL27NmD119/Hf379wcAfPXVV7j77rsxYcIE/OQnP+mxtkR6MGPNmjWasYhcTcRb3S/92MUb9Y2bPXv2YN68eVi8eDGefvppZaxjdMNIXe9CHSQwGm/oxz96cpqvfE/ehtsy+GEUXPC1Ooi3gAWgDXQAvgMmnaHPDsnPzwfQEdRZuXKl1+27EkTp7OoswfB2vn2t+MLsRKLQC6bvivH5U4oasnM3m834zne+43c7o/mlanLA0K9fP7hcLkydOlUZCEglJSUenXN6eroyf89utyvzJH/+859rvvznzp2rdAwlJSXo06ePsg85B9Uf/TxW+fjcuXNKgVX18/KxL6WlpWhubtacL6OCUsHs81IUyPtnwCNy9XRKMH8PqDd58cUXUVhYqAQyAKB///545plnMHny5B4NZnSnNWvWKP/WZ38KIWA2m/HZZ5/h0KFDyg0Wt9uN5uZmNDc3G2ZfGNX5UgskkCGnrMj95eTkwGKxYP369ZqxjsxOldM/pkyZAuC/F+YAkJqaCgAoKiqCy+VCfn4+qqurldVGfJFBDABKYES91CkANDU1obW1Vcni0PepcoxkNA6R+5KrjRQXFyuvlX/LY4W6MGZWVpZyDFk01GQyKXU99PTjs+LiYp/9glG/IYNH6vfp61idGZ95O9/q9+trOyLqeQxmkAez2QwhBJKTk5UVP86fP49nn31WMxBITU1FZWWl8tjhcGDYsGE4evQoAKCtrU2T2tjS0gKbzaZ0TrJjqK6uRkNDA8xmM+Li4jTLjfmaG6rvTORjdeRf/by/Tkd9J0OuH++to+3Jjiycc1O9CeT9+xtQRML76K16MhjX2wN/1PucO3cOX3zxBdLT0zXPnzlzRlnyM9rpgwrqx+p/X3bZZbjsssuU5ywWC+Lj42G32wMKTHSGfkl6AKipqfFob1tbm/LYarWipKTEY6WU8vJyHDlyBO3t7Zrn1GJiYjT7UhNCwGKxwGq1eox91LUfRo4cqRm3qPtHb3f+5c0amS2rv7mkvwAPVLBZFOrxgLft9OOz5uZmn/2CUb+hrjUHQKkDkp6ejjlz5hi2J1iBnrPOnlsiCj0GM3ox/dQSST6n7vxlZXApJSVFMx/UarWiubkZxcXFykCitbUVRUVFKCkpwcWLF9Ha2gqr1YqioiIUFxcrgY6qqirlGOola/VZFQA0/w70Aj/QTkdmmphMJo9iXEade0+lFnrLQAnXxWGgQQh/A4pwv49o0F0Bn54MxoXqWAx+UbS44447cPfdd+PFF1/EjTfeCAD48MMP8dhjj2H27NlB7Wvfvn14/vnncfDgQdTW1uK9997DrFmzuqHV3cNbwKK7AhmBkBkiaq2trWhqatIsiSqpAxlGTCYTbDYbXC4XrFYr4uPjlTEP0DG2WblypbKEpvrmjNFYY8OGDWhpaUFTUxMKCgoAQBMAUX8PpqamoqGhATExMcjNzQ3J96S3LFdvAhlj6bfR36jSM+o3cnNzNc/JMWhFRYUmmNFTgQb2SUSRgcGMXswokBEIh8PhcZejra0Nra2tsNlsmudNJpOm/oRczgyAkmaovqPR3t6OhoYGFBUVAYAmmwPoSJlsaWlRUisDvcAPpNPRd57r1q1TAhvBdu7eGLXDX9u8ZaAY3cHpiQ410CCEvwEF0zT9f3bdFfAJ5WDPX2X5UB2LwS+KFq+++ioeffRR/PCHP1T6u5iYGCxatAjPP/98UPtqbGzEtddei7vvvltzwUadZxScaGtrQ0xMjEeQRX3TR65IohcTE4O0tDRUVFRg1KhRmDNnjrJyGwAkJycD0H6HyWBFcXExqqurUVFRgeTkZNTW1moyLPRTMuQ+ZJBDZsfKcVQovidlgESutNIdfbS/fsHo5/rn0tPTlcyMrurMSi7yM2SfRBRevicnRqCNGzdixIgRsNlsyMzMxPvvv+9z+7179yIzMxM2mw1XXnmlss44dZ46OKFms9k85kumpaXBZrPBarUqf0tNTU1oaGiAxWKB3W7XBEJMJhNKSjrWEo+Li1M6sbi4OGUJNPW65fp1zIGOFMQnn3wSTz/9NPLz85UOvrS0FOvXr8e6des0a79nZf137XOZuik7Krn/9PR0j+MEQ59pon6uoKDAcC16dbvkY/USckb77Cx5V8ioHZLRue4M/fvqbHsC2SYSGLXT32cXqnPdGYGeV/n/VAYXu0s4zwVRMOLj47Fx40bU1dWhrKwMH330Eerr67Fx40YkJCQEta+pU6fimWeeCTijo7m5GefOndP8ocDImzKSnHIr6QMZ6hXgysvLIYSA0+kEAM1Yp6amBk8++aQydkpISFAKh8rMAiEEampqvK5uJ/uJ7Oxs5bmSkhJN8CU/Px/nz5+H3W4P6ntS/10v30N7e7sy/uos9b599Smd6cfnzJmDadOmwel0Bvw6b8cJZhwlt1XXQQlWtIxbiKJBVAUztm/fjmXLlmH16tUoKyvDhAkTMHXqVJw+fdpw+5MnT2LatGmYMGECysrKsGrVKixduhR5eXk93PJLnyyipe8InE4n4uLi0Nrairi4OEyePNljbm1ubq6ycoqcitLW1qbMLb548aLyhS8vaHJycjwu8NWPS0tLlXmt6gwUGQTwd/GVnZ2tpJumpqYq+58zZ47fC3BfjC7I9IOTQKg73lBe5AXSoQcShAgVX+2Rg4GioiI0NDSguLg4ogcHRu/F32fXk+daL9DBXXZ2thKMDHWgQT3gC+e5IOqMhIQEXHPNNbj22muDDmJ01tq1a+FwOJQ/soAlBc/tditZFXrqLA11AKKlpQWlpaXIycnRBDTU1Jmtubm5SE9PV27Q6CUnJyvTWOT3q81mg81mQ2pqKoQQmuO43W7ExsZqxkL++kX9d312djZsNhvsdjtSU1O71K+q9y3/XVRUFHRg34iscxbM67wdJ5hxlNw2Nze3031SKG9CEfV2URXMeOmll7Bo0SLce++9GDNmDDZs2IDU1FRs2rTJcPtXX30VV1xxBTZs2IAxY8bg3nvvxT333IMXXnihh1veO7S2tnpkbWRnZ2s6iaysLEybNk3J0rBYLCguLtbc8ZB3R2QQorW1Ffn5+Vi3bh0AeAQtvEXZ9Ww2m5LV4O/iKysrSxl8yrsUgfI1eDC6IJPnJJiAhP6cyvodXb2YD1VgJFR3HXy1Rw4G5CBQffcqEhm9l0jOTgn0dyErKwsrVqzAypUrQxJoUL8XDviIgvP444+joaFB+RNs/0VaNTU1hkEJ9ZhF3qCxWq1oampSaoXJzA1fiouLAQB9+/bVFB+Vx6ytrdVkqKozVp1OJ1wuF+Lj4zVtVE9BDeRiX/9dn5WVhZUrV2LFihVwOp2dCjLI73D1vuW/5fTjYAL7RtR1zgIp8C6vGbo6xglFYJ2ZhkShYxLeFr6OMC0tLYiPj8f//M//4I477lCef/jhh3Ho0CHs3bvX4zW33HILrr/+erz88svKc++99x6+//3v4+LFi4YdlFwyTDp37pwyf7Cra7QDCMm67NFETj0x+tJft26dUoXbbrcr2R1yWTdvRbcyMjLgdDqRkJCg3OHwtia4nNMo148PtE6FfhsAfreXr9GvTx+sztbAkEvf6t9rOPTEGuyRuMpLdwjkXAayTajOT3efZ/V7CaRKPhkLZo12ig4mkynoAqCh/D3obeOXQNhsNrS1tSk1K+x2O4COabTqJWJ9kZmZcnur1YpBgwahtrYWycnJqK+vV+p3AdCsgFJZWQmTyYS0tDQ4nU5llTm5Wos8/vTp0wEENo5RKy0tVZZ31Rcr9cZffxSOvshfm3pizEJEgQum74qazIwvv/wS7e3tGDJkiOb5IUOG4PPPPzd8zeeff264fVtbG7788kvD1zBFM7TUUzlkrYpf/vKXePbZZ5VABtARRHK5XMqgwFf18PLycjQ0NGhSNfXRbRk5T0tLg8lkQnp6uqazC3Y6RXFxMRoaGlBYWOhz/mcwdwq86ezd6M5MV+kuPXHXwaiWSLRNQwhVjZJAtglVlkN3Z0sYZR1F02dKRL2HfonXpqYmpR6GPpChrh+mlpycDJvNplkJThYCbWxsVGqFyaKfsbGxcLlcqKiogMlkQlNTEyoqKpCdna1kazQ1NWmKz5aUlGDXrl2aAuuBUNcqC3Qqp/wO9zZFRf293pmsQvkaAAH3D/76SGZKEEWvqAlmSPpK0/LCMZjtjZ6XmKIZeufPn8eTTz6JgoICNDU1eRTZAv5b18JXopC3+admsxnV1dWGHaLT6dQU5ZL8dbaA9kJTtktOpVHP/9SnStrtdsTFxfk4I77JfTQ3NwfVwXdmukp36epFaKinTURqIVFfc4glfcFXI4Gc71AN1rp70McABhFFGm9TRtLS0mA2BzaUlnXBXC6X5jUyYKGmrpOhvlFRUFCA1NRU5bEcA8tVT/RTVRwOB6xWKxoaGpSgi7+lZvXkd77RVE6jflN+hxtNUdFv35ngeGde469fYb9DFL2iJpgxcOBAWCwWjyyMM2fOeGRfSEOHDjXcPiYmBgMGDDB8TVxcHPr27av5Q11jFKgwWv7MF5vNhlWrVmH69OlKYU4Z3HC73aioqNB0bv7mR8qOq7KyEg0NDcjPz/e4UJTZGPn5+Rg5ciQcDgcyMjI85n/qayDIuyby4jPQC2T13Qb1PoLRmQ450DaG6mI/kOrmob77H8j+wlGfwdcc4kDbFujnYvS70ZnPlIM+op5z4cIFHDp0CIcOHQLQUdj80KFDXgufU/fwNkYpLy/v1DL3brcbVqsVJpMJqampHlnAbrdbs5KbXH5UCIGqqirlxoW6gGhrayuampqU4p1AR9ZrWlqa8hjwHpjxRn7n5+bmeox5fPVNRmMko2Kj/m4sSaGse0FEl46YcDcgULGxscjMzMTu3bs1NTN2796NmTNnGr5m/Pjx+Otf/6p5rrCwEFlZWV7v8lPP0GdgJCYmKsVDHQ6HMu1EMplMyMvLg9PpxJQpU5T0RPXcUafTiezsbE3dinPnznnUj1DPs1TfodCvFa5uo9PpNJxHaXRBJ++ip6amKu0IZB1y/Qol6nodoZCXl6esyT5nzhzD4/pqYyjWr9fvB4DhPkP9/n3tT/4+yMFkTw6Q5EBVX59Fz1f7u/K5hOozJaLuUVpaikmTJimPly9fDgBYuHAhtm7dGqZW9T6tra0wm83KlIuuslqtcLvdEELg6NGjXgMMTz31FNLT0zXZpUIIpe8A/ju9ta2tDVarFe3t7ZopK06nE7GxsUr9sNzc3JDUrSgtLUVLSwtiYmKUVVz0xc31+9b3ZXIbWbPCV18k+ysArGtBRIqoCWYAHZ34ggULkJWVhfHjx2Pz5s04ffo0Fi9eDKBjishnn32GN954AwCwePFi/PrXv8by5ctx33334YMPPsDvf/97vP322+F8G4SO4ITJZFLuaMjClUIIpKamoqqqSrN9U1OTstSqnO8pi1Kpi2BVV1cr2wHQ1I9QBwkaGhpQUFCgZHkA8LgzkpubqxwjmAtcdeccTP0MdRCkO4oeyrXsKyoqNMGMQAMHoQow6PdjtE+jQVBX+NpfJAyQ/L1f9c/1g9BgPpeuvJaIet6tt97qc/ol9Ry32420tDTNGKOz1FNtvRUJlf1SeXk57Ha7MkZKSkrCb3/7W9TU1CAlJQXZ2dnKjRP9vmJiYjR9rRxfyOLoMngQSHBDH/yWK6vIuiHeAhH6ffu6CeSvNhT7KyLSi6pgxty5c1FXV4ennnoKtbW1yMjIQEFBAYYNGwagYwkrderliBEjUFBQgEceeQS/+c1vkJKSgldeeUVzIUfhIed66p8DoAQmYmJiNIW1JJfLhfz8fOWxvEgHYDjIaG5uRlFRkTJtQ93x67MvAul0A9XZlRgqKyuVtobqgr60tFQZCCUnJ2t+Fuj7DFWAQb+fcGcEdNcAqbsqtquDcUBwn4t+MBrqoBER0aUs0ECGxWIJuDZFSkqKpqA50DGtRWasWiwWTTaIetuamhrcd999KCsr89gHACQkJGi+4+XYx2azaaZqBJKlp+8r1TdgZGaskUD2HUhfxP6KiIxEVTADAB588EE8+OCDhj8zSrmcOHEiPvroo25uFXWFvMiWWlpa4Ha7YbPZ4Ha7/c5HTU5ORk1NDcxms2ZbucyZy+WC3W5X5mWWlJQoaZsXL15Utj9//rwSJOlsIMFbMCSQ5c1khy/bql4rvqsXxSUlJcq5qa2t9UgH7c26MkDy9dmEagpHcXGxUkleZlMEM31JLRLubHVmqT8iomgSaCAjIyPDIxMV6MgCOXfuHADtlBQAmrFOSkoKSktLNYEMu92uLNWqr1ehLoKv7rcC6RuMbkQEupJIuPsdIrp0RV0wg6KXfmqJZLPZNHcd5M/VNTN8qa+vB6CtcZGSkoKvvvpKk3KpnmICdEwpWL9+vbKNul2hWFZVvxSsfD/qn6kvhr1lchjtM9gAR3Z2tpKd0pmL4K4IVYZCT+tM2q1adnY2iouLlVVpOvve5e+1/Fs9TzrY39Nw3tmS51NdD4e1OoioN6uqqtKMXWSgQmamyhoX6im0cXFxiI2NVTIiiouLlddnZGQo2ccycFxYWAiLxYKRI0cCgDLFRJ/dp14qNZT9NTMqiKg7Rc1qJhT9hBCGWRZdLabV3t4Om82mWSGlpqZGE8iQHbe6CrZcclVfDNZutyud+vr167Fu3bqAV3vIzs6GzWZTimHJgUFCQgKAjjss3ip7e1uCM5CK4JK31SmysrKwcuVKTJ8+vcergIdjlZBQCKTdvpYp1a9sYySQ1URkBfnc3FzNvtUrioRjWdlgyfMplxy02+28U0dEUcNoRRO5GlVKSkrA+1FvKwuYS263Gw6HA1OmTIHD4VAKmMsCoDKjTb30qRACDocD06dP9yjurV4OVhYyz8nJUTJi9X1TsGMLIqJwY2YGRT2LxYK4uDiPTA6r1YqEhASkpqaivLwcQgiUl5cjIyMDwH/njppMJiXVs62tDUII5OXlaebG+iuSpU6dN5lMSmEtoKOIl0wXjY+P97lqh9Gd/kAqgqvb6Wtqg9G+ujvtP1pTTNXzgb3dqfJ3x8nfe9fXv5DPqY/l7xjq1XsiOdOhszVkiIgigdHNGLlimv7iX2ZIrF+/XnPDJiMjQxOYyM3N1bxWPf3D21hBPp+amopz585h5MiRhrXgsrOzNfXF1CuIAMbZfZ0dWxARhQuDGRS1rFYrhBBwuVxISkpSpo9I6gvzqqoqZUBRUVGhSe2UK3zIAUlDQ4NHka+WlhZlaVOji0b1NBKjwlreCmTpUzuNlgfVB1B8TX8INHCg3oe3KTCh0tUU03BNUwl0yThf7fOXuquvfwEYL1Xri5wHHeyqO8Hq6ufAVGMiikYOh8NjfCHJ7+7U1FRlm4yMDAwbNkyZxmq1WmGxWJCbm+vRhwMdRcplhqi3JVqNvj+dTieEEJplW/WvATwD5Pr9hWq1ESKicGAwgyKSviioEfU0EqMq3mVlZSgoKEB6erpmOTWTyaS8Vs5LFUKgsLBQs1Sr/LnVatUsDQt4LuMqa1IYZTd4u4DTBxRkBodMK5WMlkPzdsEbyAWj/k6+uu2ROFAJRc0Qb3ztR/7MKMDkr32BbmN0h0z+O9D32FMZD529MxetNVOI6NJjtHII4HvM4S2QIV8nv4OlyspKzSprcsxhtIrUhg0b4HK54HA4lGMVFRUF/d3vTSBjglCtNkJEFA4MZlBE8hfICIQcsJSXlysDBQDKCik2mw25ubnKahGtra0ea7S3t7dj9OjRmoEJ0DFYUd9p13f06gthmZHha7ULeVfHKOvD23JoXSlSqr6THwmDFJn1kp6e7pEua/R+Q5Xy6ms/+mKxRkpLS9Hc3Oy3/oO3z0z9eyJ/Lo8VTEaIt/aFUjC/d0aBOqYnE1G4mEwm5UZBaWkpdu3apRTYlEuW19fXK7W3mpqaYLVa0dbWhuTkZDQ2NnpkWOoDDoWFhWhtbUV7e7syXrDZbBg5cqTXpUuNppq2tLQE9J0Zqr6bWRdEFM0YzIhSdru9y4UzLwUmkwkxMTEeQQg1WXRTPXcU6FgtpbCwEKNGjfIIVkgyhVNOQUlNTUVlZSVcLpdS2NHojr48lrwY9rbahfqCT74fo6koch/V1dVegyP+yIBBcnKycvxAskZ64gJUnv+KigqPYIbRgE3OFdZnyATL1yAukAGenKLjcDg6NehUZ+ToA1mB1tuIxCCBum3hGCgzG4SI1EuY9u3bV1OHSB2sBoDGxkasXLkSQODfH/osO7lPi8WC+Pj4gL5/9H2D0TSU7hYJNzSIiDqLwYwIZrVaMXnyZABAcXGxUvVaXsyWlZUZpkuqmc1mCCGQnp6OYcOGKXcO1BwOh3IxFWqBTBfxRk4B8cVmsynvx+hYJpNJOYc2m82jSGhrayucTifS09M96mTI16sHJCUlJcqdG/kz/cBHnW5qtVqVQY2efgDha/Ckv+gtLi72ub3RYEwGDGpra/GLX/wCeXl5yjQcfQAh1BfK/gaH6enpSmZGIPzNFQ6Ur0Gc+mfeMke6eqGuLjKqv3PX1eKioRbM74R+6ktPD5QjOdBDRIGTBTPVmQtyeWeZbSGflzcb1NM91X2P0VLonfnu9SZUU/4YXCAiCpxJdMcV7CXk3LlzSvGnvn37dnl/Tz75pNefyWkPwXSG6lU3jIIS06dPN9yP0cXZs88+6zPDQR7DaP6oUaaIxWJBe3u73/cggzbqQItcarW5uVlTQVwW1pKZD3I79bFNJhPS09M1GQwyZd9msyEuLg4JCQlKICglJQW1tbXKeZPFuuSAqLq6GhUVFUoGSExMDNrb25VzJ/ftcDiwbNkyzeogaWlpnc6kUNNPW5GDOXlMPX2bAM/P/KmnnlLes/73xNeqLcEGUYza09U756G48x7MPtT/b9esWdOp40W7UGY7dHfmBDMzQt93UXQK5e+Br/GLL+pxQ0pKChobGzXBCdkfq2/gdKaf8ceoXyQiosgTTN/FzIwIYDabERsbq1S6DqZznjNnjsdd9UA6eqPXjRo1yjA7QZIDAHUARWZDpKWl4ZNPPkFbWxvMZjPi4uKUAINcrlS/rJlMAY2JifFI+1RnUMhjOBwOJXggnx85ciQqKytht9vRv39/1NbW+q27oC681dDQoAlkyICS+o68fK+tra1wOBxoaWlBW1sbqqqqDFcg0Rf3kktvGk0RCSY4YDS3Vv5bvx+ju/b6z1ydjVJcXOy10rmavzve3n5u1PbO3jkP1YVqMG2QWUL6ArHh1NMX7KG8WxiKzIlAVpEhop4jC2YLITT9qC++MuN8vaarmRNERHRpiJyReS8l1yIPpc529P5S9uUAYM6cORg2bJhSqKqpqQlOpxNTpkzRzJFvamoyTAM1KqAl9y9X2ZCMimfpa03IDIXGxkZl6oHRUmPe6i7I4l5GF0UVFRXKv+VnJffd3NyMhoYGtLS0IDY21us5k+/JKHBQVFQEl8uFoqIiwwwImYEhp67ISudxcXEe01/UF4eB/A7MmTNHqf/hL0Er0JU9vA0W9e3pyqAyVFMIgmmD+nc7UkTzVIpQXFRE8/snupTob8hEKgY5iYguPQxmhFFMTEzIAxldoS6SaTabkZiYiIsXL6K1tRV2u92jvoOctuFwOAznxuvvmhoNIvTZAAC8LnGq3s7oOOp/B3KhI+suNDY2ek05VddykJ+VPL4+qFFQUOD1PZWUlKCxsVFZBlaS9TdcLhdKS0s9ghMANCtlGFU6Ly0tRUtLC2w2W9ArTainNfmibo+v9NxAB4vq96l+HIhQ3V3z11ajgFgkUc/5Vq+sEw1CcT55l5Uo/LrjhgwREVGgWDPDj+6smRGJc+/1F3CdrZ0QboG0L1TvobS0VMm+kDU5jPZpdDz1a9XzeH09rw/2yKks6qXnfL1fmU0TzLzh7qiX0Jl2hLqdvrbvyfnVXTm/nAdORlgzg4DurZnBQAYREXUH1sygTtPfMVU/jvQ71WqBtC9U70GdZRDo+vBGr9VXU1c/r592ol4KVD2VxddxjbJpgnmPwWQx+NKVdvgTzNQDdcDIaPuevPPflSkTkZahEOmBTiIKDQYyiIgo3BjMoIBxjrp3+qknRheW3s6ft0CBUSFRu93uEQDwFhDRC9WycUaC+d3oznYEc2FfUlICIYSyxK5eTwbsuhKQiLTAIr8niIiIiKgnMJhBAYu0O8CRyNeFZVfOn78AQE9mohiRxVRlgdDuboe3u//B7Ls7gyrBirSARFfwe4KIiIiIegKDGRSwS+mCKxy6cv4i/dzLYqr+VsQJlVDc/Y/0cxqteF6JiIiIqCeYw90AIop+2dnZIa9/EUnHu9SVlpZiw4YNKC0tDXdTiCgKyJW4iIiIwomZGUTUZT19N553/0MrFJkuLPxJ1HvYbLZwN4GIiIiZGUREvV0oMl3UAREiurTl5OSEuwlERETMzCAi6u1ZBaHIdGHhT6Leozd+TxIRUeRhMIOIej0uJ9p1nPpDRERERD2J00yIqNdjQVEiIiIioujCzAwi6vWYVUBEREREFF2YmdHDzGaz5m8iIiIiIiIiCg6vqHtYXFyc5m8iIiKiSGe1WjV/ExERhRuDGT0sJycHDoeDy5oRERFR1Jg8eTIcDgcmT54c7qYQEREBYM2MHse5+URERBRtOH4hIqJIw8wMIiIiIiIiIooqDGYQERERERERUVRhMIOIiIiIiIiIogqDGUREREREREQUVRjMICIiIopQGzduxIgRI2Cz2ZCZmYn3338/3E0iIiKKCAxmEBEREUWg7du3Y9myZVi9ejXKysowYcIETJ06FadPnw5304iIiMKOwQwiIiKiCPTSSy9h0aJFuPfeezFmzBhs2LABqamp2LRpU7ibRkREFHYMZhARERFFmJaWFhw8eBCTJ0/WPD958mTs37/f8DXNzc04d+6c5g8REdGlisEMIiIiogjz5Zdfor29HUOGDNE8P2TIEHz++eeGr1m7di0cDofyJzU1tSeaSkREFBYMZhARERFFKJPJpHkshPB4Tnr88cfR0NCg/HE6nT3RRCIiorCIimDGqVOnsGjRIowYMQJ2ux1paWlYs2YNWlpafL7urrvugslk0vy58cYbe6jVRERERJ0zcOBAWCwWjyyMM2fOeGRrSHFxcejbt6/mDxER0aUqJtwNCMQnn3wCt9uN1157DSNHjkR5eTnuu+8+NDY24oUXXvD52ttvvx1btmxRHsfGxnZ3c4mIiIi6JDY2FpmZmdi9ezfuuOMO5fndu3dj5syZYWwZERFRZIiKYMbtt9+O22+/XXl85ZVX4ujRo9i0aZPfYEZcXByGDh3a3U0kIiIiCqnly5djwYIFyMrKwvjx47F582acPn0aixcvDnfTiIiIwi4qghlGGhoakJSU5He7PXv2YPDgwejXrx8mTpyIX/7ylxg8eLDX7Zubm9Hc3Kw5DgBWBCcioqgh+ywhRJhbQl0xd+5c1NXV4amnnkJtbS0yMjJQUFCAYcOGBfR6+flzDENERNEimDGMSUThSKeqqgrjxo3Diy++iHvvvdfrdtu3b0efPn0wbNgwnDx5Ej//+c/R1taGgwcPIi4uzvA1TzzxBJ588snuajoREVGPcTqduPzyy8PdDAqTTz/9lCuaEBFRVApkDBPWYEYggYMDBw4gKytLeVxTU4OJEydi4sSJ+N3vfhfU8WprazFs2DBs27YNs2fPNtxGn5nhdrtRX1+PAQMGeK0eHoxz584hNTUVTqeThbkM8Pz4x3PkH8+RfzxHvkX7+RFC4Pz580hJSYHZHBW1vqkbuN1u1NTUIDExsctjmGj/P9ETeI784znyj+fIN54f/6L9HAUzhgnrNJOHHnoIP/jBD3xuM3z4cOXfNTU1mDRpkjJvNFjJyckYNmwYjh8/7nWbuLg4j6yNfv36BX0sf1hl3DeeH/94jvzjOfKP58i3aD4/Docj3E2gMDObzSHPzInm/xM9hefIP54j/3iOfOP58S+az1GgY5iwBjMGDhyIgQMHBrTtZ599hkmTJiEzMxNbtmzp1J2muro6OJ1OJCcnB/1aIiIiIiIiIooMUZF7WlNTg1tvvRWpqal44YUXcPbsWXz++ecea6+PHj0a7733HgDgwoULePTRR/HBBx/g1KlT2LNnD2bMmIGBAwdqljgjIiIiIiIiougSFauZFBYWorKyEpWVlR7pkuqSH0ePHlVWH7FYLDh8+DDeeOMNfP3110hOTsakSZOwfft2JCYm9mj71eLi4rBmzRqvBUh7O54f/3iO/OM58o/nyDeeHyIt/p/wj+fIP54j/3iOfOP58a83naOoXM2EiIiIiIiIiHqvqJhmQkREREREREQkMZhBRERERERERFGFwQwiIiIiIiIiiioMZhARERERERFRVGEwowdt3LgRI0aMgM1mQ2ZmJt5///1wNyli7Nu3DzNmzEBKSgpMJhP+9Kc/hbtJEWft2rX45je/icTERAwePBizZs3C0aNHw92siLFp0yZcc8016Nu3L/r27Yvx48djx44d4W5WRFu7di1MJhOWLVsW7qZEjCeeeAImk0nzZ+jQoeFuFlHYcQzjHccwvnH84h/HMMHjGMZTbxzDMJjRQ7Zv345ly5Zh9erVKCsrw4QJEzB16lScPn063E2LCI2Njbj22mvx61//OtxNiVh79+7FkiVL8OGHH2L37t1oa2vD5MmT0djYGO6mRYTLL78czz33HEpLS1FaWoqcnBzMnDkTFRUV4W5aRDpw4AA2b96Ma665JtxNiTjp6emora1V/hw+fDjcTSIKK45hfOMYxjeOX/zjGCY4HMN419vGMFyatYfccMMNGDduHDZt2qQ8N2bMGMyaNQtr164NY8sij8lkwnvvvYdZs2aFuykR7ezZsxg8eDD27t2LW265JdzNiUhJSUl4/vnnsWjRonA3JaJcuHAB48aNw8aNG/HMM8/guuuuw4YNG8LdrIjwxBNP4E9/+hMOHToU7qYQRQyOYQLHMYx/HL8EhmMYYxzDeNcbxzDMzOgBLS0tOHjwICZPnqx5fvLkydi/f3+YWkXRrqGhAUBHZ0da7e3t2LZtGxobGzF+/PhwNyfiLFmyBNOnT8dtt90W7qZEpOPHjyMlJQUjRozAD37wA5w4cSLcTSIKG45hKNQ4fvGNYxjfOIbxrbeNYWLC3YDe4Msvv0R7ezuGDBmieX7IkCH4/PPPw9QqimZCCCxfvhzZ2dnIyMgId3MixuHDhzF+/Hi4XC706dMH7733Hr7xjW+Eu1kRZdu2bfjoo49w4MCBcDclIt1www144403cPXVV+OLL77AM888g5tuugkVFRUYMGBAuJtH1OM4hqFQ4vjFO45h/OMYxrfeOIZhMKMHmUwmzWMhhMdzRIF46KGH8O9//xslJSXhbkpEGTVqFA4dOoSvv/4aeXl5WLhwIfbu3cvBwH84nU48/PDDKCwshM1mC3dzItLUqVOVf48dOxbjx49HWloaXn/9dSxfvjyMLSMKL45hKBQ4fvGOYxjfOIbxrzeOYRjM6AEDBw6ExWLxuINx5swZjzsdRP78+Mc/xl/+8hfs27cPl19+ebibE1FiY2MxcuRIAEBWVhYOHDiAl19+Ga+99lqYWxYZDh48iDNnziAzM1N5rr29Hfv27cOvf/1rNDc3w2KxhLGFkSchIQFjx47F8ePHw90UorDgGIZCheMX3ziG8Y1jmOD1hjEMa2b0gNjYWGRmZmL37t2a53fv3o2bbropTK2iaCOEwEMPPYR3330XxcXFGDFiRLibFPGEEGhubg53MyJGbm4uDh8+jEOHDil/srKycOedd+LQoUMcBBhobm7GkSNHkJycHO6mEIUFxzDUVRy/dA7HMFocwwSvN4xhmJnRQ5YvX44FCxYgKysL48ePx+bNm3H69GksXrw43E2LCBcuXEBlZaXy+OTJkzh06BCSkpJwxRVXhLFlkWPJkiV466238Oc//xmJiYnKXTKHwwG73R7m1oXfqlWrMHXqVKSmpuL8+fPYtm0b9uzZg507d4a7aREjMTHRY45yQkICBgwYwLnL//Hoo49ixowZuOKKK3DmzBk888wzOHfuHBYuXBjuphGFDccwvnEM4xvHL/5xDOMfxzD+9cYxDIMZPWTu3Lmoq6vDU089hdraWmRkZKCgoADDhg0Ld9MiQmlpKSZNmqQ8lvO6Fi5ciK1bt4apVZFFLol36623ap7fsmUL7rrrrp5vUIT54osvsGDBAtTW1sLhcOCaa67Bzp078e1vfzvcTaMo8umnn2LevHn48ssvMWjQINx444348MMP+V1NvRrHML5xDOMbxy/+cQxDodAbxzAmIYQIdyOIiIiIiIiIiALFmhlEREREREREFFUYzCAiIiIiIiKiqMJgBhERERERERFFFQYziIiIiIiIiCiqMJhBRERERERERFGFwQwiIiIiIiIiiioMZhARERERERFRVGEwg4iIiIiIiIiiCoMZRBQxtm7din79+oW7GX794x//wNixY2G1WjFr1qxwN4eIiIjCjGMYop5nEkKIcDeCiAgAmpqacP78eQwePDjcTfHphhtuwNVXX421a9eiT58+UTF4ISIiou7DMQxRz2NmBhEFrKWlpVv3b7fbI34QAABVVVXIycnB5ZdfzkEAERFRFOAYpgPHMHQpYTCD6D9uvfVW/PjHP8ayZcvQv39/DBkyBJs3b0ZjYyPuvvtuJCYmIi0tDTt27NC87uOPP8a0adPQp08fDBkyBAsWLMCXX36p/Hznzp3Izs5Gv379MGDAAHznO99BVVWV8vNTp07BZDLh3XffxaRJkxAfH49rr70WH3zwgc/2NjQ04P7778fgwYPRt29f5OTk4F//+hcA4OzZsxg6dCieffZZZfv/+7//Q2xsLAoLCwEATzzxBK677jq89tprSE1NRXx8PL73ve/h66+/Vl5z1113YdasWVi7di1SUlJw9dVXAwA+++wzzJ07F/3798eAAQMwc+ZMnDp1Snndnj178K1vfQsJCQno168fbr75ZlRXVwMA/vWvf2HSpElITExE3759kZmZidLSUgDGKZqbNm1CWloaYmNjMWrUKPzhD3/Q/NxkMuF3v/sd7rjjDsTHx+Oqq67CX/7yF+XnX331Fe68804MGjQIdrsdV111FbZs2eL1vDY3N2Pp0qUYPHgwbDYbsrOzceDAAc1nVVdXh3vuuQcmkwlbt2413E9tbS2mT58Ou92OESNG4K233sLw4cOxYcOGgD5D9Wf0hz/8AcOHD4fD4cAPfvADnD9/XtlGCIH169fjyiuvhN1ux7XXXos//vGPnX7/REQUfTiG4RgG4BiGeiFBREIIISZOnCgSExPF008/LY4dOyaefvppYTabxdSpU8XmzZvFsWPHxAMPPCAGDBggGhsbhRBC1NTUiIEDB4rHH39cHDlyRHz00Ufi29/+tpg0aZKy3z/+8Y8iLy9PHDt2TJSVlYkZM2aIsWPHivb2diGEECdPnhQAxOjRo8Xf/vY3cfToUfHd735XDBs2TLS2thq21e12i5tvvlnMmDFDHDhwQBw7dkz85Cc/EQMGDBB1dXVCCCHy8/OF1WoVBw4cEOfPnxcjR44UDz/8sLKPNWvWiISEBJGTkyPKysrE3r17xciRI8X8+fOVbRYuXCj69OkjFixYIMrLy8Xhw4dFY2OjuOqqq8Q999wj/v3vf4uPP/5YzJ8/X4waNUo0NzeL1tZW4XA4xKOPPioqKyvFxx9/LLZu3Sqqq6uFEEKkp6eLH/7wh+LIkSPi2LFj4p133hGHDh0SQgixZcsW4XA4lOO/++67wmq1it/85jfi6NGj4sUXXxQWi0UUFxcr2wAQl19+uXjrrbfE8ePHxdKlS0WfPn2U87BkyRJx3XXXiQMHDoiTJ0+K3bt3i7/85S9efw+WLl0qUlJSREFBgaioqBALFy4U/fv3F3V1daKtrU3U1taKvn37ig0bNoja2lpx8eJFw/3cdttt4rrrrhMffvihOHjwoJg4caKw2+3iV7/6VcCf4Zo1a0SfPn3E7NmzxeHDh8W+ffvE0KFDxapVq5TjrFq1SowePVrs3LlTVFVViS1btoi4uDixZ8+eTr1/IiKKPhzDcAwjBMcw1PswmEH0HxMnThTZ2dnK47a2NpGQkCAWLFigPFdbWysAiA8++EAIIcTPf/5zMXnyZM1+nE6nACCOHj1qeJwzZ84IAOLw4cNCiP8OBH73u98p21RUVAgA4siRI4b7KCoqEn379hUul0vzfFpamnjttdeUxw8++KC4+uqrxZ133ikyMjJEU1OT8rM1a9YIi8UinE6n8tyOHTuE2WwWtbW1QoiOgcCQIUNEc3Ozss3vf/97MWrUKOF2u5Xnmpubhd1uF7t27RJ1dXUCgNIR6SUmJoqtW7ca/kw/ELjpppvEfffdp9nme9/7npg2bZryGID42c9+pjy+cOGCMJlMYseOHUIIIWbMmCHuvvtuw+PpXbhwQVitVvHmm28qz7W0tIiUlBSxfv165TmHwyG2bNnidT9HjhwRAMSBAweU544fPy4AKAOBQD7DNWvWiPj4eHHu3Dnl54899pi44YYblPbabDaxf/9+zT4WLVok5s2bF/T7JyKi6MQxDMcwHMNQb8RpJkQq11xzjfJvi8WCAQMGYOzYscpzQ4YMAQCcOXMGAHDw4EH8/e9/R58+fZQ/o0ePBgAlDbOqqgrz58/HlVdeib59+2LEiBEAgNOnT3s9dnJysuY4egcPHsSFCxcwYMAAzbFPnjypSf984YUX0NbWhnfeeQdvvvkmbDabZj9XXHEFLr/8cuXx+PHj4Xa7cfToUeW5sWPHIjY2VnPsyspKJCYmKsdNSkqCy+VCVVUVkpKScNddd2HKlCmYMWMGXn75ZdTW1iqvX758Oe69917cdttteO655zTt1Tty5AhuvvlmzXM333wzjhw54vXcJSQkIDExUTl3DzzwALZt24brrrsOK1aswP79+70er6qqCq2trZpjWq1WfOtb3/I4pi9Hjx5FTEwMxo0bpzw3cuRI9O/fX3kc6Gc4fPhwJCYmKo+Tk5OV9/bxxx/D5XLh29/+tmYfb7zxhrKPYN4/ERFFL45hOIbhGIZ6m5hwN4AoklitVs1jk8mkec5kMgEA3G638veMGTOwbt06j33JznzGjBlITU3Fb3/7W6SkpMDtdiMjI8OjEJWv4+i53W4kJydjz549Hj9Tz9c8ceIEampq4Ha7UV1drekwjcjjyr+Bjo5Vf+zMzEy8+eabHq8fNGgQAGDLli1YunQpdu7cie3bt+NnP/sZdu/ejRtvvBFPPPEE5s+fj/z8fOzYsQNr1qzBtm3bcMcdd/hskySE8HjO6HOT527q1Kmorq5Gfn4+/vd//xe5ublYsmQJXnjhBY9jif8s7hTIMX0RXhaJUj8f6Gfo673Jv/Pz83HZZZdptouLiwMQ3PsnIqLoxTEMxzCBHtMXjmEomjCYQdQF48aNQ15eHoYPH46YGM//TnV1dThy5Ahee+01TJgwAQBQUlISkuN+/vnniImJwfDhww23aWlpwZ133om5c+di9OjRWLRoEQ4fPqzcmQE67qzU1NQgJSUFAPDBBx/AbDYrRbK8HXv79u1KwSdvrr/+elx//fV4/PHHMX78eLz11lu48cYbAQBXX301rr76ajzyyCOYN28etmzZYjgQGDNmDEpKSvCjH/1IeW7//v0YM2aMz/OjN2jQINx111246667MGHCBDz22GOGHeHIkSMRGxuLkpISzJ8/HwDQ2tqK0tJSLFu2LODjjR49Gm1tbSgrK0NmZiYAoLKyUlOYLJDP0J9vfOMbiIuLw+nTpzFx4kSv2wX6/omIqPfgGIZjGCMcw1A04TQToi5YsmQJ6uvrMW/ePPzzn//EiRMnUFhYiHvuuQft7e1KpezNmzejsrISxcXFWL58eZePe9ttt2H8+PGYNWsWdu3ahVOnTmH//v342c9+plTVXr16NRoaGvDKK69gxYoVGDNmDBYtWqTZj81mw8KFC/Gvf/0L77//PpYuXYrvf//7GDp0qNdj33nnnRg4cCBmzpyJ999/HydPnsTevXvx8MMP49NPP8XJkyfx+OOP44MPPkB1dTUKCwtx7NgxjBkzBk1NTXjooYewZ88eVFdX4x//+AcOHDjgtWN/7LHHsHXrVrz66qs4fvw4XnrpJbz77rt49NFHAz5Xv/jFL/DnP/8ZlZWVqKiowN/+9jevx0tISMADDzyAxx57DDt37sTHH3+M++67DxcvXvQ4d76MHj0at912G+6//37885//RFlZGe6//37Y7Xbl7kggn6E/iYmJePTRR/HII4/g9ddfR1VVFcrKyvCb3/wGr7/+etDvn4iIeg+OYTiGMcIxDEUTZmYQdUFKSgr+8Y9/YOXKlZgyZQqam5sxbNgw3H777TCbzTCZTNi2bRuWLl2KjIwMjBo1Cq+88gpuvfXWLh3XZDKhoKAAq1evxj333KMsY3bLLbdgyJAh2LNnDzZs2IC///3vyp2HP/zhD7jmmmuwadMmPPDAAwA6ovizZ8/GtGnTUF9fj2nTpmHjxo0+jx0fH499+/Zh5cqVmD17Ns6fP4/LLrsMubm56Nu3L5qamvDJJ5/g9ddfR11dHZKTk/HQQw/h//2//4e2tjbU1dXhRz/6Eb744gsMHDgQs2fPxpNPPml4rFmzZuHll1/G888/j6VLl2LEiBHYsmVLUOcvNjYWjz/+OE6dOgW73Y4JEyZg27ZtXrd/7rnn4Ha7sWDBApw/fx5ZWVnYtWuXZq5oIN544w0sWrQIt9xyC4YOHYq1a9eioqJCmfPr7zMM1NNPP43Bgwdj7dq1OHHiBPr164dx48Zh1apVnXr/RETUO3AMwzGMNxzDULQwCW8To4jokvbEE0/gT3/6Ew4dOhTupvQKn376KVJTU5U5n0RERNQ5HMP0LI5hKFIxM4OIqBsUFxfjwoULGDt2LGpra7FixQoMHz4ct9xyS7ibRkREROQVxzAULRjMICLqBq2trVi1ahVOnDiBxMRE3HTTTXjzzTc9KnsTERERRRKOYShacJoJEREREREREUUVrmZCRERERERERFGFwQwiIiIiIiIiiioMZhARERERERFRVGEwg4iIiIiIiIiiCoMZRERERERERBRVGMwgIiIiIiIioqjCYAYRERERERERRRUGM4iIiIiIiIgoqvx/o2UrKOXjf+wAAAAASUVORK5CYII=",
"text/plain": [
"
"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sc.pp.highly_variable_genes(adata, min_mean=0.0125, max_mean=3, min_disp=0.5)\n",
"sc.pl.highly_variable_genes(adata)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Filter the genes to only those marked as highly-variable"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"adata = adata[:, adata.var.highly_variable]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Regress out effects of total counts per cell and the percentage of mitochondrial genes expressed. Scale the data to unit variance."
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [],
"source": [
"sc.pp.regress_out(adata, [\"total_counts\", \"pct_counts_mt\"])\n",
"sc.pp.scale(adata, max_value=10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Run PCA"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAisAAAHFCAYAAAAzCLlHAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAASYBJREFUeJzt3Xt8z/Xj///7i80YNmU0bAwRMpsQplpShJLUW9S7Axo5vTvLIe8oNYfePiSiE717Fx2mUjo79VaOIXtHiSE5hLrYZmPHx+8Pvz2/Xl4bz63Xa6/nttv1cnld2vP5fLyez8drT3rdPU5PlzHGCAAAwKEq+bsCAAAA50NYAQAAjkZYAQAAjkZYAQAAjkZYAQAAjkZYAQAAjkZYAQAAjkZYAQAAjkZYAQAAjkZYAVCoRYsWyeVyad++ff6uit8dOnRIkyZN0rZt2zyOTZo0SS6Xq/QrBVQghBUAherdu7fWrVunevXq+bsqfnfo0CFNnjy50LBy//33a926daVfKaACCfB3BQA4y6lTp1S1alXVqVNHderU8Xd1fKLgM3qjRSQiIkIRERFeqBWAotCyApQRH374oVwul1asWOFx7KWXXpLL5dL27dslSZs3b9aAAQMUFRWlatWqKSoqSgMHDtT+/fvd3lfQ1fPll19q8ODBqlOnjoKDg5WVlVVoN9BXX32lW265RREREapataouvfRSDRs2TMePH3c7b0HXyI8//qiBAwcqNDRUl1xyiQYPHqzU1FS3svn5+ZozZ45iY2NVrVo11apVS506ddKyZcvcyr3zzjvq3Lmzqlevrho1aqhHjx7aunXrBX9v5/uMu3fv1qBBg9SsWTMFBwerQYMGuvnmm5WcnGy9f/Xq1erQoYMkadCgQXK5XHK5XJo0aZLbZz33M02fPl0tWrRQUFCQ6tatq3vuuUe//fbbBesLwBNhBSgjbrrpJtWtW1cLFy70OLZo0SJdccUVatOmjSRp3759uuyyyzRr1ix98cUXmjZtmg4fPqwOHTp4BAtJGjx4sAIDA/Xmm2/q/fffV2BgYKF12LNnjzp37qyXXnpJX375pf75z39qw4YNuuqqq5STk+NR/rbbblPz5s2VlJSksWPH6u2339bDDz/sVua+++7Tgw8+qA4dOuidd97RkiVL1KdPH7eQ9Nxzz2ngwIFq1aqV3n33Xb355ptKT0/X1VdfrR07dtj6/RX2GQ8dOqTatWtr6tSp+vzzzzV37lwFBASoY8eO+vnnnyVJV1xxhfU7f/LJJ7Vu3TqtW7dO999/f5HXGj58uJ544gndcMMNWrZsmZ555hl9/vnniouLK/T3D+ACDIAy45FHHjHVqlUzJ06csPbt2LHDSDJz5swp8n25ubnm5MmTpnr16mb27NnW/oULFxpJ5p577vF4T8GxvXv3FnrO/Px8k5OTY/bv328kmY8++sg69tRTTxlJZvr06W7vGTFihKlatarJz883xhjzzTffGElmwoQJRdb9119/NQEBAWb06NFu+9PT0014eLjp379/ke+90Gc8V25ursnOzjbNmjUzDz/8sLV/06ZNRpJZuHChx3sKPmuBnTt3GklmxIgRbuU2bNhgJJnx48dfsB4A3NGyApQhgwcP1qlTp/TOO+9Y+xYuXKigoCDdeeed1r6TJ0/qiSee0KWXXqqAgAAFBASoRo0aysjI0M6dOz3Oe9ttt9m6/tGjR/XAAw8oMjJSAQEBCgwMVKNGjSSp0PP26dPHbbtNmzY6ffq0jh49Kkn67LPPJEkjR44s8ppffPGFcnNzdc899yg3N9d6Va1aVfHx8Vq9erWtuhf2GXNzc/Xcc8+pVatWqlKligICAlSlShX98ssvhX4eO1atWiXpTIvR2a688kq1bNmy0G48AOfHAFugDLn88svVoUMHLVy4UEOHDlVeXp7+85//6JZbbtHFF19slbvzzju1YsUKTZw4UR06dFBISIhcLpd69eqlU6dOeZzXzoyf/Px8de/eXYcOHdLEiRMVHR2t6tWrKz8/X506dSr0vLVr13bbDgoKkiSr7LFjx1S5cmWFh4cXed3ff/9dkqxxI+eqVMnev7kK+4yPPPKI5s6dqyeeeELx8fG66KKLVKlSJd1///2Ffh47/vjjjyKvV79+fY9xQwAujLAClDGDBg3SiBEjtHPnTqWkpOjw4cMaNGiQdTw1NVWffPKJnnrqKY0dO9ban5WVpT///LPQc9qZFfO///1PP/zwgxYtWqR7773X2r979+4Sf5Y6deooLy9PR44cKTIwhYWFSZLef/99qxWnJAr7jP/5z390zz336LnnnnPbf/z4cdWqVatE1ykIaIcPH/aYJXTo0CHr8wCwj24goIwZOHCgqlatqkWLFmnRokVq0KCBunfvbh13uVwyxlitGAVeffVV5eXllfi6BV/25553wYIFJT5nz549JZ2ZzVSUHj16KCAgQHv27FH79u0LfZWUy+Xy+DzLly/XwYMH3fad2yJ0Ptddd52kM0HobJs2bdLOnTvVrVu3EtcXqKhoWQHKmFq1aunWW2/VokWLdOLECT322GNuXSEhISG65pprNGPGDIWFhSkqKkpr1qzRa6+9VuLWAklq0aKFmjZtqrFjx8oYo4svvlgff/yxvvrqqxKf8+qrr9bdd9+tKVOm6Pfff9dNN92koKAgbd26VcHBwRo9erSioqL09NNPa8KECUpJSdGNN96oiy66SL///rs2btyo6tWra/LkySW6/k033aRFixapRYsWatOmjb7//nvNmDHDo0WkadOmqlatmt566y21bNlSNWrUUP369VW/fn2Pc1522WUaOnSo5syZo0qVKqlnz57at2+fJk6cqMjISI/ZUAAurFy1rMyfP9/fVQBKxaBBg3T06FFlZ2d7DOSUpLfffltdu3bVmDFj1K9fP23evFlfffWVQkNDS3zNwMBAffzxx2revLmGDRumgQMH6ujRo/r666//wic5M+165syZ+u6773T77berf//++uijj9S4cWOrzLhx4/T+++9r165duvfee9WjRw+NGTNG+/fv1zXXXFPia8+ePVt///vflZiYqJtvvlnLli3T0qVL1bRpU7dywcHBev311/XHH3+oe/fu6tChg15++eUiz/vSSy9p6tSp+vTTT3XTTTdpwoQJ6t69u7777juPcTwALsxljDH+roS3NGzYUL/++qu/qwEAALyozHUD9e/fv9D9xpgiBw8CAICyq8yFleXLl2vWrFmqUqWK235jjLW+AQAAKD/KXFiJjY1V27ZtC50BMHHiRD/UCAAA+FKZG7OyefNm1atXTw0aNPA4lpKSoiZNmvihVgAAwFfKXFgBAAAVS5mbujx16lTNmTPHY/+0adM0ffp0P9QIAAD4kqNbVvLz83Xo0CHVrFnTWj2zS5cu+vrrr1WtWjW3spmZmerevbvWrl3rj6oCAIBiMsYoPT1d9evXP+9zvhwdVn777TdFRkb6uxoAAMCHDhw44LFy9NkcPRuoZs2aks58iJCQEElSTEyMtm3b5vFQsvz8fMXGxmr79u2lXk8AAFB8aWlpioyMtL7vi+LosFIQSEJCQqyw0rVrV73wwgse05SfffZZde3a1SoHAADKhgs9+d3R3UBpaWkKDQ1VamqqFUKOHTum+Ph4hYSEqHPnzpKk9evX68SJE1qzZo3q1q3rzyoDAACbCvueL0yZmw1Up04dLVmyRP369dOuXbuUnZ2thIQEbdmyhaACAEA5VObCyrx58xQfH6+kpCStXbtW3bp10+DBgz1mBwEAgPKhTIaV5ORkbdiwQWvXrtXMmTP9XSUAAOBDZS6sBAYGWtOboqOjlZGR4ecaAQAAX3L0bKDCZGVlaefOnSoYF3zudqtWrfxZPQAA4GVlbjZQVFRUkVOcXC6XUlJSSrOKAACghOzOBipzLSv79u3zdxUAAEApKnNjVgAAQMVCWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI5GWAEAAI7m07DSp08fNWzYUFWrVlW9evV0991369ChQ768JAAAKGd8Gla6du2qd999Vz///LOSkpK0Z88e3X777b68JAAAKGdcxhhTWhdbtmyZ+vbtq6ysLAUGBl6wfFpamkJDQ5WamqqQkJBSqCEAACgtdr/nS23Myp9//qm33npLcXFxtoIKAACAVAph5YknnlD16tVVu3Zt/frrr/roo4+KLJuVlaW0tDS3FwAAqNiKHVYmTZokl8t13tfmzZut8o8//ri2bt2qL7/8UpUrV9Y999yjonqeEhMTFRoaar0iIyNL/skAAEC5UOwxK8ePH9fx48fPWyYqKkpVq1b12P/bb78pMjJS3333nTp37uxxPCsrS1lZWdZ2WlqaIiMjGbMCAEA5ZHfMSkBxTxwWFqawsLASVaogF50dSM4WFBSkoKCgEp0bAACUT8UOK3Zt3LhRGzdu1FVXXaWLLrpIKSkp+uc//6mmTZsW2qoCAABQGJ8NsK1WrZqWLl2qbt266bLLLtPgwYPVunVrrVmzhtYTAABgm89aVqKjo7Vy5UpfnR4AAFQQPBsIAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4GmEFAAA4WrkOK+vWrdPbb7+tI0eOuO1/4403/FQjAABQXOU2rLz44ou65557tHjxYkVHR2vp0qXWsdmzZ/uxZgAAoDh89tRlf3vllVf0/fffKyQkRDt27FDfvn2VmZmpv//97zLG+Lt6AADApnIbVowxCgkJkSS1atVKK1eu1A033KC8vDy5XC4/1w4AANhVbsNKQECAjh49qrp160qSIiIitGLFCl1//fU6cOCAn2sHAADsKrdjVsaOHas9e/a47atfv75WrFihgQMH+qlWAACguFzGwQM40tLSFBoaqtTUVKtLBwAAlA92v+fLbcvK1KlTNWfOHI/906ZN0/Tp0/1QIwAAUBLltmUlJiZG69evV7Vq1dz2Z2ZmKi4uTtu2bfNiTQEAQHFV+JYVSR5BRZKCg4M9pi5/8skn+vTTTyVJa9eu1YMPPqjXXnutVOoIAADOr9zOBjp58qSMMR7TlPPz85Wenm5tT5w4UV9++aWys7O1cuVKbdy4Ub1799aiRYt06NAhTZw4sbSrDgAAzlJuu4EGDx6sxo0be4SNZ599Vrt379bChQslSdHR0dq6dasyMzMVHh6uAwcOqHbt2kpPT1eXLl20fft2r30eAADw/9j9ni+3LSvTpk1TfHy8li9frs6dO0uS1q9frxMnTmjNmjVWucDAQAUEBCgkJETNmjVT7dq1JUk1a9ZU5cqV/VJ3AADw/5TbMSt16tTRkiVL1K9fP+3atUvZ2dlKSEjQli1brIXiJCkvL8/6ef78+dbPxhhlZ2eXap0BAICnchtW5s2bp/j4eCUlJWnt2rXq1q2bBg8e7DHodvLkycrIyJAkqwVGkn755RfdeeedpVpnAADgqdyOWWndurU+//xzRUREKDk5WcOHD9fatWt9VFMAAFBcFX7qcmBgoCIiIiSdGURb0HpyLhaPAwDA2crtANusrCzt3LnTWlPl3O1WrVpJkhYvXqz169d7vH/06NGKi4vTmDFjSq/SAADAQ7kNK5mZmerVq5fbvoJtl8ullJQUa7/dxeMAAEDpK7dhZd++fbbK2V08DgAA+Ee5HbNiV3x8vKZMmeKxPzExUfHx8UW+r2fPnr6sFgAA+P+V29lAdh07dkzx8fEKCQkpdPG4unXrqn///h7v++yzz6zA8u677/qkbgAAlGcVfgVbuwoWj/v888+1Zs0aRUVFKSEhQQMHDrTGsnzzzTfq3bu3rrnmGklnFoxbtWqVevfu7c+qAwBQIVT4lpV58+ZpwoQJat68uX766SctXLhQ/fr1cytz7NgxjRw5UuHh4Zo6daqCg4PVuHFj7d271yd1AgCgIqjw66zYNW/ePCUnJ2vDhg1au3atZs6c6VGmTp06evfdd9WxY0fFx8dr1apVHgNyAQCAb1T4biC7i8dJ0l133aWuXbtq6NChSktLK60qAgBQoVX4sGJ38bgC9evX1yeffFLq9QQAoKKq8GNWoqKiiuzSKVg8burUqapevbpGjx7tdnzatGlyuVyscgsAQAnY/Z6v8GHFjpiYGK1fv95jpdvMzEzFxcVp27Zt/qkYAABlGANsvYwl+QEA8A/Cig0FS/KfiyX5AQDwPcKKDSVdkh8AAPx1jFmxwc6S/AAAoHgYs+JFBUvy9+vXT7t27VJ2drYSEhK0ZcsWggoAAD5GWLFh3rx5io+PV1JSktauXatu3bpp8ODBhQ66BQAA3kVYscHOkvwAAMA3CCs2FGdJfgAA4F0Vfrl9O4q7JD8AAPAeZgPZYGdJfgAAUDx2v+dpWbFh3759/q4CAAAVFmNWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAo5VKWMnKylJsbKxcLpe2bdtWGpcEAADlRKmElTFjxqh+/fqlcSkAAFDO+DysfPbZZ/ryyy/1/PPP+/pSAACgHArw5cl///13JSQk6MMPP1RwcPAFy2dlZSkrK8vaTktL82X1AABAGeCzlhVjjO677z498MADat++va33JCYmKjQ01HpFRkb6qnoAAKCMKHZYmTRpklwu13lfmzdv1pw5c5SWlqZx48bZPve4ceOUmppqvQ4cOFDc6gEAgHLGZYwxxXnD8ePHdfz48fOWiYqK0oABA/Txxx/L5XJZ+/Py8lS5cmXdddddeuONNy54rbS0NIWGhio1NVUhISHFqSYAAHA4u9/zxQ4rdv36669uY04OHTqkHj166P3331fHjh0VERFxwXMQVgAAKL/sfs/7bIBtw4YN3bZr1KghSWratKmtoAIAACCxgi0AAHA4n05dPltUVJR81OMEAADKMVpWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAoxFWAACAo/k0rERFRcnlcrm9xo4d68tLAgCAcibA1xd4+umnlZCQYG3XqFHD15cEAADliM/DSs2aNRUeHu7rywAAgHLK52NWpk2bptq1ays2NlbPPvussrOziyyblZWltLQ0txcAAKjYfNqy8uCDD+qKK67QRRddpI0bN2rcuHHau3evXn311ULLJyYmavLkyb6sEgAAKGNcxhhTnDdMmjTpgoFi06ZNat++vcf+pKQk3X777Tp+/Lhq167tcTwrK0tZWVnWdlpamiIjI5WamqqQkJDiVBMAADhcWlqaQkNDL/g9X+yWlVGjRmnAgAHnLRMVFVXo/k6dOkmSdu/eXWhYCQoKUlBQUHGrBAAAyrFih5WwsDCFhYWV6GJbt26VJNWrV69E7y9rNm7cqG+//VZt27bVtdde6+/qAABQJvlsgO26dev0f//3f9q2bZv27t2rd999V8OGDVOfPn3UsGFDX13Wr7p162b9XNDl9fPPP2vo0KFasGCBH2sGAEDZ5bMBtkFBQXrnnXc0efJkZWVlqVGjRkpISNCYMWN8dUm/+/PPP62fZ82apS+//FItWrTQ77//rh49emjYsGF+rB0AAGWTz8LKFVdcofXr1/vq9I7kcrmsnzMyMtSiRQtJ0iWXXKJKlXiyAQAAJeHzReEqkr1796p///4yxujgwYM6ffq0qlatKknnXV8GAAAUjbDiRbNmzbJ+vummm5SZmamqVavq0KFDuuWWW/xXMQAAyrBir7NSmuzOvwYAAGWP3e95BlJ40dSpUzVnzhyP/dOmTdP06dP9UCMAAMo+Wla8KCYmRuvXr1e1atXc9mdmZiouLk7btm3zT8UAAHAgWlb85NygIknBwcFycCYEAMDRCCtedPLkyUJDSX5+vtLT0/1QIwAAyj7CihfFx8drypQpHvsTExMVHx/vhxoBAFD2MWbFi44dO6b4+HiFhISoc+fOkqT169frxIkTWrNmjerWrevnGgIA4ByMWfGDOnXqaMmSJerXr5927dql7OxsJSQkaMuWLW5BZcOGDUpLS5MknT59WmPHjtV1112nhx9+mO4iAADOQVjxonnz5ik+Pl5JSUlau3atunXrpsGDB3sMuh00aJC1su0TTzyhw4cP6/HHH1d+fr4eeOABf1QdAADHYgVbL5o3b56Sk5MVERGh5ORkDR8+XP369fMoV6lSJVWpUkWS9N///lebN29WpUqV1LNnT8XExJR2tQEAcDRaVrwoMDBQERERkqTo6GhlZGQUWi44OFi7du2SJNWuXVsnT56UdOb5QTk5OaVTWQAAyghaVrwoKytLO3futKYvn7vdqlUrSdKMGTPUo0cP3X333YqJiVG3bt3Uu3dvrVq1SoMGDfJb/QEAcCJmA3lRVFSUXC5XocdcLpdSUlKs7UOHDumll17Sjh07lJOTo0aNGunOO++0ZhEBAFDe2f2eJ6wAAAC/YOqyg/HAQwAA7KNlxQ944CEAALSsOB4PPAQAwB7Cih/wwEMAAOwjrPgBDzwEAMA+xqz4AQ88BACAMSuOZveBhwAAgLDiF3YfeAgAAAgrflHwwMMNGzZo7dq1mjlzpr+rBACAYxFW/MDuAw8BAAAPMvQLuw88BAAAzAbyi+I88BAAgPLK7vc8LSt+sG/fPn9XAQCAMoMxKwAAwNEIKwAAwNEIKw722GOP8QRmAECFR1hxsFdeeUXdu3fXFVdcoblz5+rEiRP+rhIAAKWOsOJgTZo00cGDBzV27Fh9/PHHioyM1F133aVVq1b5u2oAAJQawoqDuVwuBQYGqn///vr888+1Y8cOtWjRQvfff7+aNm3q7+oBAFAqCCsOdu4SOJGRkZo4caL27NmjV155xU+1AgCgdBFWHGzAgAFFHrvuuutKsSYAAPgPK9gCAAC/sPs9T8uKg02dOlVz5szx2D9t2jRNnz7dDzUCAKD00bLiYDExMVq/fr2qVavmtj8zM1NxcXGswQIAKNNoWSknzg0qkhQcHOwx+BYAgPKKsOJgJ0+eLDSU5OfnKz093Q81AgCg9BFWHCw+Pl5Tpkzx2J+YmKj4+Hg/1AgAgNLHmBUHO3bsmOLj4xUSEqLOnTtLktavX68TJ05ozZo1qlu3rp9rCABAyTFmpRyoU6eOlixZon79+mnXrl3Kzs5WQkKCtmzZQlABAFQYhBUHmzdvnuLj45WUlKS1a9eqW7duGjx4cKGDbgEAKK8IKw42b948JScna8OGDVq7dq1mzpzp7yoBAFDqCCsOFhgYqIiICElSdHS0MjIy/FwjAABKX4C/K4CiZWVlaefOndb05XO3W7Vq5c/qAQBQKpgN5GBRUVFyuVyFHnO5XEpJSdGff/6p6dOnKyIiQqNGjdLkyZO1evVqtW7dWs8884xq1apVupUGAMAmu9/zhJUyrm/fvqpbt65Onjxp/b7uuusuffLJJ/rjjz/0zjvv+LuKAAAUyu73PN1AZVxKSoo+/PBD5ebmql69ejp8+LACAgLUs2dPtWnTxt/VAwDgL2OAbTlhjFFubq41nsXlcik/P9/PtQIA4K+jG6iMGzRokDIyMpSZmamaNWuqUqVKuvXWW/XFF1/oxIkTeu+994p1vvz8fL322mtavHixfv31V0lSw4YNNWDAAA0ZMkSVK1f2xccAAFRAdANVEPPnz9eCBQskScOGDdPXX3+t+fPnq0mTJpo+fXqxzzdixAgdPnxYjz/+uKKiomSM0f79+zV//nx9//331rUAACgttKzATbNmzfTLL7947DfGqFmzZtq9e7cfagUAKI94NlAFMXXqVM2ZM8dj/7Rp0zxaVpKSktS3b1/FxsaqQ4cOuvfee5WcnOxWxuVy6fjx4x7nO378eJHTqAEA8CW6gcq4xYsXa/369R77R48erbi4OI0ZM0aSNH78eG3ZskXdu3fXH3/8oS5duuiiiy5Snz59NGvWLN1yyy2SpCeeeEIxMTHq16+fGjVqJJfLpb179+rDDz/U5MmTS/WzAQAg0Q1U5sXExOiHH3644LGWLVvqf//7nypXrqz09HT17dtXK1as0N69e3Xrrbdq27Zt1vv27dun999/322A7W233abGjRv7/PMAACoOBthWECdPnpQxxqOLJj8/X+np6dZ2YGCgNZOnatWq+vPPPyVJjRs3Vm5urtt7o6Ki9Nhjj/m45gAA2MOYlTIuPj5eU6ZM8difmJio+Ph4a7t9+/YaMmSI3n//fd11113q2rWrJOnUqVPKycmxda358+fbKjdhwgRb5QAAsINuoDLu2LFjio+PV0hIiDp37ixJWr9+vU6cOKE1a9aobt26kqTMzEw999xz+uGHH9ShQwc98cQTCgoKUnp6ulJSUhQTE3PBazVs2NDqGvJGOQBAxcazgSqQ7du36/PPP9eaNWsUFRWldu3aaeDAgapWrVqxz9W/f/9C9xtj9Nlnn+nkyZOSpCuvvLLIcsnJyTp9+nSxrw0AqFgYs1JBzJs3TxMmTFDz5s31008/aciQIerXr59HualTp6p69eoaPXq02/5p06bJ5XJZs4aWL1+uWbNmqUqVKm7ljDFatWqVtb1r1y4tXrxYwcHBHuXuuOMOb308AAAIK2XdvHnzlJycrIiICCUnJ2v48OGFhhW7U5xjY2PVtm1btW/f3qPsxIkTrZ/btm2r0NBQxcXFeZQ7N+gAAPBX+HyA7fLly9WxY0dVq1ZNYWFhhX6RouQCAwMVEREhSYqOjlZGRkaRZQvrFgoODtbZPYGzZ89WvXr1Cn3/mjVrrJ8XLVqkli1bFlpuz549tuoOAIAdPm1ZSUpKUkJCgp577jldd9111ngGeE9WVpZ27txpBY5zt1u1aiXJ/hTnwlpUCjRp0sT6uVGjRkWWo2UFAOBNPhtgm5ubq6ioKE2ePFlDhgwp0TkYYHthUVFRRS6D73K5lJKSIkkaPHiwGjdu7NaVI0nPPvusdu/erYULF0qyP7blfObPn68HHnigJB8HAFCB+H020MaNG9WxY0e9/vrreuGFF3TkyBHFxsbq+eef1+WXX17oe7KyspSVlWVtp6WlKTIykrDiBXanOMfExGj9+vUeXUaZmZmKi4tzW+m2KExdBgDY4ffZQAX/op80aZJmzpypqKgo/etf/1J8fLx27dqliy++2OM9iYmJPH/GR+rUqaMlS5a4TXFOSEgodIqznbEt55viXLA6LgAA3lDsAbaTJk2Sy+U672vz5s3Kz8+XdGY109tuu03t2rXTwoUL5XK59N577xV67nHjxik1NdV6HThw4K99OljmzZun+Ph4JSUlae3aterWrZsGDx7sEUwKxrac69yxLcuXL9cNN9yg3r17e7yqVq1qlVuyZImee+45jxaZxMRE735AAEC5VeyWlVGjRmnAgAHnLRMVFWV9sRUM8JSkoKAgNWnSpMgugqCgIAUFBRW3SrDB7hTnguX7zx3bcu7y/XamOD/55JNas2aNYmNj1bNnT40dO1YPPvigJOm9997TuHHjvPkRAQDlVLHDSlhYmMLCwi5Yrl27dgoKCtLPP/+sq666SpKUk5Ojffv2nXcmCXzD7hTnadOmKT4+XsuXLy90bEsBO1Ocly1bps2bN6tKlSqaMGGC+vTpo4yMDI0fP77Q1hsAAArjszErISEheuCBB/TUU08pMjJSjRo10owZMyRJf/vb33x1WRTB7hRnu2Nb7Exxzs/Pt6Yxh4eH6+uvv1bPnj2Vl5fnNoMpJydHCxYsUKVKlTR06FAtXbpUb731ltq0aaOJEycyFRoAKjifPhsoJydH48aN05tvvqlTp06pY8eOmjVrVpGzgc7F1GXvsTvF+dzl+xcuXFji5fs7deqkJUuWKCoqyjp+8uRJ9ejRQ5s3b7Zmfg0fPlxHjx7VqVOnFBoaqqysLN1xxx1aunSpwsPDNXv2bC/9FgAATuL3qcveQFgpfa1bt9bnn3/uNrZl7dq1HuXsTHFetWqVLr74Yo8nOmdkZGj27NkaP368pDPdUgUPP6xbt66OHDmi4OBgZWdnq127diwkCADllN+nLqNs8uby/V27di30fdWrV7eCSsE1Jalq1apq2rSp9XDEKlWqKCCAP6IAUNH5/NlAKFsKxrLs2LFDO3bs8NguYHeKc1Hmz59v/WyMUV5eniTpww8/tPbn5uYqNzf3vOeJjo6+4LUAAGUb3UBw4+3l+4ty9iq369evV2xsrNv6LNKZhQVXr16twYMHS5KuvPJKj/Ns375dbdq0kXRm1WQAQNnBmBX4lJ3l+8+3yu1nn32mkydPFuuarVq1UpcuXXT33XfLGCNjjAYOHKglS5ZIkts6MAAA5yOswOe2b9/uNsW5Xbt2blOcq1evrlmzZnlMPTbG6LHHHtPx48cveI2zH4qYlZWlCRMmaPfu3Zo/f77Cw8PVpEkTq7WnQFZWltvigu+//77++9//qm3btrrvvvv+4qcGAHiL7e9542CpqalGkklNTfV3VXCOuXPnmlq1apkrr7zShISEmKSkJI8ycXFxZtOmTYW+PyIiwtZ1IiMjPfZ98803pn379mbhwoWmcePGHsfbtm1r/Txv3jzTpk0bk5iYaLp06WKefvppW9cFAPie3e95WlZQInamOG/evFn16tVTgwYNPN6fkpJiLR5Xku6izMxMPf7441q7dq1++OEHt2Nt27bV1q1bJUkdO3bU0qVL1aBBA6WnpysuLo6p0ADgEExdhk/ZmeJsZ5Vb6cxDEYvqLlq1alWh7w8ODtbcuXMLPXb2AOG8vDwrLNWsWZOp0ABQBvF/bpSIneX77axyK9l7KOKFnD225eeff9aVV14pY4x++eUXpaenq2bNmjLGKCcnp8SfGQDgH3QDoUTsTHG2s8qtZL+76HzOngp99gMXpTMP1axRo4aOHDmipKQkjRw5UtKZNV2uvvpq1a5dW8ePH9fIkSO1adMmxcTEaO7cuapfv/4FrwsAKDlmA8HvYmJiPMaT2DlWFG9PhW7ZsqV+/PFHVapUSffdd58iIiJ011136bPPPtPKlSv1ySef6Pbbb9ff//533XzzzapcuXKxzg8AOD+73/OsYAufsbvK7dSpUzVnzhyPctOmTdP06dOt7eXLl+uGG25Q7969PV7nLihXlLNXznW5XKpU6cxfgR9++EFTpkxRy5Yt9cgjj+jgwYOSzrTS/POf/1SDBg30+OOP66effir0vCkpKRowYIAee+wxZWZmatCgQWrcuLFuvvlmq8UHAFAyhBX4THx8vKZMmeKxPzEx0W0Bt8WLF+v+++/3KDd69Gi9/fbb1nbB2JZ7773X7XXfffcV+pyiwjz33HPWz5dcconWrVsnSYqIiNCRI0ckSenp6dbYloiICG3fvl3Lli1Tenq6OnfurC5duuj11193G1SckJCg6OhoVa5cWfHx8apbt66WL1+uLl26aNiwYVa53NxcPf/884qNjVWtWrVUp04dxcfH69NPP7VVfwCoiOgGgs/YWeVWst9d5O2p0MnJybr11lsVFxenoKAgffnll+ratas2bdqkcePG6e9//7uuuOIKbdmyxTrHqVOn9P777+v111/Xli1blJqa6lZPY4waNGigQ4cOWe+JjY21xuckJCTIGKNevXppyZIlat26tS677DJNnTpVI0aMUEJCgqQzrU+vvfaaFi9ebLXMNGzYUAMGDNCQIUMu2CU1YcIEPfvss+ctAwD+xpgVOMKFVrmVpKZNm2r37t0eA3bz8/N16aWXeqxQeyHFWTn31KlTevvtt7Vjxw7l5OSoUaNG6t+/vyIjIyW5r9lyrrMDUuvWrbV9+3ZlZmYqPDxchw8fVs2aNZWbm6vo6Gjt3LlT0plZUgUPhMzOztb111+vb775RsePH9c111xjHXvggQd0+PBhPfDAA4qKipIxRvv377dW7l2wYMF5fwdnDziWpCVLliglJUW9evVSbGystT8xMVHjxo2z82sFAK9jnRX43bx58zRhwgQ1b95cP/30k4YMGaJ+/fp5lCvoLjp3mvK53UW+mApdrVo1DRkypMjP8MgjjxR57OxZSr169VKXLl2UlZWloUOHasCAAeratatWrFihq6++2ipXuXJl5eXlqXLlysrOzrZaecLCwtxaS1asWKFffvnF7XqtWrXSjTfeqGbNmkkq/MGO0plQdvToUWv7ySef1Jo1axQbG6uePXtq7NixevDBByVJ7733nhVWcnJytGDBAlWqVElDhw7V0qVL9dZbb6lNmzaaOHGiR/g7W3R0NIvtAfAZWlbgM3ZWuZWK111UWlOhz163xa5ly5ZJkm6++WZt2bJFixYtUpMmTTRy5Ejri/6xxx7T//73P11//fX66KOP1KNHDz355JNKTU1Vp06drBaY5s2b67vvvlNYWJjbNY4dO6a4uDj98ssvqlWrlhYvXqzg4GC3MsYY3XHHHfr9998lSW3atNHmzZtVpUoVHTlyRH369FHfvn01fvx4t5aj4cOH6+jRozp16pRCQ0OVlZWlO+64Q0uXLlV4eLhmz54tyd7Tr+0+n8nu9HG7s7JSUlI0fvx4RURE6Omnn9bIkSO1evVqtW7dWnPnzlXDhg0vcBcBlCaeDQS/i42NPe/22X744Qczbdo006tXLzNixAjz2muvmczMTLcybdq0KfL95ztWEoU9k+hcL730UrHPm5+fb1555RUzatQo88Ybb1j7T58+bX7//Xdr+9VXXzX169c3o0aNMjNmzDDPP/+8GTlypGnQoIF59dVXjTHGXHvttebbb78t9DpnP3vp8ssvdzuWmppq4uLizNNPP+32HKXWrVsbY4w5deqUqVmzpsnIyDDGGJOVlWUdM8aYli1bmvvvv9+sWbPGrF692qxatcqEh4eb1atXm9WrVxtj7D+fqUWLFiYvL88YY8y9995rJkyYYHbs2GH+9a9/md69e1vlwsLCTHR0tLnkkkvMY489Znbu3Fno577uuuvMlClTzJgxY0z79u3NmDFjzI8//mgSExPNjTfeaIwxJicnx8yYMcPExMSY0NBQExYWZq655hqzfPlyt3Pl5eWZl19+2XTt2tU0bdrUNG3a1HTt2tUsWLDA5ObmFnr9c40fP95WOaCisvs9T1iBz7Rs2dLs2LHD/Pjjj+bHH3/02C5g56GIxhjTpEkTk5+f77E/Ly/P7YGGiYmJ5oUXXvAoN3XqVDNt2jRr+29/+1uhr9tvv91Ur179gp/PTqAxpmShxhhjUlJSzIwZM8zo0aPN6NGjzYwZM0xKSop1fN++febPP/8s9L1ZWVnWzx07djR79+51O56enm7i4uJMlSpVrH1nB4zzBc3Tp0+bRx991Nxyyy3m8OHDxhjj8UDJs8tfeeWV5rfffjPGGJOWluYRfOxcs+DnDRs2mGHDhplatWqZuLg489prr5mTJ09a5QpCa35+vqlXr57b+WJiYowxxtx///1myJAhJikpyfztb38zkydPNkuWLDGxsbHm5ZdftsoPGzbM9OnTx3z66afWn9tPP/3U9OnTxwwdOtTYce6fkcWLF5tnn33WbN261W3/c889Z4wxJjs728yZM8fMnTvX5OTkmHfeecf06dPHPPnkk273FCgvCCvwu0aNGpmoqKhCX2d/uV1++eXmwIEDxhhjtm/fbrp06VLo+QYNGlToU5OnTJli7rvvPmu7TZs2Hq0yxhiTkZFhfWEZY0xwcLB5+eWXzaJFi9xeCxcuNLVr1zbG/PVAY4zvQ82FrFy50mzbts1j/8mTJ82zzz5rbcfGxlotBvv27bP25+TkuAWMAud7+vXZwaddu3Zux84OIddee6357rvvjDHG3HTTTVb4SUtLc2sROvt8xhiTmZlp/v3vf5trr73WhISEWPsvv/xyk5eXZ9LT00316tVNWlqa9RlatGhhjHEPSFlZWebqq682xhhz7Ngxt2OXXnqpx2c25kwQatq0qbXdoUOHQl/t27c3QUFBVrkJEyaYq666yowaNcqEh4ebWbNmeXy+Bx54wPTr18/07NnTDBgwwNx6661myZIlpn///uYf//hHofUpUNg9On36tNv2e++9Z/7xj3+YhQsXuu3/4IMPzPHjx63fQ//+/U3jxo1N3759zcGDB61yt912m/nggw9stywBF0JYQZlht7vo6NGjpmXLlqZjx47moYceMg899JDp1KmTadGihVsXit3uori4OLNp06ZCyxV0o9gJNMb4JtR899135q233jJHjhxx279o0SKPMgVf8oWVsWvdunXm1KlTHvv37NljXnvttULfk5GRYUaMGOHxOw8ODra+sENCQqzQkJ+f7xZCtm/fbpo2bWruvvtuc//995uGDRuae++917Rq1cq8+eabVrnzdSHu2bPH+vnxxx83nTp1Mm3btjUPP/yw6dWrl5kxY4a58cYbTUJCgjHmzJd6wZdtenp6oV1hxhjTrFkzc+zYMY/rHT161C3IhIaGmk8//dTqBit4rVq1ytStW9cqFx0dbbWOHD582HTo0MEKiwWfz25XXGHhKCgoyPq5gD+64/bs2WPuuOMO8+ijj5qMjAxz3333maioKHPTTTeZ/fv3W+XojoMx9r/nmQ0Ev7PzUERJqlOnjpYsWeI2FTohIcFjKnTByrmFTYU+e+Xc2bNnq169eoXWqeD5QnZnFtl9cvT51oD5888/re0XX3xRs2fPVosWLfTggw9qwYIF1kyq2bNn695777VVpsC6deu0d+9eXXfddQoPD7f2v/HGG1a5Tp06FVq3Jk2aFDkouainX5+7yF3Bvfj99981fPhwa3/BLKKC6eO33HKLGjVqpGeeecaaPi7Zn5U1ffr0Qgc6d+/e3XomVI8ePdS7d29rkHPB7yw1NVW5ubnWuZ544gnFxMSoX79+atSokVwul/bu3asPP/xQkydPtsq1bdtWoaGhiouL86jb2X8e8vPzre3w8HB9/fXX6tmzp/Ly8qzfT2BgoCSpatWqatq0qTV4ukqVKm5PDD958qS6dOmiu+++W+bMPzo1cOBAzZgxw+365qz5E4sWLdKnn36qBg0aaOTIkYqLi7P+DJ+7mvOiRYsknXkkxZtvvmmdIyIiQlu3btXGjRv1+uuvq3PnzmrVqpWGDBmiO+64Q9WrV1dCQoKuu+46paWlKT4+Xtddd52WL1+uZcuWadiwYfrss88knRnQbYzRP//5T7c1hyZMmKCDBw9aaw6NGDFChw8f1uOPP+4xjf/777+/4DR+SXrzzTetdYfsTOP39sw4uwPOcR4+jUx/ES0rFYPd7iK7Y1vsdhfZsWnTJmu8xbnO/he9nVYaY+y31LRp08b6c//jjz+aZs2aWS0NBf8Kt1PGGGPmzJljLr30UnPTTTeZsLAwt9/bud0rdltqvNmiU9rsDnI25sLjhozx7tghu11xdsYNGeOf7jg744aM8U93nJ2uOGPsd8d5s4XL211x3mzh8kbrVlHoBkK5Y3dsi93uIrsDce3wdqiJjo52O3bgwAHTokULs2jRIut/fHbKGOP9UGO3nLeDT1kPUnbGDhW3K+5844aM8U93nJ1xQ8b4pzvOTlfc2df3xsy4c89d1IBzf8yMM8begHNvDDYvCmEF5Y4vpkJfaCCuNwONMfZDTdu2bT3+lX/w4EHTsmVLU6NGDdtljPF+qLFTztvBx19B6uOPP7b+hfnf//7X/OMf/7CmjpeknC8UNW7IGOPxxZ2enm6MOfNF/eKLL7qVzczMNK+++qp55JFHzOjRo83zzz9vfv31V7cy//73vy9YHzvjhowx5tFHHzU9evQwM2bMMFdddZV55plnjDHGnDhxwi3U2JnGb4y9qfx2p/F7c2bcuecrqoXLHzPjzr1uUS1cdlu3SoKwgnLH21Oh7QzEtTuzyNuh5p133rGa5c926NAh63/4dsoY4/1QY6ecN4OPt8vZDTRPPvmkufLKK01sbKx59NFHzdVXX22mTp1qrrrqKrduRrvljLEXavwZfLzlo48+Mh999JHJz883mzdvNqNGjTIzZ8506xYrTnfc3r17vdIdZ3cavzdnxhljr4XLHzPjjLHXwmW3daskCCsod7w9FdrOui12Zxb5K9TY4e1QY6ecN4OPt8vZDT6tW7c2OTk5JjU11VSrVs0aT5CWluZ2Hbvl7IQabwcff5ZzIrvT+L05M84Yey1c/pgZZ4y9Fi67rVslQVhBhWW3u8jOQFy7C9H5I9T4ozXHbjlvBh9vl7MbfM7tIjvb2X+m7JazE2q8GXz8Wa6wQFPYl3xJg8/ZYy68Uc4pvNUVV8CbLVx2WrdKgrCCCstud5Gdgbh2Zxb5I9SUh9Ycf5SzG3zOvh9nnzM/P9+0atWq2OXshBpvBh9/lfN28Cls/aIaNWpYPxe33Lm8HXzKa5DyFdZZQYWVmZmpXr16ue0r2Ha5XEpJSZFkb92WadOmKT4+XsuXLy/0IYsF7D452u4aMJI8HtgonVnXxJy1doadMosXL9b69es9yo0ePVpxcXHW06rtPtXaTrmi1pOpV6+eXn75ZWvbH+XGjh2rPXv2WA/IlKT69etrxYoVeuqpp6x9kydPVkZGhqpXr27de0n65ZdfdOeddxa7XF5envXz/PnzrZ+NMcrOzrZdRjqzHktAQIBCQkLUrFkz1a5dW5JUs2ZNtwc9+qPchx9+qK1btyozM1Ph4eE6cOCAateurREjRqhLly7W3xG75b755hv17t1b11xzjfW7WLVqlXr37q2z2SlX2J+PtWvXWvvfffddv5Y7V8+ePa11ac6ntMtt3LhR3377rdq2batrr732gufzBsIKyp19+/bZKjdv3jxNmDBBzZs3108//aQhQ4ZYC4QVsLsQnT9CjbeDj91QY6ecN4OPt8vZDT59+/YttFzz5s01YcKEYpezE2q8GXz8Vc7bASk5OVkjR47U999/r6lTpyo4OFiTJ092W/TQbjlvBh9vl3NykOrWrZtWrFghSUpKStLDDz+sXr166aWXXtKjjz6qYcOGeZzD63zSruMldAPBl+wMxLU7s8gYe9Ol7a4BY6f7yZ9dVBcqZ7eLyh/l7HaLebucN33wwQdu01QL/Pzzz2bKlCl+Left7rMC//nPf0z79u3NypUrC51xY6fc0aNHzd/+9jczevRoa+2UqKgoj3P4o9wll1xiBg8e7LZQZFhYmLXtz3Jnd0FeddVV1touR44ccfv7VxKMWQEuwM5AXLszi/wRarwZfIyxH2q8PYuqtMv5K0j5Y9C0P3g7IJ3t4MGDpnfv3m6rPZeknDeCj7fLOTlIFTXwvLDt4iKsABdgZyCu3ZlF/go1pd2aY7ecN4OPt8v5M0iV9qDp8tCK5CveCj7eLufEIFWrVi3rwax169Z1m9Z97kJ7xUVYAS7AzrotdmcW+SPU+KuLypuzqPxRzl9Bym6rz4XKFPxMNxtBqrTKnfsssz/++MMq/1efaE1YAbzA7kJ0/gg1/mrNsVPOm8HH2+X8FaTshBp/jC/ydjmCVPkOUt5GWAFKkT9Cjb+6qOyW81bw8XY5fwUpfwyarmjdbBcqR5AqWTlfhi3CCuBA3gw1/uqi8uYsKn+V80eQ8seg6YrUzUaQ8l05u+cqCcIKUIbZCTX+6qLy5iwqf5TzV0AypvQHTVekbjaClG+DlJ1zlQRhBYAxxvuhxpuzqPxRzl9Byh+Dpv1RjiBV/oKU3XOVhN3veZcxZy1h6TBpaWkKDQ1VamqqQkJC/F0doFyLioryWA23wNmPKbBTrlWrVkpKSrJWyL399tvdtlu1amX9t7TLtW3bVlu3brXqfO52Ufv/arnWrVvr888/V0REhJKTkzV8+HCtXbu22GUkz9WXFy5c6LH6sj/Lbd++3W3V53bt2nms+uyPcseOHVN8fLxCQkIKXWm64BEM/ig3ePBgNW7c2GN162effVa7d+/WwoULJckv5eyeqyTsfs8TVgB4nTeDj7fL+StI2Qk13gw+/ipHkCp/QcruuUrC9vf8X2q/8TG6gQB4m91uMW+X88eg6YrUzebk8UpOHvhtt5zdcxUXY1YAwEH8MWjaH+UIUiUr5+QgVZyxVMVl93uepy4DQCmw8zRwu08Mz8zMVK9evdz2FWyf3S3mj3JZWVnauXOn1Q127nZBt5g/ygUGBioiIkKSFB0drYyMjEJ/v/4oN2/ePCUnJ7t1sRXV3VXa5eyey5cIKwBQxtgNNf4oR5Aqf0HK7rl8ibACAPAaglT5C1J2z+VLzAYCAMDPnDyDzu65SoKpywAAwNHsfs9XKsU6AQAAFBthBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOBphBQAAOJrPwsrq1avlcrkKfW3atMlXlwUAAOVMgK9OHBcXp8OHD7vtmzhxor7++mu1b9/eV5cFAADljM/CSpUqVRQeHm5t5+TkaNmyZRo1apRcLpevLgsAAMoZn4WVcy1btkzHjx/XfffdV2SZrKwsZWVlWdtpaWmlUDMAAOBkpTbA9rXXXlOPHj0UGRlZZJnExESFhoZar/OVBQAAFUOxw8qkSZOKHDhb8Nq8ebPbe3777Td98cUXGjJkyHnPPW7cOKWmplqvAwcOFLd6AACgnCl2N9CoUaM0YMCA85aJiopy2164cKFq166tPn36nPd9QUFBCgoKKm6VAABAOVbssBIWFqawsDDb5Y0xWrhwoe655x4FBgYW61rGGEmMXQEAoDwq+H4v+L4vis8H2K5cuVJ79+69YBdQYdLT0yWJsSsAAJRj6enpCg0NLfK4y1wozvxFd955p/bv369vv/222O/Nz8/XoUOHVLNmzSKnO6elpSkyMlIHDhxQSEjIX60uSoj74H/cA2fgPvgf98AZ7NwHY4zS09NVv359VapU9DBan7esvP322yV+b6VKlRQREWGrbEhICH8oHYD74H/cA2fgPvgf98AZLnQfzteiUoBnAwEAAEcjrAAAAEcr82ElKChITz31FFOe/Yz74H/cA2fgPvgf98AZvHkffD7AFgAA4K8o8y0rAACgfCOsAAAARyOsAAAARyOsAAAARyvTYWXevHlq3Lixqlatqnbt2um///2vv6tUrn3zzTe6+eabVb9+fblcLn344Ydux40xmjRpkurXr69q1arp2muv1Y8//uifypZTiYmJ6tChg2rWrKm6deuqb9+++vnnn93KcB9876WXXlKbNm2sxa46d+6szz77zDrOPSh9iYmJcrlceuihh6x93AffmzRpklwul9srPDzcOu6te1Bmw8o777yjhx56SBMmTNDWrVt19dVXq2fPnvr111/9XbVyKyMjQzExMXrxxRcLPT59+nTNnDlTL774ojZt2qTw8HDdcMMN1jOe8NetWbNGI0eO1Pr16/XVV18pNzdX3bt3V0ZGhlWG++B7ERERmjp1qjZv3qzNmzfruuuu0y233GL9T5h7ULo2bdqkl19+WW3atHHbz30oHZdffrkOHz5svZKTk61jXrsHpoy68sorzQMPPOC2r0WLFmbs2LF+qlHFIsl88MEH1nZ+fr4JDw83U6dOtfadPn3ahIaGmvnz5/uhhhXD0aNHjSSzZs0aYwz3wZ8uuugi8+qrr3IPSll6erpp1qyZ+eqrr0x8fLx58MEHjTH8XSgtTz31lImJiSn0mDfvQZlsWcnOztb333+v7t27u+3v3r27vvvuOz/VqmLbu3evjhw54nZPgoKCFB8fzz3xodTUVEnSxRdfLIn74A95eXlasmSJMjIy1LlzZ+5BKRs5cqR69+6t66+/3m0/96H0/PLLL6pfv74aN26sAQMGKCUlRZJ374HPH2ToC8ePH1deXp4uueQSt/2XXHKJjhw54qdaVWwFv/fC7sn+/fv9UaVyzxijRx55RFdddZVat24tiftQmpKTk9W5c2edPn1aNWrU0AcffKBWrVpZ/xPmHvjekiVLtGXLFm3atMnjGH8XSkfHjh3173//W82bN9fvv/+uKVOmKC4uTj/++KNX70GZDCsFXC6X27YxxmMfShf3pPSMGjVK27dv19q1az2OcR9877LLLtO2bdt04sQJJSUl6d5779WaNWus49wD3zpw4IAefPBBffnll6patWqR5bgPvtWzZ0/r5+joaHXu3FlNmzbVG2+8oU6dOknyzj0ok91AYWFhqly5skcrytGjRz0SHEpHwehv7knpGD16tJYtW6ZVq1YpIiLC2s99KD1VqlTRpZdeqvbt2ysxMVExMTGaPXs296CUfP/99zp69KjatWungIAABQQEaM2aNXrhhRcUEBBg/a65D6WrevXqio6O1i+//OLVvwtlMqxUqVJF7dq101dffeW2/6uvvlJcXJyfalWxNW7cWOHh4W73JDs7W2vWrOGeeJExRqNGjdLSpUu1cuVKNW7c2O0498F/jDHKysriHpSSbt26KTk5Wdu2bbNe7du311133aVt27apSZMm3Ac/yMrK0s6dO1WvXj3v/l0oweBfR1iyZIkJDAw0r732mtmxY4d56KGHTPXq1c2+ffv8XbVyKz093WzdutVs3brVSDIzZ840W7duNfv37zfGGDN16lQTGhpqli5dapKTk83AgQNNvXr1TFpamp9rXn4MHz7chIaGmtWrV5vDhw9br8zMTKsM98H3xo0bZ7755huzd+9es337djN+/HhTqVIl8+WXXxpjuAf+cvZsIGO4D6Xh0UcfNatXrzYpKSlm/fr15qabbjI1a9a0vou9dQ/KbFgxxpi5c+eaRo0amSpVqpgrrrjCmr4J31i1apWR5PG69957jTFnpqk99dRTJjw83AQFBZlrrrnGJCcn+7fS5Uxhv39JZuHChVYZ7oPvDR482Pp/T506dUy3bt2soGIM98Bfzg0r3Affu+OOO0y9evVMYGCgqV+/vunXr5/58ccfrePeugcuY4zxQssPAACAT5TJMSsAAKDiIKwAAABHI6wAAABHI6wAAABHI6wAAABHI6wAAABHI6wAAABHI6wA8Lv77rtPffv2LfL4pEmTFBsbW2r1AeAshBUAjvfYY49pxYoV/q4GAD8J8HcFAJRd2dnZqlKlis+vU6NGDdWoUcPn1wHgTLSsALDt2muv1ahRo/TII48oLCxMN9xwg2bOnKno6GhVr15dkZGRGjFihE6ePGm9Z9GiRapVq5a++OILtWzZUjVq1NCNN96ow4cPF3md77//XnXr1tWzzz4rybMbqKDb6Pnnn1e9evVUu3ZtjRw5Ujk5OVaZw4cPq3fv3qpWrZoaN26st99+W1FRUZo1a5bXfy8AfIuwAqBY3njjDQUEBOjbb7/VggULVKlSJb3wwgv63//+pzfeeEMrV67UmDFj3N6TmZmp559/Xm+++aa++eYb/frrr3rssccKPf/q1avVrVs3TZ48WRMmTCiyHqtWrdKePXu0atUqvfHGG1q0aJEWLVpkHb/nnnt06NAhrV69WklJSXr55Zd19OhRr/wOAJQuuoEAFMull16q6dOnW9stWrSwfm7cuLGeeeYZDR8+XPPmzbP25+TkaP78+WratKkkadSoUXr66ac9zv3RRx/p7rvv1oIFCzRw4MDz1uOiiy7Siy++qMqVK6tFixbq3bu3VqxYoYSEBP3000/6+uuvtWnTJrVv316S9Oqrr6pZs2Z/6bMD8A/CCoBiKfjyL7Bq1So999xz2rFjh9LS0pSbm6vTp08rIyND1atXlyQFBwdbQUWS6tWr59HKsWHDBn3yySd67733dOutt16wHpdffrkqV67sds7k5GRJ0s8//6yAgABdccUV1vFLL71UF110UfE/MAC/oxsIQLEUBBBJ2r9/v3r16qXWrVsrKSlJ33//vebOnStJbuNHAgMD3c7hcrlkjHHb17RpU7Vo0UKvv/66srOzL1iPws6Zn58vSR7nLlDUfgDORlgBUGKbN29Wbm6u/vWvf6lTp05q3ry5Dh06VKJzhYWFaeXKldqzZ4/uuOMOt7BTXC1atFBubq62bt1q7du9e7dOnDhR4nMC8B/CCoASa9q0qXJzczVnzhylpKTozTff1Pz580t8vrp162rlypX66aefNHDgQOXm5pboPC1atND111+voUOHauPGjdq6dauGDh2qatWqyeVylbh+APyDsAKgxGJjYzVz5kxNmzZNrVu31ltvvaXExMS/dM7w8HCtXLlSycnJuuuuu5SXl1ei8/z73//WJZdcomuuuUa33nqrEhISVLNmTVWtWvUv1Q9A6XMZOnEBVAC//fabIiMj9fXXX6tbt27+rg6AYiCsACiXVq5cqZMnTyo6OlqHDx/WmDFjdPDgQe3atctjcC4AZ2PqMoByKScnR+PHj1dKSopq1qypuLg4vfXWWwQVoAyiZQUAADgaA2wBAICjEVYAAICjEVYAAICjEVYAAICjEVYAAICjEVYAAICjEVYAAICjEVYAAICjEVYAAICj/X9y2cxs7NpTGwAAAABJRU5ErkJggg==",
"text/plain": [
"
"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sc.pl.umap(adata, color=[\"IGHM\", \"JCHAIN\"])"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Save AnnData\n",
"\n",
"We can save this `AnnData` object for now."
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"adata.write(\"adata2.h5ad\", compression=\"gzip\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Save DandelionPolars"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"To save the vdj object, we have two options - either save the `.data` and `.metadata` slots with pandas' functions.\n",
"\n",
"\n",
"\n",
"
\n",
"\n",
"Polars backend update\n",
"\n",
"Since `vdj.data` is a Polars DataFrame in the polars backend, the pandas `.to_csv()` method is no longer available.\n",
"Use Polars' `.write_csv()` instead (note the parameter is `separator=`, not `sep=`).\n",
"\n",
"
\n"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"vdj.data.collect().write_csv(\"filtered_vdj_table2.tsv\", separator=\"\\t\")"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"vdj.write_h5ddl(\"dandelion_results2.h5ddl\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Concatenating multiple bcr objects\n",
"\n",
"It is quite common that one might be trying to analyse data from multiple samples. In that case, `dandelion` has a `concat` function to merge the data.\n",
"\n",
"We will simulate a second object but reading in the same file."
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [],
"source": [
"vdj1 = ddl.read_10x_vdj(folder_location, filename_prefix=\"filtered\")\n",
"vdj2 = ddl.read_10x_vdj(folder_location, filename_prefix=\"filtered\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Before you merge the objects, make sure that the \"cell_id\" and \"sequence_id\" are distinct so that you can distinguish them later"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [],
"source": [
"vdj1.add_cell_prefix(\"run1_\")\n",
"vdj2.add_cell_prefix(\"run2_\")"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Lazy Dandelion object with n_obs = 1988 and n_contigs = 5202\n",
" data: cell_id, is_cell_10x, sequence_id, high_confidence_10x, sequence_length_10x, locus, v_call, d_call, j_call, c_call, complete_vdj, productive, junction_aa, junction, consensus_count, umi_count, clone_id, raw_consensus_id_10x, sequence, rearrangement_status\n",
" metadata: cell_id, clone_id, clone_id_rank, locus_VDJ, locus_VJ, v_call_VDJ, v_call_VJ, d_call_VDJ, j_call_VDJ, j_call_VJ, c_call_VDJ, c_call_VJ, junction_aa_VDJ, junction_aa_VJ, junction_VDJ, junction_VJ, umi_count_VDJ, umi_count_VJ, locus_VDJ_main, locus_VJ_main, v_call_VDJ_main, v_call_VJ_main, d_call_VDJ_main, j_call_VDJ_main, j_call_VJ_main, c_call_VDJ_main, c_call_VJ_main, junction_aa_VDJ_main, junction_aa_VJ_main, junction_VDJ_main, junction_VJ_main, umi_count_VDJ_main, umi_count_VJ_main, isotype, isotype_main, isotype_status, locus_status, chain_status, rearrangement_status_VDJ, rearrangement_status_VJ"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"vdj_merged = ddl.tl.concat([vdj1, vdj2])\n",
"vdj_merged"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"