diff --git a/workflow/notebooks/iqtree_stats.py.ipynb b/workflow/notebooks/iqtree_stats.py.ipynb index f345762d..1eb1610a 100644 --- a/workflow/notebooks/iqtree_stats.py.ipynb +++ b/workflow/notebooks/iqtree_stats.py.ipynb @@ -1605,7 +1605,7 @@ }, { "cell_type": "markdown", - "id": "numeric-encoding", + "id": "advanced-printing", "metadata": {}, "source": [ "# Without Bootstrapping" @@ -1614,7 +1614,7 @@ { "cell_type": "code", "execution_count": 21, - "id": "broad-martial", + "id": "anticipated-intranet", "metadata": {}, "outputs": [ { @@ -1731,7 +1731,7 @@ }, { "cell_type": "markdown", - "id": "professional-salon", + "id": "recorded-injection", "metadata": {}, "source": [ "## Manual Bootstrapping" @@ -3274,7 +3274,7 @@ { "cell_type": "code", "execution_count": 142, - "id": "boolean-forty", + "id": "departmental-michael", "metadata": {}, "outputs": [ { @@ -3547,7 +3547,7 @@ }, { "cell_type": "markdown", - "id": "moderate-integer", + "id": "neither-trademark", "metadata": {}, "source": [ "---\n", @@ -3556,25 +3556,13 @@ }, { "cell_type": "code", - "execution_count": 188, - "id": "global-humanitarian", + "execution_count": 194, + "id": "minor-scientist", "metadata": {}, "outputs": [ - { - "ename": "TypeError", - "evalue": "__init__() got an unexpected keyword argument 'eedgecolor'", - "output_type": "error", - "traceback": [ - "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[0;31mTypeError\u001b[0m Traceback (most recent call last)", - "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[1;32m 38\u001b[0m )\n\u001b[1;32m 39\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 40\u001b[0;31m \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlegend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0meedgecolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"black\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 41\u001b[0m \u001b[0mframe\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlegend\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_frame\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_linewidth\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0.5\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 42\u001b[0m \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_xlim\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m1200\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", - "\u001b[0;32m~/miniconda3/envs/plague-phylogeography/lib/python3.7/site-packages/matplotlib/axes/_axes.py\u001b[0m in \u001b[0;36mlegend\u001b[0;34m(self, *args, **kwargs)\u001b[0m\n\u001b[1;32m 415\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mextra_args\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 416\u001b[0m \u001b[0;32mraise\u001b[0m \u001b[0mTypeError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'legend only accepts two non-keyword arguments'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 417\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlegend_\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmlegend\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mLegend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandles\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabels\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 418\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlegend_\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_remove_method\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_remove_legend\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 419\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlegend_\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", - "\u001b[0;31mTypeError\u001b[0m: __init__() got an unexpected keyword argument 'eedgecolor'" - ] - }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAosAAAFsCAYAAACkQ/OwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAB7CAAAewgFu0HU+AAAyY0lEQVR4nO3deXwddbnH8c+TNElXurcpbai0bC1lp+yyyKYCpRS8XBW5oF6vKFdxu1e87ldEEQsqIupFZBGpLAUqRaAsZVdaERHZKW0DbWkp3Zc06e/+kRMIaSY5SbOcJp/36zWv5Mz8zsxzMoR+MzPPTKSUkCRJkhpT1NkFSJIkqXAZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMiy2QEQcHhEzIuL1iEgRMbmdt/dqbjsNp5+353YlSZLqGBZbpg/wFHBuB21vIjCi3nRsbv6NHbR9SZLUzRkWWyCldGdK6esppVsaWx4RpRFxUUS8FhFrI+LPEXHkVmxvaUppcd0EnAi8DMxu7TolSZJawrDYtq4CDgX+FdiT2iOAf4qInbd2xRFRCpwB/CallLZ2fZIkSfkwLLaRiBgLfBj4UErpoZTSyymli4GHgbPbYBOTgQHAb9tgXZIkSXnp0dkFdCH7AgG8EBH155cBbwJExHuAec2s5+cppcauifwEcGdK6fWtL1WSJCk/hsW2UwTUAPvlvta3Jvf1NWBcM+t5q+GMiBgNHANM2coaJUmSWsSw2HaeBIqBYSmlhxobkFLaBDzXinWfDbwB3NH68iRJklrOsNgCEdEX2KnerB0jYm9geUrphYj4HXBNRHyJ2vA4BHgf8HRKaWYrt1lEbVi8OqVUvVUfQJIkqYXCxtr85W6Dc38ji65OKZ0VESXA14EzgZHUXqv4GPCtlNLTrdzmccBdwK4ppRdasw5JkqTWMixKkiQpk7fOkSRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRl8j6LTYiIMmCP3MulbPlkFkmSpEJRDAzNff90SmljW6zUsNi0PYAnOrsISZKkFpoIzGmLFXkaWpIkSZk8sti0pXXf/OUvf2HEiBGdWYskSVKmRYsWccABB9S9XNrU2JYwLDbt7WsUR4wYwahRozqzFkmSpHy1WZ+Fp6ElSZKUybAoSZKkTJ6G3gqbNm1i9uzZ3Hvvfbz55gpKy0rYffxunHTSSZ6yliRJXYJhsZXuu+8+fvSjn1FZuYLeffagT58dqK7eyIOz7+TKK//A8ccfxte+9lX69OnT2aVKkiS1mmGxFWbOnMk3vvEj+vQ9nMOPPIv+/Xd4e1l19Ubmv3oft8+4gkWLzuOyy35C7969O7FaSZKk1vOaxRZasGAB3/3ujxk0+CQOPewb7wqKAD16lDF2pw9w2Hsv4YknXuNnP/tZJ1UqSZK09QyLLXTrrbeybn0/9tv/s0RE5riBA8cwdqczuH3GvaxcubIDK5QkSWo7hsUWqKqq4tZb76Ki4gSKipo/gz9m7PtZuTJx9913d0B1kiRJbc+w2AJvvvkmK1asY/jwvfMa37Nnf/r0GcuCBQvatzBJkqR2YlhsgZqa2puhR7Tkx1b09vskSZK2NYbFFhg0aBBlZcW8teKVvMZXV29k/fqFDB06tJ0rkyRJah+GxRbo3bs3xx13GK/Om0FKqdnxC+Y/QGnJGo499tgOqE6SJKntGRZb6JRTJpPSQl566Y4mx23cuIrnn/8dRxyxn09zkSRJ2yzDYgvttddenPmxSTz7z5/w7LM3sXlz9RZjVq6czwP3f4nhw9byxS9+oROqlCRJahs+waWFIoLPf/7zlJaWctVVV/DyS79n5Mjj6dtve2pqqnj99UdYu/opdtmlnKlTf8LIkSM7u2RJkqRWi3yuveuuImIUsBBg4cKFW5xOnj9/PrfffjszZ97H8uUrKCkpYcKEnTnttFM44ogjKCkp6YyyJUlSN1RZWUlFRUXdy4qUUmVbrNew2ITmwmJ9KaUmn+giSZLUntorLHrNYhsxKEqSpK7IsChJkqRMBREWI+L8iHgiIlZHxBsRcWtE7JrH+46IiLkRsSEiXomIT3dEvZIkSd1FoXRDHwH8HHiC2pouAO6OiPEppbWNvSEidgRmAr8GzgAOBS6PiKUppZvz2WjumsSmlOdZvyRJUpdUEGExpfT++q8j4mzgDWA/4MGMt30aWJBSOi/3+tmI2B/4MpBXWCTXvCJJkqTGFcRp6Eb0z31d3sSYg4G7G8y7C9g/IrxnjSRJUhsoiCOL9UVtW/FU4OGU0j+aGFoOLGkwbwm1n2kIsCiPzVU0s7yc2lPjkiRJ3VLBhUXgMmBP4LA8xja8SWRkzG/8zc3cf8jb4UiSpO6uoMJiRPwMmAQcnseNJBezZQPKMKAaeLMdypMkSep2CuKaxah1GTAFeF9KaV4eb3sMOLbBvOOAOSmlTW1doyRJUndUEGGR2tvmnAF8BFgdEeW5qVfdgIi4MCKuqfeeK4DRETE1IsZFxMeBTwAXd2jlkiRJXVihhMVzqO2AfoDaxpS66fR6Y0YAO9S9yB19/CBwJPA34BvA5/K9x6IkSZKaVxDXLKaUmu0kSSmd1ci82cC+7VGTJEmSCufIoiRJkgqQYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSpoIIixFxeETMiIjXIyJFxORmxh+ZG9dw2q2DSpYkSeoWenR2ATl9gKeAq4CbW/C+XYFV9V4vbcuiJEmSuruCCIsppTuBOwEioiVvfSOltKK1242IUc0MKW/tuiVJkrqCggiLW+HJiOgJ/BP4Xkrp/ha+f2E71CRJktRlFMQ1i62wCPgUcCowBXgeuDciDu/UqiRJkrqYbfLIYkrpeWoDYp3HIqIC+DLwYAtWVdHM8nLgiRaWJ0mS1GVsk2Exw+PAGS15Q0qpsqnlLbx+UpIkqcvZVk9DN2Yfak9PS5IkqY0UxJHFiOgL7FRv1o4RsTewPKW0ICIuBEamlM7MjT8PeBV4Biil9ojiqblJkiRJbaQgwiKwP1C/k3lq7uvVwFnACGCHestLgYuBkcB6akPjCSmlme1eqSRJUjdSEGExpfQAkHmBYErprAavLwIuat+qJEmS1JWuWZQkSVIbMyxKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJylQQYTEiDo+IGRHxekSkiJicx3uOiIi5EbEhIl6JiE93QKmSJEndSkGERaAP8BRwbj6DI2JHYCbwELAP8H3gpxFxartVKEmS1A316OwCAFJKdwJ3AkREPm/5NLAgpXRe7vWzEbE/8GXg5ny3GxGjmhlSnu+6JEmSuqKCCIutcDBwd4N5dwGfiIiSlNKmPNezsG3LkiRJ6loK5TR0S5UDSxrMW0Jt+B3S8eVIkiR1TdvqkUWA1OB1ZMxvSkUzy8uBJ1qwPkmSpC5lWw2Li9nyesJhQDXwZr4rSSlVNrU8z+snJUmSuqxt9TT0Y8CxDeYdB8xpwfWKkiRJakZBhMWI6BsRe0fE3rlZO+Ze75BbfmFEXFPvLVcAoyNiakSMi4iPA58ALu7YyiVJkrq2QjkNvT9wf73XU3NfrwbOAkYAO9QtTCnNi4gPApcAnwVeBz6XUsr7tjmSJElqXkGExZTSA7zToNLY8rMamTcb2Lf9qpIkSVJBnIaWJElSYTIsSpIkKZNhUZIkSZkMi5IkScpkWJQkSVImw6IkSZIyGRYlSZKUybAoSZKkTIZFSZKkbVxlZSXXXHNN8wNboSCe4CJJkqSWW7VqFRde8EPuu/NRqpaXtMs2DIuSJEnboNWrV/OZ//gcLz26gsMHnM+wUTvxt9fGtvl2DIuSJEnboIt/9GNeemwFp29/GYN6jmL5hsp22Y5hUVKH2LhxI2+99RYAAwcOpKysrJMrkqRt15IlS7h7xoMc1O8LDOo5ql23ZViU1K5efPFFpk+fzr33zGTDhnUA9OzZm2OPP5HJkyez0047dXKFkrTtmTlzJtVv9WL3HY5p920ZFiW1i5QS119/Pb/65U8ZvF0Rk4/bjV12Kgfg+RcXcc/9NzHj9hv59Dnn8eEPf7iTq5Wkbcu8efMYynhKi3u1+7YMi5LaxU033cSvfnEJp504ntMnH0Rx8Tt36tp7wmhOm3QAN0x/nCt+/mPKysqYMmVKJ1YrSduWmpoailJxh2zL+yxKanMrV67kil/8lBOP2ZmPnHrIu4JineLiIj562iF84KixXP7zS1m1alUnVCpJ26ahQ4eyIl4lpdTu2zIsttCaNWu4+eab+cqXvsQ5n/okXzzvPK677rq3L9yXBHfeeSepZi0fOvmAZsf+y+QDqdm0mj/96U8dUJkkdQ3HHXccVX2W8Mrqv7T7tgyLeUopce2113LapA9yxUXfpsf8xxlTNY++i+Zy3WU/5F9O/iA/++lPqamp6exSpU5376y7OWjfUfTr2/y1NP23682B+4zi3ll3d0BlktQ1jBs3jj0m7szDy66kqmZDu27LaxbzdNVVV3H/jJv414mj+NAhExnU751/BFev38gfn3iJK6/7FcuWLeVb3/4ORUXmcHVfy99cyl47D8p7/PblA3nlr8vasSJJ2vZs2rSJlStXvj2tWLHiXV979yvjxU0z+M6T+zG49D2sq17RLnUYFvN05/RpnP/B3Zl0wC5bLOvXq4wPH747Owzdjm/ccivT99qbU089tROqlApDaVkZG6uq8x6/sWoTpaXed1FS15FSYs2aNY2GvIZfs5atW7cuc/29e/emf//+9OrVi+XrFrFq/WJ6M7hdPothMU/jhpQ2GhTrO3RcBUfvsoBb/nADU6ZMISI6qDqpsOy62wTmPnUv/3Z6avb3IKXE3KcWMm6vYzuoOklq3saNG/MKdFnLVq5cyebNmxtdd48ePejfvz8DBgx419fy8vK3Xzdc1vD7kpJ3ngO9bt06Zs2axTW/vZ6bp7/Y5j8Lw2KePrD3jnmNm3zgLsy69s/89a9/Zb/99mvnqqTCdPLJJ3PerD/y9LOV7Dm+osmxTz2zgEVvbOD8yZM7pjhJXd7mzZtZtWrVVoW9DRuyrwPs27fvFoFu5MiRjB8/fouQ11jo6927d5seUOrduzeTJk1i33335ebp09psvXUMi3kaXzE0z3FD6JE2sWDBAsOiuq29996bXcbtyc/+7z6+/z+nMHTwdo2OW/rmKi678n7GTdiXPffcs4OrlFSIUkps2LAh71O1jY1ZvXp15i1lSktLGw1yFRUVmUfx6gfD7bbbjuLijrm/YaEwLOapuDi/vwAiguKiIrui1a1FBBdc8APO/eyn+e/v3Mzpp+zP4QfvSq+epQCs31DF7Eef4w+3zaVXvwr+938v8LINqYuorq7eIsC19Ajfpk2bGl13RLDddtttEfbe8573vCvQNXVkr2fPnh38E9n2GRbz9OobK9lt1JBmx1UuW8XGzTBs2LAOqEoqXMOGDeMXV/yaqVN/zK9+dz9XT/szFdv3h4CFr61kY3UJ7z38WL7whS8yaFD+ndOS2k8+TRnNLVu7dm3m+nv37r1FoBsyZAhjx4591/zGjugNGDCAvn37ereRTmBYzNPdT83j/fuObXbc7X95gX5DRnDQQQd1QFVSYRs8eDAXXPB9lixZwl133cXrr78OwOHHbM/73/9+/6iS2lhdU8bWhL18mzLqpvLy8iabMbKaMrTtMCzm6c+vruCpeUvYa8fhmWNeXbKCO55+nUn/9llKS0s7sDqpsA0fPpwzzzyzs8uQClr9pozWhr2mmjL69eu3RXhrrimj/vdt3ZShbUfBhMWI+AzwFWAE8AxwXkrpoYyxRwL3N7JoXErpufaob/x+h3L+tMf56kn78N7xFe/6hUkp8bd5S/juzX+hfJe9+djHPtYeJUiSClTDpozWhL2WNmX079+fioqKJpsx6l53x6YMtZ3oiAdQN1tExOnAtcBngEeA/wA+CYxPKS1oZPyR1IbFXYFV9RYtTSm1WWdJRIwCFgK88MIL/PY3v+HR+/7EyL7BcXuMYlDfXqxat5H7/vkaLy/fxO77H8T3vv8DBgwY0FYlSJI6QF1TRmuaMVrblJFPM0bd1KtX84/OlCorK6moePt2ZRUppcq2WG+hhMU/A39NKZ1Tb96zwK0ppfMbGX8ktWFxYEppxVZsd1QzQ8qBJwAWLlzIyJEjeeaZZ7h1+nQef/gB1q5ZQ6/evdln4kGcMuVU9ttvPw/RS1IHa6opI9+w11xTRlONF1nNGHXz+/XrZ1OGOkR7hcVOPw0dEaXAfsAPGiy6Gzikmbc/GRE9gX8C30spNXZquikLWzI4IpgwYQITJkwAvvGu/0H17NmTlJp/WoUk6d2aasrIN+w115TR8Khd/aaM5o7s2ZSh7q7TwyIwBCgGljSYv4TaI3uNWQR8CpgLlAEfA+6NiCNTSg+2V6F1ampqePjhh7n+2mt4YNYslm/YwKYUlPXsyb57TOBzn/0sxx13nNeHSOryNm/ezOrVq1vdebtixYq8mjLqh7qGTRlNndq1KUPaeoUQFus0PB8ejcyrHZjS88Dz9WY9FhEVwJeBloTFpp9DVu80dJ21a9fy9a+dz50338RbNZvpMaycwcccQtngwaxZt4HH//EUj55zLofuNYFfX345I0eObEE5ktRx6poyWtt525KmjPpBrmFTRlbosylDKgyFEBaXATVseRRxGFsebWzK48AZLdlwc+fyG/41Wl1dzTe+9jVm3TSNjb37scPk0xl1+JEUl5a9M+aDJ/Dcs8/zyE3Xc8YnPsm0a672XnKS2kV1dTWrVq3aqg7cljZljB49mr322iuva/d8UobUNXR6WEwpVUXEXOBYYHq9RccCt7VgVftQe3q63dx77708PHMGG8t6Mfz0Mxn13iO3GNOjqIjx43cjzj6Hp6/+JRdPncpFP2h4Oaak7q7+Nc/5PvO24es1a9Zkrr9Xr15bhLfBgwczZsyYvLpybcqQVKfTw2LOVODaiJgDPEbt9Yg7AFcARMSFwMiU0pm51+cBr1J7P8ZSao8onpqb2s2tN91I6cZ1pDHjGHnYEZnjiiLYfthgVhxyBDNn38MXFi1ixIgR7VmapA5WVVXV6maMFStWsGrVqsxnyBcXFzca6IYPH95s523d6VsfDCCprRREWEwpTYuIwcA3qb0p9z+AD6aU5ueGjKA2PNYpBS4GRgLrqQ2NJ6SUZrZXjW+88QbPzHmC5UU9GHzoEc1eMD2grJS+4/dk+YP3cM899/j0CqmAZDVltCTsNdWU0bdv3y0C3fbbb8+4ceOa7Lz1SRmSClFBhEWAlNLlwOUZy85q8Poi4KIOKOttq1atYtPGDWzuUULf7Zu7PWPt9T5lPctYN2wEixcv7oAKpe6hqaaMfE/nrlq1qsmmjMYCXNaTMhoe6bMpQ1JXUzBhsdD16NGDKIrMf2AaszlBJLzuR6qnrilja8JeVVVVo+uua8poGPZGjx7NnnvuuUUzRmNH9mzKkKR3Myzmqby8nO2GDGPZyy+x4pWXGbDTLk2O31Bdw4a1a4mlixg1qvkjkdK2IN+mjKYemZZPU0b9QJfVlNHYtXs2ZUhS2zMs5qlnz558cMpp/PoH32PFo7PZfPSxFBVn//iWrd/AxmefZmSvMo4//vhm179x40buv/9+Zt11F8tef52S0lLGjBvHpJNPZvz48V6/pDZRvymjNbdZaa4po7Fbqey8885NBr36X23KkKTCY1hsgVNOOYU/XPNbls+fx7w7/8iYE05uNMSt3FjF4tcXs/mx2Zx8wvEMHDiwyfXOnj2bH3/3u6x+9VX2r6lh/5ISqlNizkMPcfd117HbIYfwne9/3/s1dnMNmzJaE/aaa8poGN5GjBjRbFOGT8qQpK7NsNgCFRUVXPDjS/jPf/8ElbdMY93q1ez8gRPpOaA2DFbV1LB09ToW/P3vVN17B8fsMpbPfe5zTa5z1qxZXPClL3Hk6tX8e3k525e9c4PvlBJzVq9m6j33cO6SJfz8N79h6NCh7foZ1T4aa8poadhrqimjpKSk0VupjBo1KjPc1X+93Xbb0aOH/zuQJG0pWtKw0d1ExChgIcDChQvfvvZw7ty5fPUrX2bOP55hc+++9NxzH4oHDaWqqopNzz9D/6r1TDnuWL7zzW/Sv3//zPUvXbqUj5x4IscsW8Z/VVRkHpV5c9MmPrNwITtMmcKPLrmk7T+omlVTU9PkdXj5hL6WNmU09Ri0hst69uzpUb12UlNTw6OPPsr999/P8rdW0LOslPHjx3PCCScwePDgzi5Pkt5WWVlJRcXbTzKuaO5JdfkyLDYhKyxC7qjfnDlceeWV/P0f/2BDdQ2DBw/i2COO4NRTT2Xs2LHNrv83v/kNN//v/3JjRQW9m7nVxr1vvcX3Nm/mtzNmMHr06K37YN1MSom1a9fm1XyRFfZa8qSMfB6D5pMytg0PPvggF13yU156bRmbh+xEj/7lbN60kerKvzOwpJpTTzyOL5x3HmX1zghIUmdpr7DoeadWiggmTpzIxIkTW/X+lBJ33HgjxxYXNxsUAQ7v358B8+czc+ZMzjnnnFZtc1tVVVXVquvzWvqkjPoBrn5TRnNH9mzK6Jr+9Kc/8dXv/ICqHQ7iPR//Nn1HvPMHYPWGtSz56z383y1XsaDyNS798cX+dyCpyzIsdpL169ezfNEi9uzbN6/xJUVFjEuJyoUL27mytlXXlNHasLdy5UrWr1+fuf76TRl1U2NNGVmnb/v06ePpW21h/vz5fOOCH5F2ez/jTv7cFv+N9OjZh5GHTKbfqF2453fn86tf/Ypzzz23k6qVpPZlWOwkdaf/WxJTIlp2U/CtVb8po7VhryVNGXVBrrGmjMZCn00Zai/Tp0/nrejHnid9tsk/JrbbYTz9DziVP9x6Kx//+Mfp3bt3B1YpSR3Df2k7Se/evek3aBAvvPYaRzVzax2AzSnxEnBIeXne26h7UkZT4a650NeaJ2Xke+2eTRkqRBs3buSWO+5iwN6Tm7yXap0R+3+A5x67gVmzZjFp0qQOqFCSOpZhsZNEBO8/7TRmXnQRZ2/eTGmDBoeUEmtralhZU8OK6moeWrGCZ1PiwJoarrjiiry6cPNtyqgLb4MGDWLHHXfMqyvXpgx1VYsXL+at1espH7tPXuPL+g+lePBoXn311fYtTJI6iWGxndQ1ZdQFuMaO3r322mv8ZfVqDpw7lwElJaysrmZldTUrcl8ba8n40pe+RFFRUWZTRnOnbuu+92J8qXE1NTWkBFHUfOPZ24qKqK6ubr+iJKkTGRbzNGvWLEpKSvK+r15zTRlvH9UbOpQFlZUUV1ezb9++DCkpoX+PHgzITes3b+a2FStYOXo0P/jJT5gwYYJNGVI7GjRoEGU9grVvzGe7HcY3O75m00ZqVixmyJAj2784SeoEhsU8nX322cA7TRkNj9iNHDmyyXvp1X3fWFPGLbfcwi9+8AMqly5l56IiRpeVUQM8vn49f+/Rg+EHHsgll1zCmDFjOuGTS93LgAEDeN+hB3D73JmM2P8DzY5f+vRs+rKeY445pgOqk6SOZ1jM05w5cxg/fny7NGVMmTKFI488kpkzZ3LPjBk8uGgRJSUljNl9d741ZQqHHXaYXb9SBzp1yinMnP1V3nh6NsP2OCJz3Kb1q1ny8DROPnQi22+/fQdWKEkdxye4NKGpJ7hI6rpSSlxwwQVcPeN+hh//nwzf51iiQUPX+uWv8+IfLmB0LOOqX13u/x8kdTqf4CJJHSQi+OpXv0pxcTG/v+1Sljx0PQP2Oo5eg7anpmo9bz33ONXz57BbxRAuvfhSg6KkLs0ji03wyKLUvaWUeP7555k+fTp3zHqAdes3EgF77LYz/3LqKRx99NE+F1pSwWivI4uGxSYYFiXVSSmxceNGSkpKKM7jee6S1NE8DS1JnSgi6NmzZ2eXIUkdzkdwSJIkKZNhUZIkSZkMi5IkScpkWJQkSVImw6IkSZIyGRYlSZKUybAoSZKkTAUTFiPiMxExLyI2RMTciHhvM+OPyI3bEBGvRMSnO6pWSZKk7qIgwmJEnA5cClwA7AM8BNwZETtkjN8RmJkbtw/wfeCnEXFqhxQsSZLUTRTKE1y+CFyZUvq/3OvzIuJ44Bzg/EbGfxpYkFI6L/f62YjYH/gycHO+G809zq8p5fmuS5IkqSvq9COLEVEK7Afc3WDR3cAhGW87uJHxdwH7R0RJCza/sJnpibqB37/wQu677z6qq6tbsHpJkqRtW6eHRWAIUAwsaTB/CdlH9sozxvfIra/NPbZ6CZ//7jc4+UOn8cgjj7THJiRJkgpOIYTFOqnB62hkXnPjG5vflIpmpol1Aw/7xuc54pcXsHKX7Tnv6+dz//33t2AzkiRJ26ZCCIvLgBq2PIo4jC2PHtZZnDG+Gngz3w2nlCqbmnLbeVv/HUZy+LfOo++R+/ON73+PN954I99NSZIkbZM6PSymlKqAucCxDRYdCzya8bbHGhl/HDAnpbSpbSt8t4hg4rlnsbI4cccdd+T1nldeeYW77rqLP/7xjzz66KNUVVW1Z4mSJEltplC6oacC10bEHGqD4KeAHYArACLiQmBkSunM3PgrgHMjYirwa2obXj4BfLgjii3p3Yvyow/m5hm3c9ZZZxERjY578MEH+d0N03jsr8+wdlOCCEoCRg3djg+dfAIf+chH6NevX0eULEmS1CoFERZTStMiYjDwTWAE8A/ggyml+bkhI6gNj3Xj50XEB4FLgM8CrwOfSynlfducrTV8z3E8f8dDrF27lr59+75rWUqJK6+8kkt/fR3V2+9F+aSvs9MuE4keJaxf9hqL5szk4qtv497ZD3PZpT9m2LBhHVW2JElSixREWARIKV0OXJ6x7KxG5s0G9m3nsjJFcREJ2Lx58xbLbrvtNqb++jq2O/yTjDrs3fcJ7z10FGM/8CnWH3giT179Vb7w5f/mqv/7JaWlpR1UuSRJUv46/ZrFbdVbL82nX69eWxxV3LRpE7+48rcU73bsFkGxvl6Dtmfsh7/DE8/O44EHHmjnaiVJklrHsNgKm6urWfinBzn5+A9QVPTuH+GDDz7Iq4tXMOqw05pdT9/yHSmu2Ic/3HxLe5UqSZK0VQyLrfDMtBmUrFjLpEmTtlg2Z84cGDKGPsMafaz1Fobs+T7++vRzrF27tq3LlCRJ2mqGxRao3rCRp66+ifnX3c5/nv0JxowZs8WYdevWEb36573Okj7bUZMwLEqSpIJUMA0uhe6pa27irbn/oNf6ar7yyU/z0Y9+tNFxffr0Ia1fkPd6N61dRXHUvk+SJKnQeGQxT72efInPT/kwt/9+GmeccUbmvRUnTpwIy+axdsn8Rpc3tOzv97H/XuMNi5IkqSAZFvN05RW/5JOf/CTl5Q2fMvhu733vexlTPpDKR25qdp1rFr1MzcInOe2UyW1UpSRJUtsyLLaxHj168JlPfZzNz81iwexppJQaHbduWSUv3fBtDtx9DEcddVQHVylJkpQfr1lsByeddBLLli1j6i9/yzMvPcHwAyYxqP4TXJ6Yydpn7mHi2OFccvFFlJSUdHbJkiRJjTIstpOzzz6bXXbZhd/dMI1H7riQ126tfQxgSXEwetgAPvTxUzn99NN9NrQkSSpohsV2dOihh3LooYcyf/58XnrpJaqqqhg0aBD77ruvRxMlSdI2wbDYAUaPHs3o0aM7uwxJkqQWs8FFkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpU0GExYgYGBHXRsTK3HRtRAxo5j2/jYjUYHq8g0qWJEnqFnp0dgE51wOjgPfnXv8KuBY4qZn3/Qk4u97rqrYvTZIkqfvq9LAYEeOoDYkHpZT+nJv378BjEbFrSun5Jt6+MaW0eCu2PaqZIeWtXbckSVJX0OlhETgYWFkXFAFSSo9HxErgEKCpsHhkRLwBrABmA/+TUnqjBdte2Ip6JUmSuo1CuGaxHGgs4L1B00f27gQ+CrwP+BIwEbgvIsravEJJkqRuqt2OLEbEt4FvNTNsYu5ramwVGfNr35DStHov/xERc4D5wAnALXmWWdHM8nLgiTzXJUmS1OW052noy4AbmhnzKrAnMLyRZUOBJfluLKW0KCLmAzu34D2VTS2PiHxXJUmS1CW1W1hMKS0DljU3LiIeA/pHxAEppb/k5h0I9AcezXd7ETGY2iOFi1pXsSRJkhrq9GsWU0rPUnsLnF9HxEERcRDwa+CP9TuhI+K5iDgl933fiLg4Ig6OiPdExJHADGrD6fQO/xCSJEldVKeHxZyPAk8Dd+emvwMfazBmV2qPNgLUAHsAtwEvAFfnvh6cUlrdEQVLkiR1B4Vw6xxSSsuBM5oZE/W+Xw8c3951SZIkdXeFcmRRkiRJBciwKEmSpEyGRUmSJGUqiGsW1byUEs899xwvvPACmzZtYuDAgRx88MH07t27s0uTJEldmGFxG3Dvvffy+xuu558vPsPmqKZHSQ+qq6rp17M/Hzz+RM4880wGDx7c2WVKkqQuyLBYwFJK/PKXv+SaG65ip/0r+Ph3TmXc/jtRVFTEimWrePRPc5gx42YeeexhLp36E0aNGtXZJUuSpC4mUsp8/HK3FxGjgIUACxcu7PAwNn36dC76yYWc9B9HcdTkQxods2LZKn5+/jX0qRnIVVf+1tPSkiR1U5WVlVRUVNS9rGjuscb5ssGlQFVXV3P1tVex33HjM4MiwIAh2/Gp73yEBYvncc8993RghZIkqTswLBaohx9+mMXLF/G+KdlBsc7Q7Qcx7qCx3HLrzXikWJIktSXDYoGaO3cuw0cPYsTo4XmN3/+oPXhp3ousXLmynSuTJEndiWGxQK1fv55e/XrmPb53v94kEuvWrWvHqiRJUndjWCxQffr0Yc2K/IPf6rdWU0QRffr0aceqJElSd2NYLFAHHXQQb1au5NXn82tk+vOsp9h9twn079+/nSuTJEndiWGxQB144IGMGr4Ds258qNmmlcqXF/HSXxcw5ZRTO6g6SZLUXRgWC1RRURGf/tQ5PPvIq9z2m7szA+OiBW/w62/9nvFjJ3DUUUd1cJWSJKmr8wkuBezoo49m+fLlXPrzqbzw13kceuL+7HHgrpT2LGVJ5TIeufMJnrzvWcaO3IUfXXQxpaWlnV2yJEnqYnyCSxM6+wkudZ588kmmTbuBR/78MJtSFZAIiigfvD2TJ53CaaedZmOLJEndXHs9wcUji9uAffbZh3322YfFixfz8ssvU1VVxcCBA9ljjz0oLi7u7PIkSVIXZljchpSXl1NeXt7ZZUiSpG7EBhdJkiRlMixKkiQpk2FRkiRJmQyLkiRJymRYlCRJUia7oTtJdXU1jzzyCI8++ihr1qyhV69e7L333hxzzDH07Nmzs8uTJEkCvCl3k9rrptx33XUXl//6CirfXMyAncvpOagfm9asZ9k/Kxncuz8f+5eP8LGPfYyIaJPtSZKkrs+bcncR06ZN48dX/ISKo/fgXz/yUYaNGfn2shWL32TuLQ9w6VWXs3jJYr7y5a8YGCVJUqcyLHagv/3tb0y94qdM+OgRHP6Jk7YIggPKB3P0Z05l+E6jmHbhNHbbdTcmTZrUSdVKkiQVSINLRPxPRDwaEesiYkWe74mI+HZEvB4R6yPigYjYvZ1L3Sp/uPEP9B4zpNGgWN+E4w5k9FET+P2NN+BlApIkqTMVRFgESoEbgV+04D3/BXwROBeYCCwG7omIfm1f3tZbunQp9z/2EHuf8t68Ti3vO/lwXqp8lSeffLIDqpMkSWpcQZyGTil9CyAizspnfNSmrfOAC1JKt+Tm/RuwBPgI8Ms819Ncx0qbPYh53rx5VKVqxhwwLq/xo/YYCz2Lefnll9l3333bqgxJkqQWKYiw2Ao7Uhvk7q6bkVLaGBGzgUPIMyyS63TuCNXV1SQSxSX5/cgjgqKSHlRXV7dzZZIkSdkK5TR0S9Ud8VvSYP4S2vBoYFsaNGgQxRTx5oKGJTduzfJVbFqzgUGDBrVzZZIkSdnaLSzmmk9SM9P+W7mZht0f0ci8plQ0M03cyvretuuuu7LTqPfwtxmP5DX+qTseoX9pHw499NC2KkGSJKnF2vM09GXADc2MebWV616c+1oOLKo3fxhbHm3M1NzNKtvyHocRwZRJp/DDX13Kog/NZ8SuozPHrlr6Fn+/5WFOPeYD9O3bt81qkCRJaql2O7KYUlqWUnqumWlDK1c/j9rAeGzdjIgoBY4AHm2D8tvFpEmT2G/nPZj+X79kwVMvNjrmzQWLufGLlzGydBAfP/vjHVyhJEnSuxVEg0tE7AAMAnYAiiNi79yil1JKa3JjngPOTylNTymliLgU+FpEvAi8CHwNWAdc39H156usrIyLf/gjvvb1/+G2837BoPGjGH/8RPoOGcCG1Wt5/v4nef0vL7Lz9u/hx1MvZvDgwZ1dsiRJ6uYKIiwC3wX+rd7rupsLHgU8kPt+V6B/vTEXAb2Ay4GBwJ+B41JKq9u10q3Ur18/Lp16CY899hjTb5vOo5fMYDObKSLYfZdxfOYr3+R973sfZWVlnV2qJEkS4RNCsuXuw7gQYOHChYwa1dxtGVtu48aNrF27ll69etGrV682X78kSeoeKisrqaioqHtZ0VxvRr4K5chit1VWVuZRREmSVLC21fssSpIkqQMYFiVJkpTJ09BNK677ZtGiRU2NkyRJ6lQNskpx1riWssGlCbknzDzR2XVIkiS10MSU0py2WJGnoSVJkpTJI4tNiIgyYI/cy6VATSeW05WU884R24m88/hGFRb307bB/bRtcD9tG7b1/VQMDM19/3RKaWNbrNRrFpuQ+yG3ySFcvaPBM7cXt9V9oNS23E/bBvfTtsH9tG3oIvtpfluv0NPQkiRJymRYlCRJUibDoiRJkjIZFiVJkpTJsChJkqRMhkVJkiRlMixKkiQpkzflliRJUiaPLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDItqMxFxfkQ8ERGrI+KNiLg1InZtMCYi4tsR8XpErI+IByJi9wZjyiLiZxGxLCLWRsTtETGqYz9N95DbZykiLq03z31UICJiZERcFxFvRsS6iPhbROxXb7n7qpNFRI+I+F5EzMvtg1ci4psRUVRvjPupg0XE4RExI/czTxExucHyNtknETEwIq6NiJW56dqIGND+n7BjGRbVlo4Afg4cBBwL9ADujog+9cb8F/BF4FxgIrAYuCci+tUbcylwCvCvwGFAX+CPEVHc3h+gO4mIicCngL83WOQ+KgARMRB4BNgEfAAYD3wJWFFvmPuq8/038Glq98E4avfJV4D/rDfG/dTx+gBPUfszb0xb7ZPrgb2B9+emvYFr2+gzFI6UkpNTu0zAUCABh+deB7AI+O96Y8qo/cfvP3Kv+wNVwOn1xmwP1ADHd/Zn6ioTtf/TewE4BngAuNR9VFgT8APgoSaWu68KYAL+CFzZYN7NwLXup8KYcv8OTa73uk32CbV/HCTgwHpjDsrN27WzP3dbTh5ZVHvqn/u6PPd1R6AcuLtuQEppIzAbOCQ3az+gpMGY14F/1Bujrfdz4I6U0qwG891HhWMSMCcibsxd1vFkRPx7veXuq8LwMHB0ROwCEBF7UXsUamZuufup8LTVPjkYWJlS+nO9MY8DK+li+61HZxegrikiApgKPJxS+kdudnnu65IGw5cAo+uNqUopvdXImHK01SLiX4F9qT310pD7qHCMAc6h9vfo+8ABwE8jYmNK6RrcV4Xih9T+YfxcRNQAxcD/pJR+n1vufio8bbVPyoE3Gln/G3Sx/WZYVHu5DNiT2r+wG0oNXkcj8xrKZ4yaEREVwE+A41JKG5oY6j7qfEXAnJTS13Kvn8xdgH8OcE29ce6rznU6cAbwEeAZaq9ZuzQiXk8pXV1vnPup8LTFPmlsfJfbb56GVpuLiJ9RewrtqJRSZb1Fi3NfG/7FNYx3/sJbDJTmLu7PGqPW24/an+XciKiOiGpqG5M+l/u+7mfsPup8i4B/Npj3LLBD7nt/nwrDj4AfpJRuSCk9nVK6FrgEOD+33P1UeNpqnywGhjey/qF0sf1mWFSbyd2K4DJgCvC+lNK8BkPmUfvLdWy995RSG1Yezc2aS233Z/0xI4AJ9cao9e4F9qD26EfdNAf4Xe77V3AfFYpHgF0bzNsFmJ/73t+nwtAb2NxgXg3v/Pvqfio8bbVPHgP6R8QB9cYcSO1lCV1rv3V2h41T15mAy6ntJjuC2r/Y6qZe9cb8d27MKdT+0l0PvA70qzfmF8BC4GhgH2oDzt+A4s7+jF1xol43tPuocCZqryndBHwN2Ina05xrgY+6rwpnAn4LVAInAO/J7YulwA/dT526X/ryzh/ECfhC7vsd2nKfAHdSe4ueg3LT34EZnf352/zn2dkFOHWdKfcL2dh0Vr0xAXyb2lNsG6jtPpvQYD09gZ8BbwLrgBlARWd/vq46sWVYdB8VyAScCDyd2w/PAv/eYLn7qvP3UT9q78c3H1gPvAx8Dyh1P3Xqfjky49+j37blPgEGAdcBq3LTdcCAzv78bT1F7sNKkiRJW/CaRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSJsOiJEmSMhkWJUmSlMmwKEmSpEyGRUmSJGUyLEqSJCmTYVGSJEmZDIuSJEnKZFiUJElSpv8HrhEMuPSJobwAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAGFCAYAAAD+TDXqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAB7CAAAewgFu0HU+AABarUlEQVR4nO3dd3xUVfrH8c+ThITepVdBKRaqBV0BC9gF7L0rYlvBdRW7rqx1lbWtP921rm1VmgooqNgoSlMsgCJIL0oNECDJ8/vjTnCYTJJJMpMJ5Pt+ve4rc88999xn7lCenDn3HHN3REREREQqupRkByAiIiIiUh4oMRYRERERQYmxiIiIiAigxFhEREREBFBiLCIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEREREBFBiLCIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEREREBFBivMcxs55m9q6ZLTczN7P+Cb7eotB1IrenEnldERERkXhTYrznqQZ8A1xbRtc7CGgctvUJlb9VRtcXERERiQslxnsYdx/n7re7+4hox80s3cweMrNlZrbZzKaZWe9SXG+Nu6/M24CTgAXApyVtU0RERCQZlBhXPC8AhwNnAwcS9OyON7N9StuwmaUD5wPPu7uXtj0RERGRsqTEuAIxszbAOcAZ7v65uy9w90eAL4BL4nCJ/kBt4MU4tCUiIiJSptKSHYCUqa6AAfPNLLw8A/gdwMxaAQuLaOcpd482hvkyYJy7Ly99qCIiIiJlS4lxxZIC5ADdQj/DZYZ+LgM6FNHOusgCM2sJHAOcWsoYRURERJJCiXHFMgtIBRq4++fRKrj7DmBuCdq+BFgNvF/y8ERERESSR4nxHsbMqgNtw4pam1lnYK27zzezV4GXzexGgkS5PnAUMMfdx5bwmikEifFL7p5dqjcgIiIikiSmyQP2LKGp1z6Jcugld7/YzCoBtwMXAk0JxhZPAe5y9zklvGZf4AOgnbvPL0kbIiIiIsmmxFhEREREBE3XJiIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQAt87DHMrB5wLLAIyEpuNCIiIiKFqgy0Aj5w99+THMtOSoz3HMcCryY7CBEREZFiOA94LdlB5FFivOdYBPDf//6XDh06JDkUERERkYL9+OOPnH/++RDKX8oLJcZ7jiyADh060LVr12THIiIiIhKLcjX8Uw/fiYiIiIigxFhEREREBNBQConBsmXLGDNmDJ9/Po0NGzKpXr0qPXp0pV+/frRu3TrZ4YmIiIjEhRJjKVB2djbDhw/njTfeZ9v26jRs2IvKleuwdv0mXnjxM155ZTQnndSLW28dSkZGRrLDFRGRQmRmZrJ8+XJyc3OTHYrs4VJSUqhduzb169cnJWX3GpygxFiiys3N5Z577mXkqMm0a3ctbdoeT1paRtjxQSxa9DHvjPgn69ffwiOPPESlSpWSGLGIiBRk6tSpDBkyhO3btyc7FKlAunXrxp133knTpk2THUrMlBhLVOPHj2f06C/o2vVumjU/PN/xlJQ09t67L9WrNeLjj2/izTffzJt2RUREypHMzEyGDBlC9+7dueKKK9SJIQmXk5PD0qVLeeqppzjnnHOYOHEi6enpyQ4rJkqMJR93580336FmrUOiJsXhGjQ8kAYN+/LWW2M455xzSE1NLaMoRUQkFsuXL2f79u1cccUVHHDAAckORyqI/fbbj4YNG3L55ZezePFi2rZtm+yQYrJ7DfyQMvHTTz/x7ZxfaNv2lJjq77PvKSxcuIbp06cnODIRESmuvDHF6imWsla5cmUgeGZpd6HEWPJZvnw5O7ZD/b06xlS/Tp225OSksmLFigRHJiIiIpI4SoxFRERERFBiLFE0adKESumwZs33MdVft+5nUlNzaNy4cYIjExEREUkcJcaSzz777MOBB+zNgp/fjan+T/NH07r1XnTv3j3BkYmIiBTs7rvvxswwM1JSUmjSpAnnnXceS5Ys2VlnzJgxHHPMMXTr1o0ePXqwcOHCYl1jzpw59OrViypVqtC0aVPuvfde3L3Qc9atW8cFF1xArVq1qFWrFhdccAHr16/fefz333/nuOOOo0mTJmRkZNC8eXOuvfZaNm7cWKzYpPSUGEs+ZsZZZ53Gxo3TWLLki0Lrrl71LatXTeCMM07RjBQiIpJ0++23HytWrGDp0qW8+eabzJkzhzPPPHPn8WOPPZaJEycyY8YM9ttvP8aPHx9z2xs3bqRPnz40adKEr7/+mieeeIJHHnmERx99tNDzzj33XGbPns348eMZP348s2fP5oILLth5PCUlhX79+jFmzBjmz5/Piy++yMSJE7nqqquKfwOkVDRdm0R1/PHHM23aV7wz4m9s2XwVbdqeELHARzYLF07k+zmPc/TRB3D22WcnMVoREZFAWloajRo1AoKhgVdccQXXX389GzdupGbNmjtXah03bhxLlizhySefjLntV199laysLF588UUyMjLYf//9mT9/Po8++ihDhgzBzPKd8+OPPzJ+/HimTp3KIYccAsBzzz1Hjx49mDdvHu3ataNOnToMGjRo5zktW7bk6quv5uGHHy7NrZASUGIsUZkZd9xxOzVqDOeNN5/ip/kvsVeDnlSpUpesbRtYveozKqVt4LTTgiWh09L0R0lERMqXlStXMmLECFJTU3d+q5mbm8v999/PwoULGTVq1M4pxSD4v++FF17g4osvjtrelClT6NWr187kGoIe6KFDh7Jo0SJat24d9ZxatWrtTIoBDj30UGrVqsXkyZNp165dvnOWL1/OiBEj6NWrV0nfupSQhlJIgdLS0vjLX/7CmNGvcs01J9Gwwffg46lf9xsuvaQ377zzAvfcc/cu/0CIiIgk05w5c6hevTpVq1alcePGTJo0iWuuuYZq1aoB8M9//pP77ruP7777jiOPPJInnnhi57nt2rWjVq1aBba9cuVKGjZsuEtZ3v7KlSsLPKdBgwb5yhs0aJDvnHPOOYeqVavStGlTatasyb///e/Y3rTEjbr5pEhNmjRh4MCBDBw4MNmhiIiIFKpdu3aMGTOGbdu2MXr0aN566y2GDRu28/jgwYMZPHhw1HPnzp1bZPuRwyXyHryLNoyioHPyzossf+yxx7jrrruYN28et956K0OGDOHpp58uMiaJHyXGIiIissdIT0/fufzwfvvtx08//cSgQYN45ZVXSt12o0aN8vXyrl69GiBfT3L4OatWrcpXvmbNmnznNGrUiEaNGtG+fXvq1avHEUccwR133KHpUMuQhlKIiIjIHuuOO+7g9ddfZ+bMmaVuq0ePHnz22Wds3759Z9mHH35IkyZNaNWqVYHnbNiwga+++mpn2bRp09iwYQOHHXZYgdfK64netm1bqeOW2CkxjpGZDTWzr81sk5mtNrNRZpZ/xHz+83qZ2QwzyzKzX8xMc6+IiMhuZd68eTz88CMMuvJarrx0EHfcfidffPEFOTk5yQ6tSHvvvTf9+vXjzjvvLLJu+/btGTlyZIHHzz33XDIyMrj44ov57rvvGDlyJH//+993mZHiq6++on379ixbtgyADh06cNxxx3HFFVcwdepUpk6dyhVXXMFJJ52088G7sWPH8sILL/Ddd9+xaNEixo4dy6BBgzj88MMLTLglMTSUIna9gKeArwnu2zDgQzPr6O6bo51gZq2BscBzwPnA4cDTZrbG3d+J9cJm1iyGavlH9ouIiJTC6tWrueuOe5nxxfdUymxAk9TOpFgav+T8wtj/3UHL9g24/e6b6dq1a7JDLdSNN97I4YcfzrRp03aZHSLSvHnz2LBhQ4HHa9WqxYQJE7jmmmvo3r07derUYciQIQwZMmRnnS1btjBv3jx27Nixs+zVV1/l+uuvp2/fvgCccsopu0wTV6VKFZ577jkGDx7Mtm3baN68Oaeeeiq33HJLad62lIAVtVqLRGdmewGrgV7u/lkBdR4ETnH3DmFlzwCd3L1HMa4V84c0Y8aMcv8PlIiIlJ25c+dy/vnn89///pf27dvHfN7q1au58tJr+G1OGr3rXUebmoeSYn980bxyy3w+XfV/bKz/DY8+/XcOPvjgRIQvu7HC/uzNnDmTbt26AXRz99KPc4kTDaUoubz5XNYWUqcH8GFE2QdAdzOrlJCoRERE4uC+e/7OmjmpnNP8SfapddguSTFAo6r7cnrLB6m79hBuu+keNm3alKRIReJHiXEJWDCQ6FHgC3f/rpCqjYDIR1FXEQzFqF+MSzaPYTu+GO2JiIgU6Oeff2bqp9/Qq+7VVK9Ur8B6qSlpHNfkL6xdvL1YSyuLlFcaY1wyTwIHAn+KoW7kMAgroLzgBtyXFlXHzDTGWERE4uK9994jbVM99ml9eJF1q1WqQyvrxci33uWMM84og+hEEkc9xsVkZk8ApwBHxpCwriToNQ7XAMgGfk9AeCIiIqW26JfFNEo5kBRLjal+i2qd+XXBUvTckuzulBjHyAJPAqcCR7n7whhOmwL0iSjrC0x39x1R6ouIiCRdTnYOVowUwSyFXHdyc3MTGJVI4ikxjt1TBFOunQtsMrNGoa1KXgUzu9/MXg475xmgpZk9amYdzOxS4DLgkTKNXEREpBgaN2vIb7k/x9wDvHrrzzRoWI/U1Nh6mPcUU6ZM4aijjqJatWrUrl2b3r17s3Xr1gLrf/bZZ5x88sk0adIEM2PUqFG7HN+xYwc333wzBxxwANWqVaNJkyZceOGFLF++PC7xvvPOO3Ts2JGMjAw6duyYb87mouIL17t377jEVN4oMY7dIIKZKCYBK8K2s8LqNAZa5O2EepVPAHoDs4E7gOuLM4exiIhIWTv22GPZXOVXlm6eU2Td7TlZzNvxASefdlwZRFZ+TJkyheOOO46+ffvy1Vdf8fXXX3PttdeSklJwarV582Y6deq0yxzG4bZs2cLMmTO54447mDlzJiNGjGD+/PmccsopcYn3rLPO4oILLuCbb77hggsu4Mwzz2TatGkxx/fee+/lW0HwjTfeYP78+aWOr7zQw3cxcneLoc7FUco+BTSxsIiI7Da6du3Kvgc257OvnuWsqo+SlpJeYN0vVr5IpbpbOemkk8owwpJr1aoVl112GT/++CNjxoyhZs2aDB06lOuuu65Y7QwePJjrr79+l0U49tlnn0LPOf744zn++IInkcpbQCTcE088wcEHH8zixYtp0SLoe9uwYQM33XQTo0aNIisri+7du/PYY4/RqVOnAtsePnw4ffr0YejQoQAMHTqUTz/9lOHDh/P666/HFN/ee+/N4MGDOfTQQ1m/fj1nnnkm9erV49hjjy30fe9O1GMsIiIiuzAzbr/7FrY3/om3f72FjdtX56uzLWcLHy17mh/S3mTw0EE0ahT5rHn59fDDD3PggQcyc+ZMhg4dyuDBg3dJSC+++OJChwqsXr2aadOm0aBBAw477DAaNmxIr169+OKLL+Ie64YNGzAzateuDYC7c+KJJ7Jy5UrGjh27c2Gvo48+mrVrC15aYcqUKTtX3stz7LHHMnny5Jhj6dixIx988AE//fQT33zzDX369OFf//oXderUKdF7K4/UYywiIiL5dOzYkcf/72FuufEOXlhwNs28By2qdSXV0li9dQE/Z08kvf42bhl6Laeddlqywy2Www8/fGdP77777suXX37JY489Rp8+wfPyjRs3LvRBwl9++QWAu+++m0ceeYTOnTvz8ssvc/TRR/Pdd98V2XMcq6ysLG655RbOPfdcatasCcAnn3zCnDlzWL16NRkZGQA88sgjjBo1irfffpsrr7wyalsrV66kYcOGu5Q1bNiQlStXxhzPvHnzuOGGGzjooIPo1KkTEyZMYNasWQwbNmyPSY6VGIuIiEhUBx54IO+MeYOPPvqIkW+P4bt5z5KTk0OD1vW5esAZnHTSSey1117JDrPYevTokW9/+PDhO/fvv//+Qs/PS5oHDhzIJZdcAkCXLl346KOPeP7554s8PxY7duzg7LPPJjc3l6effnpn+YwZM8jMzKRevV0XXtm6dSsLFixg8eLFdOzYcWf5rbfeyq233goE3wSEc/d8ZYWZP38+w4YNo2vXrnz22Wf873//4/XXX2fNmjVKjEVERGTPV6VKFU466aSdY4iLm0ztLorznho3bgywSwIK0KFDBxYvXlzqWHbs2MGZZ57JwoUL+fjjj3f2FkOQlDdu3JhJkyblO6927drUrl2b2bNn7yyrW7cuAI0aNcrXO7x69ep8vciFOfnkk/OVnXPOOTGfvztQYiwiIiIx2xOS4qlTp+bbb9++fcznt2rViiZNmjBv3rxdyufPn1/ow2uxyEuKf/rpJz755JN8PcNdu3Zl5cqVpKWl0apVq6httG3bNl9Zjx49mDBhAoMHD95Z9uGHH3LYYYeVKM5oifmeQImxiIiIVChffvklDz30EP3792fChAm89dZbvP/++zuPDx06lGXLlvHyyy9HPd/MuOmmm7jrrrvo1KkTnTt35qWXXmLu3Lm8/fbbO+sdffTRDBgwgGuvvRaAzMxMfv75553HFy5cyOzZs6lbty4tWrQgOzub008/nZkzZ/Lee++Rk5Ozs5e3bt26pKenc8wxx9CjRw/69+/Pgw8+SLt27Vi+fDljx46lf//+dO/ePWrMf/7zn+nZsycPPvgg/fr1Y/To0UycOHGXBwaLiq8iUGIsIiIiFcqNN97IjBkzuOeee6hRowb/+Mc/dplybMWKFUUOibjhhhvIyspi8ODBrF27dufDaG3atNlZZ8GCBfz2228796dPn86RRx65c3/IkCEAXHTRRbz44ossXbqUMWPGANC5c+ddrvfJJ5/Qu3dvzIyxY8dy2223cemll7JmzRoaNWpEz549Cx0Wcdhhh/HGG29w++23c8cdd9CmTRvefPNNDjnkkJjjqwiUGIuIiEiFUrNmTd58880Cj8eaBN5yyy27zGMcadGiRbvs9+7du9DVBFu1ahXTaoM1atTg8ccf5/HHH48pzjynn346p59+eoHHi4qvItA8xiIiIiIiKDEWEREREQE0lEJEREQqkMjhDSLh1GMsIiIiIoISYxERkT1aSkrwX31OTk6SI5GKZseOHcAffwZ3B7tPpCIiIlJstWvXBmDp0qXJDUQqnFmzZgF/rBS4O9AYYxERkT1Y/fr16dq1K0899RQNGzakcuXKyQ5J9nA7duxg1qxZPPHEEwwYMIAaNWokO6SYKTEWERHZg6WkpHDXXXdxzjnncPnllyc7HKlABgwYwNChQ5MdRrEoMRYREdnDNW3alIkTJ7J48WKys7OTHY7s4VJSUmjcuPFu1VOcR4mxiIhIBZCenk7btm2THYZIuaaH70REREREUGIsIiIiIgIoMRYRERERAZQYi4iIiIgASoxFRERERAAlxiIiIiIigBJjERERERFAibGIiIiICKDEWEREREQEUGIsIiIiIgIoMY6ZmfU0s3fNbLmZuZn1L6J+71C9yK19GYUsIiIiIsWQluwAdiPVgG+AF4B3inFeO2Bj2P6a4l7YzJrFUK1BcdsVERERkT8oMY6Ru48DxgGYWXFOXe3u60t5+SWlPF9EREREiqChFIk3y8xWmNlHZnZksoMRERERkejUY5w4K4ArgRlABnAB8JGZ9Xb3z4rZVvMY6uxPqEdbRERERIpPiXGCuPs8YF5Y0RQzaw78BShWYuzuS4uqY2YaYywiIiJSChpKUbamAvskOwgRERERyU+JcdnqQjDEQkRERETKGQ2liJGZVQfahhW1NrPOwFp3X2xm9wNN3f3CUP0bgEXA90A6cD5wWmgTERERkXJGiXHsugOfhO0/Gvr5EnAx0BhoEXY8HXgEaApsJUiQT3T3sQmPVERERESKTYlxjNx9ElDgBMbufnHE/kPAQ4mNSkRERETiRWOMRURERERQYiwiIiIiAigxFhEREREBlBiLiIiIiABKjEVEREREgAQnxmZW28wuN7P7zaxuqKyrmTUtYXsZ8Y1QRERERCSQsOnazOxAYCKwAWgFPAesBQYALYELY2jjWOAc4AiCOYJTzGwLMBP4EHjB3ZcnIn4RERERqVgS2WP8KPCiu+8DZIWVjwN6FnaimfU3s3kEi2fkAg8DpwLHApcBnwLHAL+Y2TNmtlcC4hcRERGRCiSRC3wcBAyMUr4MaFTEubcCfwHed/fcKMf/BxAakvFngt7nf5Q8VBERERGp6BKZGGcBNaOUtwPWFHaiux8cywXcfRnw1+KHJiIiIiKyq0QOpRgN3GlmlUL7bmYtgAeAd4rbmJmlm1k7M9My1iIiIiISd4lMjP8C7AWsBqoQjAv+GdgE3BZrI2ZW1cz+A2wBvid4CA8ze9zMbol30CIiIiJSMSWs99XdNwJ/MrOjgK4ESfhMd59YzKbuBzoBvYHxYeUTgXsIeqBFREREREolkdO1XQi86e4fAx+HlacDZ7v7yzE21R84y92nmpmHlf8AtIlXvCIiIiJSsSVyKMULQK0o5TVCx2KVNxwjUjXAo5SLiIiIiBRbIhNjI3ri2oxg0Y9YfQ2cGLaf1+YVwJSShSYiIiIisqu4D6Uws1kEyasDH5lZdtjhVKA1u44VLspQYLyZdSSI989mth/QA+gVn6hFREREpKJLxBjjUaGfnYEPgMywY9uBRRRjujZ3n2xmhxPMcrEA6EuwJHQPd59T+nBFRERERBKQGLv7PQBmtojg4busws+Iqc05wEWlbUdEREREpCAJG2Ps7i/FIyk2sxwzaxClvJ6Z5ZS2fRERERERSOx0banAYOBMgkU50sOPu3vdWJsqoDyDYGiGiIiIiEipJXJ55buAy4FHgb8Bw4BWBPMS31vUyWZ2feilA5ebWfhY5VSgJzA3fuGKiIiISEWWyMT4POAKd3/fzO4CXnf3BWb2LXAo8HgR5w8O/TTgKiB82ETeQ3xXxTdkEREREamoEpkYNwLyZo3I5I/FPt4j6EEulLu3BjCzT4BT3X1dIoIUEREREYHELvCxFGgcev0zwTRrAAcB22JtxN2PVFIsIiIiIomWyB7jkcDRwDTgn8DrZnYZwYN4jxV2opk9Ctzh7ptDrwvk7kPiFK+IiIiIVGAJS4zd/Zaw12+b2RLgcOBndx9TxOldgEphrwu8TOmiFBEREREJJLLHeBfuPo2g9xgzO8jdvy6k7pHRXouIiIiIJErCxhibWXUzqxJR1tnM3gWmlqLdmmbW38zalzrI4l23p5m9a2bLzczNrH8M5/QysxlmlmVmv5iZZtEQERERKafinhibWTMz+xLYAGwws0fNrKqZvQx8TfDg3Z+K0d7/zOza0OsqwHTgf8AcMzst3vEXohrwDXBtLJXNrDUwFvicYDjI34HHyzhmEREREYlRIoZSPABUB/4MnBb62YsgqdzX3RcWs72eBIuDAAwgmNe4NnARcDvwTulDLpq7jwPGAZgVtBjfLq4CFrv7DaH9H82sO/AXihmzmTWLoVq+ZbNFREREJHaJSIyPBM509y/N7G1gOfCWuz9QwvZqAWtDr48D3nH3LWb2PvBw6cNNmB7AhxFlHwCXmVkld99RjLaWxC8sEREREYkmEWOMGwELANx9JbAVGF2K9pYAPcysGkFinJds1gGyStFuojUCVkWUrSL4ZaR+2YcjIiIiIoVJ1KwU4cs351K6BHY48CrB6nm/ApNC5T35Y2W98ipyOjkroLwozWOosz+hoR4iIiIiUnyJSIwN+MjMskP7VYB3zWx7eCV37xpLY+7+tJl9RZAcTnD33NChXwjGGJdXKwl6jcM1ALKB34vTkLsvLaqOmWmMsYiIiEgpJCIxvidiv0TDKMxsCjAKGO3u0wlmo9jJ3d8vUXRlZwpwckRZX2B6MccXi4iIiEgZiHti7O6RiXFJPQOcAtxuZisIEuwxwBfuXuYr3plZdaBtWFFrM+sMrHX3xWZ2P9DU3S8MHX8GuDa0pPVzBA/jXQacU4Zhi4iIiEiMErbAR2m5+0vufhrBg2o3ADWBN4HVZvaimQ0ws6plGFJ3YFZoA3g09Pre0H5joEVe5dC0dCcAvYHZwB3A9e5eJtPLiYiIiEjxlNmS0CXl7tsIFsoYCww0s0MIepLvBV41s4+B+939ywTHMYk/Hp6LdvziKGWfAjGNpRYRERGR5Cq3PcYFcfdp7n6bux8AHAB8RNBbKyIiIiJSYuW+x9jMmgOeNzODmR0MnAv84O7PAo8lMz4RERER2TPsDj3GrxGspoeZNQImAgcDfzezO5MZmIiIiIjsORLaY2xmRwNHE8zfu0sS7u6XxtjM/sBXoddnAnPc/XAz60sw88O9BZ4pIiIiIhKjhCXGZnYXcCfB/MMrKP5qb3kqAdtCr48hmLINYC4aWywiIiIicZLIHuOrgIvd/ZVStvM9cJWZvQ/0IZj2DKAJxVxBTkRERESkIIkcY5wOTI5DOzcDA4FJwOvu/k2o/BT+GGIhIiIiIlIqiewx/jfB7BF/K00j7j7JzOoDNd19XdihZ4EtpWlbRERERCRPIhPjysCVZnYM8C2wI/yguw+JtSF3zwHWRZQtikOMIiIiIuXWmjVreO+99/h80hQ2rs+kWrUqdDu0E/369aNly5bJDm+Pk8jE+ECCpZAhmFkiXMwP4plZQ+AR/pjdYpfV59w9teQhioiIiJQ/ubm5PPvss7zy7FtsX5tBq5QjqFlpL7JyNvH65Em8+txIju3fi1tvu4XKlSsnO9w9RsISY3c/Mk5NvQi0IBiSUZrZLURERETKPXfnkUf+wZv/N47uGZfRrcUAMlKr7jyek3sdP67/mHGvDGfd7zfzj+EPk56ensSI9xxlsvKdmTUjWL1uWQlO/xNwhLvPjm9UIiIiIuXPp59+yv+eH8tR1Yeyf92++Y6npqSxf92+1MloyojxQ3j11Ve55JJLkhDpnidhs1KYWYqZ3WlmG4BfgcVmtt7M7jCz4lx3CRHDJ0RESmLVqlV8++23zJkzh99/12yPIlI+vfO/key1rXPUpDhc02r70T71RN55fQzZ2dllFN2eLZE9xsOAy4BbgC8JktvDgbsJHsy7LcZ2bgAeMLOBeuBORIrL3fniiy8YMeIdZk6fgns2BlhqOocd3pvTTjudbt26JTtMEREAli5dyleff8uRte6KqX6X+v14beFIJk+eTM+ePRMc3Z4vkYnxRcDl7j4mrOwbM1sGPE3sifGbQFVggZltIf/sFnXjEayI7HlycnJ48MEH+WDsO7RrU5s/X3YobVo3xN35Yf5yxn30JUNumMh5F1zBFVdcgZm+nBKR5Fq2bBnZW6Fpw/1iql+/cktSs6uxfPnyBEdWMSQyMa5LsGxzpLmhY7G6IS7RiEiF8/jjjzNh3NvccEVPeh7WfpdjzZvWo2/v/RkzfiYvvvwM1atX59xzz01SpCIiAfdgjgErxihSw3aeJ6WTyMT4G+Ba4PqI8mtDx2Li7i/FMygRqRiWLFnCqBFvcOlZB+VLivOYGf2O78ba9Zt54fn/4+STT6ZGjRplHKmIyB8aN25MWmVYvuVH2qUfUWT9tVlL2ZGaSePGjcsguj1fIhPjvwLvhxb4mEIwzdphQHPghOI0ZGapQH+gQ6idH4AxoYU/JAHWr1/PuHHj+GTih6z9bTXp6Rns02E/+vUfQJcuXfSVs5R7o0ePpnrVXI496sAi6556Ynfen/gj48eP54wzziiD6EREomvZsiVdDuvINxPG0K520YnxrN/H0LBVTQ477LAyiG7Pl7BZKdz9U2BfYCRQm2D4xAignbt/Hms7ZtYW+BF4GTgVOB34L/C9mbWJc9gVnrvz9ttvc2b/k3hh+DBaZC3kxJYp9Nori1+njOXGQZcy8PJLWb16dbJDFSnUZ59+RM9D2lCpUtFrANWqWZVuBzbhs88mJT4wEZEinHZGf1ZVms6P6yYVWm/llp/4YccYBpx5ouYxjpOEzmPs7suJ/SG7gjwOLAAOdfe1AGZWjyA5fhw4sZTtS5jXXnuNfz/+MKd3bsQFvXtQs2rGzmNX9nW+WbiKB8fM4LpBA3n62X9Tr169JEYrUrBNGzdQv17svzvvVbcGPyxaV3RFEZFScncyMzNZu3Yt69atY+3atfleUzOT//vxTKpbI4wUtmSvZ3P2Wq7bbwQd6hzFTxu+ZOK6h+h6TGsuvPDCZL+lPUZcE2MzOxD4zt1zQ68L5O7fxthsL8KS4tC5v5tZ3jRwEicLFizg308O58JDmnPJ0Z3yHTczOu/diCcuPZJB//6YRx95hGH335+ESEWKlpFRmS1bt8Vcf0vWdipX1iQ3IhK77du3s27duqjJbeTPyLJo8w6bGbVr16ZOnTrUrl2bBs3q8NuqlaTlVKVhlX2oWWkvft44lcnr/8326qvodfrB3H3PnVoSOo7i3WM8G2gErA69dqIvzuFA0d9vBrYB0Z6GqQ5sL3aEUqBRI0dSP30HF/Y+oNB69WtW5aIj2vHoJx+xcuVKGjVqVEYRisTugE5dmTL9U845tUeRY+K3b89m+jdLOf7ko8ooOhEpL9ydjRs3FprcFlS2efPmqG1WqVKFOnXqUKdOHerVq0edOnXo0KFDvrK6devuUlazZk1SU3dNj5YtW8aYMWP4dOJkNm7IxKp+z0mHdaF//360bx/9wWIpuXgnxq2BNWGv4+E94Fkzuwz4KlR2CPAMMKbAs6RYtm/fzoSxYzi3SwtSU4seen5Mp9b866Mf+OCDD7jooovKIEKR4hkw4FQGfzKOb39YQqf9WhRa9/Op88jcAv369Suj6EQk3rKysvIlsAX13oYnt+vXrycnJ/+z/CkpKbskr3Xr1qVx48bst99++cojE9549uA2bdqUQYMGMWjQoLi1KQWLa2Ls7r+G7bYEJrv7Lt8VmFkawewU4XULcz3wEsHMFnmLe6QRJMV/LlXAstO6devYtmUzHZu3jal+5fQ09q6vCcWl/OrSpQsdD+jOP5/9iPuG9qdJozpR6/20YCX/eW0KvY86gebNm5dxlCISLicnhw0bNsQ8HCG8bOvWrVHbrFat2i5JbN26dWnatGm+sshkt0aNGqSkJGyOAimnEvnw3SdAY4JhFeFqhY7FNJTC3dcD/cxsH6A9wdCMH9z95/iFKiWjKduk/DIzhg27nz9ffy1/vXcEJ/fZjz6996duneoArFy9gfEff8u4j+bSbr+DueWWoUmOWGTPsXXr1mIPS1i3bh3r16+PulBFampqvkS2efPmdOrUqdDktk6dOpqtQYolkYmxEYwljlQPiD4opxDu/hPwU2mDkujq1KlDRtVqfL/4N7q3bVJk/a3bdvDLb5l0b1J0XZFkqVu3Lk//6/949tlnGT3+Xf737jfUrVUFx1m7Povqtfbi9LOv4JJLLiEjI6PoBkUqkJycHNavXx9zUhv+c9u26A++1qhRI18i27Jly0KT27p161K9enXNny9lIu6JsZmNCL104EUzC//bkQocCEwuoo1HgTvcfXPodYHcfUhp4pVAeno6fU7sx7ujX+G8XvuTVsQ444++XcTWlCocd9xxZRShSMnUqFGDG2+8kYEDBzJp0iRWrVpFSkoKTZs2pVevXkqIZY/m7mzZsqXIhDba8Q0bNkRts1KlSvkS2datW9O1a9dCk9vatWtTqVKlMr4DIsWTiB7jvL9JBmwCwgf9bAemAs8V0UYXoFLYaykDAwYM4L133uClj7/lsj6dC6y3ev1mXvxsLocfeRINGzYsuwBFSqF69eqcdNJJyQ5DpER27NgR9WGyWB422749+gROtWrVypfE7r333oUmt3Xq1KFatWrqvZU9VtwTY3e/BMDMFgGPuHtJhk0cGe11eWBmVwM3EYyf/h64oaCV/MysN8F46kgd3H1uomIsqb333puB19/I/z32IFu27eCC3gdQu/ofT9a6OzMXrOShd2eQ1qANNwxWZ72ISKzcnU2bNhVrzG3e68zMzKhtZmRk5Etk99lnn12S2WiJbu3atfNNCyYiCRxj7O73xKMdM3se+LO7b4oorwY84e6XxuM6McZyFjAcuJpgcZGBwDgz6+juiws5tR2wMWx/TUEVk+3ss88mPT2dZx5/lDHfjOdPberRvH4Ntu3IYcrPq1myKZd2B3bnb3+/X6veiUiFlLeoQ0keLos2LVjeog7hCWzDhg1p3759kb23VapUScIdENlzWbSnP+PSsNlCoj98B4C77x1jOzlAY3dfHVFeH1jp7gld1jrimtOAme4+KKzsR2CUu+d7pD2sx7hOaHaNkl63WQzV9gfGzZgxg65du5b0Ujtt3LiR8ePH8/GED1j3+xrSMzLYp/1+9Os/gAMPPFBfo4nIbi03N5dNmzaVKLktbFGHohLZaGW1atXStGBS4cycOZNu3boBdHP3mcmOJ08ik8rhEfuVCMYLHwc8XNTJZlaTYJyyATXMLCvscCpwAvmngksYM0sHugEPRBz6kGBe5sLMMrPKwA/Afe4ebXhFYZYUs36p1axZkzPPPJMzzzyT3Nxcfvnll51LWG7ZsoVq1aqVdUgiIvlkZWXFNOdttLG4ubm5+dorbFGHwpLbeC/qICLJkcihFP+MVm5m1wDdY2hiPUGPswPzo10CuKuk8ZVAfYKEfFVE+SqCZbCjWQFcCcwAMoALgI/MrLe7f5aoQOMlMzOT0aNH8+Jz/8fcufPYmON4Sgqpqakc0K4d1109iAEDBpCWVmad9iKyBypqUYdoCW+sizqEJ7HNmjUrsjdXizqIVGzJyGjGAfcDlxRR70iC3uKPgdOAtWHHtgO/unsyll2LHB5S0HzNuPs8YF5Y0RQzaw78BShOYhzLclz7E9zbuFixYgU3XHsNkz//jO1VqlG1U3f2PegQUqvXYMOGDXw3YzqX/fUWXnz1VV587jn22muveF1aRHZD7r5zUYfiPly2YcOGqIs6pKWl5Uti8xZ1KKr3Vos6iEhJJCMxPp1dk9yo3P1TADNrDSz2RA2Gjt1vQA75e4cbkL8XuTBTgfOLc2F3X1pUHTNrUJw2C5OZmcmN113L9C8+I63l3rS7ZCC12+yz83gzoEOvI5n74zw+f+0FLrr8ct545RVq1qwZrxBEJEmys7NZv359sZPbdevWFbioQ82aNaPOexue2NapU4d69ertUqZFHUSkrCUsMTazWezak2oESeVeBLM6xKol0LKgfxzLakiCu283sxlAH2Bk2KE+wOhiNNWFYIhFuTV69Gi+mzYZ36sh+w66gRpN8z/7l2JGhw7tyLn4Kr56/ileeeUVrrnmmiREKyKR3J3NmzeXKLktbFGHyKm/2rRpU2Ryq0UdRGR3ksge49HsmhjnEkxTNqmYc/hOilIW3m5ZTsT4KPCKmU0HphCMH24BPANgZvcDTd39wtD+DcAigvmO0wl6ik8LbeVSbm4uo//3JrnuVDvkT1GT4jxmRtOmTVjb5RBeGzmKyy+/XKuIicRR5KIOBS3iEC3h3bFjR9Q2wxd1yEtk27ZtW+DMCXn1tKiDiFQEiXz47u44NVUnYj9vdou/AbfF6Roxcfc3zawecCfBAh/fASe4+6+hKo0JEuU86cAjQFOCFQC/B05097FlF3XxzJ8/n4XzfiSrclX2PuTwIuvXqZxOjc4HsWj6l0ydOpVevXqVQZQiuw93JzMzs0TTgm3atClqm7Eu6hCZ5NaqVUsPy4qIFCKRQykKmn+4HrDa3WPq6XX3aN/rTTCzbcBjBFOolRl3fxp4uoBjF0fsPwQ8VAZhxc2mTZvYsX0buelVqNqwoMk2/mBmZNSry5bUNH7//fcyiFAkObZt21bobAkFJbd50xxGMrOdvbF5PbcNGjSgffv2BQ5L0KIOIiKJlciug4K+c8sgmFWitNYQrCgncZSRkYGlpODueE7+/8yjyclxyM1VT5SUe7m5uWzcuDGmeW4jj23ZsiVqm1WrVs2XuHbs2DFqchter2bNmpoWTESknIl7JmNm14deOnC5mYUv8J4K9ARiHmNsZgdGFhEMWbgF+KYUoUoUbdu2pXb9Bqxdtpx18+fRoEvhHfKbtu9g66+/UD0tldatW5dRlFLRZWVllWhasPXr1xe4qEP4eNq6devStGlTDjjggF3Koj1kpnH1IiJ7jkR08Q0O/TTgKoIpzvJsJ3gY7apitDebIMmO7IGeClxaogilQFWrVuWUs87hqQeGsfKLSezVuWuhD9ys2ZLF9m9mcETHDuy///5R6+Tm5vLVV18x6p13mDdrFtu3baNWvXr0PuEETj75ZBo3bpyotyPlWEGLOsSS8Ba0qEP16tXzJbLNmzfPVxaZBGtRBxERgQQkxu7eGsDMPgFOdfd1pWwyshsyF1jj7lnRKkvpnX766bz50gssnD2DJZ9+RIvex0Stt2rzVlbOnEHln3/g/PvujZpAr1y5ktv+8hd++fpr9tm2jQFVqlA5JYXlK1YwZvZs3vjXvzj3mmu47LLL9MT7bqiwRR2KSm7Xr18f86IOLVu2pHPnzoUmt1rUQURESiuRs1IcGaemegFvuvsuM8ebWTpwtru/HKfrSEjz5s0Z/syzXHbBefz83xdYt3IVbY8+hmp7NQQgc/sOlq9cxYopX1Jp+mSuu/A8jj/++Hzt/Pbbb/z5yitJ/eYbnmzYkI5Vq+6S/F6Vk8Pba9bwnwcfZPv27Vx9dXGmt5Z4ylvUoTiJbSyLOkSOvW3VqlWRMydoUQcREUkWi+eCcmb2KHCHu28OvS6Quw+Jsc24zG6xpzOzrsCMGTNm0LVr17i0+f333zP05pv54uuv2ZFRhfR9O2A1a5GdmQmLFtCybm1uvHoQF1xwQdRE5s7bb+fHV17hX82aUb+QnryRa9bwODD81Vfp1KlTXGKviPIWdYjlQbLIso0bN0ZtMz09vcAZEgrqva1Xrx61a9fWw5glsHXrVj766CPGfTiBVb+tJaNSJfZr15Z+/frRsWNH/cIgInuMmTNn0q1bN4Bu7j4z2fHkiff/XF0I5hkG6MquC3GUlBXQTjMg+hJNEhf77bcfY957j59//pnnn3+eOXPmkLNtEy32bs7x11xB3759C5w2as2aNXzx/vvcULNmoUkxQP/69Rm1aBGjRoxQYkzBizrEkvBGmxYM/ljUIW+LXNQhWrJbp04dqkb08kviTJgwgWEPP8aydVtJa9mdyvW6kZuTzVcTp/Pa6A85rEtHhv3tHvbaa69khyoisseKa2IcPnzC3XuXpq2wJaUd+MjMwv/HTyUYezy+NNeQ2LRt25a///3vxTrnww8/pPL69RzTvHmRdc2MftWr8/T775N5001Ur169pKGWG+7Opk2bSjRzQmZmZtQ2wxd1yEtc991330KT27xFHVJT9cVKefb+++9z630Pk9v2KNqddyGVazfceczdWffTdD5+/3GuGHQt//m/p6lXr14SoxUR2XMlcoGP54E/u/umiPJqwBPuXtSMEqNCPzsDHwDh2ULe7BbvxCNWib9Vq1bRDKgaY0K2b5Uq5GzZwu+//16uEuO8RR1KMnNCTk5OvvbMjNq1a++SyDZs2JD27dsXOOZWizrs2ZYvX87dD/wD63gi+558Tb4eejOj7r4HUbXBP/ju+cE8+NDDPPTgA0mKVkRkz5bIQYAXEcw1HLmmaRXgQoqYas3d7wEws0UED99pFordSEpKChTjK/i8mWUT8bV95KIOxUluN2/eHLXNghZ1KCq5rVWrlqYFk12MGTOGdblVOeD4Kwr981+5dgMa9b6ICROGs3z5cpo0aVKGUYqIVAyJWOCjJsG4YANqmFl4QpsKnACsjnZuNO7+UnwjlLLQtGlTxrqzMTubmjE8hDVn82bSa9UqdPxk3qIOxU1uC1rUITU1NV8Cm7eoQ2HJrRZ1kHjJzc3lrdHvU+OAPqRWKvrP1F4H9GLVxOcYO3Ysl19+eRlEKCJSsSSix3g9f4wNnh/luAN3FdaAma0F9nX338xsHYU8xOfudUseqiRK3759ee7hhxm3di1nNWiwszzHnQ3Z2azdsYN12dmszc7mtx07eHTlSuo3a8btt99eYMKblRX9S4Pq1avnS2RbtGhRaHKbt6iDHiyTZNq4cSO/r99E7Z7RF8eJlFopg9QGbVm2bFmCIxMRqZgSkRgfSdBb/DFwGrA27Nh24Fd3X15EG4P5YwjGYOIzu4XESbRFHaIlsquysrj2l18Y9uuvbM7JYW12Nhuys6N+mAbUnTmTXxcv3pm8tmzZki5duhSa3NauXVuLOsjurzi/oJlF/QZERERKLxEr330KYGatgcVegomSw4dPuPuL8YtOCvP999/zzTffxDQt2Pbt26O2UbNmzT8S2bp12bjXXmxct47DatTgwOrVqVepEnUrVaJuWhoOfLJuHZMqV+aae+7h7LPPLts3LJJkNWrUoHaNqmQum0/9Dj2KrJ+bk03Oml9o1Oi4MohORKTiSeTDdx2A5sAXAGZ2DXAF8ANwTaxLRYcWrtjh7nNC+/2AS0Lt3O3u0TM0KbY33niD++67b+eiDuG9s23atInaYxteVqtWLSpVqrRLm5s2beJvd93F1x9+yPbMTA7NzmZbbi4zNm9mmhnpzZvzl5tuol+/fkl61yLJk5qaSv8T+vLUiPHkHnkeKamF/5P82w9fUjV7Q9SVJkVEpPTiuvLdLg2bzQFudvexZnYAMB34B3AU8KO7XxJjO18DD7j7O2a2N0FCPAI4CHjf3W9IyBvYzcRj5bvMzEzMLCGLOixYsIAxY8Ywd/ZstmVlUbt+fY7s25c+ffpQtWrVuF5LZHeyaNEiBpx7MandzqbVMRcVWG/H5g388J8hHLt/A5785/CyC1BEJAEqysp34VoTJLEQjDV+191vDSVwY4vRzr7A7NDrM4BP3f1cMzsceAO4IT7hSiLnD27Tpg2DBw9OWPsiu6tWrVpx07VXMuyfz7JgexYtjjyXSlVq7FJn0/KfWTjiIfauuoVbb/5rkiIVEdnzJTIx3g7kdQUeA7wcer0WqFmMdgzIm/j1GOC90OslQP1SxigiknTnnHMOlSpV4pEn/4/vvxlLlXY9qVq/ObnZ29n489ekrJlP57bNeeSBxzV/sYhIAiUyMf4CeNTMvgQOBs4Kle8LLC1GO9OB281sItALGBQqbw2silOsIiJJY2acccYZHH300YwdO5b3xk9g9ffTSKuURs8ObTn11mEceuihWtpbRCTBEpkYXws8DZwODHL3vIk3jwfGF6OdG4BXgf7AMHf/OVR+OjA5LpGKiJQDdevW5fzzz+f8889PdigiIhVSwhJjd18MnBSlvFgDTd39W+CAKIduAnJKFp2IiIiIyK4S2WOMmaUCAwimbnNgLjDK3bNL0Fb38HbcfXo8YxURERGRii1hibGZ7Q+MARoC80LF+wJrzOyUvHmJY2inGfA6cDjBctMAtc1sMnCOuy+Ja+AiIiIiUiGlFF2lxP4NfAc0c/eu7t6VYMGPb4Fni9HO80AloIO713X3ugQ9xwb8J84xi4iIiEgFlcihFJ2A7uEr3Ln7OjO7Dfi6GO0cARzm7nm9zrj7PDO7DvgybtGKiIiISIWWyB7jeQTDKCI1AH6OUl6QxQQ9xpHSgGVRykVEREREii2uibGZ1czbgFuBx83sdDNrFtpOB4YDNxej2b8CT5hZdwutUxx6EO+fwF/iGb+IiIiIVFzxHkqxnmDWiDwG/C+szEI/3wUKnKnezNZFtFMNmAZkh3LjNCCbYPzxqNKHLSIiIiIVXbwT4yPj1M4NcWon7szsaoI5lBsD3wM3uPvnhdTvBTwK7AcsBx5y92fKIlYRERERiV1cE2N3/zRO7bwUj3bizczOIhgKcjXBg38DgXFm1jG0oElk/dbAWOA54HyCKeeeNrM17v5OmQUuIiIiIkVK5DzGPQs77u6fxdhOiyLayZeQJtAQ4D/u/u/Q/g1mdiwwCBgapf5VwGJ3vyG0/2NofPRfgJgT49BczkVpEGt7IiIiIpJfIqdrmxSlLHzccIFjjCMsijgvUqztlIqZpQPdgAciDn0IHFbAaT1Cx8N9AFxmZpXcfUeMl0/YIia5ublMnTqVkaNHM2/hL2Rn76Bhvfqc0PdYjjvuOGrUqJGoS4uIiIiUK4mcrq1OxNYAOI5gDuO+xWinC9A1bDuEoCd2PnBGHOMtSn2CJHxVRPkqoFEB5zQqoH5aqL2k+vXXXzn3ogu46ta/8unapeQc2YW0Ew7jl0bVuOeZJzjx9FN5//33kx2miIiISJlIWI+xu2+IUjzBzLYBjxH0vsbSzjdRiqeb2XKCh+BGlDzKEonsvbYoZUXVj1ZemOYx1NkfGBdrg0uXLmXg9dexpmY6PYbfTv12bXY5vnXdBr59+W3uePgBsrOz6devXzHCFREREdn9JHIoRUHWAO3i0M584KA4tBOr34Ac8vcONyB/r3CelQXUzwZ+j/XC7r60qDpmVqwxxsMeeIDV1dM46h+3kVGjer7jVerU4uDrL2VGejr3//MxDj74YBo3blycS4iIiIjsVhL58N2BkUUEU5zdAkTrBS6onZoFtHM38FMpQiwWd99uZjOAPsDIsEN9gNEFnDYFODmirC8wvRjji+Pup59+Ytp337D/HVdHTYrzmBmdLz2LcR9PYcyYMQwcODBqvaysLD777DNWrFgBQJMmTejZsycZGRkJiV9EREQkERLZYzybYLiARZRPBS4tRjvriT4cYQlwdgljK6lHgVfMbDpB0nsl0AJ4BsDM7geauvuFofrPANea2aMEU7b1AC4DzinjuHcxfvx4cuvWpOkhXYqsm5aRTpNjDmfk2Pe58sorCS2wAsCWLVt44YUX+N+o91mxfitUrQ04bNlAs3r/5Ix+J3LJJZdQuXLlxL0ZERERkThJZGLcOmI/F1jj7lnFbCdy0ZBcguEYP7t7dkmDKwl3f9PM6gF3EvRafwec4O6/hqo0JkiU8+ovNLMTCMZUX0OwwMf1yZ7DeNmyZVRr24KU1Ngm9KjXbm/mj/qYrVu3UrVqVQA2bdrE9TcM4YsfllKz28ns2/14KtcJRo1sXbuCFV+/z6Mvj2LGrG/452P/oFq1agl7PyIiIiLxkMiH734tulZM7cRl0ZB4cfengacLOHZxlLJPCWbT2G25Bx32eb3F7s4dd93N53NX0vbif1C90a6/A1Wp25i9j72cTQf0ZtLLN3PP3+7joQfuL/O4RURERIoj7tO1mdkhZnZ8RNmFZrbQzFab2bNmFvPgUzO7yMxODNt/yMzWm9lkM2sZz9grimbNmrH5p1/JzY6tw/23uT9Tv3adnUMi5s6dy0eTZ9DspBvyJcXhajRpS5MTrmP8pCksWLAgLrGLiIiIJEoi5jG+G9j54J2ZHQD8B5hIsDjGyURfJa4gtwJbQ231AK4F/kowS8RjcYm4gjnhhBNIXbeJJVNmFlk3O2sbKyZM5rSTTt7ZYzxy5EiyqjSgXvtDizy//n5/Ymt6HUaNGlXasEVEREQSKhGJcWfgo7D9s4Fp7n6Fuz8KXA+cWYz2mgM/h173B95292cJkusjSh1tBbT33ntzWKeufP+fN8lav7HAeu7OjGdfpcb2XE4++Y/JNb6aNYfq7Q7HUor+45OSmkbVfQ/nq1nfxiV2ERERkURJRGJch13n9e0FjA/b/5rYFqzIkwnUC73uS9DzDJAFVClhjBXebUOH0nRHKh/feB+rvv1x5zjiPJtX/86UR/6PdeO+4I6//JWGDRvuPLYlK4u0yrE/TJdWuSqbt2yNW+wiIiIiiZCIh+9WEcxIscTM0gkePLsr7HgNoDhz+E4A/m1ms4B9gbw1ivcDFpU62gqqUaNGPPvkU9xy+21Mv/khUps3Yq9u+5OSlsrGxStY9/W3NKxSnQduu5M+ffrscm6dmjVYueG3mK+1bcNv1K0dOR21iIiISPmSiMR4PPCAmd1MMPRhC/B52PEDgeI8iXUNcB9BL/Np7p63Ylw34PVSR1uBNWnShJf+8zyzZs1i1OjRzJvxMzuyd9Cy3l4cP/ivHHPMMTunZwt3TO8jmPHCSHK2DyQ1vfA5irOzNrPtpy85etC5iXobIiIiInGRiMT4dmAE8CnBMIiL3H172PFLgQ8La8DMWrj7YgB3X0/wwN0u3P2usPpN3X1Z6UOveMyMrl270rVr7DPKnXzyyfzrxddYPu1dmh9xRqF1l00ZTa20HZx00kmlDVVEREQkoeI+xtjd17j7EQRjjeu4+8iIKmcA9xTRzNdm9pyZHVxQBTOrZWZXmNl3wKmli1qKo1GjRlx27ums+/QFVs2aWGC9FdPHs3Hyfxl40TnUq1evwHoiIiIi5UEiF/jYUED52hhO70AwTdt4M9sBTCdYNS6LIOHuSDDGeDpwk7uPi0vQErNBgwaRmZnJyyP+wZoZY2lw0MnUaLYv7s6mJXNZPf09Kq2Zx1Xn9OPyyy9PdrgiIiIiRbLI2QjKEzOrDJxAMC1bK4JZKH4DZgEfuPt3yYuufDGzrsCMGTNmFGtYRGm4O1988QVvvTOCz6bNYltOUJ6RCkcd1p3TTzuVHj167Jz/WERERARg5syZdOvWDaCbuxe9sEIZSViPcTy4exbBeOURyY5F8jMzjjjiCI444ghWrFjBypUrAWjcuDGNGjVKcnQiIiIixVOuE2PZfTRu3JjGjRsnOwwRERGREkvEAh8iIiIiIrsdJcYiIiIiIigxFhEREREBlBiLiIiIiABKjEVEREREACXGIiIiIiKAEmMREREREUCJsYiIiIgIoMRYRERERARQYiwiIiIiAigxFhEREREBlBiLiIiIiABKjEVEREREACXGIiIiIiKAEmMREREREUCJsYiIiIgIoMQ4ZmZWx8xeMbMNoe0VM6tdxDkvmplHbFPLKGQRERERKYa0ZAewG3kNaAYcF9p/FngFOLmI88YDl4Ttb49/aCIiIiJSWkqMY2BmHQgS4kPdfVqo7Apgipm1c/d5hZy+zd1XlvL6zWKo1qA01xARERGp6JQYx6YHsCEvKQZw96lmtgE4DCgsMe5tZquB9cCnwG3uvrqY119SzPoiIiIiUkwaYxybRkC0ZHZ16FhBxgHnAUcBNwIHAR+bWUbcIxQRERGRUqnQPcZmdjdwVxHVDgr99GhNFFAenOD+Ztjud2Y2HfgVOBEYEXukNI+hzv4EibiIiIiIlECFToyBJ4E3iqizCDgQaBjl2F7Aqlgv5u4rzOxXYJ9Yzwmdt7SoOmamMcYiIiIipVChE2N3/w34rah6ZjYFqGVmB7v7V6GyQ4BawORYr2dm9Qh6f1eULGIRERERSRSNMY6Bu/9IMO3ac2Z2qJkdCjwHvBc+I4WZzTWzAaHX1c3sETPrYWatzKw38C5BIj6yzN+EiIiIiBRKiXHszgPmAB+Gtm+BCyLqtCPoRQbIAQ4ARgPzgZdCP3u4+6ayCFhEREREYlehh1IUh7uvBc4voo6Fvd4KHJvouEREREQkPtRjLCIiIiKCEmMREREREUCJsYiIiIgIoDHGUg6sXr2aSZMm8fvvv1OpUiXatm3Ln/70J9LS9MdTREREyo4yD0maZcuW8fTTT/PZ5El4pRxq1a9Bzo4c1v93E3vVbsSZp53JeeedR0qKvtgQERGRxFNiLEmxYMEC/jz4OrKrbOXkq3tz8NGdyKiSAcCKxav5bMxU/vX8E8yfP5+7776b1NTUJEcsIiIiezolxlLmtmzZwl/+eiOV6udywwMDqVajyi7HG7dowFnXnkLH7vvywr1v0+zfzRg4cGCSohUREZGKQt9RS5n78MMPWbl2GZffeU6+pDjcAYe258gzD+Htkf9j8+bNZRihiIiIVERKjKVMuTsjRr1Dxx5tqNuwdpH1jzjpEDK3beTDDz9MfHAiIiJSoSkxljK1ZcsWflm0gM5/6hhT/dr1a9KiQ2O+//77BEcmIiIiFZ0SYylTWVlZOL7zQbtYpFetxLZt2xIYlYiIiIgSYylj1atXJ8VSWbt6fUz13Z11qzZSs2bNxAYmIiIiFZ4SYylTGRkZHH7I4Uz7YDbuXmT9X+ct5fclG+jZs2cZRCciIiIVmRJjKXMDBpzKyl9+54evfyq0nrsz/rVPadG4FQcddFAZRSciIiIVlRJjKXMHH3wwPQ/pzcsPjGDe7F+i1snJyeH1x0ezYPoyrr36Oq1+JyIiIgmnBT6kzJkZd999D7ffcTvPDH2Ntl2ac9gJ3WnUvD7ZO3L4/uv5TBk7i6y12dz21zs44ogjkh2yiIiIVABKjCUpqlSpwoMPPMjHH3/MyFEjeXXYu+SSg2FUzahG36NP4LTTTqNNmzbJDlVEREQqCCXGkjRpaWn07duXvn37smLFCtauXUtaWhrNmjWjWrVqyQ5PREREKhglxlIuNG7cmMaNGyc7DBEREanA9ESTiIiIiAhKjEVEREREACXGIiIiIiKAEmMREREREUCJsYiIiIgIoFkpJEE2bdrE+PHj+XLKl2zM3ES1KtXo1qUrJ510EvXr1092eCIiIiL5KDGWuHJ3Xn75ZV547WU25W6l6SHtqNasBus2bWHqm//huVee5/STB3DdtdeRlqY/fiIiIlJ+KDORuHF3Hn/icV4e+QadzutN99OPpFrtGjuPb9uSxTfvT+a/z77N6jVrGPa3+0hJ0WgeERERKR+UlUjcfP755/x35JscMeRUel1+yi5JMUBG1cocfMZRnHjfJXw4dRJvvfVWkiIVERERyU+JcYzM7DYzm2xmW8xsfYznmJndbWbLzWyrmU0ys/0SHGrSvPXO29Tv1IIuJ/+p0HptDtmPvft05q1R75Cbm1tG0YmIiIgUTolx7NKBt4B/FeOcvwJDgGuBg4CVwAQzq1HoWbuhJUuW8NW3M+jU/4iY6nftfwSLVi5hxowZCY5MREREJDYaYxwjd78LwMwujqW+mRlwAzDM3UeEyi4CVgHnAv8X67XNrFkM1RrE2l4iLF68mB3k0KJT25jqN27XEjJSWbx4MQcddFCCoxMREREpmhLjxGkNNAI+zCtw921m9ilwGMVIjIElcY4t7vKGRFiMD9OZGZaSoqEUIiIiUm5oKEXiNAr9XBVRvirs2B6jQYMGpJHCqp9iy+HXLl1NztbtNGiQ1I5uERERkZ0qdGIcejDOi9i6l/IyHnnZKGVFaR7Ddnzpwiydfffdl3at9mH2mC9jqj9rzOc0qFGXHj16JDgyERERkdhU9KEUTwJvFFFnUQnbXhn62QhYEVbegPy9yIVy96VF1TGzpHa9mhmn9z+Ve//5AItmzqNV13YF1l2zcDlz3/uKy/qdR3p6ehlGKSIiIlKwCt1j7O6/ufvcIrasEja/kCA57pNXYGbpQC9gchzCL3dOPPFEjux+OO/e+m9++Hh6vvHD7s7C6T/y9pCn2K9JWy644IIkRSoiIiKSX0XvMY6ZmbUA6gItgFQz6xw69LO7Z4bqzAWGuvtId3czGw7camY/AT8BtwJbgNfKOv6ykJaWxrC/DWPY34cx/m+v8+Vz79OuTzeq16vJ1o2bmf/xbDYv+o1DDujGsL/dR9WqVZMdsoiIiMhOSoxjdy9wUdj+rNDPI4FJodftgFphdR4CqgBPA3WAaUBfd9+U0EiTKCMjg3vvuZdz557LqNGj+HT052Ru3kyVKlU4okt3Blw3gC5duhDMZiciIiJSfph7cZ8Dk/LIzLoCM2bMmEHXrl2THY6IiIhIgWbOnEm3bt0Aurn7zGTHk6dCjzEWEREREcmjxFhEREREBCXGIiIiIiKAHr7bk1QG+PHHH5Mdh4iIiEihwvKVysmMI5IevtsDmFkzILa1mEVERETKj2vd/alkB5FHPcZ7nn5AkSvlSak1AMaFXh8PrE5iLBWF7nnZ0z0ve7rnZU/3vOyF3/NPkxlIJCXGe56ZsSwhLaUT6qXP853ueeLpnpc93fOyp3te9nTPy17EPV+frDii0cN3IiIiIiIoMRYRERERAZQYi4iIiIgASoxFRERERAAlxiIiIiIigBJjERERERFAibGIiIiICKCV70REREREAPUYi4iIiIgASoxFRERERAAlxiIiIiIigBJjERERERFAibGIiIiICKDEWEREREQEUGIsIiIiIgIoMRYRERERAZQYi4iIiIgASoxFRERERAAlxiIiIiIigBLjPYKZXW1mC80sy8xmmNkRyY5pd2RmQ83sazPbZGarzWyUmbWLqGNmdreZLTezrWY2ycz2i6iTYWZPmNlvZrbZzMaYWbOyfTe7p9Bn4GY2PKxM9zzOzKypmf3XzH43sy1mNtvMuoUd1z2PIzNLM7P7Qv9ObzWzX8zsTjNLCauje14KZtbTzN4N3T83s/4Rx+Nyf82sjpm9YmYbQtsrZlY78e+w/CnsnptZJTN70MzmhO7lcjN72cyaRLRR7u65EuPdnJmdBQwHhgFdgM+BcWbWIplx7aZ6AU8BhwJ9gDTgQzOrFlbnr8AQ4FrgIGAlMMHMaoTVGQ4MAM4G/gRUB94zs9REv4HdmZkdBFwJfBtxSPc8jsysDvAlsAM4HugI3AisD6umex5fNwNXEdzPDgT39ybgurA6uuelUw34huD+RROv+/sa0Bk4LrR1Bl6J03vY3RR2z6sCXYG/hX6eCuwLjImoN5zyds/dXdtuvAHTgH9FlP0I3J/s2Hb3DdgLcKBnaN+AFcDNYXUyCBKKgaH9WsB24KywOk2AHODYZL+n8roR/GM4HzgGmAQM1z1P2L1+APi8kOO65/G/5+8B/4koewd4Rfc8Iffbgf5h+3G5vwS/1DhwSFidQ0Nl7ZL9vsvTPS+gzkGhei3K8z1Xj/FuzMzSgW7AhxGHPgQOK/uI9ji1Qj/Xhn62BhoRdr/dfRvwKX/c725ApYg6y4Hv0GdSmKeA9919YkS57nn8nQJMN7O3LBgyNMvMrgg7rnsef18AR5vZvgBm1omgd2xs6LjueWLF6/72ADa4+7SwOlOBDegziEUtgoR2fWi/XN7ztEQ0KmWmPpAKrIooX0Xwj4CUkJkZ8Cjwhbt/FyrOu6fR7nfLsDrb3X1dlDr6TKIws7MJvmo7KMph3fP42xsYRPDn++/AwcDjZrbN3V9G9zwRHiRICuaaWQ7Bv9u3ufvroeO654kVr/vbCFgdpf3V6DMolJlVJvi26jV33xgqLpf3XInxnsEj9i1KmRTPk8CBBL06kUpyv/WZRGFmzYF/An3dPauQqrrn8ZMCTHf3W0P7s0IPIQ0CXg6rp3seP2cB5wPnAt8TjJEcbmbL3f2lsHq654kVj/sbrb4+g0KYWSXgDYJ/e66O5RSSeM81lGL39hvBWJzI35oakP83Y4mRmT1B8HXzke6+NOzQytDPwu73SiA99IBTQXXkD90I7s0MM8s2s2yChyCvD73Ou2e65/GzAvghouxHIO+BXf05j7+HgQfc/Q13n+PurwCPAUNDx3XPEyte93cl0DBK+3uhzyCqUFL8P4LhLH3CeouhnN5zJca7MXffDswgmEEhXB9gctlHtHsLTefzJMHTs0e5+8KIKgsJ/pL2CTsnnSCRy7vfMwie9g+v0xjYH30m0XwEHEDQg5a3TQdeDb3+Bd3zePsSaBdRti/wa+i1/pzHX1UgN6Ishz/+D9Y9T6x43d8pQC0zOzisziEEw2T0GUQIS4r3AY5x998jqpTPe57sJxm1lW4j+IpuO3ApwdObjwGZQMtkx7a7bcDTBA8F9CLoWcjbqoTVuTlUZ0DoL+9rwHKgRlidfwFLgKMJptD7CJgNpCb7Pe4OG2GzUuieJ+T+HkTwn9GtQFuCr/c3A+fpnifsnr8ILAVOBFqF7usa4EHd87jd4+r88cu1A4NDr/NmQIjL/QXGEUxRdmho+xZ4N9nvv7zdc4KhuqND97MTu/6fml6e73nSb6y2OHyIwZidRcA2gt/AeiY7pt1xC/3FjrZdHFbHgLsJvo7OIniqef+IdioDTwC/A1uAd4HmyX5/u8tG/sRY9zz+9/gkYE7ofv4IXBFxXPc8vve7BsF8rb8CW4EFwH0RCYLueenuce8C/v1+MZ73F6gL/BfYGNr+C9RO9vsvb/ec4BfAgv5P7V2e77mFLioiIiIiUqFpjLGIiIiICEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEREREBFBiLCIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEpFjMzM2sf+h1q9B+5+RGtWtcUjxmdpSZzTWzcvl/opllmNliM+uW7FhE9nTl8h8BEam4zKyBmf1fKBHYZmYrzewDM+uR7NiiWAI0Br5L1AXM7O5Q0lvY1ioUx7hExVFIfC8WEFPbso6lFB4Chrl7rpndYWYrzKxueAUz62Rm282sX1kH5+7bgEeAB8v62iIVjbl7smMQEdnJzD4HKgFDgV+AhsDRwLfu/n4yY4OgZxYY4O6jyuh61YHqYUVfA88Cz4WVrXH3nLKIJ5KZvUjwGV0ScShfTGaW7u7byyq2WJjZYQS/UDR09ywzSwMmAwvc/ZxQnUrAV8A37n5xAmIo8r6YWT1gOdDZ3X+MdwwiElCPsYiUG2ZWG/gTcLO7f+Luv7r7V+5+f3hSbGZDzGyOmW02syVm9nQogcw7frGZrTezk8xsnpltMbO3zayamV1kZovMbJ2ZPWFmqWHnLQr1GL5mZplmttzMrisk3l2GUphZ79D+0WY2PXTdyWbWLuK8281stZltMrN/m9kDZjY72jXcPdPdV+ZtQA6wKbzM3XMKGOJxppl9bmZbzexrM9vXzA4KxZZpZuPNbK+I2C4xsx/NLCs0vODqGD66bRHx5MU0ycyeNLNHzew3YEKCP790M3vIzJaF2p5mZr2LiP1s4EN3zwrd72zgQqCfmZ0eqnMbUBe43sxqmdmzoc9vo5l9bGadwmJoY2ajzWxV6B5/bWbHRNzjRaE/Ay+a2QbguVDsT4Z6q7NCdYbmnePuvxMk7OfE8HmISAkpMRaR8iQztPU3s4xC6uUC1wP7AxcBRxF8HR6uaqjO2cBxQG9gBHBCaLsAuBI4PeK8m4Bvga7A/cBjZtanmO9jGHAj0B3IBp7PO2Bm5xEkWjcD3YDFwKBith+re4D7CN5LNvA6wX36M3AE0Aa4Nyy2K0Kx3wZ0AG4F/mZmF5UihotC1z4cGBgqS9Tn90LoOmcDBwJvAePNbJ9C4usJTA8vcPe5BO/9X2Z2LMG3F5cAm4D3gUahGLoBM4GP7I+hF9WBscAxQBfgA+BdM2sRcd2bCIbgdAP+FnqvpwBnAu2A84FFEed8RfC5iUiiuLs2bdq0lZsNOA1YC2wFvgT+DhxYxDlnAL+F7V8MONAmrOwZYDNQPaxsPPBM2P4iYFxE228AY8P2Hegfet0qtN85tN87tH90WP0TQmWVQ/tTgScjrvEFMDvG+7MIuCFKebS4Lgs7fnao7KiwsluAuWH7i4FzItq9HZhcSDwvEiS+mWHbW6Fjk4BZMbynUn9+BEl+LtAkou2JwN8LufZ64IIo5QZ8QtBDPzxUdhSwAciIqPszcGUh1/geuDbiMxwZUedx4CNCQxwLaOd6YGG8/q5p06Yt/6YeYxEpV9z9HaAJQe/ZBwTJ5kwzuzivjpkdaWYTQl+ZbwJeBuqZWbWwpra4+4Kw/VXAInfPjChrEBHClCj7HYr5Nr4Ne70i9DPvOu0Iev7CRe7HS3gcq0I/50SUNQAIDaloDvwnNAQg08wyCRLjNkVc5xOgc9h2fdix6ZGVE/T5dSVIZudHxN+riPirAFmRhe7uBL3nKQS97hD07lYHfo+4Ruu8a4SGezxkZj+EhoNkAu2ByB7jyPvyIsG9m2dmj5tZ3yixbiXoSReRBElLdgAiIpE8GO85IbTda2b/JhgW8KKZtST4qvoZ4A6C3uU/Af8heGgvz47IZgsoi6WDoLhPKYdfJ+/clChleayY7ZcmjsiyvLjyfl4BTItop6gH+za7+88FHQvfSeDnlxKKs1uUeDMp2G9AnQKOZUf8TCH4Rad3lLrrQz8fBo4F/kLQk7wVeBtIj6i/y31x95lm1ho4nmAYxv/MbKK7hw8VqQusKeS9iEgpKTEWkd3BD0D/0OvuBP923ejuuQBmdmYcr3VolP25cWx/HnAw8EpYWfc4tl8i7r7KzJYBe7v7qwm8VKI+v1lAKtDA3T8v5nkdY6w7k2B8cba7LyqgzhHAi+4+EnbOKtIqlsbdfSPwJvCmmb1NMD66rruvDVXZPxSviCSIEmMRKTcsmJLqLYKH1b4leNipO/BXYHSo2gKCf7uuM7N3CR62uiqOYRxuZn8FRgF9CMa/nhjH9p8gmIVgOsEsA2cRPCj2SxyvUVJ3A4+b2UaCKcwyCO5/HXd/NE7XSMjn5+7zzexV4GUzu5EggaxPMC54jruPLeDUDwgeAIzFRIKhNaPM7GaCX3KaEIwjH+Xu0wl6iU8NvTcneLCuyG8lzGwwQW/0bIKx0mcAK/mjJxqCpPuOGGMVkRLQGGMRKU8yCb7GHwx8RvDU/t8I5uy9FsDdZwNDCGZ1+A44j2DWgHj5B8HX8bMIkpAb3f2DeDUe6o29n2DBhpkE41NfJMo417Lm7v8GLid4+G0O8Gno9cI4XmM2ifv8LiEYr/wPgqR1DHAIwUIsBfkv0NEiptSLJjTu+ASCP5vPA/MJHs5sxR9juAcD6wh+6XmXIPGeGUPsmQT3ZDrBXNWtgBPCetV7ALUIhmWISIJogQ8RkRAzW0QwA8HwMr7uBGClu19QlteVgJk9BNRy94FFVk4SM3uLYIaPvyc7FpE9mYZSiIiUITOrSjB04AOCh8TOIXjYqrhzJUv8DAOuMbNUT9IKgoUJzen9DfBYsmMR2dOpx1hEJKQseozNrArBV+xdCcbwzgPuc/cRibqmiIjERomxiIiIiAh6+E5EREREBFBiLCIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEREREBFBiLCIiIiICKDEWEREREQGUGIuIiIiIAEqMRUREREQAJcYiIiIiIoASYxERERERQImxiIiIiAigxFhEREREBFBiLCIiIiICwP8DVYe6zajMCnwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -3625,17 +3613,23 @@ " ) \n", ")\n", "\n", - "ax.legend(eedgecolor=\"black\")\n", - "frame = legend.get_frame().set_linewidth(0.5) \n", + "legend = ax.legend(edgecolor=\"black\")\n", + "frame = legend.get_frame().set_linewidth(0.50) \n", "ax.set_xlim(0, 1200)\n", "ax.set_xlabel(\"Sampling Time Frame (Years)\")\n", - "ax.set_ylabel(\"Substitution Rate\\n(subsitutions/site)\")\n" + "ax.set_ylabel(\"Substitution Rate\\n(subsitutions/site)\")\n", + "for spine in ax.spines:\n", + " ax.spines[spine].set_linewidth(0.5) \n", + " \n", + "out_path = os.path.join(out_dir, \"time-dependency_clades\")\n", + "plt.savefig(out_path + \".png\", bbox_inches=\"tight\")\n", + "plt.savefig(out_path + \".svg\", bbox_inches=\"tight\")" ] }, { "cell_type": "code", "execution_count": null, - "id": "stuck-lecture", + "id": "regulation-location", "metadata": {}, "outputs": [], "source": []