{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "a540c86a-3bed-4543-b952-342c796d7507", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import seaborn as sns\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "dc1ec4fb-4073-4d32-887f-e488ef0a9437", "metadata": {}, "outputs": [], "source": [ "t = np.arange(0,45)" ] }, { "cell_type": "code", "execution_count": 3, "id": "74ff2c1e-6c77-4aef-a26f-e86e5b1ce2d8", "metadata": {}, "outputs": [], "source": [ "# create signal\n", "sig = -3*np.ones((6,45))\n", "for rownum, start in enumerate([15,12,16,2,14,100]):\n", " sig[rownum,start:] = 0" ] }, { "cell_type": "code", "execution_count": 4, "id": "fdc289c2-b05f-4851-9fb2-01a7bda5dd5f", "metadata": {}, "outputs": [], "source": [ "# create measurement\n", "meas = np.ones((6,45))\n", "data2 = {'well':[], 'cycle': [], 'fluorescence': []}\n", "data = {}\n", "for row in range(6):\n", " value = sig[row,0]\n", " for col in range(45):\n", " diff = sig[row,col] - value\n", " value += 0.4*diff\n", " meas[row,col] = value\n", " data2['well'].append('well {}'.format(row+1))\n", " data2['cycle'].append(col)\n", " data2['fluorescence'].append(value)\n", " data[f'well {row}'] = meas[row]\n", "data[f'cycle'] = t\n", "df = pd.DataFrame(data2)" ] }, { "cell_type": "code", "execution_count": 5, "id": "6e5b1731-6d98-41dd-8a38-f8c168aeaace", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGwCAYAAACq12GxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB/OElEQVR4nO3deXxU9dX48c+dNfsCCVkgEFYhsi8KWkhwA7S4tValj0pbKbYgIFXR+qhYF1yoRaVF7Y8Hpa5tKVZaF2g1UQSULYCIbAYIhEAC2Ulmvb8/JjMhJpmZJDfMHXLeL+eVzOTmzjdMTE7OPd9zFFVVVYQQQgghOilDqBcghBBCCBFKEgwJIYQQolOTYEgIIYQQnZoEQ0IIIYTo1CQYEkIIIUSnJsGQEEIIITo1CYaEEEII0amZQr0AvXO73RQVFREbG4uiKKFejhBCCCGCoKoqVVVVpKenYzD4z/1IMBRAUVERGRkZoV6GEEIIIdqgsLCQHj16+D1GgqEAYmNjAc8/ZlxcXIhXI4QQQohgVFZWkpGR4fs97o8EQwF4L43FxcVJMCSEEEKEmWBKXKSAWgghhBCdmgRDQgghhOjUJBgSQgghRKcmwZAQQgghOjUJhoQQQgjRqUkwJIQQQohOTYIhIYQQQnRqEgwJIYQQolOTYEgIIYQQnZoEQ0IIIYTo1MIuGPrTn/5E7969iYiIYNSoUXz++ed+j8/Ly2PUqFFERETQp08fXn755XO0UiGEEEKEg7AKht59913mzZvHQw89xPbt2xk/fjxTpkzhyJEjzR5fUFDA1Vdfzfjx49m+fTu//e1vmTNnDqtWrTrHKxdCCCGEXimqqqqhXkSwLr74YkaOHMmyZct8jw0aNIjrr7+eRYsWNTl+wYIFvP/+++zZs8f32F133cWOHTvYuHFjUM9ZWVlJfHw8FRUVMqhVdEqqqoL35naDqjZ+rP7m+UlSf9/ziQ0373kaTtr82++/3+Tzml1hc4sO7rHW0OOPykb/xm4a/v2/9zp8/3HVXf+xRier/7eufy3dKqr3fd+pVM8pANXtxq2CW1Vxq56Pqb7XuuF8vmPcKm5UVLeKW/Ws1V1/MrV+Pb5zeZ+3fl3u761XrX9Hrf+c7/2jNPrcs7lb+jds8lAzjzX3qc2dr1kh/N4J8qmDDQU66itJ6pFBvyEjND1na35/h83UervdztatW3nggQcaPX7VVVexYcOGZj9n48aNXHXVVY0emzRpEsuXL8fhcGA2m5t8js1mw2az+e5XVlZqsHohPFRVRa2txX3mDO6aGs9b763G81a11eG22VDtDlS7HdVm87x12HHb7ag2e/19B6rLCQ4nqsvV+H2nE1xOVKcL1eUClwvV7Qa3G9XtApf3/cZvcbsbAp36wEeXQUAYUFFwGS24jNZmb06jFbfBjNtgwm0woSoN73tu9fcVE6rBhKoYcCtGVMWAevZbg/f9+vsoqIoBFAUVw/feV0AxoKLUP6YAiu9x/VC+976e1iY6QqRjNf2WaxsMtUbYBEOlpaW4XC5SUlIaPZ6SkkJxcXGzn1NcXNzs8U6nk9LSUtLS0pp8zqJFi3jssce0W7g4b6lOJ87SUlxlZThPn8ZVVo6rrMxzv+ys+6dP4yov9wQ/tbUSXLREUZp/P9DHmntMVRse8v1716c0aOl+y9wYsFvicFhisJtjcNTf7OYYHKaz3jdH4zBF4zJG4DJadBZg6IjabI7G+0Hfe0qwmcAOyFco8v/pudX8i33OhE0w5KV874eeqqpNHgt0fHOPez344IPMnz/fd7+yspKMjIy2LleEOVdFBfbCoziOFmIvLMThff/oMRxFReB0tu3EioIhKsp3U6Lr34+MwhARgWKxoFitKBYzisWCwWpFMVsaP242o5jMKCYjGI2+9xWTCYymhvcNRhSj4ay3nptiNIJiaHjM+743S6CAYvB+TPH8P+O9GQyA4gk4zn68/qac9XX6bt779W/9/X8b+IVxQEUhlB0663a44f268ladTlXBpsZQ6Uqhwt2DSnpQ6U6n0pVMpT2JKns8ahuzEyoqLsWB0+DAbrThMNqwmc7gNNpwGhy4DA5cihOXof7me7/hcbfiqr+5cSsuUMBgMGJQjBiNJowGEyajuf6tEaPRiNFgwlz/mMVoxmTy3DebzJiNZkxGMyaDAZPBhMlkxmgwYjGaMBlNmExGzAYzRsWAyWTCZDBiNBgwKob651M8jylGTAYDBqMRgwJGxYjBYMCoKBgUBYNirP92M2JUDBjqX3cDCgaDAQXP94pSf7znredxz3+e7xGFhu8X32Pt+f4ROjQ5pM8eNsFQUlISRqOxSRbo5MmTTbI/Xqmpqc0ebzKZ6Nq1a7OfY7VasVqt2ixahA3V4aBu7z5qd+RTu2MHtgMHcBQexV1V5f8TTSaMiQmYErtgTEzEmJiIqUsixoTE791PwBAT0xD8REbKD/NgVZ2Ao5vh6FdQlA9lBVBxNEB2AbDGQ1QXiOoK0Umet1FdcFmTOFWXSnFZIsUlEZSdMlBZ4cJe5/98bqDOoFKDyhmDSq2iUqvgeWusoy7yJLaIIuxRR3FHnMRuKcFprMNpcDS+6nOWaFMM8dYEEiMS6RqRSHJEIokRXUmMSCTRmkhC/cfiLHFEmaOIMkcRaYrEbGh6iV8I0XZhEwxZLBZGjRrFunXruOGGG3yPr1u3juuuu67Zzxk3bhxr1qxp9NjatWsZPXp0s/VCovNwnDxJbb4n8KndsYO6r3ej1tU1e6wxOQlLjwzMGT3q32ZgyeiBOSMDU3KyJ3sitOG0Q/HO+uBnMxRuhormd4tiioDEzMa3hF717/cCSzQA1WU2ThRUUFxQyYkdFZw8UoXL4Q18XPW3+lNGm3BEGDiluDlcZ+eE20m5QaXC4KbaU14DipO4uJMkdjmOIfIwtUoBNe6SZpdoNVrpFdOXHrE9PLeYhrfpMelEmaO0+FcTQrRT2ARDAPPnz+e2225j9OjRjBs3jldffZUjR45w1113AZ5LXMeOHWPlypWAZ+fY0qVLmT9/PjNmzGDjxo0sX76ct99+O5RfhggBZ2kplR9/zJktW6jdsQNn0fEmxxji4ogcNozIYcOIyBqEJSMDc/fuGKLkF1aHcbvg0Oewb60n+Dm+A1y27x2kQLcs6DHac0sa4Al4YlKarR+qqbBx8IuTFO0v4ERBBdVl3z8fWKNMJGXGUhVt4KDdzq6KGnaW1+CAhtjIBFaTgcHdYxmZegyH9RtO2vdxuGofdred01C/rcpz6aZfYj+GJQ9jSNIQesf3pkdMD5IikyQDKEQYCKtg6Oabb+bUqVP87ne/4/jx4wwePJgPPviAXr16AXD8+PFGPYd69+7NBx98wD333MMf//hH0tPTefHFF/nRj34Uqi9BnENum43qTz6h4r1/Ur1+PbgaMgAYDFj79/cEP8OHEzl8GJbMTMnynCsle2HH27Dzr1B5rPHHIrtAjzGQMcbzNn0kRPjfFuu0uyjYUcq3m4op/OZUoxp1RYGuPWJI6R1Pt8xYis0qHxSU8OHXxdQUuRqdp2eXKEb0TGBERgLdupazpyqXDw/9m32VJxodl2hNZGjyUIYmD2VY8jAGJw0m2hzdnn8RIUQIhVWfoVCQPkPhRVVVardvp+K9f1L54YeNan4ihg4l9rKJRA4fTsTgIRhj5JfXOVVzCr5eBTvegqLtDY9HxMOgayHzB57gp0uf5neMfY/qVjl+sJxvNxVzcOtJ7HUNgU1qnzgyhyaR2jue5F6xFFbVsXrbMVZvP8ax8lrfcRldIrl6SBqje3VheEYCRnMNHxZ8yJqDa9h9arfvuFhLLFf1uopRKaMYljyMjNgMyfgIoXPnZZ8hIfyxFxZS8c/3qfjnP3EUFvoeN6WnEX/ttcRfex3WPr1DuMJOymmDfR/Djndg/8fgrt99ZzBBvyth2C0wYDKYI4I+ZfmJM+z9spi9XxZTdaqhziu2SwQXjE3lgotTSUiJoqzGzr92FrHqP9+QX1jecFyEiR8OTedHI7szqlciDreDvKN5/G7z+6w/uh6n6lmjSTHxgx4/4Nq+15LdIxuL0aLJP4kQQn8kGBJhrWbTl5QsfYnaLVt9jxmiooidNIn4664j6qIxcukrFOoq4bNnYdtfGm9zTxsOw26FwT+CmORWnbJwz2m+WlNA8XcVvsfMEUb6jerGwLGppPVNQDEoFJTW8Ns3t7LumxM4XJ7Et9GgkDMgmRtH9uDyQd2IMBupslfx+y2/Z/WB1VTaG5qrXtj1Qqb2ncqU3lPoEtGlPf8KQogwIcGQCEvO0lJOPPssle/X7xY0GIi+5BLir7uW2Msvl6LnUFFVz6Wwjx+C6vq2FrHpMPQnnixQt0GtPmVNhY0v/raf/VtOAp4raBlZXRk4NpXMYUmYLUYA6hwulv33IMvyDmJ3enaLDe4ex40jejB1WDrJsdb6Jap88N0HPLflOUprSwFIiUrhh31+yNS+U+mb0Le9/wpCiDAjwZAIK6rLRdm771LyhyWeeiBFIfHWW+g6cybmFvpNiXOk9AB88Bv4Ltdzv0tfmPQU9L8SDMZWn87tcrMr7xhfvv8djjoXigKDc3owanIvouMb9wLL21fCI//8msOnzgAwYUAyD04ZyKC0xnUCBRUFPPnlk3x5/EsAMuMyuW/MfVyafinGNqxRCHF+kGBIhI3aXV9T/Nhj1H39NQARF15I6sJHiRwyJMQr6+QctfD58/DFEnDZwWiFCffCpXPB1LYGpicKKsl961tKC6sB6JYZR860C0juGdvouOKKOn73r918sMuThUqJs/Lo1AuZMji1UYFznbOOP+/6Myu+XoHD7cBqtDJjyAx+NvhnUgskhJBgSOifq7KSkiUvUPb226CqGGJiSL5nHom33OIZKSFCZ/86+OBezwgMgH5XwNXPeXaEtUFdjYNN7x1k9/oiUD39gMZe35esH6RjMDQEN06Xm9c2HOIP6/ZRY3dhNChMvySTe64cQIy18Y+1z45+xlNfPsWxas8W/vHdx/PgxQ+SEStjdoQQHhIMCd1SVZXKf/2LE888i6vUU9sRN3UqKfffhym5dcW3QmMVx+CjB2DP+577sekw5WnPFvk2bDlXVZW9m4rZ8I8D1FY5ABg4NpVxN/YjKq5x5mbr4dM8tPprvi32tE0Y2TOBJ64fQlZ640tixTXFPP3V0/z3yH8BT13QAxc9wOU9L5dt8UKIRiQYErpkP3KE4488yplNmwCw9O5N6qOPED12bIhX1smpKmz6E3zyJDhqQDHC2F9BzgNgjQ38+c0oP3GGT9/4lqL95QAkpkWTM20A6f0TGx1XVefgyX/v4Z3NntYJCVFmHpg8kJ+MzmiUNXKrbv7yzV/4Y/4fqXXWYlSM3JZ1G78a9isZfyGEaJYEQ0J3bN8VcPj223GVlqJYrST96ld0+fnPMFiktiOkVBU+uA82/9lzP+NiuOZ5SB3c5lOeKqrmvee3U1ftwGQxMOaa3gy7PAOjqXE7hKo6B7ct/8rXL+gno3vwwJRBdIlu/D2hqipPffkU7+59F4AR3Ubwv2P/lwGJA9q8RiHE+U+CIaEr9sOHOTJ9Oq7SUqwDB9Jj6UtYevQI9bKEqsLHv60PhBSY8gyMmQHt6OFUVlzDP5fkU1ftILlnLJNnDiaua2ST42psTn62YjP5heUkRJl55X9GcXGfrs0sUeWZzc/w7t53UVB44KIHuGXgLRgU6TMlhPBPgiGhG/bCQg7fMR3nyZNY+/en54r/w5SYGPgTRcdSVVj3iOfyGMC1L8LI29t1yvITZ3jvD9uprbSTlBHDtXOHExFtbnLcGbuTn722mS2Hy4iLMPHGLy5mcPf4Zpao8vstv+fNPW8C8Nglj3FD/xvatUYhROchfzIJXXAcO8bhO+7AWVyMpW9fer62QgIhPVBV+ORx2PCi5/41z7c7EKoo8QRCZyrsdO0e3WIgVGt3cefrW/iq4DSxVhN/8RMIvbDtBV7/5nUAHhn3iARCQohWkWBIhJzj+HFPRqjoOJbMTE9GqGvTyyAiBHKfhs9/73l/ynMw5hftOl1laS3v/WE7NeU2EtOiuXbuCCJjmtaC1Tlc/PIvW9hw8BQxVhOv/+IihmUkNHvOP+b/keVfLwfgtxf/lpsG3NSuNQohOh8JhkRIOU6c4PD06TiOHsXcsyc9X38Nc7duoV6WAPjsOch72vP+pKfg4l+263RVp+t47w/bqT5tIyEliuvmDW+ybR7A5nRx1xtb+Xx/KVEWI6/9bAwjezafJXx5x8u8svMVABaMWcCtA29t1xqFEJ2TBEMiZBwnT3Lkjuk4Dh/B3KMHvV5/TUZq6MX6JfDJE573r3gMxs1q1+mqy2y894ftVJ2qIz45kuvvGdFkpAaA3enm129sI3dvCRFmA/83fQyjM5sflvr/dv0//pj/RwB+M+o3/E/W/7RrjUKIzkuCIRESztJSjvzs59gPHcKUnkbP117DnJYW6mUJgI1/hP886nn/sv+FH8xr1+lqKmy894dtVJbUEpcUwXX3jCA6oWkg5HC5ufvtbfz325NYTQb+744xjG1m1xjA67tf54VtLwAwd+Rcpg+e3q41CiE6NwmGxDnnLCvzBEIHD2JKTaXX669j6dE91MsSAF++4tlCD5D9AEy4r12nO1Np559/2E7FyVpiu3gCodguEU2Oc7rczHsnn493n8BiMvDn20dzSb+kZs/5xjdvsHjLYgBmDZ/FnUPubNcahRBCgiFxTrnKyznys59j278fU3IyvV5bgSVDZkTpwubl8OH9nvfH/8bTVbodaqvs/HPJdsqKzxCTaOX6+SOa7SPkcqvM/+sO/r3rOBajgVf+ZxQTBjQ/buWdb9/hmc3PADBz6EzuGnZXu9YohBAgwZA4h1SnkyMzZ2L79luMSUn0fP11LJmZoV6WAPjmffj3fM/7l8yByx5u04wxL5fLzZqXdnC6qIboeAvX3TOCuKSmgRDA4//6hvd3FGEyKPzppyOZOLD5Avq1h9by5JdPAvCLwb9g1vD21TEJIYSXBEPinCn761+p27ETQ1wcvV5bgbVP71AvSQDYqhsyQmNmwJW/a1cgBLDr06OUHKkiItrMdfeMIKFb8zPBdh4t5/WNhwB46dYRXJHVfAF9jaOGp7/y7Gz76aCfMnfkXBm2KoTQjARD4pxwlpVR8oKncV/yvLlY+/UL8YqEz+eLoeo4JPSCq55odyB0ptLO5n8VADDuxr4kpkY3e5zbrfLo+7tRVbhhRHemDGm5gP6Vna9QUltCRmwG80fNl0BICKEpCYbEOVGy5AXcFRVYL7iAxJ/8JNTLEV6nDsKGpZ73Jy8Cc9Pi5tbauPoA9joX3XrFMmhcywHOP7YfY/uRcqItRh6YMrDF4woqCvjLN38BPL2ELEYZ2CuE0JYEQ6LD1X3zDeV//SsAqf/7EIpJRuLpxkcPgtsBfS+HC65u9+mKCyr4dmMxAONvHoBiaD6DU1Xn4OkPvwXg7sv7kxLXfBDmHb7qdDsZ33082RnZ7V6jEEJ8nwRDokOpqkrxE0+CqhJ39dVEjRkT6iUJr30fw/6PwWCCyU+3+/KY6lb5/J19AAwcm0pqn6ZzxLxe/O9+Sqtt9E6K5meXZrZ4XN7RPL449gUmg4n7x9zfrvUJIURLJBgSHapyzRpqt21DiYyk2/3t61kjNOS0wUf1W+fH/gqSB7T7lHs2Hufk4SrMEUbG3tC3xeMOnKxixReHAHhkahZWk7HZ42wuG8985dlGf3vW7WTGZ7Z7jUII0RwJhkSHcVXXcPI5T3O8pLvuwpyaGuIVCZ+Nf4TT30FMCkxof8bFdsbBpvcOAjDmmt7NjtoAT6bwsTXf4HSrXDGoGxMvaHkO3crdKzlafZTkyGR+ObR9c9GEEMIfCYZEhzn18jKcJSWYe/aky8+mh3o5wquyCD7zBKlc8RhExLX7lJv/dYjaKgeJqVEMndijxePWfnOCz/eXYjEa+N9rslo8rrimmD/v+jMA80fPJ9rc/I40IYTQggRDokPYCgo49fpKAFIefACDRXYA6ca6R8BRAz0ugqE3t/t0p4qq2Zl7FIAf/KQ/RlPzP1bqHC4e/9c3AMyY0JvMpJYDnOe3PE+ts5YR3UZwTe9r2r1GIYTwR4IhoTlVVTmxaBE4HERPGE9MTk6olyS8Dm+EXX8DFLj6WTC070eAqqqs/+t+VLdK72FJ9MxqfrAqwKuffcfRslpS4yL4dU7Lfaa2FG/hw0MfoqDw4EUPSk8hIUSHk2BIaK7601xqPvsczGZSHpRfZrrhdsGH9UXsI2+H9BHtPuV320s4+m0ZRpOBH9zUv8Xjjpad4U+5BwD47TWDiLY2317B6Xay6KtFANw04CYGdR3U7jUKIUQgEgwJTbltNk9WCOg6/Q6svWXkhm5sXQHFuyAiHi5/pN2nc9hdrP/7fgBGXNWzxdljAIs++JY6h5uLendh6tCWGzH+bd/f2Fe2jzhLHHePuLvdaxRCiGBIMCQ0dXrFazgKCzElJ9N1pkwU140zp+GTJzzvT3wIopPafcrtHx+m+rSNmEQrIyf3avG4DQdK+feu4xgUWDj1whYzhWV1ZSzd7umGffeIu0mISGj3GoUQIhgSDAnNOIqLKX3lFQC63X8fxhjZAaQbnzwBtWXQ7UIY/Yt2n66ytJZta48AcOmP+2O2NN8ryOlys3DNbgD+Z2wvstJb3rn20vaXqLRXMiBxAD8e8ON2r1EIIYIlwZDQzMlnn0OtrSVy5EjifvjDUC9HeB3f6blEBjDlGTC2fxzKF6sO4HK46X5BAn1HJrd43F82HWbfiWoSo8zMv7Llxo57Tu3h7/v+DsCDFz2IySAjW4QQ544EQ0ITNV99ReUHH4CieOaPSdG0PqgqfLgAVDdceAP0Ht/uUxbuOc1320tQDArjfzKgxde6tNrG8+s84znunXQBCVHNt1dQVZVFXy1CRWVK5hRGp45u9xqFEKI1JBgS7aY6nZx48ikAEm7+CRFZLTfTE+fY16vgyAYwRcJVT7T7dC6Xm8/f9QQ4Q7K707V7TIvHLv54L1V1Ti5Mj+OWMT1bPO7fBf9m+8ntRJoimT96frvXKIQQrSXBkGi36s8+w7Z3L4b4eJLnzg31coSXqjZ0mh7/G4hvuTN0sArySykrPkNEjJkxP2x5p+DJqjr+uqUQgMeuvRBjC9PrVVXlzzs9naZnDJlBarSMbBFCnHsSDIl2q/zXvwBIuOEGTImJIV6N8DmxG0r2gNECF83Q5JT7vioGIOvSdCKizS0e968dx3GrMLJnAqMzu7R43N6yvXxX8R0Wg4VbBt6iyRqFEKK1JBgS7eKqrqHqk08BpGhab3b91fO2/1UQmdDu09XVODi8+xQAAy5K8XvsP/OPAXDd8O5+j/v3d/8GIDsjm1hLbLvXKIQQbSHBkGiX6v/+B7WuDktmJhEXSq2QbrjdsGuV5/0hN2lyyu+2l+B2qnRJj/ZbK1RQWsOOoxUYDQrX+Gmw6FbdfFDwAQBX975akzUKIURbSDAk2qXiX56/7OOm/lB2kOlJ4SaoPArWOBgwSZNT7tvsuUQWKCv0fn4RAD/ol0RSjLXF47ae2MrJMyeJNccyvkf7d7kJIURbSTAk2sx56hQ1GzYAEH+NTBbXlZ31l8gGTQVzy2MyglVdZuPYvnIA+o9uORhSVfWsS2Tpfs/pvUR2Ra8rsBpbDpqEEKKjSTAk2qzyw4/A5SJiyBAsmZmhXo7wctrhm/c87w/RppPz/i0nQIW0vvF+Z5B9fayS70priDAbuOrClneG2V121h1eB8A1fSSQFkKElgRDos0q16wBIH6qFE7rysFPPKM3ortB72xNTrl/8wkg8CWy9+qzQlcMSiGmhcn0AF8c+4JKeyXJkcmMTpEmi0KI0JJgSLSJvbCQ2h07wGAgbsqUUC9HnM27i2zwj8DQ/Myw1igrrqHkSBUGg0LfUd1aPM7lVlmzw1MvFHAXWYHnEtnk3pMxarBGIYRoj7AJhsrKyrjtttuIj48nPj6e2267jfLycr+fM336dBRFaXQbO3bsuVnwea7y355fZtFjx2JKbnk2lTjHbNXwrWeHFkO12UW27ytPVijjwi5ExjQ/UgPgy+9OcbLKRnykmewBLX9P1DhqyC3MBeQSmRBCH8JmGuK0adM4evQoH330EQC//OUvue2221hTf6mmJZMnT2bFihW++xZLyz/MRXBUVaVijafRovQW0pm9H4CzFrr0gfSR7T6dqqrs814iGxPcJbKrh6RhMbX8d9YnRz7B5rKRGZdJVhdpxyCECL2wCIb27NnDRx99xKZNm7j44osB+POf/8y4cePYu3cvF1xwQYufa7VaSU2VFv9asn37LfaDB1EsFmKvvCLUyxFn8+4iG3ITaNDq4MShSipLajFZDGQOTWrxuDqHiw+/9my9D3YX2dV9rpZ2DEIIXQiLy2QbN24kPj7eFwgBjB07lvj4eDbUb+1uSW5uLt26dWPAgAHMmDGDkydP+j3eZrNRWVnZ6CYa82aFYiZOxBgrXYN1o6bUUzwNmjVa3F9/iaz3sGQsES3/7ZS79yRVdU7S4iO4yM/4jdLaUjYd3wRIo0UhhH6ERTBUXFxMt25NCze7detGcXFxi583ZcoU3nzzTT755BN+//vfs3nzZi677DJsNluLn7No0SJfXVJ8fDwZGRmafA3nC9Xt9tULxf1Q6j10ZfdqUF2QNhyS+rf7dG6Xm/1bPX88BB6/4SmcvnZYOoYWhrICrD20FpfqYnDXwfSK69XuNQohhBZCGgwtXLiwSYHz929btmwBaDadrqqq3zT7zTffzDXXXMPgwYOZOnUqH374Ifv27ePf9b/Mm/Pggw9SUVHhuxUWFrb/Cz2PnNmyBeeJExhiY4mZMCHUyxFn2/U3z1uNskJH95ZRW2knItpMRlbL2Z7KOgf//dYTNF0b6BJZ/S4yKZwWQuhJSGuGZs+ezS23+J9UnZmZyc6dOzlx4kSTj5WUlJCS4v8v1rOlpaXRq1cv9u/f3+IxVqsVq1W64baksv4SWeykqzDIv5N+lB2Cwi8BxbOlXgPeS2T9RnXDaGz576aPvi7G7nTTv1sMWWlxLR5XWFXIzpKdGBQDk3tP1mSNQgihhZAGQ0lJSSQltVyU6TVu3DgqKir46quvuOiiiwD48ssvqaio4JJLLgn6+U6dOkVhYSFpaS0PjxQtc9vtVK5dC0C87CLTl11/97ztPR7i2v/97bS7OJhfAkD/IGeRXTc83W+m9oPvPFv+L069mKTIwP/fCyHEuRIWNUODBg1i8uTJzJgxg02bNrFp0yZmzJjBD3/4w0Y7yQYOHMjq1asBqK6u5t5772Xjxo0cOnSI3Nxcpk6dSlJSEjfccEOovpSwVrN+Pe6KCkzJyUSNGRPq5QgvVdX8EtmhXadw1LmI6WIlrU98i8edrKxjw8FSAK4d1nKjRVVVfZfIru4jhdNCCH0Ji2AI4M0332TIkCFcddVVXHXVVQwdOpS//OUvjY7Zu3cvFRUVABiNRnbt2sV1113HgAEDuOOOOxgwYAAbN24kVnZAtUlFfU+nuGuuQTFK12DdOLEbSr4FowUGXavJKfd9VT+hfkwqip+C6DU7j+NWYWTPBHp2jWrxuL1leymoKMBisHB5z8s1WaMQQmglLPoMAXTp0oU33njD7zGqqvrej4yM5OOPP+7oZXUaruoaqj/5FJBGi7rjHb/R/yqITGj36epqHBzefQoIvIvsfd+E+gDjN+p7C2VnZBNrkT9GhBD6EjaZIRFaVf9Zh2qzYendm4gLpWuwbrjdsGuV5/2hP9HklN9tL8HtVOnaPZqu3WNaPK6gtIYdRyswGhSuGdpynZJbdfNBgade6JresotMCKE/EgyJoFT+q6G3kHQN1pHCTVB5FKxxnsyQBvZt9lwi6x9g/MY/67NCP+iXRFJMyzsLt57YyskzJ4k1xzK+x3hN1iiEEFqSYEgE5CwtpWbjRkB2kemOd/zGoKlgjmz36arLbBzbVw74D4ZUVW20i8wf7yWyKzOvxGKU2YBCCP2RYEgEVPnhR+ByETF0KJZe0jVYN5x2+OY9z/tajd/YcgJUSOsXT1zXloOrXccq+K60hgizgasubHn2n91lZ93hdYCM3xBC6JcEQyKgyn95Gi3Gy/gNfTn4CdSWQXQ36K1NN/D9QU6o947fuGJQCjHWlvdhfHHsCyrtlXSL7MbolNGarFEIIbQmwZDwy37kCLU7doDBQNyUKaFejjibdxfZ4B+Bof2tDsqKayg5UoXBoNB3VNNZgF4ut8qaHd5LZAF2kdX3FprcezJGDdYohBAdQYIh4Zd3KGv02LGYkpNDvBrhY6uGbz07tBiqzSWyffXjNzIu7EJkTMu1PZu+O8XJKhvxkWayB7T8PVHjqCG3MBeQWWRCCH2TYEi0SFVVKupnkUlvIZ3Z+wE4a6FLH0gf2e7TqarKvqAvkXl2kV09JA2LqeUfIZ8c+QSby0ZmXCaDugxq9xqFEKKjSDAkWmTbswf7d9+hWCzEXnVlqJcjzubdRTbkJ6BBq4MThyqpLKnFZDHQe1jL2Z46h4sPd3m23l8f5C6yq/tcLe0YhBC6JsGQaFHVfz8BICY7G2NMy833xDlmr4HvPN3AtZpQX7DDM1+s99AkzNaWa3u+LDhNlc1JalwEYzK7tHjcGccZNh3fBMguMiGE/kkwJFpUu3MnAFFjLw7xSkQjx3eA2wkxqZDUX5NTnijwzPTrfkGi3+Pyj5QDMLZPFwx+ZpbtPrUbl+oiJSqFXnHSjkEIoW8SDIlmqapKXX0wFDl0WIhXIxo5ttXztsdoTS6Rud0qJw9XAZDSu+UJ9QA7jpYDMDwjwe9xX5d+DcDQ5KHtXp8QQnQ0CYZEsxyFhbgqKlDMZiIuGBDq5YizeYOh7u0vnAYoLz6Do86FyWKgS1rLk+dVVWVHYTkAwwIEQ7tKdwEwOGmwJmsUQoiOJMGQaFbtTs8vM2vWIBSLjFDQFV8wNEqT0504VAlAcs9YDMaWfyQcLavlVI0ds1FhUFqc33N6M0NDkoZoskYhhOhIEgyJZtXu3AFA5BC5zKEr1SVQfgRQIH2EJqc8WR8MpWT6D3Dy67NCg9LiiDC3XGRdWlvK8ZrjKChkdc3SZI1CCNGRJBgSzaqrzwxFDpNgSFe8WaGkARDhv74nWN7MUMB6ofpgKFC90K4Sz/dO34S+RJuj270+IYToaBIMiSZUh4O6b74BIHKIXObQFY0vkTntLk4drQagW2as32O9xdPDeiT4Pc5bLySXyIQQ4UKCIdFE3b59qHY7hrg4zDKlXl80Lp4uPVqN260SGWsmtktEi8c5XW52HfNsvw9UPO2tF5LiaSFEuJBgSDTh21I/ZIh0DtYTVdW+eLqgoV7I32u990QVdQ43sVYTfZJavvTlVt18fUqKp4UQ4UWCIdGEdydZxFD5ZaYrp7+DunIwWiFFm6yLt16oW4Di6R2FnqzQ0Ix4v80WD1cepspehdVopV9iP03WKIQQHU2CIdFE7S5vs0UpntYVb1YobSiYtGl30FA8HSgYKgeCb7Y4qMsgzAZzu9cnhBDnggRDohFXdTX2g98BUjytOxpfIqurdlBZUgtAt14BgqHWFk8ny/eOECJ8SDAkGqn7ejeoKub0dExJSaFejjib1vVChz1ZofhukUREt5zFqbE52XfCM64j2MyQ1AsJIcKJBEOiEe9w1gi5RKYvTjsc97w2WgVDwTZb3HWsArcKafERdItreceZ3WXn29PfArKTTAgRXiQYEo3USb2QPp34Glw2iEiALn20OWWBtvVCe0/vxeF2kGhNpEdMDy2WKIQQ54QEQ6IR706ySNlJpi9nXyLToN2BqqrB7yTz1gu1YjirtGQQQoQTCYaEj+PECZwnToDBQESWzJTSlWPbPG81ukRWdaqOumoHBqNCUo8Yv8d6t9UHKp6WeiEhRLiSYEj4eOuFrP37Y4iKCvFqRCMdNKk+qUcMJj9DV09W1XGsvBZFgSE9/M8uOzszJIQQ4USCIeHjG84q9UL6UlcBpfs872scDAUqnvZmhQZ0iyXGamrxuEp7JYcqDwESDAkhwk/LP91Ep1O7SzpP61LRdkCFhJ4Qk6zJKU/WF093C7J4eliG/6yQ9xJZj5geJEYktn+BQpwHVFXF6XTicrlCvZTzktFoxGQyaVKjKMGQAEB1u6nbJZkhXdL4EpnL5abkiKdvUMDMUJDF0756IWm2KAQAdrud48ePc+bMmVAv5bwWFRVFWloaFkv7uvJLMCQAsBcU4K6pQYmMxNq3b6iXI86mcfH06aIanA43lggjCd1arg1zu9WGzFCwnaeleFoI3G43BQUFGI1G0tPTsVgsssNSY6qqYrfbKSkpoaCggP79+2MwtL3yR4IhAUDtjvpmixdmoZjk20JXNM4MnTxrS73iZ+hqwakaKuucWE0GLkiNbfE4VVXZVSLBkBBedrsdt9tNRkYGUbIZpcNERkZiNps5fPgwdrudiIiWm8IGIgXUAjh7OOuwEK9ENFJZBFXHQTFCmjavTfDF0+UADOkej9nY8o+KE2dOcKruFEbFyMAuAzVZoxDng/ZkKkRwtPo3lldKAGfvJJO/7HXl6BbP225ZYInW5JTeztMBmy36iqcT/B63s8QTSA9IHECEqe1/mQkhRKhIMCRw22zU7d0LyKR63fFdIhupyensdU5OH68BAo/hyD9a32xRhrMKIc5zEgwJ6r75BpxOjF27YkpPD/VyxNk0rhcqOVIFKsQkWomOt7Z4nM3pYk+RJ4M0PMjiaekvJMT5KzMzkyVLlvjuK4rCe++9F7L1aE2CIdFoS73seNARtwuK8j3v9xitySmDrRfac7wKu8tNl2gLGV0iWzzO5Xax+9RuQDJDQojwJcGQkOGselW6H+xVYI6GZG0Kk08GO5zVt6U+3m+A/F3Fd9Q6a4kyRdE7vrcmaxRCiHNNgiHh20kWMUSaLerKsfri6fThYGh5flhreIung91J1ppJ9UaN1iiEaL81a9aQkJCA2+0GID8/H0VRuO+++3zHzJw5k1tvvRWADRs2MGHCBCIjI8nIyGDOnDnU1NSEZO2hIMFQJ+csK8Nx+AgAkUOk5kNXNC6erqmwUV1mAwWSe7XcNwggvw3BkBBCPyZMmEBVVRXbt28HIC8vj6SkJPLy8nzH5Obmkp2dza5du5g0aRI33ngjO3fu5N1332X9+vXMnj07VMs/5yQY6uTqvvbsBLJkZmKM9z9/SpxjvmBIm3oh7yWyLmnRWCJabqxZccbBd6WevwgDFU/LTjIh9Ck+Pp7hw4eTm5sLeAKfe+65hx07dlBVVUVxcTH79u0jJyeH5557jmnTpjFv3jz69+/PJZdcwosvvsjKlSupq6sL7Rdyjkgw1MnV7qy/RCb1QvriqIUTnsLkcz2pfuexcgB6dY0iMbrleT+1zlr2l+0HJDMkhB7l5OSQm5uLqqp8/vnnXHfddQwePJj169fz6aefkpKSwsCBA9m6dSuvvfYaMTExvtukSZN8Y0U6A5m70Mn5mi1KvZC+FO8CtxOiu0F8D01O2fri6QS/x317+ltcqovkyGRSolK0WKIQQkM5OTksX76cHTt2YDAYyMrKIjs7m7y8PMrKysjOzgY8s9RmzpzJnDlzmpyjZ8+e53rZISHBUCemqiq1u2QnmS55O093HwUatDtQ3SonDgU3qT6/MLhmi97O00OShkhLBiF0yFs3tGTJErKzs1EUhezsbBYtWkRZWRlz584FYOTIkezevZt+/fqFeMWhEzaXyZ588kkuueQSoqKiSEhICOpzVFVl4cKFpKenExkZSU5ODrt37+7YhYYRx7FjuE6fBrMZ66BBoV6OOJu3XqiHNpfIyk+ewV7rxGQ20KV7y2M9VFX1FU8PD7bzdLIE0kLokbdu6I033iAnJwfwBEjbtm3z1QsBLFiwgI0bNzJr1izy8/PZv38/77//PnfffXfoFn+OhU0wZLfbuemmm/jVr34V9Oc8++yzPP/88yxdupTNmzeTmprKlVdeSVVVVQeuNHzUeeuFBg7EYGm5NkSEQAdNqk/uGYvRz9DVooo6SqttmAwKF6b7zyDJTjIh9G/ixIm4XC5f4JOYmEhWVhbJyckMqv8jeOjQoeTl5bF//37Gjx/PiBEjePjhh0lLSwvhys+tsLlM9thjjwHw2muvBXW8qqosWbKEhx56iBtvvBGA119/nZSUFN566y1mzpzZUUsNG75mizKPTF/OnIay+qLF9BGanNJ7iSzYeqGBabFEmFvuG3S67jTHqo8BcGHXCzVZoxBCe4sXL2bx4sWNHsvPz29y3JgxY1i7dm2L5zl06FCj+6qqarE83QibzFBrFRQUUFxczFVXXeV7zGq1kp2dzYYNG1r8PJvNRmVlZaPb+cpbLyQ7yXTGmxXq2g8iEzU55YkCTx1Q0M0Wg9xS3zu+N7EW/z2LhBBC787bYKi4uBiAlJTGu1xSUlJ8H2vOokWLiI+P990yMjI6dJ2hojoc1NXXT0UOlZ1kuqLxJTKXw03p0WogcGZoeyubLUp/ISHE+SCkwdDChQtRFMXvbcuWLe16ju/vclFV1e/OlwcffJCKigrfrbCwsF3Pr1e2AwdQ6+owxMZiycwM9XLE2TRutlh6tBq3SyUixkxcUkSLxzldbnYd9WSQRkgwJIToREJaMzR79mxuueUWv8dktvEXdWpqKuDJEJ1dBHby5Mkm2aKzWa1WrFZrm54znDTUCw1GMZy3CcLwo6qaZ4bObrbo7w+BAyXV1DpcxFhN9EmO8bNEVTpPCyHOKyENhpKSkkhKSuqQc/fu3ZvU1FTWrVvHiBGeIlS73U5eXh7PPPNMhzxnOJHhrDpVfhjOnAKDGVK12aXV2maLQ7rHYzS0HDQdrTpKha0Cs8HMgMQBmqxRCCFCKWxSAkeOHCE/P58jR47gcrnIz88nPz+f6upq3zEDBw5k9erVgOfy2Lx583jqqadYvXo1X3/9NdOnTycqKopp06aF6svQDV/naSme1hdvs8XUIWDSJkMZ7BiOoJstlnoC6UFdBmE2mtu/QCGECLGw2Vr/yCOP8Prrr/vue7M9n376qa9/wt69e6moqPAdc//991NbW8uvf/1rysrKuPjii1m7di2xsZ1794u7pgbbgQMARMi2en05ts3ztoc29UJ1NQ7KT5wBggmGygFptiiE6HzCJhh67bXXAvYY+n7fA0VRWLhwIQsXLuy4hYWh2t27we3GlJaGuVu3UC9HnE3jeqGSw57+QnHJkUTEtJzFOWN3su+E59hAwZA0WxRCnG/C5jKZ0E7dLmm2qEsuBxzf4Xn/HE+q311UicutkhJnJTW+5R1nDreDPaf2AFI8LYQ4f0gw1AnVSr2QPp3cA85asMZDl76anDLYYCjYZov7y/Zjd9uJs8TRM7ZzTLMWQrRdZmYmS5Ys8d1XFIX33nsvZOtpiQRDnZC9vq26tX//0C5ENHbKU8dFt4GgUbuDipOeeqEu6S0PZwX4rrQGgAtS/dfTHa48DEC/hH4yqV4IobnPPvuMqVOnkp6efk4DJwmGOhlVVXEc88yUMvfoEeLViEbKj3jeJvTS5HSqqlJ5qg7Ab7NFgKNltQBkJEb5Pc47j6x7THcNViiEEI3V1NQwbNgwli5dek6fV4KhTsZdUYG7vh2Bubv8QtOVck/WhURtgqEzlXZcDjeKAjGJAYKh054MUo8ukX6P8wVDsfK9I0RrqKrKGbszJLdgh6quWbOGhIQE3G434BnoqigK9913n++YmTNncuutt/rub9iwgQkTJhAZGUlGRgZz5syhpqamzf9OU6ZM4YknnvANWD9XwmY3mdCGvT4rZExKwhDh/xekOMd8mSFtanGq6rNC0QlWjKaW/+5xu1WOlgeXGSqqLgIgPTpdkzUK0VnUOlxkPfJxSJ77m99NIsoS+Nf9hAkTqKqqYvv27YwaNYq8vDySkpLIy8vzHZObm8s999wDwK5du5g0aRKPP/44y5cvp6SkhNmzZzN79mxWrFjRYV9PR5DMUCfjOOoJhiySFdKfsvrMkEaXySpPeQKcuCT/2Z6Saht2pxujQSHNz04yaMgM9YiVS6xCnG/i4+MZPnw4ubm5QEPgs2PHDqqqqiguLmbfvn2+3n7PPfcc06ZNY968efTv359LLrmEF198kZUrV1JXVxe6L6QNJDPUyfjqhSQY0hdVhYr6ocAaZ4ZiuwaqF/JcIkuNi8Bk9JNBUt0NmaEYyQwJ0RqRZiPf/G5SyJ47WDk5OeTm5jJ//nw+//xznnjiCVatWsX69espLy8nJSWFgQMHArB161YOHDjAm2++6ft8VVVxu90UFBQwaNAgzb+WjtKmYKi8vJy///3vHDx4kPvuu48uXbqwbds2UlJS6C6/ZHXNcfQoIMXTulN9Apx1oBggXpvXprK0vng6QDBUeLr+ElmAeqGSMyU43A6MipGUqJaHHQshmlIUJahLVaGWk5PD8uXL2bFjBwaDgaysLLKzs8nLy6OsrIzs7GzfsW63m5kzZzJnzpwm5+nZM7xab7T6ldm5cydXXHEF8fHxHDp0iBkzZtClSxdWr17N4cOHWblyZUesU2hEMkM65a0XiusOGs37qqq/TBbb1X+Q480M9QhUL1TjyQqlRqdiMuj/h7oQovW8dUNLliwhOzsbRVHIzs5m0aJFlJWVMXfuXN+xI0eOZPfu3fTr1y+EK9ZGq2uG5s+fz/Tp09m/fz8RZxXgTpkyhc8++0zTxQnt2Y95M0MSDOmKxvVCcFZmKMC2el9mKEAwdLTK870j2+qFOH9564beeOMNX23QhAkT2LZtW6N6IYAFCxawceNGZs2aRX5+Pvv37+f999/n7rvvbvPzV1dX+waxAxQUFPiGtHekVgdDmzdvZubMmU0e7969O8XFxZosSnQMT48hz1/3UkCtM95t9RrVC6lularTQdYMlXszQ/4zSFIvJETnMHHiRFwuly/wSUxMJCsri+Tk5EZ1QEOHDiUvL4/9+/czfvx4RowYwcMPP0xaWlqbn3vLli2MGDHCN4x9/vz5jBgxgkceeaRdX1Mgrc51R0REUFlZ2eTxvXv3kpycrMmiRMdwnT6NWlsLioIpXX6h6YrGPYZqKmy4XSoGgxKwx1BDzZA0XBRCwOLFi1m8eHGjx7yZmu8bM2YMa9eubfFch+onHngF6nmUk5MTdF8kLbU6M3Tdddfxu9/9DofDAXiKwo4cOcIDDzzAj370I80XKLTjrRcypaRgsFhCvBrRiMbdp72dp2O6WDEYWh6b4XKrFJUHV0DtzQxJMCSEON+0OhhavHgxJSUldOvWjdraWrKzs+nXrx+xsbE8+eSTHbFGoRHfTjK5RKY/ZdpeJqsqDa54uriyDqdbxWxU6BYb4HJatdQMCSHOT62+TBYXF8f69ev55JNP2LZtG263m5EjR3LFFVd0xPqEhuzehotSPK0vbhdUeAINrS6TBTuTrLB+DEf3hEiMfjJITreT4hpPTaAEQ0KI802b98dedtllXHbZZVquRXQw2VavU1XHwe0Agxli2154eDZfMBSw4aIngxRoW/3JMydxqS7MBjPJUVIbKIQ4v7T6MtmcOXN48cUXmzy+dOlS5s2bp8WaRAdpuEwmDRd1xXuJLL4HGILvFOtPsJfJvJmhQPVC3uLp9Jh0DIpM8RFCnF9a/VNt1apVXHrppU0ev+SSS/j73/+uyaJEx5DMkE5pPKAVtM8M+YIhGdAqhDgPtToYOnXqFPHx8U0ej4uLo7S0VJNFCe2pbndDMCSjOPRF4231bpeb6jIbEHhIa2FZcD2GfNvqYyWQFkKcf1odDPXr14+PPvqoyeMffvghffr00WRRQnvOklJUhwOMRsypMldKVzTODFWX2VDdKkaTgag4/y0UjgWZGZJt9UKI81mrC6jnz5/P7NmzKSkp8RVQ//e//+X3v/89S5Ys0Xp9QiMO7xiO1FQUk8yV0hXftvpMTU5Xeda0esXPDjGHy83xiuB6DMkoDiHE+azVmaGf//zn/P73v2f58uVMnDiRiRMn8sYbb7Bs2TJmzJjREWsUGpB6IR3TODPUMKDVf73Q8fI63CpYTQaSY6x+j/UOaZVRHEKI1sjMzGyUKFEUhffeey9k62lJm7aF/OpXv+Lo0aOcOHGCyspKvvvuO26//Xat1yY05NtJJvVC+uJyQKXGPYZKgyuePrteSFH8ZZAcnKg5AUhmSAjRsRYtWsSYMWOIjY2lW7duXH/99ezdu7fDn7dde2STk5OJiYnRai2iA9l9mSH5y15XKo+B6gajFaK7aXLKqlNBDmj1BUP+64WKa4pRUYkwRtA1oqsmaxRCiObk5eUxa9YsNm3axLp163A6nVx11VXU1NR06PO2Ohg6ceIEt912G+np6ZhMJoxGY6Ob0CeHr/u0ZIZ05ewxHAZt+vdU1l8mC7iT7HSQ9UL1YzjSY9L9ZpCEEH6oKthrQnMLcvDpmjVrSEhIwO12A57hrIqicN999/mOmTlzJrfeeqvv/oYNG5gwYQKRkZFkZGQwZ86cdgUuH330EdOnT+fCCy9k2LBhrFixgiNHjrB169Y2nzMYra6knT59OkeOHOHhhx8mLS1NfjiGCakZ0qkO6DGkdWbIu5NM6oWEaAfHGXgqRP8P/bYILNEBD5swYQJVVVVs376dUaNGkZeXR1JSEnl5eb5jcnNzueeeewDYtWsXkyZN4vHHH2f58uWUlJQwe/ZsZs+ezYoVKzRZekVFBQBdunTR5HwtaXUwtH79ej7//HOGDx/eAcsRHUF1uXAcPw5IzZDuaNxjyOVwU11e32MoUPfp+m31GUE2XJR6ISHOb/Hx8QwfPpzc3FxGjRrlC3wee+wxqqqqqKmpYd++feTk5ADw3HPPMW3aNN/0if79+/Piiy+SnZ3NsmXLiIjw/wdZIKqqMn/+fH7wgx8wePDgdn51/rU6GMrIyEANMuUm9MFZXAxOJ5jNmJJlrpSuaD2t/nQdqGCyGIiMNfs91juKI1DDRZlWL4QGzFGeDE2onjtIOTk55ObmMn/+fD7//HOeeOIJVq1axfr16ykvLyclJYWBAwcCsHXrVg4cOMCbb77p+3xVVXG73RQUFDBo0KB2LXv27Nns3LmT9evXt+s8wWh1MLRkyRIeeOABXnnlFTIzMztgSUJrvuLp9DQUqevSF99lMm0yQw2XyPzvEKtzuDhZ5ckgZXSRhotCdDhFCepSVajl5OSwfPlyduzYgcFgICsri+zsbPLy8igrKyM7O9t3rNvtZubMmcyZM6fJeXr2bN8feHfffTfvv/8+n332GT3OwRWNVgdDN998M2fOnKFv375ERUVhNjf+6/P06dOaLU5ow1c8LQNa9Ufjy2QNxdP+09PHyj3HRVuMJEb5zyDJKA4hOg9v3dCSJUvIzs5GURSys7NZtGgRZWVlzJ0713fsyJEj2b17N/369dPs+VVV5e6772b16tXk5ubSu3dvzc7tT5syQyK8SPG0TjltUOWp5dIqM+Qb0Nol+AGtfjNIzjpKaz0zB7tHy/ePEOc7b93QG2+8wQsvvAB4AqSbbroJh8PhqxcCWLBgAWPHjmXWrFnMmDGD6Oho9uzZw7p163jppZfa9PyzZs3irbfe4p///CexsbEUFxf71hUZ6f+Sfnu0Ohi64447OmIdogNJw0WdKi/0vDVHQ5Q2/XuqSuu7TwfcVu+pFwq0rd7beTraHE28temAZiHE+WfixIls27bNF/gkJiaSlZVFUVFRozqgoUOHkpeXx0MPPcT48eNRVZW+ffty8803t/m5ly1bBtAo6AJYsWIF06dPb/N5A2nTkKqDBw+yYsUKDh48yAsvvEC3bt346KOPyMjI4MILL9R6jaKdJDOkU+VnFU9r1KLClxkKuK2+dQNapceQEJ3H4sWLWbx4caPH8vPzmz12zJgxrF27tsVzHTp0qNH9QBuwQrVBq9Vd3vLy8hgyZAhffvkl//jHP6iurgZg586dPProo5ovULSft4Da0kOCIV3RuF4IzgqGAmWGyoLbSXasSrbVCyHOf60Ohh544AGeeOIJ1q1bh8Vi8T0+ceJENm7cqOniRPupdjvOE565UpIZ0hmNGy467S5qK+1AMA0Xg8sMHauRYEgIcf5rdTC0a9cubrjhhiaPJycnc+rUKU0WJbTjKC4GtxvFasWYlBTq5Yiz+XoMaVs8bYkwYo3yfwX8aJA1Q5IZEkJ0Bq0OhhISEjhe3834bNu3b6e7ZB505+x6Ian50BmNM0PB9hg6Y3dyqsaTQZJRHEII0YZgaNq0aSxYsIDi4mIURcHtdvPFF19w7733cvvtt3fEGkU72H07ySRQ1R2tewyVBtdjyHuJLC7CRHxkcD2GesTITkQhxPmr1cHQk08+Sc+ePenevTvV1dVkZWUxYcIELrnkEv73f/+3I9Yo2kF2kumU/QzUlHje1zwzpM2A1jOOM5TZygDJDAkhzm+t3lpvNpt58803efzxx9m2bRtut5sRI0bQv3//jlifaCdf92npMaQv3ktk1niITNTklL7u04EGtJ6uH9AaqF6oPisUZ4kj1hKrwQqFEEKf2tRnCKBPnz706dNHy7WIDiCZIZ3SuF4ItM8MyUwyIURn0erLZD/+8Y95+umnmzz+3HPPcdNNN2myKKEdX/dpmUumLx3RY6g0yB5D3syQTKsXQgigjU0Xr7nmmiaPT548mc8++0yTRQltuOvqcJZ46lKkgFpnyg553mqUGbLXOamrcQCBu08XBpkZ8l4mk3ohIURbZWZmNpppqigK7733XsjW05JWB0PV1dWNmi16mc1mKisrNVmU0IajyNMCwRAVhTEhIbSLEY35LpNpkxnyXiKzRpuwRAboMVTmrRmSy2RCCH1ZtmwZQ4cOJS4ujri4OMaNG8eHH37Y4c/b6mBo8ODBvPvuu00ef+edd8jKytJkUc158sknueSSS4iKiiIhyF/s06dPR1GURrexY8d22Br1xnHMe4lMegzpztlzyTTg21YfoHi6ss5BRa0ngxRwFEe1NFwUQpxbPXr04Omnn2bLli1s2bKFyy67jOuuu47du3d36PO2uoD64Ycf5kc/+hEHDx7ksssuA+C///0vb7/9Nn/72980X6CX3W7npptuYty4cSxfvjzoz5s8eTIrVqzw3W8uq3W+8hVPy04y/fFmhrTqMRTsgNb6eqEu0Rairf7/95dgSAjtqKpKrbM2JM8dafLfiNVrzZo13HbbbZw+fRqDwUB+fj4jRozg3nvv5bnnngNg5syZVFZW8vbbbwOwYcMGHnjgATZv3kxSUhI33HADixYtIjo6uk1rnTp1aqP7Tz75JMuWLWPTpk0dOgi+1cHQtddey3vvvcdTTz3F3//+dyIjIxk6dCj/+c9/yM7O7og1AvDYY48B8Nprr7Xq86xWK6mpqR2wIv3zFU9LMKQvdZVQ6+nfo1mPofri6dggB7QGKp6utFdSZa8CpGZICC3UOmu5+K2LQ/LcX077kiiz/8viABMmTKCqqort27czatQo8vLySEpKIi8vz3dMbm4u99xzD+AZzzVp0iQef/xxli9fTklJCbNnz2b27NmNkhBt5XK5+Nvf/kZNTQ3jxo1r9/n8adPW+muuuabZImo9ys3NpVu3biQkJJCdnc2TTz5Jt27dWjzeZrNhs9l898O5Dsru21Yvv8x0xZsViuwCVm369zT0GNJmQKu3XqhLRJegfogKIcJffHw8w4cPJzc3l1GjRvkCn8cee4yqqipqamrYt28fOTk5gGcX+bRp05g3bx4A/fv358UXXyQ7O5tly5YREeH/51FLdu3axbhx46irqyMmJobVq1d3aBkOtCEYKiwsRFEUetRnG7766iveeustsrKy+OUvf6n5AttjypQp3HTTTfTq1YuCggIefvhhLrvsMrZu3YrVam32cxYtWuTLQoU7abioUx2xrT7IHkOF9QNae8iAViHOqUhTJF9O+zJkzx2snJwccnNzmT9/Pp9//jlPPPEEq1atYv369ZSXl5OSksLAgQMB2Lp1KwcOHODNN9/0fb6qqrjdbgoKChg0aFCb1nvBBReQn59PeXk5q1at4o477iAvL69DA6JWB0PTpk3jl7/8JbfddhvFxcVcccUVDB48mDfeeIPi4mIeeeSRoM+1cOHCgIHH5s2bGT16dGuXCcDNN9/se3/w4MGMHj2aXr168e9//5sbb7yx2c958MEHmT9/vu9+ZWUlGRkZbXr+UJOGizrVgQ0XAxVQB5sZkm31QmhLUZSwyLLm5OSwfPlyduzYgcFgICsri+zsbPLy8igrK2tUDuN2u5k5cyZz5sxpcp6ePdv+881isdCvXz8ARo8ezebNm3nhhRd45ZVX2nzOQFodDH399ddcdNFFAPz1r39lyJAhfPHFF6xdu5a77rqrVcHQ7NmzueWWW/wek5mZ2doltigtLY1evXqxf//+Fo+xWq0tZo3CibumBtfp04AEQ7pT5t1Jpk1mqK7Ggb3WCUBswCGtwdUMSfG0EJ2Tt25oyZIlZGdnoygK2dnZLFq0iLKyMubOnes7duTIkezevdsXuHQUVVUbla90hFYHQw6Hwxcs/Oc//+Haa68FYODAgRw/frxV50pKSiIpKam1S2izU6dOUVhYSFpa2jl7zlBxFHlqPgxxcRjj4kK8GtGIxpkhb1YoMtaM2WJs8ThVVVtdMyTBkBCdi7du6I033uCFF14APAHSTTfdhMPh8NULASxYsICxY8cya9YsZsyYQXR0NHv27GHdunW89NJLbXr+3/72t0yZMoWMjAyqqqp45513yM3N5aOPPtLiy2tRq/sMXXjhhbz88st8/vnnrFu3jsmTJwNQVFRE165dNV+g15EjR8jPz+fIkSO4XC7y8/PJz8+nurrad8zAgQNZvXo14GkOee+997Jx40YOHTpEbm4uU6dO9W39O9/ZfTvJ5JeZ7vhqhjI1OZ2veDrATrLyMw6qbZ4MUqAeQzKKQ4jOa+LEibhcLl/gk5iYSFZWFsnJyY3qgIYOHUpeXh779+9n/PjxjBgxgocffrhdCYcTJ05w2223ccEFF3D55Zfz5Zdf8tFHH3HllVe298vyq9WZoWeeeYYbbriB5557jjvuuINhw4YB8P777/sun3WERx55hNdff913f8SIEQB8+umnvhds7969VFRUAGA0Gtm1axcrV66kvLyctLQ0Jk6cyLvvvkts7Pk/gdtxzPOXvUUukemLqnZYZijwgFZP0JQcayXC7D+D5M0MSc2QEJ3P4sWLWbx4caPH8vPzmz12zJgxrF27tsVzHTp0qNF9VVX9Pndr+ghqqdXBUE5ODqWlpVRWVpKYmOh7/Je//CVRUR1XHPbaa68F7DF09j9yZGQkH3/8cYetR+9kQKtO1ZaBrb5dg2bdp4Mrng62x1C5rZwzTs+xEgwJITqDVl8mA0/QsXXrVl555RWqqjyN2SwWS4cGQ6J1zh7FIXTEe4ksuhuYg9/u6o/3MlnQ2+qD3EmWHJmM1Rj+mwmEECKQVmeGDh8+zOTJkzly5Ag2m40rr7yS2NhYnn32Werq6nj55Zc7Yp2ilXwNF6VmSF80HsMBZ22rD7iTzDugVXaSCSHE2VqdGZo7dy6jR4+mrKyMyMiGH6o33HAD//3vfzVdnGg7b8NFyQzpTJm2A1pVVQ16SKv3Mpn0GBJCiMZanRlav349X3zxRZOBp7169eJYfTZChJarshJ3/RgRKaDWGV/xtEY9hqodOO1uAGK7BJkZkm31QgjRSKszQ263G5fL1eTxo0ePdopdWuHA23namJiIoY2Tg0UHKdc2M+Qtno6Ot2A0t/y/s6fHkDczJNvqhRDibK0Ohq688kqWLFniu68oCtXV1Tz66KNcffXVWq5NtJFvDIfMJNMfjWuGgu0xVFptp87hRlEgPcH/sb7MUKwEQ0KIzqHVl8mef/55LrvsMrKysqirq2PatGns37+fpKQk3n777Y5Yo2glX8NFuUSmL416DGkTDAXbY8hbL5QaF4HF5D+D5AuGouX7RwjRObQ6GOrevTv5+fm88847bN26FbfbzS9+8Qt++tOfNiqoFqHja7goO8n0paYUHGcABeK1ydpV+naSBTegNVC90Km6U9hcNgyKgdSYVE3WKIQQeteqy2QOh4M+ffpQUFDAz372M5YuXcqf/vQn7rzzTgmEdMTXcFEuk+mLt14oLh1M2vTvqSptZY+hANvqj1Z5vndSolIwG8warFAI0ZllZmY2Ka157733QraelrQqGDKbzdhsNhRF6aj1CA34aobkMpm+aFw8DWdlhoIcxRHsgFbZVi+ECLVFixahKArz5s3r8OdqdQH13XffzTPPPIPT6eyI9Yh2UlVVRnHola/HkDb1QqpbPavhYqDLZMGN4pCGi0IIPdi8eTOvvvoqQ4cOPSfP1+pg6Msvv+Qf//gHPXv2ZNKkSdx4442NbiK0XOXluM94fvGZu8tf97qi8YDWM1V2XE7PDrHoRP+X3YLNDEkwJETHUFUV95kzIbkFGo7qtWbNGhISEnC7Pb3L8vPzURSF++67z3fMzJkzufXWW333N2zYwIQJE4iMjCQjI4M5c+ZQU1PTrn+r6upqfvrTn/LnP/+50QzUjtTqAuqEhAR+9KMfdcRahAa8nadNyckYrDJXSle8l8m02lZf32MoJjECo7Hlv2vcbpVjMopDiJBSa2vZO3JUSJ77gm1bUYKYHTphwgSqqqrYvn07o0aNIi8vj6SkJPLy8nzH5Obmcs899wCwa9cuJk2axOOPP87y5cspKSlh9uzZzJ49mxUrVrR5vbNmzeKaa67hiiuu4IknnmjzeVqj1cFQe75A0fGkXkjHNM4MVQU5oPVklQ27y43RoJAa5/9YqRkSovOKj49n+PDh5ObmMmrUKF/g89hjj1FVVUVNTQ379u0jJycHgOeee45p06b5anr69+/Piy++SHZ2NsuWLSMiwv/Pm+a88847bNu2jc2bN2v4lQXW6mDIq6SkhL1796IoCgMGDCA5OVnLdYk28k2rl51k+uJ2a95jyJsZCjSg1dtjKD0hApOfDJLL7aKoxhMM9YiR7x8htKRERnLBtq0he+5g5eTkkJuby/z58/n888954oknWLVqFevXr6e8vJyUlBQGDhwIwNatWzlw4ABvvvmm7/NVVcXtdlNQUMCgQYNatc7CwkLmzp3L2rVr2xRItUerg6GamhruvvtuVq5c6buuaDQauf3223nppZeICiIVJzqONFzUqepicNlBMUKcNq9NpS8zFGBAq3dbfYL//zdLaktwup2YFBPdorppskYhhIeiKEFdqgq1nJwcli9fzo4dOzAYDGRlZZGdnU1eXh5lZWVkZ2f7jnW73cycOZM5c+Y0OU/Pnq3PgG/dupWTJ08yalTD5USXy8Vnn33G0qVLsdlsGI3Gtn1hAbQ6GJo/fz55eXmsWbOGSy+9FPAMb50zZw6/+c1vWLZsmeaLFMFrGMUhwZCueLNC8d3B2OaEbCMNO8mCHNAaZL1QanQqRkPH/MARQuibt25oyZIlZGdnoygK2dnZLFq0iLKyMubOnes7duTIkezevZt+/fpp8tyXX345u3btavTYz372MwYOHMiCBQs6LBCCNgRDq1at4u9//7vvmiHA1VdfTWRkJD/5yU8kGAoxbwG1TKvXGY231QNU1jdcDNRjyJcZkp1kQogAvHVDb7zxBi+88ALgCZBuuukmHA5Ho9/9CxYsYOzYscyaNYsZM2YQHR3Nnj17WLduHS+99FKrnzs2NpbBgwc3eiw6OpquXbs2eVxrrd5af+bMGVJSUpo83q1bN87Ub+kWoaGqqgxp1SuN64XcbpXq0zYg8GWy1maGZECrEJ3bxIkTcblcvsAnMTGRrKwskpOTG9UBDR06lLy8PPbv38/48eMZMWIEDz/8MGlpaSFaedu1OjM0btw4Hn30UVauXOkrcKqtreWxxx5j3Lhxmi9QBM9VWopqs4GiYE6VuVK6Un7I81ajnWQ15TbcbhWDUSE6wX8LBW8BdcDMUJUnGEqPlp1kQnRmixcvZvHixY0ey8/Pb/bYMWPGsHbt2hbPdejQoUb3g+155JWbm9uq49uq1cHQCy+8wOTJk+nRowfDhg1DURTy8/OJiIjg448/7og1iiB5s0Km1FQUiyXEqxGNeDNDGvUY8m6rj+kSgcHQ8ngcp8vN8QpPbVGgIa3enWSSGRJCdDatDoYGDx7M/v37eeONN/j2229RVZVbbrlFptbrgP2ot8eQ/GWvO2XaziXzbasPUC90vKIOl1vFYjTQLdZ/BsmbGZKaISFEZ9OmbS2RkZHMmDFD67WIdvJmhiwyk0xfXE6o9Lw2mvUYauWA1u6Jkf4zSG4nJ86c8BwrwZAQopNpdTC0cuVKvx+//fbb27wY0T4O6TGkT1VF4HaCwQyx2tRyVZUG2WPIVy/k/7jimmJcqguLwUJSZJImaxRCiHDR6mDo7B4DAA6HgzNnzmCxWIiKipJgKIRkJ5lO+XaSZYBG/XsqW9ljKFDx9NljOAxKqzeZCiFEWGv1T72ysrJGt+rqavbu3csPfvAD3n777Y5YowiS3TeKQzJDutIRPYbqC6jjkgJsq6/vMSQDWoUQomWa/AnYv39/nn766SZZI3HuqKqKs+g4AOZ0+YWmKxWeIJWEDE1O53ar1JTV9xjq4j8zdKy8vmYowX8wdLzG872TFhN+/UGEEKK9NMuHG41GioqKtDqdaCV3TQ2qwwGAqWuXEK9GNHKm1PM2WpthxrYzDrytOiJizX6PPV1jByApxv9OstN1pwHoEiHfO0KIzqfVNUPvv/9+o/uqqnL8+HGWLl3qm1Umzj1XeTkASkQEBmlxoC9nPIEGUV01OV1dtSfotUSaMPqZQg9QdsZzbGKU/75TFbYKz3HWRA1WKIQQ4aXVwdD111/f6L6iKCQnJ3PZZZfx+9//Xqt1iVZylZUBYExICO1CRFNnTnneRmqTdfEGQxHR/v/3VVWV8jOezFBitP8MUpnN8/0Tb43XYIVCCOGRmZnJvHnzmDdvHuCJGVavXt0klgi1VgdDbre7I9Yh2smbGTImyl/2ulOrcWaopj4YivGf7amyOXG6PdfTgs4MRcj3jxAidBYuXMhjjz3W6LGUlBSKi4s79Hnb1HRR6E9DZkj+stcd32UybTJDtb7MUIBsT329UKTZSITZ/5b+sjrP90+CNaH9CxRCiHa48MIL+c9//uO7bzRq05LEn6CCofnz5wd9wueff77NixFt580MmSQzpD8aB0PezFBkTIBgqL5eqEt04Dl13syQBENCdAxVVXHaQ3NlxWQxoCgtd6D3WrNmDbfddhunT5/GYDCQn5/PiBEjuPfee3nuuecAmDlzJpWVlb5WOhs2bOCBBx5g8+bNJCUlccMNN7Bo0SKio6Pbvl6TidRzPGw8qGBoxYoVDB48GJPJhKIoLU6dDeYfW3QMp9QM6ZOjDhw1nvc1rxkKLjOUEOX/uFpnLXUuTxNHCYaE6BhOu5tX5+aF5Ll/+UI2Zmvg7MqECROoqqpi+/btjBo1iry8PJKSksjLa1h3bm4u99xzDwC7du1i0qRJPP744yxfvpySkhJmz57N7NmzWbFiRZvXu3//ftLT07FarVx88cU89dRT9OnTp83nC0ZQwVBFRQWrVq2iW7du9OnTh82bN9O1qzb1D0IbvpqhBMkM6Yq3XkgxQoQ2lzAbaoYCZYY8wVCgzJA3K2QymIg2t/2vOSFEeIuPj2f48OHk5uYyatQoX+Dz2GOPUVVVRU1NDfv27SMnJweA5557jmnTpvmKo/v378+LL75IdnY2y5YtIyLCfx+05lx88cWsXLmSAQMGcOLECZ544gkuueQSdu/e3aFxR1DBUGJiIgUFBXTr1o1Dhw5JEbUOucrKAckM6Y53J1lUF9Aoc+rLDAUIhk77MkP+g6Gz64UkuytExzBZDPzyheyQPXewcnJyyM3NZf78+Xz++ec88cQTrFq1ivXr11NeXk5KSgoDBw4EYOvWrRw4cIA333zT9/mqquJ2uykoKGDQoEGtXuuUKVN87w8ZMoRx48bRt29fXn/99VaV7LRWUMHQj370IyZMmEB6ejqKojB69OgWC5q+++47TRcoguMroJaaIX3xBUPa/UUT9GUyb2YowGUy77Z6uUQmRMdRFCWoS1WhlpOTw/Lly9mxYwcGg4GsrCyys7PJy8ujrKyM7OyGgM7tdjNz5kzmzJnT5Dw9e/bUZD3R0dEMGTKE/fv3a3K+lgQVDL366qvceOONHDhwgDlz5jBjxgxiY2M7dGGidRoukyWEdB3ie7zF0xrVC0FrLpN5jguUGZLiaSGEl7duaMmSJWRnZ6MoCtnZ2SxatIiysrJGY7dGjhzJ7t276devX4etx2azsWfPHsaPH99hzwGt2Fo/efJkwJMWmzt3rgRDOiOZIZ06+zKZRrxb6wPuJqsJrmbIe5lMegwJIbx1Q2+88QYvvPAC4AmQbrrpJhwOh69eCGDBggWMHTuWWbNmMWPGDKKjo9mzZw/r1q3jpZdeatPz33vvvUydOpWePXty8uRJnnjiCSorK7njjju0+PJa1OrZZCtWrJBASGdUVT1ra31CSNcivqfWE2hodZlMdavYWllAnRhkAbVkhoQQABMnTsTlcvkCn8TERLKyskhOTm5UBzR06FDy8vLYv38/48ePZ8SIETz88MOkpbV94PPRo0e59dZbueCCC7jxxhuxWCxs2rSJXr16tffL8kuaLp4H1DNnUO2eX3xymUxnNM4M2WqdDUNaA26t984lk5ohIUTwFi9ezOLFixs9lp+f3+yxY8aMYe3atS2e69ChQ43ut9Sax+udd94Jao1a02xqvQgd35BWiwUlKiq0ixGNddCQVnOEEaMp0JDW+sxQgJqhcls5IMGQEKLzkmDoPOA8a1u9bI3WGa2HtNYEt5NMVdWgL5OV15UDkBCR0O71CSFEOJJg6DwgQ1p1TOshrUEWT9fYXThcnnR0F8kMCSGEXxIMnQdcMopDvzSuGWrtkFaryUCkxX9vEwmGhBCdXVgEQ4cOHeIXv/gFvXv3JjIykr59+/Loo49iry8abomqqixcuJD09HQiIyPJyclh9+7d52jV545khnTsjLa7ybQexQENwVCiVb5/hBCdU1gEQ99++y1ut5tXXnmF3bt384c//IGXX36Z3/72t34/79lnn+X5559n6dKlbN68mdTUVK688kqqqqrO0crPjYbMkDazr4RGnHaw13+vRWoTaATbfTrYURw2l41aZy0A8RrNThNCiHATFlvrJ0+e7Gv6CNCnTx/27t3LsmXLmmz/81JVlSVLlvDQQw9x4403AvD666+TkpLCW2+9xcyZM8/J2s8FyQzplG9IqwE0Kk4ONjNUXt99ukuAoMlbPG1UjMSapX+YEKJzCovMUHMqKiro0qXlOoyCggKKi4u56qqrfI9ZrVays7PZsGFDi59ns9morKxsdNM7V7knM2SSmiF98e0kSwSDNv+raZ0Z8l4ii7fGy05EIUSnFZbB0MGDB3nppZe46667WjymuLgYgJSUlEaPp6Sk+D7WnEWLFhEfH++7ZWRkaLPoDiSZIZ3SuMcQtCYz5B3SGlwwJPVCQojOLKTB0MKFC1EUxe9ty5YtjT6nqKiIyZMnc9NNN3HnnXcGfI7v/7Wrqqrfv4AffPBBKioqfLfCwsK2fXHn0Nl9hoSOaNxjCM7aTRYgGDrta7gYXPfpeKvUCwkhtJeZmcmSJUt89xVF4b333gvZeloS0pqh2bNnc8stt/g9JjMz0/d+UVEREydOZNy4cbz66qt+Py81NRXwZIjOnpNy8uTJJtmis1mtVqxWaxCr1w/JDOmUxj2GoCEzFHBIa33NUMC5ZHWeuWQypFUIoRfHjh1jwYIFfPjhh9TW1jJgwACWL1/OqFGjOuw5QxoMJSUlkZSUFNSxx44dY+LEiYwaNYoVK1ZgCFCD0bt3b1JTU1m3bh0jRowAwG63k5eXxzPPPNPuteuJ9BnSKV+PIW0CDVVVsbWyz1CgURwyl0wIoSdlZWVceumlTJw4kQ8//JBu3bpx8OBBEjr491tY1AwVFRWRk5NDRkYGixcvpqSkhOLi4ia1PwMHDmT16tWAJxU3b948nnrqKVavXs3XX3/N9OnTiYqKYtq0aaH4MjqEu7YWta4OkMyQ7mjcY8he58Lt9nSVDhgMBZsZkon1QpwTqqriqKsLyS3QcFSvNWvWkJCQgNvtBjzDWRVF4b777vMdM3PmTG699Vbf/Q0bNjBhwgQiIyPJyMhgzpw51NTUtPnf6ZlnniEjI4MVK1Zw0UUXkZmZyeWXX07fvn3bfM5ghMXW+rVr13LgwAEOHDhAjx49Gn3s7Bd57969VFRU+O7ff//91NbW8utf/5qysjIuvvhi1q5dS2zs+bOF2HuJDJMJQ3R0SNcivkfruWTVnmyPyWLAFKCrdENmSCbWC6EHTpuNF+/4cUiee87rf8ccERHwuAkTJlBVVcX27dsZNWoUeXl5JCUlkZeX5zsmNzeXe+65B4Bdu3YxadIkHn/8cZYvX05JSQmzZ89m9uzZrFixok1rff/995k0aRI33XQTeXl5dO/enV//+tfMmDGjTecLVlhkhqZPn46qqs3ezqaqKtOnT/fdVxSFhQsXcvz4cerq6sjLy2Pw4MHnePUdq6FeSIa06o7mc8mcQODi6UZDWoOdSyZDWoXo9OLj4xk+fDi5ublAQ+CzY8cOqqqqKC4uZt++feTk5ADw3HPPMW3aNObNm0f//v255JJLePHFF1m5ciV19VcsWuu7775j2bJl9O/fn48//pi77rqLOXPmsHLlSo2+yuaFRWZItMxbLyQ9hnRI87lkngAn0CWyWocLm9OT5g56Yr1khoToUCarlTmv/z1kzx2snJwccnNzmT9/Pp9//jlPPPEEq1atYv369ZSXl5OSksLAgQMB2Lp1KwcOHODNN9/0fb6qqrjdbgoKChg0aFCr1+p2uxk9ejRPPfUUACNGjGD37t0sW7aM22+/vdXnC5YEQ2HOlxlKkHoh3dG4z5CtlTvJLEYD0TKkVQhdUBQlqEtVoZaTk8Py5cvZsWMHBoOBrKwssrOzycvLo6ysjOzsbN+xbrebmTNnMmfOnCbn6dmzZ5uePy0tjaysrEaPDRo0iFWrVrXpfMGSYCjMOWUnmX55gyGNaoZaO7E+Icoc8NKpBENCiLN564aWLFlCdnY2iqKQnZ3NokWLKCsrY+7cub5jR44cye7du+nXr59mz3/ppZeyd+/eRo/t27ePXr16afYczQmLmiHRMukxpFMuB9Tv1NJ+Yn2A7fJBTqx3uBzUODy7PqTPkBACGuqG3njjDV9t0IQJE9i2bVujeiGABQsWsHHjRmbNmkV+fj779+/n/fff5+67727z899zzz1s2rSJp556igMHDvDWW2/x6quvMmvWrHZ+Zf5JMBTmXNJ9Wp9q67fVo0BkgianbJhL5j+he/qszJA/3qyQQTEQazl/dlgKIdpn4sSJuFwuX+CTmJhIVlYWycnJjeqAhg4dSl5eHvv372f8+PGMGDGChx9+uFGj49YaM2YMq1ev5u2332bw4ME8/vjjLFmyhJ/+9Kft/bL8kstkYc7XcFEyQ/ri21afAAb/dTvB8gVDgTJDNcFlhnyjOCzxGBT5u0gI4bF48WIWL17c6LH8/Pxmjx0zZgxr165t8VyHDh1qdD+Ynkc//OEP+eEPfxjwOC3JT8Aw11BAnRDSdYjv0bheCM6+TOb/bxhvAXWgifXehosyl0wI0dlJMBTmGjJDCaFdiGjMt61eu7lk3gLqyEAZnyAn1pfVeb53pF5ICNHZSTAU5ryZIZNcJtOXDhzSGqjpYrCjOGQnmRBCeEgwFOaccplMnzRuuKiq6lk1Q8EOaQ2ugFqCISFEZyfBUBhz22yoZ84AUkCtO76Gi9oEQw6bC7erfkhrwMxQfTAUbGZIRnEI0SG8A09Fx9Hq31h2k4Ux35BWoxHDeTR89rygcQG1NytkNBswBz2kVUZxCBEKFosFg8FAUVERycnJWCwWmR2pMVVVsdvtlJSUYDAYsFj8/7wLRIKhMHb2TjL5H01ntB7SGuQoDmioGQpUQO3NDCVaJasohJYMBgO9e/fm+PHjFBUVhXo557WoqCh69uyJwdC+C10SDIUxl4zi0C/Nh7R6AhxrgFEcdQ4XtQ4XAAkBjvUGQ7K1XgjtWSwWevbsidPpxOVyhXo55yWj0YjJZNIkGSDBUBhrGMWRENJ1iGZoPKTVe5ks8JBWzyUyk0Eh1ur/f29fZki21gvRIRRFwWw2YzYHzuiK0JIC6jAmmSEd83Wg1rZmKNCQ1oZRHIFrFLw1Q5IZEkJ0dhIMhTHpMaRTLifUddSQ1gCXvrz1QgGCJofbQZWjCpCaISGEkGAojDklM6RPdeVA/fydSG0CjbZkhvzxjuJQUIizxGmwQiGECF8SDIWxht1k8pe9rnjrhSLiwahNWV7wmaHgRnF4g6E4axxGjQbJCiFEuJJgKIy5ysoByQzpjsb1QtCwmyxwZsg7iiNAoXX9XDLpMSSEEBIMhbWG3WSSGdKVDpxLFuxuskANF72ZIQmGhBBCgqGwJrvJdErjHkNA8HPJggyGymySGRJCCC8JhsKY9BnSKY17DKmq2lAzFOjyl0ysF0KIVpNgKEypdjvu6mpAMkO6o3HNkNPuxuXwDCPUbGK9zCUTQggfCYbClKuivo+NwYAxTrZG60qtthPrvVkhg0nBbA0wpFUm1gshRKtJMBSmfD2G4uJQjLI1WlfOaBwMnbWTLFBX6aAn1stlMiGE8JFgKEzJTjIdC9FcMpvTRY3dMxBSJtYLIUTwJBgKU9JjSMc0rhmqrc/2BCqe9o7iMCgQGxHckFaZSyaEEBIMhS3JDOmYxn2G6qqdQOu21RsMAYa0ysR6IYTwkWAoTEmPIZ1yu6HW89poVzMUXGaoYS6Z/+NcbheVtkpAMkNCCAESDIUt6TGkU3XloHq2wWt1mayuJrjMUMPEev/1QpX2StT6QbISDAkhhARDYUsyQzrlLZ62xILJf1ASrNZnhoLrPh1rjsVs8H9OIYToDCQYClPOcs8vNJPUDOlLB4ziqA1yN5l3W33AnWTehovSY0gIIQAJhsKWFFDrVAcOaY2ICZDxaeUoDtlWL4QQHhIMhSnZWq9THTmkNeBcsiBHcci2eiGEaESCoTAlmSGd0rjhIpydGQoyGAo2MyTb6oUQApBgKCypTifuSs/WaMkM6YzmQ1pdOO2tHdIaXM2QZIaEEMJDgqEw5BvSqigypFVvtG646B3SalCwRAQa0urdWh/cZTKpGRJCCA8JhsKQd1u9IS4OxeR/7II4x3yXybQJNLw7yawxwQ9pDXZrvWSGhBDCQ4KhMOStFzLJJTL90XpIa01w2+odLjdVNk9zxkBb6ytsnsyi1AwJIYSHBENhyCkNF/VL45qh1u4kUxSIiwxwbJ3n+yfBmtD+BQohxHlAgqEwJDvJdEzzIa3B7STzjuJIiDRjDDCk1ZsZkmBICCE8JBgKQ9JjSKdU9azLZFrNJQsuM3Q6yJ1kbtVNhV2CISGEOJsEQ2FIMkM6VVcBqsvzvtaXyQJmhoLrMVRlr8JdP0hWgiEhhPCQYCgMyZBWnfLWC5mjwRyhySlrg6wZOl2fQQrUfdpbLxRtjsZslCGtQggBYRIMHTp0iF/84hf07t2byMhI+vbty6OPPordbvf7edOnT0dRlEa3sWPHnqNVd5yGzFBCSNchvqfWE2ho2X3aFuRusoZRHMF1n5askBBCNAiLJjXffvstbrebV155hX79+vH1118zY8YMampqWLx4sd/PnTx5MitWrPDdt1j8/7IIB5IZ0infXDLtLl8GmxnydZ8OchSHBENCCNEgLIKhyZMnM3nyZN/9Pn36sHfvXpYtWxYwGLJaraSmpnb0Es8pX58hqRnSl5DOJfNeJgsyGIpIaPfahBDifBEWl8maU1FRQZcugYtUc3Nz6datGwMGDGDGjBmcPHnS7/E2m43KyspGN72RzJBOadxjCDpgYn39XDLJDAkhRIOwDIYOHjzISy+9xF133eX3uClTpvDmm2/yySef8Pvf/57Nmzdz2WWXYbPZWvycRYsWER8f77tlZGRovfx2UV0uXN4hrZIZ0heNewy5HG4cNs/uNM0n1stcMiGE8AlpMLRw4cImBc7fv23ZsqXR5xQVFTF58mRuuukm7rzzTr/nv/nmm7nmmmsYPHgwU6dO5cMPP2Tfvn38+9//bvFzHnzwQSoqKny3wsJCTb5WrbgqKz39bABjvMyW0hVfzZC2PYYUBayR/q9oBz2xvj4YkrlkQgjRIKQ1Q7Nnz+aWW27xe0xmZqbv/aKiIiZOnMi4ceN49dVXW/18aWlp9OrVi/3797d4jNVqxWq1tvrc54q34aIhNhbFLFujdaWD5pJFxJhRAnSVlon1QgjRdiENhpKSkkhKSgrq2GPHjjFx4kRGjRrFihUrMBhan9Q6deoUhYWFpKWltfpz9cJVLvVCuuUNhiK1nVgfqF7I6XJTUVs/jiPQxPr6PkPxEZIZEkIIr7CoGSoqKiInJ4eMjAwWL15MSUkJxcXFFBcXNzpu4MCBrF69GoDq6mruvfdeNm7cyKFDh8jNzWXq1KkkJSVxww03hOLL0IR0n9axEM0l8wZC4JlN5vdY78R6yQwJIYRPWGytX7t2LQcOHODAgQP06NGj0cfU+voZgL1791JR4flhbzQa2bVrFytXrqS8vJy0tDQmTpzIu+++S2xs7Dldv5YadpLJX/a600E1Q8HuJIuLMGEy+v/7pswmE+uFEOL7wiIYmj59OtOnTw943NmBUWRkJB9//HEHrio0pMeQTjUa0npuM0MN9UL+L5GpqioT64UQohlhcZlMNGjIDEkwpCu2KnDXX67SeEhroFEcp4PsPl3lqMJVP0hWmi4KIUQDCYbCjFPmkumTt17IFAmWKE1O6b1MZg20QyzIuWQVdZ6sUKQpEqtRvzsmhRDiXJNgKMx4t9ZLZkhnfPVC2o3iqA06MxTcKA5vvZAUTwshRGMSDIUZGcWhU756Ie0CjbpqT8ZHs1Ec0nBRCCGaJcFQmJGt9TrVoUNaA2R8WjmxPjFCvneEEOJsEgyFGckM6VSHDmkNMIojyJohX8NFyQwJIUQjEgyFEdXtxuXtoyQF1Pqi9ZBWlxt7nWfnV2SgzFCQozik4aIQQjRPgqEw4q6sBLcbAJNkhvRF64aL9VkhFLBEBTekNeAoDmm4KIQQzZJgKIx464UM0dEoFv+/+MQ51lFDWqPMGAIOafUEQ4GaLvoaLkqPISGEaESCoTDilHoh/dK4ZijY7tMut0q5b0hrgF1ndZIZEkKI5kgwFEZkJ5mO1XoCjXM9l6yy1oF3Ck2gAmrvbjIJhoQQojEJhsJIQ8PFhJCuQzSjg2qGAmWGTtdfIou1mjAHGNIqwZAQQjRPgqEwIpkhneqIIa01wQVDvlEcQQxplT5DQgjRPAmGwoj0GNIpew24bJ73NaoZqq0O7jJZwygO/8fVOGpwup2A9BkSQojvk2AojLhkSKs+eXsMGa1gidbklLYg55KVBZkZ8maFIowRRJoi279AIYQ4j0gwFEZc5ZIZ0qWz64UU/9vgg1UbZAG1bxRHkMXTkhUSQoimJBgKI94CapPUDOlLR8wlC7KA2tt9OthgSOqFhBCiKQmGwohTMkP65A2GIrWcWN/azFBwPYYkMySEEE1JMBRGXOXeuWTy172uaDyXDILfTRZszZDMJRNCiJZJMBQmVFVtKKCWzJC+aNxjyO1yYzvj2fkVMDMU7MR6m2SGhBCiJRIMhQl3dTU4Pb8gJRjSGY1rhryBEEBEdIAhrd6aoWAn1kvNkBBCNCHBUJjw9hhSIiMxRESEeDWiEY3nknl7DFmjTBgCdJUOdjeZzCUTQoiWSTAUJqTHkI5pXDMU7Fwy91lDWoOeWC/BkBBCNCHBUJiQ7tM6FqK5ZFV1Tlxuz5TWgBPrbZIZEkKIlkgwFCa8mSFTgtR86M6ZDppYH+ROsmiLEavJ6PdY35DWiIR2r08IIc43EgyFCac3MyTb6vXHlxnS6DKZdxRHoLlkrRnSWlfuOVa21gshRBMSDIUJ2VavU/Yz4Kz1vK9RAbU3GLIGO7E+QPF0rbMWu9tzrFwmE0KIpiQYChPeURySGdIZb/G0wQzWWG1OWRPckFbfxPogh7RaDBYZ0iqEEM2QYChMSGZIp3w9hrQb0hrsKI6GzFCA47z1QtYEFI3WKIQQ5xMJhsKEbzeZbK3XF43rhSD43WSng51YX18vJMXTQgjRPAmGwoRsrdcpjRsuQvB9hlo7ikPqhYQQonkSDIUJ39Z6qRnSl1ptt9VD8Jmhshpvw8XgL5MJIYRoSoKhMKCqKk6pGdInrYe0ulXqzgSXGfJurU8IdJlMgiEhhPBLgqEw4K45Aw7PL0jZTaYzGg9ptZ9xgqepdMDMkLeAOtAoDqkZEkII/yQYCgOu8vohrVYrhkjZGq0rmg9p9QQ4lggjxgBDWr1b6wON4pDMkBBC+CfBUBiQHkM6pvmQVicQOCukqmrwmSEJhoQQwi8JhsKANzMk9UI6pPmQVk+AE6heqMrmxFk/pDXg1noJhoQQwi8JhsKAr+Gi9BjSH9+QVq0yQ96dZAECHG+XarORCHNwQ1oTIySzKIQQzZFgKAxIjyEd89UMaRNo1Pq21Zv8Hnc6yO7T0FBAHW+Nb9/ihBDiPCXBUBiQHkM65agDR43nfY0yQzZvxidAHVBZkBPra5211LnqPMfKxHohhGiWBENhwCmZIX3yFk8rRojQJusSbGaoLMhRHBW2CgBMioloc7QGKxRCiPOPBENhoGFIq/xlrysdOaQ1QM1Q2ZnWTaxPiJAhrUII0RIJhsJAw9b6hJCuQ3xPKOeS1QRXM1RWJ3PJhBAiEAmGwoBkhnRK4x5D0Iq5ZEEOafVeJpNgSAghWibBUBiQ3WQ6pXGPIWjLxPoAx8nEeiGECChsgqFrr72Wnj17EhERQVpaGrfddhtFRUV+P0dVVRYuXEh6ejqRkZHk5OSwe/fuc7RibaiqelafIckM6coZbSfWq27V14E6MsiJ9a2pGRJCCNG8sAmGJk6cyF//+lf27t3LqlWrOHjwID/+8Y/9fs6zzz7L888/z9KlS9m8eTOpqalceeWVVFVVnaNVt59aW4tqswGSGdIdjWuGbLVO1Pqu0sFnhoIc0iqZISGEaJH//bs6cs899/je79WrFw888ADXX389DocDs7npLw5VVVmyZAkPPfQQN954IwCvv/46KSkpvPXWW8ycOfOcrb05BXt2UbR/X8DjlMoqzBFdwGRk046NgOwI0ovoQ8eJcKdSXWahdvvOdp/PVulCVR0YLApbC7/1e2xZTSFmgwO78zhHTlW3eNyJ8iJMToV4JQZHXV271yiEEB3BZLWGdMeroqqqGrJnb6PTp0/zq1/9imPHjrF+/fpmj/nuu+/o27cv27ZtY8SIEb7Hr7vuOhISEnj99deb/TybzYatPhMDUFlZSUZGBhUVFcTFxWn2NSyf+QDl5V9rdj4hhBAiXM15/e+YIyI0PWdlZSXx8fFB/f4Om8tkAAsWLCA6OpquXbty5MgR/vnPf7Z4bHFxMQApKSmNHk9JSfF9rDmLFi0iPj7ed8vIyNBm8d8jLV+EEEIIfQhpZmjhwoU89thjfo/ZvHkzo0ePBqC0tJTTp09z+PBhHnvsMeLj4/nXv/7VbGptw4YNXHrppRQVFZGWluZ7fMaMGRQWFvLRRx81+3znKjOkqirOs55HCCGE6Kw64jJZazJDIa0Zmj17NrfccovfYzIzM33vJyUlkZSUxIABAxg0aBAZGRls2rSJcePGNfm81NRUwJMhOjsYOnnyZJNs0dmsVitWq7WVX0nrKYqieUpQCCGEEK0X0mDIG9y0hTehZWshu9K7d29SU1NZt26dr2bIbreTl5fHM88807YFCyGEEOK8ExY1Q1999RVLly4lPz+fw4cP8+mnnzJt2jT69u3bKCs0cOBAVq9eDXgyL/PmzeOpp55i9erVfP3110yfPp2oqCimTZsWqi9FCCGEEDoTFlvrIyMj+cc//sGjjz5KTU0NaWlpTJ48mXfeeafRJa29e/dSUVHhu3///fdTW1vLr3/9a8rKyrj44otZu3YtsbGxofgyhBBCCKFDYbm1/lxqTQGWEEIIIfThvN1aL4QQQgihNQmGhBBCCNGpSTAkhBBCiE5NgiEhhBBCdGoSDAkhhBCiU5NgSAghhBCdmgRDQgghhOjUJBgSQgghRKcmwZAQQgghOrWwGMcRSt4G3ZWVlSFeiRBCCCGC5f29HcygDQmGAqiqqgIgIyMjxCsRQgghRGtVVVURHx/v9xiZTRaA2+2mqKiI2NhYFEXR9NyVlZVkZGRQWFgoc890RF4X/ZLXRp/kddGnzv66qKpKVVUV6enpGAz+q4IkMxSAwWCgR48eHfoccXFxnfIbVe/kddEveW30SV4XferMr0ugjJCXFFALIYQQolOTYEgIIYQQnZoEQyFktVp59NFHsVqtoV6KOIu8Lvolr40+yeuiT/K6BE8KqIUQQgjRqUlmSAghhBCdmgRDQgghhOjUJBgSQgghRKcmwZAQQgghOjUJhkLkT3/6E7179yYiIoJRo0bx+eefh3pJnc5nn33G1KlTSU9PR1EU3nvvvUYfV1WVhQsXkp6eTmRkJDk5OezevTs0i+1EFi1axJgxY4iNjaVbt25cf/317N27t9Ex8tqce8uWLWPo0KG+Bn7jxo3jww8/9H1cXhN9WLRoEYqiMG/ePN9j8toEJsFQCLz77rvMmzePhx56iO3btzN+/HimTJnCkSNHQr20TqWmpoZhw4axdOnSZj/+7LPP8vzzz7N06VI2b95MamoqV155pW9enegYeXl5zJo1i02bNrFu3TqcTidXXXUVNTU1vmPktTn3evTowdNPP82WLVvYsmULl112Gdddd53vl6q8JqG3efNmXn31VYYOHdrocXltgqCKc+6iiy5S77rrrkaPDRw4UH3ggQdCtCIBqKtXr/bdd7vdampqqvr000/7Hqurq1Pj4+PVl19+OQQr7LxOnjypAmpeXp6qqvLa6EliYqL6//7f/5PXRAeqqqrU/v37q+vWrVOzs7PVuXPnqqoq/78ESzJD55jdbmfr1q1cddVVjR6/6qqr2LBhQ4hWJb6voKCA4uLiRq+T1WolOztbXqdzrKKiAoAuXboA8trogcvl4p133qGmpoZx48bJa6IDs2bN4pprruGKK65o9Li8NsGRQa3nWGlpKS6Xi5SUlEaPp6SkUFxcHKJVie/zvhbNvU6HDx8OxZI6JVVVmT9/Pj/4wQ8YPHgwIK9NKO3atYtx48ZRV1dHTEwMq1evJisry/dLVV6T0HjnnXfYtm0bmzdvbvIx+f8lOBIMhYiiKI3uq6ra5DERevI6hdbs2bPZuXMn69evb/IxeW3OvQsuuID8/HzKy8tZtWoVd9xxB3l5eb6Py2ty7hUWFjJ37lzWrl1LREREi8fJa+OfXCY7x5KSkjAajU2yQCdPnmwSuYvQSU1NBZDXKYTuvvtu3n//fT799FN69Ojhe1xem9CxWCz069eP0aNHs2jRIoYNG8YLL7wgr0kIbd26lZMnTzJq1ChMJhMmk4m8vDxefPFFTCaT799fXhv/JBg6xywWC6NGjWLdunWNHl+3bh2XXHJJiFYlvq93796kpqY2ep3sdjt5eXnyOnUwVVWZPXs2//jHP/jkk0/o3bt3o4/La6Mfqqpis9nkNQmhyy+/nF27dpGfn++7jR49mp/+9Kfk5+fTp08feW2CIJfJQmD+/PncdtttjB49mnHjxvHqq69y5MgR7rrrrlAvrVOprq7mwIEDvvsFBQXk5+fTpUsXevbsybx583jqqafo378//fv356mnniIqKopp06aFcNXnv1mzZvHWW2/xz3/+k9jYWN9ftPHx8URGRvp6qMhrc2799re/ZcqUKWRkZFBVVcU777xDbm4uH330kbwmIRQbG+urp/OKjo6ma9euvsfltQlC6DaydW5//OMf1V69eqkWi0UdOXKkb9uwOHc+/fRTFWhyu+OOO1RV9WxJffTRR9XU1FTVarWqEyZMUHft2hXaRXcCzb0mgLpixQrfMfLanHs///nPfT+zkpOT1csvv1xdu3at7+PymujH2VvrVVVem2AoqqqqIYrDhBBCCCFCTmqGhBBCCNGpSTAkhBBCiE5NgiEhhBBCdGoSDAkhhBCiU5NgSAghhBCdmgRDQgghhOjUJBgSQgghRKcmwZAQQgghOjUJhoQQohm5ubkoikJ5eXmolyKE6GASDAkhhBCiU5NgSAghhBCdmgRDQoiw5na7eeaZZ+jXrx9Wq5WePXvy5JNPctlllzF79uxGx546dQqr1conn3wCgM1m4/777ycjIwOr1Ur//v1Zvnx5i8+1YcMGJkyYQGRkJBkZGcyZM4eampoO/fqEEB1PgiEhRFh78MEHeeaZZ3j44Yf55ptveOutt0hJSeHOO+/krbfewmaz+Y598803SU9PZ+LEiQDcfvvtvPPOO7z44ovs2bOHl19+mZiYmGafZ9euXUyaNIkbb7yRnTt38u6777J+/fomAZcQIvzI1HohRNiqqqoiOTmZpUuXcueddzb6mM1mIz09nWXLlvGTn/wEgBEjRnD99dfz6KOPsm/fPi644ALWrVvHFVdc0eTcubm5TJw4kbKyMhISErj99tuJjIzklVde8R2zfv16srOzqampISIiomO/WCFEh5HMkBAibO3Zswebzcbll1/e5GNWq5X/+Z//4f/+7/8AyM/PZ8eOHUyfPt1332g0kp2dHdRzbd26lddee42YmBjfbdKkSbjdbgoKCjT7moQQ554p1AsQQoi2ioyM9PvxO++8k+HDh3P06FH+7//+j8svv5xevXoF9bnf53a7mTlzJnPmzGnysZ49e7bqXEIIfZHMkBAibPXv35/IyEj++9//NvvxIUOGMHr0aP785z/z1ltv8fOf/7zRx9xuN3l5eUE918iRI9m9ezf9+vVrcrNYLJp8PUKI0JBgSAgRtiIiIliwYAH3338/K1eu5ODBg2zatKnRjrA777yTp59+GpfLxQ033OB7PDMzkzvuuIOf//znvPfeexQUFJCbm8tf//rXZp9rwYIFbNy4kVmzZpGfn8/+/ft5//33ufvuuzv86xRCdCwJhoQQYe3hhx/mN7/5DY888giDBg3i5ptv5uTJk76P33rrrZhMJqZNm9akyHnZsmX8+Mc/5te//jUDBw5kxowZLW6VHzp0KHl5eezfv5/x48czYsQIHn74YdLS0jr06xNCdDzZTSaEOK8VFhaSmZnJ5s2bGTlyZKiXI4TQIQmGhBDnJYfDwfHjx3nggQc4fPgwX3zxRaiXJITQKblMJoQ4L33xxRf06tWLrVu38vLLL4d6OUIIHZPMkBBCCCE6NckMCSGEEKJTk2BICCGEEJ2aBENCCCGE6NQkGBJCCCFEpybBkBBCCCE6NQmGhBBCCNGpSTAkhBBCiE5NgiEhhBBCdGr/H1jScaN08IatAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sns.lineplot(data=df, x=\"cycle\", y=\"fluorescence\", hue=\"well\");\n", "plt.savefig(\"pcr_results.png\")" ] }, { "cell_type": "code", "execution_count": 9, "id": "841e114f-dd6d-4136-a093-8070a3069643", "metadata": {}, "outputs": [], "source": [ "df.to_csv('pcr_sample_1.csv', index=False)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.10" } }, "nbformat": 4, "nbformat_minor": 5 }