InducibleHIV-Fractionation / Figures / Figure7.ipynb
Figure7.ipynb
Raw
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Setting matplotlib display setting as inline plots\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Importing necessary packages\n",
    "import os\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "plt.rcParams['font.family'] = 'Arial'\n",
    "import matplotlib\n",
    "import matplotlib.image as mpimg\n",
    "from scipy import stats\n",
    "import matplotlib.ticker as ticker\n",
    "import matplotlib.patches as mpatches\n",
    "from collections import OrderedDict\n",
    "from math import sqrt, log10"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Setting working folder\n",
    "os.chdir('C:/Users/ooma1/OneDrive - UC San Diego/Guatelli Lab/Cell Fractionation/FractPaperFigures/')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### STRING Diagram for Figure 6\n",
    "[Peroxisomal proteins from WT centroid distance analysis, labeled by STRING](https://version-11-0b.string-db.org/cgi/network?networkId=bIBUqTqaoPI4)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAksAAALDCAYAAAAFXkS3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd5xU5fXH8c/3Ttm+Sy8qLFh3sUQDRo0ajbEBURNN1V+MmthjSdXEhhqNKdZYoibGGkuiiQWwJBG7YgEbYBcQRUDq1in3/P64d3Fddoe2sLvkvF+veTFz597nPjO7w5w9z3PPIzPDOeecc861L+jqDjjnnHPOdWceLDnnnHPOFeDBknPOOedcAR4sOeecc84V4MGSc84551wBHiw555xzzhWQXNkOkgx4Dci32vyCmf1Q0lRgTzNbvK462E5/egGT2mzeFvi5mV2yvvrhnHPOuf8NWlmdpThY6m9mC9ZPl1aPpJOAbwB7m1m2q/vjnHPOuQ3LWg3DSTJJ/SQlJF0i6W1JL0q6WtKkeJ9Jkr7R6pjljyU1S7pL0huSRkmqlfRw3MZUSUet5PybA2cC3/NAyTnnnHPrwkqH4WKPSmo9DLevmc1r9fiHwEhgGyAE7l/FdtPA/Wb2LUlJ4GWiwOclSVXAM5KmmdmzHRx/AfBHM5u1iudzzjnnnFstqxosfXklw3BjgJvNrAlA0rXAyavY9hPxv1sCmwE3SGp5rgTYAVghWJI0BNiPKFBzzjnnnFsnVjVYWpkcoFaPW2ehrM1z6TbH1sX/JoAlZrZ9yxOSBgJLOjjnN4B/mtmyNeqxc84559wq6KzSAeOB/5NUFA+nHUEUJAHMB0YBSBoBbNdBG28AjZL+L953CNFVeCM72H8P4D+d0nvnnHPOuQ50VmbpRmArYApRpug9oCF+7tfATZLGAjOAx9trwMwykg4CLpf0CyAFnGVmT3Vwzi2A9zup/84555xz7Vpp6YBVakTaFxhgZrfGjy8HmszstLVu3DnnnHOuC3VWsLQxUXZpINHco5eB482so/lGzjnnnHM9QqcES84555xzGypfG84555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5acc8455wrwYMk555xzrgAPlpxzzjnnCvBgyTnnnHOuAA+WnHPOOecK8GDJOeecc64AD5Z6AEkpSR9JmthqW4WkiZJKurJvznU1ScMk5SVNjW+vSnpO0q7+OXHOdQYPlnqGg4GpwChJtQBmtgy4HTi/KzvmXDfRaGbbx7dtgYuBG/1z4pzrDB4s9QzHA/cCdwKntNp+F3CYpIFd0ivnuq++wEfxff+cOOfWigdL3ZykEcAuwN+Bm4DDJfUFMLMm4HlgTNf10LluoaTVMNxM4HLgN+CfE+fc2vNgqfs7HnjAzD4xs+eB94BjWj3/HrBVl/TMue6j9TBcNVFg9HdJw+Pn/XPinFtjHix1Y5LKgO8Bu0l6X9L7wGDgR5JS8W5ZIN9FXXSuWzKzfwNvA1+IN/nnxDm3xjxY6t4OAz4BNjKzYWY2DNgUKAe+Ge8zHJjRNd1zrnuStCUwDJgSb/LPiXNujcnMuroPrgOSpgB/NrOr2mw/H9gf2JVoeGFrM1vcBV10rstJGga8A7zaanMA/NrM7pKUxj8nzrm14MFSDybpCKIvgJ93dV+c6678c+KcW1vJru6AWzOSyoFDga93dV+c6678c+Kc6wyeWXLOOeecK8AneDvnnFtt7SwzM1XSy5KOWk/nnyqp1/o4VzvnHiaproPnzpA0S9IN67gPO0v6RFLQatvtkpolVbTadrWk30r6R6ufk8XLAk2V9GgH7W8s6T5F+ku6R9IrkqZJ+n3r88b7HyXp/gL9NUn92tl+oKQr4vtjJZ23ktc9JO7Leo1ffBjOOefcmmo0s+1bHkjaGHhN0gtm9sq6PHHr83YzPwAONbMn1/F5JgMGbAdMlZQEvgxMIroA6O/xfnsBx5jZ4y0HSjLgy2a2oED71wPnmJlJuhSYZmYHSyoGHgaOAG6Q1Ae4kOjq7cdW90WY2X3AffHDHYE+K9l/tqSpwAnAlat7vjXlmSXnnHOdwszmAG8BWwJI+oGkFyVNkfRvSTXx9hvjrMXrcdajStKtkl6LMx6/k5SUVCtpkaTt4+NulvSX+L5J6idpkKSHJb0U35avAyjprDgT8kqcWRkUb58k6WJJT0l6S9LP48cvSJouadt4v50lPa5oYeZZLefuiKQ7gU2Av0j6dnyee+I+nCRpE0n3x6/xNUk/j48bJuldSdfGfZgaZ1zGS3pH0p1tMylmFhIFLXvGm3YDXiFa3ufAuN2NgYHA06vzc5S0EzAgLoQM8E/iwCSuiP8aUB0/9y3gQ+Bnq9D0ufHvw/uSTozPdYSkB+JzHgd8W9IFhX6uwJ+BXyq60nX9MDO/+c1vfvOb31brRlTHqq7Ntl2AhcAQYA/gcaA0fm5fYHp8/0bg362Ou4loiRoBRcBDwOnxc0cDLxNlbF4GSuLtBvQDzgL+FG8rA+4AqoAjiYKEsvi5ccCD8f1JwN3x/Z3itg6IH18KXBffvx3YM75fDswHRrb32lu9lveBUa3O85dWzz0G/CS+XxW/nu/E7RlwYPzcNUTlLiqBYqJg5IvtnOt7wL/i+xcDJxMVLl4AJIDDgTvaOc6AfgV+tn8AxnXw3A7AYmD7NtuPIFptoqM2DfhpqzaagFTr4+Kf0ZXx/XZ/rq3ae4UoO7Zeft99GM4559yaKomHRCCa1rEAOMyioZKTgM2BpyW17N87HrYBaD1MNRrY1aJvwWZJfwJOBS4ys+sl7Qf8EficmTW26cODwARJQ4F/EwVZSySNBv5qZvXxfpcDZ7TKRtwT//tOq3ZaHu8Z3/8+MEbSr4AaoIQoaPpkVd8g4AlYviLDrkRBI3Efb4xf+7NEVeZb5vy8AzxtZkvjYz+k/eGpB4FL46zTAcD+ZvaRotUeRhENwY1fjb62qCEKTj4j/jncCpxkZlNXOGrl/hb/O5UoKK4ssG+7P9dWz7csYdTunKvO5sGSc865NfWZOUttJIBbzOw0gPgLfSNgUfx86wnSAVHmofXjVHxcEVHQtQTYnmiYbzkze17RGoB7EwUHk+NAKdFOm0mi7BVAc5t2su28hseJMhgPEg1v7dTq+FXV8jqDdo5d/jqBTBwstmivP59hZvMlvQccDOTM7N34qfFEw3J7AGtSX8xoM01H0k+A04HvWrSc0JrIxv22OIDu8L3s6OdqZi+2amu9LWHkc5acc86tCw8B35U0OH58HPCfAvv+SJEiosXCH4mf+z3RHJl9gT9Kqm59oKSLgLPM7F/AKcDrwDZEAc5RcUYHoiGqx83sM0FSRxRdabcjcJqZ3UM0F2lzoiBstZnZMqIMUstcnSqiYbJHCh23CiYSDVk90GrbA0RDdHPNbP4atPkGsFnLg3h+0YnAzmsRKK2KHJ8GyR39XFus1yWMPFhyzjnX6czsYeC3wCOSXiEqDnpwm+xJi5OBAURL1rxK9GV9gaSxRAVFf2RmrxLNJ7pd0ZVfLS4Dtpf0GvAC0fDMHcBfiIZvJkuaDnye6IqtVe3/YuA3wEtx26cDTxEFTGvqMOArkl4luprtHqL5W2tjItEVca2DpReAQcCENWzzH0RX1BEPW15ENHfqHn1afuCMNe9yh/4L7Cfpj3T8c0XSQKLfl6fWQR/a5UUpnXPOOfcZkh4iyuxM7uq+tCVpHDDf2qybui55Zsk555xzbR0LnK1Ws/O7A0lDiLKEf1qv5/XMknPOOedcxzyz5JxzzjlXgAdLzjnnnHMFeLDknHPOOVeAB0vOOed6hP3339+ICib6rfvdNmgeLDnnnOsRFixY0NVdcP+jPFhyzjnnnCvAgyXnnHPOuQI8WHLOOeecK8CDJeecc865ApIr38U5t5aK6+sXbRUEqW1k+YoQRHQDyAemLLJkCEok9KGZvV5cXDUbyHZhn51zzsU8WHJu3RIsLobEdldfffVPPvzgw/6/vuC8xWedfW4fDGprt2z69ne+teSWm2+rqqyoDA/5xsEvm/F3YCGwhP+BS3Kdc66787XhnFu3tGzZsr6JhB269daf++1rr75iN910Y+ob3zxEJcXFAHbQQYcE9957d/jQw4/YvvvuPTOVSN7c0JS7oW/fvh8BYRf33zkk7Q8MNbPrurgf/oW1HqxhXNCtFtztbJ5Zcq4TnXbaoLFBkNg9n7fxv/vdh08AmJmB5c3MsHDe97//vWQun68sKipJJFPJ8JF/TwylhEaNHGm9ew/YbOnSTxr79CnK4Vkl102Y2YNd3QfnupJnlpzrBBdfPNz22ONMRo36wfJtDz/8K5599tbM2Se9PKAhqa8ondwnn80KSKWTyQEokUJhNp/NZxKJIJ23UIS8owQTiourngLq8YDJrQJJRwAHAZVAP+A8YAFwAZAH3gGOBQ4DjiK6uOccYBBwKtAMvAUcE992NbNDJd0EPAc0ADXAOOAuoAooAX5hZpMkHdZOO4cBB8T7DQYuj/u4DfAzM7tX0jeBn8R9fNLMTl/J6/TPw3rgmaUVebDk3Fq64oqt7ZBD/kIiEdDy/4UE+bwB4h//OIqTj361upGlG5tl88opkUqny5rz+YSZ5RJm2byUSiLlc7awrFdiFlQsxIfg3CqKg6XDgP2A/sBkogBkZzObJ+l8YBbRRQNfN7ODJPUlCoR2MLNlki4F3jGzKyX9C1gMFJnZd+P2a4BbgL8CewMDgC3jNlZoB6gDDjWzfSV9B/gxsDOwJ3AKUdD2JDDKzBok3QLcbGaPFHid/oW1HniwtCIfhnNuLVx11eesX79+fPLJVJLJBCAkMIMwDAnDkD59enP5dSNmnnLKtM/8Z+IfPtfJHjOzEPhYUj2wBXCXJIiyOw8TBTFvxPtvCrxuZsvix48D+8b3LwKeAUa2PoGZvS7pKuB2IAVcUaCd54Ap8bbFwHQzM0mLgGJgc6LAbkLcx4q4Lee6Hf//2rk1dNFFI86srh5IZWUlc+c+Qnn5xsT/6QNgFlJX9yFVVVWEoSeJ3Do3EkDSQKJg5G3gIDNbIulAokzPUD7NWL4HjJBUZmb1wB7Am5LSwGVEw3bXSNq95QSStgUqzGyspMHA08CO7bUTH1IoRfEeMBvYx8yycfZq6lq/C86tAx4sObeGBg/udX6vXr0oLS0lmUwCS5GiOq9mhgRlZeWkUhkArr56x2dOOOH5Xbqwy27DNkjSf4jmE51AFBSNV/RLuRQ4nChYAsDMFkg6B3hUUkgUXJ0O/BZ4wMyuk7QRUZbplfiwt4BzJB0OZICzC7TznUKdNbP5ki4BHpOUAN4nmg/lXLfjc5acW0N3372fDRw4kOLiYhKJBJLaZJYMMyOXy1FXV0dd3QIOOOCRDXpc33WNljlFK5sg3RNJapl0Dm2GBd264XOWVuSZJefWUGlpKcXFxZSWlhYMlrLZLPl8nkympAt761zPFNd2ug58grfrOh4sObeGSktLKC0tpaysjFQqRRB8dqlFMyMMQ7LZLGEY8uGHwa1d1FW3gTOzG7u6D85tyHwhXefW0E9+8tjP0+n08oCprKyMioqK5bfy8nLKysooLi4mlUpxyCH3XtbVfXY9i6Qj4snZq3PMUEkHxPcvkzR0ZcespL0tJP1mLY4/VdJFrR4fIul5SZMl/TDeNkjSlWvTT+fWJQ+WnFtDL7304btSnjAMkUQqlSKVSpFOp0mn0ySTSYIgIJ/P89ZbizL9h+x0RvWIsZd2db9dz2FmN5rZfat52F7ArvHxp5rZrLXsxh+Ai1f3IEklkm4FTmy1LUE0YXxvYBfg55L6mdlcYJmkPdayr86tEz7B27k1JGkAcMYbb/z45LKyCkpKSkgmk0jCzMhkMjQ2NrJw4eL89ttfeSFw/tCtRl+PFCbMLn13xoRXu/o1uDUjKQXcAGwGJIBLzOxOSZOA+UBvomrVNwIbEV0i/yUz2ygOCM6JmyolukotQ1S7aHbc5mQzO17SOGAu8DFRIUeATeL99gauBYYAfYGJRBW2X4/b/RFRdezj4jZuJarwnQTONLP/SnoFeAzYjugy/4PMbEmr17kV8AczO0DSsA76+CPgG23eosOBZURlBTai1eRzSUkzy8WfnyeB7cysKS5LcK6ZHVzgffcvrPXAJ3ivyIMl59aQotncOwPfveGGbx2w225DhqXTRcuzSZlMM7/+9dP/ueWWZ18Crjaz9wE222yfqlw6fSlmwcwZE47owpfg1lAcIGxmZj+WVAG8RJQp+QdwuZn9U9IpwMZm9gtJNUSFGxOSTgD+ZWYfSvoV0ZfMbcCLwDCipUXeBXYiDnTM7E/xeYcRXV7/zfi4vc3sz5KKgQ/MrF/rK+Pi4O044IfAbDO7XNLGREHKZvF5DjWzpyXdBtxvZne0ep3HANVmdkZ87hX6GGeFCr1Xy/vTatvBwFXAeOBYM8vHWacPzWxggbb8C2s98GBpRT7B27k1FC2QyzOSZh511F1Tga2BCui9Jyx6iOgL5SVgopktbTnunXceWQIcNaxmzE+ra8feEJp9PHvGhF92xWtwa6wW+DdAvMTHNKLgAz6tkF0LPBjvM0PS/Hj7HOAKSXXAxsBT8fa3W6pgS/qIqLDkcpIGEQVjR5rZTEmVwI6SvkxUR6loJf29Le7LHElLiapnw6dVtme3PSfROnMft3q8Qh87yiwVGv4zs3viJVVuJMpC/TUOmHKSgrgS+QpGjhzJCy+8UOBlOrdueLDk3FqKMwR/BdJVg2u2LikaOGju+4+dAywzs2xHx70/Y8LFANU1Y2+qrh17g/L8/v03x09fX/12a2U6sDvwzziztC1RRWr4tEL2a0TZpn9J2owo8AD4M7BpHGTdxKd/kXf457ykXsC/gJ+YWcvw7RHAYjM7VtLmwDFxtjNkxfmoLf2dEmeWegOfrOy8wDyiYbQWK+xrZlcCqzQ5Ow7w7gf2NbPmeFmWMH5OQK6jQMm5ruQTvJ3rBBZp7tNrs7J0UWnOzBYWCpRamzlj/PdnDaw/xhL8fGjNmJvXdV9dp7gO6CvpSWAS0VybeW32+QswTNLjRHOJmuLttwDPSXqKaD20jVi5C+L9zpE0SdLDwH+AMZKeBq4hqq69EfAqcFC8eG2LC4G94r78CzjGzHKrcN5JRMOBnSLOsN4GPB6/d0Y0lwqigPOZzjqXc53J5yy5rhCYLenVVBduqyCxaSg2DcP8wIT0TggLEgQhhJbHwmTSFoSZVH0Y5HeSUW6hvaNkcknSgvfSZWVvAY0U/st4vRpaM/ZYycbOnD5htS73Xn78iDHflmk/gJnTxx/VSd0K6urq+iWTbGv57EYhilb8jW5hYJYlERCGljJTfZrEG+ny8rfpZu9tTyPpi0C5mT0saQvgQTPbbGXHdTeS7gd+aGYfr3TntTvP74jmTD3R0T6jRo0yH4brtnzOknOdTHV1wYCE8vv/99FHx+y2627FCxbOKx00cNCyp59+Nsjm8qSLkuEuO+28zPLpN03ZWelU+vPZXPOgkuLShrv/eW8wZsx+94dNTXcUFxfPBPJd/YJamNkOkjWs6fGzpk24E7izesRXr6quHXtDgJ59b/oD161ltxJpaXA2nz1g0qOPf3nM2DG5xyY9VtLQ2Bhsu83WTQMHDWwKgoTy+Wzyrjvv5uBDvv73sGnxncXFvWbRjd7bHuhd4PZ43bQUrS6h72F+QXRV3Wnr6gTxfKzKQoGSc13JgyXXFRQE+eJQqY333nuvIUuXLC275557E/M+nrfRX264MWhubmbex3OaQ1OuuCiobGpkeFFx+a4XXHh+7s0Zb+qCC85ZIrRtGDY/AsVrW0OmUymwCkwfrm07M6c9cCJAdc3YG6trx+yixobj339/UtPKjutIPsiXBqGq99tv3+oLfn1h+Z+uvT4x7fWXc5L4y1/+qq222MJ23mVnvn7wgcskbdPcnH2kuJjZa/s6/pfFV4l9uav7sbbMbDrrMFCKzzGX6Ko957oln7Pk1rlf/OIXY376059+pe32APJKJJaWlJbNC8Ow+Xe//31zMpkwgGQiuSCdTtaB9TGC4ffc8/clixctzt94082JPn375UpKivLZbBa6WepXpsACe6uz2ps5Y/wRsvBcKym7urp27G2rcswvfznoqyedtNGurbeZJY1AWSlYcsqpp8wNw3w+kUhmIcideOLJib332T/5+z9cTCqZzKVSyTCZTBrd7L11zrmu4nOW3Dpx+umn7zlgwIBHTzzxRIqKoiua6+vrueaaa1iwYEH1uHG/TJkFh1g+t5WhpmQq3TebyfQqKysxFDTU19fPwygtLSvtk8vlyWWb6kMjSUhZWUX5ovr6umcSWe4r6d17Np9efdTlqmtG32Fw26wZE+/v7LaH1I75dYA2Cs3mzZ4x4TOry59++sYXbrzxwF8ed9xTpFLRgr11dfO45ppdGDbsoqID999/8zDIfSNvVAs1l5QU95eC8lwmZ825TEYhKLCikrLypY2NDY9L6ftKSkrm0I3eW+d8zlK3tkH/ceXBkut0F1544Tt77LHHpsOHD19hcdkwDHnmmad56623J/z4xyf+PJ9PlCeTlsnnwzJJ6Vwup4RZLsyrTukgFYZhWTJJmEgEjblcLiEFRYlEUVM+r4+Kioo+ktTcRS+zXdU1Y+40Ur+fNePedfY/enXNmL8hNRHaX2a+MeGpP/xhMzvkkD9TVFT6mffbDMxCxo8/lWXLio465pgHXgiCfEkymcw2N2fKJaWTyWRolsiQyymvfDqRCJuh5IPi4uK53e29dc6DpW5tgw6WfM6S61S//OUv+26++eablpaWsnDhQhKJBFH5lKgqbD6fZ/DgjVi4cNGYXXf98tdeeOGFVbq8HiCVav9+d5MIcvNXvteamzljwqHV2321N1ku/vo3d7r0wDEDqKt7m4aGBEEg4rcbMwjDPIMGDSMM37yhvLx8+X9m6XTZuuyic85tUDxYcp1q2LBhC3r16sXcuXOpqKhYIVjKZDLU19fTu3dvfvCDH2TY4P4aMXt/+sR1Pul85isPLAKOuv76Ha1Xrz7MmTORysohSEGr9xuamhaSSDTSu3dvrr56OzvhhFc2sPfbOefWPQ+WXKfq3bs3VVVVFBUVLV9QtoWZkUwmqaioIJlMksutSk28nkZm63Fsu3fv3lRWVpJOpzFbzKflk6L3O502EolyEokE8YR455xzq8mDJddpLr/88qE1NTWUl5eTTqcJgk+zHC1ahuI23GBp/fnzn3c+Y+jQTSgvLyeVShV8vz1Ycs65NefBkus0/fv3bygpKaGkpISioqKVfnk3Na1x2SAHFBdXZkpLSykpKWk3ODUzzIwwDAmCwN9v55xbQx4suU5z6KGHLpg8eTJlZWUUFxcv//JuPWepdWbpf+DLO1i2bFmfpNmWRr4fyaQSoDx5kcfMwiwJCMOgKB0EzZawN4qLq2YCq5QC+r//e/j3L7xwzO/Kyso+E5y2934nEgkaGhqQpPU5TOiccxsCD5Zcp5FUNnXqVIqKiigtLSWZTLZbOiCbzRKGIaWlpUgqMrMN9RL1RFFRsEkulzu4pLh0BzNShimttJozTRZmE9lUMkUml00l0ukPMtnsHQsWLPi4X79+OTpYk616y69uaonwV4JcJttY8eDdIyguLqa4uJh0Ot1uZimTyZDP50lFlxAG+BImnUJSApgAlAEHmNmidvb5ErDYzF6R9HXgOTNb6wrvrdq/DLjEzGa12lYD/MnM9uykc/zIzK4s8Pwk4Dgzm7EKbd0I3BE/HGpmq7WUz4svvrhCtron8r9Xeh6v4O06U3777be/JJ/PI4lkMkk6naaoqIiioiJSqRSJRIIgCMjn83zuc5/7NRt20UPl82FlKpXewiz8/NnnnLvLPvuO+WJDY9OO2Uz4ha1GbL3nHXfduduihYtHnXveeV+CcGjfvkXFtLlCcHjN/ltV14y5sbp2zJ8ssNODvC6eOX3CcXPfmXT8//3fv25pyRwlEonPvN/pdHr51Yi5XI6ddvrLmfjCuJ1pMNDPzHZrL1CKHQVsFN8/BajszA6Y2amtA6V15MzObtDMHlzdQMm5ruTBkus0ZtYEvHjGGWc8snjxYurr62lsbKS5uZnm5mYaGxupr69n4cKFjB49+m/AZDPbYGYdf+tbSudDrRj85XMoKGo441e/+MTC0CyM4pUFCz7Rffc9EFRWVSiVSpWYsWldHX0BhtSM2XZozdjrqmvH/ClU4jQsee7M6ROOmzVj/DHvvzl+etzystde+/C5hx6a/vaSJUtoaGigqanpM+93XV0dixYtYostfn8VMMnMNuTgdH27DthC0rWSNpF0v6RHJL0k6WuSRgL7A7+T9H1ge+BmScdI+j1E2SlJr0oqamlUUqWkuyQ9HLd1fLx9J0nPSnpO0j2SSiRNklQjabCk/0p6FDinbUclDYuP+4ekFyT9TNKNkl6WdGG8z7aSHo3bvFtSlaQzgD6Sru6oX7Hz4vNPlNQ/bu/i+JzPSTqlTX+OkHRRfP/MuE9TJR3beT8e5zqPV+iB8Z0AACAASURBVPB2nSr+T/94YPMZM2acmEqlSKVSmBm5XI58Ps+WW255JfAq8NeeHiydeuomOw0apGe/+907GDp0Z7LZBqZMuYUnnvgDP/vZu0EmU/+5XCY81JQbKlOGRDAgMKsAhSH5poQSMlGkkObJL81YfNlVd+iVaW9/ksnkEskge+47rz9SMGsgqRw4KZ1OD5869fiji4tLSCaTmBnZbBazDFtscdnVwJPAXWbmQ3CdRNIw4A4z21nS3kDOzCZJ+iJwrpnt0zLsZGYPtgxXAXOAl4AaomBqPzM7uVW7nweGmdk9kjYCHjOzLSS9DHzHzKZLOgF4FrgkbvOHwBtmdr2kbwPHtx6Gi/v6PLA5UAK8B2wMNAAzzWygpGeBo8xsmqQfAJua2RmS5prZoAL9mkQ07HdH3K/Ngf8SZdUOIZru8WTcx58SDcMNil//ncCVwJeAIuA3wKkdzauTtEF8YW2g37s9f3y0AJ+z5DqVmTVLuhrYq6amZilQUVJSUtLU1JQzsxywCPgP8GT8uMc699yhU7797Z9uv/nm+yCJhQvfAIzNNtuD4cN357LLRoSnnjqtfyaz5JFAKjOCTIBV5c1KEyhUEDS+8NL0IU899+qBmBIfzZ1fOnzYRufeffslk1nF4TIzq5N0eSaT2W/EiMsXAaWVfQbtWr/UqvK5jycC84EHgec9q7ROfQScGQcZBnRYY97Mlkl6DNgPOBI4r80uc4FTJR0MLG3V1kAzmx63cTXQev7O1sAt8f2niP5gaetdM1sSL2PzsZktjNto+V2rBa6O20wBb65ivwAej/99Ghgbvx9PxEFPNg7ERrTTp62IMsx5osDtlHb2ca7LebDkOp2ZZYAHJf0HqGpsbKwg+gJZAizdELIbZ5+92derqwdvX1nZjwULpnzmKjQw8vmQLbbYkauu+tz8E098+TN/cQ2p/eo2svzJkkJQ2iz57Vkz7l3jSb9m1gD8U9J4oFdR+dBflFSlR3z83sdnsYG83z3A+cD1ZjZR0pHAEfH2kE+nO7S+fz1wGtGcp1fatPUz4Bkzu0bSl4mCD4APJW1hZm9JOo3PBjMzgF2Al4EdO+jjygLwN4DDzWyWpF2J5mTBpxmDjvoF8AXgX8DuwGvAdKJA8FJJKeCLwE3A6DbnnAEcLykAWibMf3UDvujD9VAeLLl1Jh5iWxDfNiibbNL7nl69ejFz5t306bNlmyt0jGy2EbNF9OnTh4svHvHhFddveoop3E+QC1AqDBNnz3rzgTmd2ac4SJ03rGbMR6GoLjDp2HW+vwNXSJoLzAb6xdufAy6S9B5R1uVmSfua2XOSNgeuaqet+4FrJB0GfALk4uHtY4EboiCbj4DL+DQTcxZwp6TvEA2xrYnj4/4l4sc/iP+dJulW4C8d9Avga5JOJco4fd/MFknaU9IzQJpoCPiltleymdlUSQ8SZcMC4BoPlFx35HOWXE8nWNKrsVG1YEMIVRRim2A2FGl2gGZh1kyCgDwJiK6bl7QkmeStdLpilqTG1TnhuHG79Nl228pPKisrOyy+GYYhYRjS0NDAggULOOuCyqvenzH+R533sjtWXTP6JJO+NGv6hG+uj/O51RdnUp4imq+0tKv7051JOgY4Jn44siv70lk20O9dn7PkXDcWLFlC32Qyv3dJSdmXJMobGhrKSoqLSpsyzQ3FxSV1+Vwu35zJBGBBUToN+ZwVFaXnN9U3PtgULn6gpKT3B6xGCYONN07+vLy8/DPLurTVUjk7kUiQyWRYX4ESAML/Mu/GJA0H/glc64HSysUlBq6DDWeCt+t5PFhyPV5xcZAKMzagub5h6+envNT3i7vsoksuvSI49tijw+bmZguCwCSRSqV15FFHJ66+5spsfd0nS0tKSuckE6kngQ9ZjWCptLSsrLS0lNLS0uWFINWmGCRAS72p0tLSzn7JBRkeLHVnZvYeURkB51wP4XWWXE8XKKsUCaX7D9pk0Pbb7xD26TsgOW/+/OyyuvocBGEqXRqGIdbcnAlDCy2XydG7T2+JZJEUpFnN9PGyZck/lpaWUlZWRnl5ORUVFcv/bbm1ZJ5KS0spLi5eRy+9fR4suc4iqa+ka9fi+K9L+lubbYm43tP+8eMSSTep7Vi2c92IZ5Zct3fyyScXbbLJJk1HHnkk/fv3B2DRokVcf/31JBKJ/ieddJwsw8JPFn74cBgqPW/uB425kEwiZSUWhulspiFUEAYQBNdfH/2/n2nONZpspjWridTqVbU+/vjxy15++USKi4spKYnqGqnNmmxhGJLJZAjDkGQyGUoK1tel+4m8ZcKEf++4TvFr2p+EvlKSLicqjzC11bbNiK6KGwL8GcDMGiU9DRweP+dct+PBkuvWzjrrrNNGjhx50X777UcYhsybNw+IApLDDz+cyZMnz7/qymsnHnf80b/J53kOUJBINslyOWu2YkQylwhDKSlaZVIN8lLyo1Rp7iNWf8mV8NFHZ7/+jW/03bql6GYy+elHqaUgZBiG5PN59tx7wqRNavY9Hzhjrd+QVelcQI+uX/W/RtIRwEFES6H0A84zs7slvUZUHqCZ6Eq1W+N9ksCZZvZfSV/l04rdU4gKVO4OXEB0LcM7RFfRDQduJFqkOUcUmGSIikIGRDWTjjOzV1v1qxLY0cxaKoi/RTQpfSvgY6KCk18H2s7H+4WZTSa6+u9f8flblANHE5VNaO0uonpgHiy5bsmDJddtjRs3Lhg6dOhFtbW1fPLJJ+0uyltdXc38+fNHn3PO+T/83e9+90TLc0VFKzTXmT459dT7btt9900uMDOVlpaSSqUIgmB5Vqm5uZmlS5ey885X/zXTUH9jZVD1bHXt6HPMlJ01Y8KF67JzoQXZwBNLPU05sA/QH5gs6d542/lmNkXSH4BHzOxySRsDT0raiqj69RfMbJ6ks4kyNtcDu8Xbzieq+ZQGXgR+QhRM9QaqiWqfHUpUMLLtunU7E9VearEpsJeZzZb0FFEg9Q/gH+29IDO7U9Kebba9DKxw9WhcaqCfpCozW7Jqb5lz648HS67bGjx4cL5Xr14sWLCAxsZGitpEQNlslmXLllFVVcWmm246h/V06aqZ5SXdMHLk1ZWHH/6FHc84Y/evtCwUbGbx8Fszm29+6ZVEQxBPv/76XTng3KFbjN6sunb02WAvzpz+4Ph10b9EIsxumFcmb9Aei4dpP5a0iChogk+DlVrgNgAzmyNpKVHRyEVmNi/efp6kAfH2u+KApAR4mCjTdBpR9mYJ8CtgIrAFcC9RxunXbfrUjyiD1GKBmc2O788GiiV9g44zS6vrY6BP3L92jRw5khdeeGENmnZu7Xiw5LqtyspKqqqqKC4uJggCwvCzo2VBEFBeXk4ymSSfX79Fqs3sY0l/uOe+uTfffPPFVwDBQQdts9miRc3Njz/+1mygkWhZl0dbL+sy662J7wDnVdfsf3J17eizh89ovHCSTerUYTOFljefK9vTjASQNJAowzMv3t7ySz+dKCM0Jc4s9SYqTNlLUh8zWyjpCqKhug+Ag+KlTQ4E6oiG+Z4ws3MlfZcocLoF+MjM9pW0C3Ah8OVWfZoH9Gr1eIUQvFBmaQ30Ilqex7lux4Ml121VVlZSVlZGUVERiUSivdQ9+Xy+S4Kl+PyfVNeOnly3eNYlwMB7732tF9GX2zxgXlxRu10zZzx4BUB17ehzhtaM0awZE8Z1Vr/CUDkSK9/PdSuDWpYHAk6Is5etn7+QqHr3N4iyRceYWSZeuHa8pDzRnKXniap6j48LXy4lmp9UAdwqKUf0O/pjYCZR1e9TieY3tV2j7lngt+vm5X6WpF7AYjOrWx/nc251ebDkuq2WWkYdVckGlgdLmUyHcck6M7RmzLmbDmo8f+Z0ywHL1qSNmdMnnhu3NQ4smDVj4tlr268gQc5XzO1xHjOz01tvMLNhre4vBL7W9iAzm0g0nNbaw/GttXlEa8e1tXdHHYoXaX5e0g5mNsXMBrV67jsdHdemjUnApHa2H9Fm06HA1avSpnNdwYMl122VlJRQVlZGSUnJ8mCp9eX5rSdTNzd3QWmhACZN6pwhtFkzJozrV3NQRXXt6LMNls6aPvGyNW3LSGZY8YI42dKlfZoC+4oF2iYIeBcScy0b1pPEwjAsL0omlbOcQaIp29iYT6VSiYDk/GXNzbP79u27jJUvxOo2PGcTzXc6el2dQFIJsCvwvXV1DufWlgdLrluSlJ4yZcryWkYtV5u1rWWUzWYxM1KpFJLKzKx+ffRv2Igx586aNuGcle+56hbMuHcZcN7GI0bvUF07+myTPTpr2oNPrPTANhJhPpsPVlgWIpktTgwhm98d2DXMW0NJaemyhtziRvJYaVmv4oa6xUWl5b1y2Wx9c5BK5orKqqyhYemrFRXF9wKvE00Cdp3MzG7s6j50JJ48vs4CpfgcjcBh6/Iczq0tD5Zcd2U77LDD7959991fBEGwvJZR62Apl8sRhiFmxl577XUl0byL9dO5dVhgcs60iVOAKdVbjzm6unb0l8O68t/Pnv33VV7sN1wxs6TFixeXFSW19cIli3caMmT4Ds1Ndc2NDXVheUWfEoD6+kVL33z73fTcuXOD/fbdJ5csLQmbm5flMplMb8leSaXsrXi+i2eXnHP/c3y5E9ctmVkWeHHcuHH/Xrp0KfX19TQ1NS0fcmtqaqKhoYHFixdzxx13vLJkyZIXWU/LfAwdMebcWTMq2l5m3elmvj7h+pnTJ54XlNf/rLp29CpnsYLe7y9qbEq0zgIFUnFfC2yb4cO32HHw4MEmoXyYTe68005hVVWVmSnYeKPBCmRBU1NzUrIkJFJTp0ypTEopX4rCOfe/zDNLrju79+abb9745ptvfnv69OnHFRcXL6+UncvlaGpqYtSoUdfX19e/Ddxltn6qCwWYzO5ab1msmdMnnA/RlXPR42hSeFtDh57fVFu7SdEjtz1G//51PPvsVYc++eTFbL31VcVf+cqXynNNCpYsnDfRRKqpuan8rbfeLnr6mScWNzU12eSnn+szcqeRuQEDByVzuWxdQ0OYffaZZ0t33mWnZ/PoA6JUlWeVnHP/k7Sevl+cWyOSUsBuwFeICta1DH+JqHjdY8AkM1t/WaVOnqu0OiRpaM3+Z4ehwtlvREEUQG3tZXbXXYciBUgwcODHQDQJftq0+5kyZeLzRx/9wFlSUAp5wlCliQTkctaYTqfzYTZbHipMAqGkRkhkgWQul/1Ilak3K6hYyOovC+Ncpxo1apR5Ucpua4POPntmyXVr8XDco5KeICpaVxE/tYyoLst6XQdN62kx3I7E2bNzN9lu9CbVtaPPxvR6r6Kx/xg1agumTPl0SZgvfOElzISZ0a/fRlRUZHa88MIdFl1wwbsPtdvwul0exjnnejQPllyPEAdFC+JblxhWO/r8WTOGti3c1yU+eGXiB8B5W27126V9+/ZlyRJj2rQly8srbLXVtOXB0oIF0+jVqxdDhmSfYwP/688559YFD5bcBsvMipYtW1adIL+lkkpHFbVVCwyXghlhNjdXSTUkk0WBstkSgFwqbEgq9X4qVfaupAZazdMJLRGaXdutLp/v329wRWVlJaWlpcyZk10eLOVyH5DPR3FReXk5qVSqS6qcO+fchsCDJbehEiwsUjK9TTKRPiSZDPo1NTY3Ab2CQH3DPCOLSkrqmjPNmTDMB4lEkFYyQVEQLGtqan7OLPOPoqKid4ivwa+uHXPW7DcmjuvSV9SOqqpKKioqlq+f1zIMV1paShgGy5eESafTmBl/+tMufzjuuGd+1sXdds65HsWDJbcBSyZTMPBLu3/5gP/895HEqFE7F49/4L7mzbfYqmTxovnZRD4R9u0zoKhP3752zNE/DK+97rrgwzmzlpaUlCQymcxzUDST5QWLTOvrarvV0XpJmEQiQRBEAVJJScnyYKmleGc2m6WysqzXylt1zjnXmgdLbgNWSchSKcDKyiqDl16cnBswcKOSz3/+8yHI8rmstt5m63DrESPYc88vhdded13Q3JxJTH15au/ttvtcEXEdsmG1Y86fNbCx3cv1u1pJSQmlpaXLl4RpqXJeUVGxfBgun8+TyWTI5XIsWDDixC7u8gZP0v7AUDO7bj2c6wigpu26cuvoXMOAO8xs53V9ro4sW/Yq//nvZl11+hV8Za93uroLbj3xYMn1aOPGbTJ3yJCagd///v0kkyXMnfsqt932TebPbzjytNNe/mdxcTD/8ccffai+fslGqSCpjz+e05xUkMuG+dJUKqWnnnoskwgSAUqkZs58Vwpo2Hbb7V6A8CMgJ0lDthodWietAdeZJAVHHvkMxcXFlJaWfmZJmJZgqaXSeSKRoKGhgVNO+WPlySdfMb+r+74hM7MHu7oPzrnO5cGS67EuuaTGTjhhImAsWvQ2AIkEfO97dyEl/nrrrbv99YjvPlYb9Eo1JRIqy+bzSgVBJp/Ph4mkijLZZgVBqimXySZC5YoAArOmokR6Zioz713olR2y1X7ndMe5SjH76193Oeu00945P5VKUVxcTCKRQNLyYbh8Pk82myWTyfDqq/Pe6+oOd7U4E3MQUAn0A84zs7sl7UG0YGweeAc4lmi9sqOIMoznAIOAU4kqxb8FHBPfdjWzQyXdBDwHNAA1wDjgLqAKKAF+YWaTJB3WTjuHAQfE+w0GLo/7uQ3wMzO7V9KPgIOBFFGNsYMLvMaVtfVN4Cfx633SzE6XNA7YPH5f+gBXA4cAWwLfB+YC/SXdBwwAxpvZ+ZJuBPrGtwOA3wJD4scTzeyseJ9mYFjcpyPM7KX2+lHgx+dcl/FgyfVIV165jfXv35cFC14gCBLL14wDMAvJ5fJUVfXmT3/d5ZWf//yt9Gr/ppcOAyAIkt1yrhJENZckvXDv/SfP+NqBf6ypqKigqKiIIAjIZDLk8yKbzVJfX09d3RK+971/3AAs7Op+dwPlwD5Af2By/OV/PbCbmc2TdD5wBNHCwYvM7CBJfYkCoR3MbJmkS4FjzexKSfvEwUDazK6OgxWAzYgCrL2Jgost43bObdsOUAdUmNm+kr4D/BjYGdgTOEXS/UTBx95mFkp6CNixwGss1NYTcR9GmVmDpFsk7RMf12hm+0s6HRhjZgdIOhL4DnBZ/N59L+7vE/F7B/BfM7s0Hqp71sx+KKkY+AA4K95nppkdK+lo4BhJv2qvH2b2yEp+fs6tdx4suR7nwgtrzqyuHkhlZSVz5jxMRcUmbYIlY/Hit6mqqiIMw9SanmfoiNHnzZo28exO6fQ6MmSrfXaZPefV+RddtOmjp546/fiyslKSySRLly4lmzUaGxu5++6XX/nlLx95GLjBzLx+ADwWL4T8saRFwEZE2Y674t+jEuBhogzTG/ExmwKvm9my+PHjwL7x/YuAZ4CRrU9iZq9Lugq4nSgbdEWBdp4DpsTbFgPT42B4EVAcB0gZ4HZJdcAmcZsd6bAtouxRf2BC/Hor4n4BvNTquGnx/ZbjAF42syUAkiYTZZ1o9T4tBHaU9GVgKZ8td9rSp9nArivph3PdigdLrsfZZJPe51dVVVFWVkYymSSfn7c8WGpJApWVlZNOZzEzrrrq8y+eeOJLIwu12VY0V2m/br28R3Xt6LNFsqJ+yay9gb2uuvaLW2QbF83eeONtSx95ZFnquedmLsjnrQl4HfiHmXlWKTISQNJAouG4D+LbQWa2RNKBRJmToXy6xMt7wAhJZWZWD+wBvCkpTZRxORa4RtLuLSeRtC1RhmespMHA00TZoBXaiQ/pMIMpaTvga2a2k6RS4EUKFxgtlA19jyhg2cfMsnEmbCrwtZUcB1ArqRxoAnYCrgPG8un7dARRZf1jJW1OlEFq6Wfbtjvqh3PdjgdLrscpLa2gvLyc0tLS5XN02maWwjBcXluof/+Gz6/uOaprR497vwvXgFuZ4TX7H0sQ9Jo5fcJPACQ9NGDjz+825+3/3DJnzquD5swhIMoIzAYWdtehxC4ySNJ/iOYSnWBmeUmnAOMlBUQZkcOJgiUAzGyBpHOIlt4JgbeB04nm5zxgZtdJ2ogoy/RKfNhbwDmSDgcywNkF2vnOSvr8NlAv6QWiuT8fEWXEVpuZzZd0CfCYpATwPtHcqlWxELiTKCN0p5lNa/3ZA/4D3BEHjfVE70G7/VzLfji3XvlCuq7HmTDha9avXz9KS0uXz1cKgk//ww5Di+ct5aivr+eTT+bx1a+OX61lPobVjjn//ekTzlr5nutfde2YszB6z5wRBUoA0rcSQ7esP3PmG+O7ZYmD7mJ9XmrvOoeklon0DBiQHPm324eu5Ij1x0sHfMYGvZSSZ5Zcj1NeHhVibBmGa6la3VrLVWBhGGJWvFqXyg8dMebcWd02UBp9Dqaq1oESwNCa+jNnTvdAyW144npV1wFstVWR/3XvuoQHS67H+eCD5ucHDy7asaSkhOLi4jizFAVMZrZ8iY9MJkMmk+HAA+8buzoJVHXTSdDDa0efY6jPrBkTTunqvvRUZnZjV/fBOdfzrPgnuXPd3GGH3XNtGIYkEgla6gu13FoCqFQqFQdRBuRWOT08dMTo82ZOn3jeOuz+GqmuHX12aOo1a/qKgdKwmq/uKQWTu6JfrnNIKpb0w67uR1txbaf2tu8ez/Na3fY+L2mOpEnx7duSSiTdpDaTn5zrTjxYcj3RMzU1l1y1aNEi6uvraWpqIpPJkM1maW5uprGxkbq6OhYtWkhNzeWXATNWtWGZut1nYmjN/uMM9Zs5Y8KP23veZHu8P+3+ieu7X65TDQK6XbAEnNl2QxzUjAOuWYP2Pg9cYmZ7xrc7zayR6ErBw9eqp86tQz4M53qiGcCErba6WK++etIJLQvJSiKfz9Pc3Mwbb8xb+NWv3nQbcIuZLV2VRofWjBk3a8aEFb4culJ1zdgzRVA1c/qEkzvcyT5Ty8b1TGcQlRQ4m+iP2C8SFYD8AVEQMYqoDtF0MzsyrrY9nKjYZTXwYzN7SNIFwF5xG7eb2WWSJhF9ZmqIJuF+28zmSroY2C0+/9/M7PI21bjHA30kXW1mJ7Tq6z7ANDPLSNoTOI3oar/hRFfIXSDp163abrEvUdmGrSQdRHSl3Klxzam7gAeBm9bqXXRuHfFgyfU4cYG+icDMbbf944dElyZrxx2rBzz//MwPiL4o3gTuNbPZq9qupG41eXR47ZizJPq/XyBQGlYz9gw19e92w4ZutV0AbGtm58WB0HQzO0VSJVEV8X3isgavS9o4PqbZzEbH1bd/CjxEFFh9CfiQqOZRi6fN7DhJJwC/kvQwUXCzM9H3wJOS/hvv+18zuxRA0kltAiWIKoG/0upxNbAdUQHKD4ELzKzdPzriQpZ/NrMXJZ1BtIzMz8xskaR+kqpail461514sOR6pLhu0OuSpgO9gf7PPz+ziKi2y8dA3erUFho6YvR5s6Z3n2rd1bVjzgpEn/enFcgoASaSM9+/sWl99cutNy0VsRuBAZJuJyqUWc6nlbtbV8RuqbD9HeA3RMN6rYdmWwKhp4nWiZsNPBF/RrKSngVGtDl3R/oBz7Z6/KqZ5YgWnm4EKJBZ+qeZLY4f/xP4Y6vnPyZak86DJdfteLDkerR42YpP4tsa605zlaprx54jrM970yYWnEA7fMTYakP+xbJhCPnsHNKWitijgSFm9m1J/YGv82k9m8/8MSCpCPgm8N14n9cl3RE/PZKoSvmuRBXdpwNHApdKShEN+90Un6915fr2Jl3PA3q1erzCHyUFMktPxNmqycBXiCqRt+gFFCzzUVGxLV/Z64VCuzi3Tniw5P7ndae5SsNq/5+9Mw9vqtra+LvOydx0oi1zkyJCkyIigoL3OsB16KRerp/igAoqgnIRUcEREQQVBUFxQEARFCdwuIodcMQRFRQZSlLGpmWmdG6TNDlnfX+ckxJKC2UoLXp+z9Onyc4+a6+TTqtrr/2ujOcIbChwHTlQAgBJ4uGF+VmtUg9K45jZB8BARM9CySaF+A3A42rmxw9gGxpXxPYTUQmUliGlUPrbFaovDyOi+6FkXm9h5gNENICIVgIwAFjCzH80cCBtIxEtZuabw8ZWQAna3jqO+7wbwMtqn7s9UMUmiSgGSpuUquOwqaHR7GjBksbfntZSq2RPyZwgEPTbN+aMbdIFhFapB6Vx7DCzD8A5DYzvgdJPrj4/hc1xQ6kjAjM/CeCQGjY1AHpEnRdue1wD6w2r93xgA2t/A2ACERmYeQWU4Ck0v30D88Pt/QEli1WfmwC8eqRrNTRaEi1Y0jhRqKKioo1OJw8UQP1loA0xYpk4lpgOgFAqM1cTi/4gZCbmVVYrfUMUXYqjN+1sduzO9Kc8rpzHWtoPmyNjksAUsz0/u0mBki05c7DA9L/m9ktDoz7qAYvJAEZBaSJ8QhCRGcr24C0naktDo7loNXUaGqcrk0iv11uIxR7miJgH9u7dc43MuNTp7HWxDFxeVlZ63WefLbtNbxBvtnW23y8I1MPnM1jRCvoI0aRJApHQ4tkZe0rmBBIQXZCf1bSMEgAS4CzYlKV1aNc4KqqeUZO1xppo81tmPuFASbXlZeYhav2hhkarRMssaZwEfILEMANA//4XRy94Y7783YovJZZly8B/pYkmk5GvvPJKMwCBSEpi5gQoXdNb9Jej7YNVjxe4slr0BFyXHplPEMttmlKjFKJtjwFWMyJaPCunoaGh8XdBC5Y0msw0u50FUcTS6+IhgOCt8Gavm7PuSuaaoChKZQCwdesmz5Aht8S++95bO64ZdF2XDRv+3HfgQDGxLFF5eXmiQW8x+/y+aMAiAgi01L0MHjxY5BY+AdfFmfE4g48pUAIAk2Qd73F//kRz+aWhoaGhcSh0DFI0Gn9TZnTtyiN/+AGQlURQtVWJsRmMwdMH47mbJsT17tr/3GCQexGzVRCoRjmZI0TIkPUAWIAgyLLMol63iwL4yRgZuQUtSwahBgAAIABJREFUGCzZnZlPeFxZk1tu/fSJJKNNQX4Ti7nDsDkyJhe6s7VgSeNvR9++fXn1ak06oJXS4qUVzYmWWdI4Ii84nZxy2WXY+9tvEIhAgoAtZ5oBmSGzjEcHP4qJH8848OXjXxKAr5pk1HT0Kc0JUV+9PbldsKXWt6dkTiBwfEF+7hEFJxuiszPtGpbFT5rDryNAVVVVCUDtpYKg78MyW0AAWFZ+OZIgQ5ZlCJCJyUfEJUSwMJMfgvBjZWXhH23b9qhGKyjo19DQ0DgetGBJo1Ge6d59XOd27RDYsQNl69eDoBxDdkVGgxhgZrh3u9GmTRtc+uSl/PXEr0+L/yxsjnaPFLiyWqRFSBdn+hNgiilwH3ugBAAixJ6eTVkfn2y/jgLpdFIkZH3v+W8sGHrnncPIaDBhf/EBoXPnpNhAIIC89WtKZWZ57H3jLZ99+uE2izW2R5ekLrV5eWsft1q7bAFQAy1YAgAQ0TAADmZ++BStV6Cud1SldyJ6AUqj28ImzH0KQCqUdiUrGnjdBqAXMy8jop4AYpn5+2P1X0OjNaCdhtNolHaxsdOjo6NhtVpR5XKh2u1GdX4+tpVtw9bSrdhauhWiQUR0dDTaxLZpaXebBHVLNxKReAqWEvbsWRtRWVmZwFVVbQHED0y789UuSR07e9xLn62qqmrHXJnAzAnMHAVAxFHS2IMHDxbB1BJ5OaJaEiVJsixdujR6zD0PtAkEgrE52bkxTocDAJBos0Wfd/4/4gwGnREk2CrKSn233nqLn2S5rSAEWjiXqNFUmHlsUwIllesBDGwoUFL5FxRJAAD4Pxxsp6KhcdqhZZY0GiUqOhpRUVEwm80QBCEkbgez2YyQjKMkSdDpdJBlGVc+fWXO549+nt6CLh8Vmw4PF2xs9uJogbkyLtqadKlMfFE1yeZfv/3p3CkThht6905eV1MdeIFYDlRXU0AvVgclSd4oysg2RkVtA9Do9uAv66oeL3RnP9JcTt91112dEhMTd4S+1rIso6Ki4txp06atZYMhSEFfxTXXDDpw98g7Jb1Ohzlz57bbvWc3ARD79f+nXFZWvCsurm3nCy640LR16za9I7k7P/jgAz5BSSidFlnHUw0RPQCln1sQwPfM/JDaSLcLgLZQmtTex8zLiehKKIKT5VAUutcx86R69q6E0pwWUHrH3aU+nkNEXdTH/wEgAXgdSouReADzmXkOEa1Qr7mhIR/C1pkIoDOALCJKB/AigEQAcVB60k0C8DAAi9q/cRiAWiL6A8DLzHy+aucDADOYedVxvoUaGqcELVjSaJSIiAhYLBaYTKYGgyVmBjNDp9NBkiREWCN6trDLR4SIyObIOBVbQUIgQO1l4otliTOmPftsxwsuuEAakHpphc5gvV4URQQCvgOAzIHaYHDPrqLEtu06bjAqrSkaDZaaU2l8+vTpPHPmzMPGRVH847nnnsP48ePjvLJv9YgRwyWJYarx+ejrr75Q/GKwDMiBQJD37dsjoC4wknyQ5D+Dsq8Kk2YxJk1qLvdPS9StqcFQFK2DAD5Sgx0A8DNzOhFdDuABIvoKwGwAFzDzXiJ6pwF7OgAvAzifmfeFBTQA8AYz/0hECwFcDmALgPeZ+WMi6gjgOwBz6pk8xAcAdcESMz9JRLdDaY7bHsAvzDyciEwAdjDz40Q0Dcr23yI1UNvDzL8RkZeIUqC0O+miBUoapwNasKTRKGazGRaL5ZDMEhHBarXWBUuyLEOv10OSJBhE47ct7fMRoJTe10zd8PuHL1dXl/WWJE40GMQ2gLjbLwdLLYBUyywLgtEqy4EuOh1VSFLtHp3OWCNJtaUmU+w+tZ9VUwIWgSTZKgpCQklJSfunn56mj42N0Xk8BdaK8gP+qOg4Y0JCuzbbt22pve22O/Rvvf1mlAghEkfYFu92VnpXYt3ek/ZuhPHiiy/yJZdcgj179hwSFIeC4cGDB+Oll146MGbMGAJwfPVSWqDUEA4oQUYAUJrMAuihvrZG/VwE5UhEAoAKZg59D/wAoD0RjQZwrTo2BEApM+8D6lqfhL6eoYa1ewBY1M9jiegaABUA9A34V9+HxigBcB4RDVRtGY9y3/OhZJoKASwmojOhZLkA4G1mfuMo12tonHK0YEmjQYiIVt9+O0wmEywWC3Q6HQRB+VteP1jy+Xzw+XxYNmmLbH85cw7ABJmYzPpHCv78X1kL3woAYOHChcaePc+M8VWX32Q2mfr6/P6OLCOKubbMKOqrJZBsEIWA5K/VizqhgxSUa/V6SwmYi2VZ+L66el+21dpuP8KCpc6dH1tIZB3aqVM0gkEZe/ZUAKgcvGPHM5/4pWDljyvXdbhx8DWG7t26ef9cu3p7bGxbx9df5e6Mj49rX7B982pRpzN/9fU3Z/3yy6+1/fv3P+LpPL9Ewwrdy05609wpU6YsTEpKwq5duyCKIkRRREhOJLQVt3fvXsTHx2P27NmsBkwaJwc3lKyRDsq22MVQmtP2wuFB+T4AkUSUwMz7AfQHUMDML0PJJoGIBAAxRNSGmUuIaDaAxer19e2NA7BS3XobCCCzAf+amskcBqUJ7kg18BlBSoQm4+A/AOGPP1TXPwDgOmYugdrb7mis37seXWd0baJbx8fWcVub1b7G6YkWLGk0SubSpU/9Mnr0Y2azGQaDAaIogogQGRkJYkCWZQQCAUiSBFOECeXFK6eW7efNAJDkzLxGrg08Z3dmSswQSMAWz8as6S10K/TS68snrP7h3azff//tkbPPTunx4UefWIv376dzeveWbrnlVuOuXbvh9VZVGywRvPKnnwxn9zqbH5/4KJ55ZurvRPQ7qq0E60GDZ5wxnVeufBDhpTjMgNEoLklOnoGYjn9O7Xu2Y+/OnZ5RgKznYDBYXLzHQgSdx7PNJzAqAJh37yyIkIEDomTciiNswYGEZlE779ix49Do6GiYzWZUVFRAFA+tfZdlGQaDAdHR0QgGW0xt4S8JM68noiVQmuIKAH4E8D8owVL9ubKaRcomonJ1/uYG5oyCUkckQckMNbbFtQxKHdMQKEFLkIiOlhFqjK8BvE9EFwGoVv3qCGA9gMfUOqXfAUwnIhczf0tE3wNIUAMlDY1WjxYsaTQIMzMRufaVlAQNBoOOmesCJlmWwZIMSZLg9XpRUVGBxy587HkoqX0AQIEr62OEbdnYnZkTbSmZ84ghMUOQCdt2urOf41OgikpEOPfCG7m6ujRy6dIPOzqSHzeVl5WKJBAu/Gd/EgQl4Fm2bJn56quulO7+72jdqy+/FJRlWQoEAgWiiDxjgqUS6n/avXq9wrNmDcKqVfvVbau69wyyzBg58mIsWGCY8MqLJycLk5SceZ8+KDTL1kRsbCyio6NhNBoP2YILUX+r9YUXXogZO3Zsq8gWnq4w88KwxzMB1C8WmxT2uhsHsy7nALhQEXylxVC2x+rbzoFSYB1OUtjr4XIFjgbcC63VmA/ha4Xs5gFoqF5xJ4DksOdZYY91ULbjNDROC7RgSeNIfN5v7twZH9544y3nJCZ2Cm3HVVWZwJIMv9+Pam81HrvwsVcBfMLMlY0Z8oTpGhFNEuzJqyfaHBmv2p2ZAGAUiJdv35j9wYk4O6pduzPsJtPWHoMGIemii7AvLw+/zJ2Lq7uc/emnP7x7XU1Nee+JEycsmzRp0tXDh9/hKzlwQAzWBmvvvvuuaEmSKTMtrWzy5KmRP//84463Fi1qe8Vll9YA7KqtFXeYTOQDwGefPZcTEuIxd64L/fp1qNuaBA7W+KxatQfx8fFwOJ7b5nY/eMaJ3BMAsICorVuXNfU49zFhtVphNpsPK+IPR5Ik6PV6BINBWK3W7gB+aw5fNI5KJYBfiKgGQAGAE/p5aSmI6AsAu5j5m5b2RUOjqWjtTjSOCBHFAxgK4IynUlP/kXHmmb0+u74d1ZRUe5+/9vl3gsFgKZQahFXHmyVKShpmYvP+aQCboFRDkUjCO9s2fv5DU208k5Qk3ZGVJYg6XV0hOqAEMCDCgowMjFu1KrbGIDmIdF1lmS10aGRAAJiI/AB7ISt6RjqBturN8kYgugwAX3LJOxwbG1tX9B4eLAHKtlUwGITP50NxcTF++mnoCWWXiIgSHRkPFrqynj0RO43xww8/cJs2bRoNlkKZJb/fj/Lycqxfvz5mxIgR5c3hi4bG0TAmGrnzvZ2PPvEE0GqWjpu/dD2jFixpHBVVxLEDlFqKzlBOzpQC2AjAzczek7me7az0rpDERwjKKSEwzJJMs3Zs+nxtQ/Of7dq1sG3bton9Ro+uq6sixXGwLEOWJPz66quoKC7GPfn5x/0D3b//i93atu266WjbVpIkwe/3o6SkBNnZg07oF4jdkfm0x5316InYaAwi0v3666+B+Ph4WCwWiKJ4WPAXup9Q8Ne3b19zU5SgNY5OSNNI3eYKH48FMANANyhipUUARjLzEYNUInoIwFgox/FP2teIiEYAeDN0aq+B1xdCkSHIbYKtSVC2638BcHXoxF5T0YKlVs1fOljStuE0jgozSwB2qB/NTuGGnK0Ahoeed+me2gui7i5ly44JIF9tUH5q9+ac/QCQkJCQGBsbix3LliH6jDMgUHjRNcN74ABiYmJARBg8eLC4ZMkS6Vh9ogEDdOdR2nkh7amjBUuiKKK2tvZ434KwheE/cSONovvqq69+vvnmm/9hMpnqatLCkWUZtbW1kCQJUGq2TAC0YKl5eQ/AXGb+BACI6D4Ac6EIRR6JIQDeV+ctPIn+PArllN5Ja3zNzH8C+PNk2dPQaG60YEmj1bN90/K1AO4OPe+SnPlvg0jT7c5M73D93v87q1MCrFYrDDodAjuUeI6IDh6Bl2VYrcpRtku3bvWj3ve9PfnqLhL5zyWI/wSzjUE+Ij7k9FlnmGiv8L3lLPMARERENClY8vlOLKZIcqQOJDIuP/rM48b/2GOPLb3uuuv+4fV66+4llF0KbSn6/X5UVlaib9++0wBUNaM/pxVEFIWGVbD7QVG0JihFzrcB+ANAd2aWiOhZAKsbsWkH0D4UKKnMBpSzmOrptbEA/FBOnY1g5gARDQCwFcBrUOQCFjZgewuAn6FkrL4BEA3gfAD5zHwLESUCmIeDAfEIHBSdfJ+I/g9K0Fan1M3MITmLUUQ0HsrP1h3MvKUhdfIwXwZAyardQER3QPn5FgF82oAq+QjVF+hitD9ZGi2D9p2ncdqxPT/rUwCfAsDSgQPvioqKgtVqrdtGonqZpVAAIwgCav1+MdGR/tZBa0wMgQiCzCSX6QV5NtdY1hUUHK4PRUQRhounXmM2m2E2mw/RngqtFQowiAgGgwFEpGPm4zpzz6S72OP6bPLxXNsk+8qJx6Xdu3dPWr9+/b1+vx9ms7luKzMUKNXU1MDpdL4K4MPjvZe/KGeiYRXseQBuYGaXepS/KxRZgFQiWg4gHcDjAP7bgM2OALaHD6iZ3XIiigMwGUBvZq4kolkARkLRWRoO4HVmziciPxH1Y+Zf69lOgtKvbTcUIcl+AO4BsI2IYqBs/c1m5hwiuhTANGYeQkSPQwl6ElFPqVu9DwD4mZmnEVEGgOeI6Ak0rk5eBxG1hdIW5WwoAeDzRGRl5rqgnJnnqe8pjIlGrW5Eo0XQgiWN05pQ4BL6I1+/5gY4mO2RZRlmsxlF7mW3Hs9anbsOXOT17avV69sbQttW9RWvJUlCbW0tamtrsWbNa//DCezjy6CGVJVPNrsAzO3Zs6fXZDLF/fzzz7dbrVaRmVFTU8OXXXbZogMHDpRCkYFYcxRbfzcaU8Fux8wuAGDmVwGAiOYDGANFH+krZq4Nfd8Q0etQAq/9ULJGhxTlEJEewHVQMkl5YadOvwdwhVrjlAGgLRHdAyVjNBpA/WDpQKhJLhFVM/NG9XE5lGxSTwCPqrVPBKD+PvKRlLq/Vz//DGA6jqxOHs4ZADaE1T3e18AcDY0WRwuWNE5vjMZSi8USa7FYoNfrG8wshbI9zIwIixl2Z+Y7CMpPeTbnbGzKEjZH+icg1JCo3/nWovRnxo51PREKysJVr0PrVFdXo6KiElsLViQ2VhR7NOwpmeMeHj1sWmXlgRSdTt+HWY4DBGKWk5nJIAi0izlYoSP9VgniLxaLZTcUleRjQj3B6CKiyT6f7+xzzz33dyhbSoBSxL8ewJ9HkoX4G9OYCvYuIurGzJvVwGMTM39CRC8CuAPAhHAjzDw8/DkRFRPRv5n5U3XoXijbZaMApBBRBDNXA7gEwCYAN0Pp/TZevd4CYHuY2nfdUke5HzeUprY/E5FDtQ8cVN8ehoaVuqH69zOAiwBswJHVycPZCsBBREZVP+pDAPcy886j+KqhcUrRgiWN05rrPv8884/Ro3+OiIioy/SEZ5dkWQYz1xUp59V4lwJcDpEmJDkzZABgZoF0/GHBhtxD+p7ZHGnvgUgGobLQlXMrABgMEb1em3fRtXfe8X0Pq9VXJyEAoE4y4MCBwpr58wfMaZv0zyK7I+19iPSRJy9n6bHclyiI1mHDMiKJgwPnzHntrilTnnHu2lmQ/8GSj9oOuek6qWOnM+KZWd65s2C2HPTlA5a9OI5gKYR6euo3AL+F/gCeCsHQ05zGVLBHAlhARDKULa8X1PnvQGnvkXcUu7cAeIWIxgEwQAko7mTmcnV761vV9hYoW1i/qtcAAJi5hog+AnAngKeP4X7GqfdjAmCGEqQBSh+6bCjbhg0pdQNAfyL6BkpAdjsze5qiTs7M+9Uaru9IaRS9TAuUNFojmnSAxmkNEV3kGjv2+zZxcYiIiIBer6/TWWJmMIBAbS28Xi9KSkqQ/Pzz1zPzknAbZzjThkssZBChBoAsAZcCvEsAVntcOXeHrSUkJqe+VpS/fBqAYYJgSh80aHZyYuI/I2U5iG3bvtz49ddPCzKkFbW+8ieZeTcA2B1p7zPYUOhefk1T7snuTOugF01D8n5/6zNZ5lG+at+1OqMpFiyLM2bN0keYzbhn9OigyajfVOnzLdILxg/MZvMuKP/Ba7RSiOhBAMXMvKClfTld0aQDWjV/aekALVjSOK0hogQAE/IfeGBMqMdZqEA5tDXm8/lQUVGB7rMXuOPaJn23z7PqroZs2ZLT3oAAA4B4AcKB0DgzRAIXM8lWM5vGuN2fVhKRvp29/3t7Pb+8BSAhKr7rbQFv1ftSoOaLdmdc+FChK/uOcNt2R+oSAPC4lw+uv+4jjzwyxmw2T2ZmvdfrffedT9buLMrPmewtLbXJBuHWyMjYKT5v+T4iEkTREJjz2tyIEXfe7tXpTatqaqpeiYiI/h5Azcl8X+shMFfGeb24iAg9WGI9BLaAuSOYY0kUt8tBaQ+JZASTiSGXA8JvwSB+i46OLkXTG7L+ZVG1iOIBXMPMJ0FT4u+JFiy1arRgSUOjtUJE1KZDrxXe8t0d8+66+UwyGmE0GuuCpUAgAH9NDbrPmrUIQG7n5CsyBBJMHldOXdBic6a9DpCZQYYiV/Z19dfo2PEqiz4y+CsTPAJR6JSckcEGjyvn30REic70l0Hy1sK83Jn25Iw3dMHA/Vu3fnmIiGBi8hVzBBJidRBmbHVnr544ceLKiy++uP8ll1xyUOaACHv27MH777+PcePGWX2+qt6yHPwXs5AoAAxCJcushwCDAHmrIOuzjJGRm3CkJrwnjuDz+ZJk2T8iOjr+odLS/UWRkbGJNTWlpZddmh7588pfdFWVBwqskXFJcW3a8KD//KfypZdmLhJFYZ7BYHU3s28afyO0YKlVowVLGhqtlURnxlQCdSlyZz8nknGYJJqv+E8PuyczJaV77qZC84e///QxQN42Hc7uUbJ77VXMLNtTUt9jRgAsVoI4hhlikTunUcE/myNtHglY4tmY+xUAJCUNMLHRlA2iTkSkdHVnTpQBU6E7p18XR1qyzMLDnvzs2+rbcjj+HemFf/6o2/91/fnnnYeuXbse0p4FUOqsAOCTTz5BaWmpOGnSJAHIE/LygB49ejAA5OXlUY8ePWQodUrHXavURER/RUXXoIgxm/LdQ87p3dsaFdVGd889/+WxY8fICfHx8Hq90Ol0/N33P9CwobcL+Ztdn4iSNMMc1WY1TqKYoUbrQJUxeJqZRx7n9Q4otVbtmNlHRE9CkWE44qGLvn378urVDUpUabQ8f+lgSSvw1jhtSUzOGEQCdSl0ZQ0BAFv31MjS4oIRn6xdS5+sXWuJb99zNIBZAO+wRrWfHxHV8WEAT0ssFBFwFRNbdrq9XZlXNJr5sDnS5jHRn4Ubc74KjRUUrPDZnekFHlfOv0JjSSnptwsyJiY5MxaDBBChfxdn6hPBqorniopW1rWDcbs/rZwyZcoz7dq2vb60tBQHDhw4TDWbmVFRUYG4uDjU1tbuANAR6IEeYQeve/Ro6BT2iTN16lQOnSgUBAGbN2+Oe/3118vYYAiKcm2V11dTkZDQIdrn86Fzp04HIiwRlSUlxdb4hE4JFRUlO2686ZZOQ268oUIUaKfIhipoW3B/VaYCeOV4LlTFPJ8HDlGnnwmlAD6zwYs0NFoYLVjSOC1p3/WatgYDXRcKlACASZQqS/J/Cj23OdKu5d3rtqlPh9pSUn+3OTNSiKErcuc4O56V8U+b0/SGzZH+YaE7Z1n9NWyO1PkAbShy5bx6NH8KNuYssDkz/hmqVbJ3Tx8BUZggWmO7JTkzwIry5X5J8j/3xIM3/hkdHQ2LxYKSkpKQeKVyD2En96KiohAIBDqc6HvVFJ599lm+8847cf/99x8yTkQHZs6cifvuuy/e6/X9ds455wYKC7eaAUGGLNUwQTKazIlVlSVtQCgs9GytJGa/TME/TBGxO6EVnR83RDQMwL8BREGpd3qSmT8iog1QJAP8UJSvF6tzdAAmMPM3qgDkE6qpNQDugnKs/ykoX5OtUE7tdYGi9h2Asl16KxR9pQ+gnGLTQ1HaXh/mVxSA85j5bvX5Ziin3pIB7AXwfwD+A0XrKZwHAayCIjD5KFRhWQBg5jIi8hHR2cy87rjfNA2NZkILljROO3r0GGA1GCJmkZ7qfhl3cmY+JRLtPXQmEwDYnBm3ApxKDCuIawvdOTcDwK4N2T8B+MmekrbI7kzN9LiW1xV+2xxp80BUVujKfbH++l0caTOJjrz15dmUMy/JkXF+gTt7OAD07TtSX1xdOJ0E/Yzo6GhER0fDZDI12DJFr9dDlmUYDIbQlhyhGTM0s2bN4oEDB6KsrKzBRrqDBg3CSy+9VDxmzBiCIk6pceqwArgcQAIUWYdP1bEpzLyGiGYA+JKZXySiTgB+JKJkKKre5zPzPiKaCEV9ez6AC9WxKVB0kwwAfgdwP5RgKhaAHUA5gJsApEAJxMLpDyA/7PkZAP7FzEVE9BOUQOpDAB/WvxlSGulmMfPa+t/3ANYBGKB+1tBoVRwud6yh0cqplCPmgoXvPes+Lw2NCTLHF7qzJtab2sbmzFgMcCqBx3ncy5MbsufZmDsURGZ7Sto7AGBLTltAxPpCV+74hubLRDHn94w6eq0GoSTJkTEXAFavnhsocOWMLXTn3hgZGYlQQ96GPkKvRUREwGq1Ys6cOVc09b05VqZOnTouLi4OO3bsQFlZGcrKylBeXo7y8vK651u2bEFcXBxmzpxZ3Vx+aDTKd8wsM/NeKCKhCep4KFhxQlXPVvWJKgB0AFDKzPvU8ScBeNXxJUS0AkrPNxuANwAUA8iFkgkKAsiB0rrlUwBP4vCauHgoGaQQxcxcpD4uAmAiomuJaEW9j/OhCGjeofrQHsAXYXZ2Q+k5p6HR6tAySxqtHdHnK0tkFi4EWB5174yhg9Iv2jZt6pjVNeXF50mCaFq7dnN86uXnt/N6Ky4IBiX999//GTvztfdusxjMZweCwdFbN36eHTLGwPt2R8ZCjzt7WPgino25Q23OtMW25PStIP7W41o+/DBPwliyZMlRt5cKXNkPJjkyXg89T0rOSK+q9V/+0cKJMJvNdc14G8rmSJIEnU4Hv98PvV6/qqlv1rHSvn376SHJhbKyMuj1BzusMHNdn7uoqCgEg0FLc/mh0Sh9AICI2kHJ8OxTx0MBjAtKRmiNmlmKhRJ0xBBRG2YuIaLZULbqdgD4typueTWUpsj/BvADM08mohsBPATgbQC7mfkKIroAirDlwDCf9kFpHhzisKxnY5klKG1doN5TAZSgLURs2P1paLQqtGBJo1XDxcUW2WLoKYq4as2fmy+85aZ038BLzquoqaoaD51IVkuU2WAynDPtyf+6amp8j2/eUtQrLiEqMH3ymHenPv9G0doN284Mt1foyvnC7ky/wd7j0hRP3teHnryRUULgIITGe7LZUjImQ0YjrT+YaPBg0ba2qr9EGC0Qy0TUL8mR/rtMcEOAMcKs32IwGGA2m2GxWBrsZxcKlnw+H3w+H+68887K4cOPGLsdN7GxsYiMjKxTIq+/NRLSqgptCc6ePTtqzJgxFc3ijEZDtCeir6H0exvFzFK9r9HTUNTCr4Wiuj1C7Ts3CkAWEUlQapZWQVHkziIiAUoG6lYAkQAWE1EQSgB2HwAPgA+IaCyU+qYn6/n0C4Bnm+Fe+0GpZdLQaHVowZJGa0b0W3VxkNHXZIoeXFpaXNuxY2f9geJdsU9OfcqyZs1a8X//+7D6zddfi9y8eXPnFd99Rz/9+F2gz7nnlIs6c5fePR0FLpfnjPpGCyP2j7RVJ8wFcHtozJ6c9iCIzR738uTE5NQ0e0raIpaNjxS6P90Vfi0xJ3ryc24HAFtKxvUkcwYIJgYFAFzYeUPlBxAgiaDNckBeWLhl+ZYkR8brhS6ldomI9FVVD9zTqVMnS0P97Ji5Th+KmWEwGAClRqUUzYDVakVERESj9VMhf/R6PYLBICwWiwNKWxSNU8N3zPxw+AAzJ4U9LgEwqP5FzJwDZTstnC9w6LYXoGRyLmhg3csac4iZq4hoFRH1ZuY1zNw+7LVGJTgasJOjtd19AAAgAElEQVQUekxEbQDomdnd1Os1NE4lWrCk0ZrRc4Dtv6/deLvdnsSxsQkGIkCSBcOjjz4qfbhkKaqrfZE//PijMHv2LHn37l0455xeAOmMAJnTrjjvz3eWLj/sKDKvXh2wOzIEe3L6U578nMcSkzMGEeHMQvfyOwGgKH957gAa8NV2J8+39bhikyzr1wnMQ0CQCHSBWgdFBNLJjDy/6J25L29Fld2Z/gJk2u5xZx9SFM5MeltK2kW1ejk/MfmKp/4z+DbX6pVf9TGbzXXBUii7JMsyJEmCJEnwer3o1avXc1DqTZoFs9mMiIgIGI3GOuXzQ31XslyiKKK2thYWi2Vzc/micVoxEcrJujtPkr37oGWVNFoxmiilRovTpctMTks7A//3fz1RVubF+PHZCAZ97sLCx89/7MnZH/bt1c2ReukFW0kUTADMz8+Yab3m2htizCYzv7lwQcRDD963iYgCOp1OBrPxmWeejRj30PiPI8xRb3c/5+pxm9cuO0wcEgDszvQFJMvPyUT3FbpzR3bunnGnKOJyZgShnHaTZZL/CUAvsDiLZLknINd68pc33C7FkTaSSHR6XFljD38tfRPAKwBpUlH+V3EARrhcrtFRUVEwGo3Q6ZT/WyRJgt/vR1VVFbp37/4alJNDn5+kt/oQiEi3cuXKQEJCQqNbgqHgze/3o6SkBH369DEwsyYyqdEiaKKUrRpNlFJDozmw2Sa+cc01/W/Pzh6qqlgDiYmxyMkZDmbZ0bv3SxXde3o+f+Lh2++TJJF0JOuff+W9+0Rzh+LPs39Y3++CXj/ExCUNIqK3mMkXCEisIzKNGz/eLwRpCxHtCATlQwQnu3VLj/Lr6TYwzgPQRxbEXwn4zObIeFcQIbHMvxfpCmdxXl5d/y6bM+0dINgbghgIlxeoD7HhC6baC8PHEpPT5hKBQFRkDEY+tnnz0v1EtLtNh15n9j73wjX/HTXMd9ddd9VtgzAz3G73/quvvnopgBVQTik1F/q8vLyNV1xxRYrZbIbBYGhwS7C2thayLMNisQDNuCWooaGh0VrRMksaLUL37o87Y2LsG++55x8QRQFE4X+kZciyjOzszVi3zoP16+8hmyP9ExBqiNjl2Zg7FQASk9PnSdA9uit/WXF9+4mJqW1gEW8UGONY5J8BCvtGl00MYYUAqQ+xEAwaMHnHupwdR/LXnpL2DCThWk9+drcjzUt0ZrxZ5Mq+Lalb6k2ySAOYWCpyLb87qUdmew7KT3vyc263OzNnMLhtkTvnPgADLNGdJ9aU7/zabE243Fu172MA2wB8A2AHN+MPqFroe5/L5ZoRExODUA1VeLAUDAbh9XpRVlaG5OTkZwFM1BrBarQUWmapVaNlljQ0TjaRkZ03xsTE4KOPCuBwxEEQCKGfNaVORkZFBSE+Ph5O5zPVBFrocWX/t54ZyWAWg3ZH2jNMQqIiQqkERRQhigT4QfKOona+obziYEsTGjxYTFxXMcfjXn4b9e2rr1/s3RDMJAhC8Fd7StoiCoqzCjZl/dnIRMHmSJ0PgYRCd+4doeGCvKw9SY4MuUv3K3uxyLGFruxb1Zc+SnRc8e+a8h0Px3Y6J77GvfyJBu02A8wsE9HHTqfT7nK57gkEAjCZTHXtV0LbbzU1NUhOTn4VwKdaoKShofF3RAuWNFqE2NhYREdHw2AwYNcupQSGiBBKpMiyjMjISIiiCL+/o6WmSnje5kx7BCScB0YNyehLoBrZ538VEEAsfOJxf760/jo2R9q88EAJABLXVcxpa7X/Fwgr9nZmPO5xZU9pzF8CEgpcX9wMAPaUtEX2lNSdno3LDylItTtSFxKJ/Uhfc4Fn3Q+HbVXJhA9IDCwtcn9xiDimAIGZuTYpOUMgmiQwT2ruxrh1MPN2IprjdDprAETn5eXdFQqW9Ho9Lrroord27dpVAkW5+5dT5ZeGhoZGa0ILljROORdcMHdGx45JiIiIOOzofIjQKSxBEOD3+7GjVH4SJFjA+JEl4QMI8qOF7sMyTQ3A1KPHYENe3pJaQAmeJMKbq1fPrStS9rizh9md6Qua6r9nY+5Qe0raokRn+ltFrpxb7c60qcxoC4E8YEgNBUoqNwBCRfjWWpIjbQAzRMVTBO0pq84G0HDWqplgZhcRTQXQr0ePHlugiAMSlJYXGwD8xsyHbXVqnFyISASQDSACwFXMfNj3ERFdDKCMmdcR0X8A/MrMu+rPOwEfXgAwk5kLw8YcAF5j5gFhY0kA3mfm/k20OwLAmyd6OOD334HDu6ScGFolikZT0IIljVOOxWLpYrFY0FBRcQhZVuqWBEGA1+tFoSv75tBrid2v/I8sk6FJizEJ1YHSrgBcNkfaPAA/73Tlrqw/TQD8dmf6Ao8r57DtOHty+lMkHtoQVlX8zrYnp29mkr8pdH8xAgBsjQRdNmfGG2DoRRbm2JMz3vDkKw13GTyKSGkoSoRykqkXTnGwBCjaOQC+JqJvcLANktycNVMah9EBQDwz9znCnNsBvA+lf9q9UBrknrRgiZkPO8l5kngUwFtQGvZqaJx2aL3hNE45UVG+R8xmM6xW6yEfkZGRiIyMPGQsIiICZrMZSY6Mm+oMkJTO1f6Hj7BEHQK4iknXx+ZInU9MOwrduQsbmrfdlXP3EYx08LgiDzkFZ3OkzgfzLtbx9yQIptA4A3TmWYMSw+fanRlXAkyF7uxbt+VnvUEU3h6CvIX5F7wHAMzyfmapc1Puq7lgBUn90AKlU8s8AN2IaC4RdSaiZUT0JRH9QUSDiKgPgDQAzxHRUADnAHiLiEYQ0XRAyU4R0XoiMoaMElEUES0hoi9UW3er4/2I6Bci+pWIPiYis9rDzUFEHYjoGyL6FkBjdXQJRPSZauNx1WYiEeUQ0bfq50QiugNKH7j3Vf9eJ6LlRLRabeirodHq0TJLGqec//1vRNmoUX/CbDbDZDJBp9MdklkKtdgIBoOQZRl6vR4MusrmzMyUARYhnM8mUxsAB462lgTxWwg8jUBfeNw59ds2HIpMi+2OtEUed+7Q0FBS0gATzGYwK73gbI601wHI0AnTCjfkbAWApB6pt9lT0hbpoJtKILk26OsNpaGocj/AfwrDM1aM6iRnxuwCV/YY5X6VGiWRxQIZct+mvIetFIuvqvwfEtEgAvsYvBcggwAOymCDSKJBkjhIwG4WOEAkRAuCsMFojPiViKrQQI+xvxmjoGxtjSSiywA8z8wriOgfACYz8+VElKvOySWi26BklnYC+IOIHoYSTH3LzP4wu2eq13xMRB2hNMmdAyU4u0Hdhh0FpSlviAcAvMfM84noegAN/TNhBXALlB5zPxDRZ1AySLOZOYeILgUwjZmHqMHUDQASAfzCzMOJyASlX93jJ+G909BoVrRgSaMlEHy+4gpBSIgyGo11W3HhBAIB+P1+EBEWv3v5wmBt9R3MLNuSM94gFn4hPU+xOTMlAGBw8R3unPsmMR9WGC0I0o0yYCx0LR9zNKc8+dnf2J3pN3c5K7XX9g3L1wIAzJb5YHm/zZn2CBhdiLHJk5/7XPh1BXnL3wTwpj0lbRFY6MSCeAmAzwBl+00QacYh8/Oz7w012GXwwRs3iOu4tvE2E60cwev1xsmE8zt06HTXrp1FnjffXGj69bdV0W+++fqes87q1WXfvv3Cwjfn777s8su3EIgjIqIvrq4uv7+mpmZDRERENbRgKZzdACaoWRkG0Gi/QmauJKLvAKQCuA2H93LbA2AsEV0DpSdcyFY7ZnapNl4FEP5PSw8oDXUB4Cc0HCytZeZy9brfAHQH0BPAo0T0EJS6t/qnJ0sAnEdEA1VfjOr1n0MJvtYz8z2N3auGRkuhbcNptAQHFiy47MXy8ipUVlbC6/WitrYWgUAAgUAAPp8PNTU1KC8vx/ffz/guqk1SYmL31LcBgIFgwabPby1yZ91Q6MoaYgzG3k1A7ILk9MU2Z+Y7dmfGEpsj414AsDtT32SZqonxVVMdK3RH3ilLwr2h5ww5HsRmMJ9R6M4dUT9QCsezMXcoUbBSkPg6AEh0pr0NRrBgQ5ar/lxmMtsdad8Royo0Fijbv5mIT8ufSWbWA0G7IOCc4uI9FbIUTLRYzAlz5rxKPm9NoizLdPlll8JoNMcTJAcg94iMjJSZ5S6yLLeB9ruoPlMAvMXMtwD4Fgc1bGSE1ZSFPZ4PYDiAtsy8rp6tcQBWMvPNAJaG2dpFRN0AgIgeUgvGQ7hxsGfceY346CQiKxHpoDTBzVOve0gtBh8J4MN6vg6DUqA+BMDzACxERMx8JTMP0AIljdaKllnSOOUwc4CIFr7wQnLkrbcuu7l9e2e82WyGKIp122/V1dVYsODSLeXlu5YBeDkiutPrtuS0nwDaEm5r8+bFoe7pAIBEx5WpBL7LlpwxhqBjFvCnwLLYdN+WSPaUNNHuSHuGgTYky2d6Nn2R3tTr5Zjom+hA9Ua7I2MNEa0pdGePbGieJz9riC05favAUl0tVFHRSq89OUNqaH5r4axbz3Ia4gwbdXodWGBUl1ffu3HOxpfKy8vNRh3Z33nv/TNvuPGGsm7JPSK7dEmS+/btXZ2Q0Db4xMTHrHv37RNZlgI//vizbsCAgbXBYBCiaIgXBCmOmUUiCh7Vgb8PSwHMJqI9ULZ049XxXwFMI6LtAH6GUrN0BTP/SkRnAnilAVvLAMwhoiFQtq6Dak3TSAALiEiGksl6AUrROKBsjX1ARDcA2N6IjyUAPgCQAOADZt5IROPUtUwAzGH2foBy0u+/UGqXLgJQDWAzgI5QthI1NFotmoK3RotBRO2h1DF0jY21x5zb5/qMkgP5a9es+XQTgOrYdo6+lpjO3+1wfzkRAGzJGYtJkGuZuabQvXx0Y3Y7O9JeEcDGQvfy4YmOKx8QmO9mAb+CmZhQ0bVdzegV9bSXwrE5094giS+EKKxmYFWhK+eFY7mvxJSMhZDRTxDkXz0bc4c1uk5y+laCboQnf9nXobHwk3KtjX4P9eNvp3yLQDBQl5vQkQ63vXQbyg6UvPzJ40tyBIh9mFkvQDRJJFkAChLBRxIMssAGAHXl7QJBYhK262ThS4PVmg/tpNRxo6qx/wQglZkrWtqfk4kqOzBCeWbrA3hOqn3tT+BJ4y+t4K0FSxotivpLPgZAYlz7XjMqS7bNrK2t3AJgJzPX2JwZbxe6sm+hAQN0nXebXy3Kzxlhd6S+B1BteCF2CJsjYzIgdyp05w4/OJY2r9CdO0J9/RUQxQAAGIJAvKrAlT0TAOyO9KcZHA8ilwgYZBnDPfk5R2xv0hD2lNR8K0X3rOKKiwHcEhRME3Zu+F9R/XmJjrSFIoRggTu7ztckR8br4c9bC/+a9C+eMHgCdKIOoIO1LSwrxfhTPpyC0r2lWX+8/MeVLezq3w4i6gLgEwBzmXlOS/vTnBD1ZeDktjvR/gSeNP7SwZK2DafRorBSlF0CoMSWnLbXQInf7dq1rKbudRmbbM7MdzrL5iCDFgGAx738xsTkzDS7I20REdYWuHJnAkCiI/U2AnUID5TqEy5kmdQj8xxJwiOdHRk5ItOZTIIIWRpauCnnh3jHvyMjUHtHuKBlU0h0pr8Jogr1mq969Bj8fZVcMd/mTHMXunKfqZuXnH4eEYsAgh26pSfs3pyz/whmDeW+8kS9LKcwU7woCgRJZhZ01RwMBlhgqyiKekmSCQATiSWCIK83maKLcHiB7TFz6eRLOT4+Hvn78pVWKAQQCAwGGJBkCe3btYfJYMo80bU0jh1m3g5FRkBDQ6OZ0IoqNVoPxCyKpkP+zyvKz54CsEyguB352d8cHM/K9bhzh8qM3rbkKxZ3cqZdQKALQhmkQ8wSN3iSyKcLbCbI5QLLuzz5Wd3ANB2iOMbmzHwngvzrQaiplCvfbKr7tpS0V8AQiIW6Atu8vCW1no25QwF22lNS3wuNC4TxIGGDKAjPGgSadgSzAlDdRi/jMkAYkfV59nhZ5tGfZeU8VFxcfD8I9y5e/O6jRoPxv8s++/yh5cu/eGD//j2Dg0HxDKXo+sT+27to/EVnh1rT5O3Lw/o967G5eDM2FW/C5uLNyNuXh7y9ebBGWBETE4NLJ1/aqmuu/goQ0QtEZKs35iCiFUe5brT6OU3d2jpRPxYSkfU4r7UQ0Z9ElKY+zyCiI/Zn1NBoSbRgSaM1waK4/7CkODN/y5AdDV1Q6M69hQQxoJNokeDzNSgPwCB9+15XtA0fS3SkLTT4hVmCLH9TmK80vC3M/3xOoSvrukJX1hCW6XeA8gVwf7szfXWiM+MDmyPt6sYcTzpnUAxYMBW5c4YCrB8wYMAhWdtC1/JbmUmyOUOn+mRfpOB9ZevGzzcfKlJ5GFRby3HMcm9moXdGZqZNICHlqqsyz9i/Z++5I+++55+33HJzp/XrNzozM9Lt7yx+r0u7hLYSc8BHRCe8wWBtY/0lOjq6ThzUYrFANIh1HyaTCRaLBVarFVFRUYiJidF+pzQzzDw2vB3JMTBBvT6XmeediA+q9tLvqvL78fAKwqQimDkbwHVEFH0ifmloNBfaLzaNVgOBWK83H/4HnqkfCcIvNmfG2w1cBmYO6A3ypWwyzbU50l463C4KTD76PwCwOdLutTnS5olE5YXu3BEFm754/7D5BCIBFR539vUe9/IzibFGZHoNEG+0OTPfsTkz30l0ZC7q0Ds9IXSN7K99/ox2NSMBQAbEgj2mC+vbLXTl3kwkBO0paYsIJOflragCAJmFDUmOzEWqpcO0otjPLAhisLKiHKJIEiCyXmeUR9w9Sndu73OYiGR3votFvc4XCAZ9zMJ6nY4LAPhxgtpFce3izCEldYvF0uhHREREXcCkcewQkZ6I3iain1VF7euJSKc+TyWidkSUpypiH1Flm4iuVRW0V6gf8UT0GIA2RPQqEQ0jUrKZRPQAEa0iopVE9Kw6NomIFqkK3BuJKLUBl++B0nYF6hovENFXRPQbEdmJ6Myw9UMfI9T546Cc5Ftbz2Y2gMPqEDU0WgNazZJGa+KwP+yDabCI7kyFrpybbc6Mt4kmCSHFa0Ap3maRH96yYXkJgKF2R9p79VW4ZaLtIJxvc6TNI8DgcTd+Qg0AEpPT58qg1YDS/sPuzAiApWGF7twb6+Y4M6bqvMJLNmemJAD9GLRn506zCCAoAOslga4AsKK+bc/GnNs6paReryPxqfa9rmm7Z+3H+2rLE9/V622z3O7iWwvcuQBwp8u1H5dfPh9ExXEbN07eTxL/HBkVGW00Gtt7q2t0sizI33/3LRiSYDJFBDMzM2qIIX388ZKy2tpgvtnE5VC0bU4Is9FcFxDpdDqlZqkeoabHRASz2XyiS/5dGQmgmJlvIaJIAH8A+BrATQA+h3K0fxwzF4UJRzamst0dQCYz1xDRXCgn5J4ionuYeRQRDQMAIuoJYDCAfwAIAviIiEIF+n5mTieiy9V1locWJSIzABszh9fZ/cbMY4noKQA3MvM0AAPq36Sq6t1NVSn/Z72XQ/3uZh/TO6ehcQrQgiWN1gMxGwxRhwRMK7tVLRRJCDWnfSvR8dsiKC0WYHOkzZP19PKO9bklofked+6Niclpw+2OtEUysKHInTsdEtIh0Fk+wdt3n5rNOaIbBF2hK7tOr8bjyh5lr9cgt8iVPQEA7Cnpi2QWfiWG4NcJC2zOTBCTTThCb3SdJFhB8ipjoGZ6UtKDZ2Zmpv5j1Kj+qKmRsGbNPgBKAPLpp7fj2We/O3DuubMWblo7bkKV5N/p83ljCSSAWAYEEEisqipjvaCrCjJkyednnUTbyRLnxUlQxPbt8I03djBOj4iIgNFoPKzpcag1TSAQADMjIiriRJf8u+IEFPFUVZF7I4Cuqn7Sj1AEInPrXdOYyvY+AItIaSHjAHBY42gVB5TWIwEAIKIfVJsAsEb9XATAVO+6WADF9cbC57dXNZ9erzfnXSgBlF2tr3IAOJeI9jDzn1ACwrhGfNXQaFG0YEmjVWE2bzv0DzyhpsCd9S0AFLqyv7Q5M4baHZkTGJINjKwd63PqqxWjKD/3dQCv27tnbLF3Tx8hEo2RIV/XlECpb8erLIhuoChaxjt2R8ZCjzt7WGioc/eMWwWRgoWurEO2DuxnZfwfgvy8zZn5jnILrIese8iT/9l2AGDwQJLFxVbxvC8j2nb2p6R0xpo1ByCKAg7WYzMkSUafPl1QUFA5zJY8YVRR0czvjub/yWbxw4tXPrniSRgMBphMJuj1ehARiAjMXJdVqq2tRW1tLWZkzhg9naefajf/CrgAXATgEzWz1BPAdiLqD+AsAN9DyfCEt84JqWyvhaqyrdb8TAYQKgD/Ege/qep/X7sBPKAqcEsALgbwFoBeOHKgfQBAZL2xQ+Yz8xY0kFmC0o8Oqq8LofSs+1MdioUS6DVKnz7A6pOrHKCh0SS0miWNVgOD+PczzqjbOursyLiXcOhJtkJX9s3MPJRAmwrzcz9tyI4tOe1+myNtHkT5CxY5h0m+nsBNSnnsjQ4uAPiwX8ee/JyvQSzbnVf2CY0JIi4pdB0uIOnZkP2RTPxNoStrSKErawhAeUzBp2zOzHdsjox3iYRzPZuysozGjv6YmBj8+ON+rF5dgnXryrF+fRnWry/DunXlWLXqAFauLEZMTAzi47vV1F/nFLF7xnUz5vt8PkiSBGYGEUEQhLp+frIsw+/3o8ZbAyjZAY1jZx6AODWLtAJKwOMH8AaA2wGMB3ALEYU3Wn4cwFVqliZ0+KACSpbpDyiq2V4oCtkAsJGIFocuZub1AJao838DUADgf0dzVG3Su4eI2h5t7jHSD8rWo4ZGq0PLLGm0HmQwlqbU/Yeq1wtnbVu3bGxtbclZXq/oi4qK8v4r447nbJ0TSvV6/T+Yazp5vV6YAcBsxuuLshLe/SBnfJekDsKK3AUjiagGSi0G7N3TN9mS098tzM+56UguCECNx5XbUMsIdGnnHbF9r3kegNttzow3BIEea9zSwZNoHlfW5PBX7I70/ERH6tvJne9ATEwMjEYjvF4BPp98UOyRGcwiIiMjodfrIUktdiLfU7G/YtX0wdMj7nv3vpusViv0ej0EQQAzIxgIosZbg7LSsuCUK6a8AOCbo1rUOAxmrkXDxc09wh73Uj8PCBtrSNtqcCNrDGxgbCaAmfWGJ4W97kbDGaJXAQwBMEvtAxea/1pDazfiz7B6QxloxHcNjZZGC5Y0Wg/EDExmYBIuSLzOHJVo1fuqK9KZhVsEQY74bc26rjOeuj+Y3L2DO8+1vZ/bvSPHbm/r8Qsi/7k6v3fPlDM5+6MX17KA0pqaykEREVFfQEnryyzwCiKOqV/8fcjyPQYbbEeoNVqxYkXQ5szQ2RxpeQTKLcjL3nOst9j57PTOAvi3C84+8+5AIPLmxmqBgIOF0zqdDsFgEPWL208FzCwR0aI92/cEHh84rdttMwa3s51jswmkBEsyZEy4cMLHlsj27QC8yMxlp9I/jRbjfSh96awnIB9QBxFlAvjor9aqReOvgxYsabQmmNX+O7siql7r2yX5DSLq+vwLL/S86667Ol4+4CKd3WaTV678sf2W/I1CXHx8u86dEpy1sh+PPzpe+P/27jzOqrr+4/j7c+7sMww4CAI6M5AhM25loqi/VCoVGHPLtM0WU9NMy6zULHNPLbPMXHAXxRZt0WLRLNFyDTU0ZXBlZgQEiW32mXvP5/fHuTMMA3wB2UZ6PR8PHs7c7Zw7DsNrvuec7/fPD96/VMrs94PzL2r58aWXPu7eFJuVZEd43OpmP3R8ZdW4Ayqrxt9l0ltza6df1HPjlfGKGzLu/wjtoJm/rVgH1M2Z9p11vJU1RlfU7heYafCz/3nz5wd8uFCFhYXrjKVUKqX29nZJF22VhRncvaOiqmbEogVPfH7iGRPHSBohqUjJWm4LJD1TNnTP45pXLHh7a+wftrzs39MvbsLXm7KpXgvYHIgl9Bnu6j7W5LK230++cFZbmz7Q3p4eWF9fHyW3y9zj3CgVRfLY4jj2kuKSDslyYo8zkeV1Xnrphf/qjNOv5WtAo7pPPE1KpK72oSclPVlRNW5Sxajxk+vnTP9Cj21GDbXTV7nqradho8ZunxMV7WDmE4dXj7977uzpq/1jUb7bcWWebjk8MhtTUVVzt7tWOdco8vgTsVKPdDY3/66wsPDkkpISFRQUBGOpvb1d7e3tct/4q9vei9GjR+dKQ7xt+fLXJb1uyY6mJMXZ5WpkZv+urK65rC57lSAAbEuIJWxVxx9/fP8xY8YsKysr0+RbfqCJEyee2tTUpFSO7pG2y1z18+s/NXRYeVPVLiPrjjvuU4N+cc01S2oOP+oDd91114rzzjs3M3vOa3knffXLtR//+MGD3HMW3XLbLf1OOunEl3NzC942sx7rormPqBo/6q3a6XMkqb72oS9VVo2/s7Jq/F0WZ67M9Zw3lfLgBQ85UdGV0Q6DvxG/s2h3d326fFTNHS7Flj0vKrsdMylfZs3m0fT6OX+Z3PM1KqvG39VQO+00MyvMyytUXl6eCgsLlZOTs0owxXEsd1dnZ2ey7ZxcmVmBu7dtmq/8+lvUNPgH9bV/uajr8+yoQrrnY9zdK3c9PBo9enTuzJkzO7f0PgLA5mTOksvYSq688ko/+eST1zh3j7vr1ltv1bOzFt5/z20X/iknL7ffTTf/4Utm5scdNfYP/bfr3xbHGb32WkPhn6Y8ftR3v/X5SWaWcvfOKE69Vtiv3ywlVwa5JJVXHXaDeTSjfs703/Xej8qq8Xe562Az+1n9kNYbKxbklbulTnTXkOTiL4vNfKyZNWZizTKpU1EcR7KhdbOnHbOm91ZePeFGi/Vi/Zxp3avAm5lVjBp3Z3t+/L28Vr9i0YIX9j37zJm7DxgwQHl5ecrNXXnhn7srnU6rvb1dy5cv1xVXjPiJpAuyJwJvMcOrD/9c7D64vnbqtevzeEaXsDmNHj3aZw5+EAgAACAASURBVDJ3QF+1UetQ9nWMLGGruOqqqxqHDBmixYsXK5VKrRJLcRwrjmMNHjxYB+ydc9Auex1fb7J+snhi/eyH7+r5Oh/6UH8d+blz9r71nin7NLwyZa2Lg0ZR9JRi7SfpdzvuOmGvVEYHx9KHJItNUZsiL7NM/I2dFhTuGpuKFevF2PXHea9Ne6G86vCfuenJ+tlTT+z5mpXVE+6o3HX8IXWvTH9kDZt8WJEdK6k7lsp3GX+axTowrz26LN2Ue2Z703+PefLJa79xwAHf2r/nid5dX4OOjg41NTXp0Ucve1TSC1s6lCTJpV3qa6devO5HJkxqLq+ecHDD7GlbfE4oANhciCVsFYMHDy4ZMGCAFi5cqKKiotVGlpqbm7XddttJ0uBMJs6f/9pf1xpCda9M/VzXunE2dmxO5TsDdoqt4ySPNcjdUmYexxkVW66N2WmXmtLI1Olx9E5Ovn5S9/I+c4477hV75qXG3Eg5E2WZ082sbu6r066WpMrdaj5upgG9Q0mS6mZPOzE7s/dqsZRqzX80Luz4lCSZHZ8qH7Vionk01iP7eX3tlOuS2+0Pjz32iyHPPHNP2emnPzkqPz9XOTnJX8l0Oq22tk79/Oe73CjpdUl/fM9f7Peooqrm4vyM/3hDnjN39tQrhlfXXCqJWHofMLODJC1z9xezM2kP2cjXm6tkZu6zJP3d3Z9dy+N+I+lLPX8BMLPxkj67hikFgK2OWMIW96tf/WruTjvtpJKSEknKXum1UtfEh133f/+bR35J0hld91fsXDNaOb6PK9rP5WmT0pHFwytG1by2oxfPSFtHYeT2XGyaNm/OlCclabfddstrzOx0ff2rD53Se3+GV9fcam7/eOvVKc9IeqaiatzdFaPG31M/Z/oJHusLa5p4cuXOWlRZPeGmutnTTut5c13dA8vLq8e3VlSNm1Q+ytrijN2aSsW5dbXTuxf6dffW4VWHqmnFwreuuWaX+yX1XHHdJC2V9Kikx919lXOENrfsSd167bW/tK/70avyVM704dWHnzd39pQrN8e+YZP6qpJpAFabCX9jZNeGC93/2U25PWBzI5awxQ0ZMqSytLS0e3HWrkNPPcVxrHQ6rSiKNHjwoH7lo2omSpK7UpajtMvejaPMtQuKF7/k2ROKK6pr7ra0P9rw6tR7e7/eyy+/3FFRNX6NJ+i5PKqrnXZH1+f1tQ99sbJq/BWVo2peM3nwEFRd7dSv9F43TpIqqg670d0OiqSr5tZOv0OSKqvGf73342LL3XXx/BcnmFmupFKtXEaiUdLyLR1JXd5t3uH8+tq/XLjuR66u7j8PPjG8umb8pt4nvHfZ76+bJI1UsnLDD5V8j41Xsj7bK5LyzexeJUul/FfSp5VMEXGbVq7Z9k13f8nM6pQslzLb3c9aw/buVBJhQyQdpeR7e3tJl7j773uMQI2QdLuk5uyfpZv8zQObALGELa6goGCVlex7XzIvJbHUNWt1UVGxGuZMPXVdr9tQO+1L5VUTJilZsHO9DBt2RFFu6epXwblF5YpshsXxoZVV446oq33oM4GXebpy15o7616Z+pWKURN+IvMBLqs30xNzax+6Q5KG7zZ2iJS/ytItFdUTJsnsOUnKLmb63+yfrWrEqCOPkvTGRr1Ia8vllbse/uO6V6acv2n2ChvpZEmL3f0kMxuoZLRyNzObrmR9tnozK5F0vrvPzS6hspekYyX9zd1vNLORku6Q9FFJ5ZI+4u7r8/1aIulQSYMkPWtmPZcpulTSj9z9r2Z2rpIFhYE+h1jCFldSUqKioiIVFxevdsm8tHIl+3Q6GVTJy8tTRfX42xXbiw2vTr+2a+LK3tzdy6sPb6qsPvw3dbOnrNcwf27/9O1StMp6VJW7Hv49l3c2vDL1FEmqqDpsUmXV+Lvi5n6nNTTc19r7NepmT7u5cpcJp1RWT/iTRfE7dbMf+pokVVTX3Na9b5mCH7nppa7Ph1WP/0gkS7/9yvpdZbaluHvhMUcfdPRlF56yuKVpWZUie8fM8t2Vdo9zI1m+5GlTtMhdnR55v8j1VoFy/q3i4kWSYkmaO3dG2/Dqmhwzs7X9/+oht6NpyahOSx2rjG9vUfSGKY48ijyOYzOzAkkeRTZPsrnumUEWW0GUGz2Xn9/vDSWTYyJsD0kHmtmY7Oc52WjqaYm7z81+/I6SUaU9JH3czLp+Wdgu+9/F6xlKkvRYdj6uhWa2VEk0ddlNybp0UrJGHbGEPolYwhY3adKkb5177rnXFhUVrXH26p4r2Xd2duraa6/5rMX+bOypM8pHjbulomp8dskPL/T2nFsa3pryeNdzG2ZP+XpFdc3dH9x9XPnr/3moYT12p6Vu9tRbuj4ZNvqIohz3qobZ07pP6K6vffhLO3/oU4PTarypctT4dN2c6d3nMJVXT/isuX/cI72SUtQ5d/ZDPc5dcrPRo3N95sxON5U21O53laRc9/aR3zv9c9edeuJRf2tpWXG5u7aXvEhx/KrkHWZRkZvFFmuR5ViTZ1SWijJv5hVGT5j1X6rwivAbww478vSrfn3HxXMOOujj33/8sb/Pf+utN/XJT35q5Jw5L70xdOiIDzQ2NuZMnz6lYcy++9SncqLMR/Yas9/z/3r27tZU3FIoLZPUPQ/U3NlTz6ms/uRlSg75hBR1SHvE6fTYVE5qJ4vUOmzYiD2+993vvJNfkKuJE28Z/NJL/37L3V/1jL8eRTY8v7i0vbWpqcnzvMHMiKV1q5X0trv/2MwKJf1AySGvWCsXVF/T91WtpHvc/d7swrknZ2/fkGV39pYkM9tByeG4Rb1ef39J0yXts64Xysyfr+UXr/fFmd36X/iejigD3YglbHF33HHHG+ecc45ycnKUn5+v3Nzc1WKpo6Oje2LGZ56b95XigTse2PDKtDN6vk75Lp/c1/I7T6moGvdFyVymHLm/q0zRRR1qv1zSl1bZsCmn5/pqO1YfsX+Or3oILmrsvK4pnTm79z6/MesPiyR9uWLU+HsqqsZN8k7/ieVG3zS51dcmJ41XVI+/p3LXmp/XvTL128kbsYKK5ortJS0wV8b9onj+fBX84vrrrj1k7EcqY7cv7zi0vHz+/LcWjxq1R9mc2v+8PmRoxQfz8vKjRYvmv9Xa2rrQYrWV9i/7WPPy/17S3u4vFRRomTZTLN18858Ld955WElhYd4eD/zp/syKFSs+cNutd0Uvvvjv9va2jspUKpVz2203e35eztAoUunMmc/njBmzr3mknZXODHL3VHZ27+79i+Qdw6sn7Dd39rSn17LZVFvb8kGRUh/+4KhdDmpubo6OOurITH393Ba5D5gydWruLiNHRZl0ujyKUoWxMh/41hlnV9xy621vKhXv1NTUVNivX7+WzfU12YZMlHSLmT2mJFhucPfYzJ6RdKWZvbWW510u6TYz+1r2eRe9h20PMbO/KbmA4fTseoNd950u6bdm9j1J76pHbAN9CbGEreHl6urq619//fVvSMk5TF1zLXWNKrW1tWn58uXaddddfyXpquLtyi+sqBp/cyR7em7ttNslqeHVvzyrlUP4Kt9tXJll7GKlWr/vcfR/5bvUvBlF8XRPZe6tf/mv/5SUV7bvMyVKJqtUjqfPLC4qeaSpacnukacO+sOfHz31Z5ef2XZEzf+d09K0LImoSG5u7ZK3xfK82TN//Z9rfvWbMxVp/AdG7PSTLxx/xK1Krlzz+tnTT1jlZO/Il7i1H6nkHypJ0tKlSucW5EWjRlYsOPTQwz7S2NhokeWUPfvMUx5ZNHzB/Pkd7q5MOh6Wk1JJHMc69NBDOmJ5heJooKR52rDf6tdoj1P2uCVVkDo6juOmjraOPWtvq226ZfLvLvjXY3c8+ObcN8/acciwlkGDh/VvbGyMrrjikmWf//xX+/3r2adWlJVtp1mzXky5zA844P/a991ndGNuXtF/W1oam80WuDRslWh5a/aUS7JTCTy9x1f3ODgqjH6rSJ5ekT7q5btefnbu3Lm5Awf2H5abm1M+b15DfXtbU3FbW3tuXq4yAwfuOLC4pNjr6t5c0d7WpvyCqCAy5Z1++tc6XvzPrJyRI4bnmKWCs64j4e7t6v3LQ3L7RK38/hzS4/aeh7GPXsPz1jjFgLsPz374FUkys68oOQx33loet0DSQet+B8DWxQze2OKyow/HSzqwtrb2G4WFybIfURR1T8bY2tqqqqqqG5TM13Nf13kvFdXj7pZba+zxH9+e8/C00HYqd625U4rvUewnuKwjlu0fefxfMz3viq82y720bvbUU5ubVxz2xhsNPz/33O/s/PnPHr/sm2edXdZQ/9bC73//B/3+8pephSN3+WB66pQH3x77sUNHPPHEk/buu++8NWjQkBGTfv3AS/fcO/Xthx686ZPKBkzlqAlXy3zH+jkPfeF753x/9nYDSnZpbGxp/PPDL/31pZkPHLv7Psfd+eOLT3tpj1HlRw8ZOqh66LARZU3NzdbasnzBrFkvFj7wwIMFX/3qV1qHD/9A3Ny4zKNUpE8d+9l+f/jD7/6WSkW/ys8vfUxadb25DTHmvDF+33fu0w4Dduiebzc/J18H//Bgdbw4+JoZv7ltRkZ+oLkVuHmB3ItiV2ymVnPPdbP8lYsEp9zMmyPpuY60/b20tPQt9VoGRZI+/M295p16+NeGnXzIyYqTpeSUspR+9OsfacbzM/78t8sfusrdDzPXdrGUozjup8hkHrfFiixyL/DII1Mq+9Mq9ji2hZar6SUF/Z8ys+b1OC8KW0E2lqp6x9J7tdewYT7j1HVe67EaDsNtEdv0DN7EErYKM4uU/EZ5uKR+11133cG77757+VNPPfXa+eef/7SSy5qnKrlqJ171uRdF5bs8eYssymQsunre7KmvrmkbO+36yZrI48/Xz556giQNrxp3Yka2Z0r+tsv2kWzvfqXF/zhi/IEjjj364H1GfbAi97777s+96OJLrbm5Waeddqr/9a+P2J577KnrrvuFH3DAgRr7sbF+5BE1mWt/eb0mT57025NP/3HZcy/ULlmWWf71RS/PaJKkq6/+qZ955jdXmRIhiiJdc801uuePT9331CN3/NRdNWaqkGTm1uKxt1tkxW6eL2nVM0mkWLJ5ZppWUNDveTPb4LmPJOngCw72K754hVJRShatfkL9hb+7UE1Lm77x5NVP3i7NVUNDylKplA0dOtST87QbbMGCVPcTh0nyoRk3K09LF2Wki1Yb8Trg/AP8gs9coO2Kt1NkUfeP064lbZ6c86Tun3F//OQ37i3W8OGSpPnz50eSNGzYsLihocFSqdRqP4SHDRsWSzPS0tiMOAT3P4NY6tOIJWBzyI4w5UuqlDRcUrGSSKqXVLc+i8ZWjBp/m0yZjBV+f97sP652dU5F9YTJJqv9dE3100OG7PBwJhMrJydHT/1rzrMzZ83/dl3tlGe/fd7Vj5z4+QnlL744c6e//Hlq6tEZM6KXX57Vfukll+dFqVTmsssubS/Iz2m96667i+68c1LeY4/PmHfUkUcP+s1vJ/82x/36/zv0675oydKz2jvS7T/49lGnHHHEEd0nrXedm9F1eHHq1Kmqr68/+Morr3x8tTezGR180cE+qGyQaj5Sk4SSSSaTy6VYysQZPfLSI1q2bJmmnz99k/zQ2+87+/1hyI5Djtlr571Uvn35al+P2GM9/+bzaljYoEVvLTr86eufnroptvu/ZE0jN2uaHXsTb/MMd//VhowaZa9oPMHdb80+b4m7P7ih2yaW+jRiCejLdqw+ZGDKU1fKLFU/e/pXe9//05/+1M8666zueZskycw0adIkPfr4Cy+UlO7w2LVXf+ttdxtmHtujT8z6aGtLW1lLU2tZyYCS+ZErdrMoNy+n/cO77fzPsrLSdzNueYozc6Iof0ZRUdECSfH111/vZWVlqqqqWusVfm+88YaWLFmiU089dYv+YDnm58d4WVmZ8vLylJebp+K84u770nFaTW1NSqfTamxs1DsL3ml79JJHCzd2mzVX1vj222+voqIipVIplRaWynr8PF3WskyZTEbNzc1avHjxJou0/yWb+jDXem7zHXcfsoGxNFzJfE77bcy2iaU+bZv++8sJ3njfmzf7kf9KOqWiqmZ0RdX4m+WWUz9n2lcl6YYbbvBjjz1W8+bN6x7Z6PoF4ROf+ITS6fReCxcu/Hdh4YBvDxn+ic+mcqOPxfJ/L3jtkdXWohv2wUNPMbOPSl4gSWbRsJRUPKh0+B2f+9wuZUOHDlVpaamWLFmiKIqUm7tyDsp0Oq1MJqN+/frJ3XXNNdfMPPvss0dvia/PuEvGPbRjxY4qKSlZeW6YrTxiZm4qThUrjmOlUil1dHQUbIrtDhw4UP379++eHiK2eJWRpZKSEmUyGeXm5iqON/qc9f9l+2evNuu6Wu16JbNjD1Ey+3aukkOV33T3WdnZtXeWVCDpakkvKJkeYF9JMrPfZm+/Tck5g3tmn3+UkmWHyszsBiUXV+xnZg8rmTvpRne/2cwOVnIVXUbJ5KanKpmqYFcz+5GSA8zvSLpZ0i8l7SspT9KF7t5zwkpl9+drkr4mSeX9+/e+G9giiCVsM+prp86UNLN8VM3RFVXjbz7vm0eeMnDgQC1durT7arsuXefpDBo0SLm5uScO+cAhbXn5ufn1PeZQ6m3+63+9RVL3nEyV1eOHuvsPFzXNvaFfyeiTBgwY0B0kPQ85SVJubq7cvXsJl3Q6vffm+Sqsrl9pv71KSkpUXFzcPU1D71nTu0a+oihSW1ubKnetOWdjtjmgwoYe8OUKrW0ura5txnFyWLSzs1MXXXRRdNFFq5/3hHVqVnLu3yBJz2jl2W5XS/qluz9gZh9WMgXAxyR9TNJoJQF0mLu/amatZrarkogZ4e7/MrNSSb929zPNbLKkCe5+uZmd6e6nZ0eWOiWNU3IofaqZdf0d+ai7LzKzS5VcGXe5pD3c/RIzuyi7f0dJ2t7d9zWzIUpCbLVYcveblYSV9ho2jEMh2CqIJWxzGuZM/dMVV1yxXVlZ2Sn9+vXTsmXLFEVR99IqXWEQx7FKSkpkZvrh2Ud/5YwzzijakO3UzZ6+QNI3JOm+++47uWtm8ry8PEla5QTvrjmj4jhWFEVqaXnPF7RtsMKiwmUFBQWDCgsLlZubq1Qqtdpjur4mklRYWKi6V6b+ZGO2+c17TijNLSo6a12x5O7dgUYovWf/zF4NuMjMlitZ/01KZsN+XJLc/d9mVu7ujWZ2hpL4KJV0T/axtyiJmvoet0nJqJMkNSgZierteXd3M+ua8XuQpKGSfpf9/10o6eG17PcoSU9l9+8drXvyUmCrIZawTRoxYsTtpaWlqxx6WtPIUn5+vlKplNra2golaejICYNSqXh3yfaMIh+o2IZIXmAWZdzjtGRr/M32xp+docLCwrWO3nSFQSaTkbursHCjTwlab0/e/OQhJ1x1Qt3avhaSVllepqVt40Puui9Obvr+tO+ruLhYBQUFiqJotQWT4zjujsfW1laZWdT7ykesl30kKTs6UyJpcfb22ZIOlPRgdmTpHTMbKmlvdz8me9J1g5ndLel+Sd9VsjbhcT1ee03f7xa4f7GktyUd5e7LzexISU3qfX3nyv07Lrvv/SX9zt3Hrf/bBrYcYgnbpOLiYoUOPXWdt9R16Km1tVWV1ePuyE0pJWmF3JZaHL9rsr83qeC5xbMfaAxt7+qrr/5HZeXwSaElXNxdnZ2diuNY/fr122KHE9587s2SVCql/Px8FRYWdo+w9b5Sr2t5mfvPu//ynx3xs43dbHF+Yb7y8vJUVFS02v+Drljt7OxUOp1W3Bl3Kjm35j1Ni/A/rtDM/q4klE5Vcq6RlMTPLWb2XSVf25OUHGYbYmYvKImYq909LSltZo9LGuTuS9axvVfM7B5Jj/S+Izsr+LckTclOD7JCyWSYKyTlmdlVkrrWV3xQ0iFm9k8l/xZt+DomwBbC1XDYJj366KNeUlKyyshSb3EcK5PJqKWlRe+++64OPfTQ93w1h5kd8MorrzxRVla2xiDpGkVpb2/X0qVLNXLkyC+5+93v/R1u0L7tMGDwgCsumHLBif3791deXp5SqZSiKFolWlpaWrS8cbnOP+D8E9x98kZuM0/SJVc9ddW5AwYMUEFBQfeiyVISqZlMRq2trVqyZInO/+j5F0u6hJGlrSd70vb97v73rb0va8PVcH3aNn01HEsFYJu0cOHC2wsKClRUVKTi4uLukabef4qKilRQUKB+/fpt7CYX3Xjjjf/ouaZd16GnnoHQ0dGh5uZmSVrXb++b0rvLFi2b9cRvnnhu2bJlam5uVltbW/eflpYWNTY2aunypTr/gPOvkzRjYzeYnePn8fMOOO+mpUuXqqmpSS0tLd3bbG1tVVNTU1coXS/pr4TS1pO9oq2oL4cSsDUxsoRtkpl94uWXX35k4MCB3SM9vU+4zmQyam9v17JlyzRy5MiT3P32wEuua3s5ks6+6aabTj7ssMNGdp3o3TV60zVy09jY6HvsscfPJF3m7ss3wVtd3/0rlnSapJ0vfuTir+cXJedqxXGsWLEevPLBh5964Kn/SLpL0kubYvkQM0spOSfloG/f++3PDa4cPKDr5PJ0Jq36WfXzbjr9pgckTZM0lVjCujCy1Kdt0yNLxBK2SWa2y84773z+I4888uXS0tLuE7l7Xg3XFTBvv/12esyYMUe6e3CtufXYZqmSc0bKX3311TPdvTtIzEwjR46cqGREaaK7122K97mB+5eS9CFJY5VcsZSn5MTb5ZL+LenRTR1w2Vnah0s6RNIIJVdMScn6dnOUnPcyn7Xd3l/MbKCkH7v7BpVL9kTue5RciZcn6Wx3f8rMLlEyaeUroeePHj3aZ86c+V53G5sXsQS832RHek7eb7/9Drn77ruP7Tn/kSR1dHSopaVF9fX1LWPHjr1WyQ/+pk2w3ZSk3ZWse7ejkhNr05IWSXpayaXWrWt/hc0vGzC5klJKrmbq2NyjOtltRkr+gZSkDkkxkfT+ZGY3KpmE8sUNfN7Fkpa6+y/MbJSSeZw+YmYDJE1298NDzyeW+rRtOpa4Gg7bJHdPm9mdTz/9dHrkyJGLJk+efPjee+9dYWbKZDJavnx55/777z9RScRM3BShlN1uRtIsSbOywRYpCZJ0XwmD7H5slnXD1rHNjFZeCYV1yE76eJSSUZjtlZwA/3sz+4+kV5VcOfh1rRypyZH0Q3f/u5l9UlLXsacXlByCPVCrz6w9QtKdSiaXTCu5cq1D0m+VfO/mSjrN3V/qsV+lkvZx969nP39N0hNK5k1aKOlYSccomWSyp3Mk/Vwrr3jMkdQmSe6+zMzazGzPDQ0wYEsglrDNcvc2M7tN0l+/8IUv/FPJPwyFSv4xmKdkpGdO9tLpzbH9zfK6+J9SIulQJZM9PmtmD2Rvu9TdXzCzq5WcHH+tme0o6Z/ZEZtfSdo3O4v2jySVa80za+dJek7S2Upiajsls3Evl/R5SbsqCbGe9lNyCLXLByR93N0bzOwJJSF1v5K5m9YoOyfUPZLO6nHzi0oOERNL6HOIJWzTsiMadZLqsvO+dN3OycR4P3gs+7260MyWKokmaWWsVEuaLEnuPs/MVig5H22puy/K3n6JmQ3WmmfWvlzSuZKmKwmk85WccD9SydIjnZIu67VP2ysZQeqy2N0bsh83SCows09rDSNL7v6sme0h6TeSvuvuj/W4f4GSQ9dAn0Ms4X8GgYT3ob2lZK4sJSM8i7K3d30vd83S/UJ2ZGk7JdExwMzK3H2Jmf1SySjOmmbWPkrSP9z9YjP7nJJwulvSAnc/zMz2l/RjJevJdVkkaUCPz1c7vLy2kaXs+nP3SfqMu8/qdfd2Pd4f0KcQSwDQdw0xs79J6i/pdHfP9Fqq5seSbs+O5BRK+pq7d5jZ6Upm0c4oOWfpX5LWNLN2P0n3mFlaSYB9W8lI7G/N7Cwl5zdd0mufnpZ01Xt8P1coWWPu2uz7WO7uR2XvG6NkZAvoc7gaDgD6oOwJ3lXuft7W3pfezOwmJRdGvLDOB6/f65VJusvdjwg9jqvh+rRt+mo4YgkA+qA+HkuDJV3u7qdsote7VMlCui+t43GNWvXk8r5me61cyLiv2dz7ttjdx2/G19+qiCUAwPuCmc1099Fbez/Wpi/vX1/et/cD1oYDAAAIIJYAAAACiCUAwPvFzVt7B9ahL+9fX963Po9zlgAAAAIYWQIAAAgglgAAfZqZRWZ2k5k9ZWYzzOyDW3ufJMnMxpjZjOzHHzSzf5rZP8zsxp7LK22F/co1s7uz+/KsmR3Zl/bv/YgvFgCgrztaUoG77y/pPEk/28r7IzM7R9KtSmYkl6RrJP3Q3Q9UMkHjUWt77hZwgqT/ZvdlgpKFlfvS/r3vEEsAgL7uo0oW+5W7Py2pL8wX9IakT/X4fG9JXQsDT5N0yBbfo5Xuk3RBj8/T6lv7975DLAEA+rpSSct7fJ4xs626tqm7/15SZ4+bzFdeMdWoZD2/rcLdm9y90cz6KVnQ+Id9af/ej4glAEBft0LJor9dIndPb62dWYu4x8f9JC3bWjsiSWZWLulRSXe7+73qY/v3fkMsAQD6uick1UiSme0nKbiG3FbygpmNzX48QdI/ttaOmNkOkh6WdK673569uc/s3/sR8ywBAPq07JVbN0jaU8nJySe6e+3W3SvJzIZL+o2772dmu0i6RVKepNmSTnH3zFbar2slfUZSz6/RtyT9Un1g/96PiCUAAIAADsMBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAA9GJm72zEc48xs3t7fL6fmT1jZk+Y2YU9bh9oZhOzH8/ILno7o8efvF63P2Fm95jZwOxz7jSzF7P3PWZm/zGzE7P31ZjZV9/7VwA9bdXp4gEA2JZk5zgaJ+nfPW6+SdKxkt6UNMXMPuLuz0u6TNL1hIJi8gAACyhJREFUPR73pd7zR5nZKreb2Rck3Zx9PUk6x92nZ+8rk/Symd3p7lPNbJqZ/d7dey4Vg/eAkSUAANbCzPYys39mR24eMrOK7O0XmNnz2dv+0WN27Cclfb3H80sl5bv7G9m12R6S9Ins7fu4+4sbsj/uPlnS3mZWsIa7h0hq67EG3FRJX96Q18eaEUsAAKzdLZLOcPeDlcwifo2ZfUjJkiH7SDpa0tCuB7v7byX1nO25VMnadl26FrHdT9KcXtua1OMQ3EmBfVoqaUD2459kY61e0jWSjuvxuBcljV2vd4kgDsMBACDJzC6T9NHsp6nsf4e5e9chtcclXSmpWtKz2eVCWs1sZuBley8C3LWI7faSFvZ67GqH4dawj6ZkBGlR9qZz3H26mdVIukrSGz0evkDSwNDrYf0wsgQAgCR3/6G7j3X3sZK61k2bb2Z7Zj8+WNKrkl6WtI+ZRWaWL2mvwGuukNRhZjtnQ2eckkVsF2nl6NCGOEnS39w97rWdqZL+pOR8pi7baWVUYSMwsgQAwNqdIulX2dBJSzrJ3d80s6mSnpa0WFJn9s/anCZpspLRqofd/RkzK1EyErQ+JplZc/bjeZK+sZbHXSrpeTM73N2nSBoj6W/ruQ0EsJAuAAAbwMwGS/q0u9+QHVl6WdLH3b1+A1/nJkkT3f2FzbSf0yUdnx3dwkbgMBwAABtmsZLDcP9Sckjt1g0NpawfSTp9k+5ZlpkdLun3hNKmwcgSAABAACNLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQQCwBAAAEEEsAAAABxBIAAEAAsQQAABBALAEAAAQQSwAAAAHEEgAAQACxBAAAEEAsAQAABBBLAAAAAcQSAABAALEEAAAQ8P+eOdxxdel2VQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 612x792 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Generating figure\n",
    "fig = plt.figure(figsize=(8.5, 11))\n",
    "ax1 = plt.subplot2grid((3, 4), (0, 0), colspan=3)\n",
    "ax2 = plt.subplot2grid((3, 4), (0, 3))\n",
    "ax3 = plt.subplot2grid((3, 4), (1, 0), rowspan=2, colspan=4)\n",
    "\n",
    "# STRING diagram\n",
    "diagram = mpimg.imread('Fig7_stringperoxhits_STRING.png')\n",
    "\n",
    "ax1.imshow(diagram)\n",
    "ax1.axis('off')\n",
    "ax1.annotate('A)', (-0.15, 1.1), xycoords='axes fraction', fontsize=12)\n",
    "\n",
    "pos1 = ax1.get_position()\n",
    "points1 = pos1.get_points()\n",
    "new_points1 = points1 - [[0.1, 0], [0.1, 0]]\n",
    "pos1.set_points(new_points1)\n",
    "ax1.set_position(pos1)\n",
    "\n",
    "# Horizontal bar plots for G.O. terms\n",
    "perox = [1.68e-23, 2.83e-13, 9.87e-10, 3.78e-7, 6.05e-6, 0.00026]\n",
    "perox_log = -1*np.log10(perox)\n",
    "perox_terms = ['peroxisome\\norganization (n=13)', 'fatty acid metabolic\\nprocess (n=11)', 'peroxisomal membrane\\ntransport (n=5)', \n",
    "               'long-chain fatty-\\nacyl-CoA metabolic\\nprocess (n=4)', 'fatty acid beta-\\noxidation (n=4)', \n",
    "               'ether lipid\\nbiosynthetic\\nprocess (n=2)']\n",
    "\n",
    "ax2.barh(range(6,0,-1), perox_log, 0.5, tick_label = perox_terms, \n",
    "         color=['black', 'black', 'tab:olive', 'forestgreen', 'blue', 'lightcoral'])\n",
    "ax2.set_xlabel('-log10(FDR)')\n",
    "ax2.set_title('Peroxisomal from WT (21 hits)')\n",
    "ax2.set_xlim(0, 26)\n",
    "ax2.annotate('B)', (-1.2, 1.1), xycoords='axes fraction', fontsize=12)\n",
    "\n",
    "right_side = ax2.spines['right']\n",
    "right_side.set_visible(False)\n",
    "top_side = ax2.spines['top']\n",
    "top_side.set_visible(False)\n",
    "\n",
    "pos2 = ax2.get_position()\n",
    "points2 = pos2.get_points()\n",
    "new_points2 = points2 - [[0.08, 0], [0.08, 0]]\n",
    "pos2.set_points(new_points2)\n",
    "ax2.set_position(pos2)\n",
    "\n",
    "ax3.axis('off')\n",
    "\n",
    "fig.suptitle('Figure 7', x=0, y=1, fontsize=12)\n",
    "\n",
    "plt.savefig('Figure7.pdf', dpi=500, bbox_inches='tight')\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}