{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Decay curve\n", "\n", "`.diagonal_mean` method in `Cool` and `DotHiC` can use for calculate the decay curve from a contact matrix.\n", "By using this auxiliary function, it's easy to produce the decay curve plot within specific region." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "import coolbox\n", "from coolbox.api import *" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'0.3.0'" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "coolbox.__version__" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot a decay curve of chr1 and chr2:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAFVCAYAAACU4eUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXhV1dn38e/KPBKGMCUCEgYBJzAJolZB1ApaKmqtrVKxUGi1aoXa1/KgYEGxtoI8KE8LWEuhtqjVWhHr0IpIFWUQZFaUMYnVQJQxCSS53z92EnPISchJzsn4+1zXvpKz9n32uXfCBTdr7bWWMzNEREREpGUKa+gERERERKThqBgUERERacFUDIqIiIi0YCoGRURERFowFYMiIiIiLZiKQcA5t9A5d9w5d6TCMayh8xIREREJtYiGTqARmW9mdzR0EiIiIiL1ST2DIiIiIi1YkyoGnXOTnHPPOed2OufMObe7mtgw59wE59x251yBc26fc26mcy6+irfc7JzLc85tc85Nds6p11RERESaPdeUdiBxzhmQB3wApAOHzOz0KmL/F7gL+DvwT6AvcCewErjczEoqxJ4HZAH7gfOAvwJLzOz+kN2MiIiISCPQ1IrBNDPbWfr9ZiDBXzHonDsT2AT83cyur9B+JzAHuNnM/lLN59wE/MrMegX5FkREREQalSY1TFxWCNbA9wEHzD6pfQFwDBh1iveXlL5fREREpFlrUsVgADLxCrrVFRvNrADYUHq+nHPuRudckvOcA0wFnquvZEVEREQaSnMtBlOA/WZW6OdcNpDsnIuq0HY7sBs4DDyP98zglFAnKSIiItLQmuuM2TjAXyEIUFAh5jiAmQ2u6YWdc+OB8QDx8fHpffr0qUOaIiIiIvVj3bp1+82s/cntzbUYPAZ0qOJcTIWYgJnZfGA+QEZGhq1du7Y2lxERERGpV865Pf7am+swcQ7eUHC0n3OpeEPIx2t7cefcCOfc/IMHD9Y6QREREZHGoLkWg2vw7m1gxUbnXAzQH6hTd56ZLTWz8UlJSXW5jIiIiEiDa67F4DOAAXef1D4O71nBp+tycfUMioiISHPRpJ4ZdM79AOhW+rI9EOWcu6/09R4zWwxgZpucc3OBO5xzLwCv4O1AchewAqhywemaMLOlwNKMjIxxdbmOiIiISENrUsUgMBY4eebv9NKvK4DFFdrvxlsuZjxwNd5Wc48DUypuRSciIiLSkjWpYtDMhgQQWwzMLD2Cyjk3AhjRs2fPYF9aREREpF4112cGQ0oTSERERKS5aFI9gyIizU1hYSF5eXkcPnyY4uLihk5HRJqQ8PBwEhMTadu2LdHR/lbTqxkVg7WgYWIRCYbCwkL27t1LmzZtOP3004mMjMQ519BpiUgTYGacOHGCQ4cOsXfvXrp27VrrglDDxLWgYWIRCYa8vDzatGlDcnIyUVFRKgRFpMacc0RFRZGcnEybNm3Iy8ur9bVUDIqINJDDhw/TqlWrhk5DRJq4Vq1acfjw4Vq/X8WgiEgDKS4uJjIysqHTEJEmLjIysk7PHKsYrAXtQCIiwaKhYRGpq7r+PaJisBb0zKCIiIg0FyoGRUSk0XjggQdwzrF79+6GTkWkxVAxKCIizdYHH3zANddcQ7t27YiJieHMM89k9uzZdXq+6t5778U5R0JCgt/zhYWFTJkyhe7duxMdHU2PHj148MEHOXHihN/4RYsWMWDAAGJjY+nYsSM/+tGPyM3N9Rv7/vvvc/nll5OYmEirVq0YNmwYGzZs8Bubk5PDLbfcQvv27YmNjSUjI4Pnnnuu3nP2Z/fu3Tjn/B5nnXVWs7n3UOUcdGamI8ADGAHM79mzp4mI1NbWrVsbOoVGZ+rUqQbYrl276nytFStWWHR0tHXo0MEeeOABmzt3ro0cOdIAGzduXK2uuX79eouIiLCEhASLj4/3G3PNNdcYYGPGjLEFCxbYmDFjDLDRo0dXip01a5YBNnjwYJs3b57df//9Fh8fb/369bMjR474xK5atcqio6MtLS3NZs2aZbNmzbK0tDRLSEiwjRs3+sQeOHDAunfvbvHx8Xb//ffbvHnzbPDgwQbYU089VW85V2XXrl0G2LXXXmuLFy/2OZYuXVopvineeyhz9qcmf58Aa81fXeOvUUfNjvT09FP+4EVEqqJisLJgFIOHDh0yM7Nzzz3XYmNj7dNPP/U5P378eANs5cqVAV23qKjIMjIybMSIETZ48GC/xeCyZcsMsIkTJ/q0T5w40QB75513yttyc3MtLi7OMjMzraioqLz9pZdeMsAeeughn2tkZmZaYmKiZWVllbdlZWVZYmKiXXHFFT6xv/jFLwywl156ySf/zMxMa9u2rR0+fLhecq5KWTE4derUGsU3xXsPVc5VUTGoYlBEmqCWVgwWFhbaI488Ul6ktWrVytLT0+3xxx8vjykrBrdv326TJk2y1NRUi4qKsnPOOceWLVvmc72KBcWSJUvsvPPOs5iYGBs9erTl5eUZYFdeeWWlPN59993ynqCKcnJybNu2bXb06FG/+c+aNcvi4uJs9+7dVRaDN998swG2d+9en/a9e/caYLfddlt524IFCwywRYsWVbpOWlqa9e3bt/z1jh07/OZsZjZmzBhzztlnn31W3paammo9evSoFLto0SID7Jlnngl5ztWp+LvLz8+v8mdu1jTvPZQ5V6UuxaCeGRQRkZA7fvw4V155Jffeey8dO3Zk2rRpPPTQQ6Snp/PCCy9Uih89ejQrV67knnvuYfr06eTm5jJy5Ei/E0tefPFFbrvtNoYNG8acOXMYPnw4hYWFAMTFxVWKL2t77733fNonTZpE3759Wb16daX37Nmzh/vvv5+pU6fSrVu3Ku9zzZo1pKam0qVLF5/2Ll26kJKSwpo1a3xiAS644IJK1xk0aBDbt2/nyJEjNYo1M9atWwfAZ599RnZ2NoMGDfIbW/F6ocy5JmbOnElcXBzx8fF06dKFKVOmlP/uavqZjfHeQ5lzKGhvYhERCbnZs2fz1ltvMWnSJGbMmOFzrqSkpFJ8cnIyS5cuLV8/7dJLL2XgwIHMmzePhx9+2Cd2y5YtbNy4kb59+5a3mRnJycm899575OfnExsbW35u+fLlAOzbt6/G+d92222kpaUxceLEauNycnLo16+f33OpqalkZWX5xJa1+4s1M3Jycujdu/cpYwGys7NrdN2KsaHMuTphYWEMHTqUkSNH0q1bN3Jzc3n22WeZPn06q1at4tVXXyU8PDzg+2ks9x7KnENBxaCISCPzq6Vb2JpzqKHT8NEvpRVTR5xZ6/c//fTTtGnThilTplQ6FxZWeZDqZz/7mc9CupmZmSQkJLBjx45KsVdffbVPIQjeIrwTJkxg8uTJXHfddUybNo3k5GT+9a9/MXXqVCIiIjh27JjPexYuXMjChQsrXf+vf/0rr776Kv/5z3+IiKj+n81jx44RHR3t91xMTIzPZ5Z97y8+JibGJyZUsaHMuTpdu3bl3//+t0/b2LFjGT9+PAsWLGDJkiXcfPPNAX9mY7n3UOYcChomrgXtQCIiEpgdO3bQp0+f8n/cTiUtLa1SW7t27Thw4ECl9qp6oX75y18yefJk3nrrLQYOHFjes/foo4/Spk2bGu0LnZeXx913383YsWO58MILTxkfFxdXaZizTEFBgc+wddn3/uILCgp8YkIVG8qc8/Ly+O9//+tznMrkyZMBWLZsWa0+s7HceyhzDgX1DNaCmS0FlmZkZIxr6FxEpPmpSw9cc1E2RHgy7xl4X1X9QxkWFsaDDz7IpEmT2LRpE2bGueeeS0lJCT/+8Y/9PqN1sl/96lccPXqUcePG8cknn5S35+fnY2Z88sknREdHlz9zlpKSUuWQXnZ2ts9QYEpKSnl7z549K8U658pjKsb6uy58PaQYSGwoc77uuutYsWKFT4y/319FXbp0ITw8nP379/v9zFPdT2O591DmHArqGRQRkZDr3bs327dvr7IXJpTi4+MZNGgQF1xwAXFxcbz66quYGVddddUp37tnzx6OHj3K+eefT69evcqP1atXc+zYMXr16sXw4cPL4zMzM8nOzq70POK+ffvIyckhIyPDJxZg1apVlT73vffe44wzzihf2PpUsc450tPTAejcuTOpqamVJsiUxQKV8ghFzjNnzuSNN97wOU5l586dFBcX07Fjxxp/ZmO891DmHBL+phjr0NIyIhJ6LWlpmd/85jcG2H333VfpXElJSfn31a0z2K1bNxs8eHD560DXqjMz279/v/Xq1cuSk5MtNzfX55y/pWXeffdde+655yod/fr1s5iYGHvuuefs9ddfL49/+eWXq123ruLahl988YXFxsbawIED/a5bN336dJ9rZGRkWGJiomVnZ5e3ZWdnW2Jiol122WU+sffcc0+V69a1bt26fC3GUOdclf3791dqKy4uthtvvNHvUipN8d5DlXNVtM6gikERaYJaUjFYWFhoF198cfnafzNnzrS5c+fa7bff7vMPYzCLwWXLltmQIUPsoYcesj/84Q82ZcoUS0lJsdjYWHvjjTcqxY8ePdoAW758+Snvp6p1Bs3MvvWtbxlgY8eOtSeffNLGjh1rgI0aNapS7KOPPmqADRkyxObNm2dTpkyx+Ph469OnT6WFht955x2LioqytLQ0e+yxx+yxxx6ztLQ0i4+Ptw0bNvjE7t+/37p162YJCQk2ZcoUmzdvng0ZMsQAe/LJJ+st56pce+21dtlll9n9999v8+fPtxkzZlh6eroBds0111hxcXGTv/dQ5uyPikEVgyLSBLWkYtDMLD8/3x588EHr16+fRUdHW1JSkmVkZNjcuXPLY4JZDG7ZssWGDRtmnTp1ssjISOvcubONGjXKtm/f7jc+WMVgfn6+TZ482bp162ZRUVHWvXt3mzZtmh0/ftxv/B//+Ec755xzLDo62tq3b28//OEP7fPPP/cb++6779rQoUMtPj7eEhIS7Jvf/KatW7fOb2xWVpaNGjXK2rVrZ9HR0TZgwABbsmRJvefsz5NPPmmDBw+2jh07WmRkpCUkJNj5559vc+fOrVQINuV7D1XO/tSlGHTeOamNjIwMW7t2bUOnISJN1LZt2yotiSIiUhs1+fvEObfOzCo9gKgJJCIiIiItmIrBWtA6gyIiItJcqBisBTNbambjk5KSGjoVERERkTpRMSgiIiLSgqkYFBEREWnBVAyKiIiItGAqBks552Kdc5845440dC4iIiIi9UXF4NemAXsaOgkRERGR+qRiEHDOpQPDgEcaOhcRERGR+tRkikHn3CTn3HPOuZ3OOXPO7a4mNsw5N8E5t905V+Cc2+ecm+mci/cTGwEsAH4KHA/dHYiIiIg0Pk2mGARmAEOBT4EvTxH7GDAL2ArcCTwH3AUsdc6dfM+/ANab2dvBTVdERESk8Yto6AQC0MPMdgI45zYDCf6CnHNn4hWAL5jZ9RXadwFzgO8Bfylt6wn8BBgQ2tRFREREGqcm0zNYVgjWwPcBB8w+qX0BcAwYVaHtG0BH4GPn3H7gH0C8c26/c+6SOqYsIiIBeuCBB3DOsXv37oZORaTFaDLFYAAygRJgdcVGMysANpSeL/Ms0BPoX3r8CK9g7A+8Xx/JiohI6HzwwQdcc801tGvXjpiYGM4880xmz55NcXFxja9x66234pzze/ztb3+rFF9YWMiUKVPo3r070dHR9OjRgwcffJATJ074vf6iRYsYMGAAsbGxdOzYkR/96Efk5ub6jX3//fe5/PLLSUxMpFWrVgwbNowNGzb4jc3JyeGWW26hffv2xMbGkpGRwXPPPec3NpQ5VycrK4vx48fTtWtXoqOj6dSpE8OHD2fr1q2VYl955RUuvPBC4uPjadu2LTfccAO7du3ye92PPvqIkSNH0qZNG+Lj47n44ot58803/cYePHiQO++8k9TU1PI/I7/73e8ws0qxJSUlPPbYY/Tp04eYmBi6dOnCz3/+c44ePer32qHKOejMrMkdwGZgdxXnNgGfV3HuWcCAqCrODwGOnOKzxwNrgbVdu3Y1EZHa2rp1a0On0OhMnTrVANu1a1edr7VixQqLjo62Dh062AMPPGBz5861kSNHGmDjxo2r8XVGjx5tgC1evLjSsWfPnkrx11xzjQE2ZswYW7BggY0ZM8YAGz16dKXYWbNmGWCDBw+2efPm2f3332/x8fHWr18/O3LkiE/sqlWrLDo62tLS0mzWrFk2a9YsS0tLs4SEBNu4caNP7IEDB6x79+4WHx9v999/v82bN88GDx5sgD311FP1lnN1PvjgA2vbtq1169bNpk6dan/4wx/st7/9rd1444321ltv+cQ+//zz5pyz/v3729y5c23GjBnWoUMH69y5s2VnZ/vEfvLJJ9a2bVvr0KGDzZgxw+bOnWv9+/e3iIgIe+ONN3xiCwsLLTMz0yIiImzChAk2f/58u/baaw2wqVOnVsr5rrvuMsCuvfZamz9/vk2YMMEiIiLs0ksvteLi4nrJuSo1+fsEWGv+aht/jY39OEUx+Cmwt4pzi0qLwdbByCM9Pf2UP3gRkaqoGKwsGMXgoUOHzMzs3HPPtdjYWPv00099zo8fP94AW7lyZY2uV1YM1sSyZcsMsIkTJ/q0T5w40QB75513yttyc3MtLi7OMjMzraioqLz9pZdeMsAeeughn2tkZmZaYmKiZWVllbdlZWVZYmKiXXHFFT6xv/jFLwywl156qbytqKjIMjMzrW3btnb48OF6ybkq+fn51rNnT+vfv78dPHiw2tjjx49bSkqKde3a1Sfv9evXW1hYWKXC/oYbbrCwsDBbv359edvhw4eta9eu1rt3byspKSlvnzt3rgE2Z84cn2tcd911FhkZabt37y5v27x5sznn7LrrrvOJnTNnjgH29NNP10vOVVEx6Huu1j2DAXz+CGB+z549T/mDFxGpSksrBgsLC+2RRx4pL9JatWpl6enp9vjjj5fHlBWD27dvt0mTJllqaqpFRUXZOeecY8uWLfO53q5du8p7cJYsWWLnnXeexcTE2OjRoy0vL88Au/LKKyvl8e6775b3glWUk5Nj27Zts6NHj/q0lxWDJSUldvDgwUo9QBXdfPPNBtjevXt92vfu3WuA3XbbbeVtCxYsMMAWLVpU6TppaWnWt2/f8tc7duzwm7OZ2ZgxY8w5Z5999ll5W2pqqvXo0aNS7KJFiwywZ555JuQ5V+dPf/qTT7FaUFBgBQUFfmPfeOMNA2zatGmVzg0dOtRatWplx48fNzOzI0eOWHR0tA0dOrRS7LRp0wyw999/v7ztoosusri4OMvPz/eJffvttw2wRx55pLxt8uTJBtjbb7/tE5ufn29xcXE2fPjwesm5KnUpBpvjM4M5QLJzLtrPuVRgv5nVaT1BM1tqZuOTkpLqchkRkRbj+PHjXHnlldx777107NiRadOm8dBDD5Gens4LL7xQKX706NGsXLmSe+65h+nTp5Obm8vIkSP9Tix58cUXue222xg2bBhz5sxh+PDhFBYWAhAXF1cpvqztvffe82mfNGkSffv2ZfXq1ZXeA5CUlERSUhKxsbFcccUVvP9+5UfL16xZQ2pqKl26dPFp79KlCykpKaxZs8YnFuCCCy6odJ1Bgwaxfft2jhw5UqNYM2PdunUAfPbZZ2RnZzNo0CC/sRWvF8qcq/PKK68A0Lp1ay655BJiY2OJiYlhwIABvPbaaz6xp/rMQ4cO8fHHHwOwceNGCgsLq4yteL2SkhI++OADBgwYQExMjE/swIEDcc5VuvewsDAGDhzoExsTE0P//v0D+jnVNudQaUpLy9TUGuCbwEBgZVmjcy4Gb2KI1hMUEalns2fP5q233mLSpEnMmDHD51xJSUml+OTkZJYuXYpzDoBLL72UgQMHMm/ePB5++GGf2C1btrBx40b69u1b3mZmJCcn895775Gfn09sbGz5ueXLlwOwb9++GuXeqVMnJkyYQHp6OvHx8Xz44YfMnj2biy++mFdeeYXLL7+8PDYnJ4d+/fr5vU5qaipZWVk+sWXt/mLNjJycHHr37n3KWIDs7OwaXbdibChzrs5HH30EwPXXX8/555/PkiVLyMvL46GHHuKqq67itddeK/+51vR+zjzzzIDu/csvvyQ/P99vbHR0NMnJyZV+TsnJyURHV+5rSk1N5d133+X48eNERUWFLOdQaY7F4DPA/wB3U6EYBMYBccDTdf0A59wIYETPnj3reikRkcr++Uv476aGzsJXp7Nh+K9r/fann36aNm3aMGXKlErnwsIqD1L97Gc/Ky8EATIzM0lISGDHjh2VYq+++mqfQhDAOceECROYPHky1113HdOmTSM5OZl//etfTJ06lYiICI4dO+bznoULF7Jw4cJK1//1r33ve+TIkdx0003079+f2267zSenY8eO+S0WwOtBqviZZd/7iy/rqSqLCVVsKHOuzuHDhwHo06cPL730Uvnv+rLLLqNfv35Mnjy5vBhsiJ9TWXwgP6eymKioqJD+vkKhyRSDzrkfAN1KX7YHopxz95W+3mNmiwHMbJNzbi5wh3PuBeAVoC/eDiQrKF1wui7MbCmwNCMjY1xdryUi0hLs2LGD/v37VxqOq0paWlqltnbt2nHgwIFK7VX1Qv3yl7/k2LFjzJw5s3xoLyEhgVmzZjF58mSKiooCuANfvXr14rvf/S4LFy7k448/Ls8hLi6ufIj6ZAUFBT7D1mXfFxYW+vRclsVWjKkY6++6tY0NZc55eXkcP+77VFanTp0Ayt97yy23+BT9vXr14sILL2TlypUcPXqU+Pj4kN17dbFV3fsXX3xRZWwo8vD3+wqFJlMMAmOBwSe1TS/9ugJYXKH9bmA33jIwVwP7gceBKWZWeTwiQOoZFJGQqkMPXHMRHh7ut917Bt5XVf9QhoWF8eCDDzJp0iQ2bdqEmXHuuedSUlLCj3/8Y7/P1AXi9NNPB2D//v3lxWBKSkqVQ3rZ2dk+Q4EpKSnl7Sf/e5KdnY1zrjymYqy/68LXQ4qBxIYy5+uuu44VK1b4xJT9/k477TQ2b95cXhxW1LlzZ8yMgwcPEh8f7/OZJ/cA1+Xe27RpQ2xsrN/YwsJC9u/fz+DBX5cdKSkpbN26lcLCwkq9eNnZ2SQnJxMVFVUpj2DmHCpNZgKJmQ0xM1fFMeSk2GIzm2lmZ5hZtJmlmtlEMzv1U601y0UTSEREAtC7d2+2b99eZS9MKMXHxzNo0CAuuOAC4uLiePXVVzEzrrrqqjpdt2x4uGPHjuVtmZmZZGdnV3oecd++feTk5JCRkeETC7Bq1apK137vvfc444wzSEhIqFGsc4709HTAK6ZSU1MrTZApiwUq5RGKnGfOnMkbb7zhc5Qp66mt+DximaysLCIiImjbtm2NPrNVq1blxfjZZ59NdHR0lbEV7z0sLIzzzjuP9evXV/pzuXr1asys0r2XlJRUmmBUUFDAhg0bAvo51TbnkPE3xVhHzQ6tMygiddGSlpb5zW9+Y4Ddd999lc5VXEOtunUGu3XrZoMHDy5/XXFpmZrav3+/9erVy5KTky03N9fnnL+lZY4cOVJp2REzb8HkqKioSkupvPzyy9Wu2VdxbcMvvvjCYmNjbeDAgX7X7Js+fbrPNTIyMiwxMdFnweLs7GxLTEy0yy67zCf2nnvuqXKdwdatW5evxRjqnKuyfft2Cw8Pt/T0dDtx4kR5+4YNGywsLMxn3cTjx49b586dK63ZVxY7duxYn2t/5zvfsbCwMNuwYUN5W9mafb169fL58/bEE09Uuc5gRESEz5/DjRs3VrvO4OLFi+sl56q0uHUGG/pA6wyKSBC0pGKwsLDQLr744vK1/2bOnGlz586122+/3aeQCWYxuGzZMhsyZIg99NBD9oc//MGmTJliKSkpFhsb63dXh7L1BJcvX17etn79euvUqZP95Cc/sZkzZ9rvf/97u+222yw6OtpiY2P9Llz9rW99ywAbO3asPfnkkzZ27FgDbNSoUZViH330UQNsyJAhNm/ePJsyZYrFx8dbnz59fIoIM7N33nnHoqKiLC0tzR577DF77LHHLC0tzeLj432KCDOv6O3WrZslJCTYlClTbN68eTZkyBAD7Mknn6y3nKszadIkA+zCCy+0OXPm2AMPPGBt27a1hIQE+/DDD31in332WZ/dPB5++GHr0KGDdezY0WcRbjNvTcY2bdpYhw4d7OGHHy7fzSM8PNxeffVVn9jCwkJLT0+3iIgImzhxoi1YsKB8BxJ//3G54447yncgWbBggU2cONEiIiJs8ODBldafDFXOVVEx2ECHegZFpC5aUjFo5i3O++CDD1q/fv0sOjrakpKSLCMjw+bOnVseE8xicMuWLTZs2DDr1KmTRUZGWufOnW3UqFG2fft2v/H+isHPPvvMRo0aZWeccYYlJiZaRESEdenSxW655Rbbtm1blfc5efJk69atm0VFRVn37t1t2rRp5YsMn+yPf/yjnXPOORYdHW3t27e3H/7wh/b555/7jX333Xdt6NChFh8fbwkJCfbNb37T1q1b5zc2KyvLRo0aZe3atbPo6GgbMGCALVmypN5zrs68efPKr9O6dWu79tprbfPmzX5jly5daueff77FxsZa69at7frrr7dPPvnEb+zWrVvt29/+tiUlJVlsbKxddNFFVW7r9uWXX9pPf/pT69y5c3lv7+OPP+63N66oqMgeffRR6927t0VFRVlKSopNmDChyiI4VDlX9f5TqaoYdN45qY2MjAxbu3ZtQ6chIk3Utm3bKj1cLiJSGzX5+8Q5t87MKj2A2GQmkIiIiIhI8KkYrAXn3Ajn3PyDBw82dCoiIiIidaJisBZMS8uIiIhIM6FiUERERKQFUzEoIiIi0oKpGKwFPTMoIiIizYWKwVrQM4MiIiLSXKgYFBFpQFrrVUTqqq5/j6gYFBFpIFFRUeTn5zd0GiLSxOXn5xMdHV3r96sYFBFpIMnJyWRlZZGXl8eJEyfUSygiNWZmnDhxgry8PLKysmjXrl2trxURxLxERCQASUlJREdHk5uby4EDBygqKmrolESkCYmIiCAmJoauXbsSExNT++sEMacWwzk3AhjRs2fPhk5FRBphtU0AACAASURBVJq4mJgYunTp0tBpiEgLpmHiWtBsYhEREWkuVAyKiIiItGA1HiZ2zoUD0WZ2rEJba2As0BZYYmabgp+iiIiIiIRKIM8MzgMGAWcBOOcigf8A/UrPT3TOXWBmG4KbooiIiIiESiDDxN8AXqrw+jt4heBPgQuBz4FfBi81EREREQm1QHoGOwO7Kry+GthiZr8DcM7NB34cxNxEREREJMQC6Rl0QHiF10OA5RVefwZ0CEJOIiIiIlJPAikGdwFXAjjnLsLrKaxYDKYAB4OXWuPlnBvhnJt/8GCLuF0RERFpxgIpBv8IXOOc2wy8DHwBvFbh/PnA9iDm1mhpnUERERFpLgIpBv8XmAoUAuuBa8uWmXHOtcObafxK0DMUERERkZCp8QQS83ZQn156nHzuAHpeUERERKTJqXHPoHPuTefcZdWcv9Q592Zw0hIRERGR+hDIMPEQoGM15zsAg+uUTQNxzv2fc26fc+6Qcy7bOTfbORfV0HmJiIiIhFow9yZujfc8YVP0BNDHzFoB55Ye/9OwKYmIiIiEXrXPDDrnzgH6V2i62Dnn7z1tgduBrUHMrd6YWcW8HVAC9GqgdERERETqzakmkFyLN4MYwPB2GKlql5HDwF1Byssv59wk4DwgHegO7DGz06uIDQN+hpfv6UAu8CwwxcyO+on/JXAfEA8cQFvriYiISAtwqmJwIfAWXm/Zm8AM4I2TYgw4Amw1s4Ig53eyGUAe8AHesHR1HsMrTv8OzAT6lr4e4Jy73MxKKgab2a+BXzvn+gI34+2oIiIiItKsVVsMmtkeYA+Ac+6HwNtmtqu694RYDzPbWZrPZiDBX5Bz7kzgTuAFM7u+QvsuYA7wPeAv/t5rZtuccx8Ci4FLg5u+iIiISOMSyASSp/GGT/1yzrWq4nnCoCkrBGvg+3i9mbNPal8AHANGneL9kUDvwLITERERaXoCKQZnAmurOb8GeKRu6QRNJt4kkNUVG0uHsTeUngfAOZfknLvVOdfaec7Be3aw4lZ7IiIiIs1SIMXglcDz1Zx/Hhhet3SCJgXYb2b+lrrJBpIrrCNoeD2FO/EmwbyIt63enfWRqIiIiEhDCmRYtwvwaTXnd5bGNAZxVL3mYUGFmONmdgi4vKYXds6NB8YDdO3atS45ioiIiDS4QHoGjwOdqznfCW9otjE4BkRXcS6mQkzAzGy+mWWYWUb79u1rlZyIiIhIYxFIMbgB+K6/bdqcc5HAjcDGYCVWRzl4Q8H+CsJUvCHk47W9uHNuhHNu/sGDB2udoIiIiEhjEEgx+ARwJrDMOZfhnItyzkU65zKAZUC/0pjGYA3evQ2s2Oici8HbUaW6iTCnZGZLzWx8UlJSXS4jIiIi0uBqXAya2fPAw8BlwPt4w6zHSr+/HPiNmT0TiiRr4Rm8iSF3n9Q+Du9ZwafrcnH1DIqIiEhz4cwssDc4l4k3+7ZnadPHwF/MbE2Qc/P32T8AupW+vBOIwlvyBryt6RZXiH0cuANvB5JX+HoHkneAoSfvQFIbGRkZtnZtnToZRUREROqFc26dmWWc3B7wItGlRV/IC78qjAUGn9Q2vfTrCrxdQ8rcDezGm/l7NbAfeBxvb+LGMtFFREREpEGFdMeQYDOzIQHEFuP1Gs48VWygnHMjgBE9e/Y8ZayIiIhIYxZQMVi63dxI4HygDZWfOTQzGxuk3BotM1sKLM3IyBjX0LmIiIiI1EWNi0HnXFtgOXAW3r6/VvqVCt8b3lCuiIiIiDQBgSwt8yDQB/gR0AOv+LsSb2LGX/GeI2wX7AQbI80mFhERkeYikGLwamCRmf0ROFTaVmxmH5nZKCAfb+mZZk/rDIqIiEhzEUgx2ImvZxEXlX6NqXD+ReDbwUhKREREROpHIMVgHhBf+v1h4ATQpcL5E3iTSpo9DROLiIhIcxFIMfgx3pZzlK7Ttx641TkX7ZyLA24BdgY/xcZHw8QiIiLSXARSDL4OfMc5F136ehbeEjN5wBdABvBYcNMTERERkVAKZJ3BGcCjZlYIYGbPOueK8LamKwb+1oj2JhYRERGRGqhxMWjeJsaFJ7W9ALwQ7KREREREpH7UeJjYOfemc+6yas5f6px7MzhpNW6aQCIiIiLNRSDPDA4BOlZzvgMwuE7ZNBGaQCIiIiLNRSDF4Km05qRhZBERERFp3Kp9ZtA5dw7Qv0LTxc45f+9pC9wObA1ibiIiIiISYqeaQHItMLX0ewN+XHr4cxi4K0h5iYiIiEg9OFUxuBB4C3DAm3jLy7xxUowBR4CtZlYQ5PwaJefcCGBEz549GzoVERERkTpx3ooxNQh0bjSwwsx2hzSjJiQjI8PWrl3b0GmIiIiInJJzbp2ZZZzcHsg6g38KbkoiIiIi0tAC2YEE51w8cBPQC2iHN3xckZnZ2CDlJiIiIiIhVuNi0Dk3EHgZSK4mzAAVgyIiIiJNRCDrDM4CooDvAslmFubnCA9NmiIiIiISCoEME6cDM8zsb6FKRkRERETqVyA9g4eAA6FKRERERETqXyDF4AvAlaFKpClxzo1wzs0/ePBgQ6ciIiIiUieBFIP3Ah2cc48753o4506eSdximNlSMxuflJTU0KmIiIiI1Ekgzwx+hTdbeCDePsT4qQfNzAJarqZJy89r6AxERERE6iSQwm0RXjEoZb7cC+v/DANGNXQmIiIiIrUSyA4kt4Ywj6YpOgH+8VMoKYL0Wxs6GxEREZGABfLMYLPknIt2zi1wzu10zh12zn3snLuzRm9u2wN6Xg5LfwZr/hDiTEVERESCL+Dn+5xzlwLXAmmlTTuBv5vZ8mAmVo8igP8C38S7l3OA15xzn5vZs9W+0zn43l/g2Vtg2UQ4/WJo3zv0GYuIiIgESY17Bp1zYc65xcC/gDuAYaXHHcC/nHOLmuIMYzM7amb3m9knZlZiZhuAl4BvnOq9xSUGEdEw7GGvYc87oU1WREREJMgCGSb+OXAz8DegPxBbevQHni09NzHYCVbknJvknHuudEjXnHO7q4kNc85NcM5td84VOOf2OedmOufiT/EZkcDFwMZT5bP7wFGOFhZBm+4Q1w6y1gZ8TyIiIiINKZBi8FbgdTO70cw2mtmJ0mOjmX0feAMYE5IsvzYDGAp8Cnx5itjH8PZT3grcCTwH3AUsdc5Vd99PAIfxZk9X69jxYsYvXktBUQmclglZa2pwCyIiIiKNRyDFYBqwtJrzS/n6OcJQ6WFm7czsCiCnqiDn3Jl4BeALZnadmS0ws4l4PZeXAt+r4n2zgAuA4WZ2/FTJnNYmlnc+OcCdf11PcUo67P8I8k9Vo4qIiIg0HoEUg0eBjtWc71QaEzJmtrOGod8HHDD7pPYFwDGg0sKAzrnZwBXAZWa2vyYf0iYuimnXnMkbWz/n95+29Rqz19UwRREREZGGF0gxuBK4o7TXzYdzrh/wU+DtYCVWR5lACbC6YqOZFQAbSs+Xc87NAS4HhppZbiAfdMsFp3P35b34vx2tMZyeGxQREZEmJZBicAoQDawvncTxq9Ljb8B6IBKYGookayEF2G9mhX7OZQPJzrkoAOdcN7wh5Z7ALufckdLjnzX9sJsGduUosXyZ0FPPDYqIiEiTEsgOJJucc4OB/wWuLz3KvAv8zMw2BTm/2ooD/BWCAAUVYo6b2R68IeUacc6NB8YDdO3aFYAOrWLonBTD9vAzuDDrP1BSAmEtfj1vERERaQICqljMbK2ZXYT37OCg0qOjmX3DzBrTw3LH8Hox/YmpEBMwM5tvZhlmltG+ffvy9nNOS+I/Bd2h4Cs48EltLi0iIiJS72rVfWVmuWa2uvQI6Bm7epKDNxTsryBMxRtCPuVs4ao450Y45+YfPHiwvO3cLq15/VAX74WGikVERKSJCGQHkhudc1Wuveec+5Nz7jvBSavO1uDd28CKjc65GLxFsus0y8PMlprZ+KSkpPK2c09rzaeWQlFkoopBERERaTIC6Rm8A2+GblWK8SZiNAbPAAbcfVL7OLxnBZ+uy8X99QyefVoSRhg5CWdqRrGIiIg0GTWeQAL0xduKrirrgRF1S6d6zrkfAN1KX7YHopxz95W+3mNmi6F8sstcvKVwXgBewcv/LmAF8Je65GFmS4GlGRkZ48raWsVEktY+ng+tF12/eBoKj0B0Ql0+RkRERCTkAikG4/F6/6piQGLd0jmlscDgk9qml35dASyu0H43sBtv5u/VwH7gcWCKmVXXw1lr/U9rzb8+7soIK4Gc9dD94lB8jIiIiEjQBDJMvAv4RjXnvwHsrVs61TOzIWbmqjiGnBRbbGYzzewMM4s2s1Qzm2hmR+qah79hYvBmFK846i03Q9ZqP+8UERERaVwCKQb/DtzgnBt78gnn3BjgBuCFYCXWmPmbQALejOKvSORowul6blBERESahECGiX8NXAPMd85NwNvWDeBcoB/wETAjuOk1LX07tyIizLEz9kzO3vMOFB2HiKiGTktERESkSjXuGTSzw8BFwDygM3BT6ZEC/A640MwOhSLJxqaqYeKYyHD6dm7FP0sGeYtP73i9gTIUERERqZlAdyA5aGa3A8l4u5B0BJLN7A4z+yoUCTZGVQ0Tg/fc4F/298QSOsKGOk1aFhEREQm52u5AYqW7kOSamQU7qabs3C6t+arQ+KrXdbDjNTjSGDdoEREREfHUqhiUqp17WmsA1iYNg5Ii2PRsA2ckIiIiUjUVg7VQ1TODAD07JBAXFc47h9pDynmw/mlQ56mIiIg0UioGa6G6ZwbDwxxnpSaxft9XMOBm+GIL/HdjA2QpIiIicmoqBkNgUFo7NmV9RU6XqyE8ShNJREREpNFSMRgC3804DQP+8uEh6HM1bHzWW3NQREREpJFRMVgL1T0zCHBamziGntGBJWv2ceLs70N+njezWERERKSRCWQHEpxz8XgLTfcC2gHupBAzs0rb1TU3ZrYUWJqRkTGuqphRg7rx74VreL3gbK5OTIF35kCfb4E7+UcmIiIi0nBqXAw65wYCL+MtOF0VA5p9MVgTl/Ruz2ltYlm8Opurh94H/7gdNj4D536voVMTERERKRfIMPEsIAr4Lt6uI2F+jvDQpNn0hIc5bjq/K+/tzOOTlG9Bajq8MQUKWsSOfSIiItJEBFIMpgMzzexvZpYXqoSak+9mdCEy3PHn97Ng+G/hyOfw9m8bOi0RERGRcoEUg4eAA6FKpDlKTohm+Fmdef6DLI51OBcGjIL3fgf7dzR0aiIiIiJAYMXgC8CVoUqkKTnVbOKKRg3qxuGCIpZ+mAOXTYXIWPjnvdqVRERERBqFQIrBe4EOzrnHnXM9nGu502Kr24HkZJmnt6FPp0Tmv72T4rj2MOSX8Om/YfPz9ZCpiIiISPUCKQa/AgYCtwMfA0XOueKTjqKQZNmEOee4c2gvPs09yrJNn8HA8XDaQFj6M8j9qKHTExERkRYukHUGF+EtHSMBGn5WJ3p3TODxf+/gW2d3JuyGhTDvEnjmBzDuTYhOaOgURUREpIWqcTFoZreGMI9mLSzM6x2886/reWXzZ3zrnFT4zlOweCQsvQuu/4MWoxYREZEGoe3o6slVZ3emZ4cE5vx7ByUlBmmDYej93rOD789r6PRERESkhapVMeicS3DOneac63ryEewEm4vwMMddl/Xi48+P8M/N//UaL7obzrgKXp8Mu//TsAmKiIhIixRQMeic+55zbjNwENgD7PJzSBWuPrszPdrHf907GBYG1/4e2nT3nh/8cndDpygiIiItTI2LQefcSOAveM8ZzgMc8FfgOeAEsA6YFoIcG51A1hmsqKx38KPPD/Pcun1eY0wS3PQMWDH89ftQeDgEGYuIiIj4F0jP4D3ANqA/MKW07Skz+x6QAZwBbAhueo1TIOsMnuxb56RwQVo7pr60he3/Ld2nuF0PuGGht9TMCz+GkpLgJiwiIiJShUCKwXOAP5lZAVBWrYQDmNlmYD4wKbjpNT/hYY7//X5/EmMiuf3PH3C44IR3osdQuHIGfLQM3rhfO5SIiIhIvQikGAzn672J80u/Vuwa+wg4KxhJNXcdEmN4/PsD2H3gKL98YRNWVvid/2PIHAernoC/jYET+dVfSERERKSOAikGs4BuAGaWD3wBpFc4fwZwNHip1R/n3Hedc/9xzh1xzu2uj88clNaOe648g2UbP2PRqj1licBVv4XLfwVb/g5/vAoO/7c+0hEREZEWKpBi8F3g8gqvXwLuds5Ncc49APwUeCt4qdWrL4EngMn1+aE/uaQHl/XpwPSXt7Jo1W6vh9A5+Mbd8L2nvWcIFwyFXSvrMy0RERFpQQIpBv8PeMs5F1v6ejLe0PADeBNKPsWbZNLkmNkbZrYEb7mcehMW5pj9vf4M7t2eKf/Ywv/720YKThR7J/tcDWNeBRcOf/oWLL4OclrE/BwRERGpRzUuBs1sjZn9T+kQMWaWa2b98WYXnw2ca2b7QpQnAM65Sc6555xzO51zVt2QrnMuzDk3wTm33TlX4Jzb55yb6ZyLD2WOgUqMiWTBLRncNbQnz63L4rvzVpHzVemzgp3PgTtWwxXTIecDmD8Ynh0N+3c0bNIiIiLSbNR5Ozoz22hmW8ysPtZDmQEMxeuF/PIUsY8Bs4CtwJ146yHeBSx1zjWqbfjCwhwTv3kG836Qzs7cowyb/TbPrNnrDRtHxsJFd8HPPoRL/h/seAPmng//uAMOZjV06iIiItLEBVwUOecucc496Jxb4JzrU9qWUNreOvgp+uhhZu3M7Aogp5ocz8QrAF8ws+vMbIGZTQQmApcC3wtxnrVy5ZmdWHrnN+jTuRX3Pr+Jmxa8z+79pXNyYpJg6GSvKBw4HjY+A3POg9cmQ8Ghhk1cREREmqxAdiAJd849AywH/gcYA6SUni4CXgRuD3qGFZjZzhqGfh9vh5TZJ7UvAI4Bo4KZVzB1T45nybhBzLj2bDZnH+TK2W/zxJs7KCwqfZYwoT0M/zXcuQ7OvgFWzYUnMmHz81qbUERERAIWSM/gvcD1eL1rffGKLQBKF6L+O3BVULOrvUy8hbFXV2wszXND6flypYVuDBDpvXQxzrno+kr2ZGFhjpvO78q/fj6YoX068OjrHzNs9kre/jj366DWXWHkXPjRvyGxo7cu4eJrYf8nDZW2iIiINEGBFIO3AIvM7H+B/X7ObwN6BCWruksB9ptZoZ9z2UCycy6qQtsP8BbSfhboWvr9RyHP8hQ6torhd6PSWfjDTMyMW55azW1/Xvf10DHAaekwbjlc9Shkr4PfXQBvPqQFq0VERKRGAikGTwdWVXP+K6BNnbIJnjjAXyEIUFAhBgAzW2hm7qTjdH9vds6Nd86tdc6tzc3N9RcSdEPO6MBrEy7hnm/2ZvlHX3D5rBXc9+ImvjhUeith4TBwHNyxFs68Ft7+jTfJ5OPX6yU/ERERaboCKQYPA22rOd8TqJ/q6NSOAVUN88ZUiAmYmc03swwzy2jfvn2tkquN6Ihw7hjai7d/cSnfH9iVJav3cclvl/Pb17ZztLDIC0rsCNfNh9FLISIa/nIDPDUMti+DkvqY7C0iIiJNTSDF4H+AUc45d/IJ51wbvAkly4OVWB3l4A0F+ysIU/GGkI/X9uLOuRHOufkHDx6sdYK11aFVDNNHnsW/fz6YK8/sxNzln3L5rBW8sumzr/c47n4J/OQdGPYIHMyGJTfB3ExYvQC+CulSkCIiItLEOKvhDFTnXAZeQbgKWAj8Efg5Xg/bL4EOQKaZbQ1JppXz2Qwk+BvOdc49iLdDyiVmtrJCewxwAHjbzIbXNYeMjAxbu3ZtXS9TJ+v25HH/i1vY+tkhLu6VzPRrzuL05ArrahcXwbZ/wDtz4LPSHUza9YS0S+Hs70DXQQ2TuIiIiNQr59w6M8uo1F7TYrD0IlcDTwIdS5sMb1bxF8AtZlZvD6mdohg8G/gQ+LuZXV+h/U5gDvADM/tzXXNoDMUgQFFxCU+/v5dHX/8IDGbd2J8r+nX0DTKDL7bBzuXw6XLY8w6cOAaZ4+CKX0FUo9qYRURERIIsKMVg6YWigSv4enmZHcBrZlarZ/AC/OwfAN1KX94JRAEzS1/vMbPFFWIfB+7AW/LmldJ87wLeAYbWZccU59wIYETPnj3H7djReLaGy/4qn58sXsem7IPcfXkv7hrai7CwSqP6nuPH4M0H4b3/gzanw7W/Vy+hiIhIMxa0YrAhOefeAgZXcXqFmQ2pEBsO3A2Mx5sJvR94BphiZkeCkU9j6RmsqOBEMf/z90288EE2l/XpwKV9OnCiuIQTxSW0ionkO+mnERFe4VHR3f+BF2/zniXsc7V39LoS4ts13E2IiIhI0DWLYrCxaKw9g2XMjEWr9jD95a0Ulfj+fi/ulczcm8+jVUzk142Fh2HFI7DpeTicAy4MUtMhti2ER0JYBMS3hz5XwekXe20iIiLSpNSqGHTOvRng55iZXRZock1VY+wZrOhg/gkKi4qJCg8jMjyMpR/mcN+Lm+meHM9Tt2bSpW2c7xvMvEkm21+B3Svh+FEoKYaSE96s5BNHIbaN13s4cDx0PrdhbkxEREQCVttisAQ4AdR0GRYzs1a1S7HpaezFoD/vfrKfn/x5HRHhYfzu5vM4P62Gw8En8uHTN2HrP+Cjf0Lxcbj+Seg7IrQJi4iISFDUthgsKwKX4S0l83JdJl40F419mPhUduYeYczCNew+cIyrzu7ExCvOoGeHhJpf4Oh++MuN3vZ3V/3W2/1EREREGrXaFoPt8fYkvhXoh7eEzCLgKTNr8L17G1pT7Bksc6SwiAVv7+TJlTvJP1HMDeldGNTj6w1mYiMjuKxvByLDq1iX/PgxeH4sfPQKXPQzOC3TKw6z10HebmiVAm26eTOVy47W3SCxM4QFsta5iIiIBEOdJ5A45wbi7TJyI9AKWA38AVgSrNm5TU1TLgbLHDhSyNzln/Ln9/ZwvNi30/eS3u353c3nER8d4f/NxUXwyj2w7o/e67AI6HiWt6j1kc/hy91wMAtvOcpS4dHQuqtvkZjSH1IGVL/W4YkCOJoLrVJVTIqIiNRCMNcZjAGuB34IXIq3A8ltwVjEualpDsVgma+OHefLYycAb/HIlZ/s54GXttCvcyueujWT9olVbPVsBrtWQGQ8dDobImN8zxcdh4P7vMLwy93w1Z6vv/9yNxSUbunnwqFjP0g+A0qKvGcSiwrg2AFv8sqx/V5cl0Hemohtu1d/Qwc+9XJr1wMq76AoIiLS4gR9aRnn3EXAVOBy4AEzm1a3FJuOpv7MYE39e9vn/PQvH9AhMYY/jRlI9+QQ7FJydD9kfwBZayB7LeTthPAoiIj2ehFjW3u9gUmneT2P/3kMrASGPQwDfuBb6Jl5k1xWPeF9BW/282mZXx+p6RDTYuY4iYiIlAvWdnSdgdF4zxD2AnLwniH8vZntC06qTUdz6hmsyvq9XzJm4RqOHi9mQJfWDEprx6C0dnRrF4dz4HCEOWgTH1X184XB9NU+b5Hs3SshbYg3JG0GGOx9H77YAgkdvaVv4tt7RWbWGsjdXnoBBx36Qo+hMPQ+iIwNfc4iIiKNQK2LQedcJHAN3rDwN4Fi4CW82cWvteTZxS2hGATYe+AYi1bt5v1deWzJOUiJnz8yYQ46toohpXUsp7eL59v9U7i4Z3LV2+HVRUkJvP87ePdxbzgZAOdNWjn/J3D2d7yexYryv/Imt2StgX2rvZ7D078B318C0QHMpBYREWmiajubeA5wE9AG2AQ8BfzZzPJClWhT0lKKwYoO5p9g3Z48vjhUiOF1yhWbkXuogOyvCsj+6hjbPjvMwfwTpLaO5cbMLtyQcRqdkxpZD9yHz8CLP/GGjm9+DmKSvLUUVy+A1fPhtAwY/ltIaF/9dT7fCjnrvcktR3Oh8BAMuAW6ZNbPfYiIiNRQXRadzgf+DnxQg88xM3us1lk2ES3lmcHaKiwq5vUtn7NkzV7e+eQAYQ4uPaMD3xvYlUvPaE9EeBiHCk6wLecQew4co19KK/p1bhWaXsTqbHnRWx6n41nQ/2b4zyw4/Bl0Od8r8KISvHUUz7q+8iSUY3nw72mwbiHls6UjYiEs3CsqL/kFXHKPtu4TEZFGoy7FYCDMzMIDTa6paok9g4Hac+Aoz6zZx3Prssg9XEiHxGhiIsPZm3fMJy45IYpv9EzmnNNaU1BUzNHCIo4UFNE9OZ4bMrpUWt7mUMEJPtjzJSmtY+nWLo7oiFr+sfv4NXjmB1BcCKcNhMunesPHX2yHf/zUm9TS60rocam3TmKbbpC1Fv71gDcT+vyfQOZY7znFqHivZ/CV/wcbl3iTVb79BEQnQv6XUPCVNxmmXY9a/jRFRERqr7bF4OBAP8jMVgT6nqZKxWDNnSguYfn2L3jhg2zCw1x5b2CXtnF8uO8r3t6Ry8od+8k76j0DGB7miIsK53BBEa3jIhl9wenceuHp7DpwlL++v5elG3MoOFFSHtu1bRxndEwkvVsb0k9vw1kpSURF1HBCS/YHkJ8HPS7z7QEsKYZVc+HtR6HwoO97ul4IVz8KHc/0f83NL8DLE7wCsCIXBueNhkv/BxI61Cw/ERGRIAj60jKiYjDYSkqMA0ePkxAdQUxkGM451u35kt+99Sn/2vY54WGO4hIjLiqca/qnMPyszuQdPc4nXxzh09wjbM45yL68fID/z955h0dVpQ3896aTCiEFCCVA6L2DCIILKKKAWBAVu7KKurb1c3fdtW+x94KKIhaKvSKC9CK9d0IChATSe5053x/nDkz6JCQkIef3PPPczLnn3nvumynvvBVvDzf6tW3KhVEhXBAVQu+IINJyC9kbn8He+Ayy823cMLQtYQE+lawKHRiZkwypsbpOopcfdBpXef3C9DjY+4POWG7STMcl7v8ZNn6gXcojHoIul+kEFi9/nfSSk6ILdmed1PPbXVADkjUYDAaDwSiDtYJRBs8d+iXyEQAAIABJREFUB05mMn/jMaLC/LmiTyv8y+mKciojj82xqWyKTWXd4WT2xGcA4OXuVqzDigg08XTnrpEduHNEB/y8PYhOzOKbrXH8ujuB7i0DmTk6ik7hAaePiU3O5sPVRziemsu1A1sztnsL3KsT55h0EH77l1YMK6PfjTD++Yq7sxgMBoPB4AJGGawFjDJY/0nOymdddDLbjqbRIsiH7i0D6doykPTcQl74dR8/70wgxN+b1s2asO1YGm4CA9sFs+tEOrmFNsb3bMHkvhF8t/0Ev+yMx8PNjWA/LxIy8mgT3IRbLmjP1f1bE+RbPFFEKcX+k5kcOpXFxV3D8PUqQ3mN2wKpRyA/CwqyoSgXfJvr+EP/MNj3M6x6CUI6wdUfQYuert10xgldsNuveeVzj67XNRj73WTa/BkMBsN5jlEGaxCTTXz+sOVoKi/+up+0nEIm9W3F5H4RhAf6kJJdwOzVR5izNobM/CICvD24YWg7bh0eSXM/L37bc5IPVh9hc2wqbgJ92jRlRKdQBrZrxubYVH7aGc+hU7pld7CfF7cNj2T6sEiCmlQxuzh6OXx9l66TOOgOaD0AWvaFZu3BXghpR7X7OukAHN8AxzZCxnHdHnDy29BjculzKgVHVsLKF3TxboBBd+rM6fJc39nJELNSX2vALbozjMFgMBgaFEYZrAWMZfD8Jz23kA1HUhjaIZgAn9KK3I7jaSzZc5KVB5PYcTwNu9L61JD2wUzo3Yp2wb58tOYIy/Yn4u/twbju4YQF+hAa4E1ogDcjokJo5udV8SKyEuHHB+Dg4jNFtj2a6N7NOL1/g9qcabu3+xutHF74kO604uYOhXk6hnHDLL3PvwUMv19bEte9CcPuhXHPnlEIs07puQcXQ/yOM9dq1h6mztW9qB1knNCtAkM6a6XV9IM2GAyGeodRBmsBowwanEnLKWDbsTS6twoslZiy+0Q6766IZlNMCklZ+RTa9PsuNMCbV6f2ZXhUyOm5+UU25qyNYX9CFneMaE+3llYv5aIC7dKN34795G7y3API9o0gwyeCVJ825PuEYFcKpaBtkDuRG57SdRCjxkBwR9gxX2c3N4vUil+/6eDpoy2Fix6DP97VyuPQu2HNa7DxQ11yp+0FuvVfh1FgL4Ivb9Wlci5/FbpdDqtf1VnXDuW083htlfQNrn2hGwwGg8FljDJYCxhl0FAdlFKk5xZy8FQWf/t6J4cTs5gxsiMPj+vMsn2neO7nvcQm5+Dj6UZ+kZ2r+7fm4XFdCPbzYsWBRL7dFseSPSfJLyq/DKi7mzBjZAcebLYGz1//T1vquk2E/jdB5IjS8YFK6VI4mz8Cdy+t9PWeqotnl6yLmHUKvrxNu5i9A3VtxZ5Xw5/+CfsXwW//BN8QuOoDiBxeCxI0GAwGQ3UwymAtYJRBw9mSW2Dj6R/38MWGozT38yI5u4BOYf788/Lu9GndlDeXHWTO2ljc3MDbw5303EKC/by4rFcLokL9CWziSYCPJ37e7riLnO7issAq9N0hxI8XL2tJzzYhRGd7sj8hk2MpOUSF+dO/XbPiFky7HZb8C3JS4cIHISSq/IXbimD5fyB+O4z+my6w7eDENm09TI3RmdCD76wd4RkqpiAHEnZCm8HGbW8wGACjDNYKRhk01BQ/74znrWWHuGZAa24Y2g5P9zOWu2MpObz5+yEKbHYm9mnFhZ1Ciu0vj9UHk/jbNzs4lpKLh5tQZC/9Xm8T3IRuLQLx8nDD3U1wdxP6tG7K9UPaunSNcsnP1Ikv+3+GEY/ouMWKFBKljMJSk9iK4PNr4fBSHSYw4WXdPcdgMDRqjDJYCxhl0FDfySkoYvbqI+QW2ugcHkCXFgG0bubLgZOZbIlNZXNsKodOZWGzK2xKkV9oJyEjj05h/jw5sUexWMYqYyuCnx6ELZ/o+MTLX9UZ0DFr4NASSNyrk2OyT+mi3j5NIbiDfoR20RZFnyDXr2e3Q0GWTpY513UZC7J1HGVgRP1Qan95DP54B/pM00lDyq673gy5G9zLrtFpMBjOf4wyWAsYZdBwvqGUYsneUzzz4x6OpuQwvmcLxvdqScsgH1oE+hAe6ON6mz99Qlj2nC5jE9IZ0o7peooePhDe06qpGKrrK+akQEo0pByB9GM6O3rKLGg37Mz5Di3VfaFTosHNA9w9QdyhMEdbI1H63BPfhN7X1LR4SpMRrxNvNn2kWxb6NodW/fSj20Ro2bvyc6THnels4+kHXr5gK9T3VJCjE3PcPa2Hl1aQK1I6N38MP/wFht4Dl/4H0o/DTw/DgUVaplF/go4XQ/uRujOOwWBoNBhlsAYxdQYN5zt5hTY+WBXNW8sOk1toOz3u6S4MbBfM6K6hjOoSRqcwf8QVS9im2dpC2HowdBoLkRfqNn3lcXwTfHWHVpJGPAw9pmgl8OCvOhu6y2U6ycXx8PQD7wDwCYT9v0DsGhj9OIx8pHylKXoFbPpQWyfzMyAvQ2dPezY5o5T5hULTdtrFGtQGlA3y0vUjYRfs+kqPdZsI7YZDwnYdM3lqrzV+BYz6W9k9rNPjYNWLsGWutphWBU8/XYw8tCuEd4cWvaFlHzi1Bz6ZBO0vgusXnLECKgX7foTt83SNyfwM3Sc7rAe0HqgfbYZWHCdqMBgaPEYZrAWMZdBwvpNTUERcai7x6XkkpOdxKDGLlQcS2ZeQCUDbYF+u7BfBVf1b07a5b81ePD8Tfvk/2PaZfu4dqLObh8zQfZzLoygfvr8fdsyDvjdo97SHUy3H+O2w5CkdT+cfri2W3oFakXT3gsJcyyqXpTOnU2P085J4+kG/G7QFLrh98X25qbD+HVj3NhRkaqUwvKdWWL38tbK4abZ23/a/CbpO0BbAgmz9cPfSSqmXn75Xe5G2FtoKIDtJFxlP3K8fmSfOXFfcdfb3HUvKd7HbCiFuMxxeputNHt+srZoA3SfB2GdMfKHBcJ5ilMFawCiDhsbKibRclu9P5Oed8aw5nIRSMLh9MN1aBFBkV9jsiiK7IiuviMz8QrLyisgvstPEyx0/Lw/8vN2JbO7H6K5hDGjXrOJklT3faeVl2H3apewKSsGK/+mMZ/9wXerG218rS0fXavfoiId15xVPn8rPlZ2kXdcON61PkFbqKmvhl5OiazBu+lAriA7EHfpOg5GPnr3ilZMCCTt0YfCMEzD0z9p66ip2OyQf1IXKV7+qFdTh9+uMcldjL+02HatpMBjqNUYZrAWMMmgwQFxaLt9ujeObrXEkZeXjLjor2cNN8PP2IMDHgwAfT7w93MgttJGdX0R2vo0jSdkU2OwE+HhwUedQHhjTiaiwgJpd3J7vdI/ngixtaSzI1rFyw/9y7lvq2Yr0OgqydFyj31kk59QW6ce1O37nQqeYzQvKn59yRFtvDy3R7v/+N0OncSZJpTqkHNG9yjteXNcrMZzHGGWwFjDKoMFQfbLyi1h9MInf953k190naeLpzrczh9MiqBJLnaH2iV0H396tYzZHPqrd884KXlE+rHldxzy6eege2Ad/g6yTus1hm0G6n3ZOiu5602E0jHnSdcvu+UBOCqx/W/cPD+kMYd10jGdwh9JxrIkH4OPLdFb9fZv1HIOhFjDKYDmIiAfwEjAdcAO+AmYqpfIqO9YogwZDzbA3PoOr31lLZIgfC2YMw89bKx4FRXZeW3qA1QeT6NIigJ4RQfRoFYSPpxtJWQUkZeaTmlNAgI+H7vfs70N4kDeh/t6uJbYYyic/E37+K2z/QieX9LlOW65SjsCJrdpt3n0yXPJvCIrQsYgHF+tEoZQjOrPaN1hnQe/9ETx9db3JgbedveUw+TCsfFHHWU56Syf71BdyU3VowPp3tRU4oGXxuM72F8GkN6FpW/08NQZmj9dJRHnp0O9GuPyVss9tK9TKZfJhyIjTMaSOeNKg1jqxqj7JwlDvMMpgOYjIv4CrgfFAAfA9sFEpdX9lxxpl0GCoOZbtP8XtH29kdJcwZt00kCNJ2Twwfyu74jLo26YpR1NySMkucOlc/t4edAj1o0OIHwPaNWNyvwgCfDxr+Q7OU3Ys1K0KCzJ1zGSzSN3revAduqC1KyQegF/+CtHLoXkn7aLPTtIPW76O6fQPBb8waD8CBt9VdrZ5agyseEErqO5e+tjIC2Ha/HOjBGUnw+6vdcJQyXu323Vs6NJndEJO98kw6jFtEczL0Ek/sWt1LCvAJc/pc3w0XiveN/8IG2bpjO8Hd4F/2JlzZybA51N1Rxllo1y8AqDHJJ041ap/5fGwhkaHUQbLQUSOAo8qpeZZzy8BFgLNlKroXWeUQYOhppm7PpZ/fruLEZ1C2BiTQhNPd/53VW/G9WiBUoqEjDx2xWVQZLMTGuBNiL83zfy8yMwrJDEzn8TMfOLT84hOzCI6KZtDp7KIT8/Dz8udqwe05qYLIukY6l/Xt9nwyE3VCk1Q6+oniiilYzj/eE9nSPuFaCXQw0srWdmndALMyV0Q0EoXye57vbZ+7f8Fdn6pSwuJOwy6HYY/AEdWwDczdFmf6+dXnvBit0PMSq1wpcbqa7t7a8Wzy3jdY7ssq2XcZtjwgS4lZMvXY10vh0v/C03b6FjL72ZqZbfDaBj3DLToVfYa0o7quUdWgkcTbTm96TuI6A9Jh+DNgTp5Z8wTZ+T22TW6F/iwmdA8SruRg9poObp56Ef8dtj2Oez5VlskQWe8+wZDQAsYcKu27lb0/1NKlycK6azXZTjvaPDKoIj8DegPDADaA7FKqchy5roBfwFmAJFAIrAA+JdSKttpXlMgFeimlNpnjYUCp4AopdThitZklEGDoeZ55sc9fLj6CCM6hfDiNX0IDzw768a2Y2nMWRvDjztOUGhTXNGnFY9P6HbW5zXUEjFr4Ld/agWsWXttPSzI1O7WXlfD0JkQ2PLM/B0LKlcIs5N1R5bt87R72zsIWvTUimZRvlZ204/p6418BHpdCyd3wr6fdAJS4l6tWPW5DgbcossSrXhen7vvDbBjvs6ovuQ5vb+yEAW7HTbPho2z4fKXoe3QM/sW3ASHl2vroE8gbPwQfnoIxr8AQ+6qXH4F2Vp5To3R95WTrC2KJ3dBSBftqu92RfE1phzR97D9C31c1BiY+mnFtUANDZLzQRlUQAqwBa0QZlSgDL4G3A98A/wCdAPuA1YBY5RSdmteG+Ao0FIplWCNeaLdxf2UUtsqWpNRBg2GmsduV+yJz6B7y0Dc3Gou7i8xM585a2OYtSoaL3c3HhzbmZuHtSO7wMaKA4n8tuckSZn5TO7Xiiv6tMLXy2TE1hmnrYjvaktYr2u0O7g8q9aOhfDNXdo1Ou2L4i7W1BiYe6XedhitrY1dJxRXdOx2OPCLduHGb9fWQlu+tkK2u0DXX+w9VStnDtKO6rZ/+3+CthfA5LdL15usDnFb4P3RMPZpbX1890KtLN7wVeWljMpDKdj7Pfz+rHZXB0Zoq6KyayU2/RggOtO+ZW9Y+6Z210+bd3atHW2FxsJYzzgflMEOSqlo6+9dgH9ZyqCI9AB2At8opa5yGr8PeB24QSn1uTXmsAx2VUrtt8aMZdBgOI+JTc7mX9/tZsWBRFoG+ZCYmU+RXdHcz4sgX0+iE7MJ8PZgcr8IurQIICE9j/j0PBKz8gkP8KZLiwC6tgike6tAgv28Kr+g4dyw90fdtcYvRFsIw3toi9inV2nr3/ULoO2Qis+hlM6K3v+zVsA6jdNu1opIOqSVwJqsszhnoi4oHtQakg/BPesgsNXZn9dWpC2Ah5boDjTiptcd0hl6X6uvB9qC+u3dOnHohgV6XvRy3dIwL0O3NOw0Trufy+O3J2Ddm9pi2/lS6HyJjvOMWaVd5Mc3wZ/+qRXtmqQoX587brN+JOyEEQ/pxJySHFislfkRD59J6KkOBTkw/wbdKan/9Oqf5xzQ4JVBZypRBp8F/gGMVEqtchr3AZKBFUqpy5zGjwJ/VUrNt56PA77ExAwaDOctSil+3Z3AZ38cpXurQMZ2C6df22a4CWyKTeXzP47y0854CorsuLsJYVZ8Ynx6LklZOonFw034x4Ru3HJBZJUyl3/ccYKVBxJ59NKuhPhX0EnFUHVObIXPr9Ou0pGPwKqXdNeXG7+GsK51vTrXOfy7tmYCXD0bel5V8fzaYNfXWrkOaAnZidpS6h2oLYWZ8XpOi95aziUVuh0L4es7IHKEdvMn7i2+v0nwmV7b922uuKNQVSjI0SV6TmzVz5u20+7wnBS4dxMEhJ+Zm52s4zNzU3Tdzwsf1PVHq+MaX/w4rH1DJ/Dcv6W4ZbqmUErHgy7/H/S6Spd7qgaNSRn8FRgD+Cql8kvsWwN0VkqFOo39C5gCXAYUAt8Bm0w2scHQuMnIKyQn30aIvxceTh1SkrLyOZCQyew1MSzZe5JJfVvxnym98PXyQCnFtmNpfLM1jpZBTbh+SFuCmmg3WU5BEU99v4f5m44B0DLIh3duHEDfNue4+PX5TnocfGFl3oZ0gelfn7F4NRSUgi+u00kiE16su3Xs+xnWvg6t+mnrXtthWok7uVuXEdr1lY5FHPOkTugR0Z1wPhynj7n5ez0/NUZbXO027X4O7aaTf+ZO1kk4Q+8++7UqBV/epjvpXPGqdrH7hWjL7TvDoMeVuoi6g2/v0VbSafN04s3uryGoLfSZqnuS+zY/k+TkF6Kfl+XyjtsMH4yBqLE6lrTfdH39muToH1rhPL5Bx67a8uHutRDapfi8hJ2w+1udYFWOJbkxKYM7gTClVHgZ+xYA1wDeSqkCa8wDeJkzdQa/BO5VSuVWtg6jDBoMjRe7XfH28kO89NsBuoQHcN2gNizcfJzdJzLw9nAjv8iOn5c70wa3ZVSXMJ78YTeHE7OYOSqKsd3Dmfn5Fk5l5PPUpB5MG3wWLipDafKzdDJEz6sqd/Maqk9RvnYn7/pK148c9Td4/0+6/M1dyyu3kM25Ak7ugb9s0xbcs2Hli/D7M1oxvfDB4vt+fxZWvgC3/KRjT4+sgjmXw4UPncnaPrIKFv9Dx4yWR1BbnfDTaax+XlQAsy7SBdZnrodl/9blgf68WocpuEJhLmx4X8exNu9YfJ9SsOhvOvnJv4VO/uk0Dt4aBC376ix0h1ciOxneGwkZx3XM6+A7tRxKdDpqTMrgYcBTKVXq01VEPkErfc2UUmnVvPZdwF0Abdu2HRAbG1ud0xgMhvOEFQcS+cu8raTlFNK1RQA3Dm3H5H4RxCZn8/7KaH7YEY/Nrgjx9+bVqX25sJP+cE7LKeD+edtYeSCR6wa14elJPfHyKJ4gMH/jUWatjGZw++ZcO7A1fds0NcW0DfULux1+fxpWv6LdyEX5cNsvEDGg8mOPb4YPLobR/4CLHj0zfnCJzuYecKtrbSP3/QzzpulEoynvl87mLsiBt4doq9odS7QCZy+Cu9eVrk9pK7KysK06mDnJ1t/JOqnp1G4Ydi/86QkdhrDivzoetfMl2h39ej9tFZ3+TeVZ5amxMP9G3VvcLxSmf6uz3EErgr/+A9a/pXuoj3lS91cHrTz+/Ahc/RH0nKL/B59dreMxr/5Ix7xu/0IXeh/xsHZ/WzGtjUkZrJJl8GwwlkGDwQDadXwyI4/uLQNLKWvHU3NYeSCJsd3DCQ0oHhtlsyte/m0/by07zOD2wbx74wCC/bxQSvHS4gO8uewQncL8OZaaQ16hnagwfy7r1ZIu4QF0DPMjsrkfPp41mLhgMFSXTbO18jLhZeg7zfXj5t2gE0ru36aVnd+e0JYw0CWAht6tH55NdBLLnu8hepmO8/MP19bHQ0sgpBPc+kv5MX/7f9Gu97AeWqG78SvXi6Y7KMzT7tqN70N4Lx0L2WMKXPX+mTnr34FFj51REMsjejksvFW7zsdYimVBlo5vbT0Qlj6tx4b8WbvSnT9X7DaYNUorq/du1B1vlv9bd64ZeJuek7hfn2Pfjzp2c8osCGzVqJTBKsUMVvP6VwBXREVF3Xnw4MGzOZXBYDDw3bY4/vrlDloE+vDOjf35YNURvtkax9SBbXj2yp7kFdr4aUc8CzcfZ3Ns6unjROCGIW15emLPGi3DYzBUC1tR1VsNntqnY/p6Xg2J+7SVbPAMXcpn9ctamfEOAhTkZ+i/oy7WGc6ZJ3U/bO8AuO6zyjOuv5imrWY9r4arP6z2bbL3B/juXh1DeM8f4Nf8zD5bIbw9FBCYsbK05bEwVyebLP+PzuK+7nPtHk6NhU8magWvx2TY+in0vxmueK1sC+OxDfDhWOh4MRxepuV15bvF5yoF2z7TbSU9fOCCe5GRjzQaZbCybOKVSqnxNbEOYxk0GAw1xdajqdw1dzOJmfo37CPjOjNzdFQpS2NugY3opCwOJ2az+mAiCzYdZ/rQdjw9qYdxIRsaJt/cDds/hybNdK/prhPO7IvfoZUnDy/oNgk6XFT97OP0OK1gXvSYbn94NmQn6aLlZSmgB36Fz6/VSSdD7tatG70DdaLK8v/ovtI9roSJbxSPlcyI10k1ifu0cjf5nYpLFn03UyuNod3gzqXl14RMPKCTa07uRJ7KaDTKYC9gO+XXGZyulPr0LK9vLIMGg6HGOZGWy9M/7GF8rxZM6htR6XylFP/5ZR+zVkZz54j2/P2ybkYhNDQ8Mk/q2LjBdzW8zO/yOLoeVr2sWyh6+WuXdsphHUs55imdVV0WOSk687rnVZVbWbOTdbzmBfeXTj4pSWEebJmDDP1zw1YGRWQ60M56eh/gBbxkPY9VSs11mvsGcC+6A8nP6A4k9wNrgIsdHUjOFmMZNBgMdY1Sin99t5u562O5/0+deGhsZ5eOS8rKJyO3EA83NzzcBX8fDwJ9zr5bxHfb4kjMzOeOER3O+lwGQ4MnYReseVW3/Bt+P3SbWHliSS1SXsxgQ+q3dDtwUYmxZ6ztCmCu0/gDQAw663cCkAS8ge5NXCOKoMFgMNQHRISnJvYgr9DG60sP0sTTnbtHlW8lOJaSw+tLD/L11jhsduV0HpgxsiMPje1cKqvZVZbtO8UD87ehFLQMasKE3i0rP8hgOJ9p0ROu+qCuV1EpDUYZVEqNqsJcG9pq+FJlc6uDk5u4Nk5vMBgMVcLNTfjvVb3JK7Lzv0X7aOLpxi3Di/fJPZGWyxu/H2LhpmO4uwk3D4ukT5sgCm2KIpudTbGpvLviMGsPJ/Hadf1oH1K1nrQHT2Zy/xdb6d4yEE93Nx77ege9WwfRJti38oMNBkOd0mDcxPUR4yY2GAz1iUKbnZmfbWHxnpP8d0ovrhvclsy8Qt5ZfpgPVx9BKbhucBtmjo4iPNCn1PG/7Iznsa93Umizc9/FnejSwp/wQB9aBPoQ7OdVbjxianYBk99eQ3a+je/vHY7NrrjstVVEhfuzYMYwPN2rZ2k0GAw1y3lVWqa+YJRBg8FQ38gvsnHXJ5tZeTCR6UPb8dOOeJKzC7iyXwQPj+tM62YVW+pOpOXy0IJtrI9OKTY+qkso79wwgCZexbMbC212bp69gU0xqXxx11AGtGsG6B7M936+lXtGdeTRSxtQX+BKOJmRR4CPB75eDcaxZjCcxiiDNYjJJjYYDPWZ3AIbt368gfXRKQxpH8zjE7rTq3WQy8crpYhPzyMhI4+T6Xnsjc/gzWWHGBgZzOxbBuHvrRWhkxl5PDBvG+uik3npmj5cNaB4JuhjX+1g/qZjzL5lEKO7VNKazIldcel89sdRftx+gpAAb/q2aUq/tk3p3bop7UP8Tvd7Ptcs3p3A/fO2Etncj/kzhtX4OgqK7Gw9msqaQ0msi06mQ4g/j1zSpVSxcoOhuhhlsBYwlkGDwVBfySu0sT8hk96tg2qk3MwP20/wwPxt9G4dxMe3DmbbsTQemr+NnAIbz0zuydUDSpcEyS2wMeWdtUQnZvHBzQMZ0an82m52u+KHHSf4aE0M246l4ePpxvieLcnOL2LrsbTT9RcBmvp60i7Yl/BAH/y8PWji5Y6/twdX9G5VodJbaLMTm5zDkaRsWgb50DPCdQX54zVHeOrHPXQOCyA6KYt+bZvxyW2Da6QDTFZ+ES8vPsC8jUfJKbDhJtCjVRD7EjLw8XDn4XGduXFoOzzOgbt9V1w67UP88PM+e8tnXqHNdMipZxhlsBYwyqDBYGhMLNqVwH1fbCEswIe4tFy6hAfw1g39iAoLKPeYlOwCrn9/PdFJ2cyaPoBRJSyESimW70/kf4v2sS8hk46hftw4tB1T+rUmyNfz9JwT6XnsPJ7O0ZRsYpNzOJqSQ2JmPjkFNnIKisjILUIE3ry+P2O7n+lGarcrPlkXwyfrY4lNzimWQX1Bx+bcMyqK4VHNERGOpeSw5lASe+IzaN2sCZ3CAogK8+ejNTHMXnOEsd3Def26fizek8Bf5m3jkh7hvH3DANzPovvL4t0JPPH9bhIy8riybwTjeoQzrGMIQU08OZyYxZPf72bVwSS6tQxkSPtgAn08CPDxpGOYH6O7hFWq6CulSMsppJmfV6Vr+XlnPPd8toXwQG/+Nr4bk/q2qvYPiSe/383CTcd4b/rA0/24DXWPUQZrEOMmNhgMjZVl+09x3+dbuaJPK564ortLlp/U7AJu+OAPDp3K4t3p/enWMpAjidkcTsrmh+0n2HAkhXbNfXl4XBcu79WyWq31krPyuW3OJnYeT+PZyb24fkhbEtLz+OuX21l1MInB7YMZHBlMh1A/IkP82ByTyvurojmVmU+3loHkFhQRk5wDgK+XOzkFtmLnv3V4JI9P6H5a8Zu9+ghP/7iHaYPb8u8re5ZSmg4nZvGPb3bi6+XB+zcNLKUwZuQV8teF2/l190m6tgjguSt7nY63dEYpxS+7Enh1yQHi0/LIzC86ve+VqX24sl9pi6zdrth6LJVFuxJYtDuBYym5zL9rKEM6NC8110F0YhYT31xDZIgvbiLsOJ7OwHbNeHJijypZUAEWbDzGo1/tIMDbg/wiO2/dUFxBN9QdRhnKhFelAAAgAElEQVSsBYxl0GAwNEZsdlVla1haTgHTP9zAzrj0YuNhAd7cd3EUUwe1rXZ9Qwc5BUXM/GwLy/YnMqV/BEv3nqKgyM4/L+/OtMFtSils+UU2vtkSx+cbjhLq782FnUIY0SmEjqH+pOcWcuhUFgdPZRHs58UlPVqUut7zi/bx9vLDdGsZyG3DI7miTys83d34cHU0Ly4+gIebkFNg4+GxnbnvT51OH2e3K+78ZBMrDiTyyCVduP3C9i5nXNvtisz8Im7/eCP7T2ay+MGRtAxqcnr/yYw8rn9/PYcTs/Fyd2N4VHM2xqRyWa8WPH91nzLPmVtg48q313AyI48f7x9By0AfFmw6xvO/7ic1p4DJfSN4YEwn2jU/U27oRFoui3Yl0Cncv5j7f+vRVKa+t57B7YN57bq+3PbxRnadyODla/swqW8Edrvi4Kks9iVkMDwqhBB/Ew95LjHKYC1glEGDwWBwnfScQj7bEKvdnCF+dAj1JzzQu0Zb6BXa7Pz9650s3HycPm2a8urUvlWumegqSikWbj7Oh6uOsP9kJiH+XoQF+LAnPoNx3cN59sqePPvjXn7aGc+CGcNOW/5eW3KQV5Yc4OlJPbhpWGS1rh2TlM3411YxMFLHLooIaTkFXPveOuJSc3l6Uk/G9ggn0MeThxds57c9CWx8fAzeHsUtuUopHlm4g6+3HuejWwYVc+On5xTy9vJDzFkXQ5FNMXVQG/q0acq3W+NYF52MQ30Y2TmUv1/WlWA/L654YzVeHm58P/NCmvl5kZlXyO1zNrExRicz7T6RQWaetm62a+7LZ3cMqTTDfVdcOkFNPKtcs7KgyM6MuZuYOqgNl/Y0BdDBKIO1glEGDQaDof6hlGL78XR6tAo8JzUOlVKsPZzM7NVHOHgqi4fHdWZiHx1vl5FXyGWvrQLg57+MYFNMCrfP2cSV/SJ46Zo+Z6UIz10fyz+/3cUzk3sypV8EN3zwB3viM/j41kFc0PFMnN7y/ae45aONvH/TwFLu2vkbj/J/X+2ssJXhyYw83vz9EPM2HqXQpmgb7MuU/hFc3rsVKw4k8vrSg2TmFdIi0IfUnEK+vucCurUMPH18boGNx77ewf6ETPq3a8aAts0IauLJQwu24e/twWd3Di1XYU/JLmD4f3/Hw11454YBVYo//HR9LI9/u4s2wU1Y9vCoc5KAU98xymAtYJRBg8FgMFTG5thUrn1vHSM7hbA5NpXWzXz5+p4LzjrTVinFTVaNx14RQWyKTeGdGweUcmkX2uwMfm4JIzqF8vq0fqfH03IKGPH8Mnq2CuLTO4ZU6vqPS8slOSufXhHFM9TTcgp48/dDzN94jP9e1dvlNoS7T6Qz/cMNuLsJn90xhM7hpRORXvntAK8tPUj7ED+OpeTw3JU9mTqobaXnziu0MeqF5RTa7CRnF/DytX2Y0r90fGV5KKX4eWcCkSG+dG8ZWKPW67qkPGXQqMnVQESuEJFZ6enplU82GAwGQ6NmQLtm/OVPnVi2PxE3N+G96QNqpOSKiPD81b3xcBc2xKTwv6t6lxnb6OnuxvheLfltz0lyCs4koLy7Ipqs/CKemNjdpRjQiKZN6N26aSnFqKmvF49f3p0dT46rUj/qHq2CWDBjKAJMfW8dR60EHgfZ+UXMWRfDmG7hfHfvcIZ1bM7/fbWT//6yD7u9YkPWvA1HScjI47Xr+tG1RQBvLTtULJO8MnafyGDm51uY8PpqLnl1JW8vP8SJtFyXj29oGGWwGiilflBK3RUUVLUMK4PBYDA0TmaOjmLm6I58cNPAGu3X3DKoCR/fOphZ0wdwzcA25c6b2KcVuYU2lu49BcCpjDw+XnuESX1a0bVFYLnHVYXqWM+iwgJYMGMYRTbFwwu3FVPY5m08RlpOIXeP6kigjycf3TKI64e05d0Vh/nX97soz7OZV2jjreWHGdI+mOFRzbn34igOJ2azaFeCy+taH50MwGPjuxLo48nzi/Zz0QvL+GH7iSrfY0PAKIMGg8FgMNQy7m7CXy/pysDI4Bo/94B2zRhXhkXQmUGRwYQHevO9pcy88fshimyKB8uJEzyXRIb48eTEHmyMSeWDVdGATv74YFU0g9sHn0688XB347nJPZlxUQc+XX+U/y7aV6ZC+On6WBIz83lobGdEhPE9W9Ih1I83lx0qNv9URh4HTmaWuaZ1h5PpEOLHny/qyJd3X8DKv46mb5um3D9vK5/9EVsLUqhbjDJoMBgMBsN5jrubMKFXK1bsT2RXXDpfbDjK1EFtipWLqUum9I/g0h4teGnxAfYlZPDdtjji0/O4e1THYvNEhMcu7cqNQ9vy3opo3lp2qNj+nIIi3ll+mAujQk7XVXR3E2aOimJvfAa/7ztFXqGNN38/yEUvLGfSm2vIdqrdCFBks7PhSEqxuoxtm/vyyW1DGNU5lH98s4u3SiiWDR3TadtgMBgMhkbAxL6tmL3mCLfP2Yi7m3C/U+3DukZEeO7Knlzy6koenL+dgiIb3VoGMqpz6RaGIsLTE3uSk2/jxcUHyC+y0ybYl7ScArbEppGcXVDK4jmxbyteXXqAf/+8lwKbnWMpufRv25QtR9NYcSCRy3qdiXXcE59BZn4RwzoWL9LdxMudWTcN5JGF23nh1/0s23eKApud5KwCUnMKGBgZzANjOtG/beni4fUdowxWA6cOJHW9FIPBYDAYXKJP6yDaBDfhWEouM0Z2IDzQp66XVIzm/t78d0pv7vhEV+l4fVq/cuMQ3dx08kxOgY03fj9jHfR0F6YNblOqm4unuxt3XxTF37/ZSedwfz67YwhD2gcz+N9LWbQroZgyuO6wjhcc2r60S9/T3Y1Xru1L62ZNWHUwiWA/LzqFBeDr5c5PO+OZ8vZaLuocyoNjO9O3TdOzlsm5wpSWOQtMaRmDwWAwNCTeWHqQj9bGsPShi1zqV1wXPPPjHjbHpvLln4dVWhvQblccOJWJn5cHTX098ff2KFeBVEqxOTaVvm2anj7v/325g593xrPpn2cKct/60QaOpuSw9OFRVVp3dn4Rn6yLZdbKw6TmFPLRLYMY3TWs8gPPIabOYC1glEGDwWAwNCTsdkVekQ1fr/rtGFRKnZPafsv2neLWjzfy0a2DGN0ljCKbnb5P/8bkfq14dnKvap0zK7+IiW+uBgWLHhjpcpvFHcfTmLUymtWHkujTuikjOoUwolMoncP9a0wWps6gwWAwGAyNHDc3qfeKIFSvTE11uCCqOf7eHvxqlZ3ZGZdOVn4RQzs0r+TI8vH39uCfE7oTnZTN3PVlZx4X2eycyshj94l0ft4Zz3Wz1jHxzTWs2J/IyE6hHE/N4dmf9nLJqyuZ/uGGWk9Wqf+vCIPBYDAYDIZawNvDnYu7hrF4z0meu1KxPjoF4KyUQYBRXUIZ2TmU15YcYEq/iNMu+bScAmZ+voW1h8/0dgZoEejD3y/ryrTBbQnw8QTgRFouc9bG8N7KaNZFJxdrMVhV0nMKeeL7XeXuN8qgwWAwGAyGRsulPVvw/fYTbIxJYV10Mp3C/Anx9z6rc4oIj0/oxvjXVvHqkgM8NaknJ9JyuWn2Bo4m53DXyA60btqEEH9vwgK96RXRtJQ7uVXTJjw4tjNfbj7OB6uOVFsZzM4v4paPN7D1aFq5c4wyaDAYDAaDodFyUedQvD3c+HHHCTbFpHD1ANd7GFdE5/AArh/clk//OMqQDs155sc9ZOUVMee2waXK1pSHj6c704e149UlBzl0KouoMP9Sc7LyiziVkcepzHwEGBgZfLq9YH6RjbvmbmLH8XSu7BfBq+VcxyiDBoPBYDAYGi1+3h6M7BzKgo3HKbDZz9pF7MyDYzvz7bY47vlsC6EB3syfMYzurarW/u/Goe14e/lhZq85wr+vPJPUsuVoKjPmbiYxM7/Y/IimTZg2uA1XDWjNk9/vZs2hZF66po+utVjONYwyWA1MnUGDwWAwGM4fLu3Rgt/2nATOPl7QmWA/L568ogfzNh7l5Wv7VqsvdYi/N1f1j+Crzcd5ZFwXgv28OJqcw51zNuHn7cFj47sSFuBNWIAPKTkFzNtwlBcXH+DFxQcAeOKK7lxVibXTlJY5C0xpGYPBYDAYGj5pOQUMfHYJUWH+LHpgZF0vpxQHT2Yy9pWVPDS2MzcNa8eUd9aSkl3A13dfQIfQ0q7j6MQsFmw6TkRTH6YPizw9Xl5pGWMZNBgMBoPB0Khp6uvFvRdH0aZZ1S1354JO4QGM6hLKJ+tiWHMoieMpuXx6x5AyFUGADqH+PDa+q8vnN8qgwWAwGAyGRs8DYzpXPqkOuePCDtz44R8kZaXw2nV9GVxGu7zqYpRBg8FgMBgMhnrO8KjmXNW/Nb0iApnUN6JGz93olUERuRa4H+gLJCmlIut2RQaDwWAwGAzFERFeurZPrZzbtKODVOBN4B91vRCDwWAwGAyGc02jtwwqpX4DEJHJdb0Wg8FgMBgMhnNNnVsGReRvIrJQRKJFRIlITAVz3UTkQRHZJyJ5InJMRF4SEb9zuGSDwWAwGAyG84Y6VwaBfwMXA4fRLtuKeAV4GdgD3AcsRMf7/SAixe5FROZZymV5j1E1ficGg8FgMBgMDYz64CbuqJSKBhCRXUCZRXNEpAdaAfxaKXWV0/gR4HXgOuBzp0PuBO6t4LrpZ7lug8FgMBgMhgZPnVsGHYqgC0wDBEq11nsfyAFuLHHeTKVUUgWPwrNfvcFgMBgMBkPDpj5YBl1lEGAHNjgPKqXyRGSbtb/KiIg74Gk9RER89GlVfsVHGgwGg8FgMDR86twyWAVaoesAlqWkxQEhIuJVjfNOB3KBBUBb6+/91V6lwWAwGAwGQwOiIVkGfYHyrHV5TnMKqnJSpdTHwMeuzheRu4C7rKf5VpyjoWJCgKS6XkQDwcjKNYycXMPIyXWMrFzDyMl16qOs2pU12JCUwRwgrJx9Pk5zahWl1CxgFoCIbFJKDaztazZ0jJxcx8jKNYycXMPIyXWMrFzDyMl1GpKsGpKb+ATaFexdxr4ItAu5SlZBg8FgMBgMhsZOQ1IGN6LXO9h50Er46AtsqotFGQwGg8FgMDRkGpIyOB9QwAMlxu9Exwp+ds5XZLmLDZVi5OQ6RlauYeTkGkZOrmNk5RpGTq7TYGQlSqm6XYDIdM4ENN4HeAEvWc9jlVJznea+gS4k/Q3wM9AN3YFkDXCxUsp+rtZtMBgMBoPBcD5QH5TB5cBF5exeoZQa5TTXHW0ZvAuIRGfpzAf+pZTKqtWFGgwGg8FgMJyH1LkyaDAYDAaDwWCoOxpSzGCdIyJuIvKgiOwTkTwROSYiL4mIX12vrS4Qkc4i8rSIrBeRRBHJFJFtIvKPsmQiIl1E5FsRSRWRbBFZJSIX18Xa6xoR8RWRaBFRIvJmGfsbtaxEJFhEXhSRQ9Z7LVFElonIiBLzhojIEuu1lyEii0Skb12t+1wiIv4i8ncR2Wndf5KIrBWRW0RESsxtFHISkb+JyEKn91ZMJfNdlouItBKRT6zXYq6IbBKRa2rlRmoZV+UkIj4icqeIfCciMdZ9R4vIFyLSrZxjvK3vhSMiki8ih0XkcRHxrNWbqgWq+noqcez/rGPK9FrWNzkZy2AVEJHX0DGK3wC/oGMW7wNWAWMaW8yiiPwXmAl8D6wHCoHRwLXADmCoUirXmtsR3UqwCN1fOh2d/NMTGK+UWnLOb6AOEZEXgRmAP/CWUupep32NWlYi0g5YjpbNh8ABIAjoDfyqlJpnzRtqzYsDHAr1veh6pBcopXae04WfQ0TEDVgBXADMQb//fNE93AcDzyul/s+a22jkJCIKSAG2AAOADKVUZDlzXZaLiASjK1aEAS8Dx4Hr0SFOtymlPqqF26k1XJWTiHQF9gKrgcXoEm8dgLsBP+BSpdSyEsd8C0wCZgPrgGHAbcAcpdQttXNHtUNVXk8ljuuLroCSh9az/MuYU7/kpJQyDxceQA90b+SvSozfh85yvr6u11gHMhkIBJUx/qwlk3udxhYANqCv05g/EItu/yd1fT/nUG790YreQ5ac3iyxv1HLCv3j6hjQspJ5G4AMIMJpLMIaW1zX91HLMhpmvXZeKTHuBUQDaY1RTkAHp793ATE18foBnrfkfYXTmLt1jmTAv67vvTbkBDR3/hxyGu+O7gi2qcT4ZZacXiox/pI1fkFd33ttvZ5KvC42oo0ky4GsMubUOzkZN7HrTAMEbalx5n1055Mbz/mK6hil1CalVHoZu+Zb254Alst4IrBcKbXN6fgs4AOgMzColpdbLxCdBPU+sAj4uoz9jVpWIjISuBBt2YoXEU8R8S1jXhRaDguVUnGOcevvhcAYEWlxrtZdBwRa2xPOg0oX3k8CsqHxyUkpFe3KvGrI5XrgsFLqB6e5NuANIBj95d5gcFVOSqlk588hp/E9aOWoZ4ld11vbkt+TjucN6nvSVTmV4H60snxfBXPqnZyMMug6g9CWwQ3Og0qpPGAb5/EXdDVobW1PWtvegDfaFF6S9da2scjvQaAr2h1VFo1dVo4v1aMi8gOQC2SLyAERcf6AdMigPDkJ2q1zvrIBSAMeFZFrRKStiHQVkf+g7/tJa15jl1N5uCwXEWmJthiuL2eu8/kaBVaYQkvOfMY7GATEKaWOOQ9az09wnsvJCnF5BnhKKRVbwdR6JyejDLpOK3TLu/wy9sWhW+V5neM11Tssy9c/0W7Qz63hVtY2roxDHGMRtby0OkdE2gNPAU8rpWLKmdbYZdXF2r6PtrjcjI6jKQDmisit1v5GLSelVCragpyCDiuIRcd2zQSuUkq9b01t1HKqgKrIxciwNH9GK4NzSoy3omw5YY2f73J6Bx2m8XIl8+qdnDzO9QUbML7oGImyyHOa09j7I7+Kjmf6u1JqvzXmcPOVJb+8EnPOZ96l8g+Kxi6rAGubCYy23J6OYOto4N8iMgcjJ4AstKvue2AtWnmeCXwuIpOUUr9h5FQeVZGLkaETInIB+jNsO/DvErsr+548b+UkItOAS4ELlVJFlUyvd3IyyqDr5KAzycrCx2lOo0VEnkG7P2cppf7jtMshF+8yDmsUsrNcnGOBkUqpwgqmNnZZ5VrbLxyKIGhLmIh8D9yEth42ajmJSC+0AvigUupdp/Ev0Ari+1ZWeqOWUwVURS5GhhYiMgD4Ce3KnGCFSTmTQ9lyAi2r81JOVrb5q8CHSqm1LhxS7+RklEHXOQF0FxHvMlzFEWgXcqO1CorIk8DjwEdoF4IzjiD3skzfjrHyTOYNHhHxRv+S/hlIsILX4cy9B1ljSTRyWaFLdgAklLEv3to2w8jpQfSXxkLnQaVUjoj8hP5RFomRU3lURS5GhoCI9Ad+Q5e6Gu2ceOPECcp3cUZw/srpCXSpnfedPt8BmgBijeU7xQjWOzmZmEHX2YiW12DnQRHxAfqia1A1SixF8Al0/MgdysqRd2In2iQ+rIzDh1rb81l+TYBQYAJw0Omx3Np/o/X8DoysHAlarcvY5xg7hX4/QvlyUsDmml1avcLxReJexj4Pp21jl1N5uCwXpVQ8+st5aDlz4fx+TzoUwSWcCd8oLzliIxAhIm1KHN8GHSd3vsqpHVoZ/IPin/GD0S7fg+jaxA7qn5zqsoZPQ3oAvai4zuCNdb3GOpLLv6z7/wRwq2DeQnTtvD5OY47aeQc4j2vnAZ7A1WU87rZk94v1vLORFc3Qdd6O41S7DR2sngXsdxrbaM1t5TTWyhpbUtf3UstyesV67TxaYrwp2uqQArg3ZjlReZ1Bl+UCvED5dQZTgYC6vt9alFM/dC3FozjV3Stn7gQqrp93YV3fb23ICf2joqzP+N3o0JergbH1WU6mA0kVEJE30O6Xb9Auv27omkJrgItV4+tAMhNduf8oOoO45P2fVDqI3VHXawO6S8kr6A/cO9FK9gSl1K/nat31BRGJBI5QugNJo5aViNwFvIf+IJ2NLqR8N1ohvFwptdiadwGwDK04vmEdfh8QDgxXSm0/x0s/Z1glLLaglefP0J9BwejXSSQwUyn1tjW30chJRKajrTSg79EL/QULEKuUmus012W5iEhztKWwOTrkIw5de3YU2hvyYS3dUq3gqpys19lm9GvrKeBwGaf7RimV7XTuH4DL0d2DHJ01bgc+VUpNr/m7qT2q8noq5/jlwEBVdgeS+iWnuta2G9ID/UvwYXQXiHz0B8LLNLDq8zUoj4/Rv2LKeywvMb8b8B26PloOusXRmLq+jzqUXyRldCAxslIAU9A13LLRrqnF6C/okvOGAUvRVsNM4Fegf12v/xzJqCM6NOM4+odDBrASmNJY5YQOvXDp86iqckG75ueiY3vz0Mr41Lq+59qUE1rZregzXgGRJc7tg+5CFWN9T0ajjQWedX3ftf16Kuf4Uh1I6qOcjGXQYDAYDAaDoRFjEkgMBoPBYDAYGjFGGTQYDAaDwWBoxBhl0GAwGAwGg6ERY5RBg8FgMBgMhkaMUQYNBoPBYDAYGjFGGTQYDAaDwWBoxBhl0GAwGAwGg6ERY5RBg6EBIyK3iIgSkVF1vZa6RkTGi0iRiHR1GnvSkk+ki+eIsboG1BtE5GMRcbkgrHW/H9fikhCRSOs6T9bmdRobIjJJRApEpFNdr8XQuDDKoMFQDxCRUdaXq+NhE5FUEdklInNE5FIRkRq+5pMiMrkmz1lXiIgHuk3UZ0qpfXW9HkPdU8Z7yvnxYznHXCYia0UkW0RSRGShiLQvZ24XEfnWep9mi8gqEbm4nLlBIvKGiMSJSJ6I7BaRu0u+p5VS3wE7gf+d7f0bDFXBo64XYDAYivEFuu+1AAFAF2AycBOwRESuUUqlOc2fC8wDCqpxrSfQ7cy+PasV1w+uQbfwm3aW5+mCbjXVkGkC2Gr5GrHWdYpq+To1wSxgVYmx4yUnicgU4EtgO/BXIAh4AFgjIgOVUiec5nYE1qLv/3kgHd0X+lcRGa+UWuI01wv4DeiH7oG8FxgPvI3ug/xkiaW8BswRkR5Kqd3VvGeDoUoYZdBgqF9sUUp96jwgIg+hv3AeQiuL4x37lFI2av+LvyFwD7BDKbX9bE6ilMqvofXUGUqpvHNwDYXuz9sQWFfyPVUSEfFEK2rHgBFKqSxr/BdgM1phu8vpkP8ATYEBSqlt1txPgN3AWyLSVZ3p9XoHMAi4Xyn1hjX2voh8BfxdRD5SSsU6nftr4B3gz8B91bxng6FKGDexwVDPUUrZlFIPA6uBS0XkQse+smIGRcTHcgHvF5EcEUkTkZ0i8oK1P9IpBu1mZ/eZ0zmmisj3InJURPJFJMlyifUuuT5HnJ2IdBWRn0QkU0TSReRLEWlRxvxAEXlORPZaLrNkEVktIteVmNdSRN6x1lAgIidEZJaIhJWY1wK4EG1RLQ8/EXldRBJEJFdE/hCRP5V3L9W9PxEJFpFXROSw071tFpG/VrA2lxCRUBH5xDpntogsFZH+ZcwrFTPoGBORYSKywjo+WUQ+EBH/EnPbiMhsEYm1/venLNfpzU5zSsUMWjIqzy0bU+IanURkrojEW//bGBF5QUT8zlZOZSEifiLiU8GUi4BWwAcORRDAUvSWA1MthRFrjROB5Q5F0JqbBXwAdEYrfw6uB3KA90tc81XAE5jqPGidZxVwdRVu0WA4K4xl0GBoOHyIVnomoBXD8ngLuA34BHgZ/T7vBDjimRKB6WgX8yq0G60k9wLJ1r4EoCPaMrJGRPorpQ6WmB+B/tL8Bu1i6wPMAAKBcY5JItLUWnsPtEvuHcAd7UK7HO3yRkTaAusAL+u+DwNRwN3AaNFuu3TrtBdZ2w0VyOQTtAX1f2j3+wxgkZRw6VWAS/cHLARGAu8CO9Cu1G7AKOAFF65TEYuAFLSVqgX6f7RCRIYppXa5cHxf4EfgI+Bza023A3Ysq5fo2Mvf0Pf7NnAA7S7tDYxAhxWUx3NoZciZjtZ6TzoGRGQA8DuQBrwHxKHleT8wXEQuUkoVWnO90f8vV7AppVLLGH8Nfc+IyEH0++N1J8sdnFHe1pVx/Hr0e6cz2vLXG/CuYK7jfBtExA3oj7b4l7SkbkCHJAyiNOuAS0RbGE0MrKHWMcqgwdBw2GFtO1cy70rgF6XUzWXtVEplA5+KyFwguhwX2qXWvNOIdoNtAx5Eu2WdiQKmKqUWOM23A/eISBel1H5r+N9oRXCGUqqYEmp9cTp4A2016aeUOu40ZyH6C/dBzsRadbe2h8u6X4sitPuvwDrPbGCfdZ1uFRzn8v2JSBBaaXhHKVUb7r1Y4CqHEiMiXwMbgReBS104vjcwTCn1h/X8PREJBG4VkYcsi1R3dNzk/ymlnq/K4pRSvzk/F5FmaKUmGW0dczAbiAcGKaUyneYvRbtIbwA+toanYSlyLhALRDo9LwS+R1uMT6Atf7ejLXJ9gVud5raytnFlnNcxFoFWBl2dC9AM/YOg1FylVL6IJDnNdcbxWu6Bfp0aDLWKcRMbDA2HDGsbWMm8dKCHiPSs7oUciqBoAkUkBG1R3A8MKeOQE86KksXv1raTdS434Dpgb0lF0Lqm3ZoXhLYSfg/kiUiI4wHEAIcobo0LtbYpFdzSKw5F0LrWceAzoKuIuKIMVnp/QC6QDwwRF0vZVJHnna1ZSqnNaCvemJKu3nJY56QIOvgdbRSItJ47rK2jS7rjq4LlUv0KaA9cqZQ6bI33QiulnwPeJf63q4Fsiv9vfwXGuvi4wXkNSqk1SqlJSqn3lFI/KKXeA4Za57xFRIY7Tfe1tmXFjOaVmFNTcx3zfcsYT7a21f4fGAxVwVgGDYaGg0MJzKhwls6AnAvsFJFoYBnwA/CDQ+GqDBHpBzyDdiWWjOM6UsYh0WWMOb7QmlvbELSlZFEll++C/qF6u/UoC+frORSkikrv7C1jbI+17VDO/vKu56DY/SmlCkTkAbRb8oiI7EErW98qpZZWcn5XKO8exgHt0FarinDlHnImHY8AAAZ0SURBVGJF5Dngb0C8iGwDlgILlVIbq7DWWcBo4CallHMmr0Pxfsp6lEW44w+lVDzailgjKKXsIvIf4BJ0uMUaa1eOtfUu4zCfEnNqaq5jfk4Z447XckPPbDc0EIwyaDA0HBzJG/srmqSU+s6yTF2Gjqcbg1aqVonIGGcLWVlY8Xor0UrnM9b1stFfTK8CZVmhKspormp9RMf8Tyk/Ri3X6e9EaxuMzgatDVy6P6XUuyLyHVrRuAidBHCviMxXSl1X3gnOEa7ew+OWG30COk7wDuCvIvK8Uur/KruIiPwduAV4Vik1t5zrvET5PwpOx/2JSBN0zKIr2JRSiZVPI8bahjiNOcrGRFBa6Xa4cePKmFuSknNT0a/VUnOteMgQYEUZ5wm2tq7cj8Fw1hhl0GBoODisZD9VNlEplYJWpj4VEQH+CzwKTEInOVTElWiFb6JSapnzDhFpTvkur8pIQn859qlk3iG04unlYnKHI3miE7pGXFl0K2OfI9awLItZtbGsWR8AH4iIO9pKO01EXqqida0k3TiToOCgO1rJiy09vfoopaLR8ZRvWFm4vwKPWvdwqrzjRORa4FlgPvCvMqY4Eo9sLv5vp1L9mMHycLj1TzqNOf4vw4CS6xqK/mF0wHq+E/0eGFbGuYda201w2hK5BegnIt4lShcNRivHm8o4T5S1dSUxyGA4a0zMoMFQzxERdxF5Eat8ilJqTSVzmzqPWXFmW62nwU67sko8d+CwIBWz6InInegs1mphuai/ALqLSCn3r6W0opRKRgf9TxGRoWXNE5FQpyGHZaXUXCceFF3813GO1uikhv1KqcpcxC4hIr4iUiz+y6oD6Uj8KUvWVeFRh4ys6/VHW32XOpdDORtEd8rwdB6zsmAdMmpWwbFD0ZbcP4BbSmTrOtiKVnD+LCIdyjiHh4g4y6naMYPWD5eS5/fmTOLRD067VqDd0Xc4x1+KSB90qMRCR4azJesfgFHWfsdcf7QV9SDFM9u/+P/27ua1iiuM4/j3AelOzKJduJBYW5VWVw0FQSoGaxHbRZUiCCnJpmCgKIIQCkLFP0ApuqpSQl0ofcEqSEsWpdQsVCh9ARFLCpcuiiYEQ+0LUsvj4ncik0muCSQxrfP7wIHk5szMncvk3ueec55n0LrAap1C0HKO+yhwrtsE3K4kXpktKo8Mmv23vBQRPeXn6h1IOoEhpmZlzmQ5Wut1EX3wjqJF/P1oVK76AXgFJR8MAL+iuPEc8CVax3QmIk6W7TajaedfmN/7xmGUcXs6Il5DSQOBSsssQyVvKM93GPi2ZDF/j768rkGjmx9TPtQzcyxUG3AncKjNcZehafKz6DXah7I898/jXOrWoVIv51HAcweN5vWjdZYP186F6gD2At2Z+c0c99+J7nBxEViJSsv8jUrdLJRu4MNQQeSb6AtDFwpyrs4SnFxAGeCfAm/F1Dut/ZGZX2RmRsTbaC3lT2U6+joKlp4HdqP1ioMw7zWDX0XEb6ho9GQ2cQ8aGTyRmQ8Dtsz8JyIOoMDsckScQmt0D6Kp2vdr+34P2AYMRcRxNHL4DpoOfr0WCJ9CmcvHyvKNG+ha3YWm0lvVHZeg8hWUdW32eGSmm5vbEjc0+pCV9i/K7LyORlt2tNmur/TfWn5/Ct0d4RpKDriH1kh9BKytbbsWBZi/Tx638rctKBi7i+rBXQI2olp7rdp+WqgAb7tz6qs93oHuqDKCbqM3jgKlPbV+T6PafD+jrMsJNEX3AfBire+ecqyu2uNHyuMb0LTnrbKva8D2GZ7ztHOZ6/mhJIzjqPzOBArURtA6y5W1bT9H9f3Wz+HaGCzHeQZNOY+jYP3r+vmW/gkMzvZYm+vnWVQj8Ua5Lv4sPx8FVlS2W122O1I7RrtWv2Y6y3FalWvgO3Ttrlqg/6kBVNpmDJWZmUDJVHsfsc0b6EvSXyiY/wx4rk3fF1AAPFH6DwOvtunbAZxEQek9lPjzLhAz9O0tr9nGhXp/cXObrUWmk5XM7P+trM37EfghM3tm67+USomdUeBStqkFac1V1hi2MnP3Uj8Xaw4Hg2b2RIiIHZQRzFygdYCLISJeRmvU1mWloLZZRLwJfAJsyOl3+TFbNA4GzczMzBrM2cRmZmZmDeZg0MzMzKzBHAyamZmZNZiDQTMzM7MGczBoZmZm1mAOBs3MzMwazMGgmZmZWYM5GDQzMzNrsAdtX2p22epDFAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "plt.rcParams['font.size'] = 18\n", "\n", "dot_hic_path = \"../../../tests/test_data/dothic_chr9_4000000_6000000.hic\"\n", "dhic = DotHiC(dot_hic_path, balance=False) \n", "\n", "gr1 = \"chr9:4000000-5000000\"\n", "gr2 = \"chr9:5000000-6000000\"\n", "\n", "mat1 = dhic.fetch_data(gr1) # any genome region\n", "decay1 = dhic.diagonal_mean(mat1)\n", "mat2 = dhic.fetch_data(gr2)\n", "decay2 = dhic.diagonal_mean(mat2)\n", "\n", "binsize = dhic.fetched_binsize\n", "\n", "# plot\n", "fig, ax = plt.subplots(figsize=(10, 5))\n", "ax.plot(decay1, label=gr1)\n", "ax.plot(decay2, label=gr2)\n", "plt.yscale('log') # log scale on y-axis\n", "plt.xlabel(f\"Distance(bins, binsize={binsize})\")\n", "plt.ylabel(\"Mean contacts\")\n", "plt.xlim(0, 150)\n", "plt.ylim(1e-1, 1e5)\n", "plt.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Another function `.diagonal_mean_std` can get the standard error of contacts in the specific contact distance at the same time, we can show it with a error band:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "tags": [ "nbsphinx-thumbnail" ] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAFVCAYAAACU4eUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeZicZZX///epfek13Z2VJBASJEFZTAcQdcLmgkxcYNxGFEcG/KqgwDg/h0HBQcVxRpYB+c4XUAdhVBCFGSKM2yiIIkJYZAlBluxrdzrpvaqrq87vj+oqutNLunpJL/m8rquupO7nfp463WnIyb2c29wdERERETk4BSY6ABERERGZOEoGRURERA5iSgZFREREDmJKBkVEREQOYkoGRURERA5iSgYBM7vNzLrMrK3X650THZeIiIjIeAtNdACTyC3ufuFEByEiIiJyIGlkUEREROQgNqWSQTO7zMzuNrNXzczNbMMQfQNmdomZrTOzlJltNrNrzCw5yC0fMbMmM3vBzC43M42aioiIyLRnU+kEEjNzoAl4ElgOtLj7oYP0/Tfgs8C9wP8AS4GLgIeB090916vvG4EtQCPwRuCHwJ3u/qVx+2JEREREJoGplgwucvdXe37/HFA2UDJoZkcBzwL3uvvZvdovAm4APuLuPxjic/4a+Cd3XzLGX4KIiIjIpDKlpokLieAwfBgw4Pp92m8FOoBz9nN/rud+ERERkWltSiWDJVhBPqF7rHeju6eAp3uuF5nZB82s0vKOBq4E7j5QwYqIiIhMlOmaDM4FGt09PcC1rUCtmUV6tX0a2AC0Aj8hv2bwivEOUkRERGSiTdcdswlgoEQQINWrTxeAu68c7oPN7ALgAoBkMrn8yCOPHEWYIiIiIgfGE0880ejudfu2T9dksAOYOci1WK8+JXP3W4BbAOrr633NmjUjeYyIiIjIAWVmGwdqn67TxNvITwVHB7g2j/wUctdIH25mq8zslubm5hEHKCIiIjIZTNdk8HHyX9vxvRvNLAYcC4xqOM/dV7v7BZWVlaN5jIiIiMiEm67J4F2AAxfv034++bWC3x/NwzUyKCIiItPFlFozaGYfBRb2vK0DImb2xZ73G939DgB3f9bMbgIuNLN7gAfIn0DyWeAhYNCC08Ph7quB1fX19eeP5jkiIiIiE21KJYPAecC+O3+/0vPrQ8AdvdovJl8u5gLgTPJHzd0IXNH7KDoRERGRg9mUSgbd/eQS+maBa3peY8rMVgGrFi9ePNaPFhERETmgpuuawXGlDSQiIiIyXUypkUERkekmnU7T1NREa2sr2Wx2osMRkSkkGAxSXl7OjBkziEYHqqY3PEoGR0DTxCIyFtLpNJs2baK6uppDDz2UcDiMmU10WCIyBbg7mUyGlpYWNm3axIIFC0acEGqaeAQ0TSwiY6GpqYnq6mpqa2uJRCJKBEVk2MyMSCRCbW0t1dXVNDU1jfhZSgZFRCZIa2srFRUVEx2GiExxFRUVtLa2jvh+JYMiIhMkm80SDocnOgwRmeLC4fCo1hwrGRwBnUAiImNFU8MiMlqj/f+IksER0JpBERERmS6UDIqIyKTx5S9/GTNjw4YNEx2KyEFDyaCIiExbTz75JO95z3uoqakhFotx1FFHcf31149qfdUXvvAFzIyysrIBr6fTaa644goOO+wwotEohx9+OF/96lfJZDID9r/99ts57rjjiMfjzJo1i7/927+loaFhwL5//OMfOf300ykvL6eiooJ3vvOdPP300wP23bZtGx/72Meoq6sjHo9TX1/P3XfffcBjHsqWLVu44IILimVRZs+ezRlnnMHatWv79X3ggQc46aSTSCaTzJgxg/e///2sX79+wOe++OKLvPe976W6uppkMslb3/pWfv3rXw/Yt7m5mYsuuoh58+YVf0b+/d//HXfv1zeXy3Hddddx5JFHEovFmD9/Pn/3d39He3v7gM8er5jHnLvrVeILWAXcsnjxYhcRGam1a9dOdAiTzpVXXumAr1+/ftTPeuihhzwajfrMmTP9y1/+st90003+3ve+1wE///zzR/TMp556ykOhkJeVlXkymRywz3ve8x4H/BOf+ITfeuut/olPfMIBP/fcc/v1vfbaax3wlStX+s033+xf+tKXPJlM+rJly7ytra1P3z/84Q8ejUZ90aJFfu211/q1117rixYt8rKyMn/mmWf69N29e7cfdthhnkwm/Utf+pLffPPNvnLlSgf8u9/97gGLeShPPvmkz5gxwxcuXOhXXnmlf+c73/F//dd/9Q9+8IP+4IMP9un7k5/8xM3Mjz32WL/pppv86quv9pkzZ/qcOXN869atffq+/PLLPmPGDJ85c6ZfffXVftNNN/mxxx7roVDIf/nLX/bpm06nfcWKFR4KhfySSy7xW265xd/3vvc54FdeeWW/mD/72c864O973/v8lltu8UsuucRDoZCfcsopns1mD0jMgxnO/0+ANT5QXjNQo17Dey1fvny/33gRkcEoGexvLJLBlpYWd3c/5phjPB6P+yuvvNLn+gUXXOCAP/zwwyU9t7u72+vr633VqlW+cuXKAZPB+++/3wG/9NJL+7RfeumlDvjvf//7YltDQ4MnEglfsWKFd3d3F9vvu+8+B/xrX/tan2esWLHCy8vLfcuWLcW2LVu2eHl5ub/tbW/r0/fv//7vHfD77ruvT/wrVqzwGTNmeGtr6wGJeTCdnZ2+ePFiP/bYY725uXnIvl1dXT537lxfsGBBn7ifeuopDwQC/RL797///R4IBPypp54qtrW2tvqCBQv8iCOO8FwuV2y/6aabHPAbbrihzzPOOussD4fDvmHDhmLbc88952bmZ511Vp++N9xwgwP+/e9//4DEPBglg0oGRWQKOtiSwXQ67d/4xjeKSVpFRYUvX77cb7zxxmKfQjK4bt06v+yyy3zevHkeiUT86KOP9vvvv7/P89avX18cwbnzzjv9jW98o8diMT/33HO9qanJAX/HO97RL45HHnmkOArW27Zt2/yFF17w9vb2AeO/9tprPZFI+IYNGwZNBj/ykY844Js2berTvmnTJgf8U5/6VLHt1ltvdcBvv/32fs9ZtGiRL126tPj+pZdeGjBmd/dPfOITbma+ffv2Ytu8efP88MMP79f39ttvd8DvuuuucY95KN/73vf6JKupVMpTqdSAfX/5y1864FdddVW/a6eeeqpXVFR4V1eXu7u3tbV5NBr1U089tV/fq666ygH/4x//WGx785vf7IlEwjs7O/v0/e1vf+uAf+Mb3yi2XX755Q74b3/72z59Ozs7PZFI+BlnnHFAYh7MaJJBrRkUEZFx19XVxTve8Q6+8IUvMGvWLK666iq+9rWvsXz5cu65555+/c8991wefvhhPv/5z/OVr3yFhoYG3vve9w64seS//uu/+NSnPsU73/lObrjhBs444wzS6TQAiUSiX/9C26OPPtqn/bLLLmPp0qU89thj/e7ZuHEjX/rSl7jyyitZuHDhoF/n448/zrx585g/f36f9vnz5zN37lwef/zxPn0B3vSmN/V7zoknnsi6detoa2sbVl9354knngBg+/btbN26lRNPPHHAvr2fN54xD+WBBx4AoKqqir/4i78gHo8Ti8U47rjj+PnPf96n7/4+s6WlhT//+c8APPPMM6TT6UH79n5eLpfjySef5LjjjiMWi/Xpe/zxx2Nm/b72QCDA8ccf36dvLBbj2GOPLen7NNKYx4vOJhYRkXF3/fXX8+CDD3LZZZdx9dVX97mWy+X69a+trWX16tXF+mmnnHIKxx9/PDfffDNf//rX+/R9/vnneeaZZ1i6dGmxzd2pra3l0UcfpbOzk3g8Xrz2m9/8BoDNmzcPO/5PfepTLFq0iEsvvXTIftu2bWPZsmUDXps3bx5btmzp07fQPlBfd2fbtm0cccQR++0LsHXr1mE9t3ff8Yx5KC+++CIAZ599NieccAJ33nknTU1NfO1rX+Nd73oXP//5zzn99NNL+nqOOuqokr72PXv20NnZOWDfaDRKbW1tv+9TbW3tgOf/zps3j0ceeYSuri4ikci4xTxelAyKiEwy/7T6edZua5noMPpYNreCK1cdNeL7v//971NdXc0VV1zR71og0H+S6nOf+1yfQrorVqygrKyMl156qV/fM888s08iCPkivJdccgmXX345Z511FldddRW1tbX86le/4sorryQUCtHR0dHnnttuu43bbrut3/N/+MMf8rOf/Yzf/e53hEJD/7XZ0dExYLIA+RGk3p9Z+P1A/QsjVYU+49V3PGMeSuHotCOPPJL77ruv+Gd92mmnsWzZMi6//PJiMjgR36dC/1K+T4U+kUhkXP+8xoOSwREws1XAqsWLF090KCIiU8JLL73Escce2286bjCLFi3q11ZTU8Pu3bv7tQ82CvUP//APdHR0cM011xSn9srKyrj22mu5/PLL6e7u3m8cTU1NXHzxxZx33nmcdNJJ++2fSCSKU9T7SqVSfaatC79Pp9N9Ri4LfXv36d13oOeOtO94xtzU1ERXV1efPrNnzwYo3vuxj32sT9K/ZMkSTjrpJB5++GHa29tJJpPj9rUP1Xewr33Xrl2D9h2POAb68xoPSgZHwN1XA6vr6+vPn+hYRGT6Gc0I3HQRDAYHbM+vge9rsL8oA4EAX/3qV7nssst49tlncXeOOeYYcrkcn/zkJwdcU7evf/qnf6K9vZ3zzz+fl19+udje2dmJu/Pyyy8TjUaL6+3mzp076JTe1q1b+0wFzp07t9i+7+DC1q1bMbNin959B3ouvDalWErf8Yz5rLPO4qGHHurTp/Dnd8ghh/Dcc88Vk8Pe5syZg7vT3NxMMpns85n7jgCP5muvrq4mHo8P2DedTtPY2MjKlSv7fO1r164lnU73G8XbunUrtbW1RCKRfnGMZczjRRtIRERk3B1xxBGsW7du0FGY8ZRMJjnxxBN505veRCKR4Gc/+xnuzrve9a793rtx40ba29s54YQTWLJkSfH12GOP0dHRwZIlSzjjjDOK/VesWMHWrVv7rUfcvHkz27Zto76+vk9fgD/84Q/9PvfRRx/lda97XbGw9f76mhnLly8H8snUvHnz+m2QKfQF+sUxHjFfc801/PKXv+zzKiiM1PZej1iwZcsWQqEQM2bMGNZnVlRUFEeH3/CGNxCNRgft2/trDwQCvPGNb+Spp57q93P52GOP4e79vvZcLtdvg1EqleLpp58u6fs00pjHzUBbjPVSaRkRGX8HU2mZf/mXf3HAv/jFL/a71ruG2lB1BhcuXOgrV64svu9dWma4GhsbfcmSJV5bW+sNDQ19rg1UWuaRRx7xu+++u99r2bJlHovF/O677/Zf/OIXxf4//elPh6zZ17u24a5duzwej/vxxx8/YM2+r3zlK32eUV9f7+Xl5X0KFm/dutXLy8v9tNNO69P385///KB1Bquqqoq1GMc75sGsW7fOg8GgL1++3DOZTLH96aef9kAg0KduYldXl8+ZM6dfzb5C3/POO6/Ps//qr/7KA4GAP/3008W2Qs2+JUuW9Pl5+9a3vjVoncFQKNTn5/CZZ54Zss7gHXfccUBiHozqDCoZFJEp6GBKBtPptL/1rW8t1v675ppr/KabbvJPf/rTfRKZsUwG77//fj/55JP9a1/7mn/nO9/xK664wufOnevxeHzAUx3OPfdcB/w3v/nNfr+eweoMurv/5V/+pQN+3nnn+be//W0/77zzHPBzzjmnX99vfvObDvjJJ5/sN998s19xxRWeTCb9yCOP7JNEuLv//ve/90gk4osWLfLrrrvOr7vuOl+0aJEnk8k+SYR7PulduHChl5WV+RVXXOE333yzn3zyyQ74t7/97QMW81Auu+wyB/ykk07yG264wb/85S/7jBkzvKyszP/0pz/16fujH/2oz2keX//6133mzJk+a9asPkW43fM1Gaurq33mzJn+9a9/vXiaRzAY9J/97Gd9+qbTaV++fLmHQiG/9NJL/dZbby2eQDLQP1wuvPDC4gkkt956q1966aUeCoV85cqV/U4gGa+YB6NkUMmgiExBB1My6J4vzvvVr37Vly1b5tFo1CsrK72+vt5vuummYp+xTAaff/55f+c73+mzZ8/2cDjsc+bM8XPOOcfXrVs3YP+xSgY7Ozv98ssv94ULF3okEvHDDjvMr7rqqmKR4X39x3/8hx999NEejUa9rq7O/+Zv/sZ37tw5YN9HHnnETz31VE8mk15WVuZvf/vb/Yknnhiw75YtW/ycc87xmpoaj0ajftxxx/mdd955wGMeys0331x8TlVVlb/vfe/z5557bsC+q1ev9hNOOMHj8bhXVVX52Wef7S+//PKAfdeuXevvfve7vbKy0uPxuL/5zW8e9Fi3PXv2+Gc+8xmfM2eORyIRX7p0qd94440DjsZ1d3f7N7/5TT/iiCM8Eon43Llz/ZJLLhk0CR6vmAe7f38GSwYtf01Gor6+3tesWTPRYYjIFPXCCy/0W1wuIjISw/n/iZk94e79FiBqA4mIiIjIQUzJ4AiY2Sozu6W5uXmiQxEREREZFSWDI+Duq939gsrKyokORURERGRUlAyKiIiIHMSUDIqIiIgcxJQMioiIiBzElAz2MLO4mb1sZm0THYuIiIjIgaJk8DVXARsnOggRObio1quIjNZo/z+iZBAws+XAO4FvTHQsInLwiEQidHZ2TnQYIjLFdXZ2Eo1GR3z/lEkGzewyM7vbzF41MzezDUP0DZjZJWa2zsxSZrbZzK4xs+QAfUPArcBngK7x+wpERPqqra1ly5YtNDU1kclkNEooIsPm7mQyGZqamtiyZQs1NTUjflZoDOMab1cDTcCTQNV++l4HfBa4F7gGWNrz/jgzO93dc736/j3wlLv/1sxOHvOoRUQGUVlZSTQapaGhgd27d9Pd3T3RIYnIFBIKhYjFYixYsIBYLDby54xhTOPtcHd/FcDMngPKBupkZkcBFwH3uPvZvdrXAzcAHwJ+0NO2GPg/wHHjG7qIyMBisRjz58+f6DBE5CA2ZaaJC4ngMHwYMOD6fdpvBTqAc3q1vQWYBfzZzBqB/waSZtZoZn8xypBFREREJr2pNDI4XCuAHPBY70Z3T5nZ0z3XC34E/KrX+zcBtwHHAg3jG6aIiIjIxJuOyeBcoNHd0wNc2wqcZGYRd+9y9w7yo4UAmFkD4O6+ZbCHm9kFwAUACxYsGNvIRURERA6wKTNNXIIEMFAiCJDq1acfd3/Q3Qdci9irzy3uXu/u9XV1daMIU0RERGTiTcdksAMYrNhOrFefETOzVWZ2S3Nz82geIyIiIjLhpmMyuA2oNbOBEsJ55KeQR1VP0N1Xu/sFlZWVo3mMiIiIyISbjsng4+S/ruN7N5pZjPzGkDUTEZSIiIjIZDQdk8G7AAcu3qf9fPJrBb8/2g/QNLGIiIhMF1NmN7GZfRRY2PO2DoiY2Rd73m909zsA3P1ZM7sJuNDM7gEe4LUTSB6ip+D0aLj7amB1fX39+aN9loiIiMhEmjLJIHAesHKftq/0/PoQcEev9ouBDeRLwJwJNAI3AlfscxTdiJjZKmDV4sWLR/soERERkQllOhh95Orr633NGi1BFBERkcnPzJ5w9/p926fjmkERERERGSYlgyOgDSQiIiIyXSgZHAHVGRQREZHpQsmgiIiIyEFMyaCIiIjIQUzJ4AhozaCIiIhMF0oGR0BrBkVERGS6UDIoIiIichBTMigiIiJyEFMyOAJaMygiIiLThZLBEdCaQREREZkulAyKiIiIHMSUDIqIiIgcxJQMioiIiBzElAyKiIiIHMSUDI6AdhOLiIjIdKFkcAS0m1hERESmCyWDIiIiIgex0HA7mlkQiLp7R6+2KuA8YAZwp7s/O/YhioiIiMh4GXYyCNwMnAi8HsDMwsDvgGU91y81sze5+9NjG6KIiIiIjJdSponfAtzX6/1fkU8EPwOcBOwE/mHsQhMRERGR8VbKyOAcYH2v92cCz7v7vwOY2S3AJ8cwNhEREREZZ6WMDBoQ7PX+ZOA3vd5vB2aOQUwiIiIicoCUkgyuB94BYGZvJj9S2DsZnAscFIX3VGdQREREpotSksH/AN5jZs8BPwV2AT/vdf0EYN0YxjZpqc6giIiITBelJIP/BlwJpIGngPcVysyYWQ35ncYPjHmEIiIiIjJuhr2BxN0d+ErPa99ru9F6QREREZEpZ9gjg2b2azM7bYjrp5jZr8cmLBERERE5EEqZJj4ZmDXE9ZnAylFFM0HM7P+a2WYzazGzrWZ2vZlFJjouERERkfE2lmcTV5FfTzgVfQs40t0rgGN6Xv84sSGJiIiIjL8h1wya2dHAsb2a3mpmA90zA/g0sHYMYztg3L133AbkgCUTFI6IiIjIAbO/DSTvI7+DGMDJnzAy2CkjrcBnxyiuAZnZZcAbgeXAYcBGdz90kL4B4HPk4z0UaAB+BFzh7u0D9P8H4ItAEtiNjtYTERGRg8D+ksHbgAfJj5b9Grga+OU+fRxoA9a6e2qM49vX1UAT8CT5aemhXEc+Ob0XuAZY2vP+ODM73d1zvTu7+z8D/2xmS4GPkD9RRURERGRaGzIZdPeNwEYAM/sb4Lfuvn6oe8bZ4e7+ak88zwFlA3Uys6OAi4B73P3sXu3rgRuADwE/GOhed3/BzP4E3AGcMrbhi4iIiEwupWwg+T756dMBmVnFIOsJx0whERyGD5Mfzbx+n/ZbgQ7gnP3cHwaOKC06ERERkamnlGTwGmDNENcfB74xunDGzArym0Ae693YM439dM91AMys0sw+bmZVlnc0+bWDvY/aExEREZmWSkkG3wH8ZIjrPwHOGF04Y2Yu0OjuA5W62QrU9qoj6ORHCl8lvwnmv8gfq3fRgQhUREREZCKVMq07H3hliOuv9vSZDBIMXvMw1atPl7u3AKcP98FmdgFwAcCCBQtGE6OIiIjIhCtlZLALmDPE9dnkp2Yngw4gOsi1WK8+JXP3W9y93t3r6+rqRhSciIiIyGRRSjL4NPCBgY5pM7Mw8EHgmbEKbJS2kZ8KHighnEd+CrlrpA83s1Vmdktzc/OIAxQRERGZDEpJBr8FHAXcb2b1ZhYxs7CZ1QP3A8t6+kwGj5P/2o7v3WhmMfInqgy1EWa/3H21u19QWVk5mseIiIiITLhhJ4Pu/hPg68BpwB/JT7N29Pz+dOBf3P2u8QhyBO4ivzHk4n3azye/VvD7o3m4RgZFRERkujB3L+0GsxXkd98u7mn6M/ADd398jGMb6LM/CizseXsRECFf8gbyR9Pd0avvjcCF5E8geYDXTiD5PXDqvieQjER9fb2vWTOqQUYRERGRA8LMnnD3+n3bSy4S3ZP0jXviN4jzgJX7tH2l59eHyJ8aUnAxsIH8zt8zgUbgRvJnE0+WjS4iIiIiE2pcTwwZa+5+cgl9s+RHDa/ZX99SmdkqYNXixYv321dERERkMispGew5bu69wAlANf3XHLq7nzdGsU1a7r4aWF1fX3/+RMciIiIiMhrDTgbNbAbwG+D15M/99Z5f6fV7Jz+VKyIiIiJTQCmlZb4KHAn8LXA4+eTvHeQ3ZvyQ/DrCmrEOcDLSbmIRERGZLkpJBs8Ebnf3/wBaetqy7v6iu58DdJIvPTPtqc6giIiITBelJIOzeW0XcXfPr7Fe1/8LePdYBCUiIiIiB0YpyWATkOz5fSuQAeb3up4hv6lk2tM0sYiIiEwXpSSDfyZ/5Bw9dfqeAj5uZlEzSwAfA14d+xAnH00Ti4iIyHRRSjL4C+CvzCza8/5a8iVmmoBdQD1w3diGJyIiIiLjqZQ6g1cD33T3NIC7/8jMuskfTZcFfjyJziYWERERkWEYdjLo+UOM0/u03QPcM9ZBiYiIiMiBMexpYjP7tZmdNsT1U8zs12MT1uSmDSQiIiIyXZSyZvBkYNYQ12cCK0cVzRShDSQiIiIyXZSSDO5PFftMI4uIiIjI5DbkmkEzOxo4tlfTW81soHtmAJ8G1o5hbCIiIiIyzva3geR9wJU9v3fgkz2vgbQCnx2juERERETkANhfMngb8CBgwK/Jl5f55T59HGgD1rp7aozjm5TMbBWwavHixRMdioiIiMioWL5izDA6mp0LPOTuG8Y1oimkvr7e16xZM9FhiIiIiOyXmT3h7vX7tpdSZ/B7YxuSiIiIiEy0Uk4gwcySwF8DS4Aa8tPHvbm7nzdGsU16wx1VFREREZmshp0MmtnxwE+B2iG6OXDQJINdXV10d3cTCpWUU4uIiIhMGqXUGbwWiAAfAGrdPTDAKzg+YU5eO3bsIJfLTXQYIiIiIiNSSjK4HLjG3X/s7k3jFdBUYmZ0dnaye/fuiQ5FREREZERKSQZbAGU9+0gkEuzZs4eWlpaJDkVERESkZKUkg/cA7xivQKYSM1tlZre0tLRgZiQSCRoaGrShRERERKacUpLBLwAzzexGMzvczPbdSXzQcPfV7n5BRUUFAIFAAHcnk8lMcGQiIiIipSllG+xe8ruFjyd/DjED5IPu7gfl1tpCMhiJRCY6FBEREZFhKyVxu518MigDCAaDdHR0kEwmJzoUERERkWEr5QSSj49jHFNS7zWCoVCIzs7OCYxGREREpHSlrBmclswsama3mtmrZtZqZn82s4uGc+/Wlgxd3fkag6FQiHQ6TTabHdd4RURERMZSycmgmZ1iZjeY2U97XjeY2SnjEdwBEgJ2AG8HKskX1f6imX1gfze2d2X5xwfWFxNCQJtIREREZEoZdjJoZgEzuwP4FXAh8M6e14XAr8zs9qm4w9jd2939S+7+srvn3P1p4D7gLfu7d1ZZmEc2tBQTQjMjnU6Pf9AiIiIiY6SUkcG/Az4C/Bg4Foj3vI4FftRz7dKxDrA3M7vMzO7umdJ1M9swRN+AmV1iZuvMLGVmm83sGjMbcoeHmYWBtwLP7C+eqniI/+/U+TyyoYXLHliPBfKbSERERESmilJ2E38c+IW7f3Cf9meAD5tZNfAJ4Joxim0gVwNNwJNA1X76Xgd8Fri3J6alPe+PM7PT3X2wA4W/BbSS3z29X+99fS2pTI4bHt7K87vSvH6m4+4Dld0RERERmXRKGYelj8oAACAASURBVBlcBKwe4vrqnj7j6XB3r3H3twHbButkZkcBFwH3uPtZ7n6ru19KfuTyFOBDg9x3LfAm4Ax37xpuUGccOQOAF3Z1ksvl6O7uHvYXJCIiIjKRSkkG24FZQ1yf3dNn3Lj7q8Ps+mHAgOv3ab8V6ADO2fcGM7seeBtwmrs3lhJXZTzEIZVR1u5s10kkIiIiMqWUkgw+DFzYM+rWh5ktAz4D/HasAhulFUAOeKx3o7ungKd7rheZ2Q3A6cCp7t4wkg9cNjvB8zs6CAQCqjcoIiIiU0YpyeAVQBR4qmcTxz/1vH4MPAWEgSvHI8gRmAs0uvtAW3u3ArVmFgEws4Xkp5QXA+vNrK3n9T+lfOBRs5M0tmfYk0abSERERGTKKOUEkmfNbCXwb8DZPa+CR4DPufuzYxzfSCWAwWq8pHr16XL3jeSnlIfFzC4ALgCYN29esf2oWQkA1jWkqIo4uVyOQOCgr+ktIiIik1xJ2Yq7r3H3N5NfO3hiz2uWu7/F3Z8YjwBHqIP8KOZAYr36lMzdb3H3enevnzFjRrF9cW2ccMB4YWeH1g2KiIjIlFFKaZminnV1I1pbd4BsA5aZWXSAqeJ55KeQh71beF9mtgpYtXDhwmJbJBTgiLo4z+1ox6yadDpNNDpYPioiIiIyOZRyAskHzWzQ2ntm9j0z+6uxCWvUHif/tR3fu9HMYuSLZK8ZzcPdfbW7X1BRUdGnfdnsJC/u6oRAkPb2cd1YLSIiIjImSpkmvpD8Dt3BZMlvxJgM7gIcuHif9vPJrxX8/mgebmarzOyWlpaWPu1HzU6Q6s6xuaWbtrY2TRWLiIjIpFfKNPFS8kfRDeYpYNXowhmamX0UKMzN1gERM/tiz/uN7n4HFDe73ES+FM49wAO8dgLJQ8APRhOHu68GVh999NHn925fNit/0t0LOzuZf3iC9vZ2qqr2d1CKiIiIyMQpJRlMkh/9G4wD5aMLZ7/OA1bu0/aVnl8fAu7o1X4xsIH8zt8zgUbgRuCKIY6iG5V5lRGqYiGe39nOqmXV7Nmzh8rKSh1NJyIiIpNWKcngeuAt5M/uHchbgE2jjmgI7n5yCX2z5M8kHvOzkgfaQNLTzrLZCdbu6CAYDJJKpUilUsTj8bEOQURERGRMlLJm8F7g/WZ23r4XzOwTwPuBe8YqsMlssA0kkC8+vaEpRXs6SygUorm5eQIiFBERERmeUkYG/xl4D3CLmV1C/lg3gGOAZcCLwNVjG97Us2xWAgde2NXB8kPKaGtro7u7m1BoRFV8RERERMbVsEcG3b0VeDNwMzAH+Oue11zg34GT3L1l8CdMH4PtJob8GcUAz+9oL64VVJkZERERmaxKPYGk2d0/DdSSP4VkFlDr7he6+97xCHAyGmia2N0BKI+GWFAd5fkd+QNOIpEIe/bsKV4XERERmUxGdHiu5zX0vA7aLCeXy9HR0UFHRwetra3FhO/1s5M8s72NTDZHKBQik8mQSqX28zQRERGRA29EyaDkRSIRFi5cyKJFi0gkEmSz+co7py+ppiWV5eFX85tHwuEwDQ0NGh0UERGRSUfJ4AgU1gy2trYSDocJBAIkEgm6u7sBWLGgnFllYVY/vxvIJ43pdJrW1taJDFtERESkHyWDI1BYM1hZWVlsi8VixZHBYMA4c1kNj21qZXtLFwDxeJyGhoZiwigiIiIyGSgZHCORSKTPSSNnLqsB4P61+dHBQCD/rW5qajrwwYmIiIgMQsngGAmHw8Bru4rnVEQ4fkE5P127m2wu3xaLxWhubqazs3PC4hQRERHpTcngCBTWDPY+XcTMiEajxaligHe/vpZdbRke29RS7BOJRNi1axe53LgcjywiIiJSkpKOxTCzJPlC00uAGsD26eLu3u+4uunG3VcDq+vr68/v3Z5IJGhubi6eNvKWwyqoiodY/fxu3nRofn1hOBymo6ODxsZG6urq+kwti4iIiBxow04Gzex44KfkC04PxoFpnwwOJhaLsWfPnuL7cDDAu5bO4K6nd9HUkWFGIj+VHI/H2bt3L5FIhKqqqokKV0RERKSkaeJrgQjwAfKnjgQGeAXHJ8ypobBusLe/XFZDNgc/7dlIAvnp4kQiQUNDAx0dHQcyRBEREZE+SkkGlwPXuPuP3V1bYgcQDocxsz7FpQ+dEWPF/HJ++OQu2tKvrScMBAJEo1G2b99OV1fXRIQrIiIiUlIy2ALs3m+vg1hhE8m+tQQ/9ea5NKey3LFmZ5/2UChEIBBgx44dfTaeiIiIiBwopSSD9wDvGK9AppKBdhMXxOPxfsngkTMTvP111dz19C52tvYdBYxGo3R1dem4OhEREZkQpSSDXwBmmtmNZna4HcTbYAc6gaQgHo8PWDbmk2+aC8Atf9je71oikaClpYW9e/eOfbAiIiIiQyglGdwLHA98Gvgz0G1m2X1eB/1Za4V1g/uaUxHh/cfW8bN1Tfy5of+mkcKGkra2tgMRpoiIiAhQWp3B28mXjpEhhEKh4iaSfZPCj9XP4qfP7+am323j+vce3ud6IBAgFouxY8cODjnkEGKx2IEOXURERA5Cw04G3f3j4xjHtGFmxONxMplMv1Iz5dEQHz9+Nv/22608+EozpyzuW2MwFAqRy+XYsmULFRUVVFdXD1iuRkRERGSs6Di6cVBIBgdy1htqOXJmgq//ahNbm9P9rkciEeLxOK2trWzYsIGGhgbS6bQ2l4iIiMi4GFEyaGZlZnaImS3Y9zXWAU5F0Wh00OQtHAzw1TMOBeBL/7Oeru7+m00Ko4uJRILW1lY2bdrExo0b2bNnD11dXTrXWERERMZMScmgmX3IzJ4DmoGNwPoBXge9WCxGIBAYNGmbWxnl8rctYN2uTr71u62DPsfMiMViJJNJAoEAu3fvZtOmTbz66qusX7+ebdu20dLSolFDERERGbFhJ4Nm9l7gB+TXGd4MGPBD4G4gAzwBXDUOMU46Q9UZhPxmkMrKStLp/tPABSsPr+JDx9Xx42ca+d8/7xm0X0EoFCKRSBRf4XCYrq4udu7cyc6dO1W0WkREREaklJHBzwMvAMcCV/S0fdfdPwTUA68Dnh7b8CanoeoMFpSXl+93OvfTJ83j9bOTfP1/N7F2R3tJMQQCASKRCMlkkvb2djZv3jxk8ikiIiIykFKSwaOB77l7CihkOUEAd38OuAW4bGzDm7qi0eiAR9P1FgoaX3vXoVTFQ1zy368MWH9wOOLxOACbNm1SnUIREREpSSnJYJDXzibu7Pm199DYi8DrxyKo6aKqqoqurq4h+9SVRbjxrMUkwgE+d+/LvLK7c8j+g4lEIsRiMbZv305ra+uIniEiIiIHn1KSwS3AQgB37wR2Act7XX8dUNpc5yRhZh8ws9+ZWZuZbRir5yaTSYD9bvCYUxHlxrOWEA4an7v3ZTbuSY3o84LBIPF4nB07duhoOxERERmWUpLBR4DTe72/D7jYzK4wsy8DnwEeHLvQDqg9wLeAy8fyocFgkLKysv2ODgIcUhXlhvctIefwf+7+M49ubBnRZwYCAeLxOLt27aKpqUllaERERGRIpSSD/xd40MziPe8vJz81/GXyG0peIb/JZMpx91+6+53ky+WMqcrKyiHXDfZ26IwYN79/CbXJMH/336/w7Ue3k82VXjYmEAiQSCTYvXs3r7zyClu2bGHv3r0qXi0iIiL9DDsZdPfH3f0fe6aIcfcGdz+W/O7iNwDHuPvmcYoTADO7zMzuNrNXzcyHmtI1s4CZXWJm68wsZWabzewaM0uOZ4z7isVihEKhYZd+mV8V49YPvI4zls7gu4/t4NL/foXd7QOfZjKUQCBAMpkkkUiQzWaLNQo3b95Ma2urRgxFREQEGIPj6Nz9GXd/3t0PRHZxNXAq+VHI/RXnuw64FlgLXES+HuJngdVmdsCO4TOzYW0k6S0WDnD56Qu47LQF/GlbGx+8fS0/eHInmWzp32IzIxwOE4/HSSaTuDs7duxg/fr17NmzRyOFIiIiB7lQqTeY2V8AbwdmAde4+zozKwPeCDzj7uO5c+Fwd3+1J47ngLJBYjyKfAJ4j7uf3at9PXAD8CHyBbQPiLKyMhobG3F3zGxY95gZq46q4Zi5Sf7tt1v51u+2cd/zu7n4rYdw4qEVI44lHA4TDofJ5XI0NjbS3t7OrFmzCIfDI36miIiITF2lnEASNLO7gN8A/wh8Apjbc7kb+C/g02MeYS+FRHAYPkz+hJTr92m/FegAzhnLuPYnHA6TSCTIZEqf7l1QHeOa9xzOv65aRC4Hl973Chfe8xLPbBtdPcHCNHJXVxcbN27U1LGIiMhBqpSRwS8AZwOXAj8jfxoJAO6eMrN7gXeRn8qdaCvIF8Z+rHdjT5xP91wvMrMgEO55mZnF8t19zI70qKqqYvv27UQikRHd/+bDKlkxv5x7n23k9jU7+T8/fokTFpRzwZvmsnRWYsRxxWIxstksO3bsKE4pR6NR4vE4oVCo+AoGgyP+DBEREZm8SkkGPwbc7u7/ZmY1A1x/gXwyOBnMBRoHSea2AieZWcTdCwv5Pgr8R68+neR3Fh86VgHF43ECgQC5XI5AYGRLFiOhAB88bibvfn0NP3mmkf98Yid/e9eLnH1MHZ980xySkZElbMFgsFgTMZvNkkql+p1kEgwGSSQSJJNJIpEI4XB42FPeIiIiMnmVkpUcCvxhiOt7gepRRTN2EsBgo3qpXn0AcPfb3N32eR060M1mdoGZrTGzNQ0NDcMOKBAIUFlZOSbnB8fDQc5ZPoufnHsUZx9dy0/+1MA5//kCv1/fPOpnB4NBIpEIiUSizyscDtPR0cGOHTvYtGkT27ZtK2lTzIGSy+Xo6Ohg586dOolFRERkGEpJBluBGUNcXwwMPzsaXx1AdJBrsV59Subut7h7vbvX19XVlXRveXn5mK7LS0aDXHryfP7f+48gEQny96tf5Uv/s56mjtLXJu5PIBAgGo0Wk8N0Os3GjRvZu3dvcUeyu5PJZA54PUN3J5VK0dDQwPr169m6dSvt7e1s376dXbt2aS2kiIjIEEqZJv4dcI6Z/cu+F8ysmvyGkp+NVWCjtA1YZmbRAaaK55GfQh7xsJaZrQJWLV68uKT7CucHZzKZMd29+4Y5SW778Ov4zyd2cdtjO3hsUysXvXUeZy6dMW5TubFYjFwuR0NDAy0tLZgZXV1duDvuTjAYpLy8nPLyciKRyKBT44PtsC4klt3d3cVnujuBQIBAIFBcw9je3k5zczOZTIZgMEg0Gi1+lrvT0tJCKpVi9uzZI16vKSIiMp3ZcEdwzKyefEL4B+A28mvs/o78CNs/ADOBFe6+dlwi7R/Pc0DZQNO5ZvZV8iek/IW7P9yrPQbsBn7r7meMNob6+npfs2ZNSfe0trayY8eO4hq9sbahKcU3fr2JP21rZ/khZXzg2Jm88ZCyEa8nHI5MJoOZEQwGi4ldLpejq6urOCpX2JwSCoWK93R3d5PL5QiHw8RiMaLRKMFgkPb2djo6OooJYG8DJY6F+waTTqfp7u4u1lqMRqP7vUdERGS6MbMn3L2+X3sp03lmdibwbfI1BgGcfAmXXcDH3P0XYxDrcGMZKhl8A/An4N596gxeRL7O4Efd/T9HG8NIksFcLsf69ev7jGCNtZw79z23m39/ZBut6SyhgHH0nCQnHVbBmUtrqIyXXF5y9DHlcuRyuWJyVxjhMzOy2SzZbLZ4bF9hB/NYfn/cne7u7uJIYyAQYObMmZSVDViqctg6OzuB/EipNtSIiMhkNibJYM+DosDbgKXkE8GXgJ+7+4jW4JX42R8FFva8vQiIANf0vN/o7nf06nsjcCFwL/BAT7yfBX4PnDqaE1N6TROf/9JLL5V8/549e2hsbBy30cGCru4cz25v59FNLTy2sZWXGjuJhowzl9bwoeNmckjVYMsqp79sNktnZyfV1dXU1NQQCASKaw/b2tqIxWIkk8lBE9LCEX979+7FzIjFYtTU1BCPx5UUiojIpDRmyeBEMrMHgZWDXH7I3U/u1TcIXAxcQH4ndCNwF3CFu4+uYnOPkYwMQn6UbOPGjQSDweK06YHwyu5O7nxqFz9ft4dsznnb66q56C3zqEkenKePuDudnZ3FjTHNzc1ks1mCwWDx1+rqasrLy/usQ0ylUuzcuRN3L44IZjIZurq6CIfDxanoSCRSrNGoBFFERCbatEgGJ4vRjgzC+K8dHEpje4YfPb2Lu55qIBoK8Om3zOXdR9UQOEgTlq6uLrLZbL+p+2w2W9wZ3TuZy+VyxGKxARP5bDZLJpMpTokX7guHw0QiEcrKyigrKxvX5LCwVjMWi+2/s4iIHDRGlAya2a9L/Bx399NKDW6qGunIIORHmLZs2UI2m52wXa6b9qT4l19v5smtbRwzN8lZR9exaEaM+VVRIqHxWc94MHJ3crlcMVFMJpPU1dWN+XnQ7k5zczMNDQ24O7Nnz6aiYuTnWIuIyPQy0mQwB2SA4ZZhcXc/aP72GU0yCJBKpdi8eTOJRGLCphHdnftfaOJbv9tKSyoLQNBgflWMty6q5B1HVrOoJj4hsU1XqVSKXC5HbW0tlZWVg/7ZD1Z2ZyBdXV3s2rWLzs5O4vH8n1dHRwczZsygpqZm0OcUNu+o7I6IyPQ30mSwkATeT76UzE9Hs/FiuhiLaeKCHTt2FDcs9Dx7QhLDru4cG/ekWN+UYkNTirU7O3hicytZhyW1cU47ooo3H1rJohrtmh0LuVyuuF6xrq6umMBBPrHbvXs3qVSKWbNmkUj0P3s6m83S1dVFKpWivb2dzs5OQqEQ0ehrm4Lcnfb2dioqKqiuru5zxnQ6naalpYWWlhbcnXnz5vWJYSCpVIrW1lZmzJihsjwiIlPQSJPBOvJnEn8cWEa+hMztwHfd/cXxCXXqGO3IIOTr7e3YsaM4QlP48xgoATjQmjoy/O9Le/nFuiae35nfLF6XDHPiwgre/rpqls8vn+AIp76uri4ymQzl5eVUVVXR2trK3r17i5uL0uk05eXl1NbWEgwG6ezspKWlhba2tuLIYaF+42BJemdnZ/HnqlAPMpPJFE+VyeVypNNpDjnkkAETwlwuR1NTE3v27AHyxdPnzp075tPcIiIyvka9gcTMjid/ysgHgQrgMeA7wJ1jtTt3qhmLZHBfuVyODRs2EA6HJ9XoS0NbF3/c2MofNrbw+KZW2rqyrJhfzqfePJcjZ0584jqVuTvpdJpsNksgEOhTs7BwrZD45XK54vnRIxmhLdR73HfzS3d3N+l0mvnz5xdHqbu7u4vTz4Wi3WZWjGfevHl9RiJFRGRyG8s6gzHgbOBvgFPIn0DyqbEo4jzVjEcyCNDS0sLOnTsnZKfxcKS7c9z7bCPfe3wHzakspyyuYkldnFzOyTlEgsZxh5SxbFaSYEBTymOhcJLLeBUqh9eSv/Lycjo6OopFwCORSL9RwEwmQyaTYdasWcWd1Vo+ICIyuQ2WDJZc5M7dU8D3zWwDkANOBxaNOsIpZKRnEw9XWVkZu3fvpru7+4DWIRyuaCjAh46byV8uq+GHT+3irqd28ZuX9/brVxELcvyCcpYfUs4RdQkW1cSIapfyiIxnElhQ+Fnr6Ojot/5wX+FwGDNjx44dxXsTiQRlZWXEYrFJNaotIiJDK/U4ujnAueTXEC4BtpFfQ/j/3H3zeAQ4mY3XyCBM/tHB3rpz+TOEA2YEDFrTWR7b1MqjG1t4dGMLTR35EaaAwYKqGHMqI1THQ1QnQtQkwiyujbN0ZoJkVAnEVJXL5YrnTZsZ8Xic8vJyysvLNWIoIjJJjHia2MzCwHvITwu/HcgC95HfXfzzg3l38XgmgxN1SslYy7mzrbmLlxo7ebnntau1iz2d3ezp6CaT69nYABw6I8aC6ihd3U5HJktnJkdZJMjiujhLauMsro1z2IyYaiBOcu5enEaOx+PMmjWr3zTzgZj2FhGRvka6m/gG4K+BauBZ4LvAf7p703gFOpWMZzIIE3tKyYHg7uzt7ObFhk7W7mhn7c4OtjV3EQsHSIQDxMIBmju7eXl3J+nu/M9pMACHVsdYUhfn8Jo4cyoizCrPv9zhzw0d/Lkhn3SWRYIcMy/JMXPLmFsxsg0XMjqpVAqAWbNmkUwmyWQyxR3TuVyOZDJJWVkZ8Xh8Sv+jR0RkKhhN0elO4F7gyWF8jrv7dSOOcooYyzqDQ3F3Nm7cSCAQOKj/oszmnC3NaV7qSfJeaujkpcZOGtszg94zrzJCSypLazpfSLs2GWJBVYzZFRFml0dYUhfnLxYNXvBZxk53dzepVIp4PE4qlSqWtCmc6ZzJ5P8cC+c6JxIJwuEwgUCAQCCAmRV3VXd0dNDa2oq7F89/DofD5HI5uru7yWQymBkVFRXE43GNPIqI9DKaZLAU7u4HzcKv8R4ZhOk/OjgaLaludrVl2Nnaxc7WLnIOi2vzU8rJaJCcO+t3p/jTtjae29HOtuYutrd00diewYHTllRx2WkLSET6/shu2ZumOhEiGTlofpTHnbsXN0QNdRpK4Vzn3vcVksFCQfZCWZ1Cbc5cLoeZFRNHoFhHsbKykvLy8hGX4hERmU5Gupv4lHGKR4YpmUwSCoXIZrPaobmPiliIiliIxbUDn5wRMOPw2jiH18Y56+i6Ynsmm+POpxq4+Q/beHV3iq+feRjzq6I8uaWN763ZwZrNbSQjAf5yWQ3vP6aOuZWqpTdaheLYQwkGg4P+jA90NF8gEBj0mZFIhFwuR3NzM3v27CEcDlNZWUkymSQYDJLL5XD3PklmIaEcyWhiIdnt6uqio6OD9vZ2ysrKqK6u3u9/t4WENhgMjvlIZnd3d/H52WyWcDis2pAi0k/JdQblNQdiZBCgubmZXbt2aXRwjD2+qYUrfraB7pxzaHWM53d2UJMIcfYxdWxoSvG/L+3BHU5YUEFlPEjO8xtiUhlnT2eGPR3d7OnsxoCqeIjKeIiqeIjDZsQ4cmaCI2cmmFepEanJoJCo7e+8Z3cnFouRSCSK6xhDoVAxSSuc1tLe3k5bW1ufU4MKCss6urq6CAQC1NbWDriruquri5aWFvbu3Vv87FAoVEzYotFosfh8JpMpJplmRm1t7ZBnmnd3d9PY2Ehra2ufPu5OXV3dkGdiH0i9T8YZSCqVmnQF+EUgP/vg7lPuXPcxKzotrzlQyWA2m2XDhg1Eo1GtgRpj21u6+NL/rKepo5uPLM/XTizUQmxo6+InzzTy65f2knUnYPnRxkjQmJEIU50IUR0P4cDezm5aUt3sbu9mQ1OquEu6PBrkiLo4R85KcGRdgqPnJqkrm1r/8ziYuHuf6erCFHUhSSucvlIYlew9qjiQbDZbTGgK6yQDgQCZTIbOzk7MrM9/14UTYrq7u/tNlweDweL6yK6uLmKxGLW1tX3ud3daWlpobGwsPrt3bIUzsZPJJDNnzizOOmQyGbLZ7IAFxsdDJpOhra2N5uZmzIx58+b1Wxfd1tbG9u3bCQaD1NbWUlZWpv//yaSQyWTYsmULAPPnz59Sa/qVDI6DA5UMAuzZs4fdu3dPijOLp5v9jU6UKpPN8eruFOt2dfDirg7W7erklcbOYoJ41KwEJy+uYuXhVRxSpSm7qaAwzTrSqdzeyV3vM6VH8zNXONcaXjtzuvBZ+9s8U9jlbWbF2pCF/w5isRiVlZUkEomS/pLLZrPFaWmgz5nZhXJDhbO10+k0ANFolEwmQzAY7JMQplIpNm/eXDwasZBQ19bWkkwmJ8Wopkx/7e3thMPhPqN/vRPBwkzCnDlzpszPpJLBMXSgdhP31t3dzYYNG4jFYvrX8RSUyeZ4ZXeKxza28OAre1m3qxOAJbVxTl5cxcmHV3JYTZzG9gxrNrfy+KZWWtPdnPWGOk5YqMLNsn+l1m4crH9hh3dhOrqiomLQn790Ok1bWxstLS3FJLD33ymBQIB4PF48exvySeK+o4/pdJpgMMjcuXPJ5XJs3ry5mEwW9N6VXldXV0wUx5u7FxPvdDpNZ2cnM2bM0D/Mp7k9e/awa9cuAoEAFRUVVFVVYWZs3boVoLj2tr29nZqaGmbMmDGR4Q6bksFxcCBHBgEaGxtpbm4mHh94w4RMHdtb0jz4cjMPvrKXZ7e3AzAjESqe1lIZCxIOGo3t3RxRF+ejy2exYkE56W4n3Z0j3Z3/i9x6pq6z7jS2ZWhoz9DQlqEmEeJdS2sIBZVEysgUppRjsRh1dXVEIpHi6F9XVxfNzc2k02nMbNB/pBZGVAOBwH7X/RVG/wojlYOtxSokZuXl5dTU1IzbtLa709nZSWNjY3G9aTAYLI50zp8/v99nF0ZAp9o6MumrtbWV7du3F9flFv4xY2aEQqE+f765XI6Ojg4OOeSQKfEPBCWD4+BAJ4OZTIYNGzZoI8k009CW4bev7uWZbe0sqYuzYn45S+riZHPOz9ft4T+f2MmmvemSn7ugOspn3zqPkw6t7NPu7mxr6WLdzg5e2NVBZybHB46tY2H1gRlpkamlkHztOzoYiUTGfK1U7+njoRTqTmazWWpqaqiqqhqzGZPCs3fv3k1HR8eA6yhTqRSRSIR58+YVvy+5XI5du3bR3NxMZWXluCaqMjwjWQLU3t7O1q1bSSQSfX6mCtUHBvo56+7upru7m3nz5k363fpKBsfBgU4GATZs2DBkCQ6ZfrI55/frm9nW0kUsFCAaChAN5Tct5HJODggANckwdWVhapNh1mxu5YaHt7J5b5oTF5azdFaSrc1ptu5Ns2lvuliMOxwwAgHozjrvWlrDJ06Yzazy1/7V25119qa62duZf2VzzvL55YQCGnGUiZfL5UilUgSDQerq/v/2zjw8srLK/5+TSlWqsqfTSe8L9ALd7KssytooiysuqAOD4yjKCDjMjDouP2VGRx0dxBlcARVEUUBBQQQEFGRfZV8aeu9Opzvp7Evt5/fHe2/6dnUlqXQnXUnX+TxPOzsNQwAAIABJREFUPZV673vvfe9JVd1vnXPe8zbtEDXxJ+vk26e/v590Oj0kaP0Z2/39/fT29g7VxBzpxh4MD2azWVpbW+nv76eysnIHoVpbWzulJhiMBT/fc6xCPJFI0NraSiwWo76+ftw8qapKR0cHnZ2dQ69FhObmZmpra/Pu4+fz+l7v1tZWKioqxvw/8380TabZ+vkwMTgBFEMMbt26lf7+/kn/68MoPqlMlt88187PnmhlIJVhRk2EOXUR5tZVsLSpkmUzKtm3MUpvIsN1T27hdy+0IwJLmyrpibuyOb5oDLJkeox/O3keB80yD7UxOUin00Mh62D5oGg0Sm1tLdFolFAoRG9vLx0dHUMhv9xSQ354upAf26pKf38/s2fPpqenh4GBgR3ChL5QBVdDs6KiglgsRiwWIxKJ7DAD3M+9TKVSO62e489I7+3tHVq+cU8IDT9NYGBggOrqaqLR6NB50+k0nZ2ddHV1DU2gKFQ8+QsphMPhoRqYdXV1BYlCVaWvr4+2tra8JZs6Ojpob2/fYZKRH8adNWsWNTU1O1yfn3oVPMbueLz9c1VXV9PU1DQpPcMmBieAYojBnp4etmzZYqFio2CS6SwiEA6N/Ot9c0+S655sZWN3goaYK5tTH3g0xMrZ2pfiBw+30Naf4u3Lp/GBQ5vpGEjR0pOktSfJ3PoKTllcv9OqLoaxp/ELgfv14IChcjvjFVlJp9MMDg5SXl4+Yi63Xy4onU4PCdDKykqi0Sg9PT1DeZKhUGjo79raWrLZ7NDyi+Xl5aTTacLhMI2NjTuJQj9fMZVKDYlQ3+sZCoV2EHOjXVNvby+dnZ1ks1nKysqGCpY3NDSgqmzbtg1wYtsX4bNnzx5yUviTbpLJ5JDnMBQK0d3dTUdHB7FYbOh/EAz5+zPZY7HYTkIqlUrR1tZGX18fFRUVJJNJ6urqmD59OmVlZXmFYND+QUGYSqVobW0lkUgQi8XGXVzH4/EhkVtbW5vXeeOXsfIL4O+pVZJMDE4AxRCDiUSCDRs2TIlEVWPvpD+Z4WdPtHLjs1vJBBasFECBWLiMU5fU8/bljRw0y8qAGJOD0QqO7yq+YBrrWPzajrmzpv1jJpPJHZZf9Al6Qf12lzKy/cPotwfv76FQiMbGxqFVePztfp1Jv6i5f+xc0Rws3J5buiiZTJJOp5k5cybpdJqurq68eaaqOmKxdF/M+uI3KGZ7e3spKysbmkXuT/CpqKigsrJyqPTacP8L39PZ2NhIZ2fn0BrpE4VfoD6bzRKNRqmpqRkSyP7/Pkg0GqW5uXlCJx95Kx2ZGBxviiEGs9ksq1atGtUz6Icw/A/13pqzYhSPtR1xXtkywMyaMLPqKpheFeaVLf384aUO7nu9k4FUltm1EVYsbeCt+zWwb2OMZDpLS0+SjV0JRFzIual69+rtGUYp4t+7g17PkT5Hfri5rKyMyspKUqkUyWRyaLtf+3JX7xVBkToek4v8AuyqSjabJRwO5xV6iUSCVCo1ohD0yWQyDA4O7uCZ3BP4ItcXtvmWvUwkEqTTaRobG6mrqxuqy9nb20symRwSx+FwmJqamjF5NFWV3t5e2traWLx4sYnB8aIYdQaDrFu3btRSDfF4nMrKSqqqqujs7Byq47WnanMZpc1AMsP9q7r402udPLWhl6xCQ6ycrsE0ud84tdEQS6bHePO+dZy1rJHqipG/pNMZRdFRw96GYeyMH7L27yH2Q2zy4OeZBtdM95dj9MWx//+LRqM0NjYSi8WGwuB+QfeamhqqqqqIRqOkUim2bt1KPB4nGo0yf/58E4PjTTE8gwBtbW309vYOK+xUlYGBARYsWDDkck4kEnR0dNDX1zfq6gSGMZ50DKT48+tdvLp1gJk1EebWVzCvvoJ0VnmjbZCV7YO8sqWfN9rjxMJlvG2/Bt590HQWT49RFrhRre+Mc+sL7fzxlQ76EhmmV4WZWRthRk2YcFkZWVUUqIqEeOcBjezXbKkUhmHsnfgh50gkMpSL6nt2k8nkDpOkgquozJkzx8TgeFMsMeiv2TlcqDgej1NdXU1zc/MO7b6reMuWLTstsWMYxebVrQP89vk27nmtk2RGiYSE+Q0VLGiI0h1P89SGPkJlcNKiehY0RGntdZNWtvQlyarLWRSBjoE0g6ksb5pfw3lHzuCwOXtm9qVhGMaeJp1OD+vhzVdn0cTgBFAsMZhIJFi/fn1eMeh7BRcuXDjstPZkMjk0kwoYCheUl5ebx9AoOt2DaR5c3c2ajjjrOt0D4KzljbxjeSONVSOXa+hNpLn1hXZu/FsbnYNp9m+O8c4Dp3Pa0gaqcmY598TTrN4WZ21HnNUdceKpDEfOq+Ho+bXUx0bOeUpnddh6ixM1WcEwDGN3MDE4ARRLDGazWVavXp13RvHg4CC1tbU0NTWNegx/BlgikRh6+GUNgAmZcm8Ye4pEOssdL2/jlhfaWb3NhaBPWVJPrLyMNR1x1nTEh5b/A6gMlxEqE3oTGQRYPqOS/Zorh0rrVEVCbOpOsLJtgJVtg2ztS1EbDdHkFfsOh8rY1p9iW3+KjoE08xsqWLG0gVOXNDC3Pv+sxVQmy+aeJKEyYVZtZIewuGEYxnhjYnAYRKQcuBw4D7eQw2+BT6lqfLR9iyUGAdavXw+ww4ytQryCo+Enp3Z0dDAwMGATTowpj6ry8pYBbntxG/e+3okAC6dFWTgtyr6NUfaZFmXhtBgzasKounD1o2t7eHxdzw6rtQCUCcyvj7K0Ocbcugq6BtNs7UvR3p8kkVamV7kVYOpj5bzY2j+07vTSphgNsXJEXMgmkc6yqTvB1t7U0ISaaHkZC6e58UTDZaSzSiarJDNZOgfSdAykae9PURUJ8e4DG3nXgdOpG8V7aRiGEcTE4DCIyJeB9wFnAEngNuBJVb1ktH2LKQbzTSIZHBwcKsK5u/ih6JFqQhnGVCOdVUIytrVK01mlN56mJ5GhuTpMLFx4SYrNPUn+/Honj67tIZHOkgVUIRxynsC5dRXMqYuQzsLqbYNDIet0VgmVCaEyCJeV0VBZTmNlOY1VYdZ1xnlqQx+RkPC2/adx+JxqYuEyYmG3VGEinWUglaU/mSGRdtdbHhLKy4RplWGWNsWojZqINIxSxMTgMIjIeuCzqvpr7/XbgJuBBlXdeS2uAMUUg319fbS2tg6Fiv0ClwsWLBi3moKtra3mHTSMSciqbYP85rk27nq1g0R67N/hM2siLG2KUVFexmAqSzydQRBO37+B05ZOozxkPwANY29kyotBEfk8cDhwBLAPsE5VFw7Ttwz4NPAJYCHQBtwEfFlV+wP96oFOYJmqvuq1NQFbgcWqumqkMRVTDOauRNLf309zczN1dXXjeg7zDhrG5KU/maG9P0U8lfVEXZZoeRlVkTIqIyEqysvIZJV0VklllK19SVa2DfLa1gHeaB8kk2XIo9gTT7O+K8GMmjAfPnwG71jeSDRc+ISyeCrL4+t7+OuqbpKZLCuWNnDsgloi5TYpzTAmC8OJwakUK/g60AE8A9SP0vcK4BLgVlw+4DLv9WEiskJV/XV7/FWruwL7duVsm5T4OYH++psVFRXU1taO6zkqKiqorq4mHo9P6LI9hmHsGlWR0E4zpEdi4bQoR8/P/z2hqjyytodfPL2FKx7YyDWPbeZt+zVw1vL8NRvTGWXVtkFe3NzP0xv7eGxdD/F0lpqKEKEy4b7Xu6ipCHHS4noqw2Vs6EqwoStBW1+K6VVh5tZHmFNXwf7Nlbx1v4ZRi4gn007sWojbMMafqfSpWqSqqwFE5EWgOl8nETkAuBi4RVXfG2hfA/wf8EHgBq+513uuA1q9v+tztk1KysrKiEQiZDIZkskkc+fOnRDvXUNDAxs2bDAxaBh7OSLC8fvUcfw+dTzf0sdvn2/ntpe28Zvn21kyPcaSphgDySyDqQx9iQxvbBscClE3VYU5c9k0Tlpcz6Gzq0HgqfW93PVaB/e81gnA3PoIi6fHOG5hLe39KTZ2JXi+pZ/fPt/OdU9u4ZPHzeLkxfWICAPJDA+t6eaBVd1uok1fiq5BN/N7Xn0FR86r4Yi51Rw5r2aPikNVdUsw1kaYVpl/kl4inaXCvKF7NarKdU9tobUnySeOnUXDMO+FqcSUCRMH8cVgvjCxiHwN+CJwgqo+GGiPAtuAB1T1zED7euAzqnqj9/qtwG+Y5DmDAO3t7bS1tdHQ0MDMmTMn7DwtLS0kEgkThIZRYvTE09yzspM/vtJBx0CKynCIWLiMykgZixpjHDCzigNnVTFjhPWlR5q0o6o8uq6H7z/UwpqOOAfMrGRmTYSH1nSTSCtNVWEWN8Vorg7TVBWmPCQ839LPs5v6GEg5L+QXV8znhEWjBYu2k0hnuf+NLgbTWY6cWzNs2Z8gG7ri3P1qJ3e92kFLT5IZNWGufM+SHfZNZ5Vv/2UDf3x5G8csqOWMZdM4fp86E4ZFpK0vxTWPbebExXUct3B8UqhSmSxfv3c9d7/mKhPURkNc8pa5nL5/w5RIp5ryOYNBRhGDdwMrgEpVTeRsexhYqqpNgbYvA2cDZwIp4PfAU5N9NjFsn0SycOHCcZs0ko94PM6GDRuIRCJ5S9Ykk0nC4eFvBoZhGCORzip3vtLB1Y9tJp1VTllcz4qlDRw8uypv7cV0RnlpSz/ffWAjr7UN8oFDm/in42aPmJ+4tTfJLS+0c9uL2+iKb68vObs2wpHznCisrQhRFyunvExY1T7Iyjb32NidQIAj59Vw3MJafvZkKxWhMr733sXMq4+SSGe57O61PLCqmxMX1fFS6wDt/SlqKkKctXwa5x85c9QyQJt7EqzrSHDY3OodBGQyneXOVzu4d2Unpy5p4F0HNu4V37Ut3Qkaq8ITJpb/8kYX//3n9fTEM5SXCd84ax+O36dwQTiQzHDf613MqAlzyGz3P+lLZPjCH1fz1IY+Ljh2FifsW8c379vAi639HDWvhn8/dT6zaif3yl6lJAZfAJpVdUaebTcB7wcqVDXptZUD32F7ncHfABep6uBo4yi2GMxkMiQSibzFp8ebeDzO5s2bUdWh2cV+wepwOEwmkyEWi034OAzD2HtRb33pQotvJ9NZvv9wCzc/18Z+TTFWLG2gYyBNx4ALK/clMwwkXZmdbf0psgrH71PH+w+ZTlN1hKc29PLE+l7+tqmX/mR2p+PPqnWzrg+aVcWKJQ0017gb/ettA3z61lWUh+Cbb9+XHzzcwjMb+7j0xLm8/5AmMlnl6Y293PFyB/e93klVJMT5R83gfQc37SRY+xIZrnuylZuebSOVVaojIU5aXMdpSxtY3RHnhqe30tafoiFWTudgmhMX1fH5U+fvFB5XVTZ2J3hmYx8vbO6nobKcI+fWcMjs6oInAj3f0sd3/7qJmooQh8yu4uDZ1SybUUlluGzcBOgrWwa46tEWHl/fy4zqMBceP5vTlo6fV60/keGKv27kj690sH9zJf960lwuv38Dq9rjfOsd+/KmBdtzZgdTGboGM8ys2e7MUFX+9Fon33+4hfb+FADhMuGg2VV0DqRZ3xXnC6fO54xljQBkVbn1hXZ++HAL4ZDw1TP24ch5k3fKQSmJwVVAWFXn59n2c5zoa1DVrtztBZ77AuACgPnz5x+xbt26XTnMlCSdTrN161b6+/sREUKhEE1NTUSjUdavX084HCYUKjyZ3TAMYzz466ou/uve9fQmMlSUu3qKDd6qMVWREJWRMmbURDhr2TRm1+0cFlZVBlNZehIZeuJpEmllQUPFiPmIq7YNcsktb9A5mCZUBl9asYC37T8tb78fPNTCo+t6mFUb4biFtUPFyfsSGa57agvdg2nOWDaNExfV88CqLu5/o4uBlBOnh86u5vyjZnDU/Bp+/bet/OiRzTTEyvn0CXNIpJV1na425ctbBobES32snL5EZmjJxANnVnHkPJdjuXxG1U6lg7Kq3PDMVn78SAtN1RFqKkK80T44VBA9JBANlxEtL6O8TEh7M9SzCgfOquLcI2Zw6OyqIUG1rjPOTc+28eLmfmbWRphXX8HcugqeWN/D/au6qY2GeO/BTTy8ppuVbYMcMLOSS94yhwNnVu2yKBxMZbjl+XZ++cxWeuJp/v7IGXz06FmUh4SeeJqLb3mDdZ1xvnnWviQyWe5d2cnDa9ykp8bKcg6eXc2BM6t4YFUXz2/uZ//mGJ86fg6pjPLkhl6e2tBL12CaL542P+8krA1dcf79D2tY1xnnojfP4ZxDmyalB7eUxOCYPIO7Q7E9g8VAVenqcjq6rq5uaC3jnp4eWltbqa7OO6/HMAxjQkmks6SzOq5erNFY2xHnigc2cs5hTaPmpD2xvoefPt7Kmo74DqvaHDq7mk+fMGeHGduJdJbH1vUwrTLMQbN2XIP+ta0DfOWutazvcllQoTKYW1fB4ukxDp9bw+Fzq5lfX0E8neW5Tf08vdEJmZVtTtzFwmUcOLOKJU0xFjVGmVcf5donW3lkbQ8nL67n86fOp7oiRF8iw4ub+1nZPsBA0s3kjqeyQwKzvExQXDi2azDNATMrOWt5Iw+t7uaRtT2Ey4RD51TT3p9iU3eCZMb9bz54WDMfOqyZqooQmaxy56sd/PiRFrYNpJleVc4hs6s5dE41Cxqi9CXcZKWeRJq2vhQtPUlauhO096eYWRNhcVOMxY0xEuksv362ja7BNEfNq+ETx85i+cwd7dY5kOKiW95gTYdbXKw+Ws7JS+rZZ1qUF1v7ea6ljy29zgN74XGzOXP5tDEvD9mfzPC1e9bxwKpuTlvawIqlDcysCTOzNkJNxc4/LNIZHVp//ZDZ1XnXXe8YSNGbyLCgYffq/W7uSXDFA5v45SffUjJicEw5g7t4/ncA71i8ePHHX3/99d051F6DqrJx40YymQyRyOTOmTAMwygmiXSWbf0pBtNZ9p0WHbN4HUxleL6lnxk1bhWbQoqEdw+meWZTH09t6OXFzf2s64yTzLj7f7hMuOQtczj74OljHksineUPL2/jV89spaUnSX2snLMPms7ZB08fmnGdySpb+1wOZXXFztGjgWSGe1Z28rdNffxtYx9tnoczSLS8jFm1EebURWisCtPSnWTVtsGh9cXfNL+Gj75p1k7iOUjHQIqbnm3jsDnVHDGvhvKyHa+1rS9JTUX5mOpr5pJV5bont/CTxzeTDcirWLiM+lg5ddFy6qIhuuNpVm/b/j+IhIT3HDSdc4+YQWNVmE3dCW54Zit3vLyNZEbZtzHKW5c2cNp+DcyqHdtkzle2DPCZ21eRSGd5+atnlowYHG028V9V9YzxGEcpegZHwp9oEixSnU6nh0LKhmEYxuQgnVU2diVY0xH31ufePc9TOqus3DrAoumx3ZoUoqps6k7S1pekuiJEdUU5NRUhqiL5Pb4dAykGktmCZoXvSboH02zqSdDak6S1N0lbX4rueJruwQxd8TTVkRBLm1zJppk1EW57aRt3v9pBeUg4eFY1T2/sJSTCmcunsW9jjPtWdvK8t9b5/PoKljTFWNpUycJpUbrjaTZ7XtN0Vjl6fi3H71PLtMowD6/p5v/duZb6WDnfedcijj9oUcmIwYOA54Bbc+oMXoyrM3ieqv5iN89vnsFhaGtro6ura+hDG4lESCQSRKNRE4SGYRiGMQwbuxJc+2QrT67v5bT9Gjjn0GaaqreHjjf3JLjv9S4Xvm8bpLV3e7ZbmUBTdRhV2NqXQoD9mmOsbBtkaVOMb79jEY1V4amfMygi5wELvJcXAxHc6iLglqa7PtD3SuAi3Aokf2T7CiQPA6cEViDZLcwzuDOZTIb29nZisRixWIxwOExvby+bN2+mqmrXk4MNwzAMw9hOTzzN+s4EDZXlzKiOUB4SVJU32gd5aE0Pj6zpZm59BZ85eR6V3kpFe4MYvB84cZjND6jqSYG+IeCfcbN+FwLtwI24tYn7xmtMJgYLp729nc7OTqqqhs/nMAzDMAxj4pjyaxMHxV4BfTM4r+Hlo/XdFQJh4ok4/F5JY2MjqVSKgYEBq0doGIZhGJMIWydnF1DV21X1grq68VnephQQEZqbmwmFQqRSO88UMwzDMAyjOJgYNPYYoVCIxsZGksndLvFoGIZhGMY4YWJwFxCRd4jIVd3d3cUeypTDQsSGYRiGMbkwMbgLWJh41ykvLycajVqo2DAMwzAmCSYGjT1ObW2tiUHDMAzDmCSYGNwFLEy8e8RiMaZKSSPDMAzD2NsxMbgLWJh494hEIoTDYdLpdLGHYhiGYRglj4lBoyjU1dVZqNgwDMMwJgEmBo2iUFlZSTY7LqsCGoZhGIaxG5gYNIpCJBIhFAqZIDQMwzCMImNicBewCSS7j4hQU1NjBagNwzAMo8iYGNwFbALJ+FBdXW2eQcMwDMMoMiYGjaJRUVGBiJggNAzDMIwiYmLQKBplZWU0NTUxODhY7KEYhmEYRsliYtAoKjU1NVRVVZFIJIo9FMMwDMMoSUwM7gI2gWT8EBGam5vJZrNkMpliD8cwDMMwSg4Tg7uATSAZX8rLy2lubmZwcNCWqTMMwzCMPYyJQWNSUFNTQ21tLfF4vNhDMQzDMIySwsSgMWloaGgwz6BhGIZh7GFMDBqTBhEp9hAMwzAMo+QwMWgYhmEYhlHCmBg0DMMwDMMoYUwMGoZhGIZhlDAmBncBqzNoGIZhGMbegonBXcDqDBqGYRiGsbdgYtAwDMMwDKOEMTFoGIZhGIZRwpgYNAzDMAzDKGFKXgyKyAdE5CER6RORtcUej2EYhmEYxp6k5MUg0Al8D/hisQdiGIZhGIaxpykv9gCKjareAyAi7y72WAzDMAzDMPY0RfcMisjnReRmEVktIjpSqFZEykTkUhF5VUTiIrJBRC4Xkao9OGTDMAzDMIy9hqKLQeDrwCnAKlzIdiSuAL4DvAxcDNwMXALcLiI7XIuI/NoTl8M9Thr3KzEMwzAMw5hiTIYw8SJVXQ0gIi8C1fk6icgBOAF4i6q+N9C+Bvg/4IPADYFdPg5cNMJ5bfkQwzAMwzBKnqJ7Bn0hWAAfAgT4bk771cAAcG7OcXtVtX2ER2r3R28YhmEYhjG1mQyewUI5CsgCTwQbVTUuIs9628eMiISAsPcQEYm6w2piN8drGIZhGIYx6Sm6Z3AMzAbahxFpm4DpIhLZheOeBwwCNwHzvb9f2+VRGoZhGIZhTCGmkmewEhjOWxcP9EmO5aCqei1wbaH9ReQC4ALvZcLLczRGZjrQXuxBTBHMVoVhdioMs1PhmK0Kw+xUOJPRVgvyNU4lMTgANA+zLRroM6Go6lXAVQAi8pSqHjnR55zqmJ0Kx2xVGGanwjA7FY7ZqjDMToUzlWw1lcLELbhQcEWebXNwIeQxeQUNwzAMwzBKnakkBp/EjffoYKM34eNQ4KliDMowDMMwDGMqM5XE4I2AAv+c0/5xXK7gL/f4iLxwsTEqZqfCMVsVhtmpMMxOhWO2KgyzU+FMGVuJqhZ3ACLnsT2h8WIgAlzuvV6nqtcH+l6JKyR9K/BHYBluBZKHgVNUNbunxm0YhmEYhrE3MBnE4P3AicNsfkBVTwr0DeE8gxcAC3GzdG4EvqyqfRM6UMMwDMMwjL2QootBwzAMwzAMo3hMpZzBoiMiZSJyqYi8KiJxEdkgIpeLSFWxx1YMRGSpiPyniDwmIm0i0isiz4rIF/PZRET2E5HfiUiniPSLyIMickoxxl5sRKRSRFaLiIrI9/JsL2lbicg0EfkfEXnD+6y1ichfROQtOf3eJCL3eu+9HhG5S0QOLda49yQiUi0iXxCRF7zrbxeRR0TkIyIiOX1Lwk4i8nkRuTnw2Vo7Sv+C7SIis0Xk5957cVBEnhKR90/IhUwwhdpJRKIi8nER+b2IrPWue7WI/EpElg2zT4V3X1gjIgkRWSUiXxKR8IRe1AQw1vdTzr7/7e2TN2o52exknsExICL/i8tRvBW4E5ezeDHwILCi1HIWReSbwKeA24DHgBRwMvAB4HngGFUd9Pouwi0lmMatL92Nm/xzIHCGqt67xy+giIjI/wCfAKqB76vqRYFtJW0rEVkA3I+zzU+AlUAdcDBwt6r+2ut3jNdvE+AL6otw9UiPU9UX9ujA9yAiUgY8ABwHXIf7/FXi1nA/GviWqn7O61sydhIRBTqAZ4AjgB5VXThM34LtIiLTcBUrmoHvABuBD+NSnD6qqj+bgMuZMAq1k4jsD7wCPAT8CVfibV/gQqAKOF1V/5Kzz++AdwE/BR4FjgU+Clynqh+ZmCuaGMbyfsrZ71BcBZQ4TmdV5+kzueykqvYo4AEcgFsb+bc57RfjZjl/uNhjLIJNjgTq8rR/zbPJRYG2m4AMcGigrRpYh1v+T4p9PXvQbofjhN6/eHb6Xs72krYV7sfVBmDWKP2eAHqAOYG2OV7bn4p9HRNso2O9984VOe0RYDXQVYp2AvYN/P0isHY83j/Atzx7vyPQFvKOsQ2oLva1T4SdgMbg91CgfTluRbCnctrP9Ox0eU775V77ccW+9ol6P+W8L57EOUnuB/ry9Jl0drIwceF8CBCcpybI1biVT87d4yMqMqr6lKp259l0o/d8IIAXMn4ncL+qPhvYvw+4BlgKHDXBw50UiJsEdTVwF3BLnu0lbSsROQF4M86ztVlEwiJSmaffYpwdblbVTX679/fNwAoRmbmnxl0Ear3nlmCjusL77UA/lJ6dVHV1If12wS4fBlap6u2BvhngSmAa7uY+ZSjUTqq6Lfg9FGh/GSeODszZ9GHvOfc+6b+eUvfJQu2UwyU4sXzxCH0mnZ1MDBbOUTjP4BPBRlWNA8+yF9+gd4G53vMW7/lgoALnCs/lMe+5VOx3KbA/LhyVj1K3lX9TXS8itwODQL+IrBSR4Bekb4Ph7CS4sM7eyhNAF/BZEXm/iMwXkf1F5Bu4677M61fqdhqOgu0iIrNwHsPHhukbPF5J4KUpzGL7d7xfKEmzAAAS2ElEQVTPUcAmVd0QbPRet7CX28lLcfkq8B+qum6ErpPOTiYGC2c2bsm7RJ5tm3BL5UX28JgmHZ7n6//hwqA3eM2zvedNeXbx2+ZM8NCKjojsA/wH8J+qunaYbqVuq/2856txHpfzcXk0SeB6EfkHb3tJ20lVO3Ee5A5cWsE6XG7Xp4D3qurVXteSttMIjMUuZsOd+SRODF6X0z6b/HbCa9/b7fRDXJrGd0bpN+nsVL6nTziFqcTlSOQjHuhT6usjfxeXz/QFVX3Na/PDfPnsF8/pszfzI0b/oih1W9V4z73AyV7Y00+2Xg18XUSuw+wE0IcL1d0GPIITz58CbhCRd6nqPZidhmMsdjEbBhCR43DfYc8BX8/ZPNp9cq+1k4h8CDgdeLOqpkfpPunsZGKwcAZwM8nyEQ30KVlE5Ku48OdVqvqNwCbfLhV5disJ23khztOAE1Q1NULXUrfVoPf8K18IgvOEichtwN/jvIclbScROQgnAC9V1R8F2n+FE4hXe7PSS9pOIzAWu5gNPUTkCOAOXCjzLC9NKsgA+e0EzlZ7pZ282ebfBX6iqo8UsMuks5OJwcJpAZaLSEWeUPEcXAi5ZL2CInIZ8CXgZ7gQQhA/yT2f69tvG85lPuURkQrcL+k/Aq1e8jpsv/Y6r62dErcVrmQHQGuebZu95wbMTpfibho3BxtVdUBE7sD9KFuI2Wk4xmIXsyEgIocD9+BKXZ0cnHgToIXhQ5xz2Hvt9BVcqZ2rA9/vADFAvLZEIEdw0tnJcgYL50mcvY4ONopIFDgUV4OqJPGE4Fdw+SMfU2+OfIAXcC7xY/Psfoz3vDfbLwY0AWcBrwce93vbz/VefwyzlT9Ba26ebX7bVtznEYa3kwJPj+/QJhX+jSSUZ1t54LnU7TQcBdtFVTfjbs7HDNMX9u7PpC8E72V7+sZwkyOeBOaIyLyc/efh8uT2VjstwInBx9nxO/5oXMj3dVxtYp/JZ6di1vCZSg/gIEauM3huscdYJLt82bv+nwNlI/S7GVc775BAm187byV7ce08IAy8L8/jQs92d3qvl5qtaMDVedtIoHYbLlm9D3gt0Pak13d2oG2213Zvsa9lgu10hffe+WxOez3O69ABhErZToxeZ7BguwDfZvg6g51ATbGvdwLtdBiuluJ6AnX3hul7FiPXz3tzsa93IuyE+1GR7zv+JVzqy/uA0yaznWwFkjEgIlfiwi+34kJ+y3A1hR4GTtHSW4HkU7jK/etxM4hzr3+LuiR2v67XE7hVSq7AfeF+HCeyz1LVu/fUuCcLIrIQWMPOK5CUtK1E5ALgx7gv0p/iCilfiBOEb1fVP3n9jgP+ghOOV3q7XwzMAI5X1ef28ND3GF4Ji2dw4vmXuO+gabj3yULgU6r6A69vydhJRM7DeWnAXWMEd4MFWKeq1wf6FmwXEWnEeQobcSkfm3C1Z0/CRUN+MkGXNCEUaifvffY07r31H8CqPIe7VVX7A8e+HXg7bvUgf2WNfwR+oarnjf/VTBxjeT8Ns//9wJGafwWSyWWnYqvtqfTA/RL8V9wqEAncF8J3mGLV58fRHtfifsUM97g/p/8y4Pe4+mgDuCWOVhT7Oopov4XkWYHEbKUAZ+NquPXjQlN/wt2gc/sdC9yH8xr2AncDhxd7/HvIRotwqRkbcT8ceoC/AmeXqp1wqRcFfR+N1S640Pz1uNzeOE6Mn1Psa55IO+HE7kjf8QoszDl2FLcK1VrvPrka5ywIF/u6J/r9NMz+O61AMhntZJ5BwzAMwzCMEsYmkBiGYRiGYZQwJgYNwzAMwzBKGBODhmEYhmEYJYyJQcMwDMMwjBLGxKBhGIZhGEYJY2LQMAzDMAyjhDExaBiGYRiGUcKYGDSMKYyIfEREVEROKvZYio2InCEiaRHZP9B2mWefhQUeY623asCkQUSuFZGCC8J613vtBA4JEVnoneeyiTxPqSEi7xKRpIgsKfZYjNLCxKBhTAJE5CTv5uo/MiLSKSIvish1InK6iMg4n/MyEXn3eB6zWIhIOW6ZqF+q6qvFHo9RfPJ8poKPPwyzz5ki8oiI9ItIh4jcLCL7DNN3PxH5nfc57ReRB0XklGH61onIlSKySUTiIvKSiFyY+5lW1d8DLwD/vbvXbxhjobzYAzAMYwd+hVv3WoAaYD/g3cDfA/eKyPtVtSvQ/3rg10ByF871FdxyZr/brRFPDt6PW8LvQ7t5nP1wS01NZWJAZoLPsc47T3qCzzMeXAU8mNO2MbeTiJwN/AZ4DvgMUAf8M/CwiBypqi2BvouAR3DX/y2gG7cu9N0icoaq3hvoGwHuAQ7DrYH8CnAG8APcOsiX5Qzlf4HrROQAVX1pF6/ZMMaEiUHDmFw8o6q/CDaIyL/gbjj/ghOLZ/jbVDXDxN/4pwL/BDyvqs/tzkFUNTFO4ykaqhrfA+dQ3Pq8U4FHcz9TuYhIGCfUNgBvUdU+r/1O4GmcYLsgsMs3gHrgCFV91uv7c+Al4Psisr9uX+v1Y8BRwCWqeqXXdrWI/Bb4goj8TFXXBY59C/BD4JPAxbt4zYYxJixMbBiTHFXNqOq/Ag8Bp4vIm/1t+XIGRSTqhYBfE5EBEekSkRdE5Nve9oWBHLTzg+GzwDHOEZHbRGS9iCREpN0LiR2cOz4/z05E9heRO0SkV0S6ReQ3IjIzT/9aEfkvEXnFC5ltE5GHROSDOf1micgPvTEkRaRFRK4SkeacfjOBN+M8qsNRJSL/JyKtIjIoIo+LyKnDXcuuXp+ITBORK0RkVeDanhaRz4wwtoIQkSYR+bl3zH4RuU9EDs/Tb6ecQb9NRI4VkQe8/beJyDUiUp3Td56I/FRE1nn/+61e6PT8QJ+dcgY9Gw0Xll2bc44lInK9iGz2/rdrReTbIlK1u3bKh4hUiUh0hC4nArOBa3whCOAJvfuBczzBiDfGdwL3+0LQ69sHXAMsxYk/nw8DA8DVOef8LhAGzgk2esd5EHjfGC7RMHYL8wwaxtThJzjRcxZOGA7H94GPAj8HvoP7nC8B/HymNuA8XIj5QVwYLZeLgG3etlZgEc4z8rCIHK6qr+f0n4O7ad6KC7EdAnwCqAXe6ncSkXpv7AfgQnI/BEK4ENrbcSFvRGQ+8CgQ8a57FbAYuBA4WVzYrts77Ine8xMj2OTnOA/qf+PC758A7pKckN4IFHR9wM3ACcCPgOdxodRlwEnAtws4z0jcBXTgvFQzcf+jB0TkWFV9sYD9DwX+APwMuMEb0z8CWTyvl7jcy3tw1/sDYCUuXHow8BZcWsFw/BdODAVZ5I13i98gIkcAfwa6gB8Dm3D2vAQ4XkROVNWU17cC9/8qhIyqduZp/1/cNSMir+M+H/8X8NzBdvH2aJ79H8N9dpbiPH8HAxUj9PWP94SIlAGH4zz+uZ7UJ3ApCUexM48CbxPnYbQcWGPCMTFoGFOH573npaP0ew9wp6qen2+jqvYDvxCR64HVw4TQTvf6DSEuDPYscCkuLBtkMXCOqt4U6J8F/klE9lPV17zmr+OE4CdUdQcR6t04fa7EeU0OU9WNgT434264l7I912q597wq3/V6pHHhv6R3nJ8Cr3rnWTbCfgVfn4jU4UTDD1V1IsJ764D3+iJGRG4BngT+Bzi9gP0PBo5V1ce91z8WkVrgH0TkXzyP1HJc3uTnVPVbYxmcqt4TfC0iDThRsw3nHfP5KbAZOEpVewP978OFSP8OuNZr/hCekCuAdcDCwOsUcBvOY9yC8/z9I84jdyjwD4G+s73nTXmO67fNwYnBQvsCNOB+EOzUV1UTItIe6BvEfy8fgHufGsaEYmFiw5g69HjPtaP06wYOEJEDd/VEvhAUR62ITMd5FF8D3pRnl5agUPL4s/e8xDtWGfBB4JVcIeidM+v1q8N5CW8D4iIy3X8Aa4E32NEb1+Q9d4xwSVf4QtA710bgl8D+IlKIGBz1+oBBIAG8SQosZTNGvhX0Zqnq0zgv3orcUO8wPBoQgj5/xjkFFnqvfW/rybnh+LHghVR/C+wDvEdVV3ntB+FE6Q1ARc7/9iGgnx3/t3cDpxX4+LvgGFT1YVV9l6r+WFVvV9UfA8d4x/yIiBwf6F7pPefLGY3n9Bmvvn7/yjzt27znXf4fGMZYMM+gYUwdfBHYM2IvNwPyeuAFEVkN/AW4HbjdF1yjISKHAV/FhRJz87jW5NlldZ42/4bW6D1Px3lK7hrl9Pvhfqj+o/fIR/B8vkAaqfTOK3naXvae9x1m+3Dn89nh+lQ1KSL/jAtLrhGRl3Fi63eqet8oxy+E4a7hrcACnNdqJAq5hnUi8l/A54HNIvIscB9ws6o+OYaxXgWcDPy9qgZn8vrC+z+8Rz5m+H+o6macF3FcUNWsiHwDeBsu3eJhb9OA91yRZ7doTp/x6uv3H8jT7r+Xp/rMdmOKYGLQMKYO/uSN10bqpKq/9zxTZ+Ly6VbgRNWDIrIi6CHLh5ev91ec6Pyqd75+3I3pu0A+L9RIM5rHWh/R7/8Lhs9RGwz83eY9T8PNBp0ICro+Vf2RiPweJzROxE0CuEhEblTVDw53gD1EodfwJS+MfhYuT/BjwGdE5Fuq+rnRTiIiXwA+AnxNVa8f5jyXM/yPgqG8PxGJ4XIWCyGjqm2jd2Ot9zw90OaXjZnDzqLbD+NuytM3l9y+nbj36k59vXzI6cADeY4zzXsu5HoMY7cxMWgYUwffS3bHaB1VtQMnpn4hIgJ8E/gs8C7cJIeReA9O8L1TVf8S3CAijQwf8hqNdtzN8ZBR+r2BE56RAid3+JMnluBqxOVjWZ5tfq5hPo/ZLuN5s64BrhGREM5L+yERuXyM3rVclrF9goLPcpzIW7dz911HVVfj8imv9Gbh3g181ruGrcPtJyIfAL4G3Ah8OU8Xf+JRpsD/7Tnses7gcPhh/S2BNv//ciyQO65jcD+MVnqvX8B9Bo7Nc+xjvOenYMgT+QxwmIhU5JQuOhonjp/Kc5zF3nMhE4MMY7exnEHDmOSISEhE/gevfIqqPjxK3/pgm5dn9jfv5bTApr6c1z6+B2kHj56IfBw3i3WX8ELUvwKWi8hO4V9PtKKq23BJ/2eLyDH5+olIU6DJ96zs1DfApeKK//rHmIub1PCaqo4WIi4IEakUkR3yv7w6kP7En3y2Hguf9W3kne9wnNf3vmA5lN1B3EoZ4WCbNwvWt1HDCPseg/PkPg58JGe2rs/fcALnkyKyb55jlItI0E67nDPo/XDJPX4F2yce3R7Y9AAuHP2xYP6liByCS5W42Z/h7Nn6duAkb7vftxrnRX2dHWe2/wqXFxisUwgunSONE865HANsCUy8MowJxTyDhjG5OFxEzvX+Dq5AsgD4EzvOysxHDS7X6zbcjXcrLon/QpxXLngDfAw3+eBzwHqcbvw1cCcuj+l6Efmet9/xuLDzKnbve+NLuBm314jIW3GTBgRXWqYcV/IGb7wPAX/1ZjH/DffjdV+cd/PneDd1VW0TVxvwTODfhjlvOS5M/iucjT6Jm+V5yW5cSy5LcaVebsUJnk6cN+9CXJ7lUO6cuDqA5wMnq+r9BR5/AW6Fi9uAWbjSMoO4UjfjxcnAVeIKIr+G+8FwBE7kPD6KOPk9bgb4zcD7ZMeV1vpU9XeqqiJyHi6X8nkvHP0STiwtBs7G5SteC7udM3iXiLTgikb7s4nPxXkGr1TVIcGmqikR+TROmD0oIlfjcnQvxYVqv5Jz7M8DpwJ/EpErcJ7Dj+PCwWflCOGrcTOXv+Olb7yCe6++BxdKXxs8sCcq34KbdW0YewZVtYc97FHkB877oIFHBjez8yWct+X0Yfb7iNf/JO91BLc6whO4yQEJXI7UT4ElOfsuwQnMHv+8gW0n4MRYL64e3B3Agbhae2tzjrMWV4B3uGv6SE57PW5FlTdwy+htwwmlD+T0m46rzbcSN+uyCxei+19geU7fD3jnOiKn/TKv/QBc2LPVO9YTwGl5xrzTtRR6fbhJGFfgyu904YTaG7g8y1k5+/4WV99vvwLeG9d652nChZy34cT6n3Ov1+uvwLWjtQ3z/tkHVyPxFe990e/9/Z9AXWC/hd5+l+WcY7hH7ntmgXeetYH3wNO49+68cfpMfQ5X2qYNV2amCzeZ6kMj7PN23I+kAZyY/w2waJi+y3ACuMvr/xCwYpi+9cD3cKI0gZv4cxEgefqe79nswPH6frGHPUZ7iKpNVjIMY2rj5eY9BzyrqueO1r+YeCV2tgJ36DC1II3SxcsxXKuqZxd7LEbpYGLQMIy9AhE5Hc+DqeOUBzgRiMhRuBy1pRooqG0YIvJu4CbgAN15lR/DmDBMDBqGYRiGYZQwNpvYMAzDMAyjhDExaBiGYRiGUcKYGDQMwzAMwyhhTAwahmEYhmGUMCYGDcMwDMMwShgTg4ZhGIZhGCWMiUHDMAzDMIwSxsSgYRiGYRhGCfP/AejVdHi62ImuAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "\n", "gr = \"chr9:4000000-6000000\"\n", "mat1 = dhic.fetch_data(gr) # any genome region\n", "decay, std = dhic.diagonal_mean_std(mat1)\n", "\n", "# plot\n", "fig, ax = plt.subplots(figsize=(10, 5))\n", "ax.plot(decay, label=gr)\n", "ax.fill_between(np.arange(decay.shape[0]), decay - std, decay + std,\n", " color='gray', alpha=0.2)\n", "\n", "plt.yscale('log') # log scale on y-axis\n", "plt.xlabel(f\"Distance(bins, binsize={binsize})\")\n", "plt.ylabel(\"Mean contacts\")\n", "plt.xlim(0, 150)\n", "plt.ylim(1e-1, 1e5)\n", "plt.legend()" ] } ], "metadata": { "@webio": { "lastCommId": null, "lastKernelId": null }, "celltoolbar": "编辑元数据", "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.8.8" } }, "nbformat": 4, "nbformat_minor": 4 }