{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Getting started\n", "\n", "The following tutorial demonstrates how to use the package to calculate various codon usage bias indices given the coding sequences of genes." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import codonbias as cb" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Load sequences\n", "\n", "Loading sequences is only a matter of converting them to a list of Python strings. For example, the following code loads *E. coli* genes from a FASTA file using Biopython. You may use any other package, or any other source for the sequences, as long as they are converted to a list of strings." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from Bio.SeqIO import parse\n", "\n", "fasta = parse('Escherichia_coli_str_k_12_substr_mg1655.ASM584v2.cds.all.fa', 'fasta')\n", "genes = [str(record.seq) for record in fasta]" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Initialize a codon bias model\n", "\n", "Next, every codon model (also called *score* in the package) is a class that can be initialized with various parameters, or used with its default parameters (see the documentation of each score). For example, here we initialize the Effective Number of Codons (ENC) score. By default, the variant proposed by Sun, Yang & Xia (MBE, 2013) is used." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "enc = cb.scores.EffectiveNumberOfCodons()" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, we could add background-correction to the model, as proposed by (Novembre, MBE 2002)." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "enc = cb.scores.EffectiveNumberOfCodons(bg_correction=True)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "The rest of the models in the `codonbias.scores` module are initialized in a similar way." ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Compute gene scores\n", "\n", "Once the model has been initialized, we can provide ENC with the *E. coli* sequences and retrieve a score for each gene." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "gene_scores = enc.get_score(genes)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'background-corrected ENC scores of E. coli genes')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHFCAYAAAAUpjivAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABA1klEQVR4nO3deVgW9f7/8dctyw0oELjAjSKauSWoHS2VSlwxPFpppWmZnmx1OZH6tdRMKpSy0jxWZpum5pInszIjNZcWpWOamaZmJyxM0DRlMUWFz++PftynWxYBWcfn47rmurhnPvfMe+Yehhcz85nbZowxAgAAsKgalV0AAABAeSLsAAAASyPsAAAASyPsAAAASyPsAAAASyPsAAAASyPsAAAASyPsAAAASyPsAAAASyPsVDNxcXGy2Ww6evToJbHcqmLYsGFq1KhRZZdR4Q4dOqS4uDjt2LGjzOc9f/582Ww2HThwoMh2efteYcNf35837umnny50eV9//XW+aZ9//rkGDBig+vXry9PTU/7+/oqMjNScOXN08uTJi13VS86ZM2f0wAMPyOFwyM3NTW3bti207bBhw4r8fCtbQfvppXo8qM7cK7sAAFXXoUOH9MQTT6hRo0ZF/sGqCImJifL398833uFw5Bv39NNP67777lNgYOAF5ztlyhQ9+eSTioyM1FNPPaUmTZrojz/+0ObNmxUXF6cffvhBM2fOLJN1uFTMmTNHc+fO1ezZs9WuXTvVqlWryPbe3t5av359BVV38SZPnqyHHnqosstACRB2UGWcPXtWNptN7u7slkX5448/5OPjk298Tk6Ozp07J7vdXglVlb927dqpTp06F2zXo0cPbdy4UVOnTtXzzz9fZNvly5frySef1PDhw/Xaa6+5nEmIiYnR+PHjtWXLlouuvbwVtk9Ull27dsnb21ujRo0qVvsaNWqoY8eO5VxV2WnSpElll4AS4jJWNZWSkqL+/fvLz89P/v7+uvPOO/Xbb7+5tFm2bJmio6PlcDjk7e2tli1b6tFHHy3wtPxXX32lvn37qnbt2vLy8lKTJk0UGxtbZA179+7V5Zdfrg4dOujIkSOSJGOMpk2bprCwMHl5eal9+/Zau3atunTpoi5dujjfu3HjRtlsNi1cuFBjx45V/fr1Zbfb9eOPP0qS3nzzTbVp00ZeXl4KDAxUv379tGfPHpflnz/PPOefYj5w4IBsNpuee+45zZgxQ40bN1atWrXUqVMnJSUl5Xv//Pnz1bx5c9ntdrVs2VILFiwocjucLzc3V7Nnz1bbtm3l7e2tyy67TB07dtQHH3zg0mb69Olq0aKF7Ha76tWrp7vuuksHDx7Mt47h4eH67LPPFBkZKR8fH919993OdZo+fbri4+PVuHFj2e12bdiwQZL09ddf68Ybb1RgYKC8vLx01VVX6Z133slX66+//qr77rtPoaGh8vT0VEhIiG699VYdPnxYGzdu1NVXXy1J+sc//uG8rBAXF+d8f3GXk5SUpGuvvVZeXl4KCQnRhAkTdPbs2RJt1+Jq3ry5hg8frpdeekk///xzkW2ffPJJBQQE6F//+leBl0x8fX0VHR1d5Dy++eYb9enTR/Xq1ZPdbldISIj+/ve/u3yWFbFPSFJGRobGjRunxo0by9PTU/Xr11dsbGy+3/nly5erQ4cO8vf3l4+Pjy6//HLnPIpy+vRpTZgwwWX+I0eO1IkTJ5xtbDabXn/9dZ06dcq5z8yfP/+C8y5riYmJ6t69u3MdW7ZsqYSEBJc2H3zwgTp16iQfHx/5+vqqZ8+exQq3xb2MVdzjoVT8z85ms2nUqFFauHChWrZsKR8fH7Vp00arVq3Kt/z9+/dr8ODBzn2zZcuWeumll1za5ObmKj4+Xs2bN3fum61bt9asWbMuuH7VikG1MmXKFCPJhIWFmf/7v/8zn3zyiZkxY4apWbOmueqqq8yZM2ecbZ966ikzc+ZM89FHH5mNGzeaV155xTRu3Nh07drVZZ6JiYnGw8PDtG7d2syfP9+sX7/evPnmm+b222/Pt9zffvvNGGPMxo0bTUBAgLnpppvMyZMnne0mTJhgJJn77rvPJCYmmtdee800bNjQOBwOExUV5Wy3YcMGI8nUr1/f3HrrreaDDz4wq1atMseOHTPTpk0zksygQYPMRx99ZBYsWGAuv/xy4+/vb3744QfnPKKiolzmmWfo0KEmLCzM+To5OdlIMo0aNTI33HCDWblypVm5cqWJiIgwAQEB5sSJE8628+bNM5LMTTfdZD788EOzaNEic8UVV5jQ0FCXeRZlyJAhxmazmXvuuce8//775uOPPzZTp041s2bNcra57777jCQzatQok5iYaF555RVTt25dExoa6tzGeesYGBhoQkNDzezZs82GDRvMpk2bnOtUv35907VrV/Pvf//brFmzxiQnJ5v169cbT09Pc/3115tly5aZxMREM2zYMCPJzJs3zznvgwcPGofDYerUqWNmzJhh1q1bZ5YtW2buvvtus2fPHpOenu7cHo899pjZsmWL2bJli0lJSTHGmGIvZ/fu3cbHx8dceeWVZsmSJeb99983vXr1Mg0bNjSSTHJycpHbM2/fS0tLM2fPnnUZzp0759JWkhk5cqRJTU01Pj4+ZsiQIfk+261btxpjjDl06JCRZAYOHFisz7UgWVlZpnbt2qZ9+/bmnXfeMZs2bTLLli0zDzzwgPn++++d7Spinzh58qRp27aty+c5a9Ys4+/vb7p162Zyc3ONMcZs3rzZ2Gw2c/vtt5vVq1eb9evXm3nz5rlsq4Lk5uaaXr16GXd3dzN58mSzZs0a89xzzzmPPadPnzbGGLNlyxbTu3dv4+3t7dxnjhw5Uuh8hw4damrWrJnvsz179qzJyckp1efy+uuvG5vNZrp06WIWL15s1q1bZ15++WUzYsQIZ5u3337bSDLR0dFm5cqVZtmyZaZdu3bG09PTfP755852efvNX/fT848xhSnu8bC4n50xxnksu+aaa8w777xjVq9ebbp06WLc3d3Nf//7X2e73bt3G39/fxMREWEWLFhg1qxZY8aOHWtq1Khh4uLinO0SEhKMm5ubmTJlivn0009NYmKieeGFF1zaWAFhp5rJO/A//PDDLuPzfnEXLVpU4Ptyc3PN2bNnzaZNm4wk8+233zqnNWnSxDRp0sScOnXqgsv97bffzMKFC42np6f55z//6XIw+v33343dbs/3x2PLli1GUoFhp3Pnzi5tjx8/bry9vU3v3r1dxv/yyy/GbrebwYMHO8eVNOxERES4/HH8z3/+YySZJUuWGGOMycnJMSEhIeZvf/uby8HlwIEDxsPDo1gHt88++8xIMpMmTSq0zZ49e4wklwOvMcZ89dVXRpKZOHGiyzpKMp9++qlL27x1atKkiUvANcaYFi1amKuuusqcPXvWZXyfPn2Mw+FwfmZ333238fDwcPmjfL6tW7fmCy8lXc7AgQONt7e3SUtLc7Y5d+6cadGiRYnCTkFDkyZNXNrmhR1jjJk0aZKpUaOGc18/P+wkJSUZSebRRx8tcvlF+frrr40ks3LlykLbVNQ+kZCQYGrUqOFcvzz//ve/jSSzevVqY4wxzz33nJHkEvKLIzEx0Ugy06dPdxm/bNkyI8m8+uqrznF5AaY4hg4dWujn27179xLVaIwxmZmZxs/Pz1x33XUuv8d/lfe7HhER4XIMy8zMNPXq1TORkZHOcaUNOyU5Hhb3szPmz308KCjIZGRkOMelpaWZGjVqmISEBOe4Xr16mQYNGpj09HSXeY4aNcp4eXmZ33//3Rjz5+9r27Zti1wXK+AyVjV1xx13uLweMGCA3N3dnZcxJOmnn37S4MGDFRwcLDc3N3l4eCgqKkqSnJeEfvjhB/33v//V8OHD5eXldcHlTp06VcOGDdPTTz+tWbNmqUaN/+1CSUlJys7O1oABA1ze07Fjx0JP+d5yyy0ur7ds2aJTp05p2LBhLuNDQ0PVrVs3ffrppxessTB///vf5ebm5nzdunVrSXJe6ti3b58OHTqkwYMHu1zSCAsLU2RkpMu88u6PyRtyc3MlSR9//LEkaeTIkYXWkfcZnb+O11xzjVq2bJlvHQMCAtStW7cC53XjjTfKw8PD+frHH3/U3r17nfvHX2vs3bu3UlNTtW/fPmetXbt2VcuWLQuttTAlWc6GDRvUvXt3BQUFOd/v5uamgQMHlmiZ69at09atW12GlStXFtp+/PjxCgwM1COPPFLi9SuuK664QgEBAXrkkUf0yiuv6Pvvv8/XpqL2iVWrVik8PFxt27Z1+Tx69eolm82mjRs3SpLz0uSAAQP0zjvv6Ndffy3WuubdQHx+jbfddptq1qx5Ub+b3t7e+T7brVu36uWXXy7xvDZv3qyMjAyNGDGi0N5ceb/rQ4YMcTmG1apVS7fccouSkpL0xx9/lHp9pJIdD4v72eXp2rWrfH19na+DgoJUr14957Hs9OnT+vTTT9WvXz/5+Pjk+/08ffq08xL+Nddco2+//VYjRozQJ598ooyMjIta76qKsFNNBQcHu7x2d3dX7dq1dezYMUlSVlaWrr/+en311VeKj4/Xxo0btXXrVq1YsUKSdOrUKUly3ufToEGDYi130aJFql+/vm6//fZ80/KW/dc/ankKGifl70mTN4+CetiEhIQ4p5dG7dq1XV7n3cibty3y5n3+ti1oXPfu3eXh4eEc8u53+O233+Tm5lbgPPKUdB0LalfYtMOHD0uSxo0b51Kfh4eHRowYIUnOxwf89ttvxf7cz1eS5Rw7dqxY2/RC2rRpo/bt27sM4eHhhbb38/PTY489psTERJd/AvI0bNhQkpScnFyiOv7K399fmzZtUtu2bTVx4kS1atVKISEhmjJlivOepIraJw4fPqydO3fm+zx8fX1ljHF+Hp07d9bKlSt17tw53XXXXWrQoIHCw8O1ZMmSItf12LFjcnd3V926dV3G22w2BQcHX9TvZo0aNfJ9tu3bt1ezZs1KPK/iHNMutL1zc3N1/PjxEi+7oGUU53hY3M8uz/nHMunP49lfj2Xnzp3T7Nmz882zd+/ekv73+zlhwgQ999xzSkpKUkxMjGrXrq3u3bsX+IiG6oxuL9VUWlqa6tev73x97tw5HTt2zPlLsH79eh06dEgbN250ns2R5HIjoSTngev8myALk5iYqIEDB+r666/Xp59+qrCwMOe0vGXn/SE8v96Czu6c/59X3jxSU1PztT106JBLbxwvLy+lp6fna1faZwHlLTstLS3ftPPHzZ07V5mZmc7XeXXVrVtXOTk5SktLKzSk/HUdzz8gn7+OUv5tVNS0vPdOmDBB/fv3L/A9zZs3d9Za3M/9fCVZTu3atYu1TcvDgw8+qFmzZumRRx7Rgw8+6DLN4XAoIiJCa9asuajeTBEREVq6dKmMMdq5c6fmz5+vJ598Ut7e3nr00UcrbJ+oU6eOvL299eabbxa4jL/O46abbtJNN92k7OxsJSUlKSEhQYMHD1ajRo3UqVOnQms8d+6cfvvtN5fAY4xRWlqa84xRZSvOMe1Cx5kaNWooICDgouooyfGwJJ9dcQQEBMjNzU1Dhgwp9Ixi48aNJf35j/KYMWM0ZswYnThxQuvWrdPEiRPVq1cvpaSkVKlefheDMzvV1Ntvv+3y+p133tG5c+ecd/jnHQzP74Y8d+5cl9fNmjVTkyZN9Oabbyo7O/uCyw0LC9Pnn38uu92u66+/Xvv373dO69Chg+x2u5YtW+bynqSkpAv2isnTqVMneXt7a9GiRS7jDx48qPXr16t79+7OcY0aNdIPP/zgUvexY8e0efPmYi3rfM2bN5fD4dCSJUtkjHGO//nnn/PNs3nz5i7/geYduGJiYiT9+ZyRwuRdfjh/Hbdu3ao9e/a4rGNp1qFp06b69ttvC/xPuX379s7T3zExMdqwYYPzclNBzj/7VZrldO3aVZ9++qnLQT8nJyffflIePD09FR8fr61bt2r58uX5pk+ePFnHjx/XP//5T5fPPE9WVpbWrFlTrGXZbDa1adNGM2fO1GWXXabt27dLqrh9ok+fPvrvf/+r2rVrF/h5FPTPht1uV1RUlJ555hlJf/YsK0xeDefX+O677+rkyZMXtd+WpcjISPn7++uVV14p8DOV/tx/69evr8WLF7u0OXnypN59911nD62LUZLjYWk+u6L4+Pioa9eu+uabb9S6desC51nQ2aHLLrtMt956q0aOHKnff//9gg/8rE44s1NNrVixQu7u7urZs6d2796tyZMnq02bNs7rw5GRkQoICNADDzygKVOmyMPDQ2+//ba+/fbbfPN66aWX1LdvX3Xs2FEPP/ywGjZsqF9++UWffPJJvlAl/fkf8aZNm9SrVy917txZa9euVXh4uAIDAzVmzBglJCQoICBA/fr108GDB/XEE0/I4XC4XBsvzGWXXabJkydr4sSJuuuuuzRo0CAdO3ZMTzzxhLy8vDRlyhRn2yFDhmju3Lm68847de+99+rYsWOaPn26/Pz8SrVNa9Sooaeeekr33HOP+vXrp3vvvVcnTpxQXFxcsS+5XH/99RoyZIji4+N1+PBh9enTR3a7Xd988418fHw0evRoNW/eXPfdd59mz56tGjVqKCYmRgcOHNDkyZMVGhqqhx9+uFT155k7d65iYmLUq1cvDRs2TPXr19fvv/+uPXv2aPv27c4/+k8++aQ+/vhjde7cWRMnTlRERIROnDihxMREjRkzRi1atFCTJk3k7e2tt99+Wy1btlStWrUUEhKikJCQYi/nscce0wcffKBu3brp8ccfl4+Pj1566aUSP5l427ZtBT5U8MorryzyMx80aJCee+45570zf3Xbbbdp8uTJeuqpp7R3714NHz7c+VDBr776SnPnztXAgQML7X6+atUqvfzyy7r55pt1+eWXyxijFStW6MSJE+rZs6ekitsnYmNj9e6776pz5856+OGH1bp1a+Xm5uqXX37RmjVrNHbsWHXo0EGPP/64Dh48qO7du6tBgwY6ceKEZs2a5XJPX0F69uypXr166ZFHHlFGRoauvfZa7dy5U1OmTNFVV12lIUOGXLDGwuTm5hb4GAhJuuqqq5yh293dXVFRUUXeH1SrVi09//zzuueee9SjRw/de++9CgoK0o8//qhvv/1WL774omrUqKHp06frjjvuUJ8+fXT//fcrOztbzz77rE6cOFHgE7hLqiTHw+J+diUxa9YsXXfddbr++uv14IMPqlGjRsrMzNSPP/6oDz/80HkPVt++fRUeHq727durbt26+vnnn/XCCy8oLCxMTZs2vejtUGVU2q3RKJW8ninbtm0zffv2NbVq1TK+vr5m0KBB5vDhwy5tN2/ebDp16mR8fHxM3bp1zT333GO2b99eYO+aLVu2mJiYGOPv72/sdrtp0qSJS4+v87ueG2PMiRMnzLXXXmsCAwOdvQhyc3NNfHy8adCggfH09DStW7c2q1atMm3atDH9+vVzvjevN9by5csLXM/XX3/dtG7d2nh6ehp/f39z0003md27d+dr99Zbb5mWLVsaLy8vc+WVV5ply5YV2hvr2Wefzfd+SWbKlCn5lt20aVPj6elpmjVrZt58881idzU15s+eHjNnzjTh4eHO+jt16mQ+/PBDlzbPPPOMadasmfHw8DB16tQxd955p7Nbd56oqCjTqlWrfMsoap2MMebbb781AwYMMPXq1TMeHh4mODjYdOvWzbzyyisu7VJSUszdd99tgoODjYeHhwkJCTEDBgxw2ZeWLFliWrRoYTw8PPJtr+Iu58svvzQdO3Y0drvdBAcHm//7v/8zr7766kX3xpJk1q5d62yrv/TG+qs1a9Y425/f48UYYzZt2mRuvfVW43A4jIeHh/Hz8zOdOnUyzz77rEuvl/Pt3bvXDBo0yDRp0sR4e3sbf39/c80115j58+e7tKuIfcKYP7vCP/bYY6Z58+bO5URERJiHH37Y2Rtu1apVJiYmxtSvX994enqaevXqmd69e7t0ty7MqVOnzCOPPGLCwsKMh4eHcTgc5sEHHzTHjx93aVdWvbEkmf379zvb6rxeTEVZvXq1iYqKMjVr1nQ++uCZZ55xabNy5UrToUMH4+XlZWrWrGm6d+9uvvzyS5c2F9P1vLjHQ2OK99nlbYOC9vGwsDAzdOhQl3HJycnm7rvvNvXr1zceHh6mbt26JjIy0sTHxzvbPP/88yYyMtLUqVPHeHp6moYNG5rhw4ebAwcOXHD9qhObMYWc5wPKSHJyslq0aKEpU6Zo4sSJlV0OAFQajoeVg7CDMvXtt99qyZIlioyMlJ+fn/bt26fp06crIyNDu3btKrRXFgBYDcfDqoN7dlCmatasqa+//lpvvPGGTpw4IX9/f3Xp0kVTp07lFxvAJYXjYdXBmR0AAGBpdD0HAACWRtgBAACWRtgBAACWxg3K+vOBVocOHZKvr2+Rj+YHAABVhzFGmZmZCgkJKfLBtYQd/fldKKGhoZVdBgAAKIWUlJQiv/yVsCM5v8MnJSWl1F81AAAAKlZGRoZCQ0Odf8cLQ9jR/74008/Pj7ADAEA1c6FbULhBGQAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWBphBwAAWJp7ZRcAAKj6WrVuq9TU1CLbOBwO7d65o2IKAkqAsAMAuKDU1FRFT11ZZJs1k26ukFqAkuIyFgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDS+LgIAUCYyMrMUWDeoyDZ8fxYqA2EHAFAmTG4u35+FKomwAwCXuOJ8o3lGZmYFVQOUPcIOAFziivON5stHdauYYoBywA3KAADA0gg7AADA0gg7AADA0rhnBwBQYeiejspA2AEAVBi6p6MycBkLAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYWqWGnYSEBF199dXy9fVVvXr1dPPNN2vfvn0ubYYNGyabzeYydOzY0aVNdna2Ro8erTp16qhmzZq68cYbdfDgwYpcFQAAUEVVatjZtGmTRo4cqaSkJK1du1bnzp1TdHS0Tp486dLuhhtuUGpqqnNYvXq1y/TY2Fi99957Wrp0qb744gtlZWWpT58+ysnJqcjVAQAAVVClfhFoYmKiy+t58+apXr162rZtmzp37uwcb7fbFRwcXOA80tPT9cYbb2jhwoXq0aOHJGnRokUKDQ3VunXr1KtXr/JbAQAAUOVVqXt20tPTJUmBgYEu4zdu3Kh69eqpWbNmuvfee3XkyBHntG3btuns2bOKjo52jgsJCVF4eLg2b95c4HKys7OVkZHhMgAAAGuq1DM7f2WM0ZgxY3TdddcpPDzcOT4mJka33XabwsLClJycrMmTJ6tbt27atm2b7Ha70tLS5OnpqYCAAJf5BQUFKS0trcBlJSQk6IknnijX9QGAqqBV67ZKTU0tsk1GZmYFVQNUjioTdkaNGqWdO3fqiy++cBk/cOBA58/h4eFq3769wsLC9NFHH6l///6Fzs8YI5vNVuC0CRMmaMyYMc7XGRkZCg0Nvcg1AICqJzU1VdFTVxbZZvmobhVTDFBJqsRlrNGjR+uDDz7Qhg0b1KBBgyLbOhwOhYWFaf/+/ZKk4OBgnTlzRsePH3dpd+TIEQUFBRU4D7vdLj8/P5cBAABYU6WGHWOMRo0apRUrVmj9+vVq3LjxBd9z7NgxpaSkyOFwSJLatWsnDw8PrV271tkmNTVVu3btUmRkZLnVDgAAqodKvYw1cuRILV68WO+//758fX2d99j4+/vL29tbWVlZiouL0y233CKHw6EDBw5o4sSJqlOnjvr16+dsO3z4cI0dO1a1a9dWYGCgxo0bp4iICGfvLAAAcOmq1LAzZ84cSVKXLl1cxs+bN0/Dhg2Tm5ubvvvuOy1YsEAnTpyQw+FQ165dtWzZMvn6+jrbz5w5U+7u7howYIBOnTql7t27a/78+XJzc6vI1QEAAFVQpYYdY0yR0729vfXJJ59ccD5eXl6aPXu2Zs+eXValAQAAi6gSNygDAACUF8IOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNPfKLgAAUDqtWrdVampqkW0yMjMrqBqg6iLsAEA1lZqaquipK4tss3xUt4opBqjCuIwFAAAsjbADAAAsjbADAAAsjXt2AABVSkZmlgLrBhXZxuFwaPfOHRVTEKo9wg4AoEoxubkXvPF6zaSbK6QWWAOXsQAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKVVathJSEjQ1VdfLV9fX9WrV08333yz9u3b59LGGKO4uDiFhITI29tbXbp00e7du13aZGdna/To0apTp45q1qypG2+8UQcPHqzIVQEAAFVUpYadTZs2aeTIkUpKStLatWt17tw5RUdH6+TJk84206dP14wZM/Tiiy9q69atCg4OVs+ePZWZmelsExsbq/fee09Lly7VF198oaysLPXp00c5OTmVsVoAAKAKca/MhScmJrq8njdvnurVq6dt27apc+fOMsbohRde0KRJk9S/f39J0ltvvaWgoCAtXrxY999/v9LT0/XGG29o4cKF6tGjhyRp0aJFCg0N1bp169SrV68KXy8AAFB1VKl7dtLT0yVJgYGBkqTk5GSlpaUpOjra2cZutysqKkqbN2+WJG3btk1nz551aRMSEqLw8HBnm/NlZ2crIyPDZQAAANZUZcKOMUZjxozRddddp/DwcElSWlqaJCkoKMilbVBQkHNaWlqaPD09FRAQUGib8yUkJMjf3985hIaGlvXqAACAKqLKhJ1Ro0Zp586dWrJkSb5pNpvN5bUxJt+48xXVZsKECUpPT3cOKSkppS8cAABUaVUi7IwePVoffPCBNmzYoAYNGjjHBwcHS1K+MzRHjhxxnu0JDg7WmTNndPz48ULbnM9ut8vPz89lAAAA1lSpYccYo1GjRmnFihVav369Gjdu7DK9cePGCg4O1tq1a53jzpw5o02bNikyMlKS1K5dO3l4eLi0SU1N1a5du5xtAADApatSe2ONHDlSixcv1vvvvy9fX1/nGRx/f395e3vLZrMpNjZW06ZNU9OmTdW0aVNNmzZNPj4+Gjx4sLPt8OHDNXbsWNWuXVuBgYEaN26cIiIinL2zAADWkpGZpcC6BZ+9z+NwOLR7546KKQhVWqWGnTlz5kiSunTp4jJ+3rx5GjZsmCRp/PjxOnXqlEaMGKHjx4+rQ4cOWrNmjXx9fZ3tZ86cKXd3dw0YMECnTp1S9+7dNX/+fLm5uVXUqgAAKpDJzVX01JVFtlkz6eYKqQVVX6WGHWPMBdvYbDbFxcUpLi6u0DZeXl6aPXu2Zs+eXYbVAQAAK6gSNygDAACUF8IOAACwNMIOAACwNMIOAACwtEq9QRkArKZV67ZKTU0tsg1dooGKRdgBgDKUmppKl2igiuEyFgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDTCDgAAsDS6ngMALCkjM0uBdYOKbMMzjy4NhB0AgCWZ3FyeeQRJhB0AqJKK8yTmjMzMCqoGqN4IOwBQBRXnSczLR3WrmGKAao4blAEAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKURdgAAgKXxnB0AqGDF+RoDHhgIlB3CDgBUsOJ8jQEPDATKDpexAACApRF2AACApRF2AACApRF2AACApRF2AACApRF2AACApRF2AACApRF2AACApRF2AACApV102MnJydGOHTt0/PjxsqgHAACgTJU47MTGxuqNN96Q9GfQiYqK0t/+9jeFhoZq48aNZV0fAADARSlx2Pn3v/+tNm3aSJI+/PBDJScna+/evYqNjdWkSZPKvEAAAICLUeKwc/ToUQUHB0uSVq9erdtuu03NmjXT8OHD9d1335V5gQAAABejxGEnKChI33//vXJycpSYmKgePXpIkv744w+5ubmVeYEAAAAXw72kb/jHP/6hAQMGyOFwyGazqWfPnpKkr776Si1atCjzAgEAAC5GicNOXFycwsPDlZKSottuu012u12S5ObmpkcffbTMCwQAALgYJQ47knTrrbdKkk6fPu0cN3To0LKpCAAAoAyV+J6dnJwcPfXUU6pfv75q1aqln376SZI0efJkZ5d0AACAqqLEYWfq1KmaP3++pk+fLk9PT+f4iIgIvf7662VaHAAAwMUqcdhZsGCBXn31Vd1xxx0uva9at26tvXv3lmlxAAAAF6vEYefXX3/VFVdckW98bm6uzp49WyZFAQAAlJUSh51WrVrp888/zzd++fLluuqqq8qkKAAAgLJS4t5YU6ZM0ZAhQ/Trr78qNzdXK1as0L59+7RgwQKtWrWqPGoEAAAotRKf2enbt6+WLVum1atXy2az6fHHH9eePXv04YcfOh8wCAAAUFWU6jk7vXr1Uq9evcq6FgAAgDJXqrAjSWfOnNGRI0eUm5vrMr5hw4YXXRQAAEBZKXHY2b9/v+6++25t3rzZZbwxRjabTTk5OWVWHAAAwMUqcdgZNmyY3N3dtWrVKueXgQIAAFRVJb5BeceOHZo7d65iYmLUtm1btWnTxmUoic8++0x9+/ZVSEiIbDabVq5c6TJ92LBhstlsLkPHjh1d2mRnZ2v06NGqU6eOatasqRtvvFEHDx4s6WoBAACLKnHYufLKK3X06NEyWfjJkyfVpk0bvfjii4W2ueGGG5SamuocVq9e7TI9NjZW7733npYuXaovvvhCWVlZ6tOnD5fTAACApFJcxnrmmWc0fvx4TZs2TREREfLw8HCZ7ufnV+x5xcTEKCYmpsg2drtdwcHBBU5LT0/XG2+8oYULF6pHjx6SpEWLFik0NFTr1q2jxxgAACh52MkLFd27d3cZX143KG/cuFH16tXTZZddpqioKE2dOlX16tWTJG3btk1nz55VdHS0s31ISIjCw8O1efPmQsNOdna2srOzna8zMjLKtGYAAFB1lDjsbNiwoTzqKFBMTIxuu+02hYWFKTk5WZMnT1a3bt20bds22e12paWlydPTUwEBAS7vCwoKUlpaWqHzTUhI0BNPPFHe5QMAgCqgxGEnKiqqPOoo0MCBA50/h4eHq3379goLC9NHH32k/v37F/q+vLNMhZkwYYLGjBnjfJ2RkaHQ0NCyKRoAAFQpJb5BWZI+//xz3XnnnYqMjNSvv/4qSVq4cKG++OKLMi3ufA6HQ2FhYdq/f78kKTg4WGfOnNHx48dd2h05ckRBQUGFzsdut8vPz89lAAAA1lTisPPuu++qV69e8vb21vbt2533vmRmZmratGllXuBfHTt2TCkpKXI4HJKkdu3aycPDQ2vXrnW2SU1N1a5duxQZGVmutQAAgOqhxGEnPj5er7zyil577TWXnliRkZHavn17ieaVlZWlHTt2aMeOHZKk5ORk7dixQ7/88ouysrI0btw4bdmyRQcOHNDGjRvVt29f1alTR/369ZMk+fv7a/jw4Ro7dqw+/fRTffPNN7rzzjsVERHhvJEaAABc2kp8z86+ffvUuXPnfOP9/Px04sSJEs3r66+/VteuXZ2v8+6jGTp0qObMmaPvvvtOCxYs0IkTJ+RwONS1a1ctW7ZMvr6+zvfMnDlT7u7uGjBggE6dOqXu3btr/vz5cnNzK+mqAQAACypx2HE4HPrxxx/VqFEjl/FffPGFLr/88hLNq0uXLjLGFDr9k08+ueA8vLy8NHv2bM2ePbtEywYAAJeGEl/Guv/++/XQQw/pq6++ks1m06FDh/T2229r3LhxGjFiRHnUCAAAUGolPrMzfvx4paenq2vXrjp9+rQ6d+4su92ucePGadSoUeVRIwAAQKmVOOxI0tSpUzVp0iR9//33ys3N1ZVXXqlatWqVdW0AAAAXrVRhR5J8fHzUvn37sqwFAACgzJU47PTr16/ApxPbbDZ5eXnpiiuu0ODBg9W8efMyKRAAAOBilPgGZX9/f61fv17bt293hp5vvvlG69ev17lz57Rs2TK1adNGX375ZZkXCwAAUFIlPrMTHByswYMH68UXX1SNGn9mpdzcXD300EPy9fXV0qVL9cADD+iRRx4p96+PAAAAuJASn9l54403FBsb6ww6klSjRg2NHj1ar776qmw2m0aNGqVdu3aVaaEAAAClUeKwc+7cOe3duzff+L179yonJ0fSnw/6K+pbxwEAACpKiS9jDRkyRMOHD9fEiRN19dVXy2az6T//+Y+mTZumu+66S5K0adMmtWrVqsyLBQAAKKkSh52ZM2cqKChI06dP1+HDhyVJQUFBevjhh/XII49IkqKjo3XDDTeUbaUAAAClUOKw4+bmpkmTJmnSpEnKyMiQ9OeXgP5Vw4YNy6Y6AADKUUZmlgLrBhXZxuFwaPfOHRVTEMpFqR8qKOUPOQBgZa1at1VqamqRbTIyMyuoGpQFk5ur6Kkri2yzZtLNFVILys9FhR0AuJSkpqZe8A/j8lHdKqYYAMVW4t5YAAAA1QlhBwAAWBphBwAAWFqpws6oUaP0+++/l3UtAAAAZa7YYefgwYPOnxcvXqysrCxJUkREhFJSUsq+MgAAgDJQ7N5YLVq0UO3atXXttdfq9OnTSklJUcOGDXXgwAGdPXu2PGsEAAAotWKf2UlPT9fy5cvVrl075ebmqnfv3mrWrJmys7P1ySefKC0trTzrBAAAKJVih52zZ8/qmmuu0dixY+Xt7a1vvvlG8+bNk5ubm9588001adJEzZs3L89aAQAASqzYl7H8/Px01VVX6dprr9WZM2f0xx9/6Nprr5W7u7uWLVumBg0a6D//+U951goAAFBixT6zc+jQIT322GOy2+06d+6c2rdvr+uvv15nzpzR9u3bZbPZdN1115VnrQAAACVW7LBTp04d9e3bVwkJCfLx8dHWrVs1evRo2Ww2jRs3Tn5+foqKiirPWgGg3LRq3VaBdYOKHPjeK6B6KvV3Y/n7+2vAgAEaPny41q9fLx8fH23atKksawOACsP3XgHWVaqws3PnTtWvX1+SFBYWJg8PDwUHB2vgwIFlWhwAAMDFKlXYCQ0Ndf68a9euMisGAACgrPHdWAAAwNIIOwAAwNIIOwAAwNIIOwAAwNJK3fUcAKqLVq3bKjU1tcg2PEMHsC7CDgDL4xk6wKWNy1gAAMDSCDsAAMDSCDsAAMDSCDsAAMDSCDsAAMDSCDsAAMDS6HoOAEARMjKzFFg3qMg2DodDu3fuqJiCUGKEHQAAimBycy/4nKY1k26ukFpQOlzGAgAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkZvLADVWqvWbZWamlpkm4zMzAqqBkBVRNgBUK2lpqZesFvw8lHdKqYYAFUSl7EAAIClcWYHAICLxFOWqzbCDgAAF4mnLFdtlXoZ67PPPlPfvn0VEhIim82mlStXukw3xiguLk4hISHy9vZWly5dtHv3bpc22dnZGj16tOrUqaOaNWvqxhtv1MGDBytwLQAAQFVWqWHn5MmTatOmjV588cUCp0+fPl0zZszQiy++qK1btyo4OFg9e/ZU5l96VsTGxuq9997T0qVL9cUXXygrK0t9+vRRTk5ORa0GAACowir1MlZMTIxiYmIKnGaM0QsvvKBJkyapf//+kqS33npLQUFBWrx4se6//36lp6frjTfe0MKFC9WjRw9J0qJFixQaGqp169apV69eFbYuAACgaqqyvbGSk5OVlpam6Oho5zi73a6oqCht3rxZkrRt2zadPXvWpU1ISIjCw8OdbQqSnZ2tjIwMlwEAAFhTlQ07aWlpkqSgINe724OCgpzT0tLS5OnpqYCAgELbFCQhIUH+/v7OITQ0tIyrBwAAVUWVDTt5bDaby2tjTL5x57tQmwkTJig9Pd05pKSklEmtAACg6qmyYSc4OFiS8p2hOXLkiPNsT3BwsM6cOaPjx48X2qYgdrtdfn5+LgMAALCmKht2GjdurODgYK1du9Y57syZM9q0aZMiIyMlSe3atZOHh4dLm9TUVO3atcvZBgAAXNoqtTdWVlaWfvzxR+fr5ORk7dixQ4GBgWrYsKFiY2M1bdo0NW3aVE2bNtW0adPk4+OjwYMHS5L8/f01fPhwjR07VrVr11ZgYKDGjRuniIgIZ+8sAABwaavUsPP111+ra9euztdjxoyRJA0dOlTz58/X+PHjderUKY0YMULHjx9Xhw4dtGbNGvn6+jrfM3PmTLm7u2vAgAE6deqUunfvrvnz58vNza3C1wcAAFQ9lRp2unTpImNModNtNpvi4uIUFxdXaBsvLy/Nnj1bs2fPLocKAQBAdVdl79kBAAAoC4QdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgaYQdAABgae6VXQCAS1Or1m2VmppaZJs/Tp2Wj7dXkW0yMjPLsiwAFkTYAVApUlNTFT11ZZFtlo/qpugZiRdsAwBF4TIWAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwtCodduLi4mSz2VyG4OBg53RjjOLi4hQSEiJvb2916dJFu3fvrsSKAQBAVeNe2QVcSKtWrbRu3Trnazc3N+fP06dP14wZMzR//nw1a9ZM8fHx6tmzp/bt2ydfX9/KKBeApFat2yo1NbXINhmZmRVUDYBLXZUPO+7u7i5nc/IYY/TCCy9o0qRJ6t+/vyTprbfeUlBQkBYvXqz777+/oksF8P+lpqYqeurKItssH9WtYooBqoiMzCwF1g0qso3D4dDunTsqpqBLSJUPO/v371dISIjsdrs6dOigadOm6fLLL1dycrLS0tIUHR3tbGu32xUVFaXNmzcXGXays7OVnZ3tfJ2RkVGu6wAAgMnNveA/Af/+Zw8CUTmo0mGnQ4cOWrBggZo1a6bDhw8rPj5ekZGR2r17t9LS0iRJQUGuO0VQUJB+/vnnIuebkJCgJ554otzqBgCgNIoTiNZMurlCarGSKn2DckxMjG655RZFRESoR48e+uijjyT9ebkqj81mc3mPMSbfuPNNmDBB6enpziElJaXsiwcAAFVClQ4756tZs6YiIiK0f/9+5308eWd48hw5ciTf2Z7z2e12+fn5uQwAAMCaqlXYyc7O1p49e+RwONS4cWMFBwdr7dq1zulnzpzRpk2bFBkZWYlVAgCAqqRK37Mzbtw49e3bVw0bNtSRI0cUHx+vjIwMDR06VDabTbGxsZo2bZqaNm2qpk2batq0afLx8dHgwYMru3TAsuhWDqC6qdJh5+DBgxo0aJCOHj2qunXrqmPHjkpKSlJYWJgkafz48Tp16pRGjBih48ePq0OHDlqzZg3P2AHKEd3KAVQ3VTrsLF26tMjpNptNcXFxiouLq5iCAABAtVOt7tkBAAAoKcIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwNMIOAACwtCr93VgAKhbfaA7Aigg7AJz4RnMAVsRlLAAAYGmEHQAAYGlcxgIuEdyPA+BSRdgBLhHcjwPgUsVlLAAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGmEHQAAYGl8EShQxRXn28odDod279xRMQUBQDVD2AGquOJ8W/maSTdXSC0AUB0RdgALyMjMUmDdoAu0yaygagCgaiHsABZgcnMvePZn+ahuFVMMAFQxhB0AAKqR4pzJ5T4+V4QdAACqkeKcyeU+PleEHQAALIazP64IOwAAWAxnf1zxUEEAAGBphB0AAGBphB0AAGBphB0AAGBphB0AAGBp9MYCAOASdCl1TyfsAABwCbqUuqcTdoBSaNW6rVJTU4tsY5X/iACguiPs4JJSnJDyx6nT8vH2KrJNRmambv3Xp0W2scp/RABQ3RF2cElJTU0t1reDR89IvGAbAED1QNiBZRTnrE1GZmYFVQMAqCoIO7CM4p61qSjF6elQ3EtmAIDSI+ygWqiOZ22K09OBS2YAUP4IO6gWqtpZGwC4FFjlWTyEHQAAUCCrPIvHMl8X8fLLL6tx48by8vJSu3bt9Pnnn1d2SdCfl58C6wYVObRq3bayywQAWJglzuwsW7ZMsbGxevnll3Xttddq7ty5iomJ0ffff6+GDRtWdnkXZOUH1BXn8tO//9njgqdJq9r9OACAP1WHS12WCDszZszQ8OHDdc8990iSXnjhBX3yySeaM2eOEhISKrW24t5YW1YPqCurh+YVp01Z7bzFvZEXAFD1VIdLXdU+7Jw5c0bbtm3To48+6jI+OjpamzdvrqSq/qesbqwtTnL+s92Fg1NxewBdqE1l77wAABRHtQ87R48eVU5OjoKCXINAUFCQ0tLSCnxPdna2srOzna/T09MlSRkZGWVen8nN1dlTJ4tuY8wF2+Tm5KjrY29fcHnvjetTJssrVpvc3Atus7Jaf9rQhja0oU01blOMvxelkTdPY0zRDU019+uvvxpJZvPmzS7j4+PjTfPmzQt8z5QpU4wkBgYGBgYGBgsMKSkpRWaFan9mp06dOnJzc8t3FufIkSP5zvbkmTBhgsaMGeN8nZubq99//121a9eWzWYr13ovJRkZGQoNDVVKSor8/PwquxxLYduWH7Zt+WHblp9LddsaY5SZmamQkJAi21X7sOPp6al27dpp7dq16tevn3P82rVrddNNNxX4HrvdLrvd7jLusssuK88yL2l+fn6X1C9fRWLblh+2bflh25afS3Hb+vv7X7BNtQ87kjRmzBgNGTJE7du3V6dOnfTqq6/ql19+0QMPPFDZpQEAgEpmibAzcOBAHTt2TE8++aRSU1MVHh6u1atXKywsrLJLAwAAlcwSYUeSRowYoREjRlR2GfgLu92uKVOm5LtkiIvHti0/bNvyw7YtP2zbotmMuVB/LQAAgOrLMt+NBQAAUBDCDgAAsDTCDgAAsDTCDgAAsDTCDi7KnDlz1Lp1a+eDrDp16qSPP/7YOX3YsGGy2WwuQ8eOHSux4uorISFBNptNsbGxznHGGMXFxSkkJETe3t7q0qWLdu/eXXlFVlMFbVv23dKJi4vLt92Cg4Od09lnS+9C25Z9tnCW6XqOytGgQQM9/fTTuuKKKyRJb731lm666SZ98803atWqlSTphhtu0Lx585zv8fT0rJRaq7OtW7fq1VdfVevWrV3GT58+XTNmzND8+fPVrFkzxcfHq2fPntq3b598fX0rqdrqpbBtK7HvllarVq20bt0652s3Nzfnz+yzF6eobSuxzxaGMzu4KH379lXv3r3VrFkzNWvWTFOnTlWtWrWUlJTkbGO32xUcHOwcAgMDK7Hi6icrK0t33HGHXnvtNQUEBDjHG2P0wgsvaNKkSerfv7/Cw8P11ltv6Y8//tDixYsrseLqo7Btm4d9t3Tc3d1dtlvdunUlsc+WhcK2bR722YIRdlBmcnJytHTpUp08eVKdOnVyjt+4caPq1aunZs2a6d5779WRI0cqscrqZ+TIkfr73/+uHj16uIxPTk5WWlqaoqOjnePsdruioqK0efPmii6zWips2+Zh3y2d/fv3KyQkRI0bN9btt9+un376SRL7bFkobNvmYZ8tGJexcNG+++47derUSadPn1atWrX03nvv6corr5QkxcTE6LbbblNYWJiSk5M1efJkdevWTdu2beNJn8WwdOlSbd++XVu3bs03LS0tTZIUFBTkMj4oKEg///xzhdRXnRW1bSX23dLq0KGDFixYoGbNmunw4cOKj49XZGSkdu/ezT57kYratrVr12afLYoBLlJ2drbZv3+/2bp1q3n00UdNnTp1zO7duwtse+jQIePh4WHefffdCq6y+vnll19MvXr1zI4dO5zjoqKizEMPPWSMMebLL780ksyhQ4dc3nfPPfeYXr16VWSp1c6Ftm1B2HdLJysrywQFBZnnn3+efbaM/XXbFoR99n+4jIWL5unpqSuuuELt27dXQkKC2rRpo1mzZhXY1uFwKCwsTPv376/gKqufbdu26ciRI2rXrp3c3d3l7u6uTZs26V//+pfc3d2d/x3n/bec58iRI/n+c4arC23bnJycfO9h3y2dmjVrKiIiQvv373f2HGKfLRt/3bYFYZ/9H8IOypwxRtnZ2QVOO3bsmFJSUuRwOCq4quqne/fu+u6777Rjxw7n0L59e91xxx3asWOHLr/8cgUHB2vt2rXO95w5c0abNm1SZGRkJVZe9V1o257fw0Vi3y2t7Oxs7dmzRw6HQ40bN2afLUN/3bYFYZ/9i8o+tYTqbcKECeazzz4zycnJZufOnWbixImmRo0aZs2aNSYzM9OMHTvWbN682SQnJ5sNGzaYTp06mfr165uMjIzKLr1aOv9Sy9NPP238/f3NihUrzHfffWcGDRpkHA4H27cU/rpt2XdLb+zYsWbjxo3mp59+MklJSaZPnz7G19fXHDhwwBjDPnsxitq27LNF4wZlXJTDhw9ryJAhSk1Nlb+/v1q3bq3ExET17NlTp06d0nfffacFCxboxIkTcjgc6tq1q5YtW8bzNMrI+PHjderUKY0YMULHjx9Xhw4dtGbNGrbvRXJzc2PfLaWDBw9q0KBBOnr0qOrWrauOHTsqKSlJYWFhkthnL0ZR25bjbdFsxhhT2UUAAACUF+7ZAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAQAAlkbYAVCtDBs2TDabLd9www03SJIaNWokm82mpKQkl/fFxsaqS5cuLuMyMjI0adIktWjRQl5eXgoODlaPHj20YsUK8bxVwDr4uggA1c4NN9ygefPmuYyz2+3On728vPTII49o06ZNhc7jxIkTuu6665Senq74+HhdffXVzm8/Hz9+vLp166bLLrusvFYBQAUi7ACodux2u4KDgwudfv/992vOnDlavXq1evfuXWCbiRMn6sCBA/rhhx8UEhLiHN+sWTMNGjRIXl5eZV43gMrBZSwAltOoUSM98MADmjBhgnJzc/NNz83N1dKlS3XHHXe4BJ08tWrVkrs7/wsCVkHYAVDtrFq1SrVq1XIZnnrqKZc2jz32mJKTk/X222/ne//Ro0d1/PhxtWjRoqJKBlCJ+NcFQLXTtWtXzZkzx2VcYGCgy+u6detq3LhxevzxxzVw4ECXaXk3H9tstvItFECVQNgBUO3UrFlTV1xxxQXbjRkzRi+//LJefvlll/F169ZVQECA9uzZU14lAqhCuIwFwLJq1aqlyZMna+rUqcrIyHCOr1GjhgYOHKi3335bhw4dyve+kydP6ty5cxVZKoByRNgBUO1kZ2crLS3NZTh69GiBbe+77z75+/tryZIlLuOnTZum0NBQdejQQQsWLND333+v/fv3680331Tbtm2VlZVVEasCoAJwGQtAtZOYmCiHw+Eyrnnz5tq7d2++th4eHnrqqac0ePBgl/EBAQFKSkrS008/rfj4eP38888KCAhQRESEnn32Wfn7+5frOgCoODbDY0IBAICFcRkLAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABYGmEHAABY2v8DlYpRAmorHmYAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import seaborn as sns\n", "\n", "sns.histplot(gene_scores, bins=50)\n", "plt.xlabel('ENC')\n", "plt.ylabel('# genes')\n", "plt.title('background-corrected ENC scores of E. coli genes')" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Compute positional weights\n", "\n", "For many scores, it's possible to calculate a positional score for each codon in the sequence. For example, below is a meta-gene analysis of the first 100 codons of *E. coli* genes using the weights of the Codon Adaptation Index (CAI) model. The model requires a set of reference genes to calculate codon frequencies, and in this case we'll provide the entire set of genes as reference." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "cai = cb.scores.CodonAdaptationIndex(ref_seq=genes)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Slicing the genes, as well as padding the resulting vectors to have constant dimesnions, can be done automatically." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(4140, 99)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "gene_vecs = cai.get_vector(genes, slice=slice(3, 300), pad=True) # skip the START codon\n", "gene_vecs.shape" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Finally, we can visualize the result and confirm the presence of a \"ramp\" in the 5'-end of the gene." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'metagene plot of CAI weights')" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHFCAYAAAAaD0bAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAADlU0lEQVR4nOydebgcZZX/v29tXb3ffcmeEEgCSQiLIrsCAgooLoiiCAjOCM78FDd0QBFUGB2HcRnBQUFGB4UZUEAEIayC7GsgJCEhy01y96X37lrf3x9vVfVW3bfvktwLeT/Pkwdud3UtvVSdOud7vodQSik4HA6Hw+Fw9iGEmd4BDofD4XA4nL0ND4A4HA6Hw+Hsc/AAiMPhcDgczj4HD4A4HA6Hw+Hsc/AAiMPhcDgczj4HD4A4HA6Hw+Hsc/AAiMPhcDgczj4HD4A4HA6Hw+Hsc/AAiMPhcDgczj4HD4A4nBnimmuuwV133TXTuzHrOf/887Fo0aJJvfb666/HLbfcMq37AwC6ruMLX/gCuru7IYoi1qxZM+5r/vznP+OMM85AZ2cnFEVBS0sLTjzxRNx6660wDKNq+eHhYQQCARBC8MILL/iu8/zzz0ckEpnq4UyZqXxG3/3ud0EIwfDw8LjL8t8MZzrhARCHM0Pwk/meZ08FQDfccAP+67/+C5dffjmefPJJ/O53v6u5LKUUF1xwAT70oQ/Btm1cd911eOihh/Df//3fOPjgg3HJJZfg+uuvr3rd7373O+i6DgC46aabpv0YppNvf/vb+NOf/rTHt8N/M5zpRJrpHeBwOJy3G6+//jqCwSD+6Z/+adxl/+3f/g233HILrrrqKnznO98pe+6MM87AN77xDWzZsqXqdTfffDM6OjqwcOFC/OEPf8B1112HYDA4bccwney3334zvQsczoThGSAOpwQ3Hb9u3TqcddZZiMfjaGlpwVe+8hWYpolNmzbh1FNPRTQaxaJFi/CjH/2oah2pVApf+9rXsHjxYiiKgrlz5+LLX/4ystmstwwhBNlsFv/93/8NQggIIXjve98LABgaGsIll1yCAw88EJFIBB0dHTjhhBPwxBNPVG1r165d+PjHP45oNIqmpiZ8+tOfxvPPPw9CSFXm44UXXsCHPvQhtLS0QFVVHHLIIfjf//3fsmVuueUWEELw6KOP4uKLL0ZbWxtaW1vx0Y9+FL29vVXbv/3223HkkUciHA4jEonglFNOwcsvvzzu++xuZ+3atbjgggvQ0tKCcDiMM844A1u3bh339YVCAd/61rfK3uMvfvGLSCQS3jKLFi3C+vXr8fjjj3vv8XhlmkbWSwjBr3/9a+TzeW+9tbJMhmHghz/8IZYvX45vf/vbvst0dXXhmGOOKXvs2Wefxeuvv45zzz0Xn//855FMJnHnnXeO+7748Ze//AWEEDz//PPeY3feeScIITjttNPKll29ejU+9rGPeX9TSnH99ddjzZo1CAaDaG5uxsc//vGqz8ivBJZIJHDhhReipaUFkUgEp512GrZu3QpCCL773e9W7efAwAA+9alPIR6Po7OzE5/73OeQTCa95+v9ZnK5nPebU1UVLS0tOPzww/GHP/xhUu8ZZx+BcjgcjyuvvJICoMuWLaPf+9736Nq1a+k3vvENCoD+0z/9E12+fDn92c9+RteuXUsvuOACCoDeeeed3uuz2Sxds2YNbWtro9dddx196KGH6E9/+lMaj8fpCSecQG3bppRS+vTTT9NgMEg/+MEP0qeffpo+/fTTdP369ZRSSjdu3Egvvvhietttt9HHHnuM3nvvvfTCCy+kgiDQRx991NtWJpOhS5cupS0tLfQXv/gFfeCBB+ill15KFy9eTAHQ3/zmN96yjzzyCFUUhR577LH09ttvp3/961/p+eefX7Xcb37zGwqALlmyhP7zP/8zfeCBB+ivf/1r2tzcTN/3vveVvVc/+MEPKCGEfu5zn6P33nsv/eMf/0iPPPJIGg6HvWOphbud+fPn08997nP0/vvvpzfeeCPt6Oig8+fPp2NjY96y5513Hl24cKH3t23b9JRTTqGSJNFvf/vb9MEHH6Q//vGPaTgcpocccggtFAqUUkpfeuklumTJEnrIIYd47/FLL71Uc58aXe/TTz9NP/jBD9JgMOitd3Bw0HedTz31FAVAL7vssrrvRyWf//znKQC6fv16mkqlaCgUou9973urljvvvPNoOByuu650Ok1lWabXXHON99gXvvAFGgwGaTgcprquU0opHRgYoIQQev3115fthyzL9Ktf/Sr961//Sn//+9/T5cuX087OTtrf31+2H6WfkWVZ9JhjjqGqqtJ//dd/pQ8++CC96qqr6P77708B0CuvvNJbtvQ3953vfIeuXbuWXnfddTQQCNALLrjAW67eb+Yf//EfaSgUotdddx199NFH6b333kv/9V//lf785z9v7A3n7JPwAIjDKcE9Gf/7v/972eNr1qyhAOgf//hH7zHDMGh7ezv96Ec/6j127bXXUkEQ6PPPP1/2+jvuuIMCoPfdd5/3WDgcpuedd964+2SaJjUMg5544on0Ix/5iPf4L37xCwqA3n///WXL/+M//mNVYLN8+XJ6yCGHUMMwypY9/fTTaXd3N7Usi1JaDEwuueSSsuV+9KMfUQC0r6+PUkppT08PlSSJ/vM//3PZcul0mnZ1ddFPfOITdY/J3U7p8VBK6d///ncKgH7/+9/3Hqu8uP71r3+lAOiPfvSjstfefvvtFAC98cYbvccOOuggevzxx9fdl8mst5HAg1JKb7vtNgqA/vKXv2xoHyhlQXQsFqPvec97yrZHCKFbtmwpW7bR/TjmmGPoCSec4P29dOlS+vWvf50KgkAff/xxSimlt956KwVA33zzTUopCzj8fgs7d+6kwWCQfuMb3yjbj9LP6C9/+QsFQG+44Yay11577bU1A6DK9/2SSy6hqqp6Nw2U1v7NrFy5kp555pnjvg8cTim8BMbh+HD66aeX/b1ixQoQQvCBD3zAe0ySJCxduhQ7duzwHrv33nuxcuVKrFmzBqZpev9OOeUUEELw2GOPNbT9X/7ylzj00EOhqiokSYIsy3j44YexYcMGb5nHH38c0WgUp556atlrP/WpT5X9vWXLFmzcuBGf/vSnAaBsvz74wQ+ir68PmzZtKnvNhz70obK/V69eDQDesT7wwAMwTROf/exny9anqiqOP/74ho/T3SeXo446CgsXLsSjjz5a8zWPPPIIAFZ2KeWss85COBzGww8/3NC299Z6J8r//u//IpVK4XOf+5z32Oc+9zlQSvGb3/xmUus88cQT8fe//x35fB47duzAli1b8MlPfhJr1qzB2rVrAQAPPfQQFixYgP333x8A+y4TQvCZz3ym7DPu6urCwQcfXPczfvzxxwEAn/jEJ8oer/xuluL3nSsUChgcHBz3+N797nfj/vvvxze/+U089thjyOfz476Gw+EBEIfjQ0tLS9nfiqIgFApBVdWqxwuFgvf3wMAA1q1bB1mWy/5Fo1FQShtq9b3uuutw8cUX44gjjsCdd96JZ555Bs8//zxOPfXUshP7yMgIOjs7q15f+djAwAAA4Gtf+1rVfl1yySUAULVfra2tZX8HAgEA8LbvrvNd73pX1Tpvv/32ho4TYPoXv8dGRkZqvmZkZASSJKG9vb3scULIuK+tx55Y74IFCwAA27Zta/g1N910E1RVxamnnopEIoFEIoHVq1dj0aJFuOWWW2BZ1oT346STToKmaXjyySexdu1atLW14ZBDDsFJJ52Ehx56CADw8MMP46STTvJeMzAwAEopOjs7qz7jZ555pu5n7L6Xlb8jv++ry3jfuXr87Gc/w2WXXYa77roL73vf+9DS0oIzzzwTmzdvHve1nH0X3gXG4UwjbW1tCAaDuPnmm2s+Px7/8z//g/e+97244YYbyh5Pp9Nlf7e2tuK5556ren1/f7/vNr/1rW/hox/9qO82ly1bNu5++a3zjjvuwMKFCyf02lIq99V9bOnSpTVf09raCtM0MTQ0VBasUErR39+Pd73rXZPalz2x3sMPPxwtLS24++67ce2114IQUnf5N998E08++SSAYvBUyQMPPIAPfvCDE9qPI444ApFIBA899BC2b9+OE088EYQQnHjiifj3f/93PP/88+jp6SkLgNra2kAIwRNPPOEFI6X4Pebivpejo6NlQZDf5z0dhMNhXHXVVbjqqqswMDDgZYPOOOMMbNy4cY9sk/P2h2eAOJxp5PTTT8dbb72F1tZWHH744VX/SjtlAoGA790tIaTq4rJu3To8/fTTZY8df/zxSKfTuP/++8sev+2228r+XrZsGfbff3+8+uqrvvt0+OGHIxqNTug4TznlFEiShLfeeqvmOhvh1ltvLfv7qaeewo4dO7zuHj9OPPFEACxQLOXOO+9ENpv1ngdqv8dTXW+jyLKMyy67DBs3bsT3vvc932UGBwfx97//HUDR7+dXv/oVHn300bJ/9913H2RZrhlcj7cfxx13HNauXYtHHnkE73//+wEAxx57LCRJwhVXXOEFRC6nn346KKXYvXu37+e7atWqmts7/vjjAbAuwVIqv5sTpZHPs7OzE+effz4+9alPYdOmTcjlclPaJuedC88AcTjTyJe//GXceeedOO6443DppZdi9erVsG0bPT09ePDBB/HVr34VRxxxBABg1apVeOyxx/DnP/8Z3d3diEajWLZsGU4//XR873vfw5VXXonjjz8emzZtwtVXX43FixfDNE1vW+eddx7+4z/+A5/5zGfw/e9/H0uXLsX999+PBx54AAAgCMX7m//6r//CBz7wAZxyyik4//zzMXfuXIyOjmLDhg146aWX8H//938TOs5Fixbh6quvxuWXX46tW7fi1FNPRXNzMwYGBvDcc895d+Tj8cILL+Ciiy7CWWedhZ07d+Lyyy/H3LlzvdKcH+9///txyimn4LLLLkMqlcLRRx+NdevW4corr8QhhxyCc88911t21apVuO2223D77bdjyZIlUFW15oV7IuudCF//+texYcMGXHnllXjuuedwzjnnYP78+Ugmk/jb3/6GG2+8EVdddRWOOOII/Pa3v8WKFStw0UUX+a7rjDPOwD333FOVpWqEE088EV/96lcBwMv0BINBHHXUUXjwwQexevVqdHR0eMsfffTR+Id/+AdccMEFeOGFF3DcccchHA6jr68PTz75JFatWoWLL77Yd1unnnoqjj76aHz1q19FKpXCYYcdhqeffhq//e1vAZR/NydCrd/MEUccgdNPPx2rV69Gc3MzNmzYgN/97nc48sgjEQqFJrUtzj7ATCqwOZzZhtuRMjQ0VPZ4rW6b448/nh500EFlj2UyGXrFFVfQZcuWUUVRaDwep6tWraKXXnppWevwK6+8Qo8++mgaCoUoAK9bSdM0+rWvfY3OnTuXqqpKDz30UHrXXXdVddpQyrqxPvrRj9JIJEKj0Sj92Mc+Ru+77z4KgN59991ly7766qv0E5/4BO3o6KCyLNOuri56wgknlHUoud1ZlV1sjz76KAVQ1oZPKaV33XUXfd/73kdjsRgNBAJ04cKF9OMf/zh96KGH6r7P7nYefPBBeu6559KmpiavxXnz5s1V733lcefzeXrZZZfRhQsXUlmWaXd3N7344ovL2ucppXT79u305JNPptFolAKoWk8lja630e6rUu6++2562mmn0fb2dipJkmct8Mtf/pJqmkbvuusuCoD+5Cc/qbkOt1PN7cyayH68+uqrFADdf//9yx7/wQ9+QAHQr3zlK76vu/nmm+kRRxxBw+EwDQaDdL/99qOf/exn6QsvvOAt4/cZjY6O0gsuuIA2NTXRUChE3//+99NnnnmGAqA//elPveVq/ebc78i2bdu8x2r9Zr75zW/Sww8/nDY3N9NAIECXLFlCL730Ujo8PNzQe8PZNyGUUjoDcReHw9lDXHPNNbjiiivQ09ODefPmzfTu+HLLLbfgggsuwPPPP99wuYzz9uf3v/89Pv3pT+Pvf/87jjrqqJneHc4+Di+BcThvY/7zP/8TALB8+XIYhoFHHnkEP/vZz/CZz3xm1gY/nH2DP/zhD9i9ezdWrVoFQRDwzDPP4N/+7d9w3HHH8eCHMyvgARCH8zYmFArhP/7jP7B9+3ZomoYFCxbgsssuwxVXXDHTu8bZx4lGo7jtttvw/e9/H9lsFt3d3Tj//PPx/e9/f6Z3jcMBAPASGIfD4XA4nH0O3gbP4XA4HA5nn4MHQBwOh8PhcPY5eADE4XA4HA5nn4OLoH2wbRu9vb2IRqPjWtdzOBwOh8OZHVBKkU6nMWfOnHENN3kA5ENvby/mz58/07vB4XA4HA5nEuzcuXNcKxAeAPngzkXauXMnYrHYDO8Nh8PhcDicRkilUpg/f35D8w15AOSDW/aKxWI8AOJwOBwO521GI/IVLoLmcDgcDoezz8EDIA6Hw+FwOPscPADicDgcDoezz8EDIA6Hw+FwOPscPADicDgcDoezz8EDIA6Hw+FwOPscPADicDgcDoezz8EDIA6Hw+FwOPscPADicDgcDoezz8EDIA6Hw+FwOPscPADicDgcDoezz8EDIA6Hw+FwOPscMx4AXX/99Vi8eDFUVcVhhx2GJ554ou7yt956Kw4++GCEQiF0d3fjggsuwMjISNkyd955Jw488EAEAgEceOCB+NOf/rQnD4HD4XA4HM7bjBkNgG6//XZ8+ctfxuWXX46XX34Zxx57LD7wgQ+gp6fHd/knn3wSn/3sZ3HhhRdi/fr1+L//+z88//zzuOiii7xlnn76aZx99tk499xz8eqrr+Lcc8/FJz7xCTz77LN767A4HA6Hw6mLbtqwbTrTu7FPQyilM/YJHHHEETj00ENxww03eI+tWLECZ555Jq699tqq5X/84x/jhhtuwFtvveU99vOf/xw/+tGPsHPnTgDA2WefjVQqhfvvv99b5tRTT0VzczP+8Ic/NLRfqVQK8XgcyWQSsVhssofH4XA4HI4v63cn0RYNoDOmzvSuTDvJvAHbpmgOK3t92xO5fs9YBkjXdbz44os4+eSTyx4/+eST8dRTT/m+5qijjsKuXbtw3333gVKKgYEB3HHHHTjttNO8ZZ5++umqdZ5yyik118nhcDicfQvTsrGxL4W8bs3I9jXTQiJvQDPsGdn+nqYvkcemgTQKxsy8v40yYwHQ8PAwLMtCZ2dn2eOdnZ3o7+/3fc1RRx2FW2+9FWeffTYURUFXVxeamprw85//3Fumv79/QusEAE3TkEqlyv5xOBwO551JIm9gNKvP2AU6r1vI6xayujkj29+TWDbFWM7AULqAXWO5md6dusy4CJoQUvY3pbTqMZc33ngD/+///T985zvfwYsvvoi//vWv2LZtG77whS9Mep0AcO211yIej3v/5s+fP8mj4XA4HM5kGM3qGM3qe2dbGR2pgglzhjQ4Od1CwbSQ1d55AVBON1EwTMSDCnaN5pHMGTO9SzWZsQCora0NoihWZWYGBwerMjgu1157LY4++mh8/etfx+rVq3HKKafg+uuvx80334y+vj4AQFdX14TWCQDf+ta3kEwmvX+unojD4XA4e4eRjIbhdGGPb0czLQxmNFi2DWuGAqBMwYRlUWimDdN6Z5XBcroF3bLRHFJg2Da2j2Rn7H0ejxkLgBRFwWGHHYa1a9eWPb527VocddRRvq/J5XIQhPJdFkURAMvyAMCRRx5Ztc4HH3yw5joBIBAIIBaLlf3jcDgczt4jp1vI7wVNTDJnIFswIAoEpr33gw9KKcZyOsIBCYZlQ3+HBUDpggHBqbi0hgPoT+YxuBcC28kgzeTGv/KVr+Dcc8/F4YcfjiOPPBI33ngjenp6vJLWt771LezevRu//e1vAQBnnHEGPv/5z+OGG27AKaecgr6+Pnz5y1/Gu9/9bsyZMwcA8KUvfQnHHXccfvjDH+LDH/4w7r77bjz00EN48sknZ+w4ORwOh1MbSilyjiB5PMnCVBlKaxAFARalM5KZKBg2CoaFqCohXTChmzZCE2iWGs5oaArKkMQZV7BUQSnFaFZHQGKJCVkUEJQlbB/OojmkQJXFGd7DcmY0ADr77LMxMjKCq6++Gn19fVi5ciXuu+8+LFy4EADQ19dX5gl0/vnnI51O4z//8z/x1a9+FU1NTTjhhBPwwx/+0FvmqKOOwm233YYrrrgC3/72t7Hffvvh9ttvxxFHHLHXj4/D4XA442NY1MvGGBaFIu2ZACivWxjJ6oiqElJO8LG3yekmNNNGU0hBIm9MaB8KhoW3hjJY0BJCdzy4B/dycuR0CzndQlgphhZNIRm9yTx2juawf2d0Bveumhn1AZqtcB8gDofD2XtkNRPPbWOO/u9e3IpwYM/cm/cl83h1ZwJz4kEMZ3R0xgM4aE58j2yrFjtHc1jfm8LcpiB6E3kcOCeG+S2hhl47mtXxwvZRNIcVHDK/adZlgQbTBby8Ywzd8WBZFi+nm8jrFg5b1IyoKu/RfXhb+ABxOBwOhwMAhmXDtADTYv8/EQqG1VApi1KKwZQGRRRBCIEokBnJACXzOmSRBQeCQCbUip/TWefa2F7smJsI2YIJEFJVwgwpEjTLnpH3ux48AOJwOBzOjKKbNmxQ2KATFgVvHkyjPzW+yDajmRhzyl8AIAoE2l6+IFs2RSpvQnU0MoooIKs1HgCl8gYCogBRIOhN5GfdKI3RnO4dWyUEmDHbgVrwAIjD4XA4M0pp0GNYjV8kLZsinTcxmCpgPDVHImegYFqeEFckBJa9d4XQOd0s2wdZJMjrjWWwLJsimTMRkEU0hxQMZ3SM5WZPFqhgWMhqFlTZP6ygdOLZvT0ND4A4HA6HM6Popg0CliXQJ1AS0k0bhmUjkdORqWMqaNsUA6kCgnJRWyQILKjYm63wed2CblIoErv0yqIA3bIaCgzyhoWCaUKVBciO9md3Ij9u4Le3yGomCoZVs9OLEMDgJTAOh8PhcIrkdBOSIEASBOQmHABR5A0biTqOw+mCiWTB8MpfACuB2Xs5A5TVTJTKY2RRgOEYIo5HTjehW9RrMY8HZQxlNCTzs8NpOadboBSeB1AlokBQ4AEQh8PhcDhF8oYNSSSQROL5ATWCZlmwKEVIFjGQKtTUxAxnNJgm9TInACuBmZTuVV3KWF5HoHQfBAKL0obEwTnNBEp2VZVFmBZFXzK/J3Z1woxmdS+z5YckCLNu+CsPgDgcDoczY5hOd5AsCJAFAbrZ+IgKN3CIqjKSBQNpnzJYVjOxO5FHrKL92ssATUBzNBU009XIVJeINHP8oG8sZyBQEWA0BWUMpDSkCzObBdJMC+mCgWAdo0OWAbJmTckO4AEQh8PhcGYQw6IwrWIGyLTshsWymmFDAIEiCTBMioSPKLgvmUdONxFRy72FCCGg2HudSXndgmbYVQGQgPFb4XXTRlYzq14bUiQUDAsDyfpdcFnNxJv9qT02dyynWSj4HFspkkBg2ns34zYePADicDgczoyhWzYMm0ISCCSBwLAbb4XP6qbnqeOWwUqzRxkn+xMP+s+aIMBe0wDlnG4vUSjXyMgiGbcVnk1Y9w8wmoIKdo3lfYM/gPkf7RjJYmcij3Rhz0yfz+omLGpXHVspokBgWTbMvZRxawQeAHE4HA5nxjAsVvKSRAGSKMCyaUPdQu78MNcNOaJKSOXNsnJQfyKPvG4jUsdZem91gaULBgSfAEESBUdAXDswyOlM6+QXYIQDEkybYvtw1jfDM5TR0JsooKBbSOb3TNt8MmdAEuqHE5LANFfGDAygrQUPgDgcDoczY/iVuxrxAtKdUpmbAZJFARa1MeY4JGc0E7uTeTQF649e2BsZIEopEjkDqo9IWBYJDMuu2wmWKZg1u6sAoDWsoD+loa+iFKabNrYPZyEKBNGAjKGMPu3miaZlI5mvr/8B3AwQ5RkgDofD4XAAwDBpWXcTGjTMcz2ASju7QrKEgbQG07LRl8izwZx1sj8C2TvjMPKGVdMjRxYFGLZds+xHKcVYTvcNnlwkUUAkIGH7SLYsA7Z7LIfRrIHmkIJwQEK6YCCjT28ZLKtbyNXx/3HxNFezyAyRB0AcDofDaQjNtBrqWJoIecMsK+2IAkG+gYu0bjI9iVTy2khAQjpvojdRQG8ij+Ya2h8XgZBpPx4/croFzbSrurgApzRk1W6Frxc8lRIPyshpFnaM5GDbFMm8gZ7RHOJBGaJQFIpPtw4oq5mwrHKLgdpQGFwEzeFwOJy3G71jeWzqT09r2Siv22VBjCQQ5Bvwi9EtGxQoG7wpiQJsUAxnNOSN+tkfgAVbE509Nhnyjklg6b66AY/7WK0AKKdbKNQInippjwbQm8ijL1XAjpEsNLNc/ySLAobT2lQOpYqxrF72+dWH8AwQh8PhcN5+ZHXWVdWbmB7zPUop8oZVlj1oRBQM1A4YwoqEvmQezaH62R/ACYDMPZ+RYBPgi8d4y1Pb8clfPY23hjIAWDeaVqMVPq9bQEXwVAtZFKBKInYMZ9GfLKAtEih7PhwQkcjrbJ0N4HbR1cKwbCTyBoJK/eyUi0jGb/nfm/AAiMPhcDgNUTBsUBvYOpypOYLBsmnDF1jdsmHazAPIRRYJTNseVwid1UzfzENMldAeURFSamd/Nvan8PTWkb0yENWdAF+awXnqrWEYFsXru5MAmEtytkbZbyynN1heYjSFZKQ1EyFFqnpdUBaR122kGjRO7EvksX04W7NMmNMs5HVrXAG0iyQQFGaRGzQPgDgcDoczLpbNdCpNIRmaYfu2Xed1C+t7k3hu+wje6E1iOKPVDS5cE8SyDJAgNGSGmNUt38CAEFI3I6GZFq68Zz2uvW8DRnPaHh+ImtNNaCUT4HO66XVrDTrlqFpeQKZlI10wa05Y94MQgjnxIOI+3W+EEIgC8Trl6lEwLAykCmyOWo1gN6ObMKntWRGMhyjwDBCHw+Fw3mYYlg3DtiEJAtojAfQl82XlkWTOwOu9SfQmClBEEb2JAl7uSeClnjHsTuR9tR+GacO0UCaCZm7Q9TvBTKdtfCKZEZdXdyZZiQ1svMRUB6ImcnrdgCJvWNCt4gT4rUNZ77mBFAuE2FBUu6qslzMsFMzxBdATIaJIGM7q43a/jWZ1ZDQTAgESWf8AaCyrQR7H/6cUSRBg7uUBtPXgARCHw+HMQlIFAz0juRmf8+RiWDZMmzojKwTEVBnbh7NI5HQMpApYtzuBZM5Ad1xFJCChM6aiLawgWzDx2q4kBnzEt4Zlg8KGQAhyTglIIAQ26rtB65YN0yx6AE2EZ7aNeP+f08wpD0TtTebRM5qr+XzlENOtwxnv/4sZICcA8smoVWbIpkpQEZHXzbrfK9tmQ1ZVSURYkTCc0aoCUqb/MRGSy0uNybyB3z+7A4Pp6vEcouCWN2dHGYwHQBwOhzNJCobVUDlholBK0TOSw6u7Enhxxxje6E1iNDv9JnYTwbDY4FBXdxNVZRgWxebBNF7bnQS1gc6YWmbYJ4kCWiMBKDW6j3TLBqUEj2wcwNk3PoOHNwyUba8WuukEYxPIPgCsjPfctlHv74xmTmkgqu3oe1IFo2ZpJ1mh/3mrJAM06GWACIyKVnjdZF5GIpney7QoENgUSORqB0DJvIGxrIFYUEY4ICGrm0hVlMGymomCblWVG+94cSf+8PxOfPee9VVaMEkkTsmRZ4A4HA7nbc1oVsfOsdp3/5MllTcxmC6gK6YiKIvYncjjpR1jWLcrgYzPxPO9Abtrp2XdSG2RAMayBiKKhOZw7a6rcEBEMl8dJGgGy/6s28XEwK/3sv8SAHodrYhu2jVHQ9RjY3+qTM+SLphTGoiaNywUdAs5zfL9XAzLRlozESgpYW0dKmaAsjp7XWUrvGVTbBlKoz9V3sllWLaXKZsKIVnEUB191lBag0WZt48oENh2dcCU0UyYdvX8r1edz3LnWB7XP7alrJuPuUHPHjNEHgBxOBzOJGEXpPFbtidKXzIPw6RQZREhRUJXLIiWsIK+ZAHDPqWFvYHhZGtKEQWC7nhwXL+doCwip5tV3UdZ3YQkEgw52SH3v5IgIFcnANJMu9w9ukGe2Tpa9rdrCjhZTQozOLRAQZH2EQrndMfE0MkA6abtlcvcTJqbBaKUPe8OL+0ZyaM9opYFGN+++3V8/rcvYCQzNS+fcEBCRjOR8TFFLBgWBtIFxFS5bPnKgGk0q0MRy7M/qbyBbcMswyUQ4LE3h/DX9f3e8255s5FRJ3sDHgBxOBzOJCkYFhMHT+MJPV0w0J8qoClU3sUjiwJCioThjD7tAVcj6KaNBqxofHG7jyrFtHndgiQQTwvjBUAiQa5OK33BsHwHi9aDUopnHf3PotYQAPZeg05+IGpetwDCRnAMZ6s/l7xuwbKKXVI7RrKwKRBVJSxqCwOAp42SBIKs0yG2dSiL5pDsCacBVpZa35tCqmDi4Y2Dk9pfF1lknXbDmUJVWXUkqyOrmQiXlLZCiohMoagb0k0bKR//n/VOBm9+SwjnHbkIAHDj37Zi80C6bLm9NYB2PHgAxOFwOJMkr9uwLDqtos6BVAEFw/b1sXEvRNkGfXYA1va9uT895dJJTjcnrLkpJSQzMa1b/jAtG4ZNIQgEw05Gww3uZEFgZa4amZmsZkGeYADUM5pDX7IAWSQ4/oAOAECqYIKQyWeAEnmWBQkqIrI+n0tGK58Av9XJjuzXHkFnTAVQqgMSkC4YeHMg7WX+StkyWCydPbRhYMpBcFNQwZbBDDb2p7zSpG1T9CWY+Lm01CmLAixnvAbA9D953a7y/1nn+BqtmhvHRw6Zi/csaYFpU/zrXzcWs00Us2YgKg+AOBwOZxK4LsamPX0TrnM6m2MVV/0nmKuyCM20Gu4Ms22KrUMZvDmYRs/I1LRKBaPcsHCihBQROd3yyk66ZTs+N4anwdEtG6mC6bTC+3cL2Y7R4kQ7o57ZyrI/a+Y3oSPKdDXpglF3IOpwRqtp6uh69AQkgX0ullVVUkrkDARKykSu8/N+7WF0OvtQ2gk2ljVg2/D18Nk8WMyi9CULeKMv1dBx1yIckNAWUdEzksO6XQmMZnUk8gYSOSZ+riQoixhMa7BtysTjqNZgvebof1bPjYMQgi+deAC6YioG0xque2gTKGUB796Yv9YIPADicDicSWBYFJbNupGMaUrpD6YKyBkmwoHavi+SIGA001jn2a6xHHpGcogHFexO5jE6yY411wSx8ZlP1UiiANMuuhAbJgscK8tiQ2kNkkBg2P6t8LplQ7cnEwAx/c97lrR6F/h0gQ1i9bsgMy1ODgMp/1EQOcMqMzgUiIBErvj+FgwLOd1CoMTE0PUAWtIWQYeTAXK9gEIBEZGAhLaIv5h88wALntysy9o3BnyXmwiKJKC7KYh0wcS6nQn0jGZhU//Bpu40+bRmOvqf8mWSeQM7HH3TyrlxAGw47Tc/sBySQPD89jHsHMs7btA8AOJwOJy3LYbFppFbtj0tGaCCYWF3ooCIIted+xQKSBjL6ePeRQ9nNLw1lEFUlREPyqA2sHM0N6lyT6kJ4lQIyhKG0hooZcGNZdsYrgjKhtIFSE7JxfDJzOiO5koWCZJ5A685ZZd6DKU1bBnKgAB496IWRFVWXkoXTIjEfyCq7nRcuVmPSgq65ewHe09CsojRrO6V+PK6Bd20EZBYwGLZFNtGnADIJwMkEIJYsPZn75bAPvmu+QCAJ7cMT0tHmEAIOqIqFEnA7rGCb/YJYMGSbtkYzmhIFYyq8pc71mNhS6hsHfu1RzCvOQiAfQ7iLBqHwQMgDofDmQSuMaBAyLRogIbSGtIFAzG1kY6qYinJj6xmYstABpQSRJ1yWktYwUCq4AmNJ0KpCeJUCCki0gUDWZ2Jx0FIlWHeUEmHk5+4XDdtT1j8y8ffwr/86TU8tqm+KNgVP6/ojqEppHgBUKpggBBAN2mVpqZgMLfptGYi4xNoZCva3t0SnzvSIm9YZa36uxN56KYNVRYwpynoZYBcDVA9RjIaRnM6BAJ8cFU35jYFoZk2ntwyPO5rGyWqyljQEqo7Q00VRYxmNdbZVhEAuYHoqnnxqte1OBYJ7uR4w6Qz6mnlwgMgDofDmQSGRWFRZsY33liB8bBsit2JPEKKNO7Ub1EgoJSNnvDfLxtbBjNIFvSycoosClBEAdtHag+3rEWlCeJkYRom1kHkBo2VAdlQ2skIUf9xGLppw+3Gd3Uxf3pld11RsKv/ec+SFgDwWrxNmzqlzOrxDJphwXQE7pUmgAAwltcRqJhib9o20hpb1tUXubj+P4vbIk7WhWWAXC+gerzpZH8WtISgyiLef2AnAOChaSiDTYRwQEK2YMG2UaX/KRVAV9IcYt/D0ZzO3KCpPW1l46nAAyAOh/OOopE7y+mYReRenEWBoDBFUadmMr+YRqdqBxURwxn/0syOkSz6knl0RNSqYKo5rCCRM9CX8Ne11MLPBHGyyKKAkayGnG5BIkUPoIUtrDXdzQCJAvEt8WiGBQEs6+a+dutQtqYoOFMwvezEEYtbAQABSfDGaOR009edWDNtEEKhiqJXtis+x8wPK7MgAUnEiNPJlsgZnv8PUHSA3s9pf1dl0SsVjZcFctvI9++IAgDet6wDAgE29Kf3iBFnLVRZRNYwq457LKdj52gOBMDKOeNlgARY1vQ1DkwFHgBxOJwZJ10w0J+cuMGfYdkYy+roS+axZTCNl3pG8WLP6LjaiE39KU98OlncAEgSSM1OoYbX5QiCGxX2hhU2nqCyNDOa1dEzmkNzSPGd0C0QgpgqoWc0j4zGppSnCgaG0hp2J/I1j8PPBHGyhBUJiZyBrGZCEgVPA3PgnBgAeCMz2PtanSXIOu34A6kCSmOWe17t9d3es9tGYFMWYM1pYloUQoqlwZxu+Q5ETRcMSIKAcEBCMm+U+RIVdBsFw64KBEIKc7xO5IyqMpGbAVrSHvYec7NAfnPSSnH1P/t3RgCwgOKwhc0AgIc3TM0TaKLMawqhucKjytX/LGoL+3aQuQGQlwGaxs7JqcADIA6HM+MkcsxBdqLdITtGsnhh+xjW7Upi61AW6byF0axRNyAxLRtjOQM9o7kpWfIXDJbBEKdB06BbNuwJjHZQJAGGSct0QIZlY9twBraNujqOqCojp5t4fXcSz20bxfPbR/HSjlHWCp3z7xKbigliJWwYpwXDpBBJsQR2YDcLgIpmiALyRrnLNqWUeQCJBH1OwOzqeZ7ZOlKVScnpJm59rgcAcMz+bWXPuVqrrM9AVEopUgUTAUlkZTvDKnOxzhmm7ygOVRaRNywMZTRopu0ZGVJK8dawGwBFvOUb0QFRSosBkJMBAoD3r2BlsEc2Doyb0czp5rR5VQmEVGUCX6tT/gKKJbCxLAuALDp9nZNTgQdAHA5nxsloJsZyGsZqXID9sG2K4YzORKXxILrjbFwEBWWjEmqgmTYM08ZIRsPIFAaZ5nU2B0kSp65pYFPRJxZAyRUDRneN5jCU1spmR9WiI6rCMG0oooDmoII5TSHIgoBMDX+hqZogluLqYjTLgmbayDtBr5sBGsvpMCw26Z1NDi++L4ajyZFFAX1JVsZbNTeOg+fFYVPg3tf6yrb126d3YCitoTMWwJlr5pY952aAvIGodmmJy4ZeMm1eEgXPrBFgZTUB1RGhe2wFwwIt+XswrSGrMdfrBU6pD0BVJ5gf/akC0poJSSBY2Fp87eGLWhAPyhjLGXhxx1jN128eSONz//08vvK/r+wxB/HxAiA3A+T+3ghmhxkiD4A4HM6MQillZngWRX+y2pq/FlndRF43q+z4BdQvSWmmDd2yoYjipNvCXRNEyRkWOVVNg2ZYID4X1HqUDhhN5HTsGGV+P41kkWRRQFNIQUiRvLJbQBKQzBu+F8n8FE0QKwkpzFPGDXjjQRntztR4CnahlAR3XIPmZQZ1xxxRFgX0JVjWpDsexIcOngMAePCNfm/Z9b1J/MUJiP75fftXlavczFFGcweiFgNYzbDLWtgjAQljWfZeU0qRyBtl/j6lBGWRmTkK1QLoBa2hsjJnpReQH67/z+K2cNlrZVHA+5a1AwB+/9yOsiGvLjvHcvjun9cjq1nYPpLznKink9Gsjl1j+Zr6HwDeoNyxHNNH0RoC970ND4A4HM6MoltMT9ESDmAkq1cNzKxFVrOgW8WLlIsiCnW7atjwSqA5JGM0q09qsKRrgiiLhF2op6hpyBvWhCebq7KIvG5iLKdj+3AWpkURGWcoaT0CEivfVGbPXD+eqXaAlRJTmQvxqNPJ1h4NgBDida0NpTXIIoFlA6/tSuC5baN4dWcCA6kCDIu14/cm3QBIxeGLWtAdV5HVLDyycRCaaeFnD28GAJx8YCcOnt/ksw9FM0T3OF0KpgWLFiedB53SVipvQDOZP5Aq+QvWQ4qIRFYve/4tdwRGW6Rs2UYyQG6n2/6d0arnTls1B2FFxFtDWXz9jlexe6wobh9Ka/jO3euRKimTPr99tGodU8XN/ixuDyNSw8KhxSmBGRYrYRLCAyAOh8NBwWAZmbAiwrJowz416YIBwecUJosCcppVM7NT0C0QsLKGJAjYOTbxLJBrgigKTANk2VMrgeV0C/IES0wCIRAEgoFkAQPpAlrD/g7CjaJIAnSDVmXPpssEsRRCCCIBCUOOB5ArBm53/juU1kAIQWdMRXc8iIAkYCSjY8tABjZlx+6WwObEVQiE4PTV3QCAP6/rxR+e60FvsoCWsILPHb3Ydx9KvYAq51Nphg2UiL4JIRAJwWhWZxPgDbtmBiggiWgJK2XBwFuD1QJooDEN0GZP/xOpeq4rruJHHz8YHdEA+pIFfP2OV/H67iSSeQPfued1DGc0zG0K4rNHLgSwZwOg1TXKXwD7brnBuSuEng1u0DwA4nA4M4pmsGBFEgXEVNkZBlr/5GjbFCNZvar8BbiOtVZNb560ZnqlhOaQjJGMjpHsxLJAnjGgGxQQUjMDZFg2+pL5mvoLb8zEJEpMQVlCX6qAaED27fqaCKJAYDulvVLcY80bJv7yWt+4njUTwQ12qwKgkqwcIQQhRUJ7NIA5TUHMbQrCtGwva+J2dp20ohNBWcSusTzufGk3AOCL790P4RpZMTcDlMqzgah2yeeTKuhVHXkhRcRwVkcyb8CmtMzjp5KoKpd7AJUMQS1lPC8gy6be/DC/AAhg3kA/PutgLOuMIq2Z+Pbdr+OyO9dh11gebREFV3/4IJywjA1/fXMgMyGdXSO8tisBoLb+x6W5ohV+qtYR0wEPgDgczoxSMGxP/RIOiMg4s4bqkdVN5HTT1zdHFgUYFq0534l1EbFTn6vh2T2Wn1AWyDVBLJatak+ET+UN7B7L1xRmuxmWic62AlgWozmoeILeqSIQUuVw7Jog/uW1Pvzy8bfww79unDYxrRvEtEdZJqTdEXAPj5MFHExrsGwKRRK8C2tIkXDSig5vmeP2b8e7Hd8fP4rjMJhhoeYEfpRSZDTL6+ByCSkS8rqJRE6fULYukdMxmtVBACxqLc8AjecFtGss57TbC5jXHKp63qU5pOAHH1mJo/drhemYakZVCVd/aCU6oipaIwEsdYKveoLpiTKS0dCbLEAgwIE19D8uLU7r/IjTCaYZ1e7bexseAHE4nBkloxleIEEIgSqJ6Evm64qhs5oF07KrLlKAm8mwfQMOJoC2ygY5NocUDGe0CWWBKoMdAUJNd2XNtJHVzZoBkDvbajIaG4GQmhmOyeAKoUtxTRDdtvNXdibw4DQ5EA82kAHyw92X7phalmk5ffUcyCJBU1DGPxy3pO46omrlQFT2+WimDc20qoZ9su8Vm/EVaNCwEigaIM5pCvpmLOt5Abnlr/3aI+NqxAKSiG+cuhyfetd8LO2I4LtnHIT5JR1n71rEfIOmswzmlr+WtEXG1Z+VCqElt2w8w51gPADicDgzRqnfikssKGMsa/h2tbgk8zpIndMXpYTpOCrQDBuGScsCJ1kUIICgdyzfsC9QZQBUT9NQMCzkdatmgKSbtlcCnGkCsuAIoYv76pogjpRMoL/579vK2sInS2UJzG3hH08H5up/upvUssfnNAXx808eip+cvabmUE+XWIkGqHQgasGwqr4jLmGFDaIN+DxXizc9F2f/ElY9HdCWOvofPwRCcM4RC/Efn1iDAypE0+9axMaAvNyTmDYBsmuAuHKc8hdQFEKPZotmiNPhyD4VZv4Xx+Fw9ln87rZlUYBFa4uhLZtiLGf43k27SCJBRqsOoDTT8jWwawkr6E8VsHUo21AbvmuC6G2vzoTrdIFlf/wCMmB2dMO4BCQRmmGjUOLA7JogugFPPCgjp1v4z0e3VJUwehN5/OejW/DE5qFxt1UwLC/I7XBLYCUi6HrlES8DFA9WPTe3OYjWBryQSjNAgkC8gagFg02p98u4RFUJTUGlqqW+Hm4AVBmQuOaE9TrBvA6wjuoOsImyX0cEzSEZecPC+l7/sSETZTz/n1JaKjJApj3zZog8AOJwODNGrbvtmCpjIO0vhq6n/3FRRAHpgll1EWVljuoLqyQKaA0HsG04g57R8b1SXBPE4uv93aBtmyKnWaAUNcdzGCb126UZwU8IndNNrwMKAL5y0gGQRYIXd4zhkY1sDINNmUbo/932Mh5Y348fPbAJP35wU13BtFvmCsoiwgH2WboZoLxhIVvHy6k34XaAVQdAjeJqgFy3aYuyjETBMFHL9logpGzUg00pc9+uEaxRSrHJCYCWdZUHMUMZDcm8UdMLyLBsbHXKZ+4IjKkgEILDF7IsUCNlMEopnt02UvNGpFz/Ext3fd44DNcNehbMA+MBEIfDmTE004ZNq++2w4qIrGahP1k9tDOnWczFuE4ZgnWC2V5ZwyWjmZCI/+uYIFXBlsGsV2Lxo9QE0aXWhGtXcxRSxLKxFaXkDdM322DZMyMSJQRlQui8wdya3VERq+bF8al3LwAA/OrJrdjUn8aV96zHLx9/C5ppY3FbGAIBHn9zCF+67WWs7036bmco5QqgA95oBVUWvdJUPSG0lwGqKIE1imVTjGY1T3yf1y1YNut2SxXMKv1PLe57rQ//77ZXcNfLu2vuZ7pgQhYJFrcVBdC6aSOoiHUzQDtGcjBt5u3UFZvccVZSqgMa77v12JtD+P5fNuBHD2z0ff51J4vUiP4HKJkIn9VBCAEhdEqjaKYDHgBxOJwZI6+boD4OyIQQxFUZO0ZyVVqgZF6HME4XjiyyAKhSeJzRTMh1AqdwQEJAEvBmf7qmQWKpCaKL61pceUermRZ0kyKsSCiYtm+5K69Xd4C9OZDGJ/7raXzmpmfx/b+8gTtf2oUNfam9Ui5TRBHJAsv2uCaIrjllU1CGLAr46CHzsLQjgqxm4Wt3vIpXdiagSAL+4dgl+MnZa/DDj61GV0zFYFrDv/zpNfzPMzuqLriVAmiX8YTQlk29bEl3fHKBQVYzEZBFT0CedQaiGpZd1iU4Hs9sHQEA/P2tYd/n3fLXkrZI2To100JQEkFBa2qA3Ncu7YhUzd6aLAfPb4IksDlquxP1g/w/OUHdxv607w3Ba57+Z/zsD1BeAgMACgKDa4A4HM6+Sr277YgqQTNtbB/OemJJdudu1C1/ASwAMit0N4ZlQzOqu3sqaQopsGx24vdzpS41QXRhZojV840KBpvxpcoiDLPaZdm2mfFgZQfYHS/ugm7ZSBVMPLttFLc8tR3fuHMdzr3pWewcy9Xd/6kSkJiRpO4EbIZtI5VnGaEWx6lZFAi+fOL+3n4v64ziZ2cfgjMOngOBECzviuGnn1yDE5d3wKbA7S/sxBOby4MEN8Bp9zIgBQymC+MKoYcyGkybQhZJQ3PP/MgZJqKq7GUucroJizJDSt20GhI5m5aNjf0sSNkymEHGJ8O3qd+//FUclEp8vYA00/ICkEYExo0SUiRPr1OvDPba7iS2lYzNeHJzdYD3+gT0P0AxA1QwmJM2mwfGM0AcDmcfxLbpuHfbbZEA+pN5744/q5vIGyZCdQTQHoSUdTN5M8AauLi1RRRkCib6fO6STYuWmyACnht0ZcnNHbshiwSGRavMGQ3HQbrUBHEwVcCz21hm4RunLMMFRy3Ce5a0IBwQkdUtPL5pfIHxVAhIAgrOkFLXBDGRZ3ftpW7TC1vDuPrDK3HpSQfghx9bjbnN5XqckCLhyycd4A0hfezNwbLnBz0XaJYBcbuCSoXQfrifSVdFC3yjWDaFTdlnEnG0R+kCG4ia0UzoJm0oA7R1OOsFtDYFXtudqFrG0/9UCKB104YqixAJYXPZKryA7nhxF/pTzN37DMfherpwu8Ge317bD+juV3oBFD/vJ7aUB0CjWR27E2z+13j+Py5BRfRuXMayrPNupt2geQDE4XAawrb9zQUnC5u4Xf9uWxYFBBUJ24azyOkmspoJw2zMNFAg5UNRNcNq2G/HHdUwktGr7lJ1y4ZF2eDW83/zHO56Zbf7orKBmgAz2ZMFAYQQUKDq/dNNG2aFCeJ9r/fBpsCa+U04dv92fPTQebj8gwfiH45lvjbP75j+cQalSKIAy7aZQN0xQRxzBNCVGZdVc+M4YXlHXY+a9x/YCYC1X5dmSSpb4Bm0aIZYowRWrwOsEfK6hXCADYF1S2Bu5iWvW6Co7hL0442KTqpXd5VrnXTT9rIoB1RkgECY0J8FxjY6YkUvoN1jedzx4i4AwOePXYKQUq6vyevWlLRhbgC0vjfpK1LvTeS97NA3T10OUSDYNpwtyzy62Z8l7eEJzZ/zhNDeOAyeAeJwOG8Dto9k8fruVM0RExOlYDB9zHgZmaagjHTBxI6RLJI5Y1z9j4siCkiXtMJrJitHNaqnCCkicj4jClwdzks9YxjJ6iXlgfKulkpHYQI2h6x8Xew1blBWMCw8sJ6ZDJ5eced/6IJmEABbh7ITGuD6yMYB3PnSrgleNAlymumZII5m2fs4mXljC1pCWNQagmlTPL21mEko1QDZlEIAAUC8AKiWBsibATZJAXRWN9EWVhAOiF5w4ZY6a3Vz+fFGHwuAljvBzSs7E2XPbx3KwLQpmoKyJ3QG4EyzJ4gH2fgS06JeFmwgVcAv//YWTJvi0AXNOGq/cidrw7IxkM7X9cgaj664ivnNQdjU3xX6nld7QQEcvrAZy7tjWOMMki0tg43X/m5a/nq3ZscNeiyrQ/C5Ydjb8ACIw+GMS6pgYOdoDgOpAnaOTo8GpWBasFF/phLAsjGtYQW7xwoYa0D/46KIAvK67WVw8rrlXGQbQxIFmLZd1b3lntjdC/iIkx2pdINmGa6iWFoWBaR9g6liUPb4m0PIaCY6YwGvZdmlKaR47dAvNDjO4Km3hvEfD23GLU9tx998dBy1CIgCUgWzaILouGS3hienuTl2/3YA8PbBsqkXxLVHAywIFAFRJJ5jcM0S2BQyQDZlAvaWsAJFLLbfu5+xbtoNdYBRSr0A6OzD50MgwO5EvmyfN5X4/5QG3ZphIyCxbQckgXWCORmg+17rY4JyUcAXjl9SFaxrpo2mkFLXIqARjnBGhPzqia3YMVLU+mQKJh7awAJwt3R57NI2AMATm4e8IHq8AGgwrWEoXagaq1LaCj8b4AEQh8Opi21T9IzkUDBttIQU7BzNeSWRqeBOZW8EVRYhEIKxvNaY/geALBEYZrETLKM13t7sbVcSyxyQgaIJoqthGcvpsCmFJJZrGjzNkbNNWRSQ1awyryDdtOH+SSnFveuY9uK0Vd2+ZRg3KGpkntOusRx+8tBm7+9bntrecAkzILGZbHndAiHw3gNXBD1Rjt2fXUTX7UogkdMxktFgU2Yg2RxWvDIg+5tlCYYzmq9TcG9y8h1ged1CKCAhFpQhisUxIumCAUrhBKzjf0d2JVgWRhEFHDy/yTMqfNUZDAqgpv9PwbQQCUiQnBKcYdteBsgN7s46fJ5vgMe0QwIUiUypHP2xQ+dhv/YwknkDl9/1uhcEPfBGPzTTxqLWEFbPY8HNe5a0QhIIdo7lsWMkN67+x7IpBAFY0BpCRjPLslVeK/w0D2SdLDMeAF1//fVYvHgxVFXFYYcdhieeeKLmsueff77jH1D+76CDDipb7ic/+QmWLVuGYDCI+fPn49JLL0WhUG0zzuFwxmc4o6EvWUBbOIBwQIJpU2wfyU65JXsifisA0BpR0BoONNyiXNoK7xoSTnTgaFBhLeGlWiLXBNHNAFk2RTJnVGkaNMNyMhtsm4roTKkved8004LolPRe353E9pEcApKA96/o8t0fV7/xys764wzyuoVr7t+IvGHhoDkxtEUCGM5onrh1PAIyy2alCgYkQfCyXJMpgQEsW7O0IwKbAk+9NeK9d22RAATCBOKSKEASBEQV2Zu7VTm53KbU84bqbpp4BiijmWgJKZ4AOeIFQMyLKW8WS5YvbB/Fjx7Y6FtucvU/B3Sy9vaDnTJRWQDUX1sAHQ+y7YYUEaZNPQ0QAMyJq/jYofN8918zLTQHFbSEFK8zbzJEVAnf+/DKsiBo61DGC8A/vGaul30KByQctpD5Bz2xZdjT/yyuof/J6SZCioT92qNYMScGw7K8jE9LyUT42cCMBkC33347vvzlL+Pyyy/Hyy+/jGOPPRYf+MAH0NPT47v8T3/6U/T19Xn/du7ciZaWFpx11lneMrfeeiu++c1v4sorr8SGDRtw00034fbbb8e3vvWtvXVYHM47Bt20sWM0B0kg3oWhLRLAYLrgufFOBsumyGpmQx1ZLgIhVYLQ8ZZ3hcfeENQJbA9gLsV53fK0RKUmiKXljpEss/fXnWALYBkgQooZDEUSYJi0rDU/q1mQnUzPn9f1AQBOWN6BiOp/nEvaw+OOM6CU4uePbsbO0RxaQgouO2U5zjtyIQDWXdTIxUcWBaZPsiksans6qEZGTNSitJRS6QFkWjaCsghFEmCDeoFWpRnicEbzhOztE9wX1+25zdmmKBCEnWxiyg2A9OJ35PfP9eCJzcOeILkUNwByMyBrnGzJqzsToJRiLKdjMM2MFqtcnAm877G7rc4So8OL37u0ZqBOAYRVCZ1xFYY9NTF0VJXLgqCv/t+rGM7oaArJOP6A9rJljyn57LzyV43ur5xuoTWsQJEEzG0K4qA5cRDCuv6aHA0QzwABuO6663DhhRfioosuwooVK/CTn/wE8+fPxw033OC7fDweR1dXl/fvhRdewNjYGC644AJvmaeffhpHH300zjnnHCxatAgnn3wyPvWpT+GFF17YW4fF4bxj6EvmMZLRvNQ1wC4c0YCM7SPZSYsxmUEga0m3bIq/vt7vlZSmE1d4rJmso2miGSBCCESBIJljx+maIFq2XXbsI1nNG/DoukFXuk6LAoFFi510pcFUaev7aatqtz0LhHh347V8XO55tRdPbB6GKBBc9oHlaA4rOO6AdhzQGUHesPC7Z3c0dvCUtfynnUyDKgtewDAZjnHKYOt7U9jYzwIIt+XdsClCCtPEmKWt8BVCaLdE1BlTvRJhxhNr1ydvWFBlATFnBpgkEO//0wUDIUWEKrFSq00pehyt20MbBqratV39z0HdzARweXcMiiRgLGegZzTnmRjObwmVBe2uAFp13kc3Azo3ruLjh87DhUcv9kTHlVg2E4qrkojmkIKQIk1ZC1QaBLlO3x9c2V31O3n34hYoooC+ZAGPv8lsGFbNqw6A3CCzpSRT2BFTsXJuHJJIEHbei30+A6TrOl588UWcfPLJZY+ffPLJeOqppxpax0033YSTTjoJCxcu9B475phj8OKLL+K5554DAGzduhX33XcfTjvttJrr0TQNqVSq7B+Hs6dI5HT0jGRnZMzBRMhoJnpGc4gG5Co9SlSVoRnMpDCvM7+YRoaIuhQM5n8jiwKe2DyEXzy2BT9+8M3pPgRPeFxr5EYjBGUJI1mduSI7JoiVpZnRrA5JYO3jhtMJ5u86TT1NEusAsyEJxGt9P3heHAtbw6iHqwN6wScAWt+bxM1/3wYAuPDoxTjQuUALhOCiY1gb/UNvDGDbcKa4R5Ti5Z4xrH2jv6wLSpEE5PXiwNLWcGBKjsQdURUrumOgAB52Zoi5GSCbMl8ctyTUXsMMsS9Rrf9J5HRPpF2PjGaiOaR4Q3RFkSASLJbAQorkZWIGU1qZdqx0uOtIRkN/is3AWt7NyluyKHjB0Ku7EjXLX64AOuQI+WVJgCQSWBQ476hFOPOQuTX3XzdtKDJBUBGhyiLaIwHfgb8TxQ2CVnRF0RZR8EGfADykSDjcGaORN5h276Du6gAob1gIymLZvDSAlb5USfQym7NFBN14PnmaGR4ehmVZ6OzsLHu8s7MT/f39476+r68P999/P37/+9+XPf7JT34SQ0NDOOaYY0AphWmauPjii/HNb36z5rquvfZaXHXVVZM7EA5nguway6M/WUBAFstS37MJ07KxczSLrGZiblPId5m2SAB9yTwSeR0iIRCJAFFkFwNRIJDcf6KAiCqhPVK8gBYMCzZlF2bXK2VDXwp9yfyk/V38cIXHBcMCpZO7eIcUEaM53fOwMe1iW7jLSEYvcYNm3V9+rtMiEbzOGDeYUmWChzewgOD01XPG3Z9DFrBxBr1JVoac42hhMpqJHz/4JmwKHH9Ae1Ub/YruGI7dvw1PbB7Gr5/chm+fdiAe2zSEe9b1ep19kYCEI/djmZqAJGAooyGRn3wLfCXHLm3Dhr6ilYIr/iUgXtbBpnZNM8RiCzw7ZstmWT0KJ5isk+EzreJ6AZYBipdkgGxa7EjcUTEQ997X+nDSik4QQrzsz6K2cFl2Z838Jry8M4FXdia846ucAK+ZFprDSpkuTBYEJzNUPx+hWyx4cn2z2qMB7BzLj3vcjRBVZfzwY6tBgZpdmcfu346n3mJZysVtYd8ybVaz0BphGqtKVFlEtGT0iG7aCGPyGcXpYMZF0JV3FJQ25tNxyy23oKmpCWeeeWbZ44899hh+8IMf4Prrr8dLL72EP/7xj7j33nvxve99r+a6vvWtbyGZTHr/du7cOalj4XDGI6uZGHG6WzYPZupOy54JKKUYzmh4dVcCO0byddueZVFAVyyIoCRBEgTYlOlbUnkToxkdfckCekbz2DKYwas7E9jQl/KON6ebXkv6rrGiluixSbocD6QK+L8XduJLt72MK+953Wt9d4XHY1ndc1t+bNMgrrxnfcNeOrLbDq8Z0C0btvMeleKWwNwMkGZavq7Tiih474HrsqwZthdk1Cp/lBJSJG/69guOKSKlFNc/tgXDGQ3dcRVffO9S3/PoeUcugiwSrNuVxHm/eQ6/eGxLma3B5sFiZigSkNAZVb0M0GQ7wEo5ZmkbSpNw7bGAlwmVRRYEEZBxS2BuBsgtKUVVueawWYCJwgOy6JW8gPLJ7jZlQ3ZdekbYe3LYwmYoooCtQ1kvq+Ppf7rLZ2C5QujXd6fw5gB7H6tGYFhFATRQDIAamYquGRZiQcn7XONBGVFVqrJWmCyEkLqWFIcvbIYqs+9zrfEchl0eZJbCutcE76bAdRefSWYsAGpra4MoilXZnsHBwaqsUCWUUtx8880499xzoSjlP8pvf/vbOPfcc3HRRRdh1apV+MhHPoJrrrkG1157LewapkuBQACxWKzsH4ezJxjL6cgbNrrjKrKaibcGM3tlwGUjpAoGNvSl8EpPAomsgc5owPdOrhRRYCn5sNNa3BRS0BJW0BoJoCOqojOmojseRHNIwc6xHF7aMYaekaw3IRtA2VDGxzYNNlwa1E0b973Wh2/cuQ4X/fYF/PaZHdg6nMVLPQlsdy5grvC4UOLvcseLu/BSzxh++3SDWhgAAZG1wxuWDQpWIgHgteR7rfKOuZtmsuxO5V297IyZMJxJ9Ral3kW+KSiP+367vGth+TiDRzcN4YnNwxAI8NX3L/PKPJV0xlR8+GBWZsnpFjqiAVx4zGKc+x4mI9gxUgyGCGFt4m6gOFkPoFKawwpWlohnO6IBmDaFKDKRvSQy4XotM0Q3A+RmCd3MSVdMRc6sHQhkdVb+Cpd0LUkCgSoXMyqlc992OEHhQXNiOO4AlhH7y2tMpL7e1f9UiIAXt4URVSXkDcvTGy1oqc6elmaNBOf308g5wLQpooFiACeJArrjKnL63rmJUmUR71/RCYEUbQ1K0UwLikgQrSHgV2URNuDZHCRyUy/fTZUZC4AURcFhhx2GtWvXlj2+du1aHHXUUXVf+/jjj2PLli248MILq57L5XJVTrGiKIJSOus1F5x3NrZN0ZcsICiLIISgwykh9YxMj7HgVNg1lsPLO8awczSPppCMjpg65bR6KQFJxJx4CCIh2NCfwlhOR0BiJ373oiY6ZZ3SLEQ9bnu+Bzc8/hY29KVAAKyeF/cM5dygyhUeG1ZRcO0+9+imwbKBj/UIKSKSecPL3rhdTK4LsNsm7gqHCwabAVaJIgreUFS3TOLOf5pIOdTVY7y+O4ntw1n88vG3AADnvHtBVdahkk+9ewHOP2oR/uUDy3HjuYfjzDVzscLJZmwfqX4/RrwxGFPPAAFFU0QCVkZ1nbBlUfDKpy0+Zog2pZ4HkOsCbVgUqiKgLRqAKolldgUuhsVGrlS+v6LAMh6lrfAuri/OwpYwTlvFypJPbhnG7rE8tjvfmcoMkEAIDp7X5P29tD1SpjkzLDbzTa0ITsMBcdyp6OzaRasCZFdbs7dmal14zBLcetF7sLyrOkmQ1ayyAbOVSKIAQoAWJ5BOTMHNerqY0RLYV77yFfz617/GzTffjA0bNuDSSy9FT08PvvCFLwBgpanPfvazVa+76aabcMQRR2DlypVVz51xxhm44YYbcNttt2Hbtm1Yu3Ytvv3tb+NDH/oQRHFm642cfZtk3kAyZ3h3SJIooDmkYPtodo90QLlsG8rU7dbSTAs9IzkQEMxpCiIg7bnfSSwoozMahCwICAVE9CcLsClrNz/a0Z48unFwnLUwXnZGD5y+uhu/Of9d+MGZq7wL0K6yiekUupMpGEwXvG4XCuC3T29vaFuqzC4yOc0sM0F0LwSuCFcQmEFdXrd8ywnuUFTNYBoIAmDAySaVesGMx9ymILpiKkyb4vK7XkPesHBgdwwfP2z+uK9VJAEfO3QejtyvzbtAL3QyFYNprSqj4Ga3pkMDBABHL21FRzSANfObvPKiLLDSiCwy3Vg86GpzTDz0xgA2D6TRlyhAN5mQ3dUOGZaNsCIhEpDQFgmUZXFcRrM6OmNqVWmGEOIMRC2aIQJMK+SWZRe2hrC0I4JlnVGYNrMXoGAluBaf96O0hLmsIkhwBdCVTuaqLMKm9TNAhsVGxlRm9iIBCc1hxfe49wSiQGoGOHnDQnu0tlBeFghAKVpCbgZo5ktgMyaCBoCzzz4bIyMjuPrqq9HX14eVK1fivvvu87q6+vr6qjyBkskk7rzzTvz0pz/1XecVV1wBQgiuuOIK7N69G+3t7TjjjDPwgx/8YI8fD4dTD9fZtrQsElIk5HULWweziAbkmqWLyZLXLexK5JHVLcTmxHxPTiMZHWnNQFds+sTH9RAFgianrX6Xk42Z2xTECcs78LfNQ/jb5iFceMziuhmonG5i6xDLFH3s0HmeP808ZyJ5qa5IIgI0w4YcEbzHW0IKkgUDL+wYw2u7ElhVcufuhxvM6Fa5CaKbOXGF1pLjJWPY/iMVikNRbeR0E5IgYKBiKnojEEJw+KJm3LuuD6mCiZAi4ivvP2BSXW4AC0xbQgpGczp6RnNld/jeGIwpeACVElVl3Hju4d6+mhZFJChCcLJAkihAIgRNQRmJvIGfPrK57PUd0YD3WovaXkmpIxZAbyIPyy4OM83rFogALGgN+743iiR4F/SUkwHqS7IgOSiLXtB02upubFqb9ryXVnT7yyRKM0DLKvx/XAF0ZVlUkcYf0KKZFhRJhFqhKSOEoCuuYjBdQEYzawYniRz7jc+Jhyb9HamH4XQzxiu6v0qRRAGUEs9SY58ugblccskl2L59OzRNw4svvojjjjvOe+6WW27BY489VrZ8PB5HLpfD5z//ed/1SZKEK6+8Elu2bEE+n0dPTw9+8YtfoKmpaQ8eBYdTH820MJDWfOvjLWEFyYK+R+7iUgUDWc1Ef6qAMZ8TjmnZ2J3IIyCK487k2hO4mZp5zUGsmd+EpqCMVMH0sju12NCXhk2BrphaNqHc7VgrzQBFVMnzH3EfP3BODKcexNyWf/PU9obK40FZRKpggqLoY7K4LexpSEazrBMsb9jI68UA6InNQ9jYV7TWEAjrgssbrCTi6ok6J5ABAoo6IAC4+Pj9ptxRuLCVvXelOiDLpl7L8nRlgABUlYaCspMVdToHTZvimx9YjtNWdWPV3HjZhXVFSXBGUDQTbA4piKlyWWPBaE7H3HjQN1sDsE63cEUGyPX/md8S9H4TxyxtK9uHg+b4B0BdcRUHdscQD8pVGiHNshHz+f2z7wmpO4hVN21vfEYl7ZEADuiMIqubVVkVSikGUwXYoGgOBvZY00VOtxBRJUTV2gGQKBCIIvHMEGdDCWxGM0Aczr5CImcgUzB95xe5WZnxpqwncwZAUPcuqxLWnSQAlF38m0NyWRZoNKdjLKtPKPswnbgZmXnNQYgCwXEHtOOeV3vx2KZBb+yDH+t7mRtt5YXIzQD1JgpeW3NIkeD6OJZu74OruvHwxgFsHszgqbdGcPTSamFnKeGAhJHRHASw8pkiCYipElrDCnqTBYxkdSxSQ2wStmkjokroTeTxowc2oSkk47cXvJuVXQQBybzhzZ0acDVAE/wMDp7fhJNWdKA9EsB7l3VM6LV+LGwN4+WdiTIdUDJvOHYF8LJ2041JqScmJ4QgKIsY03UcNCdeFkQk8wYGUgUvULMpBYFQNmutKx7Axv4M4kEZybyBSEDEfB8hsosislIsUNQAlep/XGRRwMkHduL/HFdoPw8cl++fuRKWXa3XAVAmwvbWLbHSn2lRKJL/TUit4AlgZdeFrWEokoA3+9MYSmtoiyiwKTCQziMeVLCsM4rhjIbtI9kJnT8aJaebWNLun2VzccubbufdbCiBzXgGiMNpBNOy8UZv0lfkONuhlKI/WYAs1m4zFYgwrpCxN5nH1qGM74BIPzTTwmhGR0SR0BSSMZjSyrJAtk3RlygwD589kBZvhN1OQDK3mV2k3udcyJ/ZOlq3u+V1pxSxsuIuuzOmspEUlu07TbwYAIXQHFLwEWfi9W+f3u61ztdCFgUsbAl5pZJOR+/gloZGMkU3aHcIqiuyTuQMT4cliwIKOusEE0lRUD0RDRDA7qi/dOIBOOeIheMv3ABuYFEqync7wJpCyh79jpSWhVRZ8HRapcSDMg7ojHoaNdZlR8qsBloiAQQkgpxuIqMZWNAS8g06XNhAUnccBvt83A6wBa3lgdMHVnYjpIiY3xz0RNi1jqUy+KklgAacVnhRqNsJRkARqnMcAOuMWzk3DlUW0J8qoD9VQGdUxaq5cTSHFcRDMmyb1s00TQZ3feMFyJIglHkvTdZFfjrhARDnbYFm2hjN6ntULDwZGglGMpqJsaxe5kFSiSwQZLX6AVBWMzGQ0qo8aGqRypvI6hZCAdG7aOway3nlnkTewFBa22N39uNBKcWuhFMCc4zt9msPY15zELple6ZrlWimhc3OqIGD5pZngESBeCZ5pTogl9KSGwCcechcNAVl9CYLePCNgXH3WXKE1ADQ7mRs3NLQaFaHLDI3aEpZNmNXSYu/uz+KJEC32aDUgmkj7wS+M5WFc1nkOFBvL3Epn+4OMF9ouS4uKEuwGrhIG5YNSRLKAqBoQEJLOID+ZAGtkcC4A1NFgSCisN9lyssAse/IworMUXs0gOvPORQ/+tjBE3bE1kx/ATTgZEZE/6APcEZgEAFqA80JrZEAVs6LoyWiYFFrCAfOiXsBYEyVEXQ0h9NJTrcQUsSa7e8uskggCkJJBogHQBxOQ5g2RcGwvU6Q2UAyb2B9b3LcIGg0q6NgWnU9XmRRQN6wao6TMC3XW8bGztHcuNkKgKWYCYoC3uaQgsGU5mk6+pN5UNAJDwidLhJ5A1mN2eq7QQshxCvnPLbJvxvszf40TJvNG+ry0b0UhdDl9gKpvOFd5OY62wspEj75LtY59YfnejyX5noMVAzydPUlI44GyLAoXCeO0n1wAyBZZN5Elm17wz5bQsqMfQ4u81uCIGCBgHtxKk6BZ8eaLhjTeufOBMtCmWDcNawcD8OyEZCEsuCJEILOmIrmsIKFLaFx3ZUlgSDilcBYWdK1ZfAbSdIaCdQcVFsPzbAQUSTf/SGEIFTHC0gzLQRkAQG5se9HTJWxZl4TlnVFy75TqiyiKSQjO82+QQWDtb+P1z1KCEFAJl4pL93gDLc9CQ+AOG8LTMeBN5E3Zs0cmVTewFhWr3tCsRzvn/GmmEsi8czx/NAtVjJpjwYwktWrDOIqMZwSUOl23ZPh7kQeyZyBwZSGuDoz2R+gGBB0xtSyE/V7nUnU63YlfbNdxfKXf1fb3BoZINf/p73C4PHkg7owtymIRN7Abc+Xd536MVTRtt7qZEdGMhoEQqCZrK27ch/cYMj1JgIhky5/7QkCkuhp1NwyUNEEkR1jWjOYHcA0XUTd0pBcon1pdGCtYVEEfYKC1oiC/TsjZeL4WogC8UpL6YKJXWM52NRpLw9Nn1am0gG6kkhAqhMA2WWGjY3APHeqfxttkUDNTNNk0Uy7psi8ElVic8wkp5w601kgHgBx3hawHy1FQBKwO5Gb0ODNPcVoVkdGM+tmDTIFs257qossCjDsOgGQ4ywckNjE6h0jubqZsHTBRFazqrbrZoF6RrMomNa0t91PhKL+p7xM0RlTcdAcNjTzr+ur5wK+7gmg/YWo8xw9UanDNFBS/qooi8iigM8fywaF/nldX9loCD8GK9rW3eyImy3piAYQD8qglHrHWL0/1PEAmrgJ4p5kYUkZDCh6ABXHYLAxFYmc3lAWcjxMm0ISyzNAskggEjJuZtW0ba+7rxRZFNAdD0JoQLMkVkyEdzvAFraGpjT41Y96WiRFEnyNMwEWYMRVeVr2J6pKUETmVTUdWDaF6DiGN0JAEmHTol6ocqjw3oYHQJy3BabN7pjjQRljWWPGfziaaSFVMGDZLBNUi4xuMqO3BlLxlkVh1AhqdJONTRAdr41kvn4WKJnTYaN6+rmbaRnMzGz2B6gdkADAGc5Q0Ptf6ys7WRuWjY3OTKZa84hqlcB21Qi4ADbz6YjFLbBsihuf2Fq3LX6wogTWWlICA1hZTRIFjGR1T99Tun32mgBaw4Gqcloped3CYKqw11x+AWCR1wrvBECuB1DJGIzOmIqueBCDaW3K7vqmZUNxNDAuzA16fH0dBRBocHRILUSBIBos+gC5+h+/ERaTpZ4A2qVe+dOmNsKTKLv5EQlIiKjyuHrDRikYFlSlKCQfD0USYMNGizMOw8+aY2/CAyDO2wLTYsJSN5Bw75xnikzBRMGw0BSSMZo1ap6sR7OaVw6ph2uQVysDpJk23FtEUSAIyhJ2juR87+Qsm2IgrSEk+580W8IKRJCGT1p7Clcg7GZsSnnPEuYWnCqYZQNS3xrMQDdZS/B8n0CGrY89PpYzynxPSjvA/LjomCWQRYJXdibwzLZR32UsuzgI1dMARYoi6NIOG3d7rjh0IFXUrwUkEYok1B2DMZbXEQ/JyGgGdidyyBTMPT7Ox80AuYGApwGKKN7dflAWsV9HBBFVqrqAFQwLvYk8ehM59CZzGM5oyBRqaz1cw8FSJE8UXDvD5L4PU9VNSYKAuPP56KaNzYMsuPbT/9Tbl3qfi2bYUGsIoF0Couib9WKt/qTKAHGyEELQEQ0gbzRWwtRMq267esGwEGtA/+PC9F1FM8SxGZYz8ACI87ZAN22ITgo4prKW7vResn/3I6OZsGyKsCIhb5i+OiDNtJDMGePqf0oxTP8TacGwylL68aCMRF7HQLI6EHTLbrXS0rIooDVS27J+b1GrBAawIO+Mg1kW6O5XdnsXGFf/c9CceM39DykSWpwT7G4fDc68GoFTV1zFRw6ZBwD49RNbfYPLkawGm7KMXbOT+WkJKSBgwVGpQNgzXeyOIRwQQQH0VpTlamWAWOcPweK2MA5d2IKlHRFY1EbfHs4Iea3woznYlHolsLZwAKbN3H5liY1D2L8jAsNijta6aaM/lUdaM7CgNYhDFjRjRVccLWEFNpiZop942rRo1e9DEQVIpP6EdFY6I1MOgESBlW/cc8sGJ7tY2QFWj8G05mUF/dBMVoqulwWWJeLpAEvRTdt3BMZUiKkyJEFoqIQ5kmFZzHr6pEb1PwALOEGp95qZzuTzAIjztiBvWF45J6iIKJi256A7E4xkNajOXbxh2sj5pJSzmoW8bte98ytFIqSmoDqrWWyWjgMb4iijZyyP4Ux5KSJVMHwnkc8mdNP2sni1ApKTD+xEUBaxcyyPl3oSAGobIFbirnO302ZvWDb63e3VaY0+67B5aIsoGExr+NPLu6ued79z7dGA110niQLijmC2VKBfmnGa11StS3JdeoHqDFBGMxENSIgHZcSDMpZ2RHH4ohbMbw7uUQfd7ngQssiE3NuGs14JrzWisKGlImEXMbD3YFFbCKNZHWM5Nm/rkPnNWN4VQ0dMxYLWEA6e34R3L2YBXEar3m8btKq7iRCCgELqinUNqzg/bCqIAjueSEkWCGi8BEYphWlTx/rAf381s74AGih6AVUGfZppIyCKDbXAN0pUlRBSRGTHaYcvGBZkiRkXZgrV5yWbsm7H0ASCM7cV3tMAZXkJjMMZF3cIoks0IO3xu+Fa5HUL6YLp3ZUJREAyX30nkymYsClt2EBOFgVfQbVtU+R1qyqgiakSNMPCKz0JrO9NIZk3QCnFQKrQcNDVKD2jOVxy64u4/YWd01KG6U3kQcEmYTfVcKYNKRJOPrATAHDXK7th2RRv9BUzQPWYWzETrK9k6Gq9O1ZVFvG5oxcDAP7vxV1egOJSqf9x8XRAmdIAyBmp0Bws2Z+iLimZN6A5A1ErB3VmNRPt0UCZNkaVRbSEA7DrlIamiigQzzn5pR1jAICwwjp3TKdlXXba1AlhDsTLuqJYs6DJM9yrzMzJooDmMGvz9xPu+wUxQVmqWwIznAB/Ip1RfogCqRrw2RySPa+a8WAdWoIzLLfG/hKMa2IoOcdiVByzZlrenLTpQhIFtEXGL4Ml8wY6oiq64ypyPssWDAsBSWxYAO1uWxTg/eZ5BojDGQfbpjBM6rVOAuwuJlMwGjYFnE4ymgnNKGZ2goqIsZxR1Zk2ktUm2LpKmNi5Yj26ZUO3qwMgVs9X0RSSsTuRx8s9Y9jUn0a6YEy7vufmv2/DzrE8/ueZHbj9hZ1TXp+bCZnXVL/b5oyD50AgwCs7E3j8zUHkdAtBWcTitvoaDVfn4wZApeWv8Up/xyxtw8o5MeimjT9WZIGGKjrAcrqJkYxW0glW/D7uLM0A+QxpdYOplooBmaYzcLXZJ1ALB8SagcR0scgZAfFSDwuAWpx2csvR65S+f7IoYEk7azmv977GVAkxVS4rW9uUQgDxzVSqklBXBG1aNkKKOC1l3MqL+ET0P3nDQjggIaJKvtYAhmVDFklDNySVXkDpggFK0VA7/0RpCtd3hXb3Y05TEE2OR1VlSTivW4iqUl1/s0okgUAoGZrKRdAczjiYNoVJyzNAhBAoojgjnkCZggGK4gyvoCwir5frgPI66xKbSO3etcOvvLhppu3c8fqf7AOSiDnxIFRJxPaRLDTDbliU2Ahv9KXw4o4xb2L1rc/24M6Xdk1pnW5A4qf/KaUzpuLIJa0AgBv/thUAG2Q6XlatshOsdAbYeBBC8InDmTni428OlX0ermbHzdhkNBM53Sp6ATnfx5xuet/Nuc3BqoAMqN0Cn9FMRFTJd/ZTWJEQVPwvttOFqwNy9TBudsu07IbN+CpxDQoLVvEiykpqbBZWJYos1mwLB1gGaLqCfFkiCJf8TifSAZY3LLSEFLSGFWhWdTa6YFjjCqBdwgHJK/sl8wZyuoVlXVFfs8+p4rpC18qgj2V1tEcDaArKiKmsFFtZBtMsa0L6H4Cd4yRBwOLWMC49cX/80/v2m/QxTAc8AOLMekzbdhxjyy96klAtGtzTUEoxnNURLAkwFEmAbtnIldTUM5qJgpOtaBTmBUSrOsF004Zl2b6ToEsJByTMbQp5rsrTAaUUv3t6OwBmGHjue9jcqVue2o57Xu2d9HonEpCc6czrcjUL4+l/gKLOpy9ZgGXTkgxQYxe31fOa0BZh07Of3VYcyeHOF3Mnt5s2K8W4XS2jTgnMPb7mkIxIQCrRJOW9EuJAhaGiS84w0VlR/nIRBILWsILCNPm4+OGOxHAzMG4AZNHqjq2JEA/JCIjFmXemzb7TfoG9PE6Aa8N/2OhkUEShIgPUuP4HlCIWlBFVZQg+XVyaMxR3vN8uUOxoS+R0aJaFFXNimN8y/X5EACunxoOyNwC2FNNilhvdTSoEgRSD15IbAbc7bSLlL8DVXBGEAiIOWdCMReNkcvc0PADizHoMiwkNxYoTgeCMHdib5HQLWc2syuwIEMr8gFJ5AyBkQicvUSCwfDJAumUDEzgHTucJ85WdCbzem4IsEnzyXfPxicPn42xndMSvntiK+1/vm9R667XAV7K8O4ZlnVHv71Xj6H8AoC0agCKx+UoDqcKEAi6AfRYnLmcjOR7aUJwR5mqC2kvmdkki0OSIoN0SWGXLfVdMhUBYxsDNDLmGiqVT4E3LhkiEuvPZ4kEZlo091hJfGQC4JRgK6putaZRoQEJTSPF0boZFIdUQMsuiAIL6xzhdo0MUqdzHpjQDlNNrt/AXDObQHHFKYEwHVB6Y6paFpmBjWRJZFEApCwwP7I55juZ7ijlNQcgSwWC6UPY+J/MGWsJKmfdTU0iBKheD14LBxnOMZ/DqhyqLDQ903tPwAIgz67FsCsuqzgCJhMC26V79MWU1EwXDrtL2qDIzvrNt5gkyktUmd7dMUHXC1QwLZCIR0DRBKcXvntkBgE3Cdi+En373Anz0EJaVuf6xt/D9v7yBrUOZCa3XbU+v15FVyofXsJZ4RRKwX0dk3OUFQrx17xzLjesB5MeJK1gA9HJPAkNpDTal1SJoypyMXU3DiJcBKm+5d92J2XNsX4olsOKFxit/1RHhhgMSVFlg3lDjQClrzd81lmu4bNYSVsoubK0lLtDjZWbqQQhzkXYzCaZlV2mKXJgXkH8nmOtHFBCnJwMkCqTMUbo0ABrJ6DXL7K7+R5UFBCQRMVUuywKzDBEaLoOrsoi2iIID58S978qepD0awKq5TQgrrKHEcs6lBdPG3OZg2fk2rIhoDilexmgy+h8XVa49+HVvwwMgzqzHdDIglSdKQWBp+T0RANEa603mDYhC9b4EFaYDyhkWsrqFnOZf/spqJm7821vY4HQzVUJAoFXcRWY0c0Za2p/dNorNgxmosoCzDpvnPU4IwflHLcLHD50H4iz3pdtfwTX3bcD24ey46x11HJIFwrx3GuGo/drwqXfNx5dO2L/h98INPl7blfS2193g9gDWEr7SGcnxyKZBJHIGTJtCICwr4g7ylIiAuOpmgMpLYKUBV6UuqVgCK+5TVjfRFQvU1TiFFCbarTfV2w18epMFEMKyOmM5o6ERCKy7q7jfrWEFlLLRHY2UcurRFFKgSiyTYNoUIcV/fcwNmvh6AbnOytOVAZKEYgmsIxrwfIlMy4YisxEVfp45BcNCW6TY9dYSVqBbpa7ldEIePpGAhFXz4nt1LEpLWMHKeXG0RRT0p/IYzepoCslVwmtCCDpiAeiWBUop8/+pk6WshyqLbBbeLGB6/LU5nD2IOwesEsHJANXqZJgstk2xdTiDVN7E0s6INyvItilGsjqCPg7LAUnEiKV76X3NstAqV3dv/Onl3fjzuj48tGEQ//bx1VUdJ5IglImpKaXI+bTA72lsSnHrsyz7c8bqOVUlGUIIzjtqEU5Y3oHbnu/BE5uH8fTWETy9dQTdcdXzyBGIM3H93fNx+MIWAMXyF/OcYcdlWjZMu7auQxQIzjli4YSOwS0hPONoeDpj6oTfx5NWdOL13hQe3jCAg53RG60RFqAUDAuSSCCK8DJAGc2EZlq+povzmoN4dhsLjiilRT2RUwIzLBuSUL/8BbD3vjWsYHM2g2af57OaiUTeQFSVcGB3FB3OcQuEYPtIFh3R8d+Hha1hrHdMJ1vCAc94sJYQv1HCiojmsILRrM66ymqYhMqiAFlw3aDLvxOux9VU98VFEFiAKBBgzfwm73HdsqGIbAp7qmCWCX5tSkEIEFGLmbqoKnkGg5KjdVIlaUKZ4OlsXmiUSEDCQXPi2DKYQV+igLlNQd/vR1NQQdDxDyIECAcmNyyWBdGzIwDiGSDOrMe0KECrT3aiQGA7RmTTBaUUPaNZbB3KYiitYd3OBAZSrEae1U3kdLPmCU0grHU1kdMh+Py0LJtiraMnyRsWrr73jSofDEkkyJf4ibhT4KfrZN8oT2wexvaRHMKKiI8eMq/mcvNbQvj6Kcvx808dgqOXtgFgwuPdiTx2J/LYOZbHpoE0vnfvG7jvNaYX8mZylZS/xnIGBp0y03ThZl/cTEuj+p9Sjl7ahqAsoi9ZwKNvspEcbvnLvRCrMmtN9+aspTX0JatNHl0zxF2JPMZyBnTLdrJJ7MKa0UzEghKi6vgXlqgqAz4jGCybIpE3cEBnBIctbMaC1jBUWYQoEOzXEcHcphAG04Vxs6aLSjNAjgmiKPi3rE8EtwymmRZAak9+FwXiabgqMSxWgp5qNspFEgQsaAnhtxe8G19839KS7bAMzpx4EHmjfAxJwbA8/Y9LOCCxTLCTwdVMG9Gg2LAP2EyiyiKWd0Vx0NxYzQxUUBHRFg5gJKNNWv8DMIH7bHlHeAaIM+vRTMvXCIx1XYw/NHEi7E7ksWUwi5gqIxyQMJbV8dquJBa3hRGQmetzrdS7KokYyeiwbf+6/4s7xjCa1RFVJUQDEnqTBfzgLxvwg4+s9O78ZFFw2t7ZAFXd+X+/qdd7AkopntwyjF8/yVrOP3LoPM8ltx4LW8P45qnLMZAqYCSrO/OR2Poe3jiIhzcO4obH30J/yTys0uDAsJlbblYzGwoAGqEy4JmI/sdFlUUcs38b1r4xgAedyfRuAGTYNoIKuwj20wJawwr6kgW80ZuCaVMEJKGslFBqzuiKqVsjxW6vnG5hYWuooQtmJCBBVZj5Xul3bSynozWiYEFLqCpAkEUB+3dGYNo2BtJ5dMWCXqauEjczKTmeLTnNgiQI05KJjAdlBB3jwHqBvSoLSOX9vXXCgenzxhEFwpzVnU4uF9200RRS0BYJIByQkNUt76Kf1y3EQ3JZxpIZPsroHSsgqsrQLQvxBgXQswFJFMbtIG2PBtAzlkNYYdqnyW6H+tzQzgQ8AOLMejTTKjNBdBEIS0VPVwA0kCrgzYE0gnLRGK05rCCnm9gylGZDROsMNlVlEVnNhAWKZp8T34NvsAvoics7cOpB3fjaHa9i00AaP3t4M7528jIQwkoMmmaVBUCmRX2Pf7rZMpjBr5/c6pU+5jYF8SFnKnujdMbUqjvIlXPj6I6r+J9ne/Cnl3d7F/jKACWiyhhMFaYcALnvXeXJvHR7usl0JLUCgFJOWtGJtW8MeNkI1wTRtCiCMsv8UAovAFq3KwGABTyl63czXsMZDTtGc8663GyS7QUbjRBURERVCalcsSPRcLKFC32CHxdVFrGsKwrTptidyEMSCCgFiMCC1biqIOzM+Vo5J4aFrWEIhHjB3nRkM8JON9hwVqvbVRaUJVh2tdGp6czgmy4kgXjT50uPz7BtRFUZQUVEZ1TFtpGsFwBppu1rUNgUVLBzlJU4BUKmrVV/thAPyYipMlp9HL8bhZWN64862VvwAIgz6ykYdlULPOBMUJ+mAGg4o2FjXwqSIFR14IQUNshwJKshrta+owtIAkaybLRB5Z3yaFbH89vZhPH3H9iFuc1BfOsDy/Gde9bjb5uHMbcpiHOOWOhoCCh000ZIYSWwUtPFibKhL4VXdibw8cPm1bx7z+kmfv3kNjz0xgAoWKfVxw+dh48cMndaTuCEEJz9rgXojKn46cObvRPfXCcjYzsC20hAxJDjpTLRCy2lFOmCyZyGCdAcUhBSJHREA17nVmkGaCBdgCgQzGmg22ZFVxRzm4Kee7VrgmjaLDMniwIIoWh1LojrdrF5ZW7JyyXmzPVK5g287LgsF/U/FLI0sQtmazjg6YgA9h3riAa8/ahFSJGwojuGzqgOQuD8IxjNahhMagg7gzuv/ehq7zWWk9GaLjpiAeQNq+4sL0USfEuihExfCzzgZICEag8fgqImx818uCN5mAam+vIZUSXIIkFWY8c2kTlZbwcCkoiFrSFEJ6n/AVjAKc2AhYkfXAPEmdVQSmH4tMC7EEKmrBvJaiY29adhU3iGdpXIooCuWLBuRwchBKos+GqEHt44AJuyi6nbZrt6XhMueS9zQv3D8zvx1lCGeQHRohniVEYeUErxHw+9id8/14M/1zEt/K/Ht2KtE/wct387fvnpw/Cpdy+Y9rvX9y7rwPc+vBKRABvG6OpMDEds2h5VEVWlsnEJ42HZbNJ4XzIPSoBlXVG0RgJei3ipzshtizct2/F+YSXO8SCEeC3xQNG52aYUiiw62TmCFscLyB1WOr+lOrhys1Cv7Ew46yoGUxMd7hkJSKwRgFKvu2t+S+MltAWtIcxvCWFecwhzm4JoCio178pZx9b0fR9awgq642rdoEqRBFTG/a4OZ7oDIFEoD7Zcoz93/+JBlvVIFQwUDMvpxKt+P8KKhJAiYSynIyA35gD9dmNec8gb/jsZJEGAJBBYPADicOpjWGzScq2TOsXURNCUUuwYzSJdMKZl5k5rOFA1w8mmFGvfYOLnkw/sKnvu5AO7cPA81l20ZbDopeMGPlnNnHT5q2c054lx736l19fQrT9ZwGNvDgIAvnvGQfj6KcuqBnNOJyvnxvGrzx6OX376MK/d2DCZwV5YEdEVC5R5qYzHYLqAoCLgoLlxHL6wGYvbI4gEJBgVOqNYibeO7ohoF7aEoFt23XZylxOWdcD9GNzRBG7JUpYEiCKpCp79NEduQOa6Wrst8KZFoSrChIZeMg8aEXmdmSt2xdUJjyYoJSALoDW6c2w6veNVWCYhXDez6X7vS8XH7kiY6QyAJIFArHBx1k0bskS80R+CQNAVV6GZNrI606n5vR+i49StW0zTNp1DTN8puBPh6w273VvwAIgzq7FsxwW6zomkcgjpRBjO6OgdK5S5nk43r+9Ooi/JJrQfs39b1fNzvW4lFqwIIJ7jatZpgaeUljlNN8IzW4sjHEZzOh7dNFi1zB9f3gWbAocuaMJhC/2aqqefSEAqCxJ1y0ZQYV09zaEAJIk0lPmiThfUorYw5jWHvIxVJCDBpOz17mTz0mBEN1nGaV5zEIvaQhjNaeOWUVsjAXz5pANw4TGL0RVXPTM+xWnXlgipal/3M3ms1D11uoJqy/a1V6iHKjPzvbGcDlkkmN88tbEJqixCloQazsfMmHBvIktCVWnK1XdNJFM2HoSwY7No+XYUUSgLclrCCqIqE+rXCzRjQRkhRZw2Mf87DUIIAnJ1yXEm4AEQZ1Zj2GwuzYNv9OMPz/VUtf0SEF+TskbQTRvbR7J7XKz4oJP9Of6Adt/tuBdBt11bEglyOhNCayY74d/3Wh8+fdOzeMDpRGqEZ7YxzdGSdtbR88eXdpel+UezupeZOuuw+ZM4sumhtMstqkqIV0wNr4Vm2gjIYpUglpUt2MX6mKVtOP6AdnzyXcXj0y0bUVVyDP/C6IqpGHDGUtTjfcs6vLlkps1Ey64fjSgKZWUBAvh21FRmhdxymkknV2JqCStsblM8OKWyBMC6GANijUnzlO51LypZECALBMm8gbGcjpGMhkTOmFYPIBdFFMpupJgGTyq78QpIIrpiKiKqVLcFPBKQ0BSSJzQIeV9DlUTMgviHB0Cc2Y1lURR0Czc9uQ2/f67H61ByEQipGh7aKH3JPEYy2pTKBuORLhh46q1hAMDJB3b6LuNeBN0MkCwKyOkW6wAzWavwK05n0f88s6PmBOdShjMatgxmQABcdspyhAMidifyeLYkK3TXK7th2hQrumMNDRjdU1g29cphpaWG8SgYzG27MnAISKJX0oiqMr528jIcsqCY3TItiogj4pRFNlqjUT1Q6TpcN2JJZBfqeMkdf2dM9S3TlGaABIIywfJkAoxIQEJXTPVa7KeCIgkIyGJVAOQ6Xst1OiD3BIoksFKT8xk3hxUsbAthXnNw2geEKpJQlgFyg+RKOmIBtEUCda0hQoqI+S2hhjv69kVmwvDRDx4AcWY1hm0jXTC9u4UH3ijPgIiTnAifLhjYPsL8fvakUdmjm4ZgWBSL28JYWmOGlRsAucMxZZFltbKa6WgeBPQl2HOJvIF7140/gPRZJ/uzvCuKOU1BnLaKtbPf8dIup2PK8AaZfuKweXtk4nTDEHhaC4CNSwiUDF6shWbaaPFpxw3IAuQ6ZTRCaNn2oqqM/TsiKJhWw98ltxTjBi3MFK94Uq9lutgRVT1tS1ukZOQFpZPStTSFZKzojk1buSWmSlU3FG62a2+XwESB4JAFTXjPkhYcvqgFB89vworumFcynk4CkoBSSYoN6pvBiaoyVs2N1w1WCSFlLuecahRJAJkFb88s2AUOpzamRb0BfADw1JYRZEr+FsjE2yltm2LHSA4F3Z70hSOvW7j8rtdwy1Pb6y631gnYTj6ws2aQ4XYCubOaJIHpMDKaCQp2Vu5LFUs0f3xp17iDLV39zxFLWgEAZ6zuhiIKeHMgg9d3J3Hvuj4UDBuL28J7Tfvjh2kxgXvpHWEkIKE5xDpu6mFT/88vIDHthl9m0LIpBCJUdR+1RwJojShINqizMu3yC2RIEQFC0OTc9ZeWutxONYBd1N3SmCumtmwKUZycySAhxLcde7KEA1LVnCbLG4Ox9y8XhJC9EpxLJV1gtKIDzG+fOFNDEsle8TYbDx4AcWY1lk3L9CC6ZZeJeQUy8Vbx4YyGvmR+Sl1fL+wYxbpdSdz1yu6awchgqoDtIzkIhOl/ahEJFOcFDaY0SCIL6nK6BQpgLKtDN9nYhLlNQaQ1E3e/UrutPaOZeG0386J5z2IWADWFFK+Vu7Qt/qwZzv64XT2VF5vOmArTtqs0X8XX2RBFASGfVmQWFFSXctzXse2Vv05wAhPdshqyVTCtcs1OQBJh2RQtzliL0gxQRmP+RJUO2B2xogBaFsi0CnsniyqLIKjuvBKF6dfdzCZEkQCEHbNhUSgiQeAd2MI+W3Bb4Weamf/FcTh10EzLywC5P5cH1vd7J2hRILDtxjvBKGUOuCIRptRK+7Lj42LZFG/UmOzuGuLt3xGtm2kihHhZoIFUwXEPdrrfiIBep5W9I6ri00csAMD0O7WEwi/uGINlU8xrDpZpQz56yDwIBHi9N4W0ZqI7ruKo/aq70vYmbOCkWBUANYVkhBTJaxevJK9bCEoiQjUuUtGABMOnzZZ1gFVvD2CC4khALssw1sKm5YNbJZEABPjgym4c0BnBuxe1lG0zEiiWllY7Q1VXdDPdFdMTzY4Ag83YKnfptWyKgLx3MjEzhSQQuBOqdMuGLPl/RzjTgywSSHtZU+bHzO8Bh1MH3bS9CeuHLWyGIgrYMZrDpoE0AGceGLWr0va1SOVNjGT1KQkUKaWekR0AvLoz6bvcq7vZMgeXTJiuhSeEdpx9KVh5SBYF9Cfd6ekqjl7ahkWtIeR0C396ebfvup51pp+72R+XLuf1Lh8/bN6MD2pkc53EqourKotoiyjeZ19JwbDQFJLrjHyQfAcuGpaNUED09WcJSCK64yrSNbZZBqFlGRtZFEAAnHJQF/79rDVlbf7uNl1vog+u6sYt57/L84QybNsTU880bLCrWCZCN+2Jt+i/3RAF4t1UGaYNVZqeuWccfyRRwGx4e2fBLnA4tSmYlndBmtccxNFL2UX9wfWsfZsFQI2PwxjOaLDs2gNNG6E3USgbQbDOCXRKoZRinRMYrXaMDutR2QkmECYClkXimRl2N7HZUuccsRAA8Od1vVWaFcOy8cJ2NmbhiCUtqOTjh7KgpzMWwPuWdVQ9v7cxLBvRGhqW9qgKCv/P1qA24jVcuwGWySCodgmv1d3j0hYNQJXqC7BtSiGg/AIpO3PFfL+HhLVZuxkgQkhZ95dpUYRnScu0LAoIKeWt8JZNob7DsyEiIRAI0wGN9x3hTB3Jcd+eaWZ+DzicGlBKoRlFEXQ8qOCUg9hd8982DyGnm+zOrcGJ8JppoT9VQESZWsfMKztZgOGOtNg2lK0yKdyVyGPUMahb0TV+i3lpCQwAZIEgr7MMkFsC646zIOk9i1uwtD2CgmHjzpd2la3ntV1J5A0LzSEZB3RGq7azpD2Cn33yEPzbxw6eFXe4lKKm1qIpKCPq4wlk2RQSEWqWv4DanWAUqOvPElNltEUD3jgLP1jJCmWDPGXnjrbye+iKvBVHI+S7PtueVSMTogG57H2jQN2hpe8ERJFAEFkpvdSWgbNnYL+XmS+pvrO/1Zy3NaZNYdGiCLopKOPA7hjmNwehmTYef3MIAmEmco2UwBI5AxnNrOvh0Qiu/uf4A9qxoCUECniiY5d1zjIrumMNZZt8vYAMk7XAOyWwOfHiCIZPv4dpge5+ZTf+89EtGM6wjNQzTvnr3Ytba046X9ASqhrX0QjTNXi2dH2EoKbWQhIFdMUCVTqggmFBlf0F0C4BSYQilrsau0NXx/Mg6YypdY/VtG1IFW7EkkggidX2/u7YjVhQgm9NDmycizKLAqBQQIKN8uOYDQLtPYk3DoPSKlsGzvQjCrwLjMOpi2VTWJaNlJsBCskghHjaiQfW90NwM0DjtMJTStGfLEASSM3AoNF9coOdNfObvPLWuooA6FVHAH3wvKaG1utOBR903KAjKjO4Ewib1wUAXSWTyw9b0IwTlnXApux9+IffvYBfPbHV8/95j0/5a6oMZ9jQ0cn4LvnhdYDVufg3hwNQJOIN+wRYAFRrFpOLKBCEFKksk1Ecb1D/tNccktEUkmuOHjEsCqmiK4qNwxBgVnwPNdOGKrOxCKJQy7WczAoBtIvqlA/dUSMEmBX6pD2J6JwX3Gnvs8Wo753MnnTfb5R39reaM+PkdQvJ3MRmWLkYFhM3uxci12Plfcs7IAkEbw1lsXUoCxvjZ4DSmonRrIbYFA3jNg+kkdMtRAIS9muPeAHOqyWi6NIgqRH9D1DMAKU1EzndhEDYBTxVMJHTLRAUfWMAlgW69P0H4F8/ugoHzYnBsCjuebUXo1kdQVlsOPBqFJtSGLaNzpiKwXShZnv6RGik2yamSmgJKUjli8JkzbKqBo/6Ea0w9XOHro4XAEmigDlNQeQMy/c4TctGSCkXbgsCQUCp1gC5HWABSSjTAbm4M8UC4sxfDFxUWYQsMod106YQxdkVoO0JREIgCgR5w2ooSOZMHXUWZNlmfg8472hGshp6xrKTeq1pUZgm9YS+7qyjeFDGkfsxMfTDGwZAQMdtgx/NaNAMe8p3HW756+B5cYgCwcq5cQgE2J1gYzUAYNtwFhnNRFAWsX9HtQ7Hj6AiIuaU5gZKTA/7Eqz81RoJ+JbSDpoTx7UfWYWrPnSQ5zR93P5t067vyRRML+hrCikYzjQ+NqIWjXTbEELQGVdh2CwYYaJmUrf85RJUxLLZ5qVDV8ejNRxASBF92/AN239ulyqJVa33lk2LAZAkVJl2snIagSzNngCD7asIw6RM7yTMjAni3kRyNCkFgzVI8ABoz6PKImbaWYF/ypw9im7ayOStSU1sN22KnGF5niSlreuHOrOddifyoCj3LanEsGz0p7Rpccx129/XzGfbjwQkLGlngYdbBlvnzO1aOTc2IaFfh5Ph6U8VO8xcAbSr//GDEIJDFzTjurMOxg2fPhT/ePx+DW+zUTK6ga6YinhIxtKOCAihyDTSLl4H1gI//mfSHFI8TyDNYEFT5QBUP1gnWNHUT7dsRAONZQCDCht86edGbVP/QFqV/YTO1DnRE0R8zBkNzwNo9pyKJZG9v5ppMcH5LDGt29Mojmg+4mPLwJl+JFGAOMPzMGbPr47zjqRgWNBta1IDS03b9i5AQVksq8u7baoZzQRo/S6wRM5AKm9MeV5STjc9/6E1C5q8xyvLYBPV/7h4M8FKM0AlHkDjQQjBvObQtF9MdZOJfl3n7LZIAEvaI0jmjSnpgUxq152q7aLKItojAWQ0AwXDQiggNZQ+DzilHDfrYtn2hCZ0d8TYrK7KoIWA+GbjKh8zLBuyJHj7GlXlqgyR6bhAz6YACCiWD81Z5FG0p1FEAYRg2uaqceojz4LS6oS/1UuWLMHIyEjV44lEAkuWLJmWneK8cyiYFgyTNjTduxI2B8zR/4TKT0ruhdPNQtQbXzCYKoA4Nf6p8PruJCybojuululxSoXQhmVjfa+r/2ma0Pq7KlrhgaIAurtp6tO+J0uqYKAlpLBOJod5zWwq92Bam4IeqPa8pUraowEABFndREuoegCqHwFJgFwyE0wgExO3xlQZreHymWRMGOw/uZ09VnwvNNOGIoletkiVRdb3X0LlTLHZQlARQSnbv32lHCQ7btxcAL13aA0HsKgtPKP7MOFv9vbt22FZ1XVxTdOwe7e/My1n38S2KTSdiWdLu3gaxbDsEg8g/wAoq5kQfO7SXbKaieGM5ulrGqE3kce/PbAJT2weKnv8Za/81VT2+IHdMUgCwVBaw+NvDkEzbcSDMha2TmxqdbEVvlgCc00QSwOuvQmlFLploSuulgUdokCwX3sEzSEZvYn8hMthnvi3wYtNPCgjqkqgQMM2BrIoQJUEGKYNw2Jam4kILwWBTfXWTNsLsE2bQpb871xlUQClRUdh3bQRlIvlrYAkQCBCWbbStGhDWbC9javPMKzZ5VG0JwlIovNv3wj4Zhq3U3MmaXjr99xzj/f/DzzwAOLxYneLZVl4+OGHsWjRomndOc7bG8O2YVIbhmlPeGApABSM4hiMygAoXJIBEoCapZi8YaFg2mgJNzb4dOdoDpff9RrGcgb+tnkIO0ZyOOeIBRAIKdH/NJW9RpVFLOuKYn1vCn94rgcAsGpufMLt9h3Rci8gAOh1PYCaZiYAymoWwork23UVVESsmhdHbyKP3Yk8did0T68zHu5Q0kYduSVRQHdchW7avgLkWkRVCWM5A4rptsBP7GLeHGaBV6ZgIhaUWSAl+Gt2JJFAkog3PV03bcTU4ucWkETIEmGDXAW2HzbolFzJ9xSqzETbqbwxK9qV9waiQBCQhX3meDkTCIDOPPNMAExncN5555U9J8syFi1ahH//93+f1p3jvL0xLArTuWDUGy1QC920vcGUTTUyQDZl4tbK7hqXSl+WemwfzuKKu19HMm+gOSRjLGfg9hd2Yncij3PfsxC7xvIQiH9p6+B5TVjfm8KgMyKjnv7HpqxrrVJX4bpBu2WlrF4cBNsVm5kSWEozsLg1XPOiEFIkLO2IojOmoi9RQK/jgC0QANSdsM0mnbeVjH+oN5S0Fm2RAAqG1ZAA2iUckGDZ1OkAEyccbLjzwd4cyCAWlGFaFLIg+BoDMi8gJsiXRMCGjVBJdsdrhTfLRdSz0WTQNZIUBMIGve4DSCJBUJ74d4Tz9qXhM4ntiPcWL16M559/Hm1tMztFmjP7MS0bls0yJJnCxAOggmkho7HXNVVkIAIS60wxbcqE1jUyQIZl1zLgLWPrUAZX3P060gUTS9rD+N6HVuK5baP4xWNb8OSWYU/gvH9H1LdksXpeHL9/rvzvWgynNei2jZaKbImbAcobFlIF05s31hySZ0QnYlg2RELQFh0/exZVZUS7ZHTGVc+3iRB2w0QpxfaRHAbTBe8YDctGc1jxHUpai3BAwrIGxoqUEpBEEFAYpu3oiCZOWzSAHSNZFJyOxFjQf5iq7LlBuwaCpGyGliAQhAMSxrLMQoDNFJt9AmiAZUPCAQmpgjErpnbvDVRJRDwoz4oRDZy9w4QLcNu2bdsT+8F5B2JYzKAwKIvI66wVvtELHgueqCdArSyBEcJO0Mm8gYJhw7Sci05F2algWOMKZrcMZvDtu19HRjOxf0cEV39oJSKqhJMO7ERnXMW1923wBrJWlr9cDuiMIiAJ0JwLba2uLcsZ77GgOYRdYzkQEC+4USQBLWEFo1kdA6mCVwrrjs9M9iddMBEPyVXZt3rEg3LVZwUAYUXC+r6UFwTpVmMdYFMlIAuQJAEFx7xyMkQDElojAc+lu1YJThJZUK6bNnSLZb0qM2eRgOR9rqbjhD1bMw4xVcJQxj/b9U6kOaxMakQM5+3LpM4IDz/8MB5++GEMDg56mSGXm2++eVp2jPP2x9XlyKKAvGlCt2yoQmOZDNOmMO2iCLqyCwxgF5Nk3kDOsLyJ8JXp+pxuQa4TdCXzBq68hwU/y7ui+O4ZB5V506yaG8ePzzoYV9/7BnoTec+AsRJZFHDQnDhe6hnD6rnxmkFXqmAgFpSxf2cUAVnAW4MZAAEvCOqMBrwAqNcbgbF39T82pdAMG3nDwn4d4QllaWrRHFZwYHcMb/QmMZTWYINC3QtZLbfspBF70uJWQgi64ir6kwVYNkVQrn3aVJ1gXzdZ+3hlAFRqzuiZIM7SACOoSAjJ4j5TAuPse0w4ALrqqqtw9dVX4/DDD0d3dzc3jOLUxNXfKJKAlEa9uUgNvdZmOplaGSCgqAPK6yZsm3raCxdKKfKGVdfD5JantiFVMLGoNYSrPnSQr4B3TlMQP//UIRjL6V4Jx4+PHToXqYKBD6+ZU3OZnG5iUSsbkLqkjRkobhnMoMUJgjpjKjb0pzGY1jwX6HomiI0wltURVMS6733BsJAqGDBtCoGwz6wjGijT7UyVlrCCFd0xvNGXgq413gI/FRRnaOlkBNClMBsAGcMZra5rc0gRMZSmgMlKfJXllFJzxqIJ4uw8hwYkAUFFnLUBGoczVSYcAP3yl7/ELbfcgnPPPXdP7A/nHUTeML0ZO5ZNnVb4xsoprn7InSPmW1ZxxiHkdYsJi308Vkyb1nSxfX13Eg9tGAQAfPG9S+t2L8miUDf4AZg4+j8+sabm8zmdjcdodYIKQSBlQVCbECibCu+2wE+lBGZYNgqmhZxhYm6Tf1u+TSlGshrmNoXQFJIRUkSEFGY2ON03OK2RAA7sjmH7SHavtFcTQhBVJRRMa0oTviVRwJy4ioxm1i0JBSTmBm1T6pl1lj1fYs5oWjaiamOeRjNBPChjUWt41pboOJypMuFvtq7rOOqoo6ZtB66//nosXrwYqqrisMMOwxNPPFFz2fPPPx+EkKp/Bx10UNlyiUQCX/ziF9Hd3Q1VVbFixQrcd99907bPHIbtiD1rUTDssuBjIq3wpk2hm5anvakUQQPFDFBOL5bASjEsG6Zl+4oaDcvG9Y+/BQA45aAuLO+emLh2MqQKBtpjgbISmxsEdcZUZAqm1wk2kCoUTRCnkAHK6WxyuiqLyOn+Xj3JnIGmkIJlXVHMbwmhNcKyUXvqwtwaCeDgeU3TMpqkESIBGUFZmnLGqTUSQEtIqRsQyBIBCLND9BOul5oz1popNlsQBMI1MZx3NBM+I1x00UX4/e9/Py0bv/322/HlL38Zl19+OV5++WUce+yx+MAHPoCenh7f5X/605+ir6/P+7dz5060tLTgrLPO8pbRdR3vf//7sX37dtxxxx3YtGkTfvWrX2Hu3LnTss+cIjtGsl6Wwo/S8pMAUrMVXjdt9CbyyJcMnjQt6nWACQS+AtZwSQBE7epxGO4wR78U/l2v7MbO0RziQRnnHblwnCOdOqZlAyDo9MkiCQJBJCDBsKk3D2zHSA6jOdYtNJUMUE430R5VMCcexFiueoCpZVPkTRPzW0J79U5/b45WiKgS2iJTz7SEAxKWd0frjkqQBAHU6e5SfUpupeaMlFLuOcPhzCAN3YJ95Stf8f7ftm3ceOONeOihh7B69WrIcvnJ4Lrrrmt449dddx0uvPBCXHTRRQCAn/zkJ3jggQdwww034Nprr61aPh6Plxkw3nXXXRgbG8MFF1zgPXbzzTdjdHQUTz31lLdvCxfu+QvcvkiqYEK1bADVF2jDsmHa1Mu+yCJBVvMPgBI5HW/0phCQmdldR0yFYRXngMVU/9ZUzw1aN2H7ZIB0y4ZlV2eABlIF3Pb8TgDA545evFdm/6QKJppCsm8mC2DiWZvaXglsxGmVjgYkz/k4r1sgBA1fNN2SYDyoQJUF9CbzyOlmWakvmTcQDyromGSL+NuBlrCClmnKZIxn8iiLxDFKJDU/p6gqIZEzgFlqgsjh7Cs0FAC9/PLLZX+vWbMGAPD666+XPT6ROyxd1/Hiiy/im9/8ZtnjJ598Mp566qmG1nHTTTfhpJNOKgtw7rnnHhx55JH44he/iLvvvhvt7e0455xzcNlll0EU/U9ImqZB04rjB1KpVMPHsa/iCoxrOTCbFoVl2Qg4FwxZFJDT/FvhUwUDNqWQiIAtgxnsGs0jFpSrOsAymolwSWmmdB4YBWBVaYBsZkZTsd+/fPwt6KaNVXPjeN+y9qm9EQ3A3isT+3WEa3qMKJIAUDbzSiDM4BEAukscoMdyOiiYKLqR31petxBURERVCaosYk48iK3DGe8ibtnsM9y/M8KFrtOELAoQBUCpM1IhpEiO2Hz2doBxOPsCDQVAjz766LRveHh4GJZlobOzs+zxzs5O9Pf3j/v6vr4+3H///VXluK1bt+KRRx7Bpz/9adx3333YvHkzvvjFL8I0TXznO9/xXde1116Lq666avIHsw9i2kzEadrEGWtQfiJ3NQ6SlwHyb4W3bYrRrI6gLCKismxHXrcwmtVhOJqheFCGZloYyWiwQooniA6XzAMD/EtgqJAoPbttFC/sGIMkEFz83v32igA1q1sIByS01hnHIYtMcCwQ5njsOkqXlr8oWPYgVTB9ReFV29VMdMZVLxPR3aSWZYESOR0tYXnSBoGcaiTHDDGq+pslAsybiGJ2d4BxOPsCM377UXkB8jOz8+OWW25BU1OTN6LDxbZtdHR04MYbb8Rhhx2GT37yk7j88stxww031FzXt771LSSTSe/fzp07J3Us+xJuicswLd9J76Zlw64ogRlW9VT4nGEhp1ll5QK3HVyz3ABIgWVTRIOSF+wApRkgCwTVAVDeMKsyLn95rQ8A8OE1czG/eWLDSv0wa2TASkkVDHRG1bpuzookQBKZs3VnyeBTVwBtOcFkc0hGRjOqOt5898220Ropln6iquxpgdyuvAUtYZ6FmEZkZ05YTK1dcgtIojN5XOAlMA5nBplwG8ZHPvIR3wCFEAJVVbF06VKcc845WLZsWd31tLW1QRTFqmzP4OBgVVaoEkopbr75Zpx77rlQlPITTXd3N2RZLit3rVixAv39/dB1vWp5AAgEAggE+F3wRHAFxjalrLur4u0znWDE/a5IIpuCXdkJltNMFEwbrZHq4MBtgW8KybBsCkUQYUls9IUqi1UlsApPTuT18i60dMHAul0JAMDJB9b/jjVCumBgJKNjfkuoZmnLHScxXpZFdgzxDMtGZyyA13azx90MkDs8dG5TCMm8gVTeqKknApivT0ASq/RN3U0q+pJ59CXzaI8G0BbhXT7TiSAQxFQJqlKvVZ75Ekk1ZopxOJy9w4R/ffF4HI888gheeukl7+L28ssv45FHHoFpmrj99ttx8MEH4+9//3vd9SiKgsMOOwxr164te3zt2rXjttk//vjj2LJlCy688MKq544++mhs2bKlzKH6zTffRHd3t2/ww5kcrsCYAo6/T8Xzpg2/JEXlsqmCgVrXgESeCYHjQZn5+UgEreGAJ44uLYEJINCt4rr9TBCf3z4KmwKLWkOY0zT18RIZjY2KKO1eqySns+GdsXFKVu4F0bD8M0BumTGiSpjfEvKE3/W2G1ElhCuyTlFVRnc8CEkQML8ltFe7sfYV9u+Moq1OuTMgscxPUN5zVgMcDmd8Jnz26+rqwjnnnIOtW7fizjvvxB//+Ee89dZb+MxnPoP99tsPGzZswHnnnYfLLrts3HV95Stfwa9//WvcfPPN2LBhAy699FL09PTgC1/4AgBWmvrsZz9b9bqbbroJRxxxBFauXFn13MUXX4yRkRF86Utfwptvvom//OUvuOaaa/DFL35xoofKqYNhMYExAVDwCQB004ZQcXKvbIW3bYqRjO7bLgzA6ZRhAZBtUwQkAV1x1iFm2dTrjspoJgQBZRPhDYvCsu0yjcVTb40AAI5c4j/OYiJYNgUFK8PlDH9/HYBZAbRExh+wSAhBSBFhWnaNAIhCVQSIAkFHVEVTSPEyZLW22xnzF0vPbQ5iSXtoWl2eOUVkZ4p6LQghiAREBOtkiTgczp5nwiWwm266CX//+98hlEwIFgQB//zP/4yjjjoK11xzDf7pn/4Jxx577LjrOvvsszEyMoKrr74afX19WLlyJe677z6vq6uvr6/KEyiZTOLOO+/ET3/6U991zp8/Hw8++CAuvfRSrF69GnPnzsWXvvSlhgIyTuO4Yy5kUfDMCkvx099UtsLnDAt5RyDsRzJfLIGZNoUqCWgOKYipMjKaiYjTzaSbTG9U2pHGTBApVJntQ8Gw8HJPAgBqzvOaCO5+t4QVb2q7Hza1EQ82lnkMB0QMposZoKAsemJnw7IRVljAokgse/P67gRidnVwZVqs9OfnRMy2I2G/jmhD+8TZM7RFVPCh4xzOzDLhAMg0TWzcuBEHHHBA2eMbN26E5ZQgVLWxNl0AuOSSS3DJJZf4PnfLLbdUPRaPx5HL5equ88gjj8QzzzzT0PY5kyNvmJCcNt6sZlWJ1ytdoIHqVvh6+h+gGADFgzIsShGQRSgS8wraNJBBd1xlc5Wc7ZXqi1yRtuQE6i/uGINu2eiKqVjUGp7y8Wd0E3ObVTbvSfTvhNNMC4pTtmqEgCTCtikO6IjgmKVtOKAz4r2nFrXLPGjaIwHEgwqSeaPK4ybrlN32xrR1zuTY2wNuORxONRM+Q5577rm48MIL8S//8i9417veBUIInnvuOVxzzTVeuerxxx+vGk/BeWeR15nBoOK0t5cOOrUdsXNON/G7Z3fg5BWdWN4dq2qFr6f/oZQi4WaAggpsp20YAFqjAQRGstBNG6GAiKxmIW9YnvAacDRKtNiF9vRWp/y1X+uUdRfUmTvWFg4grIgIy6x1Xw6WH0zOyRKFGjQuVCQBhLDjvOzU5WXP/f/27jxMqupOH/h77lpL7zRLs4O4Ia6g/MCNGELixKjRGMckKFGTMSaKIdFoNDEhCy4TozNRZzQYx0mMZqJmdUMjqDFuCBFxAwWapaEFmt5qudv5/XGrbld1V3VX0dV0Qb+f5+kZqbpdfepCqC/nfBeRej7z2vF1Eby5pRVJx80a9Bm3HUwcFmV+DxFRL4oOgH7+859j5MiRuOWWW7Bjxw4Afu+eb37zm8Ex07x58/CpT32qtCulspGZYJxr0ruV2n15ZcNuLHt7B1o6Ldz4mSOgqyK41lCVXvN/Mnd0aiI6WmLJYEep0tRQFzWxsz2JqKH5AZDlwvX8btCqIlJHdH5AZLseXtu4GwAwuwT5P3HbRVhXURXWoakKaqM6NrfEeyQ6J2wXY2rCveaDZPIrggQ8KbPypzwpIdCzYmh4pYmxteHUzDAbNREduqpAytyz04iIqEvRAZCqqrj++utx/fXXBx2Tq6qyB0mOHz++NKujsuQHOH4Qk2vSuz+F3QtmT23dEweQKoV3/cAmpvhBS77jofTxl6kpCOkqJEQQAAkh0FAdwvbWOCpMDc3tScRtvyrK8TyoiprKB/Kv/+eWPYhZLuoiBg4Z1f/cl46EgxHVZhDw1UQMbNrVmXWNlP6g2Mpw4f8T0zW/MZ7jShhaVwCULoHv3jNGVxUcMboKo2vC2NISQ3N7ErbjoSKk5c3/ISIiX7/+luwe+NDQ4Lh+F+iInvrjI4GknZF/43hwXaAlVaXU3J70E3NVBRACSceFJ6Wf/5OvAiyjBN71JFQhgnwewA86qkJGEITELD+3KN39IGa50FK7KC+nqr9mTq7rUZlWLCklbM/LqqCKmn5ukuV4QZCSsP0dsUqz8DljhqpAUxU/uMwo0LRdCU3L3TRPCIG6qIHaiI49MRtb98ShKoXPCyMiGqoKCoCOO+44PPvss6itrcWxxx7baw7FG2+8UbLFUXnyK6z8tv8AoCoCMaurEsz2PHjwR1wAfsl4c3sSo2vCQSl8wnbz5v8AXSXw6SaImiKCnwf4OTCjqs1g3lLcduFKGcwDi9tO0Hzx5Q2p46+D6vv93hO25x9/ZTQYjBoawoaGmOXA0Pyjp5jloDKsIaQXnofT1Qyx20gP10PYUHvt2CyEQG3UCOamERFR7woKgM4666ygU3L30RM09FiuH+Ckd1MMVQkGlwLpfjwSuzq6ysO3tcYxuiYclMLHbTdv/g+QXQHmeB6UbgEQAERMDVHTf43OpAvPA1xXwvMkkpYfNL3T1IbWuI0KU8O00f3fsey0HNRFjazSfUURqK8wsGFn1zFY0nUxqSJaVMJ1uhfQ7g4r63Hb9VBfYBNPNtYjIipMQQHQjTfemPO/aWiyMxKMAT93JeF4QSm44/pdoHd1dn2QN+1JABNSpfCWi6SdP/8HQFYFWLADpGTvgChCIBwcgfk5QG7qiMqRfo5SuvrrhIl1JamKshwP9TnGWlSFdHgyXSHmr21vytAjhood3XaAHE/m7ZVERER7Z68+Efbs2YNf/vKXuO6667B7t3+88MYbb2Dr1q0lXRyVJ8t2IZBxHKUqcJyuqq247cJ2vKzBp9ta/URoXVVguS4SjpdVut1dayw7B0hTRY+Gf5oiEAnGYbiQ8I/b7FSOkipEVvl7MaSUPYarJh0XhubPeuouamowNQVJx+uaVbYXicghXc054oJDM4mISqvov6HffPNNzJ07F9XV1di4cSO+8pWvoK6uDo899hg2bdqEBx54YCDWSWUkbrtZwYiuCtieH/BETSBheWhLZHeHbmpNBNe2xmWv+T9AxhFYKgcoV7CkKCJrIGp6Irzt+DlKTW0JfNSehKkpOHZ8TcHvT0qJprYEFOEPWI0YKqKmho6Eg8qQnnNnJ2L4w1ljSb8f0bAKo9cALx9DVZARW/oNJoG9ei0iIsqv6H9WLlq0CAsWLMC6desQCnV1Mz399NPx/PPPl3RxVJ7itgs94zhKCAFP+jskUkokHLdHALQtoxQ+6bgI673H3l1HYH4AlCuZWFX8mUqAPxBVAv4RWCpHqbnND7rG1UYKDiCklNjRnkB1WMfRY2swZUQFdE3B7k4LMcvJO19LCIFhUQNx24Xluhi2l1PWdVWBkuoFBPjHX1qOEngiIuqfoneAXnvtNfz3f/93j8fHjBmD7du3l2RRVL48T8KyZY+E5PRQVMfzj59aU2Xs4+oi2Lw7llUKP74u0meybmtQBWb4YzByBDCaIlCRKjPvSDqA9NeXHhH6USqZeHiOnJ18mtuTqDR1TB1dhaqQjhFVwMT6KNoTNtoTTq8DRCtCOiAATSh7PYbC6NYLKJ1X1b0JIhER9U/Rf6uGQqGgAWKm9957D8OHDy/Joqh8Wa4H2/NyDDr1h6KmZ3Cly9gPHlEBQ1OCUnigsEqlzCowAD0CLsBPNK7MmAgvhF8xZdkuIBEMKa0vcDfmo/YkQoaKwxoqs8rcVUWgJmJgXF0EYSP/TlJlSENYVxFOHZntDV1VoKVmiwF+wrmhKllT7YmIqP+KDoDOOussLF68GLbtf0AJIdDY2Ihrr70W5557bskXSOXFTnWB7t6TJj0U1XL8KezpAKi+wsTo1ODHdCJ0X1xPoi3RdQQGoEcFGJA+AksnQTtQhR84+GM6RBAAFbIDtKsjCU0VOHxU5V6PkQjpKipDGmqjeq89e3rj7wApwVwzy/EQNVWWtxMRlVjRf0v/+7//Oz766COMGDEC8Xgcp556KqZMmYLKykr85Cc/GYg1UhnxK6xkj0nvRqq6q9PyGxKmx2AMixpoqA4DSJXCF6A9YSNdgFUV1iGRewdIVQQqUwFSp+UCQGoIq5+jtLMjHQD1Pnm7I+FAQuKwhkoM6+WIqxDj6iIYnXq/eytsqMEOkON5iBosgSciKrWi/2atqqrCiy++iL/97W9444034HkejjvuOMydO3cg1kdlxnE9ABJtCQevbtiFkw8ejpCuBkNRO5N+8nO6B9CwCgOja4rbAUoff6WPt1QhspKuM9VEuv4IJxwPydSsMU0V+Kij7yMw2/XQmrBxeEMlRvQRKBWitxyhQkUNDTtcf+0SgMmxFkREJVd0ABSLxRCJRHDaaafhtNNOG4g1URmzXA9SCjz8WiP+/GYTOpMuzj52TNdQVNsPkNLdjIdFza4doNbCdoC6V4CpCqDmyYEJ61owhytpu/6cMk/CUJSgE/XwPEGJlBIfdSQwuiaMsbWRYm7DgDJ1BV5Go0lWgBERlV7RAVBNTQ1mzJiBOXPmYM6cOTjxxBMRjUYHYm1UhizHgxDA5hZ/N+fDnR1dT0rAkxLSE8Ew02FRA3ErtQO0J/cO0KrGFmiKwLQx1RBCBBVg6SaIiiJ6HLml6aqCqKH6R1+2PxDV8fxmhJ7sSmDOpSVmo8LUcdDwih5J3YMpHfDYrscSeCKiAVJ0ALRixQqsWLECy5cvxy9+8QskEgkcd9xxmDNnDk499VScfvrpA7FOKhMxy4GmKEGC8ZaWrqBGVfxJ752WkxoHkWpkKP0doKyp8Cnrmzvw/T+tBQCMrQ3jjCMb0JHK56mJpMdgKHkDIENVENZVtMBG3PKnzPtJ2F0BWK7gJm65sF0PhzVUl92YCUNVoAoBy/GgKyyBJyIaCEX/zTpr1ixce+21ePLJJ9HS0oLnn38ehx12GH72s5/hjDPOGIg1UhmJ2x5UBUF+zZaWOGSqaZ+hKohbLtpSR1h1UQOKEKiLGj1K4dNe3rAr+O8tLXH81/Mf4tcvbwKQOgKT6Unwuf+o6qoSjMNIT4QHJHb20gPI9SR2xyyMr4tgRBE9gvYVQ/MDvpjlQleVYOI9ERGVzl790/fdd9/F8uXLg50g27bxmc98Bqeeemqp10dlxPUkLMdDws6e+9USs1EXNaBrClrjFtrifiL0sKgfXChCYHR1CBt3xYKp8GkrN7YAAL568mQIAfzlzSZsTR2V1VUYcFwPFWb+JGZF8cdQAEAs6R+BSZlRAp8j/2d3ZxLDKw1MqO+7IeNg0FUFmqagNWZjWIVRlmskItrfFR0AjRo1CrZt47TTTsOcOXPw3e9+F0ceeeRArI3KTLoHULrEPW1LSwx1UQOmpiBhe2hL+s/XRbsCl4bqMDbuigVT4QGgpdPC+o/8HKKTDq5HbcTAp49swD+3tOKdpjZ8cuooJBwXZo4xGGmaonQFQJYD1wOEQEYFWHYA5KVGZRQzHmNf01UFuiKCHkBERFR6Re+tjxo1Ch0dHWhsbERjYyO2bNmCjo6Ovr+R9ntWqstzS2f3AKhr0vvEYVG0xlI7QBnl57lK4Vc2+rs/U0ZUoDaVqCyEwDHjanDBCeNRFU4PQs3/x9TfAfLj+Jjlz+HKzFHqfgSWsF2EDQ1VYb3Ha5WTqKlCVfzmikREVHpFB0CrV6/Gjh07cP3118NxHHzve9/D8OHDMXPmTFx77bUDsUYqE7bjwXVl0OMnbeue7EToXZ1+8JE+AgOQsxT+9U1+ADR9Qm3enykleu2qrAoR7AB1Wk5qhlZXE8TuO0AdSQe1Ub3sA4uIoUHXFFaAERENkL3KAaqpqcGZZ56Jk046CSeeeCL++Mc/4sEHH8Trr7+Om266qdRrpDJhu36yc3p3pdLU0J50sKUllnVdZhPEtGAcRipYcj2J1akdoBm9BEAQvQdAmqIEnZI7kg7G1kagKfnHYNiulxWYlStTVxHWVZhqeQdqRET7q6L/efnYY49h4cKFOProozFixAh87WtfQ2dnJ37+85/jzTffHIg1UpmwXQ8SMsivOXpcDYDsUngA2NWRIweoJrsU/t3tbei0XFSGNBw8ojL/D5Wy1x49igJEQ6kdoKQDXVWQdDy0JfxjuMwAKGG7CGlq2R9/AYCuCoQNlTtAREQDpOgdoH/7t3/DKaecgq985SuYM2cOpk2bNhDrojKUdFwooiu/5tjxNXhx/U40tyf94CJ1rLS7s6sHT1q6FN5yPDS3J/F6qvpr+vjavAGOJyUU0fskdFURqDS7doCArgAsrKuIZkxv70w6qAxrWY+Vq7CuosLQWAJPRDRAig6AmpubB2IdtB/wh4x2zdiaMrwiOAZrao1jUn0FYpaDuO03MhwWNWE5HmzXQ9TU0FAVwqbdfin865t2A+g9/8f1JFQVve8ACRE0MuxI7foEFWCVZlYJecJxcVBlxX5RVl4Z0nF4QxWUMupQTUR0IOE/L6lgccuF60nsSY2qGF5pYmytf7SVPgZL5/9EDBVhQ8WemBU8lu7/89bWVmzcFYMAcNz4PgIgRfSRAyRQGcqeCL8z6AHUtQNlpzpQV4XK//grjcEPEdHAYQBEBbFdD7Yn0Zrwg5+QrqDC1IIhoukAqGsIqh98pDs5SymDUvin1+4AABw6qrLXfBzHk1AVpdcdoO5HYP6A055NEDuTDipMDRWh8hp7QUREg4MBEBXEdj04roc9nandnwr/eKlrB8ivBAtK4FPBhwRgaAK2K4NS+PZUrk6v1V9Aag5Y/kGogN83qDqiB9cnbC/rCCwtZrkYXmGW1dBTIiIaPPznMOXUuKszyOUB/ODCcf0ZWgAwvNLfzRmTDoBS5e1BBVi6sSEkNFWF43lBKXza9Al1va7B9SQiptJnzk6FoUEVAq6U6LScHmMw/InyQE1k/zn+IiKigcUAiHpojdv44KMOOF72FqGuKtjVkd1fZ2yNfwS2tSUOT8qsHkBS+nO5dFXxd4AyZoDVRHRMHh7tdR2O5yGk5Z8DlmbqCiKmivaEg45EVwCU3gGKWQ4ihhbkChERERUcAJ1zzjkFXffoo4/u9WKoPDTticNyZNbQ0rSPugVAI6tMaIpA0vGwsyOZ0QXa8I+wNP8IK265WaXw08fXQuljZ8fzUNC8Ll3154G1Jxx0JJ2gC3R6B6jTcjGuNsyeOkREFCg4AKqurh7IdVCZaEvY2N6WyHtc1P14SVMVNFSHsLkljq0t8aAHUF2FCceT0ITf0G93pwVFCIyvjWD9Rx04fmLvx18A4MErKGjRVQXhVA+i7a0JJFOT6usrTEgp4UmZ1ZSRiIio4ADoV7/61UCug8rEjtYEErabd1xEcyoAGpGRYDymNozNLXFsaYkHOUDDooZfxaX61WIS/hiNr39sCt7b3obZBw0rYDUCWi9NENNURSCSqgTbsKsTAFAT1mFoCmKWg5Cu8PiLiIiylORMwPM8/PnPf8bZZ59dipejQdKRdLCtNY6acO7dEk/KruOljAAonQfUuDuGllhXAJSu4gobKgA/kJkyogKfPmp0wc0IC6naUhURdHfeuNMPgNL5P5bjIaSpqTUQERH5+hUArVu3Dtdddx3Gjh2Lz3/+86VaEw2SHa0JxJNu0Fm5u9a4DduVEMgec5EuhV+7rRWeBBQB1EQMOK4ffJiq6ldpebLgtUgpoUBAV/r+I6oqXRPh0ztAmRVgDH6IiKi7oqvA4vE4fve732Hp0qV4+eWX4boufv7zn+Piiy9GRUXFQKyR9oGY5e/+VOfZ/QG68n/qoga0jO7M6WaIm1PNEGsiBlTFD3hChgJD8+d52a4HVSksGHE8CUVFwUdg6Ryg9m5DUG1PBsERERFRWsE7QK+++iq++tWvYtSoUfjFL36Bc889F5s3b4aiKJg7dy6Dn/3cjtYEOpMOomb+YOGjHPk/QFcvoLTMLtBhXU0FQAps1yt4PW4qgVorZAdICFSY2Tk+9akxGJ70CqokIyKioaXgHaDZs2fjiiuuwKuvvopDDz10INdE+1jCdrF1TxyVpt5rbk5QAdYtAKowNdRE9GBG2LD0DC7hV4mpikDIUNARd1Go9BywQnaAFEUg0i1wSzdqFOh9lhgREQ1NBX8ynHbaaVi6dCkWL16MJ598ElIWns9B5W17axztSQeVfczJ6t4DKNPYjJ5B6QoyKQE9lcQc1rWidoDSc8B6G4ORpikCFd0CoPqMQah6AUEUERENLQUHQE8//TTWrl2LQw89FF/72tfQ0NCAhQsXAkDBFT1UfjqSDhp3x1Ed6n33BwCa2xMAunZXMo1J5QEB/hGYlH6ydDpXKGpocIpIgvY8CUMTBf3ZUhWBim5l7sMrTDiuB0UR0NkAkYiIuinqk2HcuHH4/ve/jw0bNuB///d/0dzcDE3TcNZZZ+G73/0u3njjjYFaJw0AKSW27I4hZjkF9cnp3gSxNW6jIzXYdGxGHlBd0AOo6wjL0BQUEyc7noRZYOCiiK6J8IC/I1SbWoOuFFZJRkREQ8tefzJ84hOfwG9/+1ts27YNV1xxBZ544gkcf/zxpVwbDbA9MRtbW+N5mx521z0HqNOy0Zbw834yA6BhFWaQxGykdoDSHZ0LPTp1pURILyx5WVUEKjICoLqoAUUIvxO1qvAIjIiIeuj3P41ra2txxRVXYNWqVXjttddKsSbaBzxPonF3DNJDQYFGwnbR1q3EHDLd3rCrFB7I7gKdzuExNAWaKnIeg7mexK6OJHa0JbBtTxzbWuNwvcKrt1RFoDLctYOVXp/jetBVkVWyT0REBBQRAK1btw4XXHAB2traejzX2tqKL3zhC6ipqSnl2mgAfZQKOAqdkZVOgI4YKipMLVWl5f/x8aTE8AoT9RUGKkwNI6r8HSBd6Qo+DDV/KXxb3IapK5gwLILDG6pw5JhqTJ9QhxFVhe1MAUB1uGsHKH1E53iF7yIREdHQUnAZ/K233opx48ahqqqqx3PV1dUYN24cbr31Vtx9990lXSCVnuV42LSrMwhKCtE9/8fxPD+/R/oNDk1NxR3nHwvXkzA1FR0JB5Wh7EosPwDquQMUt12Mq6vApOF730vK1FSEdRVx283YAWITRCIiyq3gHaDnn38e5513Xt7nP//5z+Nvf/tbSRZFA2tHmz+1vSZS+IT07vk/jpvKr9EU2I4f1FSFddSmdpT83ZeuP15C+OMquu8AeVJCCP97+8PQlKCJY30qSJOQbIJIREQ5FRwAbdq0CSNGjMj7fH19PTZv3lySRdHASdguGnfHUWHqBQ0aTesRAKWOuMKGAivHsVauJOao0bMXUNxy/WO1PnoQ9cVQFVSngqiG6q4yfYMl8ERElEPBnw7V1dX44IMP8j6/fv36nMdjfbnrrrswadIkhEIhTJ8+HS+88ELeaxcsWAAhRI+vI444Iuf1Dz30EIQQnFKfoS1hoyNh99n0sLueO0AewoaKSlPPGQABskfycchQ0f0ArNNyUB3R+71ToykKPn/8OFx84kQcPa4mqDYrpJEiERENPQUHQKeccgr+8z//M+/z//Ef/4GTTz65qB/+8MMP46qrrsL111+PVatW4eSTT8bpp5+OxsbGnNffcccdaGpqCr42b96Murq6nEdzmzZtwre//e2i13Sg29NpQ1EElCKbVwZdoCuyh4yGDRVeztJ20aP8PFc5uuNK1BVYht8bRQEmDovis8eOhZKaPK+qbIJIRES5FfzpcN111+GJJ57A5z73Obz66qtobW1Fa2srXnnlFZx77rl46qmncN111xX1w2+77TZccskluPTSS3H44Yfj9ttvx7hx4/ImUldXV2PUqFHB1+uvv46WlhZ8+ctfzrrOdV188YtfxA9/+ENMnjy5qDUdyGzXw86OJCJ68cdNXV2g02Mu/PwaU1P9mRcZ0l2guzcgNFUVaio4Sa9H07J7+Oyt7kNTnW59iIiIiDIV/Olw7LHH4ve//z2ef/55zJo1C3V1dairq8Ps2bPxwgsv4He/+x2OO+64gn+wZVlYuXIl5s2bl/X4vHnz8NJLLxX0GkuXLsXcuXMxYcKErMcXL16M4cOH45JLLinodZLJJNra2rK+DkTtCQedlotokQGH36fHAgCMqMzOrzF1BZqWXd7evQt05vW6KoJrY5aLqKGVJABSlK6eREBGkjYDICIiyqGoT54zzjgDmzZtwpNPPon169dDSolDDjkE8+bNQyQS6fsFMuzcuROu62LkyJFZj48cORLbt2/v8/ubmprwxBNP4MEHH8x6/O9//zuWLl2K1atXF7yWJUuW4Ic//GHB1++vWmMWvNSU9WLsiVlwPAlF+F2W0/k1uqrA1BQYqf4+6WDD9SQ0pecUdr8ZogLHlYAOxG0HE6uiRa8nF1WIrI0o2/NzlErx2kREdOAp+p/e4XAYn/3sZ0u2gO7DLqWUBQ3AvP/++1FTU5OV4Nze3o4vfelLuPfee1FfX1/wGq677josWrQo+HVbWxvGjRtX8PfvD1xPYkd7cq/64jy6aisAYGRVCKri7+Boqp/j09XgsCv6cFwJLcckd1URMDUFsaQLKSU8T6I60r/y9zRNUaAqqdyf1P8P69z9ISKi3Pp/9rCX6uvroapqj92e5ubmHrtC3Ukpcd9992H+/PkwjK5eNh988AE2btyIz3zmM8Fjnucft2iahvfeew8HHXRQj9czTROm2f9E3HLWkXDQkXRQV0TvHwD40z+34k//3AYAmP///KPGdICjqwqEEKgwVTS3WcH3OJ4X7PZ0FzU1tMZsxG0XoVQVWSkoiv/lSQkVAo4rEd6LXCciIhoaBu2fyIZhYPr06Vi2bFnW48uWLcPs2bN7/d4VK1Zg/fr1PXJ8DjvsMKxZswarV68Ovs4880x87GMfw+rVqw+4XZ1itCXs1Gyswn/LX/5wF375wgYAwEWzJuLkg4cD8AMcXVGCBOPKkA7b68oBcnsZQRE1NDieRNxyURXSES5Rp2Y1VdmWrkjzpGQPICIiymtQ/4m8aNEizJ8/HzNmzMCsWbNwzz33oLGxEZdddhkA/2hq69ateOCBB7K+b+nSpZg5cyamTZuW9XgoFOrxWHo+WffHhxIpJZrbEwgV0Wvn/R3tuPXp9yABfOqIUTj3uDHBc44rURFWoKSOuExNzU5A7uX4ydD8bOWE4+GgitLtuqUDoHSFmRAMgIiIKL9BDYDOP/987Nq1C4sXL0ZTUxOmTZuGxx9/PKjqampq6tETqLW1FY888gjuuOOOwVjyfqnTctEedwqu/trelsCP/vI2LMfDceNrcdmpB2XlZdmeh7DeFbyYmgIBf/dFESLVBTr3zzI0BZB+g8JSVH+lqUJAVQQ86Qd8UvZMwiYiIkrbq08gz/Owfv16NDc3Bzk2aaecckpRr3X55Zfj8ssvz/nc/fff3+Ox6upqxGKxgl8/12sMNW1xG0nXw7ACJqMnbBc//svb2BO3Mak+iu986tAelVSulz1k1NQV6JqA5Xipoy/ZowQ+TVcFdE3A1Po//iJTcATmSb8HkNazESMREVFa0Z9AL7/8Mr7whS9g06ZNQTl0mhACruuWbHFUGh+1JwseCfHfz3+ATbtjqInouPGMqYgYuf+IZO6umJoalML7AVDPHkBpfi8gBXVRo6Q7NEL4PzPueqkkbe4AERFRfkUHQJdddhlmzJiBv/71r2hoaCioZJ0GT9xy0Rq3CzpueubtHXjmnWYoArh63qEYlidHR0qZFVyoikDE0LAnZvnHYBA9ukCnGaqCClNDbZHVaIUwVAWdngPXk9AVNkEkIqL8ig6A1q1bh9///veYMmXKQKyHSqwt4Zec14R7LzffuLMTdz/vD7v9wgnjcdTYmpzX+X12lB4jJipCKj7qkKkZXMi7AySEwEEjKhAp4DiuWLrq5wDZnocImyASEVEviv4n8syZM7F+/fqBWAsNgJZOC4oQve7UxSwHNz35birpuQbnzcjfLsDxUk0QtezXC+n+UFTH9RsR9rb7UhXSc/YI6i9dVeCm1hBiE0QiIupF0TtAV1xxBb71rW9h+/btOPLII6Hr2TsLRx11VMkWR/1jOR52dlqoyJPHA/jHWXct/wBb98QxLGpg0ScO7XVSfL4ZW+mhqP4YjME5fjI0BZ4nIYXMm7tEREQE7EUAdO655wIALr744uAxIUQwwoJJ0OWjLWEjlnSyBph29/y6nVjx/kdQBHDNpw5DdR9HZY4n/Uqu7gFQaihqwnERNQfn+Mnf6fKH07MHEBER9aboAGjDhg0DsQ4aAHtiFiSQNxixXQ8P/GMjAOBfjx+PqQ1Vfb6m43qoCvdMYE4PRW2POxhWUfoE50KoioBEqgkiE6CJiKgXRQdA6SaFVN5cT2Jnh4VoL0dBT7zVhOb2JOoiBj577Ji812Wyu/UASjNUPwBKOh7CA5DgXAg/0JPwy/AZABERUX57nSjx9ttvo7GxEZZlZT1+5pln9ntR1H/tCRsdCQd10dy7MZ1JBw+9thkA8IWZ4/PO7urOk56f79ONEAKVIQ1CIOfz+4KqCEgJ6GyCSEREfSg6APrwww/x2c9+FmvWrAlyfwAEVUbMASoPrXE7la+TeyfksVVb0Z5wMKYmjLmHjyz4dQXy59dUmDp0NX8TxIGmCgFVKH1WoRERERX9KbFw4UJMmjQJO3bsQCQSwdq1a/H8889jxowZWL58+QAskYolpcRH7cm8R1G7Oy38YfVWAMBFsyYUnLDsB7v5gwtTVxA21EELPhRFQChgE0QiIupT0TtA//jHP/C3v/0Nw4cPh6IoUBQFJ510EpYsWYIrr7wSq1atGoh1UhE6kg7aEzYqQ7kruh56rRFJx8OhIyvx/yYPK/h1Hc+f8ZXveMnUFIR1bdCOnzRFQBX+DhWbIBIRUW+K/mey67qoqKgAANTX12Pbtm0A/OTo9957r7SrO0B5nuz7on5oSziwHJkzF2drSxxPrd0OAPjyiROLGmXiz9jKv7tiairCujJoCciqIqAoImeSNhERUaaid4CmTZuGN998E5MnT8bMmTNxyy23wDAM3HPPPZg8efJArPGA0pawsXlXDIc3VEEZoF2Kne3JvEHK/768EZ4ETphYhyNGVxf1uo7nQc8xBiMtbKiYNLwC0UEKQBQhoAoxaFVoRES0/yg6ALrhhhvQ2dkJAPjxj3+MM844AyeffDKGDRuGhx9+uOQLPNC4rkSn5aDTcvIeUfVHzHKwJ2YhavYMAvbELPz9g10AgAtnFd/OwHElKsJKr4FbX40UB5KmCOiaAp1NEImIqA9FB0Cf/OQng/+ePHky3n77bezevRu1tbWcDF+guO0iZrkDEgC1xR0kHC/ntPXVm/cAACbXRzFhWLTo17Y9D2E994T4cqAoAqaaf4eKiIgoba8/KdavX4+nnnoK8XgcdXV1pVzTAc+yJdoT9oC89u5YEmqe4aerGvcAAI4dX7tXr+3kaYJYTsYPi6I6Mni7UEREtH8oOgDatWsXPv7xj+OQQw7Bv/zLv6CpqQkAcOmll+Jb3/pWyRd4IJKQ2NVhBT2USiXpuNjdYSFq9tzYk1Ji1eYWAMCx42v26vUl8vcVKhd1UWPQGjESEdH+o+hPs29+85vQdR2NjY2IRCLB4+effz6efPLJki7uQBa3/GOwUupM+q+Za5dm464YWmI2TE0paOZXLkKi7AMgIiKiQhSdA/T000/jqaeewtixY7MeP/jgg7Fp06aSLexAl3BcdCadnLs1eyvpuJDwq6G6W9Xo7/4cOaZ6r4IY15NQe6kAIyIi2p8U/WnW2dmZtfOTtnPnTphm+SbIlhshBNoTTklfM2G5yJeGviqVAL33+T+e3wRRY6I7ERHt/4oOgE455RQ88MADwa+FEPA8D7feeis+9rGPlXRxB7KQpmJ3rLR5QO1JJ+fuTsJ2sXZbK4C9z/+xHQlN5YgJIiI6MBR9/nLrrbdizpw5eP3112FZFq655hqsXbsWu3fvxt///veBWOMBKayriCUdxG0XEaP/x2BSSnQm3ZwBytptbbBdieGVJsbWhIt+bdv10JqwMHFYlAEQEREdEIr+NJs6dSrefPNNnHDCCfjEJz6Bzs5OnHPOOVi1ahUOOuiggVjjAcnUFSRsF53J0iRCJx0PtuvlnMP1Rir/59hxNUX3arJdD83tSYyvi2DKiIqSrJWIiGiw7dXWw6hRo/DDH/6w1GsZUkTq/3YmbQyv7Dt3Km65sBwvb4+bpOPBcj1U5kiqTuf/HFdk/o/jemhuT2BMTQRTRlQO2owvIiKiUturACiRSODNN99Ec3MzPM/Leu7MM88sycKGgpCuYlenhYn1fV+7qzOJ3Z0WjorU5Hw+6bhwPdkjSPmoPYnNu2NQBHD02Nzfm4vrSexoS6ChJoxDRlXA4HgJIiI6gBQdAD355JO48MILsXPnzh7PCSHguqXtbXMgC+kKOpIOEraLUB8DPD0PaE84sBwvZzCStD0gRz716lTzw4NHVKIiVNhvtycltrfFMbI6hENHVbKxIBERHXCK/mf9N77xDZx33nloamqC53lZXwx+ihPSVSRtD53JvsvhLddFwnERt3Pf47jl5hxS+kYw/qKm4HXtbE+iJmLg0FGVfQZmRERE+6OiA6Dm5mYsWrQII0eOHIj1DCmKEJCQBSVCW46HpO0hkScA6rDsHk0KXU/in0Xm/7TFbaiqwCEjK0tSnUZERFSOig6APve5z2H58uUDsJShyVRV7O60+rzOciQsx0Msx26R43qIW16PAOiDjzrQnnQQNVQcMrKyz5+RsF10Wg6mjKhAXbTnNHkiIqIDRdH/xP/FL36B8847Dy+88AKOPPJI6Hp2VdKVV15ZssUNBWFDRXvSRtJxe821sVy/y3Nbju7RSceD7Xg9cnzS4y+OGlsDNcfxWCbXk9jVmcTEYVGMri6+VxAREdH+pOgA6MEHH8RTTz2FcDiM5cuXZ/WVEUIwACqSqalojdvoTOYPgDxPwnX9YKkj6aT6/XTt9qRL4Ls3KXxtox8AFXL81dyewIjKECYNj+bMJSIiIjqQFB0A3XDDDVi8eDGuvfZaKApLo/tLVQQkgJjl5D12cqWEKz1EDC1IhM4OgHoOQd3daeG9He0AgOMn9h4AtcQsREwNU0ZWsOKLiIiGhKIjGMuycP755zP4KSUJ2I6X92nXk3A9CVNTYDsSCSs7ETqZIzH6tY27AQCHjKzAsIr8jRY9KRGzHEwcFkFVKHeTRSIiogNN0VHMRRddhIcffngg1jJkaYqA5eYPgBxPwpMSqiIg4O8WZepIOtC7BaQvf7gLAHDCpGG9/uz2hIPKkF5QN2oiIqIDRdFHYK7r4pZbbsFTTz2Fo446qkcS9G233VayxQ0VqiKQ7HMHyD/iMjQFrfGuAEhKiY6km9UcMWG7+OeWPQCA/zepLu/r+t9r4/CGKh59ERHRkFJ0ALRmzRoce+yxAIC33nor67liB22ST1EErAKOwBQBmJrfPdpxPWiq4idAOx7MjABoVWMLbFdiZJWJ8XWRvK/bmfQn0XP3h4iIhpqiA6DnnntuINYxpKlCwHYkpJQ5g0jH8yCEH2CamorWhIW47aJSVWC5HmzPQ4Xa9Vv5ygY//2fmpGG9BqWtCQtTRlSw4SEREQ05zGQuA6oi4EoPtptjmBf8HaD0M4amwHZlMBIjaXtwnK4SeNeTQQL0zF6Ov2KWg5CuYkRVqHRvhIiIaD/BAGgfeqepDTf84S08+Epj1uOqIoJjrlxcT3YbdCqRsPwjs6TjAhm7PO9ub0NbwkHUVDG1oSrvWlrjNkZWmaz8IiKiIYlnH/tQ3Hbx1zVNqDS1rOMufwdIwvE8AD2TkV1PZsY4MFQVe+IWxiPiD0HNeDJ9/HX8hDpoau74Num4UFWBBnZ8JiKiIYo7QPvQEaOroKsC7UkHO9qSweOqIuC6yLsDZDleVi6PqSnoSDhwPZkqge967tVUAHRCL8dfe2I2hleYqA5z94eIiIYmBkD7kKmpODQ1lDTdpRnwy9s96cHJFwC5HtSsAEhF0nHRkXQQt7pK4De3xLB1TxyaIjB9Qu7uz47rQQIYXRNm1R4REQ1ZDID2sWljqgEA72cEQIBf4ZVvB8h2ZNYwU131y+ZbY3bWDLD07s+RY6rzVnZ1Wi4qDA013P0hIqIhjAHQPnZkEAB1ZD0uAdh5ukEnHReZ80mFEIAA2hI2bEcGO0CvpLo/z5ycv/tz3HZRG9Xz5gcRERENBfwU3MfSO0Af7uzoEfDk2gHyUtVharcJ7bqipoagSihCoDVu493t/q7SCRPz5/940kNNJPfQVSIioqGCAdA+Nq42jKipwnYlNuzsDB5XRO5u0JlzwDKZmgLLlpCpmGndjnbI1Ovn6+ycsF2YqoqKEIv/iIhoaGMAtI8JITBpWBRAdh6Qpgi/p083XWMwugVAuoJOy4Gm+o837o4BAManXjuXuOUiGtIQNTj3i4iIhjYGQINgUr0fpGRWgqki90BUV0q4OXaADFVBzHKCIaZbWuIA/B2gfOK2i/qoweovIiIa8gY9ALrrrrswadIkhEIhTJ8+HS+88ELeaxcsWAAhRI+vI444Irjm3nvvxcknn4za2lrU1tZi7ty5ePXVV/fFWylYOgB6f3tGKbwiYKXmgWVyXQkvNQk+kxAC4+uiqEodZ21u8XeAxtXmHn7qSQkIoJLVX0RERIMbAD388MO46qqrcP3112PVqlU4+eSTcfrpp6OxsTHn9XfccQeampqCr82bN6Ourg7nnXdecM3y5ctxwQUX4LnnnsM//vEPjB8/HvPmzcPWrVv31dvqUzoA2taaQHvCBpDZDTo7AHI8L2cSdPp7hBCQUnYFQHmmvydsF2FdRYXJ/B8iIqJBDYBuu+02XHLJJbj00ktx+OGH4/bbb8e4ceNw991357y+uroao0aNCr5ef/11tLS04Mtf/nJwzW9+8xtcfvnlOOaYY3DYYYfh3nvvhed5ePbZZ/fV2+pT1NQwutofQpouh1eFgJsKdjK50u/2/N8rPsD65o4erwX4nZ07k36p/Jia3EdgMctFdVhHSGf+DxER0aAFQJZlYeXKlZg3b17W4/PmzcNLL71U0GssXboUc+fOxYQJE/JeE4vFYNs26uryl4Ynk0m0tbVlfQ20Q1IdodOJ0Joq4Hk9d4BcT2Llpt34y5om/O/Lm3K+Vnr3Z2RVKOgJ1F3ScVEXZfk7ERERMIgB0M6dO+G6LkaOHJn1+MiRI7F9+/Y+v7+pqQlPPPEELr300l6vu/baazFmzBjMnTs37zVLlixBdXV18DVu3LjC3kQ/pAOgdO8eRQi4np/zk8lxJVrjDgBgXXN7jxwhANicqgAbmycB2nE9aIqCSpa/ExERASiDJOjuFUmZU9J7c//996OmpgZnn3123mtuueUW/Pa3v8Wjjz6KUCiU97rrrrsOra2twdfmzZsLXv/eOnRkBYBU/55UlVfXRPgunpTotPwAqD3hYEd7ssdrbQ4qwHLn/8QsFxGD+T9ERERpg/aJWF9fD1VVe+z2NDc399gV6k5Kifvuuw/z58+HYeQ+1vn3f/93/PSnP8UzzzyDo446qtfXM00Tppm7eeBAmVgfDSbDN7Um/OGkkD1ygCzHQyzZ1R9o3Y52jKrKDub6qgCL2y7G1IY5/oKIiChl0D4RDcPA9OnTsWzZsqzHly1bhtmzZ/f6vStWrMD69etxySWX5Hz+1ltvxY9+9CM8+eSTmDFjRsnWXEq6quCg4f4uULofkARgdwuAko4b7AAByJkIvWV3agcoTwWY63moibD8nYiIKG1QtwQWLVqEX/7yl7jvvvvwzjvv4Jvf/CYaGxtx2WWXAfCPpi688MIe37d06VLMnDkT06ZN6/HcLbfcghtuuAH33XcfJk6ciO3bt2P79u3o6MhdQTWYgkTodD8gKXrkANmuRMzK2AHqFgB1Jh3sjlkAcucAJR0Xhqbw+IuIiCjDoH4qnn/++di1axcWL16MpqYmTJs2DY8//nhQ1dXU1NSjJ1BrayseeeQR3HHHHTlf86677oJlWfjc5z6X9fiNN96IH/zgBwPyPvbWoakAKL0DpAgBy80eh2E5HjqT2TtAnuwajZE+/qqLGojmCHJilouoqSFqMAAiIiJKG/RPxcsvvxyXX355zufuv//+Ho9VV1cjFovlfb2NGzeWaGUD75BRfgC0YWcnLMeDqmQPRE3PAevICIDitoute+JBvk+6AizfCIy47WJMTRhKjkaKREREQxWzYgfRyEoT1WEdjifx4c6OVADUdQTmehKW6wVHYOkmh+t2dB2D9VYB5noSAkA183+IiIiyMAAaREIIHBKUw/sBUCJjIrzryWBUBgAcO64GALC+uWuGWNADKEcCdNxyETE1VIUYABEREWViADTIGqr9XZ1dnVZqHEZXKbzjeWhL+MdfEUPFoakjs8xE6PQU+PE5jsA6LAf1USNvd2giIqKhip+Mg6w6NZ29NW5BUfxdH9v184A8D2hPdYGuDGk4eIQfAH34UScc10PScbGjLQGg5w6QlBKelBx/QURElMOgJ0EPdV0BkA1V8eeBZe4AtSf9I7BKU0dDTQgRQ0XMcoPqLwmgwtRQE84+5oqnpr9XhXn8RURE1B13gAZZZgCkKQoc2TUQ1fUkOlNdoCtDGhQhMGWEnzP0/o4ONKYbINaGe4wP6Ug6qIlw+jsREVEuDIAGWXrnZk/MhiKQtQPkShn0AEoPMk0fg61r7gh2gXIlQNuuh/qKfTveg4iIaH/BI7BBli5Rb0vYEEJAAsFAVMft6gFUmarkOji1A7S+uWsmWPceQAnbRUjj8RcREVE+DIAGWfoILGF7SNj+cZebcQSW7gFUEewA+QHQxl2x4LnuPYA6kw6qIhqiBo+/iIiIcuER2D4mhD/yIj3zNKyr0FU/f6c1bgOQcFLzwCzHDY7AqlIB0PBKE1UhDa4n0dTqV4B1H4KacF0Mrwj1yAsiIiIiHwOgfSykqzB1BclUw0MhBKrDfql6a9yGgAjK4DO7QFeYenD9wakZYgBgagqGV3bl+liOB11VePxFRETUCwZA+1hIVxE1VcQzJrxXh/3dnda4DVUIJJ2uAKj7DhDQdQwGAGNqw8FgVMA//qo0NVRy+jsREVFeDIAGQW3YQNLtGnoa7ADF7KyBqJYj0Wn5AVBFngCoe/5P3HEwotLk8FMiIqJeMAAaBNFUMCOln+uTLoVvTfgBUNJx4XoSXsYk+Mx5XulSeCC7AsxxPahCCQIqIiIiyo0B0CCIGBpMVYGV2gVKl8LvSe0AOa4/DiPpuEjY/jUVGUdatVED9RV+kJOZAB2zXEQMNegZRERERLnxk3IQRHQVIUNFwvZgamrWPDBVEbBcD0nbQ2tqDpgAEO2W03PJSZPxRmMLjp9YFzxmux5qowY0lXEtERFRbxgADQJFEagJ69jSEkd1WM8IgBwoqYnwSccN5oBFTQ1qt5yek6bU46Qp9VmP2a5E1GTvHyIior5wq2CQVIX1oOFhTbcdIFdKJB0vaxJ8ITx4CBuMaYmIiPrCAGiQRAzVD3Y8GfTs8QeiCngu/AAoWVwABAiYGn9LiYiI+sJPy0ESMTSEdAUJ283YAfKPvDz4R2CdqSOwypC/W2RnlM53Z7seNFXAYABERETUJ35aDhJDU1AZ0pGw3WAHyHYl4ql5YJ4HxJL+f1eaGlpiFppa43lfz3Y9GKrCHSAiIqIC8NNyENVEdCRdFyFdRUj3fyv2xGwgNSesI+MIzPEkNEUJ8oa6sxw/ADJYAUZERNQnfloOooihAdJviJiuBGtLHYNJSHSmxmVUhnR4UiKkK0GX6O7SFWAcgEpERNQ3BkCDKGqqMHUVScdDTap785541ziM4AgspEEAMHU1aJ7Yne15qAxxACoREVEhGAANorCuIqyrqTygjIGoioDtZswBMzVIAIbaNUU+F1NjDyAiIqJCMAAaREII1EYMJDJ2gNIBkOW6GZPg/Z0dXVXguD1zgDwpIQBWgBERERWIn5iDrDKswcvIAWqN21CFgO105QBVmCogAVNXAMhgiGoaK8CIiIiKw0/MQRbR/TEXFalmh+mBqJbroSPh7wBFTR2qIhA1NeiaAtvtHgBJaBoDICIiokLxE3OQRUwVIV1FRPfzd9oSfgDUkXSChOeoqUJRgApDg6mpPSrBbMdDSFM4BJWIiKhA/MQcZLqqoNLUEDb8AGhPzIKpKdBT5eyKAExNgSoETF1B1NB6JEJbrlfEuAwiIiJiAFQGKkOa3xMIfg6QECLYzUlXgKmKgCIEqsNaj1J4V3rB9xMREVHfGACVAV1VEA35vxVtCQeelBmDUP05YIoQ0BSBsKH1SILmEFQiIqLi8FOzDKiKQKXpV4G5nkRn0kF7Ij0IVYMnUztAikBIV6BmjMRwPekfj7EHEBERUcEYAJUBTfGPvKJmKg8obqM90TUHzPVk0OMnpKswtK6GiLbrQecUeCIioqLwU7MMaKoCSKA61DUPLAiATH8OmK76SdGmpiBs+OMzgPQQVJVHYEREREXgp2YZUBUBVRFBM8Q9MRsdya4jMNeTwZR3IQSqQlpQCm+7HiKmCkXhEFQiIqJCMQAqA5riJzhXZXSDbutxBNaV4xM1NTheOgCSwdEZERERFYYBUBlIJzhXml2l8B2JriowCAk1Y4cnrKtQhICUEh48hFkCT0REVBQGQGVASx2BVYYyd4C6jsAAkRUAhXQ1NRneAySC4zEiIiIqDD85y4CmKqkAKDUPrNsOkICE1i0AMjUVMcv1Z4Dp/G0kIiIqBj85y4Sh+cNOAaA1ZgVVYH4n6OwdIFURiIZUdCQdToEnIiLaC/zkLBOGqqAiIweoPVUFVmGqUEV2AAQA1WE91QNI4REYERFRkfjJWSZCuhrsADW3J2G7fqfniKFBKOgRAIVTDRErTBVCsASeiIioGAyAyoShqkE5e7rJoaYIGJrIuQNk6ipCeteuERERERWOAVCZUFWBqKkiM8ypSM0BU3IEQOHUjpGpswcQERFRsRgAlYl0M8R0JRjgV4ClB6Gq3Y65DE1BfYXJHSAiIqK9wACoTKiKgIRAdcQIHqs0/S7QuXaAAOCg4RVB7yAiIiIqHAOgMqGldnmqsnaANHiePwiVic5ERESlwwCoTGiqAkVBMBAVSAVAUsJgnx8iIqKSGvRP1rvuuguTJk1CKBTC9OnT8cILL+S9dsGCBRBC9Pg64ogjsq575JFHMHXqVJimialTp+Kxxx4b6LfRb13jMLJzgFxPQmefHyIiopIa1E/Whx9+GFdddRWuv/56rFq1CieffDJOP/10NDY25rz+jjvuQFNTU/C1efNm1NXV4bzzzguu+cc//oHzzz8f8+fPxz//+U/Mnz8fn//85/HKK6/sq7e1V1RFQFUUVJoZO0CmXwXGTs9ERESlJaSUcrB++MyZM3Hcccfh7rvvDh47/PDDcfbZZ2PJkiV9fv8f/vAHnHPOOdiwYQMmTJgAADj//PPR1taGJ554IrjuU5/6FGpra/Hb3/62oHW1tbWhuroara2tqKqqKvJd7R0pJV7+cDeeWrsd97+0EQDw9TlTcNTYakwdXYVxdZF9sg4iIqL9VTGf34O2tWBZFlauXIl58+ZlPT5v3jy89NJLBb3G0qVLMXfu3CD4AfwdoO6v+clPfrLX10wmk2hra8v62teEENBVgUozOwlaKICSowKMiIiI9t6gBUA7d+6E67oYOXJk1uMjR47E9u3b+/z+pqYmPPHEE7j00kuzHt++fXvRr7lkyRJUV1cHX+PGjSvinZSOoSnBOAzAD4AgkTUJnoiIiPpv0JNLupd3SykLKvm+//77UVNTg7PPPrvfr3ndddehtbU1+Nq8eXNhiy8xU1OCcRiAHwBJ6fcBIiIiotIZtDbC9fX1UFW1x85Mc3Nzjx2c7qSUuO+++zB//nwYhpH13KhRo4p+TdM0YZpmke+g9EK6mtXZucLU4HiSO0BEREQlNmg7QIZhYPr06Vi2bFnW48uWLcPs2bN7/d4VK1Zg/fr1uOSSS3o8N2vWrB6v+fTTT/f5muVAVQQihoqJwyJoqA6hKqRDUQRzgIiIiEpsUAdJLVq0CPPnz8eMGTMwa9Ys3HPPPWhsbMRll10GwD+a2rp1Kx544IGs71u6dClmzpyJadOm9XjNhQsX4pRTTsHNN9+Ms846C3/84x/xzDPP4MUXX9wn76k/NEWBEAI///wxkACkBFTBHCAiIqJSG9QA6Pzzz8euXbuwePFiNDU1Ydq0aXj88ceDqq6mpqYePYFaW1vxyCOP4I477sj5mrNnz8ZDDz2EG264Ad/73vdw0EEH4eGHH8bMmTMH/P30l6oISCmhpRofxi0XipJ7DhgRERHtvUHtA1SuBqMPEAC0dFp4fWMLRlSZUIRAZ9KBIyVmTqpjN2giIqI+7Bd9gKgnVRVQVMD1/JjU9SRUCB6BERERlRgDoDKiKwpUIYIAyJMSmsZJ8ERERKXGAKiMqIq/2+PJrh0gg0dfREREJcdP1zKiKQKqUDJ2gABd5e4PERFRqTEAKiOKIqBpXUdgrpQwOAmeiIio5PjpWmYMVcnKAWL1FxERUenx07XMGFpXAAQJ9gAiIiIaAAyAyoypKXDTrZkEAyAiIqKBwACozOiqAomu3pQMgIiIiEqPAVCZ0VQBgVTQI/35YERERFRa/HQtM+kdH9eTUBUBlU0QiYiISo4BUJlJ7/i4noSi+OMxiIiIqLQYAJUZVRFQFAHb9aAK7gARERENBAZAZUZTBDQhYLkeFEUwCZqIiGgAMAAqM5rabQeIARAREVHJMQAqM5qiQFUEbEdyB4iIiGiAMAAqM+mJ8JbrwdAY/BAREQ0EBkBlyNRUuJ6EwTlgREREA4KfsGUovfPDQahEREQDg5+wZcjQ/N8W7gARERENDH7CliFDVaEqgMYAiIiIaEDwE7YMqaqAoSmsACMiIhogDIDKkK4I6JoCjQEQERHRgGAAVIZURUBXFCgMgIiIiAYEA6AypCkKDO4AERERDRgGQGVIVQX01EgMIiIiKj0GQGXIUJXgi4iIiEpPG+wFUE+GpmDamGoIwR0gIiKigcAthjLF4IeIiGjgMAAiIiKiIYcBEBEREQ05DICIiIhoyGEAREREREMOAyAiIiIachgAERER0ZDDAIiIiIiGHAZARERENOQwACIiIqIhhwEQERERDTkMgIiIiGjIYQBEREREQw4DICIiIhpytMFeQDmSUgIA2traBnklREREVKj053b6c7w3DIByaG9vBwCMGzdukFdCRERExWpvb0d1dXWv1whZSJg0xHieh23btqGyshJCiL1+nba2NowbNw6bN29GVVVVCVdIufB+71u83/sW7/e+xfu9b5Xqfksp0d7ejtGjR0NRes/y4Q5QDoqiYOzYsSV7vaqqKv4PaB/i/d63eL/3Ld7vfYv3e98qxf3ua+cnjUnQRERENOQwACIiIqIhhwHQADJNEzfeeCNM0xzspQwJvN/7Fu/3vsX7vW/xfu9bg3G/mQRNREREQw53gIiIiGjIYQBEREREQw4DICIiIhpyGAARERHRkMMAaADdddddmDRpEkKhEKZPn44XXnhhsJe031uyZAmOP/54VFZWYsSIETj77LPx3nvvZV0jpcQPfvADjB49GuFwGHPmzMHatWsHacUHliVLlkAIgauuuip4jPe7tLZu3YovfelLGDZsGCKRCI455hisXLkyeJ73u3Qcx8ENN9yASZMmIRwOY/LkyVi8eDE8zwuu4f3un+effx6f+cxnMHr0aAgh8Ic//CHr+ULubzKZxBVXXIH6+npEo1GceeaZ2LJlS/8XJ2lAPPTQQ1LXdXnvvffKt99+Wy5cuFBGo1G5adOmwV7afu2Tn/yk/NWvfiXfeustuXr1avnpT39ajh8/XnZ0dATX3HTTTbKyslI+8sgjcs2aNfL888+XDQ0Nsq2tbRBXvv979dVX5cSJE+VRRx0lFy5cGDzO+106u3fvlhMmTJALFiyQr7zyitywYYN85pln5Pr164NreL9L58c//rEcNmyY/Mtf/iI3bNgg/+///k9WVFTI22+/PbiG97t/Hn/8cXn99dfLRx55RAKQjz32WNbzhdzfyy67TI4ZM0YuW7ZMvvHGG/JjH/uYPProo6XjOP1aGwOgAXLCCSfIyy67LOuxww47TF577bWDtKIDU3NzswQgV6xYIaWU0vM8OWrUKHnTTTcF1yQSCVldXS3/67/+a7CWud9rb2+XBx98sFy2bJk89dRTgwCI97u0vvOd78iTTjop7/O836X16U9/Wl588cVZj51zzjnyS1/6kpSS97vUugdAhdzfPXv2SF3X5UMPPRRcs3XrVqkoinzyySf7tR4egQ0Ay7KwcuVKzJs3L+vxefPm4aWXXhqkVR2YWltbAQB1dXUAgA0bNmD79u1Z9940TZx66qm89/3w9a9/HZ/+9Kcxd+7crMd5v0vrT3/6E2bMmIHzzjsPI0aMwLHHHot77703eJ73u7ROOukkPPvss3j//fcBAP/85z/x4osv4l/+5V8A8H4PtELu78qVK2HbdtY1o0ePxrRp0/r9e8BhqANg586dcF0XI0eOzHp85MiR2L59+yCt6sAjpcSiRYtw0kknYdq0aQAQ3N9c937Tpk37fI0HgoceeghvvPEGXnvttR7P8X6X1ocffoi7774bixYtwne/+128+uqruPLKK2GaJi688ELe7xL7zne+g9bWVhx22GFQVRWu6+InP/kJLrjgAgD88z3QCrm/27dvh2EYqK2t7XFNfz9PGQANICFE1q+llD0eo733jW98A2+++SZefPHFHs/x3pfG5s2bsXDhQjz99NMIhUJ5r+P9Lg3P8zBjxgz89Kc/BQAce+yxWLt2Le6++25ceOGFwXW836Xx8MMP49e//jUefPBBHHHEEVi9ejWuuuoqjB49GhdddFFwHe/3wNqb+1uK3wMegQ2A+vp6qKraIzptbm7uEenS3rniiivwpz/9Cc899xzGjh0bPD5q1CgA4L0vkZUrV6K5uRnTp0+HpmnQNA0rVqzAf/zHf0DTtOCe8n6XRkNDA6ZOnZr12OGHH47GxkYA/PNdaldffTWuvfZa/Ou//iuOPPJIzJ8/H9/85jexZMkSALzfA62Q+ztq1ChYloWWlpa81+wtBkADwDAMTJ8+HcuWLct6fNmyZZg9e/YgrerAIKXEN77xDTz66KP429/+hkmTJmU9P2nSJIwaNSrr3luWhRUrVvDe74WPf/zjWLNmDVavXh18zZgxA1/84hexevVqTJ48mfe7hE488cQebR3ef/99TJgwAQD/fJdaLBaDomR/DKqqGpTB834PrELu7/Tp06HretY1TU1NeOutt/r/e9CvFGrKK10Gv3TpUvn222/Lq666SkajUblx48bBXtp+7Wtf+5qsrq6Wy5cvl01NTcFXLBYLrrnppptkdXW1fPTRR+WaNWvkBRdcwLLVEsqsApOS97uUXn31ValpmvzJT34i161bJ3/zm9/ISCQif/3rXwfX8H6XzkUXXSTHjBkTlME/+uijsr6+Xl5zzTXBNbzf/dPe3i5XrVolV61aJQHI2267Ta5atSpoCVPI/b3sssvk2LFj5TPPPCPfeOMNedppp7EMvtzdeeedcsKECdIwDHnccccFpdq09wDk/PrVr34VXON5nrzxxhvlqFGjpGma8pRTTpFr1qwZvEUfYLoHQLzfpfXnP/9ZTps2TZqmKQ877DB5zz33ZD3P+106bW1tcuHChXL8+PEyFArJyZMny+uvv14mk8ngGt7v/nnuuedy/p190UUXSSkLu7/xeFx+4xvfkHV1dTIcDsszzjhDNjY29nttQkop+7eHRERERLR/YQ4QERERDTkMgIiIiGjIYQBEREREQw4DICIiIhpyGAARERHRkMMAiIiIiIYcBkBEREQ05DAAIqJBtWDBApx99tmDvYx+K+R9LF++HEII7NmzZ5+siYjy4zR4IqISuOOOO5DZV3bOnDk45phjcPvttwePzZ49G01NTaiurh6EFRJRJgZAREQlUEhQYxhGMAGbiAYXj8CIqGCe5+Hmm2/GlClTYJomxo8fj5/85CfB82vWrMFpp52GcDiMYcOG4atf/So6OjqC513XxaJFi1BTU4Nhw4bhmmuuQfdpPMlkEldeeSVGjBiBUCiEk046Ca+99lrwfPoY6dlnn8WMGTMQiUQwe/bsHlPUM23cuBFCCDz00EOYPXs2QqEQjjjiCCxfvjzruhUrVuCEE06AaZpoaGjAtddeC8dxgud///vf48gjjwze39y5c9HZ2Qkg+whswYIFWLFiBe644w4IISCEwMaNG3MegT3yyCM44ogjYJomJk6ciJ/97GdZa5o4cSJ++tOf4uKLL0ZlZSXGjx+Pe+65p/ffKCLqW7+niRHRkHHNNdfI2tpaef/998v169fLF154Qd57771SSik7Ozvl6NGj5TnnnCPXrFkjn332WTlp0qRg6KGUUt58882yurpa/v73v5dvv/22vOSSS2RlZaU866yzgmuuvPJKOXr0aPn444/LtWvXyosuukjW1tbKXbt2SSm7hivOnDlTLl++XK5du1aefPLJcvbs2XnXvWHDBglAjh07NvjZl156qaysrJQ7d+6UUkq5ZcsWGYlE5OWXXy7feecd+dhjj8n6+np54403Siml3LZtm9Q0Td52221yw4YN8s0335R33nmnbG9vl1L6k8XT72PPnj1y1qxZ8itf+YpsamqSTU1N0nGcYO0tLS1SSilff/11qSiKXLx4sXzvvffkr371KxkOh7OG+06YMEHW1dXJO++8U65bt04uWbJEKooi33nnnX78ThIRAyAiKkhbW5s0TTMIeLq75557ZG1trezo6Age++tf/yoVRZHbt2+XUkrZ0NAgb7rppuB527bl2LFjg8Cho6ND6rouf/Ob3wTXWJYlR48eLW+55RYpZVcA9Mwzz2T9HAAyHo/nXFs6AMr1s2+++WYppZTf/e535aGHHio9zwuuufPOO2VFRYV0XVeuXLlSApAbN27M+TMyAyAppTz11FPlwoULs67pHgB94QtfkJ/4xCeyrrn66qvl1KlTg19PmDBBfulLXwp+7XmeHDFihLz77rtzroOICsMjMCIqyDvvvINkMomPf/zjeZ8/+uijEY1Gg8dOPPFEeJ6H9957D62trWhqasKsWbOC5zVNw4wZM4Jff/DBB7BtGyeeeGLwmK7rOOGEE/DOO+9k/byjjjoq+O+GhgYAQHNzc6/vIdfPTr/uO++8g1mzZkEIkbX+jo4ObNmyBUcffTQ+/vGP48gjj8R5552He++9Fy0tLb3+vL688847We81/TPXrVsH13WDxzLfqxACo0aN6vO9ElHvGAARUUHC4XCvz0sps4KHTPkez/Uaua7P9dq6rvd4fc/zCvo5udaW62dkrkdVVSxbtgxPPPEEpk6div/8z//EoYceig0bNhT9MzNfP9/PzJT5XtPr2Zv3SkRdGAARUUEOPvhghMNhPPvsszmfnzp1KlavXh0kBQPA3//+dyiKgkMOOQTV1dVoaGjAyy+/HDzvOA5WrlwZ/HrKlCkwDAMvvvhi8Jht23j99ddx+OGH9/s95PrZhx12WLD+l156KSsAeemll1BZWYkxY8YA8AOPE088ET/84Q+xatUqGIaBxx57LOfPMgwjaxcnl6lTp2a91/TPPOSQQ6Cq6l69RyIqDAMgIipIKBTCd77zHVxzzTV44IEH8MEHH+Dll1/G0qVLAQBf/OIXEQqFcNFFF+Gtt97Cc889hyuuuALz58/HyJEjAQALFy7ETTfdhMceewzvvvsuLr/88qyKqGg0iq997Wu4+uqr8eSTT+Ltt9/GV77yFcRiMVxyySX9fg933nln8LO//vWvo6WlBRdffDEA4PLLL8fmzZtxxRVX4N1338Uf//hH3HjjjVi0aBEURcErr7yCn/70p3j99dfR2NiIRx99FB999FHewGzixIl45ZVXsHHjRuzcuTPnjs23vvUtPPvss/jRj36E999/H//zP/+DX/ziF/j2t7/d7/dKRL1jHyAiKtj3vvc9aJqG73//+9i2bRsaGhpw2WWXAQAikQieeuopLFy4EMcffzwikQjOPfdc3HbbbcH3f+tb30JTUxMWLFgARVFw8cUX47Of/SxaW1uDa2666SZ4nof58+ejvb0dM2bMwFNPPYXa2tp+r/+mm27CzTffjFWrVuGggw7CH//4R9TX1wMAxowZg8cffxxXX301jj76aNTV1eGSSy7BDTfcAACoqqrC888/j9tvvx1tbW2YMGECfvazn+H000/P+bO+/e1v46KLLsLUqVMRj8dzHpUdd9xx+N3vfofvf//7+NGPfoSGhgYsXrwYCxYs6Pd7JaLeCZnrwJmI6ACyceNGTJo0CatWrcIxxxwz2MshojLAIzAiIiIachgAERER0ZDDIzAiIiIacrgDREREREMOAyAiIiIachgAERER0ZDDAIiIiIiGHAZARERENOQwACIiIqIhhwEQERERDTkMgIiIiGjIYQBEREREQ87/B+53OWV/DsuJAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import pandas as pd\n", "\n", "plot_df = pd.DataFrame(gene_vecs).melt(var_name='position', value_name='CAI weight')\n", "plot_df['position'] += 1 # skipped the START codon\n", "\n", "sns.lineplot(x='position', y='CAI weight', data=plot_df, errorbar=('ci', 95))\n", "plt.xlabel('codon position')\n", "plt.ylabel('mean CAI weight')\n", "plt.title('metagene plot of CAI weights')" ] } ], "metadata": { "kernelspec": { "display_name": "base", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.8" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }