{ "metadata": { "name": "", "signature": "sha256:c056e589defcf1f16ca1c42ba2b8cdd1e6a9960d11f8c5e5a351f9c805ffec95" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "[``anhima.sfs``](http://anhima.readthedocs.org/en/latest/sfs.html) - Site frequency spectra" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import division, print_function\n", "import numpy as np\n", "np.random.seed(1)\n", "import scipy.stats\n", "import random\n", "random.seed(1)\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "import sys\n", "# import anhima\n", "# # dev imports\n", "sys.path.insert(0, '..')\n", "%reload_ext autoreload\n", "%autoreload 1\n", "%aimport anhima.sim\n", "%aimport anhima.gt\n", "%aimport anhima.af\n", "%aimport anhima.sf" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "n_samples = 100\n", "ploidy = 2\n", "m = n_samples * ploidy\n", "\n", "# simulated allele counts for population 1 - constant size\n", "pop1_derived_af = scipy.stats.reciprocal.rvs(a=.001, b=1, size=1000000)\n", "pop1_derived_ac = np.round(pop1_derived_af * m).astype('i4')\n", "pop1_biallelic_ac = np.column_stack([m - pop1_derived_ac, pop1_derived_ac])\n", "\n", "# simulated allele counts for population 2 - deficit of rare variants (population bottleneck)\n", "pop2_derived_af = scipy.stats.beta.rvs(a=.2, b=1.1, size=1000000)\n", "pop2_derived_ac = np.round(pop2_derived_af * m).astype('i4')\n", "pop2_biallelic_ac = np.column_stack([m - pop2_derived_ac, pop2_derived_ac])\n", "\n", "# simulated allele counts for population 3 - excess of rare variants\n", "pop3_derived_af = np.append(pop1_derived_af,\n", " scipy.stats.beta.rvs(a=.1, b=6, size=1000000))\n", "pop3_derived_ac = np.round(pop3_derived_af * m).astype('i4')\n", "pop3_biallelic_ac = np.column_stack([m - pop3_derived_ac, pop3_derived_ac])" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "pop1_sfs = anhima.sf.site_frequency_spectrum(pop1_derived_ac)\n", "pop2_sfs = anhima.sf.site_frequency_spectrum(pop2_derived_ac)\n", "pop3_sfs = anhima.sf.site_frequency_spectrum(pop3_derived_ac)\n", "\n", "fig, ax = plt.subplots()\n", "anhima.sf.plot_site_frequency_spectrum(pop1_sfs, ax=ax, label='pop1', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop2_sfs, ax=ax, label='pop2', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop3_sfs, ax=ax, label='pop3', plot_kwargs=dict(lw=2))\n", "ax.set_xlabel('derived allele count')\n", "ax.legend();" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAEPCAYAAACZcRnqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcFNW99/HPr7unh22AYRUQFBUxuIJGTUxuRs0lxCcR\nTXIjagxXuebmMS7JTYwLxuXRRH1exlxjrvokQcUFXOJ+TRRcxphExAVXVNArAsOiDMsMA7N09+/5\no6qHZpiBZpiuHma+79erX119uqrmVL2G+XKqTp1j7o6IiEgUYsWugIiIdB8KHRERiYxCR0REIqPQ\nERGRyCh0REQkMgodERGJTMFCx8xGmtnzZvaumb1jZueH5Vea2XIzWxC+vp6zzSVmttjM3jeziTnl\nh5vZ2+F3N+WUl5rZ/WH5PDPbK+e7qWa2KHx9v1DHKSIi+bNCPadjZnsAe7j7G2bWB3gNOAn4LlDr\n7je2WH8cMAv4PDACeAYY4+5uZvOBc919vpn9Gfituz9lZucAB7n7OWZ2CnCyu08xswHAK8Dh4e5f\nAw539/UFOVgREclLwVo67r7K3d8IlzcC7xGECYC1sslkYLa7N7n7EuBD4CgzGwaUufv8cL27CMIL\n4ERgZrj8EHB8uPw1YI67rw+DZi4wqcMOTkRE2iWSezpmtjcwHpgXFp1nZm+a2Qwz6x+WDQeW52y2\nnCCkWpZXsSW8RgDLANw9BWwws4Hb2ZeIiBRRwUMnvLT2J+CCsMVzKzAaOAxYCfy60HUQEZHOIVHI\nnZtZCcFlr3vc/VEAd/805/s/Ak+EH6uAkTmb70nQQqkKl1uWZ7cZBawwswTQz92rzawKqMjZZiTw\nXCv108BzIiLt4O6t3SbZoUL2XjNgBrDQ3f8zp3xYzmonA2+Hy48DU8wsaWajgTHAfHdfBdSY2VHh\nPs8AHsvZZmq4/B3g2XB5DjDRzPqbWTnwz8DTrdXT3fXqoNcVV1xR9Dp0lZfOpc5nZ37tikK2dI4B\nvge8ZWYLwrJLgVPN7DDAgY+Bfwdw94Vm9gCwEEgB5/iWozsHuBPoCfzZ3Z8Ky2cAd5vZYqAamBLu\na62ZXU3Qgw3gKlfPNRGRoitY6Lj732i9JfWX7WzzK+BXrZS/BhzcSnkDQRfs1vZ1B3BHvvUVEZHC\n04gE0mEqKiqKXYUuQ+eyY+l8dh4Fezh0d2Bm3p2PX0SkPcwMb2dHgoL2XhMR6eyC/knSlo7+j7lC\nR0S6PV3xaF0hAln3dEREJDIKnXfegWXLil0LEZFuQaFz8MFwxRXFroWISLeg0AFIpYpdAxGRbkGh\nAwodEekWVq1axYknnsiIESOIxWIsXbo08joodEChIyLdQiwW44QTTuChhx4qXh2K9pM7E4WOiHQy\ne++9N9dddx0HHnggAwYM4KyzzqKhoQGAP/zhD4wZM4aBAwcyefJkVq5c2bxdLBbj5ptvZt9992Xw\n4MH8/Oc/b+4SPmTIEH74wx9yxBFHFOWYQKETSKeLXQMR6YTMOu7VHrNmzWLOnDl89NFHLFq0iGuu\nuYbnnnuOSy+9lAcffJCVK1ey1157MWXKlK22e/TRR3nttdd4/fXXeeyxx7j99ts74Gx0DIUOqKUj\nIp2OmXHuuecyYsQIysvLmT59OrNnz2bWrFlMmzaNww47jGQyybXXXstLL7201f2Ziy66iP79+zNy\n5Eh+/OMfM3v27CIeydYUOqDQEZFWuXfcqz1Gjtwyr+WoUaNYsWIFK1asYNSoUc3lvXv3ZuDAgVRV\nVW13u85CoQMKHRHplHJbL0uXLmX48OEMHz6cTz75pLm8rq6O6upqRowY0eZ2ud8Vm0IHFDoi0um4\nO7fccgtVVVWsXbuWX/7yl0yZMoVTTz2VO+64gzfffJOGhgYuvfRSjj766K1aPzfccAPr169n2bJl\n/Pa3v+WUU05p/q6+vp76+vptlqOiAT9BHQlEpNMxM0477TQmTpzIihUrOOmkk7jsssvo0aMHV199\nNd/+9rdZt24dxxxzDPfdd99W206ePJnDDz+cDRs2cOaZZ3LWWWc1f9erV6/m/R9wwAGYGekI/wZq\nPh2Ao46CefOKXR0RKYJwbphiV2Mbo0ePZsaMGRx33HE7tV0sFuPDDz9kn3322eU6tHVudmU+HV1e\nA11eExGJiEIHFDoi0mV09knpdE8HFDoi0ul8/PHH7douyvsz7aGWDih0REQiotAB9V4TEYmIQgfU\n0hERiYhCBxQ6IiIRUeiAQkdEJCIKHVDoiIhERKED6kggIt3Ck08+yZe+9CXKy8sZNmwYZ599Nhs3\nboy0DgodUEtHRLqFmpoaLr/8clauXMl7771HVVUVF154YaR1UOiAQkdEOp1CTFd96qmnMnHiRHr0\n6EH//v05++yz+fvf/x7pcWlEAlDoiEir7KqOG1LGr9j5QUWz01X36tWLb37zm1xzzTUce+yxXHrp\npcydO5dx48bxs5/9jClTpvDCCy80b5edrrq2tpavfvWrjB07lmnTpm2z/xdeeIGDDjpol45rZyl0\nQKEjIp1O7nTVANOnT+e8885j5cqVzdNVA1x77bWUl5ezdOnS5jl1stNV9+/fv3m66pahM3fuXO66\n6y7mz58f6XEpdCCYSzaTgZiuNorIFu1pnXSktqarnjBhQnN57nTV2dDZ0XTV8+bN4/TTT+ehhx5i\nv/32K/BRbE1/ZRNh7qoHm4h0MoWYrnrBggVMnjyZO++8k2OPPbbAR7AthU42dHSJTUQ6kUJMV/3O\nO+8wadIkfve733HCCScU5bgUOgodEemEcqer3nfffRkzZgyXXXYZxx9/fPN01cOHD+fjjz9uc7rq\n8ePH841vfKP5fs6NN95IdXU1Z511FmVlZZSVlXHwwQdHe1yFmqbVzEYCdwFDAAd+7+6/NbMBwP3A\nXsAS4Lvuvj7c5hLgLCANnO/uc8Lyw4E7gR7An939grC8NPwZE4Bq4BR3/yT8biowPazONe5+Vyt1\ndO/fH9avh7Vroby840+EiHRqmq66bbvbdNVNwE/c/UDgaOBHZvY54GJgrrvvDzwbfsbMxgGnAOOA\nScAttmUKvFuBae4+BhhjZpPC8mlAdVj+G+D6cF8DgMuBI8PXFWbWv9VaxuPBu1o6IiIFV7DQcfdV\n7v5GuLwReA8YAZwIzAxXmwmcFC5PBma7e5O7LwE+BI4ys2FAmbtn+/XdlbNN7r4eAo4Pl78GzHH3\n9WErai5BkG1LHQlEpAvRdNWAme0NjAdeBoa6++rwq9XA0HB5ODAvZ7PlBCHVFC5nVYXlhO/LANw9\nZWYbzGxguK/lrexrW7qnIyKdUFedrrrgoWNmfQhaIRe4e21uCru7m1lRL6ZemR3s7oYbqDjpJCoq\nKopZHRGRTqeyspLKysoO2VdBQ8fMSggC5253fzQsXm1me7j7qvDS2adheRUwMmfzPQlaKFXhcsvy\n7DajgBVmlgD6uXu1mVUBFTnbjASea62OVw4aBBs2wLnnQsQPSYmI7A4qKiq2+g/5VVdd1e59Feye\nTtgJYAaw0N3/M+erx4Gp4fJU4NGc8ilmljSz0cAYYL67rwJqzOyocJ9nAI+1sq/vEHRMAJgDTDSz\n/mZWDvwz8HSrFVVHAhGRyBSypXMM8D3gLTNbEJZdAlwHPGBm0wi7TAO4+0IzewBYCKSAc3xLX71z\nCLpM9yToMv1UWD4DuNvMFhN0mZ4S7mutmV0NvBKud1W2W/Y2dE9HRCQyBXtOZ3dgZu6HHAJvvQVv\nvAGHHlrsKolIxDrrczqdwe72nM7uQS0dEZHIKHQUOiLSTTz//PMccsghlJeXM2DAACZOnMjChQsj\nrYNCR6EjIt3EgQceyF/+8hfWrVvH6tWrGT9+PGeddVakdVDoqPeaiHRChZiuesiQIc3THGQyGWKx\nGMOGDYv0uBQ6aumISFvMOu7VDtnpqj/66CMWLVrENddcw3PPPcell17Kgw8+yMqVK9lrr72YMmXK\nVttlp6t+/fXXeeyxx7j99tubv1u6dCnl5eX06tWLJ598khkzZuzSKdpZCh2NvSYinVDudNXl5eVM\nnz6d2bNnM2vWrObpqpPJJNdeey0vvfTSVhO3ZaerHjlyZPN01VmjRo1i3bp1rFmzhkMPPZQzzzwz\n0uNS6KilIyJtce+4Vzu0NV117oRtudNVb2+7lsrLy7nhhht44oknqKmpaVf92kOho9ARkU6qENNV\n52pqaiIWi1FaWlqA2rdOoaOOBCLSCRViuuqHH36YRYsWkclk+Oyzz/iP//gPTjjhBIVOpNTSEZFO\nqBDTVa9YsYJJkybRt29fJkyYQHl5OTNnzmztxxfuuLrz8A9m5j5lCtx3H8yaBaeeWuwqiUjEOusw\nOJquuqtSS0dEJDIKHYWOiHQhmq66s1PoiEgn1FWnq1ZLR73XREQio9BRS0dEJDIKHQ2DIyISGd3T\nUUtHpNvr7DffuxKFjkJHpFvrjM/odGW6vKaOBCIikVHoqKUjIhIZhY46EoiIREaho5aOiEhkFDoK\nHRGRyCh0FDoiIpFR6Kj3mohIZBQ6aumIiERmh6FjZgOjqEjRqPeaiEhk8mnpzDOzB83sBOuKY0Wo\npSMiEpl8Qmcs8Afg+8CHZnatme1f2GpFSKEjIhKZHYaOu2fcfY67TwHOBqYCr5jZC2b2xYLXsMCu\nv0EdCUREorLDAT/NbBBwOkFLZzVwLvAEcCjwJ2DvAtav4F5/Wy0dEZGo5DPK9D+Ae4DJ7r48p/xV\nM7utMNWKTgqFjohIVPIJnbHextjf7n5dB9cncs2ho95rIiIFl09Hgjlm1j/7wcwGmNnTBaxTpNTS\nERGJTj6hM9jd12c/uPtaYGg+Ozez281stZm9nVN2pZktN7MF4evrOd9dYmaLzex9M5uYU364mb0d\nfndTTnmpmd0fls8zs71yvptqZovC1/fbqqNCR0QkOvmETrrFH/O9gUye+78DmNSizIEb3X18+PpL\nuN9xwCnAuHCbW3KeC7oVmObuY4AxZpbd5zSgOiz/DXB9uK8BwOXAkeHritzW2lYHh3qviYhEJZ/Q\nmQ68aGb3mNk9wF+BS/PZubu/CKxr5avWHjKdDMx29yZ3XwJ8CBxlZsOAMnefH653F3BSuHwiMDNc\nfgg4Plz+GjDH3deHrbS5bBt+gFo6IiJRyuc5naeAw4H7gfuACWHZrjjPzN40sxk5LZDhQG7vuOXA\niFbKq8JywvdlYT1TwIZw2J629rUNdSQQEYlOvgN+JoG1QC0wzsz+aRd+5q3AaOAwYCXw613Y1y5T\nS0dEJDr5PBx6PcG9loVAbnPgr+35ge7+ac6+/0jwoCkELZiROavuSdBCqQqXW5ZntxkFrDCzBNDP\n3avNrAqoyNlmJPBca/VZyu1cCbB0KRWVlVRUVLS2mohIt1VZWUllZWWH7MvaeARnywpmi4CD3b2h\nXT8g6HjwhLsfHH4e5u4rw+WfAJ9399PCjgSzCG78jwCeAfZzdzezl4HzgfnAk8Bv3f0pMzsnrNv/\nNrMpwEnuPiXsSPAqMIHg/tFrBJcF15PDzPxQFvAG4+HQQ+GNN9pziCIi3YqZ4e7tGgA6n4dDPyK4\nvLbToWNms4GvAIPMbBlwBVBhZocR9GL7GPh3AHdfaGYPELSoUsA5OQ+lngPcCfQE/pxzT2kGcLeZ\nLQaqgSnhvtaa2dXAK+F6V7UMnCz1XhMRiU4+LZ2HCcZZe5YtwePufn6B61ZwZuYHsJD3GAdjx8L7\n7xe7SiIinV6hWzqPh69sOlnO8m5PvddERKKzw9Bx9zvNrBcwyt27XFNAvddERKKTz3TVJwILgKfC\nz+PN7PFCVywqCh0Rkejk85zOlcBRhCMLuPsCYJ8C1ilS6kggIhKdfEKnqZWeX/mOvdbpZVs6rtAR\nESm4fELnXTM7HUiY2Rgzu5lgYrcuwWPqSCAiEpV8Quc84ECC7tKzgRrgx4WsVJRiSd3TERGJyg6f\n0+nKzMyHltWxqrY33rMntmlTsaskItLpFfQ5HTN7vpVid/fj2vMDOxu1dEREopPPw6EX5iz3AL5N\nMExNlxBPqveaiEhU8nk49NUWRX8zs1daXXk3FC8JbmuZO2QyEMt3tgcREdlZ+VxeG5DzMQYcAfQt\nWI0iliw1mkhQQirowabQEREpmHwur73OlrHWUsASYFqhKhS1kpLgWZ0SUsEltpKSYldJRKTLyufy\n2t4R1KNokkkNhSMiEpV8Lq99m21Hlc52lXN3f7jDaxWhkhINhSMiEpV8Lq+dBXyRYLpnI5gG+iUg\nO+30bh06aumIiEQnn9BJAuNyppgeBsx09zMLWrOIKHRERKKTT1etkcCqnM+rgVGFqU70sh0JAI2/\nJiJSYPm0dJ4BnjazWQSX104B5ha0VhFSS0dEJDr5hM55wMnAl8PP/8/dHylclaK1VUtHoSMiUlD5\ndJl2M3sdqHX3uWbWy8zK3L02gvoVXDKp3msiIlHJZ7rqHwAPAreFRXsCjxayUlFSS0dEJDr5dCT4\nEfAlgnl0cPdFwJBCVipKW93TUUcCEZGCyid0Gty9IfvBzBJs+7DobksdCUREopNP6LxgZtOBXmb2\nzwSX2p4obLWio8trIiLRySd0LgI+A94G/h34M3BZISsVJXUkEBGJznZ7r4WX0t5x9wOA30dTpWip\npSMiEp3ttnTcPQV8YGZ7RVSfyOmejohIdPJ5OHQA8K6ZzQfqwjJ39xMLV63oaBgcEZHotBk6ZlYa\n9lq7jC1TGWR1zd5rTU3FrYyISBe3vZbOS8AE4Gx3/15E9YlcMgn19Ag+bN5c3MqIiHRx2wudUjM7\nHfiimX2LoLXj2ffdffK2rJISqKUs+FDbJUb2ERHptLYXOj8ETgf6Ad9s5fsuETrJpEJHRCQqbYaO\nu78IvGhmr7r7HyOsU6RKSmBNNnQ2bixuZUREurgdPhzalQMH1NIREYlSPiMSdGklJbCRPsEHhY6I\nSEEVNHTM7HYzW21mb+eUDTCzuWa2yMzmmFn/nO8uMbPFZva+mU3MKT/czN4Ov7spp7zUzO4Py+fl\nPsRqZlPDn7HIzL7fVh3V0hERiU4+8+nEzOwMM7s8/DzKzI7Mc/93AJNalF0MzHX3/YFnw8+Y2TiC\nqbDHhdvcYmbZ54NuBaa5+xhgjJll9zkNqA7LfwNcH+5rAHA5cGT4uiI33HIpdEREopNPS+cW4AvA\naeHnjWHZDoWdEda1KD4RmBkuzwROCpcnA7PdvcndlwAfAkeZ2TCgzN3nh+vdlbNN7r4eAo4Pl78G\nzHH39e6+HpjLtuEHtOgyrY4EIiIFlc8wOEe5+3gzWwDg7mvNrGQXfuZQd18dLq8GhobLw4F5Oest\nB0YATeFyVlVYTvi+LKxXysw2mNnAcF/LW9nXNpJJ3dMREYlKPi2dRjOLZz+Y2WAg0xE/3N2dIg+p\no4dDRUSik09L52bgEWCImf0K+A67Np/OajPbw91XhZfOPg3Lq4CROevtSdBCqQqXW5ZntxkFrAin\nYejn7tVmVgVU5GwzEniutcrccceVLGUdVwIVn3221UYiIgKVlZVUVlZ2yL4saGzsYCWzz7Hlfsmz\n7v5e3j/AbG/gCXc/OPz8fwlu/l9vZhcD/d394rAjwSyCG/8jgGeA/dzdzexl4HxgPvAk8Ft3f8rM\nzgEOdvf/bWZTgJPcfUrYkeBVgrHjDHgNmBDe38mtm7/5pnPcoWtYw2AoL4e1a/M9NBGRbsnMcPeW\nA0HnZYctHTO7293PAN5rpWxH284GvgIMMrNlBD3KrgMeMLNpwBLguwDuvtDMHgAWAingHN+SiOcA\ndwI9gT+7+1Nh+QzgbjNbDFQDU8J9rTWzq4FXwvWuahk4WVv1XlNHAhGRgtphS8fMFrj7+JzPCeAt\ndx9X6MoVmpn5Rx85++7rNJKkhBTU10NpabGrJiLSae1KS6fNjgRmdqmZ1QIHm1lt9kVwD+bxdta1\n0ykpATDqTJ0JREQKrc3QcfdfuXsZcIO7l+W8Brj7xRHWsaCSyeC9VqEjIlJw25s59AB3fx940Mwm\ntPze3V8vaM0iUhI+cbRR93VERApuex0JfgqcDfya1p+lObYgNYpYtqVTo2d1REQKbnvz6ZwdvldE\nVpsiyLZ0al2jEoiIFFo+A37+i5n1DZd/YWYPt3a5bXeVDZ0aV0tHRKTQ8hkG53J3rzGzLxE8IHo7\ncFthqxWdWAwSCQ2FIyIShXxCJx2+fwP4g7v/N7ArA352OhppWkQkGvmETpWZ/Z5grpsnzaxHntvt\nNjSnjohINPIJj+8CTwMTw6FkyoELC1qriGnKahGRaOxw7DV3ryOYIC37eSWwspCVippaOiIi0ehS\nl8naS3PqiIhEQ6GDRpoWEYmKQgddXhMRiYpCB3UkEBGJSj7TVXdpj77/KI2DR5BWS0dEpOC6fUvn\n5PtPZs3et+rymohIBLp96AB4slYdCUREIqDQATKJGmroG3zYsAEymeJWSESki1LoAJlkDU0kaexd\nDuk0rFlT7CqJiHRJCh3AkzUAbOw7LChY2aUGXBAR6TQUOkA6HnQe2NBToSMiUkgKHaApFrR0qpMK\nHRGRQlLoAA3UAM7qmEJHRKSQun3o9Ez0xHFI1lHlw4NChY6ISEF0+9DpWxp2lU7WsrRJLR0RkUJS\n6GRDp7SGjzYpdERECkmhkxM6i2oVOiIihdTtQ6esNBz+prSWD7Khs2IFuBevUiIiXVS3D51sS6ds\nUA119CHTuw80NMD69UWumYhI16PQCUOn94DgWZ2mQbrEJiJSKAqdZBA6PfsHoVPfX6EjIlIo3T50\nsvd0SvsGQ+HUlil0REQKpduHTvbyWkmfoKWzroceEBURKRSFThg68Z5B6KwpUUtHRKRQFDrZ53R6\nBKGzIjYi+Lx0aZFqJCLSdRUtdMxsiZm9ZWYLzGx+WDbAzOaa2SIzm2Nm/XPWv8TMFpvZ+2Y2Maf8\ncDN7O/zuppzyUjO7PyyfZ2Z7tVaP3GFwABbHDgg+L1zY4ccsItLdFbOl40CFu4939yPDsouBue6+\nP/Bs+BkzGwecAowDJgG3mJmF29wKTHP3McAYM5sUlk8DqsPy3wDXt1aJsmTQkSAVD1o676YPADNY\nvDh4XkdERDpMsS+vWYvPJwIzw+WZwEnh8mRgtrs3ufsS4EPgKDMbBpS5+/xwvbtytsnd10PA8a1V\nINvSyc6ps2J9L9hnH0ilYNGi9h+ZiIhso9gtnWfM7FUzOzssG+ruq8Pl1cDQcHk4sDxn2+XAiFbK\nq8JywvdlAO6eAjaY2YCWlciGTr2HE7lVAwcdFHz57rvtPTYREWlFoog/+xh3X2lmg4G5ZvZ+7pfu\n7mZW8AHQfv/r38NL8FnJCqCSNWsq4MAD4bHHFDoiIkBlZSWVlZUdsq+ihY67rwzfPzOzR4AjgdVm\ntoe7rwovnX0arl4FjMzZfE+CFk5VuNyyPLvNKGCFmSWAfu6+tmU9/s9V/4cbr70RSpxEZQXV1dA4\n5kCSAO+802HHKyKyu6qoqKCioqL581VXXdXufRXl8pqZ9TKzsnC5NzAReBt4HJgarjYVeDRcfhyY\nYmZJMxsNjAHmu/sqoMbMjgo7FpwBPJazTXZf3yHomLCNXiW9iFmMTU2bGL1vCoBPynR5TUSkEIrV\n0hkKPBJ2QEsA97r7HDN7FXjAzKYBS4DvArj7QjN7AFgIpIBz3JvnHjgHuBPoCfzZ3Z8Ky2cAd5vZ\nYqAamNJaRcyMsmQZGxo2MPqAWhZ/UM47TWMZE4/DRx/B5s3Qs2fHnwERkW6oKKHj7h8Dh7VSvhb4\nahvb/Ar4VSvlrwEHt1LeQBhaO9K3tC8bGjaw15haoJz3/qeUk/fbDz74AN5/H8aPz2c3IiKyA8Xu\nMt0pZHuwDds76MG2aBFwyCHBl6+9VqRaiYh0PQodtow0PXjPnND58peDLzuox4aIiCh0gC0tnf57\nbACCq2pke2pUVmrqahGRDqLQAfYsC3pdb7CP6dMH1q6F6j0OhIEDoaoq6FAgIiK7TKEDHDAoGOTz\ng+r32X//oGzRh7GtWzsiIrLLFDrkhs4HW0JnEQodEZEOptBhS+i8v2ZLS2er+zrPP6/7OiIiHUCh\nA4wuH01JrISlG5YyemwdAG++CYwbB8OHw4oV8PLLxa2kiEgXoNABErEE+w3YD4Ch44LpDP72N0hl\nYjAlHMhg1qxiVU9EpMtQ6ISyl9g2JD5g332hpgYWLABOOy1Y4YEHgjl2RESk3RQ6odz7OsceG5RV\nVgITJsCYMbB6dXBvR0RE2k2hE8oNndz+A5htae3MnNnqtiIikh+FTmjswLHA1qHz4ovhFbWpUyEW\nCy6xrVhRtDqKiOzuFDqh3JbOgCGbGTMGNm6EV18FRo+Gk0+Gpib4r/8qbkVFRHZjCp1Qvx79OGL4\nETSkG5j7P3P5+teD8nvvDVf46U+D91tvhbq6otRRRGR3p9DJcfIBJwPw8HsPc9ZZQdk998CmTcAX\nvgBHHw3r1sGNNxavkiIiuzGFTo5s6Dyx6AnGHdTE5z8P69fDn/4UrnDttcH7L38JixcXp5IiIrsx\nhU6Ozw3+HAcMOoC1m9fy10/+ytlnB+V/+EO4QkVF0KmgoQF++EPIZIpVVRGR3ZJCp4Vsa+eBdx9g\nyhQoKwtGJ2h+ROeGG4IpD557Dn7zm+JVVERkN6TQaeF7h3wPgLveuotNtpqf/zwov/DCsGEzaBDc\ncUdQePHFMH9+cSoqIrIbUui0MG7wOE4ceyL1qXpuevkmfvITGDYMXnsNZs8OV/rmN+GCC4KHeL7x\nDVi4sKh1FhHZXZh34yH7zcxbO/55y+fxhRlfoG9pXz758Sf86Z7+nH12cFXtrbeCgadpaIDJk+Hp\np2HoUHjhBRg7NvqDEBGJmJnh7taebdXSacXRex5Nxd4V1DTU8IMnfsCZZzoTJ0J1NZxxBqTTQGkp\nPPIIHH98MC7bccfBhx8Wu+oiIp2aQqcNt/2v2yhLlvHgwge5cd4NzJwJQ4YE/Qd++tNwTreePeHx\nx+ErXwmGx/nKV3SPR0RkOxQ6bRg7aCx3nXwXABc9cxHPfHoP990HySTcdFPwqI470KsX/Pd/bwme\nL38Zbr6swf9rAAAQIUlEQVRZ3alFRFqh0NmOkw44ieu/ej2O86+P/itrhz7EPfcEA0//4hfwgx8E\nt3bo0wfmzIEf/QgaG+H884PwUQcDEZGtKHR24OfH/JzpX55O2tOc8qdTaBh7D/feCz16wB//CEce\nCa+/TtAE+t3v4OGHg+5u//gHjB8Pl1wCH31U7MMQEekUFDp5uPrYq/nFP/2CtKc545EzeKbXNB57\ndhX77BP0ZjvySLjssrDVc/LJQQvn7LODVs9118F++wUdDp5+WpfdRKRbU5fpnTj+G1+6kYueuYhU\nJkWfZB8uPGo6qx+7gFtv7ol70GP6mmvgW98Kpt/hH/+A224LWj/ZkalHjoTvfS+4FDdiRGEOTESk\ngHaly7RCZyePf1H1In4252c8segJAPqV9uOfBn6XN+78V5b94wuAMWYM/Nu/BcO0DR1KMGrobbcF\n0yIsXRrsKJEI0mnqVDj22KAnnIjIbkCh007tCZ2suR/NZfpz03llxSvNZUPiY6h/+fvUvPQvUL0/\niYQxcSJ86Uvw1a/CERMy2D/+HvRue+ihLZfaSkuDaROOPRa+/nX4/OeD3goiIp2QQqeddiV0st79\n9F1mvjmTe966h5UbVzaX92zYi/p3JuH/cxysHA/r9mWPoTHGjw/uAU09bhmj/3Y3PPggvPlm2P86\nNHw4HHhgMGPp6NFwzDHBK6ZbcCJSfAqdduqI0MlKZVI88z/PcO/b9/LUh0+xZtOarX9WYx985WFQ\ndSSsPhRqRrBv33F86bBhHLHPWo5ueIH9lj5Lv+cfxaqqtv0BQ4bA/vvDqFEwblwQSp/7HAweDP36\nQTzeIcchIrIjCp126sjQyZXxDK+vfJ2nPnyKecvnsWDVAlbUrmh95dphsG401IyEDSPp1bAnE3vH\n+Fyds8/mGvZLLWHC0mfpW72k7R/Yo0cQQoceCgcfHIRQ375BQI0aFQSTWkki0kEUOu1UqNBpzad1\nn/L6yteZt3wei6oXsXT9Mt5Y9RZ1qZodb9zQi/0Wj2HP6t7ssw4OqqnnkE0bGFO3lvKGesoaN29/\n+9LSIHyGD4fevYORS4cPD3rPZV+9ewf3kUaMgP79O+agRaRLUui0wcwmAf8JxIE/uvv1Lb6PLHRa\nk/EMn6z/hKUblrKsZhnLNizjvRXLeG/FJ6xqWEJNUzVNmUY2s267++lbD4eshsNWwQFroHcjlNfD\nqA3Ba+AOMqmlTWVl1I0cgZf1Id6YwXr3IT5wACWDBhIfWE6qb29KYglKEqXYkCFBy6q0dNtXjx7B\naA1lZVtCTUR2ewqdVphZHPgA+CpQBbwCnOru7+WsU9TQyddndZ/x7mfvsnbzWtZuXseKtetYWb2J\nNes3s75hHWs3VbNmUzW1jTVsztSSitWSjtfiyY0A9G6AvTbAHhuhVxMM2gQjamBE7Zb30hQkMkFI\n9Uq1r56VQEUb36UNNibj1JUm2FwSp74kQSrZk8ZkD+qTRlOyhKZkD1I9SkmX9iTToyfesxfesxfW\nuzclJaX0cGhqrKOhqY704EEkhw6hrG8ZifoMXtdEyR79SJSVkqreDIkYifLe9Bs6lF79B0HPnnjP\nHpBMEo8liMfiJGIJ4hYnHotv9R6zGFbkgKysrKSioqKodehKdD471q6ETqKjK9OJHAl86O5LAMzs\nPmAy8N72NuqMBvceTEXvip3eLuMZNjbUsXJtDYuXr2f5mnXUNq1j/eaNfLahjjc/q2PtxjpSVkcq\nVkej11Fbv5He69YyrHoNSd9EQ7KB3pnN9G3aRP/Gevo3QP/NhsdTJEgzuA76NkBpGmavgWRZsFya\ngh4p6N0UfN+7Cfo1pOnXkM6pYR6XFjtYBthcEoRg76agrCEevBrDF0DcoTYJdSUGWDB0hwfL6Rhs\nTEJTPAPuxDyJkaAxnqEpZqTiJTTG0qTiGZxYuHWMzSUlNJQkSKaN0jSUpDO4QSZmpA3SMcPD98eq\nanh/VD8ylsBipcRiPSGWpCmWIm1pzIJ9mhmYkUrEqetRStydkrRTX1pCU0mSEhKUkCBBDCcFniYV\nS0MsjsVKiMVKsERJ0AqNxSBmYDEsln3FsZhhFgveY8HPw4JljxkWj0E8Z3sseG/en2GxGLF4nHgs\ngcUTxGJxYrF4+B+AGMRixML9xmNxzIyYhe/x4D8DZjFiFiNjGTI4aRwSMeIlpSQSpcQTSRLJUswg\n1VhPKl1POtVAwuHemffT4BuIx+PE4gmwGPFEglg8Dtm6JOLE4sFyE2nSpClN9qRnSS96lvQkEY8H\nh44FpwswHCM8ffE4mBGLWfN6QPN/YLb3OWYxErEEJbGS7B90Mp4h7WncnVh47DvzH6Lsf6iL/R+o\nlrpy6IwAluV8Xg4cVaS6FEXMYvTtUUbf4WWMHd4xox+k01BfD5s2QU1dE59uqKGmrp6Nm9NsmPFr\n3vn2ucTiKSyexi3Fxro06zak2LC2gXRtLcnNtcQb66B+A02Ny4k11FJaDyWNDSQaN5Fo3ERJajPJ\nps2UpuopbWqgR6qBWCZFfcJJUYp5KYM219K3cTOJTBP1CaMxYfSrT9EjnWFjMkbcoXdTmj6pJvo0\npemZcnqmnGRmS9hk9Uq13robBoCHr+1pCF8dZznww6raDt1nd/Yy8LV7Hmj39hkgY8H/O8wh5u0b\nQ6wpFvzHpim29b5wSDmkCMp9J96z9Qp2Yznfe045uBkxB/PgPeZgQMy9uTxtRipmpMxIx4LPTja0\nYpjHSL70ZrvPI3Tt0On81812Q/F4cHumd28YPLiEfRnY/N3bL5bzg2+PKWLttpVOB2PipdNQl4L1\n9WnSGzeTakjTVNKLVNpIbWokvbmRTH0j6U0NWMwhkSFTU0O6tpbGVIZUOkMq/J8nTSmSmxuIp4N/\n7bVNG9ic3kiPdAmkmkg11VLSZCQyMTKZDOlMinSmkXj9RmKNm2iIGalYglQ8jmfAMk4sk8YyjqXT\nWCbDwvfn8/C+E/B0I5nMZtKZTbg3kUjHibmF/4t1PJPBcUpSTfRq2Bz+0YjTo7GBklQjafPmF8TA\nY8TcMM9gnsYyKYxMUAfPEHMHd4wt7+bBCxxzmt+3fJf9w5VpXjYIvgvXjYX7i4XbxMiOQZjdd/AH\n0MJ/ts0/J1zLmv85O2DNfyjj7sErE7xn/5CmDdIW/NFMx4yNTWlWl8Sbt8vWMe7ZY8n548uW/cTD\nHxtjSzi0lCH/ACrJBK/C2d6fPd/BOvn9yfw4tWsH0JXv6RwNXOnuk8LPlwCZ3M4EZtY1D15EpMDU\nkaAFM0sQdCQ4HlgBzKdFRwIREYlWl7285u4pMzsXeJqgy/QMBY6ISHF12ZaOiIh0Pt1ybBQzm2Rm\n75vZYjO7qNj12R2Z2RIze8vMFpjZ/LBsgJnNNbNFZjbHzDS0QRvM7HYzW21mb+eUtXn+zOyS8Pf1\nfTObWJxad05tnMsrzWx5+Pu5wMy+nvOdzuV2mNlIM3vezN41s3fM7PywvEN+P7td6IQPjf4OmASM\nA041s88Vt1a7JQcq3H28ux8Zll0MzHX3/YFnw8/SujsIfgdztXr+zGwccArB7+sk4BYLHtSRQGvn\n0oEbw9/P8e7+F9C5zFMT8BN3PxA4GvhR+DeyQ34/u+PJbn5o1N2bgOxDo7LzWvZeORGYGS7PBE6K\ntjq7D3d/EbYZ36it8zcZmO3uTeHDzh8S/B4LbZ5L2Pb3E3Qud8jdV7n7G+HyRoIH6kfQQb+f3TF0\nWntoVPNG7zwHnjGzV83s7LBsqLuvDpdXA0OLU7XdVlvnbzjB72mWfmfzc56ZvWlmM3IuBelc7gQz\n2xsYT/B8bYf8fnbH0FHPiY5xjLuPB75O0Pz+cu6Xnn16Udolj/Onc7t9twKjgcOAlcCvt7OuzmUr\nzKwP8BBwgbtvNTzGrvx+dsfQqQJG5nweydYpLXlw95Xh+2fAIwTN6dVmtgeAmQ0DPi1eDXdLbZ2/\nlr+ze4Zl0gZ3/9RDwB/ZcrlH5zIPZlZCEDh3u/ujYXGH/H52x9B5FRhjZnubWZLgBtjjRa7TbsXM\neplZWbjcG5gIvE1wHqeGq00FHm19D9KGts7f48AUM0ua2WhgDMHDztKG8I9i1skEv5+gc7lDFowQ\nOgNY6O7/mfNVh/x+dtmHQ9uih0Y7xFDgkXD02gRwr7vPMbNXgQfMbBqwBPhu8arYuZnZbOArwCAz\nWwZcDlxHK+fP3Rea2QPAQoIxIc/ZLebkiEgr5/IKoMLMDiO4zPMx8O+gc5mnY4DvAW+Z2YKw7BI6\n6PdTD4eKiEhkuuPlNRERKRKFjoiIREahIyIikVHoiIhIZBQ6IiISGYWOiIhERqEjXV44zP1Pd3Kb\nb3bUtBdmtnEn119iZgPy3XZn919oZja1xcOZIs263cOh0i3t1MNoZhZ39yeAJ4rx81usn8+2ne1h\nu38F3iEY80xkK2rpSJdkZtPN7AMzexEYm1O+r5n9JRwd+69mNjYsv9PMbjOzecD/Df+3frOZ9TWz\nJTnb9zazpWYW386+RpvZSxZMcnfNdur4SLjtOzkjdW/vmC40s/nhyMlX7sI6k8zsNTN7w8yeCcsG\nmNmj4XYvmdnBYflWrcSwrqPCYaTeM7Pfh2VPm1kPM/sOcARwr5m9bmY9dnRc0r0odKTLMbPDCcbU\nOxQ4Afg8W1oDvwfOc/cjgAuBW3I2HQ58wd2b/8i6ew3whplVhEXfAJ5y9/R29nUT8F/ufgiwYjtV\nPSvc9vPA+WZWvp1jmgjsF06YNx44ouXI3q2sc3gr6wwO6/0tdz8M+E741VXAa+5+KHApcFf2FLSo\nSu7n/YDfuftBwHrg2+7+J4LxDU9z9wnuXr+d45duSJfXpCv6MvBw+Aev3sweh+bBSb8IPBiOGweQ\nDN8deLCNMaPuJwixSmAK8Ltw2Pe29vVFgkEmAe4Brm+jnheYWXYirJFsf6DEicDEnLGwehP80X9x\nJ9c5GnjB3T8BcPf1YfkxwLfCsufNbGB2UNft+Njd3wqXXwP2zvmutQnURBQ60iU5W//Ryy7HgHXh\nPECt2dRG+RPAr8KWyATgOaBsB/varrDldDxwtLvXm9nzwI4uRV3r7r/fxXVanputqtVKWYqtr4jk\n1rEhZznd4rvOdp9JOgldXpOu6K/ASeE9hjKCS2KEE1F9HN53wAKHtLGP5j/A4ZS9rwC/BZ4Ip2mp\n2c6+/k7QIgI4vY399yUIrXozO4CgBbI9TwNnha01zGxEeKlsZ9d5GfgnC2aEJNtLjqA1dHpYVgF8\nFp6vJQRBi5lNIJgYrS3Zc1YbHp/INhQ60uW4+wKCS2JvAn9m60tWpwPTzOwNgh5WJ+Zu2mI59/P9\nwGnh+472dQHBbKpvEdwnau1//U8BCTNbCFwLvNTW4YTHNBeYBbwU7vdBoM8O1nkgZx3CdT4DfgA8\nHNZ7dvjVlQT3gN4EfsWWeVMeAgaY2TvAj4APWtatlc93ArepI4G0RlMbiIhIZNTSERGRyCh0REQk\nMgodERGJjEJHREQio9AREZHIKHRERCQyCh0REYmMQkdERCLz/wEdmsxiScyxnQAAAABJRU5ErkJg\ngg==\n", "text": [ "" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "pop1_sfs_scaled = anhima.sf.site_frequency_spectrum_scaled(pop1_derived_ac)\n", "pop2_sfs_scaled = anhima.sf.site_frequency_spectrum_scaled(pop2_derived_ac)\n", "pop3_sfs_scaled = anhima.sf.site_frequency_spectrum_scaled(pop3_derived_ac)\n", "\n", "fig, ax = plt.subplots()\n", "anhima.sf.plot_site_frequency_spectrum(pop1_sfs_scaled, ax=ax, label='pop1', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop2_sfs_scaled, ax=ax, label='pop2', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop3_sfs_scaled, ax=ax, label='pop3', plot_kwargs=dict(lw=2))\n", "ax.set_ylim(bottom=0)\n", "ax.set_ylabel('scaled site frequency')\n", "ax.set_xlabel('derived allele count')\n", "ax.legend(bbox_to_anchor=(1, 1), loc='upper left');" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAEKCAYAAADO/lZMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXmcTeX/wN/nzmIZxizWwRjGkn3NUpKRJIrKt4hSklIJ\npSiJKQmtokQSqR/KElmzjS1rlrEvM5ZhZhhmYQyz3s/vj+fec+81gxlh0POe13ndc5/tfM65Z87n\nfD7P53keQ0TQaDQajUaTP1jyWwCNRqPRaP7LaEWs0Wg0Gk0+ohWxRqPRaDT5iFbEGo1Go9HkI1oR\nazQajUaTj2hFrNFoNBpNPuKe3wLkJ4Zh6LFbGo1Gcx2IiJHfMtwt/OctYnFzQ/r2RUT09i+3YcOG\n5bsMd8umr6W+nrfzprmx/OcVMVlZcORIfkuh0Wg0mv8oWhEDHD2a3xJoNBqN5j+KVsSgLGLtbvnX\ntGzZMr9FuGvQ1/LGoq+n5nbG+C/7+w3DEPHzg4QEiI2F0qXzWySNRqO57TEMA9HBWjeM/3TUNACV\nKilFfOSIVsQajUbzL9AjUa7OlV5etGu6UiX1qQO2NBqN5l+T3xHdt+t2NbQirlhRfWpFrNFoNJp8\nQCtiu0WsI6c1Go1Gkw9oRaxd0xqNRqPJR7Qi1opYo9FoNPmIVsTly4ObG0RHQ2pqfkuj0Wg0mjuA\nU6dO0aFDB8qWLYvFYiEqKuq629KK2MNDKWMR3U+s0Wg0mlxhsVho164dc+bM+fdt3QB57nyCg9Xn\n/fdDaKieZUuj0WjuIoKCghg1ahQ1a9bEz8+Pl156ibS0NAAmTZpElSpV8Pf3p2PHjsTGxpr1LBYL\n48aNIzg4mBIlSjBw4EBzKFLJkiXp3bs3jRo1+tfyaUUMMGQI1K0LiYnw0UfwySf5LZFGo9HcVRjG\njduuh+nTp7Ns2TIiIyM5dOgQn3zyCatWrWLw4MHMmjWL2NhYKlSoQJcuXVzqzZs3j23btrF9+3bm\nz5/PTz/9dAOuhit6ikvn858zB55+WlnEISEQFASdO0Pr1qofWaPRaDRXnOIy2zPVJe/GHT+vaqti\nxYq8//77vPLKKwAsWbKEN998k5YtW1K8eHFGjRoFQEpKCr6+vkRERBAYGIjFYmHp0qW0adMGgO+/\n/545c+awYsUKs+3MzEw8PT05duwYgYGBV5ThatOCaovYmU6dYOxYtR8WBlOmQNu2atKPYcPUVJga\njUajyTMiN267HsqXL2/uBwYGEhMTQ0xMjIvy9PLywt/fn+jo6KvWu9FoRXw5ffpARAQsXQrDhysl\nfOIEfPyxsowvXsxvCTUajUaTR5yjmqOioggICCAgIIDjx4+b6SkpKcTHx1O2bNkr1nPOu1FoRZwT\nwcHwyCOq7zgiAlauVGk7dsDzz8PMmbBrV35LqdFoNJpcICKMHz+e6OhoEhISGDFiBF26dOHZZ59l\nypQphIeHk5aWxuDBg2natKmLlfzFF1+QlJTEiRMnGDt2LJ07dzbzUlNTSbUNe3XezytaEV8LiwVa\ntYL586FIEZg7F559FurVU4FdUVEQF5ffUmo0Go3mChiGQdeuXWnTpg3BwcFUqVKFIUOG8NBDDzF8\n+HA6depEQEAAR48eZebMmS51O3bsSMOGDalfvz6PPfYYL730kplXuHBhvL29MQyDe+65By8vr+uT\nTwdr5eH816+Hb75RE38sWuTorDAM+OMP6Njx5giq0Wg0txHXE6yVn1SsWJHJkyfTqlWrPNWzWCxE\nRERQyT4D479AB2vdKJo3h1mzYMECWLYM6tRRaxiLwGuvQVJSfkuo0Wg0mjsMrYivl9atITxcTY3Z\nrBnExsKAAXoyEI1Go7lLMG7kmKurHed2dCPcKm6YG2XfPtVnnJEBb7wBb72lFHLlyv++bY1Go7nN\nuNNc07cD2jV9s6lRA379FTw94bvvlAKuUgVmz1bDnQYPhoUL81tKjUaj0dyGaIv4Rp7/unWqr/js\nWTh9WinoDh3ANmsLgwcrZV26NLz66o07rkaj0dxCtEWcd65mEWtFfDPOPz1dWcUnTjgfzLX/eMsW\nuPfeG39sjUajucloRZx3tGv6VuPpCYMGOb5366aGN3XoAE2aqLQJE/JHNo1Go9HcVmiL+Gad/6VL\nUL26WtFpzx615jHA4cNQtSoUKqQirn19b87xNRqN5iahLeK8oy3i/KBQIfjnH9i/36GEQQVxtW6t\nFPUvv+SffBqNRqO5LdCK+GZSvDgEBGRPf+019fnll0ohazQajeaOYtGiRTRv3hxfX1/KlClDr169\nuHDhwnW1ddMUsWEY5Q3DCDMMY69hGHsMw+hrS/czDGO5YRiHDMNYZhiGj1Od9w3DOGwYxgHDMNo4\npTc0DGO3Le8bp/QChmH8ZkvfZBhGBae8F2zHOGQYRvebdZ7XRceOULeumqf6q6/yWxqNRqPR5JHz\n588zdOhQYmNj2b9/P9HR0bz77rvX15iI3JQNKA3Us+0XAQ4C1YHPgIG29EHAKNt+DWAn4AEEARE4\n+rC3AI1t+4uBtrb914Hxtv3OwEzbvh8QCfjYtkjAJwcZJd9YtUotrVm4sEjjxiIPPSSydm32chcv\nilitt14+jUajuQK2Z2dOz/38FOuKVKhQQUaOHCk1atQQX19f6dGjh6SmpoqIyA8//CCVK1cWPz8/\n6dChg8TExJj1DMOQsWPHSqVKlaR48eLy7rvvivUKz+O5c+dK7dq1ryjDla6ZiNy6YC3DMOYB39q2\nB0XktGEYpYHVInKPYRjvA1YRGW0rvxQIBY4Dq0Skui29C9BSRHrbygwTkc2GYbgDsSJSwjCMZ4EW\nIvKarc4E23FmXiaT3Krzz5EnnlCrOjkzfLhafhHg4EFo2hQee0z3J2s0mtuG6wnWMj66cdNFyrC8\nPbeDgoLw9vZmyZIlFC5cmMcff5yQkBBCQkLo3Lkzy5cvp0aNGrzzzjuEh4ezZs0aQC36EBISwpw5\nc0hOTqZ169YMHDiQnj17ZjtG//79iYuLY/r06TnKkO/BWoZhBAH1gc1AKRE5bcs6DZSy7QcAJ52q\nnQTK5pAebUvH9nkCQEQygXOGYfhfpa3bi6lT4ccf1XrHw4aptFGjwN7P8MEHaiGJP/6ArKx8E1Oj\n0WjuZAzDoE+fPpQtWxZfX18++OADZsyYwfTp0+nZsyf16tXD09OTkSNHsnHjRqKiosy6gwYNwsfH\nh/Lly9O/f39mzJiRrf3ly5czbdo0Pv744+uSz/26zyyXGIZRBJgD9BORZOdJtO2m/82W4bbFxwfs\nb1atWimFvH69mhqzRg2YM0flpaTA3r1qtSeNRqO5A8mrFXujKe80eiUwMJCYmBhiYmJo0KCBme7l\n5YW/vz/R0dEEBgZesZ4zmzZtolu3bsyZM4fK17m+wE1VxIZheKCU8C8iMs+WfNowjNIicsowjDJA\nnC09GnAa50M5lCUbbdu/PN1eJxCIsbmmi4lIvGEY0UBLpzrlgVU5yRgaGmrut2zZkpYtW+ZU7Nbw\nwgtKEf/wg8MC9vBQi0ls2qQVsUajyRdWr17N6tWr81uMf4WzlRsVFUVAQAABAQEcP37cTE9JSSE+\nPp6yZcu6lK1evbq575y3Y8cOOnbsyNSpUwkJCbl+4a7UefxvN8AApgFfX5b+GTDItv8e2YO1PIGK\nqAArex/2ZqCJrc3Lg7W+t+13wTVY6wgqUMvXvp+DjFfsWM8XkpJEChVSQVwgEhAgMmyY2n/ppfyW\nTqPRaETkzgzWqlOnjpw8eVLi4+Pl/vvvlw8++EBWrFghJUqUkJ07d0pqaqr07dtXHnjgAbOeYRjS\nunVrSUxMlKioKLnnnntk0qRJIiKye/duKVmypPz++++5kuFK10zUE/6mKeLmgNWmXHfYtrY2JbkC\nOAQsc1aQwGBUtPQB4BGn9IbAblveWKf0AsDvwGFgExDklNfDln4YeOEKMub+l7xVdOumfpZKlUQi\nI0U2bVLfa9QQmTdP5OOPRS5dcpTPzBRJS8s/eTUazX+OO00RBwUFyahRo6RGjRri4+MjL774olyy\nPUcnTJggwcHB4ufnJ48//rhER0eb9QzDkHHjxkmlSpXE399f3nnnHTNqukePHuLm5iZFihQxt1q1\nal1RhqspYj3F5e12/idPwpQp0KuXWqUpLQ28vdVCEm5uymXdsqWKti5aFNq1g61bYccO1xm8NBqN\n5iZxp01xWbFiRSZPnkyrVq3yVM9isRAREUGlSpX+tQz5HjWtyQPlysGHHyolDFCgANiDCbKy1IIS\nq1eraTInTIClSyE+HkaPzjeRNRqNRnP9aEV8J3D//erzwQdV9HTFisoKfv11R5kff4TLovk0Go1G\nc/04j/K5qce5Hd0It4rb1Y2SjTNn4NdfVVS1nx9ERirlfPo0NGsGZcrA3LnQvz98/XV+S6vRaO5y\n7jTX9O3A1VzTWhHfqee/b59yTb/1Fpw7B/XrQ7FiSjknJaklFp3Gx2k0Gs2NQivivKMV8RW4q26a\nevUgPFwFcX34IezerfqSW7TIb8k0Gs1dhlbEeUcHa/0XeOYZ9fnWW7BrlxqJ/N576lOj0Wg0ty3a\nIr5bzv/wYahaNXv6H3+oxSU0Go3mBqEt4ryjLeL/AlWqKPc0qL7iESPU/v/+ByEhcOCA+p6cDJmZ\n+SOjRqPRaLKhFfHdxIsvqs9+/eCdd6BrVzAM1VfctSvs2QMVKqggruTk/JRUo9FoNDa0a/puOn+r\nFf75Bxo1AovtHSshQUVUR0WpGbrOn1fp3bqpNY5v0Tg5jUZz96Bd0xAWFka/fv04ceIEhmHQqFEj\nxowZQ40aNXIsr6Omr8B/5qaZOxc6dVL7FSuqIU4XL8KMGdClS/7KptFo7ji0Ioa4uDgyMjIoW7Ys\nGRkZDBkyhDVr1rBp06Ycy+s+4v86Tz6pAraKFoVZs2DUKJX+yy+OMhs2qPWRp0zJHxk1Go3mJhEU\nFMSoUaOoWbMmfn5+vPTSS6SlpQEwadIkqlSpgr+/Px07diQ2NtasZ7FYGDduHMHBwZQoUYKBAwfa\nF7egZMmS5pKIVqsVi8VCmTJlrks+bRH/V84/K0stIFG4sJoKs2xZKFJEua49PODpp2H2bDXX9dGj\n4H5Tl6rW3OUkpyVTwL0Anm6e+S2K5iZwXRbxjewGy+NzOygoCG9vb5YsWULhwoV5/PHHCQkJISQk\nhM6dO7N8+XJq1KjBO++8Q3h4OGvWrAGUIg4JCWHOnDkkJyfTunVrBg4cSM+ePQG1PnHdunU5f/48\nNWvWZPXq1fj5+eUog7aINWrlpsKF1X5AANxzD1y4AFu2qJm4FixQeSdPwrx5+Sen5rbnxLkTPDvn\nWQ7HH84xf8HBBQR8FUCrn/O20o1Gc7MwDIM+ffpQtmxZfH19+eCDD5gxYwbTp0+nZ8+e1KtXD09P\nT0aOHMnGjRuJiooy6w4aNAgfHx/Kly9P//79mTFjhpkXGBhIYmIiZ8+epW7duvTo0eO65NOK+L/K\nQw+pz5UrYc4cZS0XLKjSxo3LP7k0tz3jt45n5p6ZfLnxy2x5v4T/QseZHbmQfoG/T/zNqQun8kFC\nzW2JWln9xmzXQXmnZWIDAwOJiYkhJiaGwMBAM93Lywt/f3+io6OvWu9yfH19+eKLL1iwYAHn7QGx\neUAr4v8qzop42jS1P3q06kdeu1at8qTR5MC+s/sA2HFqR7a8zzZ8hiAUK1AMgDXH1txS2W4U6Vnp\nrD62GqtY81sUzQ3C2cqNiooiICCAgIAAjh8/bqanpKQQHx9v9v3mVM85z5mMjAwsFgsFChTIs2xa\nEf9XadlSDXFau1ZtBQuq1Z0ef1zlb9iQr+Jpbl8OnFWTw+w6vYtMq2NyGKtYTXf16/eqJTpXH1t9\ny+Wz83fU33yx4QuS0/I+Zn7kupGE/BzC9N3T81z34NmDtP21LfvO7MtzXc3NQUQYP3480dHRJCQk\nMGLECLp06cKzzz7LlClTCA8PJy0tjcGDB9O0aVMXK/mLL74gKSmJEydOMHbsWDp37gzA3LlzOXTo\nEFarlTNnzvD222/Trl07rYg1ecDXVy2hCFCypIqWLlZMjTkG2LkT4uLUJCF//529/tKl0L276l++\nnNOn9drIN4Df9/7OL+G/XLvgLSQ9K53IhEgAUjNTORR/yMw7ce4EaVlplC5Smg7VOgCw+vhqQleH\nUm9CPU5fOH1LZDyXeo42v7Sh+ZTmvLv8XTrP7uzywpAbVh1bBcCO2OxW/7X4dsu3/BX5F5O3T85z\nXc3NwTAMunbtSps2bQgODqZKlSoMGTKEhx56iOHDh9OpUycCAgI4evQoM2fOdKnbsWNHGjZsSP36\n9XnsscfMQK2YmBjatm2Lt7c3DRo0wNfXl59//vn6BBSRq25AX8D3WuXuxE2d/n+Yo0dFfvtNJCXF\nkbZ8ueqFadZM5OOP1X7RoiK7djnKWK0iVauqvBEjXNtMSxMpW1bEYhHp3l3k1Klbcip3GynpKeL+\nsbsYoYbEXYjLb3FM9pzeI4Ribr+G/2rmLYtYJoQiLaa0kPTMdPEa4eVS9pM1n9w0uRYdWiStfm4l\npy+clt/3/C6EIkU/LSq+o3yFUOSVP1+RjKyMXLWVZc2SIp8WEUKRp357Ks+ytJjSQghFnpn1TJ7r\n3inYnp13zDM1KChIVq5cmed6hmFIZGTkDZHhStdMRHJlEZcCthqG8bthGG0NQ0/FdNcQFKRWbbJH\nU4Njvupdu2D5crWfnAzt2sGff6pAifBwOGSzhCZPVjN62dm0Sa2FbLWqvufnnrslp5JfxCbHkmXN\nuuHt2t2+grA1Zmuu6/2+93caTGzAnrg9OeanZabxyoJXmLV31nXJtf/sfpfv22K3MXzNcGbvm21a\nx1X9quLh5kHzwOYuZafsnGKOwfw3TAufRt0JdVl5ZKWZ9sWGL1h1dBVz9s3hWNIxAHrW78mCZxfg\n6ebJD9t/4IEpD3A86fgVWnVwKP4QF9IvAOSqvDMiQvipcEB5CDSa3HBNRSwiHwBVgZ+AF4HDhmF8\nahhG8E2WTZMfFC+uxhKnpMC6dWrs3733qmFNHTtC69bg7H45ckStf1y5Mvz4owr+ArXYRNGisGIF\nbM29IrmTmH9gPgFfBfDN5m+y5Z1LPWc+kK+H7bHbzf0t0VtyVWdL9Bae/+N5dpzaQejq0BzLzDsw\nj0nbJzEkbMh1ybX/jFLE9xS/B4AJ/0xg6Oqh9Pyzp6mkq/hXAaB1pdYANCnbhHLe5YhMjGRd1Lps\nbYoIm05u4lLGpVzJ8OP2H9l1ehdt/68tU3ZMwSpWtsVuA+Bo0lFTEQf5BHF/4P0s7baUskXLsunk\nJt5c8iYAmdZMUjNTc2x/W8w2c//4ubwp4qhzUZxLOwfAifNaEd/p3Cq7M1d9xCJiBU4Bp4EswBeY\nbRjG5zdRNk1+YbeKQfUZr1kDX38NJUrAqlUwZozKa2UbJ/rppxAZCUOHwl9/qbTnnoPevdX+6NG3\nTvZbyIRtEwD4+0T2PvQ3l7xJvYn1WHV01XW17ayIN0dvvmb5hEsJPPXbU6RnpQPwx4E/TIXkzMLD\nCwGISIi4oiK6GgfiVaBW11pdAbiUqZTn+bTzzNqnrOyq/mo5ztfvfZ1v2n7D/C7z6V6nO6Cs4suZ\nvGMyzSY3470V713z+CLC7rjdgFKmry58lW0x2zifpoaMHEk8YirPCj4VAAipGMK6HuoFYO3xtVjF\nyoNTH6Tat9XMes7YlTrA2YtnSUlPuaZcdnad3mXuxybH5rlvWnNzOHr0KK1a5X1ce1ZWFpUqVboJ\nErlyTUVsGEY/wzC2AZ8BfwO1ROQ1oCHw1E2WT5Mf2AO2QA1zKlQI+vdXAVr2scYBAcoCtr8xFikC\nsbGwebOKxn7wQVXH01PNdX3oUPbj3MGcvXiW5ZHKdZ+TC9KunCf8M+G62nceGrQless1XbrfbfmO\n6ORompZryjM1n8EqVr7d8q1LmSxrFosPLwZUhPPBswcBuJB+gY4zO/Llhuzjgi/HbhG3rtSa4oWL\nA+BmuAEQlxIHOBRxYY/C9G3Sl1JFStGjvproYPru6WYboJTpp+s+BeC3vb9dc7hQdHI0SalJ+Bfy\np0WFFmRYMxi5fqSZfyTxiItFbKeib0UCigZwLu0cyyKXseHEBqLORTEtfJpZ5o/9f/D5359n80Dk\nxSoOP+3wgmRJVr6OoxYRPfzqDiE3FrEf8JSItBGR30UkA0wr+fGbKp0mf3C2iJ3fIhs0UMrXzQ1e\nfVUtIDF3LvzxB3zyiaNco0Zq3uqAALXKkwhMz/swkNuZufvnkiWqb/hyF+SljEscTTwKwPyD84m/\nGJ+nttOz0tl9ejcGBr4FfUm4lEBkYmS2chtPbKTb3G7sidvDt1uV0h3RagQD7xsIKBfuxYyLZvlN\nJzeRcCnB/L73jBorPi18Gn8e/JOBKwZmixJOvJTI3P1zuZRxCatYzaFL1UtUp2f9nlTzr8aXbRwK\n3MAg2Dd7r1Vlv8q8XP9l0rPSeWXhK6aCmLV3FkeT1LU6nXL6im74if9M5MftP5p937VK1uLRyo8C\nyvq3E5kYmaMiBuUiBxi1fpSZ9u2WbxERMrIy6D6vOwNXDDRfomqXrA3krZ/YWRFD/vUTW8VKwx8a\n0vTHpjekX15zc8mNIl4KJNq/GIbhbRhGEwAR0QPl7kbsFrG7OzzwgGtet26QmKj6hUEtJvHEE2oM\ncqFCKs1Zef/vf+pz/nz1edGhGP7NLDn5zcw9jiEOpy+c5lLGJdpPb0/P+T05GH8QQZ1XelZ6nsei\n7o3bS4Y1g6r+Vc2AJ7uC2ndmH9titiEivLH4Dabvnk6jHxoRlxJH3VJ1CQkKoWFAQ+4NuJdzaedM\nCxhg4SHllrbP/7wnbg8iwsRtEwH18H5j8RumkpyyYwpVv61Kp9878czsZ1h8eDGXMi9RukhpfAr6\nMKr1KA70OUD3ut2xGOpRUsGnAgXccx5H+dnDn1HKqxTro9Yzev1oLqRfYMS6EQCU91azF83YPYOn\nZz1Nz/k9zSC45LRkXlv0Gr0W9GJZ5DJAKUm7InbmfNp5UjJS8C7gjU9BH5c8uyJec9wxycjB+IOs\nPLqS7bHbzQAtgHLe5bg34F4gdxbxmmNr+GrjV2b/ckWfikD+9RPHpcSx49QOtsZs5XBCzlORam4f\ncqOIvwcuOH1PAa7P36a5M6hYUa3QNHEieHllzy9aNPsE7j4+0KePWkDCNuAdUEq5SBE1LnnaNBUM\nFhICy5Ypy7tBA8i8vfrR4lLieGPRG6ZVm1P+muNr8HTzxK+QH4Kw5vgaFh9ezE87f2LdcdUf6V3A\nG1B9oDlZJRczLrIschnTd093cSHa+4frl6lP47KNAYd7utXPrWjyYxPGbBpjuq/TstQqMm83e9sM\nLulSSy1v+dve38x27f3DL9R9AVAW8ebozew6vYvihYtTukhpNp7cyLTwaSyPXM5Lf77E2YtnsRgW\nFh5aSMeZHQF4rdFrLufhW8iXRgGNAIdbOid8C/nyffvvARi8ajD3fHsPe8/sJbBYIBMeU4+UsVvG\nMnvfbH7a+ROj/1axBceSjpkvNvaXhtqlalOnVB3KFHGsdhNQNMDcv9waBmhSronL98411X363dbv\nzIlHWlRoQSXfSvSs39NsIzcWce9FvRmwbABHk47ibnGnTXAbAE6eP3nNujeDqHOO2aA2nth4S49t\nGIbectiuRl6Ctez7WYDbv/upNLc9gwbBSy/lrc6oUWqCD2fXdsGC0Lat2u/RAy5dgtWr4ZFH1BCp\nnTvVkKfbiC82fMH4f8YzdPXQHPOXRizFKlZaVWxlKh57fzHApO2TAOjVoBelvEoRfjqcufvnurSx\n8cRGAr4M4JFfH6Hb3G5M+GcCiZcSeWHeC7y3UgUtNSjdwLTKtsVu48T5E5xOOU2WZPH2srcB6N+k\nP+2qtKN5YHNTsQA8XeNpABYdWsSF9AscSzrGnrg9eBfwpncjFUS3N26vqdh61OvBFw9/AcDA5QN5\nd/m7AHzwwAcs6bYEi2HBKlY61+zMkBbZI67bVFKKp0bxnBdFt/Nk9SeZ0nEKFsNCdHI0Vf2rsuL5\nFTxc6WHTgvXyUC9/w1YPY2v0VpegM7urvVbJWhiGQdvK6t4KLBZoXivIWRE3LNMQA/VALF64OGPa\njsHd4s6CgwvMQLNXG75KZN9IQluGmsFe17KIrWI1JzkBuK/8faZ7PifXtFWsDA0byozdM7Ll3Sic\nFfGmk7fu/0tEDL1debvSdcuNIj5qGEZfwzA8DMPwNAyjH3Dkxv10mrsGi8V1TLKdJ55Qn1YrVK+u\nptcEqKAedCxerGbjGj8erjFhelZSIgkTx0B6+pULJSWpRSyuk5VH1RCsZZHLcgx2WXR4EQDtKrcj\nsJiaCm/ZkWVmvj2qt37p+oS2DAXgvZXvkZGVgYiQkp5C93ndOZd2znRhjlo/ijcWv8G08GmcvXiW\nUl6l6FCtA3VK1QGUG3lvnOv83x4WDwbeP5BFXRexrsc6F5dw+WLlua/8fVzKvMSiQ4tYdEjJ3Ca4\nDbVK1sLd4s6RxCP8uutXDAx6NehF19pdaVGhBWcuniH8dDili5Rm8AODaRPcht/+9xsD7xtoKtHL\nGXj/QD5q+RED7x94zev7Yr0XWdJtCQPvG8jGnhup4l8FDzcPXqj7AgXdCzLr6Vn0bdyXTGsmYzaP\nyTH6u1bJWgA8VV3Fi9otWTsVilXIVqdogaLULFkTgJCgEEoXKc1jVR8jS7LMSOkHKzyYrY2cjg/q\npSDTmklsciwZ1gxKFC7B9le2M+vpWZTzLgfk7JpeH7We4WuH8/ri16/af3vi3Ak6z+7MR6s/umKZ\nK+FiEZ+8tRaxJu/kRhH3Bu4HooGTQFPglZsplOYuo1075bIG5e5etUop3h9+UGmLFyt39htvKEv5\nKsr4ULe2+PV+i50j+uZcIDoaypfH+txzrD2+NtfDR6bvns6bi98kNjnWDFiKS4ljR+wORq4byTvL\n3mH2vtlczLjIXxFqiFb7qu3Nvs2cJtCwBzRV9a9KREIEPqN98PrUi0aTGhGREEHtkrXZ98Y+apao\nyYnzJ5hUxjloAAAgAElEQVSxZwYF3QuyqecmYgbEUK14NUoXKY1fIT+SUpNYfkRZ3U3LNcXLw4ve\njXpTpuiVFyJ/psYzAPy08ycWHFLLXD5W5TE83Typ4lcFQci0ZvJKw1eo4l8FwzD49tFvzSjoD1t8\nSGEP9WL1vxr/Y/TDoynkUSjHYxUtUJShDw69qjzOtAluw+iHR+NXyLF269ePfM3Zd8/yaJVHea6O\nmghm9+nd2RRhYLFA0+3/WNXHWPH8CsY8MsZFEedkEYNSwIDZv9yzfk8zL9g3mLLejgn9r2QRRyRE\n0G1uN7xHevPivBdN+Sr6VqR+mfqU9CpJ+WLqvrC7pk+cO8HwNcM5deEU8w6oZUaTUpOIOhfF6PWj\naTypMYmXzFAcwo6GUXdCXX7f+zufrv80z5PGOLvTd8ftvq75tjW3jmuu/i4ip4HO1yqn0VwRX18V\nWZ2a6gj+KlkSWrRQFnS4U6Tppk1q4YmwMGVhO3P2LJWX/QPA+XUOVzC//aaCyP76C86cgQsXsM6b\nyyPVZvN5h3H0adyHPw/+SWW/ytQokd11ujduLy/Oe5EMawaHEg6Z/ZGg+v7+ifnH/F7Nvxrn0s5R\nzb8alXwrmYo4J6r5V8PDzYOv2nzF4zMeN92qB84ewN3iztQnplLQvSAftviQLnNUn+6HLT506cs0\nDIPaJWuz5vga0336XO3n6PViLzwsHle56Kqf+MOwD1kWuQzD9vdoFaWAapasyf6z+/Er5MeIViPM\nOrVL1WbCYxMIPxXOyw1evmr7NxrDMPDyVG5p+4QhB+MPUtFXeQ0al23MlugtppfAzkOV1Epidu8C\nXFkRDw8ZzkMVH+LxamrAR9vKbSlTpAyxF2JpGdTSpWw573K4GW7EJseSnpWOp5sniZcSeWDKA+aw\npEWHF5nucedj2u+LE+dPcDzpOC1/bsmxpGNsPLnRjDwHNUztq01fEZcSx1+Rf5l9+++tfI/EVKWY\n07PSOZZ0jMBigSSlJlHCq8Q1r2XUeYdFbBUrW2O20qqiXh/6diU344hLGobxgWEYkwzDmGLbfroV\nwmnuItq3h06dXNMKFnSNsH71VfD3V6tB7dpFNqZMwSNTuYp9DjsFwUyfDllZMHOm6nMG3DOtND0J\nSyKWsCN2Bx1nduSBKQ8QmxxrVotNjmV77HZ6/tmTDGsGgBmVW7+0ihy3K+GutbtS0qskB+PV2Nt2\nVdoBmJaPHbsyqFCsgqlU2ldtz6l3TpEwMIHot6OZ0nEKy55bRoMyDQBlbbav0p62ldvyzn3vZDtt\n+zAau3VVo0QNPN08rxkAUqpIKT5/WM25IwhNyjWhpFdJQFmEBgZftfkK/8L+LvVebvAy49qNM6Or\n84OiBYoSWCyQ9Kx0cynFUQ+NYuB9A/kk5JMc61zLNQ1QrGAxOt7T0XSvu1vcGdBsAIBLH7s9r5Jv\nJQRhwUHlURiwbACnLpyiSdkm+BT0ISk1yeyDDSoWZNYNKBqAgUFsciwtprYwreYlEUvM4VoAs/bN\nMsdfb41WM9BlWjPNiUHs98jB+IO8sfgNyn1dLld9vnbXtD3Y71YHbGnyRm5c0/MBb2A5sMhp02iu\ni52ndtJncR/lLmunFBp+fvDZZ9BBrdrDihWulaxWZOJE82uVU+kcPxsJIsgWNfOUbN6E2BQxwIPH\nVH+c3S2bcCmBnn/2RET4fuv3BI4JpOEPDdkcvZmyRctSxa+KWXd4yHDcLcph1LBMQ3558hdWPL/C\ndKXaVxdytoj9C/nzxD2qP7x6ieou4pf0KolvIV8CigbwYr0XCakYYua5WdxY2HUhS7otyVH52ftD\n7eRk1V+JXg170aJCCyVz1Q5m+kv1X+Lce+d4od4LuW7rVmM/T/uUkTVK1GD0w6OpW7pujuWdLVK7\nWzk3vN3sbZIGJfFw8MPZ8vo37Q/A+yvfZ1r4NKbsnEIBtwL8/MTP5u9ijxlwPr6Hmweli5RGEKLO\nRdGkbBP6Neln5pcorKxa5zm/7XOK22c9CywWaA65Onj2IIsOLyI9K50Pwz685jnZXdP2oD3n2cI0\ntx+5UcSFRGSQbTKP2bZtzk2X7Bbx7LOq21Jz6+i9sDffbf2On8N/Vi7lZ55RE4V4e6u5rMGx4ISd\ntWsxIiM5XgyivKFQJuxcPxuiozFOqeX1ZO9eMjY5pptsc9KT82nnGb91PKAmm1gSsQSf0T68vvh1\nMq2Z1CxRk1olazHtyWkMfmAwoGaEejj4YdpWbou7xZ1v2n6DxbBQu1Rttry8hdlPzzbdmM4WcfUS\n1elWuxvFChQzH4A3gtqlapv7foX8TKs2N1gMC3OemcP37b83lYqdogWK3jAZbwbOEdiF3Atd87wL\neRRi4H0DeePeN8xZv3KDYRgUK1gsx7xeDXpRxa8KhxMO88I89dIS2jKUasWrUbOECvy60gQiDwY9\nSEH3gnza6lPW9lhLaMtQMzLcHtRm98SAGraWZc0y5yivW6ou1fyrAWrsc0yyWlp0xZEV/B2Vw9Kk\nNlLSU4i/FI+nmycPBKquoJwmhNHcPlyzjxhYaBhGexG5K63gmTPVnBKvvprfkvw32Bu315w7eXvs\ndmjch4Sp3+NdwFvdjA+p/j7WrVN9xzNnqqFUC5RVO7MW1DgDgechduMyUlIDsI90tlgFz3MXSPEA\nrwy4NyqTAhlqxiY3w42fn/iZvkv7knApgYLuBfnhsR94vu7zpmwZWRlsjd5KnVJ18HTzZEanGcRf\njHexroL9ggn2c8wcVdKrJB4WDzKsGVQvXp2GAQ1Jei+HNZr/Bc4WcY0SNfI8EX3xwsXNIUt3Es6W\nf5BPUK7Oe/TDN3Zecw83D0a1HkWn3zvh6ebJZ60/o2+Tvtnks8vozP899X+kZqaaAW+ebp7MeWYO\nW6K30K9JP4aGDTXn6na3uJOSkcL+s/tNt3SdUnXMvvIlEUtc2v547cf89dxfOcpsj9Qu712eyn6V\nATX1p4jcskUMNHkjN4q4PzDYMIx0wP76JiLiffPEurUc1hPP3DTiL8azOXozjwQ/gpvFjZ92OMIL\ndpzawa7Tu2g8qTEVfSvy8xM/qz6tOnVUH/F996mZuNLTYaGajGJhVSjm5QeHErCG72T/qWk0ArIM\ncLPFWO0q70ETr2p47N5D05OwpqKKNO5Wpxtda3fl7MWzFHQvmM0i9HDz4Lv235nfi3gWoYhnkaue\nn8WwUM67HEeTjlK9ePWrlr1evAt4U6FYBY6fO37Ncbp3E5cr4vziqepPsaTbEoJ8gkzFCJgWsR37\nUDY7FsNiKmE7rSq2MoOmapeqzZboLRRyL0Sb4DbMPzifrdFbzWky65aqS7XiyiK2L+bRo14PpoVP\nY8WRFVzMuJitfXD0DwcWC8SvkB/FChTjXNo5zlw8kydviubWkZtlEIuIiEVECopIUdt21yhhUIr4\nDp1p8ban39J+tJ/enpcXvMzFjIv8susXM29P3B5+3fUraVlpHDh7gCY/NqHy2MosDLStCmSbDtP6\n3Xdw6BAXixZkYznwbazGepY9lkDyerW60cH7qpntFmrUDMtDysX961xodQQeCX4EUG7IEl4lbqhb\n1m4hO7uQbzR2q/jyvue7GedzzU9FDCq62lkJA+aYZFB9vvbgvNxSr5Sa+KZpuabcX/5+QM2gZreI\n65auS2CxQAq6FzTrPFzpYWqUqIFVrC4rPTlj7x+u4FMBwzDM+/NIYs7TPyReSuTsxbN5kl1zY8lN\n1LTFMIznDcMYavseaBhG45sv2q2hWDG17n1cXH5LcveRac00A1mm7pxKic9LcObiGWqVrEVV/6pk\nWjPNWagerfwohdwLEZkYydd+aqWmE96wrQxYbJNzbKtbgiw3KNygKQD1Ths0iFZR1FWHjTWPW/Ph\nbvDee0jTppRLhr9+gSdTg7ILGBYGERE5C3/0qJqAJBd81eYrxjwy5qYODxl4/0A6VutIt9rdbtox\nbjd8CvqY01bmtyLOiVJepfAt6Atcn3ztq7YHVNT8vWXVrGBLI5dy4vwJCrkXItg3GIthcQkibBTQ\niPplVET/jtgd7D+zn75L+tLmlzbmxB+mReytLHR7NLnz7F925u6fS9A3QbQbVo1d3s1Y//LUPJ+H\n5t+Tm2Ct8UAzoKvt+wVb2l1BFds9fqXnseb62Rq9laTUJIoXLk4BtwJczLhInVJ1mNxhsjksIyk1\nCQ+LB7OensW5984R3jucbv1+ZPKQ9vw1bSjLnm5gtjc3WPWMlKx/P3h6UiFJKJYG6aWK4/5Qa/VW\nBXg0vBdKlcJYt45Tz7THXaDWiEmubo+ICNUf3ayZGnvszE8/QaVKMG7c1U/w1Cno1o3aUan0a9ov\nxxmnbhQtKrRgXpd5uRpDejdh9wQ4D026XTAMw7SKr0cRd6jWgdPvnOa1Rq9xb8C9lPQqaQZ+1S5V\nGzeLmljF7p72LuBNsF8wDUqr/4ntsdt5ZeErjNsyjuVHlhO6JpSk1CRzAhK7q9w+3eblAVuTt0+m\n0++dyEg+z+RpCdRJ3oTfjO/Q3Hpy8+RoIiKvA5cARCQBuPpMAjYMw/jJMIzThmHsdkoLNQzjpGEY\nO2zbo0557xuGcdgwjAOGYbRxSm9oGMZuW943TukFDMP4zZa+yTCMCk55LxiGcci2db+SjHZFrPuJ\n/z3n086z8shKc1rIvyJVMEnnmp3Z8eoONr+8mZ2v7qRx2cbmwwTg/sD78fL0wsPNgzql6vBSg570\nHL6Qlzt+RL9RYRwtXZD4QjC1tJpEoUqp6vD662rIk78/nv3eVpN/jB8PQ4Y45rp2d6f0xF/VQhPr\n1kGbNlCrFmzcqOa7FoGzZ6GfY1gJIvCFmnOZefOyn+Q336hFK5KS1PGmT1crTCXnceaijAz49luI\n1NGsV2PkQyMZ3HwwHat1zG9RcsTeT3wlRZyVnsWa2m+wptPYHPNLepU0JzIJeyHMVJ72ceyAGTnd\noEwDLIbFtIhXHFXR055unqbbfEv0FnPREHuaXRFf7poeu0XJNG6JG7VtHsGKF/eQmXp7LcLyn0BE\nrroBm1GLPOywfS9h389F3QeA+sBup7RhwNs5lK0B7EQp+SAgAjBseVuAxrb9xUBb2/7rwHjbfmdg\npm3fD4gEfGxbJOCTwzFl6FC1Ft/gwaLJJfvi9klmVqZLWlpmmjSY2EAIRfot6Sfx8SIBQ5sIocif\nB/7M1sbyyOVCKEIoMmLtiKseb8/B9VJmoJsQiviO8s27wBMn2hdcVFvnziIvvuia9sUXIllZImFh\njrSCBUXS0kRSUlTeoUMi7u4qb9IkkebNHWX79s1+3MxMkTZtRJ56SsRqdc0bN07VK1FCZM+evJ/T\n7YLVqs5lyZL8liRf2BazTZr+2FR2ndplpp3ccFzWVH9VTm44LjvHrTXvkdVPfH3N9mKTY2X4muFy\nPOm4mfZ31N9ihBry3ZbvRETkXOo583+HUKTd/7WTfov7CcOQNxe/KYQiBT8pKKkZqSIisiJyhRCK\nNP+pudlmVFKUEIr4fFRYMgxDMg3krLuXCEjkov3XlFOpjmvrAL3lbsuNMn0O+BM11/SnwCHgmVwf\nQCnVyxXxgBzKvQ8Mcvq+FDWvdRlgv1N6F2CCU5kmtn134Ixt/1nge6c6E4AuORxTfvlFXYWnnxZN\nLlhwcIEQivzv9/+J1Um5fLDyA5eHg89zrwpDLWIM9ZDktORs7ZxNOWuW3Rq99ZrHHbxisBCKNJnU\n5KrlsrJEBg0S+dNZ92dmiowaJTJkiPqxfXxEgoPVfo8e5oNS7r9fbc4Kev58kYAAkaAgkQcecKS3\naSPi4SFiGCJubuozPNxVmE2bHOUXLHCkW60ideo48kqUEImKuuY1MLl0SWToUJGNG9X38+dFEhJy\nX/96uHRJ5PTpbMkZm/4RAbno5Z/9ZeM2YV3PKbL6yTG37Hgbyj0tArKm2ssS9tBwl/tp6/Dre2FJ\nzUh1+X+rPLay+f8zYfME2eETLEuDkcIjCguhSMjUELPs0cSjQigS8GWAmTbxn4lCKNJ9+EMiIJE+\nyIoK6n/i774zrymPVsQ3dstN1PSvwCBgJBADdBSR369V7xq8aRhGuGEYkw3DsK/eHYBaVMLOSaBs\nDunRtnRsnydscmYC5wzD8L9KW9m4013TGRlw4cK1y12LrdFbafdTF8KOrDXTYpJj6DK7C6uOrjLT\nft31K4BaM9Y2FCnsaBgj14/EwODtpmq6wKTKE8FihajmSFr2IUD+hf15oe4LPF71cRc33JX48MEP\nGfbgMHPKRmciI9U48Lg4NTvm6NHQq5dTrJWbmxqLPHw4VK2q3MqRkWqe64kTYc4cKFEC/v5bbRaL\nmpIT1EIUMTFw7Jhybxe0RbAuW6Yufv366uAiMGmSq2B/OY3zHDHC0Ue9bZsanuXvDw8+qPqo7Qtg\n5IavvoKPP4Ynn4Tjx6FuXbjnHrXEZG5IS4P/+z/o3x/i43NVPq1xc7ICgyA21iXr+A/qHAulxJN1\n4Pb7Jzq9PZpmk3vy4B/9OTRn97Ur/Etit57k3pNqycsqhxfjvX01APsLqa6YSz871oc+F3WO1Y9/\nyZHFB7K144xYhQLuBVzGANtjLAwM7s9oQr2kSB6JBLdkNdLAed7sct7lcLe4E5McQ5/FfXhzzmuk\nvf0ZdU5BswQ1RCzSDw6XVz2O6f+4RmOfP3meLUMXYs3MXfCi5jq4lqYGAm1bBdsWCATmVtOT3SIu\nCRi27RNgsi19HNDNqdyPQCegIbDcKf0BYIFtfzcQ4JQXAfgDA4APnNKHkLMVLu++O0xgmHh4DJNV\nq8LkZmK1KgPs0UdF0tOvXT4jQ2T8eJGtVzAYrVaRhx5SBt6xY7mT4eRJkbNnXdOm7Jgi7qGeQiji\nOayo6WZ7beFrpptr9dHVkpaZJt4jvc03cc/QwvLJmk/FZ5SPEIq8v+J9CQ8XofIS8XryHSne/XWh\n1E6ZPz93suWGw4dF+vRxPYeXX1YGxwcfiHz9tWl85Hzd+vd3FGjZ0pGekKBcrK1aiQwfLvLrr+Js\nycijj4pYLCJffilSu7YjfcAAkZ071b6vr7Ic7dx3n0o3DPW5cqVKf/VV9f2tt1QaKIv7ShblyZPK\nqhcRiY8XKVbMcXw/P9cTTkxUsmdkOOrv2aNuvG3b1FamjKNOaOg1r3l6vwFm+ehxc1zyIso/aOad\nGPHzNdu61axtN9KUL6zJIDP92IrDsu/Xbf+6/RPrjkpY+88lISJeRETCmg9xuW+yUL/9xnfniIDs\nK1TfrBvW7jMRkHTcZVPpjrKvUANZU6O3ZKap3zr1XKpEeN4jW/3bmMfa8NbvYs2yysh1I4VQ5L7J\n98nuHx2elwavqP/N1UdXu8jpbEH3bavKLg1Gljz5kQjI+EbIG10rioBsLtnepe7qev0kDOTFuh1k\n2LBhMmzYMG0R3+Dt2gVgj03h7QYOA5nA3lwf4DJFfKU84D3gPae8pUAToDSurmnT7Wwr09S27+ya\nNt3Xtu8Tgc45HF9E1PMTRGJj5aYyY4b5/3JN5ZSWptzlIFKunHquDhqkuiV371Zl/u//HO29+WbO\n7WRliRw9qp7xUVEiRYuqrs8BA1TX58GzB8UINdQ/6ZtVhFCk/FflZV/cPtPNRShS5NMi8vnfnwuh\nSO3xtaX0qy+6uKIrvvekbNiUaXZ9PvecesaDyOuv37hr+Pjjkk1/NGggprfYues3Rx2zbJmjwNUC\nA6KiHOXuuUddQLuSHTTIzHveb6H88ouTEDNmqDIJCUpxu7uLDByo8jp3Vj+st7cISNauPUrBli2r\n8n/5RaRbN9ebY9IkpcibNVMu6LffFgHZY6llPuTNbepUdSOAyA8/qPpWq2Q1aSoCYi1UyFTi5wsW\nFwFJDHny6hd87VqxGo7jhD85zJF3/rykGx5m3r6Q13JsImvcd5IR0lokLu7qx7JhnTNX0gYOkUvD\nRubNZX95O1lWOepZ1ZTvpFt5ycrIEmuWVaLdyomAhD08QqxZ6gVo68eLZV2vqblqOysjS9Y8/4Oc\np4hqp8VQST2XKnFGSRGQAwXrmseN8Kwu506cEwFJxVMyLqmXpL8Du7j+frZtTY3eYs2yqhcF+0vO\nuqOy01t1jWx463eJPh8tD097WFZErpANb88yy3XuhBQYXkAuZVxykfeRXx4RQhH/0f7yT2kV53C6\nkLusrt1HBOTtNkiD/kXN6+TMgYKqG2VjmSfMNK2Ib7EizlYBGtit2FyWv9wiLuO0/xYw3bZvD9by\nBCqiAqzswVqbbUrZIHuwll0pd8E1WOsIKlDL176fg2wiItK4sboSbduKLFokN5SkJJGePUXef9/V\nEHnqqavX695dXP4/hw1z7Ht5iXz6qVLQ9rRChRzPuZMnRb79VimiatUc9d9/37XNAQNEBvw1QCnT\nJ7oL7peEns1MxWvva+oyu4uL0h2waIhY3LKE6rOl6NuNpNKHjwoeKXLvvY6Xhx9+cHSRVqyo9NjZ\ns0rPDB6s4p4uJz1dZN8+kVOncr4mcXGOWKk2ykiQtDTVVWs3DuvVc5zfvffm0EhqqkjhwqrAVX7s\nM2dErBUqqHLffmumR0eLbP5CBeBkGm5SlHNSooRI6pffqrLNm6u3plnq4ZjZvIVc2rZX5QUEiPz9\ntwjIxYrVxd9f/c4ZAwa6/jBubiK//SYyfbrat6cXL27uN2KLfI2y7rejTjrplXccF+Cll5TAa22y\nYjHr/mF5UuqjHvJxxYKveh9e6tRVBGQv1UVAdldxKO74nxeIgKmMjhevn72BxES5ZAsCOvVs/6se\nS0Qkbc8hlxeMxMJlJGv3Xjm//6QkHD57zfoiIpKeLsd7fiT7HlHXJ9YoI1FuFURAdoxZLcdXHHK5\n3mvq9ZXYDUckFU9Vflu0bBw4V2IsAbL1k6XZmj88f6/sKtLMpY1tvq1k22crREAOF6ghm96f52i/\nRm8RETnurizOw/NUcN6hAjVFQNb2+EnWvjhZ1r06TS5RQNXpNlHWvTzVbCOs1cfmddnuE+IiT1jH\nr8xyQ0KQFlNaZJN54cGF0vbXtrLy/+a7yG1/YXj62QJiGYpcpKAImBZ+cmyyee8k4S3pKcqVpxVx\nPiti9RuwJ5flZqD6ldNRfbkvAdOAXUA4MA8o5VR+sM29fAB4xCm9oc0ijwDGOqUXAH63WeqbgCCn\nvB629MPAC1eQT0REPvvM5d6UPn3U8/pqWK0ia9aomJlXXlEev5x45RXXtmvXVoaSh4fyMoqIRESI\n/PWXI94mLs5RpmdP1/rOyhdE6tdXXlNQijYmRqRUKdcyYBVPT4flb1fqQcGp4j/aXynYspulQAER\nCsdJsaEVTaX754E/5WL6RTMimlBkxJStLu0XKeLYt+u4AweUsWf3nB48qK6rc73mzUVmz1bn/O23\noo6PiL+/SHL2+C4ZO9ZRt2hR1f727Zefq9o81TM1Z6U+erRI+/YiFy/m+JvNnKnqLuk5S/2AtnKJ\nicqD7EaGnHq4m/xQ8gPzeN+EJpiWrrRuLVJXPeBGFftUgitZJcvPX+X16iUCsqRib7PuCw12XfkH\nBsnq9aqZnlHUR3oxUUDEIEvKEC1PMFcE5EytlpJlUYr7Ys2G6j5t114E5CM+lB5MliF8LB6ky711\n0yQdd/Vwv3Ah55s3PV0uFvQRAXneS7lWTxaubGbva62s7+8KDZAM3CQDN5ELF8R6PllOl6kjJ+u2\nk5OvjzDPI9UoIHveGC/hXk1lw7DsCk5EZG/rviIga2gh61GBc3YFlIUh+wrVl6MelSXGUlaOLD2Y\nYxv7hs10uX7LGgySVfcNVu1We1n+fk1FaEZ6VDOVb4SXI3hu59g1sj6omwjIKUtpUymJKCv7iEcV\nEZDTllKyupO6Kc9TRMJsxwhr8JYkxyabbduDnzaV7qi+95kuqedSzeufcibFbH9N1wkiIFuKt5Ww\nRu+YMl2gsMs5OUc2hzV820z/rWYZWXJYBYRtHb5E/hm5zOXahDV9L8d/mHbvKm/YrqI1xP7CIiKy\nc+wal3Lh49eLiGhFfIO3axdQ/a327V2bcv0rvwW/ISdvU8QiIidOiIwY4bCuKlQQ+egjkd69RTp0\nUN2JDRooxffqqyING4rLvWwYIs8/L7JqlVIqvXopqxVUmy++qAJy9+wReeQRlT5mjLIa7c9vEPnw\nQ+WNBGWhJyQoaxeUooqOVsbWK6+o5/2OHSLr1jnqBwaqzwYNlEe0z6RpUuBDf+Gh9wUjU2q02C9D\nVg4Vj5dbCs8+LoQiltfqCVjlW5tRF9hgnxT/rLhU/6aOPNwmU1q2FLm//THxeK+MlBhaR7p2swoo\nF3dOSrBkSUd353PPqbSmTR1Bxp07OxQ2iDz7rKMb1X79//pLstGoketxwsNFfvwx+/ErVlTXDkQq\nV1a/4cmT2du7Ek2aOK6hHatVyZ2TzKCM1ZQVG9RbhC0x3dtPgjgiILLOv4PLCXYxZoqHhzKSQWRT\ns/4iL78se/+5KCm931ZvEnXryolXPxbvolbp1iJKjg0aLw/VOyugfn/7C1BllIXn7EJOdyugbjaQ\nFApJvbJxMm2ayMiRIkeOqPt9F7VUvU2bc74Qq1aJgOzjHhk/xklxJyeLZGbKmcLlRUB+fmW9bKO+\nCEjq8jVy6M2x2X6UCCq5fD/mWVmsGZmqf2T/fpHISLEmJkmyRblH53+0Q8IWX5TFlnYiIGfxM61F\n+xbW8O0cxf67cT8RkG0eTWSxXzc5uiVODs7ZLQJyxigua2v1FgFZ+fBIWRzyWTZZ17081cXiXVep\nu9n24XnqmsYZJZS7WUSi3IPErrQFlDUsIqvrvCkRnvdI4hH1hh3WQo2VDGsySA78puIKjnhUcZE9\nepPqEkkwfFXf8GWy2V3fq+s7vAv2KG0BJbcoBZqFIRm4ycHZKuYjKyNLuZ1B9hZ2fYB1mtxeCEUW\n1GqtZGz2vpK5nev1CWsxVEREK+IbvF27AISihhwNAz4AugEF81vwG3LyTorYzubNIjXUS+E1t1Kl\nVPraTF4AACAASURBVOzPm286XKY5bUOHuh7DuW/XbgVWqqSsYItFpKbyWMmkSaq83ap+441s4pp8\n/bXDi1mixj6ZsH6G/LT9J3H7yM3hVh5Q2sXFbG6NxkvLlsrCLFFCtbHhn2QZMCjV9Vw8UgS3NPP7\n55878gY44nnkf/9zyHXihKuF3t32TDt/Xg3dtTg8pjJ8uMg7NiPg8u7bfftUure3yJNPqv3vv1fX\nBFxjlp54QmTuXIdSB4c7/0qejh9+UN7cjRtdf7sjR9RvbFeYzm2CGtFkV9yTJonyuT//vMiYMfLt\nR2cd14fPXSqWIlaee84Rq+Xv7xpoVq+eetFo2jTn+y4xUeSbb9R5vfh8pulSdNm6KatuKt3l7ct0\nltUq8runcjvHf/6jSly4UPXTjBkjkp4u8d2Va3dMgYGSmiqy21BBakl/bZKU3/4UATlMsBw+mCUz\n/F4TAYl5pp/EeqlhMHaXZoQRLNM/2CPpuEsaHnIW9WMdf32UXCzmuDlSLMqFvdHzAUlLUyKtXpUl\n7/c8LTt3WOXC6Quy/ctVsv41FUgX7VZOsjKysv2W+7zVD7JuqMMatGZZ5ah7ZREcrvR/Pl8lKecz\nZXMB1feahHojDntwmJy2KLnSUf/Yx8MiRcShmJyV8+V9vfGHcnahb3hntgjIlhKPyvreyirfGODa\nR2XNskqMRd1sdos6CUdw3pZhC0VAEg0fuRivPDW7ijhukjNGcUmOTTbd4AKy0/sBsWZZZfuX6sXq\nhFsFWdvjJzM/xlJW3lr6lhCKfD5EvaScspSW9JR02VC2k0NOJ0WvFfEN1kX5LUC+nnwOilhEKaSZ\nM9VD/quv1EN95UoVlLp2rRqS+v336mXezuHDygKtVk25ij/5RLle27fP/vDPzBR5912HVfj006pb\n8TX1LBNQCsre53v+vHrIX8GTajJv2Rmp0P85R/CVbXt29rNSdLifEIoUG1lMXpr3krz94+/CI28J\nj/cSPFJk9WrVhj2gd9AgR//yuHEiS5eqWCK7sq9YUfXpNmyogoMzMpT1CcqF7MyGDcrA8/RUgWPO\nzJql3Mw9eijlsEB1O8r99yslPmqU8pzaA7969BDTcn/+eRXDBCoA2X7thg1TbaekKMXaqZMjLzhY\n6crkZMfvY28bXN3soPqZ7ft+fiJz5ji8DqCGJv9ke6Y1b+56bk+pZ5d06iTSBIeGP2SpKiCyZYs6\nZ7sit292T4P9/ihVSqRFC7X/zDPqujizY4eYFqmAchGDWG1vOU8xW/7+O/v9MrGiiig+8nhfkTVr\nJMuzgMuFulBY9UkPf1RVXuyvFPuBdyZJVC3VH/Jdxc9FRGT0o2EufbtHCJLxj/4pR6kgM5/6TaxW\nkd8+CJcFXx+Wea2+cTnhaMpIDKXN73Ofn5NdWCecLTu7q9ROZkqqqcBOHUh0yVvZwOHqzcQiiVHn\nRUTkn+UJ8l3jqbLg0e/EWemk4ikby6i3vjXdJoqIyHafliK4jrVd/ZTjfA4WrH1FuY+tOCxK8QVI\nWGMVFxDWMjRbOfvxBeXNWF1bvW3u9lJj6Pd4qZtyXc8pIiISYynrcj3XVlFj4w8UrGNa0Otenipr\nK7+ojtl8iGmRC8iOYi1k3OZxQijSa34vifBU8QAbB841r/PuHzdJOu6SiUWSjiVqRXyDt2sXgAWo\nCT0W5LSf3yfwr07+Cor4VhEfL7J8uWOkSXS04yEcEnL1updzKeOSNPtRBVp5fOwhj/76qNT9vq70\nW9JPrFarRJ+PlmURy8zZds6dc1jxrVo52rHF9pgjZHx9XUfCfP+9cskOGZJdhkWLlDUaH589b/v2\nK/ejO7efmKja9/BwVbJ2L8XSpY5+4QoVHMpq/34xn0Vz52Y/xvLlItXV80W6d3f1Stg3Ly/Hvt3S\ntm9Tp6oIdBHXEVDLl6sXJbscERGqjNXq8ARs3y5S0JImKag+hon0knvuccg2zxHXI6+9pl64Hn7Y\nkTZxojp2TEzO1y81VeRX4zmzwjw6mPtpeEiV0udN2Z2Z2EFZV2eLV/3/9u48Torq3P/452HYNxVB\nBEERQRA3BDQqUceohMQFVKIoMRoJGo3LL2oSNcl1csM1YhK9EqOJRgSMbMYtuCBgGDFEQRAEXNhk\nuOwiyDKsQ8/z+6NqenqGWXqke2qm5/t+veY1VadPVT9d1PB0nTp1jhc0C64GX+RyX8Zx8e030sYn\nvxI8TvOP0x9yB887ob/HMN9NIx//p03uHhzzH/NEfLtnT3nE3YPH6kq/99IPd/gWgnvP8zjNH3sw\n31d+XugTfjHPn/nem75rZ+UDg8zoHTTBLGqd7Z93vsBXP/T3YN9/n+0OvqxB9wO2mTtyVjy+Txud\ncsDrRZ2ttljQmWJFw27+zuDgM83qeLVvXbU1noyKmpvd3RePmh3fb+7Jt5Ubc6wg5jsITrJPmgS9\n7N/72YFfOhKbgz9u2tuXvrTIP29wfLzuzBue8aLEXLC7wPdTz2OYL2/YvcRJ++m4D+MdvrbSMt4S\n8PmUJb53x17fS9C8M7PrjfGR7o559Bh/+7Lg/Ys6k22jhccKYv5e+yt8VserffW/85SIU/xTeQUY\nCUwELgUuC+8R/y9wHnBe1B/goD58xIm4LPcH/T189OjktyksLPTvv/T9+KNHyzcvT2q7yy8PkvF/\n/lNcFouVvOL7wQ8O3G7DhgP/c02lxJ7PULJX9L59QeJOTJrHHBNs17178HnKe+JlyZKgfsuWxY9B\nnXJKcKX+0EPBLdEGDYJ+Vlu2FH9ROffcko/4Ft2Tb9CguJ/TddcFZUW3IVasCNbbtAm27d3bfRrB\nKEaDGVficbNYLGjOHzCguPUkP9998ODgCjjxi0p5HmkbJMk9NPRB9o/4wXmLi/zuu8veZsLDq0oc\n6EkM8guy93vLJvv8It7yR7nTBzV5LR7TKze/UaL++AbfL9HPa+xY9zsbPuHjudrffbOM3nYJHj5x\ntL/Od/yJ+1dXWK88H4+eUyKW/dTzVX95w9+9Krg/ndv5hgO22bcn5hss+HY0o+uwA15fNWNFiX3O\nafPdoEMYwT3hoseEFrQs2fSxd8fe+P3r/9z1QoVxz211UYn3WDn1wMcHEofFnNnlwM+xc9POeHP1\n7F8Ftwg21DuyxJX07CMucfegqbvoqj5I3mfE9/Npk+APbUa/B31/bL93f7y7k4P/+dU/xlsVnKDj\nWGlKxCnORZVWgHnJlNXGn5qYiAsLg/uhlY0WuHX3Vh/5/khf+dVKf3Dmg04O3ux/mvmC9QuSfq/8\n/LIHAkl8zOnFilsJ0+LOO4vfP3Hsih/9qLjOqFHu553nfvzxwS1N96DJuawm2ESnFj/e6eC+alXx\nWBnuwfEo6r0+eHBwpTt/fsl9xGJBjI88Ulw2PbiY8jZtgiv/sWOD9QEDgtfvuMO9O5/4T/mj12N/\nSgc5cXcfkR0kyfc5w2/+TnGCvY2RvqTszsU+94NC/yr8D3063/IfXL3HY7HgmPz+90Enxb/9rbj+\nu69u9i9p5fk09TFc5/f88MB7oatXBx0QK7NlS9A0/3UVxgp9VvsrfWHjPv7mYde4g++w5r6kWZBc\n/nX1k2Vu92aP4B7G27ccmDD37dxX4jGv3JNvC+/ZBk2/q7OCb6gz+j90wLYzu97oefWP8615Xx3w\nWqLV/87zvPpBx7V8msYH70i0c9PO+L3pGZf8ocz95J4a9C4vesZ3cbPTgwFLwtgX/GlmvO5Xn2+J\nx5476E/x8qLm8QUj33F39wmLJsSHwrx72Fn+7Nmn+4yLH/Y1/1l1wPsrEac4F1VaAT4FjktY70zC\nABu1+acmJuJkFT3b23h4YycHtxzzlz99OSX7DjvbeuPG5T/Zkk5FTbUdOrjn5no8aU6dWvm2lRk+\nvHh/ffpUXHffvqCpPBmxWHA7AYJe7kXPjI8YEbw+aVLx+2ZlBc+Xp9JDD8b8F/zOz6z/gU+bWujr\nONL3U8+vOfPzcrfZvdv9dhvpzzHEe3fdVum/9a5d7uefke+nn7jT77uvZB+JKO3aWehTjihumnfw\nj//+YZl1t6zf428+8J7vLyj7m+7q8Hljp3iShnc7Fz/Uvyaro29fu/2g4l03Z7XPP+S8EiN9lbao\nWTC4wbyHp5f5+rJXPy7xed9rf0W8GXph8zPjg5QU+XzKEp9x8e9974698bI92/aUuCKPFcb8xD+f\nWKJ/yedbyj5/lIhTnIsqrQD9gf8D3gl/ViU+41ubf2pDIl69bbXn7y35P+S8dfOCx45+Uy/+B/PQ\nuwd+Sz8Yzz2X+sFNkhWLBT2qPwz/L7377uSbaCtT1DwNQS/qVNqzp+QcEllZxaOgrV1bXH7WWal9\nX/fiK/JvfjP48nR6vbl+IVN9QiXj92dnB838C5JvSKmRduXH/PWBT/kWDvM19Y+ODzxRVUWdsRz8\n/fuDZovEgTXm5FTPH8VnExd47vcePyChJpp/yHnxuHJ73ul7tu3xGRcO97y3k7s1VZZ56+b50FeH\nxpupy/tyr0Sc4lyUVCVoDJwa/jSKOuiUffganog/3fSpN/ptI79gzAUlyvs918/Jwe+acpf/e9W/\nfeLiiSVmZpGKnX12cF942bLU77vo1sK8eQcOJtKpU/AX9+tfp+d9R40KBlJxDzq4XXKJxx8DKs+u\nXcEoYpli79Zdvmfz12/GKepx7BSPgLV56Ze+rFGPEs/u1gSzbi8eM7e8Juyv65637nFy8JwZOWW+\nrkSc2p/6VMLMmgF3EUz0MMzMuppZN3d/rbJt5eCMmj+KvbG9vL3ybeasncMZR53BMx8+w9QVU2nZ\nqCX3n3M/hzc9POowa53Jk2HLFujSJfX7NoMTTij7tWuvDSZOuuqq9LzvD39YvP7II8lt16RJ8JMp\nGh5ycB8m1vHYYCw+oH3fYwFo1fVwWu35mDScLgel9/DL2fR4G9r4Jhp27pDSffc8sicACzYuSOl+\npWyVToMIPEswROXZ4fo64H/SFpEAECuM8fyi5+Prj81+jNy8XG55/RYAHv32o0rCX1OrVulJwpUZ\nPhy2bYOTTqr+95bkNDg+SL4b6rWjaeumEUdTsUYtG7HsJ48x9/B+nHh3/5Tu+9QjTwXgow0fpXS/\nUrZkEvFx7j6CIBnj7jvTG5IAzMibwbod62jXvB1ZlsXExRO5cOyFFBQW8NMzf8qNp90YdYhSRWbQ\nsGHUUUhFWp9/MgBrDjsl4kiSc/afrqHPl29xyNGHpHS/3Q7vRqOsRqzcupJte7aldN9yoGQS8V4z\ni7f3mNlxwN70hVS35W3N456p93DP1HsAuKn3TVxxwhXEPIbj3HHGHTx80cMRRymSmbpddSrzfjeV\no6Y8E3UokWqQ1YATjzgRgIUbF0YcTear9B4xwVjTU4AOZjYO6AvckMaY6qwVW1Zw3ujzWLtjLQBZ\nlsV1p1xH84bN6d66O4N6DOKUtrXjm7pIbdX73ouiDqFG6Nm2Jx+u/5CPNn7EOcecE3U4Ga3CRGxm\n9Qjm870SODMsvtPdN6U7sLpm0cZFXDzuYtbuWEvfjn258bQbOemIkziu1XEA/Pf5/x1xhCJSlxTd\nJ56/fn7EkWQ+C3qiV1DBbJ67966meKqVmXllnz8dCr2QMQvG0Pmwzpx7zLlMWDyBYZOHsbNgJ307\n9uXNIW/SolGLao9LRKTI3HVzOf3p02lSvwkLb1lIl1bFPRzNDHe3CMPLKMkk4oeALwnGm4531HL3\nLekNLf2iSsRPz3uam167CYBOh3Yib2seAENOHsLTlz5NkwYZ9DyJiNRaQ14awrhF4zjn6HPIvSGX\nehZ0K1IiTq1kEnEeULqSu3vndAVVXaJIxLsKdtFlZBfW56+nYVZD9sX20aZpG4Z/azjDeg3DTOe2\niNQMm3dtpscTPfhi5xe8eNWLXHHCFYAScaqV22vazL4XLn7L3Y8t9VPrk3B1+2r3Vzy/8Hluf+N2\n1uevp1e7XuTdmcf4K8ez/I7l3NT7JiVhEalRDm96OENPGwqo93Q6VdRZ637gBeAfQK/qCSczvbb0\nNYZNHsaG/A3xshEXjqBdi3YMPmlwhJGJiFSs82HBddfKrSsjjiRzVZSIN5vZNKCzmU0u9Zq7+2Vp\njCtjvL70dS4dfykAp7c/nZ5H9qTnkT25sPOFEUcmIlK5Tod2Aoj3ZZHUqygRf5fgSvg54A9AYrtp\n9fdwqoV27N3Bj1//MQD3ffM+hn9reLyzg4hIbaBEnH7lJmJ33we8b2Z93f2LaowpI7g7P5v2M9Zs\nX0Of9n347fm/VRIWkVrn6EOOxjDWbF9DQayABlkNog4p41SaGZSEq27nvp0MeWkIf533V7Isi79d\n+jey6mVFHZaISJU1zGrIUS2PotALWbN9TdThZCRdoqWQuzN6wWi6Pd6N8YvH07xhc1686sX4CDUi\nIrVRUfO0OmylhxJxiuwu2M2Ql4bww1d/yNoda+ndrjezfzSbAd0HRB2aiMhB0X3i9Ko0EZtZNzN7\n28w+DtdPMbNfpT+02mP73u1c9NxF8avgsQPHMmfYHHq06RF1aCIiB+3YQ4N5mpWI0yOZK+KnCZ4p\n3heuLwKuSVtEtUz+vny++/x3mbV6Fh1bdmTWjbO47tTr1DFLRDKGrojTK5lpEJu6++yiUZ/c3c2s\nIL1h1R43v3Yzs1bPokPLDuTekBt/+F1EJFMoEadXMpdtm8wsPu2GmQ0C1qcvpNpj5qqZjFs0jsb1\nGzP9uulKwiKSkdRZK72SuSK+DXgK6GZm64CVwJC0RlULFMQKuP3N2wG4t++9dGvdLeKIRETSo2PL\njtSzeqzdvpZ9sX2VbyBVkswVcaG7XwAcAXR3976UHGWrzlm1dRXnjzmfhRsX0unQTvy878+jDklE\nJG0aZDWgQ8sOmBnrdqyLOpyMk8w0iPPd/bRSZfPcvXdaI6sGX2caxN0Fu+n2eDdWb19N+xbteemq\nl/hGh2+kKUIRkZph/Y71tG7amgZZDTQNYoqV2zRtZicAPYBDzOwKgqtgB1oCjasnvJpn7EdjWb19\nNSe2OZHcG3Jp3bR11CGJiKRduxbtog4hY1V0j7gbcClwSPi7yA5gWDqDqqlihTH+8N4fAPj1ub9W\nEhYRkYNW0aQPrwCvmNlZ7v5eNcZUI7k7o+aPYvmW5XQ+rDNX9rgy6pBERCQDVNQ0/Qt3HwFca2bX\nlnrZ3f2O9IZWc3y560sGTRrEO6veAeCes+6hfr1kOpyLiIhUrKJs8kn4ex7F8w8X3ZyvU/MR/+pf\nv+KdVe/Qqkkrfn72z7m5z81RhyQiIhmi0l7TJSqbZQHN3X1bkvVHARcDX7j7yWFZK2AicAyQB1zl\n7lvD1+4DbgRiwB3uPjUs7w2MJugk9oa73xmWNwLGAr2AzcDV7r4qfO164JdhKMPdfWwZ8VXaa/qz\nLz/jpCdOAmDxrYvp3rp7Mh9dRCRjqdd0aiUz6cM4M2tpZs0Ixpn+xMySfXD2WaB/qbJ7gWnufjzw\ndriOmfUAriboqd0feMKKxtWEJ4Gh7t4V6GpmRfscCmwOyx8FRoT7agX8F3BG+POAmR2aZMwl3Pf2\nfcQ8xo96/UhJWEREUi6ZAT1OdPftwEDgTaATcF0yO3f3d4GvShVfBowJl8eE+wUYAIx39wJ3zwOW\nA98ws3ZAC3efE9Ybm7BN4r5eBC4Il78NTHX3reHV9jQO/EJQqbytebzy2Ss0rt+YB857oKqbi4iI\nVCqZRFzfzBoQJL/J7l7Awd0jbuvuG8PljUDbcLk9sCah3hrgqDLK14blhL9XA7j7fmCbmR1ewb6q\nZNLHkwAY0G2AnqETEZG0SCYR/5XgXm5zYKaZdQKSukdcmfAGbY3t+DVh8QQABp80OOJIREQkU1X6\nDI67jwRGFq2b2Srg/IN4z41mdqS7bwibnb8Iy9cCHRPqdSC4kl0bLpcuL9rmaGCdmdUHDnH3zWa2\nFshO2KYj8K+ygsnJyYkvZ2dnk50dbLbkyyXM3zCflo1a0r9LlVu1RUQyRm5uLrm5uVGHkbGq1Gv6\na71BcAU9OaHX9MMEHaxGmNm9wKHufm/YWWscQeeqo4DpQJdw/uPZwB3AHOB1YKS7TzGzW4GT3f0W\nMxsMDHT3wWFnrbkEvamN4BGsXkW9sxNiK7fX9G9yf0POOzlcf+r1jB44OqXHRESkNlOv6dRK66gU\nZjYeOA9obWarCXoyPwRMMrOhhI8vAbj7J2Y2ieD55f3ArQlZ8laCx5eaEDy+NCUsfwZ4zsyWETy+\nNDjc1xYz+y3wQVjvN6WTcEV2FeziyblPAjDk5Do/46OIiKRR2q+Ia7Lyrogffe9R7pp6F73b9eaD\nYR9Q/BSViIjoiji1Khri8kqCjlRFsy6V4O4vpTGuyOwq2MWIWSMAyMnOURIWEZG0qqhp+lKCBHwE\ncDbFnZ3OB/4DZGQinrB4Aht3bqRP+z5c3PXiqMMREZEMV9HsSzcAmNk0oIe7rw/X21E8iEbGefmz\nlwEY1muYroZFRCTtknmOuCOwIWF9I8EjQxknf18+01ZMA+DS4y+tpLaIiMjBS6bX9HTgLTMbR3C/\n+GqCISMzztQVU9kb28uZHc7USFoiIlItkknEtwOXA+eE639195fTF1J0Xl3yKgADuw2spKaIiEhq\nJDOylpvZh8AOd59mZk3NrIW776iG+KpNQayA15a+BsDA7krEIiJSPZKZBvEm4AXgL2FRB+CVdAYV\nhemfT2fL7i2c0PoEurXuFnU4IiJSRyTTWesnwDeB7QDuvpTgkaaMMn7xeACuOemaiCMREZG6JJlE\nvNfd9xathJMrZNRwXLsLdscfW7rmZCViERGpPskk4nfM7JdAUzO7iKCZenJ6w6pery97nfx9+fRp\n34curbpEHY6IiNQhySTie4FNwCLgZuAN4FfpDKq6Tfp4EqBmaRERqX51ftKHglgBbX7fhq17trLi\njhV0Pqxz1GGJiNRomvQhtSqa9GFRBdu5u5+Shniq3dx1c9m6ZytdWnVREhYRkWpX2aQPGe+t5W8B\n0K9zv4gjERGRuqiiSR/yqjGOyEz9fCoA3+7y7YgjERGRuiiZAT3OMrMPzGynmRWYWaGZba+O4KrD\n7DWzqV+vPtmdsqMORURE6qBkek0/DlwLLAUaA0OBJ9IZVHWKeYyzOpxFy0Ytow5FRETqoGQSMe6+\nDMhy95i7Pwv0T29Y1adrq658+zg1S4uISDSSmX1pp5k1Aj4ys4cJ5ibOmG7rS29fyv7C/VGHISIi\ndVQyV8Q/COvdBuwimPThynQGVd3q10vm+4iIiEjqVTqgh5k1B3a7eyxczwIaufuuaogvrczM6/KA\nJiIiX4cG9EitZK6I3waaJKw3BaanJxwREZG6JZlE3Mjd84tW3H0HQTIWERGRg5RMIt5pZr2LVsys\nD7A7fSGJiIjUHcn0Uvp/wCQzWx+utwOuTl9IIiIidUdSsy+ZWUOgG+DAEncvSHdg1UGdtUREqk6d\ntVIrmSEurwIau/si4HJgopn1SntkIiIidUAy94h/7e7bzeybwAXAKOAv6Q1LRESkbkgmEcfC35cA\nT7v7a0CD9IUkIiJSdySTiNea2VMEHbReN7PGSW4nIiIilUhmZK1mBJM8LHT3ZWbWDjjZ3adWR4Dp\npM5aIiJVp85aqZVUr+lMpUQsIlJ1SsSppSZmERGRCCkRi4iIREiJWEREJEJKxCIiIhFSIhYREYlQ\nZInYzPLMbKGZzTezOWFZKzObZmZLzWyqmR2aUP8+M1tmZp+ZWb+E8t5mtih87bGE8kZmNjEsf9/M\njqneTygiIlK5KK+IHch299Pc/Yyw7F5gmrsfD7wdrmNmPQgGFOlB8EzzE2ZW1HX+SWCou3cFuppZ\n/7B8KLA5LH8UGFEdH0pERKQqom6aLv0c2mXAmHB5DDAwXB4AjHf3AnfPA5YD3wgHF2nh7nPCemMT\ntknc14sE42SLiIjUKFFfEU83s7lmNiwsa+vuG8PljUDbcLk9sCZh2zXAUWWUrw3LCX+vBnD3/cA2\nM2uV8k8hIiJyEOpH+N593X29mbUBppnZZ4kvurubWdqHvcrJyYkvZ2dnk52dne63FBGpVXJzc8nN\nzY06jIxVI4a4NLMHgHxgGMF94w1hs/MMd+9uZvcCuPtDYf0pwAPAqrDOCWH5NcC57n5LWCfH3d83\ns/rAendvU+p9NcSliEgVaYjL1IqkadrMmppZi3C5GdAPWAT8E7g+rHY98Eq4/E9gsJk1NLNjga7A\nHHffAGw3s2+EnbeuA15N2KZoX4MIOn+JiIjUKFE1TbcFXg47PtcHnnf3qWY2F5hkZkOBPOAqAHf/\nxMwmAZ8A+4FbEy5lbwVGA02AN9x9Slj+DPCcmS0DNgODq+ODiYiIVEWNaJqOipqmRUSqTk3TqRX1\n40siIiJ1mhKxiIhIhJSIRUREIqRELCIiEiElYhERkQgpEYuIiERIiVhERCRCSsQiIiIRUiIWERGJ\nkBKxiIhIhJSIRUREIqRELCIiEiElYhERkQgpEYuIiERIiVhERCRCSsQiIiIRUiIWERGJkBKxiIhI\nhJSIRUREIqRELCIiEiElYhERkQgpEYuIiERIiVhERCRCSsQiIiIRUiIWERGJkBKxiIhIhJSIRURE\nIqRELCIiEiElYhERkQgpEYuIiERIiVhERCRCSsQiIiIRUiIWERGJkBKxiIhIhJSIRUREIqRELCIi\nEiElYhERkQhldCI2s/5m9pmZLTOzX0Qdj4iISGkZm4jNLAt4HOgP9ACuMbMToo0qs+Xm5kYdQsbQ\nsUwtHU+pyTI2EQNnAMvdPc/dC4AJwICIY8po+s8udXQsU0vHU2qyTE7ERwGrE9bXhGUiIiI1RiYn\nYo86ABERkcqYe2bmKzM7E8hx9/7h+n1AobuPSKiTmR9eRCTN3N2ijiFTZHIirg8sAS4A1gFzgGvc\n/dNIAxMREUlQP+oA0sXd95vZbcBbQBbwjJKwiIjUNBl7RSwiIlIbZHJnrXJpoI+DZ2Z5ZrbQT9c4\ndgAABjlJREFUzOab2ZywrJWZTTOzpWY21cwOjTrOmsrMRpnZRjNblFBW7vEzs/vC8/UzM+sXTdQ1\nUznHMsfM1oTn53wz+07CazqWFTCzjmY2w8w+NrPFZnZHWK7zM03qXCLWQB8p40C2u5/m7meEZfcC\n09z9eODtcF3K9izBOZiozONnZj2AqwnO1/7AE2ZW5/52K1DWsXTgkfD8PM3d3wQdyyQVAD919xOB\nM4GfhP9H6vxMk7p4sDTQR+qU7jV5GTAmXB4DDKzecGoPd38X+KpUcXnHbwAw3t0L3D0PWE5wHgvl\nHks48PwEHctKufsGd18QLucDnxKMwaDzM03qYiLWQB+p4cB0M5trZsPCsrbuvjFc3gi0jSa0Wqu8\n49ee4DwtonM2Obeb2Udm9kxCM6qOZRWYWSfgNGA2Oj/Tpi4mYvVOS42+7n4a8B2CpqtzEl/0oBeg\njvXXlMTx07Gt2JPAsUBPYD3wxwrq6liWwcyaAy8Cd7r7jsTXdH6mVl1MxGuBjgnrHSn5bU6S4O7r\nw9+bgJcJmqI2mtmRAGbWDvgiughrpfKOX+lztkNYJuVw9y88BPyN4qZSHcskmFkDgiT8nLu/Ehbr\n/EyTupiI5wJdzayTmTUk6GTwz4hjqlXMrKmZtQiXmwH9gEUEx/H6sNr1wCtl70HKUd7x+ycw2Mwa\nmtmxQFeCAWqkHGGiKHI5wfkJOpaVMjMDngE+cff/TXhJ52eaZOyAHuXRQB8p0RZ4Ofh7pT7wvLtP\nNbO5wCQzGwrkAVdFF2LNZmbjgfOA1ma2Gvgv4CHKOH7u/omZTQI+AfYDt7oGAIgr41g+AGSbWU+C\nJtKVwM2gY5mkvsD3gYVmNj8suw+dn2mjAT1EREQiVBebpkVERGoMJWIREZEIKRGLiIhESIlYREQk\nQkrEIiIiEVIiFhERiZASsWS8cEq8u6u4zaWpmiLTzPKrWD/PzFolu21V959uZnZ9qQE1RKQCdW5A\nD6mTqvSwvJlluftkYHIU71+qfjLb1rTBAG4AFhOM8SwildAVsWQkM/ulmS0xs3eBbgnlx5nZm+Gs\nUTPNrFtYPtrM/mJm7wMPh1d1fzKzlmaWl7B9MzP7PzPLqmBfx5rZe2a20MyGVxDjy+G2ixNmsKro\nM/3MzOaEMwrlHESd/mY2z8wWmNn0sKyVmb0SbveemZ0clpdoTQhjPTocIvZTM3sqLHvLzBqb2SCg\nD/C8mX1oZo0r+1widZ0SsWQcM+tNMIb4qcB3gdMpvmp8Crjd3fsAPwOeSNi0PXCWu8cTj7tvBxaY\nWXZYdAkwxd1jFezrMeDP7n4KsK6CUG8Mtz0duMPMDqvgM/UDurj7GQTT0vUpPeNVGXV6l1GnTRj3\nFe7eExgUvvQbYJ67nwrcD4wtOgSlQklc7wI87u4nAVuBK939HwTjuV/r7r3cfU8Fn19EUNO0ZKZz\ngJfCJLDHzP4J8QkqzgZeCMfJBmgY/nbghXLGyJ1IkNhzgcHA4+EUceXt62yCiQYA/g6MKCfOO82s\naHL1jlQ8WH4/oF/C2L/NCBLhu1WscybwjruvAnD3rWF5X+CKsGyGmR1eNLFHBVa6+8JweR7QKeE1\nO7C6iJRFiVgykVMyERQt1wO+CudRLsuucsonAw+GV6y9gH8BLSrZV4XCK+wLgDPdfY+ZzQAqa8b9\nnbs/dZB1Sh+bEmGVUbafki1niTHuTViOlXqtpt23Fqmx1DQtmWgmMDC8Z9mCoDmZcHLzleF9TCxw\nSjn7iCcld88HPgBGApPDaW63V7CvWQRXzgBDytl/S4JEvsfMuhNcqVbkLeDG8KoeMzsqbGauap3Z\nwLlm1ims0yosf7co1vBLwqbweOURfPnAzHoBx1YQY9Ex2xF+PhFJghKxZBx3n0/QnPwR8AYlm3uH\nAEPNbAFBz97LEjcttZy4PhG4Nvxd2b7uBH5iZgsJ7juXdXU4BahvZp8AvwPeK+/jhJ9pGjAOeC/c\n7wtA80rqTEqoQ1hnE3AT8FIY9/jwpRyCe8ofAQ9SPO/si0ArM1sM/ARYUjq2MtZHA39RZy2R5Gga\nRBERkQjpilhERCRCSsQiIiIRUiIWERGJkBKxiIhIhJSIRUREIqRELCIiEiElYhERkQgpEYuIiETo\n/wP1no4QKUPSeAAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "pop1_sfs_folded = anhima.sf.site_frequency_spectrum_folded(pop1_biallelic_ac)\n", "pop2_sfs_folded = anhima.sf.site_frequency_spectrum_folded(pop2_biallelic_ac)\n", "pop3_sfs_folded = anhima.sf.site_frequency_spectrum_folded(pop3_biallelic_ac)\n", "\n", "fig, ax = plt.subplots()\n", "anhima.sf.plot_site_frequency_spectrum(pop1_sfs_folded, ax=ax, label='pop1', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop2_sfs_folded, ax=ax, label='pop2', plot_kwargs=dict(lw=2))\n", "anhima.sf.plot_site_frequency_spectrum(pop3_sfs_folded, ax=ax, label='pop3', plot_kwargs=dict(lw=2))\n", "ax.set_xlabel('minor allele count')\n", "ax.legend();" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAEPCAYAAACZcRnqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcVNWd///Xp6p6o2mgQUFAcEUMihFwzTLBJYT4VTGb\nYDTjqON3EhNjzCQTo/4UvzFR800m61eTSdwzonGMW9AoUXGLBkXjhggqgjSrrL13LZ/fH+dWd9F0\nQ9t0VTXd72ce99G3zl3q3GvoT5/PPfccc3dEREQKIVbsCoiISP+hoCMiIgWjoCMiIgWjoCMiIgWj\noCMiIgWjoCMiIgWTt6BjZmPM7Akze8PMXjezb0bls81spZm9HC2fzTnm+2a21MwWm9m0nPIpZvZa\ntO0XOeVlZnZXVP68me2Ts+1sM1sSLf+cr+sUEZGus3y9p2NmewF7ufs/zGwgsBA4DTgdqHX3/2y3\n/wTgDuBIYDTwV2Ccu7uZLQC+4e4LzOwh4Jfu/hczuwA41N0vMLOZwOfcfZaZDQVeAKZEp18ITHH3\nzXm5WBER6ZK8tXTcfY27/yNarwPeJAQTAOvgkBnAHHdPuvt7wNvA0WY2Eqhy9wXRfrcRghfAqcCt\n0fo9wAnR+meAR919cxRo5gHTe+ziRESkWwryTMfM9gUmAc9HRRea2StmdqOZDYnKRgErcw5bSQhS\n7ctraAteo4H3Adw9BWwxs2E7OJeIiBRR3oNOlFr7H+CiqMVzA7AfcDiwGvhpvusgIiK9QyKfJzez\nEkLa6w/ufh+Au6/L2f574MHoYw0wJufwvQktlJpovX159pixwCozSwCD3X2DmdUAU3OOGQM83kH9\nNPCciEg3uHtHj0l2Kp+91wy4EVjk7j/PKR+Zs9vngNei9QeAWWZWamb7AeOABe6+BthqZkdH5/wK\ncH/OMWdH618EHovWHwWmmdkQM6sGPg080lE93V2LO1deeWXR69BbFt0L3Qvdix0vuyKfLZ2PA2cB\nr5rZy1HZpcAZZnY44MAy4N8A3H2Rmf0RWASkgAu87eouAG4BKoCH3P0vUfmNwO1mthTYAMyKzrXR\nzH5A6MEGcJWr55qISNHlLei4+zN03JJ6eAfH/Aj4UQflC4GJHZQ3E7pgd3Sum4Gbu1pfERHJP41I\nIABMnTq12FXoNXQv2uhetNG96Bl5ezl0d2Bm3p+vX0SkO8wM72ZHgrz2XhMR6e1C/yTpTE//Ya6g\nIyL9njIeHctHQNYzHRERKRgFHYALL4SJE+Hpp4tdExGRPk1BB2D5cnj9ddi4sdg1ERHp0xR0AMrL\nw8+mpuLWQ0Skj1PQASgrCz+bm4tbDxGRPFqzZg2nnnoqo0ePJhaLsWLFioLXQUEH1NIRkX4hFotx\n0kkncc899xSvDkX75t5EQUdEepl9992Xa6+9lkMOOYShQ4dy7rnn0hxlY373u98xbtw4hg0bxowZ\nM1i9enXrcbFYjF/96lcccMAB7LnnnvzHf/xHa5fw4cOH89WvfpUjjjiiKNcECjqB0msi0gGznlu6\n44477uDRRx/lnXfeYcmSJVx99dU8/vjjXHrppdx9992sXr2affbZh1mzZm1z3H333cfChQt56aWX\nuP/++7npppt64G70DAUdUEtHRHodM+Mb3/gGo0ePprq6mssuu4w5c+Zwxx13cN5553H44YdTWlrK\nNddcw3PPPbfN85nvfe97DBkyhDFjxvCtb32LOXPmFPFKtqWgA21BRy0dEcnh3nNLd4wZ0zav5dix\nY1m1ahWrVq1i7NixreWVlZUMGzaMmpqaHR7XWyjoQFt6TS0dEelFclsvK1asYNSoUYwaNYrly5e3\nltfX17NhwwZGjx7d6XG524pNQQeUXhORXsfduf7666mpqWHjxo388Ic/ZNasWZxxxhncfPPNvPLK\nKzQ3N3PppZdyzDHHbNP6+clPfsLmzZt5//33+eUvf8nMmTNbtzU1NdEU/a7LXS8UDfgJSq+JSK9j\nZnz5y19m2rRprFq1itNOO43LL7+c8vJyfvCDH/CFL3yBTZs28fGPf5w777xzm2NnzJjBlClT2LJl\nC+eccw7nnntu67YBAwa0nv/ggw/GzEin04W7rv48umrrfDq33QZnnw1nnQW3317saolIAUVzwxS7\nGtvZb7/9uPHGGzn++OM/1HGxWIy3336b/ffff5fr0Nm92ZX5dJReA6XXREQKREEHlF4TkT6jt09K\np2c6oN5rItLrLFu2rFvHFfL5THeopQNKr4mIFIiCDii9JiJSIAo6oPSaiEiBKOiA0msiIgWioANK\nr4mIFIiCDii9JiJSIAo6oPSaiPQLc+fO5ROf+ATV1dWMHDmS888/n7q6uoLWQUEHNImbiPQLW7du\n5YorrmD16tW8+eab1NTU8N3vfregdVDQgW3Ta71wDCYR6X/yMV31GWecwbRp0ygvL2fIkCGcf/75\nPPvsswW9Lo1IABCPQ0kJJJNhKS0tdo1EpBewq3puSBm/8sP/QZudrnrAgAGccsopXH311Rx33HFc\neumlzJs3jwkTJvCd73yHWbNm8eSTT7Yel52uura2lhNPPJHx48dz3nnnbXf+J598kkMPPXSXruvD\nUtDJKisLAaepSUFHRIoud7pqgMsuu4wLL7yQ1atXt05XDXDNNddQXV3NihUrWufUyU5XPWTIkNbp\nqtsHnXnz5nHbbbexYMGCgl6Xgk5WeTnU1YWgM2hQsWsjIr1Ad1onPamz6aonT57cWp47XXU26Oxs\nuurnn3+eM888k3vuuYcDDzwwz1exLT3TydK7OiLSy+RjuuqXX36ZGTNmcMstt3Dcccfl+Qq2p6CT\npXd1RKQXycd01a+//jrTp0/n17/+NSeddFJRrktBJ0vv6ohIL5I7XfUBBxzAuHHjuPzyyznhhBNa\np6seNWoUy5Yt63S66kmTJnHyySe3Ps/5z//8TzZs2MC5555LVVUVVVVVTJw4sbDXla9pWs1sDHAb\nMBxw4L/c/ZdmNhS4C9gHeA843d03R8d8HzgXSAPfdPdHo/IpwC1AOfCQu18UlZdF3zEZ2ADMdPfl\n0bazgcui6lzt7rd1UEdvvf4jjoCFC+GFF8K6iPQLmq66c7vbdNVJ4GJ3PwQ4Bvi6mX0EuASY5+4H\nAY9FnzGzCcBMYAIwHbje2qbAuwE4z93HAePMbHpUfh6wISr/GXBddK6hwBXAUdFypZkN2WFtlV4T\nEcm7vAUdd1/j7v+I1uuAN4HRwKnArdFutwKnReszgDnunnT394C3gaPNbCRQ5e7Zfn235RyTe657\ngBOi9c8Aj7r75qgVNY8QyDqn9JqI9AGarhows32BScDfgRHuvjbatBYYEa2PAp7POWwlIUglo/Ws\nmqic6Of7AO6eMrMtZjYsOtfKDs7VOfVeE5FepK9OV533oGNmAwmtkIvcvTY3Cru7m1lRk6mzZ88O\nK+++y1Rgqlo6IiLbmD9/PvPnz++Rc+U16JhZCSHg3O7u90XFa81sL3dfE6XO1kXlNcCYnMP3JrRQ\naqL19uXZY8YCq8wsAQx29w1mVgNMzTlmDPB4R3VsDTpLlsDixUqviYi0M3XqVKZOndr6+aqrrur2\nufL2TCfqBHAjsMjdf56z6QHg7Gj9bOC+nPJZZlZqZvsB44AF7r4G2GpmR0fn/Apwfwfn+iKhYwLA\no8A0MxtiZtXAp4FHdlhhpddERPIuny2djwNnAa+a2ctR2feBa4E/mtl5RF2mAdx9kZn9EVgEpIAL\n2vozcwGhy3QFocv0X6LyG4HbzWwpocv0rOhcG83sB8AL0X5XZbtld0q910RE8i5vQcfdn6HzltSJ\nnRzzI+BHHZQvBLZ7g8ndm4mCVgfbbgZu7mp91dIREck/jUiQpS7TIiJ5p6CTpfSaiPRxTzzxBIcd\ndhjV1dUMHTqUadOmsWjRooLWQUEnS+k1EenjDjnkEB5++GE2bdrE2rVrmTRpEueee25B66Cgk6WW\njoj0IvmYrnr48OGt0xxkMhlisRgjR44s6HUp6GTpmY6ItGfWc0s3ZKerfuedd1iyZAlXX301jz/+\nOJdeeil33303q1evZp999mHWrFnbHJedrvqll17i/vvv56abbmrdtmLFCqqrqxkwYABz587lxhtv\n3KVb9GEp6GQpvSYivUjudNXV1dVcdtllzJkzhzvuuKN1uurS0lKuueYannvuuW0mbstOVz1mzJjW\n6aqzxo4dy6ZNm/jggw/46Ec/yjnnnFPQ61LQyVJ6TUTac++5pRs6m646d8K23Omqd3Rce9XV1fzk\nJz/hwQcfZOvWrd2qX3co6GQpvSYivUw+pqvOlUwmicVilGX/6C4ABZ0spddEpBfJx3TVf/rTn1iy\nZAmZTIb169fz7W9/m5NOOklBpyiUXhORXiQf01WvWrWK6dOnM2jQICZPnkx1dTW33nprR1+fv+vq\njdO0Fso201U/+yx84hNw7LHwt78Vt2IiUjCarrpzu9t01bsXpddERPJOQSdL6TUR6QM0XfXuQr3X\nRKQX6avTVaulk6X0mohI3inoZCm9JiKSdwo6WWrpiIjknZ7pZOmZjki/1dsfvvclCjpZiUQYCTaV\ngnQa4vFi10hECqA3vqPTlym9lmWmFJuISJ4p6ORSZwIRkbxS0Mml5zoiInmloJNL6TURkbxS0Mml\n9JqISF4p6ORSek1EJK8UdHIpvSYiklcKOrmUXhMRyaudBh0zG1aIivQKSq+JiORVV1o6z5vZ3WZ2\nkvX1sSKUXhMRyauuBJ3xwO+AfwbeNrNrzOyg/FarsObPh1tugYaM0msiIvm006Dj7hl3f9TdZwHn\nA2cDL5jZk2b2sbzXsACuuw7OOQc21iu9JiKSTzsd8NPM9gDOJLR01gLfAB4EPgr8D7BvHutXEJWV\n4WdzTOk1EZF86soo038D/gDMcPeVOeUvmtlv8lOtwho4MPxsRuk1EZF86krQGe+djP3t7tf2cH2K\nItvSaUItHRGRfOpKR4JHzWxI9oOZDTWzR/JYp4LLBp3GjJ7piIjkU1eCzp7uvjn7wd03AiO6cnIz\nu8nM1prZazlls81spZm9HC2fzdn2fTNbamaLzWxaTvkUM3st2vaLnPIyM7srKn/ezPbJ2Xa2mS2J\nln/eUT3bgo7SayIi+dSVoJNu98t8XyDTxfPfDExvV+bAf7r7pGh5ODrvBGAmMCE65vqc94JuAM5z\n93HAODPLnvM8YENU/jPguuhcQ4ErgKOi5crc1lp72aDTkFF6TUQkn7oSdC4DnjazP5jZH4CngEu7\ncnJ3fxrY1MGmjl4ynQHMcfeku78HvA0cbWYjgSp3XxDtdxtwWrR+KnBrtH4PcEK0/hngUXffHLXS\n5rF98GuVDTr1aaXXRETyqSvv6fwFmALcBdwJTI7KdsWFZvaKmd2Y0wIZBeT2jlsJjO6gvCYqJ/r5\nflTPFLAlGrans3N1KNt7rT6t9JqISD51dcDPUmAjUAtMMLN/2oXvvAHYDzgcWA38dBfO1SOyLZ26\npNJrIiL51JWXQ68jPGtZBKRzNj3VnS9093U55/494UVTCC2YMTm77k1oodRE6+3Ls8eMBVaZWQIY\n7O4bzKwGmJpzzBjg8Y7qM3v2bN5+O6y/uinqGa6WjohIq/nz5zN//vweOZd18gpO2w5mS4CJ7t6t\nP/+jjgcPuvvE6PNId18drV8MHOnuX446EtxBePA/GvgrcKC7u5n9HfgmsACYC/zS3f9iZhdEdfua\nmc0CTnP3WVFHgheByYTnRwsJacHN5DAzd3eeeQY++Un47vgH+PFbM+Dkk+HBBxERke2ZGe7erQGg\nu/Jy6DuE9NqHDjpmNgf4FLCHmb0PXAlMNbPDCb3YlgH/BuDui8zsj4QWVQq4IOel1AuAW4AK4KGc\nZ0o3Areb2VJgAzArOtdGM/sB8EK031XtA06ubHpta4vSayIi+dSVls6fCOOsPUZb4HF3/2ae65Z3\n2ZbOkiUwfjycMepJ7lg1NTR7nupW9lBEpM/Ld0vngWjJRifLWe8Tsr3XNjepy7SISD7tNOi4+y1m\nNgAY6+6LC1Cngsum11qDjtJrIiJ50ZXpqk8FXgb+En2eZGYP5LtihZQNOpsa9J6OiEg+deU9ndnA\n0UQjC7j7y8D+eaxTwSUSUFoKjRplWkQkr7oSdJId9Pzq6thru43KypypDdTSERHJi64EnTfM7Ewg\nYWbjzOxXhInd+pTKSk3iJiKSb10JOhcChxC6S88BtgLfymelimHgQE3iJiKSbzt9T6cvy76nA3DE\nEfDSwgwZ4mFjJgPWrW7oIiJ9Wl7f0zGzJzoodnc/vjtf2FtVVoITI5MoIZZKQksLlJUVu1oiIn1K\nV14O/W7OejnwBcIwNX1Kttt0uqQ8BJ3mZgUdEZEe1pWXQ19sV/SMmb3Q4c67sdygU9JYGzoTDBpU\n3EqJiPQxXUmvDc35GAOOAPrcb+Ns0EnF1YNNRCRfupJee4m2sdZSwHvAefmqULFkx19LJdSDTUQk\nX7qSXtu3APUoumxLpyWmF0RFRPKlK+m1L7D9qNLZrnLu7n/q8VoVQTboJE3pNRGRfOlKeu1c4GOE\n6Z6NMA30c0B22uk+FXSaY0qviYjkS1eCTikwIWeK6ZHAre5+Tl5rVmCtQUdD4YiI5E1XhsEZA6zJ\n+bwWGJuf6hRPNuho0E8RkfzpSkvnr8AjZnYHIb02E5iX11oVQbb3WpMrvSYiki9dCToXAp8DPhl9\n/q2735u/KhVHtqXTmFF6TUQkX7rSZdrN7CWg1t3nmdkAM6ty99oC1K9gskGnIaOWjohIvnRluur/\nDdwN/CYq2hu4L5+VKoZs0KlP65mOiEi+dKUjwdeBTxDm0cHdlwDD81mpYmgNOiml10RE8qUrQafZ\n3VtzTWaWYPuXRXd72aBTl1J6TUQkX7oSdJ40s8uAAWb2aUKq7cH8Vqvwsr3XalNKr4mI5EtXgs73\ngPXAa8C/AQ8Bl+ezUsXQ2tJpUXpNRCRfdth7LUqlve7uBwP/VZgqFUc8HuZsa2xWek1EJF922NJx\n9xTwlpntU6D6FFVlpUYkEBHJp668HDoUeMPMFgD1UZm7+6n5q1ZxVFZC80al10RE8qXToGNmZVGv\ntctpm8ogq8/1XoMQdOqIehRs3VrcyoiI9EE7auk8B0wGznf3swpUn6IaOBBWMzJ8WL26uJUREemD\ndhR0yszsTOBjZvZ5QmvHsz/7yuRtuSorYQWjwodVq4pbGRGRPmhHQeerwJnAYOCUDrb3yaCzTUvH\nHax9ZlFERLqr06Dj7k8DT5vZi+7++wLWqWhC77UKWiqHUFq/GTZsgD32KHa1RET6jJ2+HNpfAg7k\njL82SM91RETyoSsjEvQb2aBTW6XnOiIi+ZDXoGNmN5nZWjN7LadsqJnNM7MlZvaomQ3J2fZ9M1tq\nZovNbFpO+RQzey3a9ouc8jIzuysqfz73JVYzOzv6jiVm9s9dqW92/LUtlVHQUUtHRKRHdWU+nZiZ\nfcXMrog+jzWzo7p4/puB6e3KLgHmuftBwGPRZ8xsAmEq7AnRMdebtT7FvwE4z93HAePMLHvO84AN\nUfnPgOuicw0FrgCOipYrc4NbZ7ItnU3lUXpNLR0RkR7VlZbO9cCxwJejz3VR2U5FnRE2tSs+Fbg1\nWr8VOC1anwHMcfeku78HvA0cbWYjgSp3XxDtd1vOMbnnugc4IVr/DPCou292983APLYPftvJBp0N\npWrpiIjkQ1eGwTna3SeZ2csA7r7RzEp24TtHuPvaaH0tMCJaHwU8n7PfSmA0kIzWs2qicqKf70f1\nSpnZFjMbFp1rZQfn2qFs0FkXV0tHRCQfutLSaTGzePaDme0JZHriy93d6UVD6mSDzpqYOhKIiORD\nV1o6vwLuBYab2Y+AL7Jr8+msNbO93H1NlDpbF5XXAGNy9tub0EKpidbbl2ePGQusiqZhGOzuG8ys\nBpiac8wY4PGOKjN79uzW9Xh8KjCVVa4u0yIiWfPnz2f+/Pk9ci4LjY2d7GT2Edqelzzm7m92+QvM\n9gUedPeJ0ecfEx7+X2dmlwBD3P2SqCPBHYQH/6OBvwIHurub2d+BbwILgLnAL939L2Z2ATDR3b9m\nZrOA09x9VtSR4EXC2HEGLAQmR893cuvmudc/bx5MmwYnTW1g7vxKKC0No01rVAIRkVZmhrt36xfj\nTls6Zna7u38FeLODsp0dOwf4FLCHmb1P6FF2LfBHMzsPeA84HcDdF5nZH4FFQAq4ICciXADcAlQA\nD7n7X6LyG4HbzWwpsAGYFZ1ro5n9AHgh2u+q9gGnI9n02samATBkCGzeDBs3wrBhOztURES6YKct\nHTN72d0n5XxOAK+6+4R8Vy7f2rd0XnkFDj8cJk6EV1MT4M034dVXQ4GIiAC71tLptCOBmV1qZrXA\nRDOrzS6EZzAPdLOuvVq2pVNXB4xSt2kRkZ7WadBx9x+5exXwE3evylmGuvslBaxjwbSOvVYPjFS3\naRGRnrajmUMPdvfFwN1mNrn9dnd/Ka81K4Jtgs4odZsWEelpO+pI8O/A+cBP6fhdmuPyUqMiygad\nhgbwvUaGObqVXhMR6TE7mk/n/Ojn1ILVpsjicSgvD72kW/YYRRmopSMi0oO6MuDnl8xsULT+/5nZ\nnzpKt/UV2dZO4xC9ICoi0tO6MgzOFe6+1cw+QXhB9CbgN/mtVvFoTh0RkfzpStBJRz9PBn7n7n8G\ndmXAz14tG3S2Vua0dLowaoOIiOxcV4JOjZn9F2Gum7lmVt7F43ZLre/qZAbA4MHQ0hJGJRARkV3W\nleBxOvAIMC0aSqYa+G5ea1VE2dlDt+k2rec6IiI9YqdBx93r3f0ed18afV7t7o/mv2rFoRdERUTy\np8+mybpLL4iKiOSPgk4724y/NlLdpkVEepKCTjtq6YiI5I+CTjsdBh21dEREeoSCDvD1uV/noF8d\nxN9X/n3b3mvqSCAi0qMUdIA19WtYunEpy7csV0tHRCSPFHSA0VWjAVi5dWXnXabT6Y4PFhGRLlPQ\nAfYetDcANVtrGDQolG3aBAwYAGPGhFEJli4tXgVFRPoIBR1yWjq1K9l//1DWGmMmRwNqL1xY+IqJ\niPQxCjps29I56KBQtnQpZDLAlCmh4KU+N1GqiEjBKejQFnRWbl3JkCEwfDg0NsL779PW0lHQERHZ\nZQo6wKho7pxVtavIeIbx40P5kiVsG3QymeJUUESkj1DQASpKKhhWMYxkJsn6+vWtQeettwg92Pba\nC7ZuhWXLilpPEZHdnYJOZPSgtm7T2wQdaHuuo84EIiK7REEn0tqZoLZm+6Cj5zoiIj1CQSeyd1Vb\nZ4JsDzYFHRGRnqWgE8mm12q21rD//pBIwIoVoRfbNkHHvXiVFBHZzSnoRFq7TdeupKSEbV8SHTMG\n9tgDNmwIkUhERLpFQSeSO/4asO1zHTOl2EREeoCCTiR3VAJAz3VERPJAQSeSOyqBu2/7gigo6IiI\n9AAFncigskFUllRSn6xna/NWdZsWEckDBZ2ImW3T2skNOu6EngWDB8OaNZpJVESkmxR0cuSOSjB8\neIgxW7bAunVs25lgwYLiVVJEZDemoJMjd1QCsw46Exx/fPh5772Fr5yISB9QtKBjZu+Z2atm9rKZ\nLYjKhprZPDNbYmaPmtmQnP2/b2ZLzWyxmU3LKZ9iZq9F236RU15mZndF5c+b2T47q1PuqATA9p0J\nTj89/LzvPmhq2qXrFxHpj4rZ0nFgqrtPcvejorJLgHnufhDwWPQZM5sAzAQmANOB683MomNuAM5z\n93HAODObHpWfB2yIyn8GXLezCuWOSgBs35ngoIPg8MPDiNOPPNLNyxYR6b+KnV6zdp9PBW6N1m8F\nTovWZwBz3D3p7u8BbwNHm9lIoMrdsw9Zbss5Jvdc9wAn7KwyuaMSQAdBB2DmzPDzrrt2djoREWmn\n2C2dv5rZi2Z2flQ2wt3XRutrgRHR+ihgZc6xK4HRHZTXROVEP98HcPcUsMXMhu6oQjsclSArm2J7\n4IFoYDYREemqRBG/++PuvtrM9gTmmdni3I3u7maW99E1Z8+e3bp+2NGHAW3ptQMPhFgM3nkHamuh\nqorQdfrII+GFF+Chh+ALX8h3FUVEimr+/PnMnz+/R85VtKDj7qujn+vN7F7gKGCtme3l7mui1Nm6\naPcaYEzO4XsTWjg10Xr78uwxY4FVZpYABrv7xvb1yA06Gc9Q8mIJGxo30JhsZMCACo45Bv72N5g3\nDz7/+WjHmTND0LnrLgUdEenzpk6dytSpU1s/X3XVVd0+V1HSa2Y2wMyqovVKYBrwGvAAcHa029nA\nfdH6A8AsMys1s/2AccACd18DbDWzo6OOBV8B7s85JnuuLxI6JuxQzGKMqhoFwKra8ALoKaeEbQ8+\nmLPjl74Ufv75z1BX1/ULFxHp54r1TGcE8LSZ/QP4O/Bnd38UuBb4tJktAY6PPuPui4A/AouAh4EL\n3FsntrkA+D2wFHjb3f8Sld8IDDOzpcC3iHrC7UzuqAQAJ58cyufOhUwm2mnsWDj22PBM589/7s71\ni4j0S0VJr7n7MuDwDso3Aid2csyPgB91UL4QmNhBeTNw+oetW+6oBACHHAL77APLl4eM2tFHRzvO\nnAnPPQc33QSzZn3YrxER6ZeK3WW618m+IFpTGzoTmHWSYjvrLBg4MDzsee65AtdSRGT3pKDTTvv0\nGrSl2LbJpA0bBhddFNavvLJAtRMR2b0p6LRzwNADAHhx1YutZVOnQmUlvPJKu9mqv/1tGDQotHae\neaawFRUR2Q0p6LRz/H7HUxYv4/mVz7O2LrynWlYG06LR3ubOzdl56FD41rfCulo7IiI7paDTzsDS\ngZy4/4k4zoNL2h7idPhcB+Dii8McCI8/Dk8+WbiKiojshhR0OnDawWH4tvsW39dadtJJoVPB449D\nfX3OzkOGhMADcMUV0YxvIiLSEQWdDpxy0CkYxl/f/St1LeHlzxEj4KijoLk5DLu2jW99C6qr4amn\n4Pe/L3yFRUR2Ewo6HRgxcATHjjmW5nQzj7zdNoXBOeeEn1dfDel0zgGDB8Ovfx3WL74Y3n67cJUV\nEdmNKOjNfFEIAAAWrklEQVR0Ysb4GQDc91Zbiu2cc8JgBIsWwd13tzvgy18OL4nW14d3eFKpAtZW\nRGT3oKDTiexznT8v+TPJdBKA0lK4/PKw/aqr2rV2AK6/HvbeG/7+d/jRdoMniIj0ewo6nTho2EEc\nvMfBbG7azNMrnm4t/5d/gf32g8WL4c472x1UXQ233BLW/8//0bs7IiLtKOjswGnjt+/FVlKybWtn\nuyzaCSfAd74TmkEzZrSbAU5EpH9T0NmBbIrt3sX30pJuaS3/ylfggANg6VL4wx86OPCaa8LYORs3\nwmc/C2vXdrCTiEj/o6CzA0eOPpLxw8azcutKrn3m2tbykpLwSg500lktkQi5tyOPhGXL4H/9L827\nIyKCgs4OxSzGb0/+LQBXP3U1i9Yvat121lkhe7Z5M3zucx3ElMrKMHzB/vvDwoUwfTqsWVPA2ouI\n9D4KOjvxqX0/xfmTzyeZSfKvD/wrGQ8zucVicNttMH48vP46nHtuB4MRjBgBDz8Mo0bBs8/ClCma\nBkFE+jUFnS748ad/zMiBI3lu5XPc8MINreWDBsF990FVVXhv58c/7uDggw4KLZ1PfhJWrYJPfSp0\nrdZwOSLSDynodMGQ8iH8v5P+HwCXPHYJb29se4hz8MFtnQkuuQR++tMOTrDXXvDYY2H+nWQSvv71\nMILo6tUFqL2ISO+hoNNFn/vI5/jihC9S11LHp2//NDVba1q3nXoq/PKXYf0734Hvfa+DhkxJCfz8\n56GDwZAhYY6EQw+Fu+4q3EWIiBSZeT9O85iZf5jrr22u5cTbT2RBzQIO3uNgnvqXp9izcs/W7f/9\n3+Hl0VQqDJnz29+GWLOdmho47zx4JBrX7fjjQzfro47atQsSESkAM8PdrTvHqqXzIVSVVfHwmQ9z\n6PBDWfzBYj7zh8+wuWlz6/YzzwwjUFdUwM03w8c/DkuWdHCi0aNDB4Mbbgitnscfh6OPhs9/Ht54\no3AXJCJSYAo6H9LQiqHM+8o8Dhx6IC+veZnJv528zUjUn/0sPPFEGBj0hRdg0iT4r//qIN1mBl/9\nKrz7bngYVFEB994LEyfC6afDa68V9sJERApA6bVuXv/yzcs59c5TeXXtqwCccegZ/OwzP2PEwBFA\neH/nG98IKTeA446D664L74t2aPVq+OEP4Xe/g5Zo9IMZM0Jg+vSnIR7vVj1FRHrarqTXFHR24fqT\n6SQ/e/5nzJ4/m8ZUI1WlVVx09EV8+9hvU11RDcCcOaGz2qZN4ZgvfjHMxzN+fCcnrakJfa9/+9sw\nYxyEZtM554QW0Ec+ElpJIiJFoqDTTbsadLKWbVrGhQ9fyNylcwEYVDaIi4+5mIuOvojqimo2bYJr\nrw093JqawjHTpsHXvhaGaEskOjjpmjVw001w440hBZd14IGhBfT5z8Mxx4S3VEVECkhBp5t6Kuhk\nPbviWa6cfyWPLXsMgIGlA/naEV/j4mMuZmTVSFauhB/8IIxkkA0+o0eH+d/OOAMOP7yDRkwmA/Pn\nw+23h2F1Nmxo2zZmDMycGYLQxIlhBlMRkTxT0Ommng46WU8tf4qrn7qaee/OA6AsXsYp40/hpANP\nYvqB0ylLjuTWW+E3v9m2d9v48fClL4Wx3CZN6iAApVLwt7+FYRDuvhtWrtx2+157wYQJodvcccfB\nscdCeXmPX5+I9G8KOt2Ur6CT9ULNC1zzzDXcu/jebconj5zMSQeexGcPPInk8qP4451x7r4b1q9v\n22fMmDBowQknwNSpMHRou5NnMmEctzvvDJPFvfUWNDZuu09ZWWg+HXJIWA47LPRkUItIRHaBgk43\n5TvoZL23+T3mLpnLw28/zOPLHqcx1RYcqsurOXL0kUwaMYWSD6bw7t8O5/F79mPN6rZnNWbw0Y+G\nV3mmTAnLoYeG6bNbZTKwfDm8+mpIxz3xBLzyyvaVMQtj9xxzTAhIhx0WUnPDhuXvBohIn6Kg002F\nCjq5GpONPLn8SR5a+hBzl87l3U3vbrfPgJIB7FMxgbKth7Dl3YN5/6WDSa05GDYeCJnQ66C0NMSL\nyZNDzBg1KmTXRo4Mnd3iccIkcq+9FobBfuMNeOklePnlti7ZufbcM0zDsP/+YT7usWNhn33CcsAB\n7SKciPRnCjrdVIygk8vdWb5lOQtXLWTh6rC8vu51VtWu6nD/BGUMbPoI6VUTqV02HjYdAJv2h837\nQONQyIQxdyoqwqOdiRNh331DI2aPPUJAOuTAZvZY+Q9YsCAEpFdfDT8bGjqvaCIRWkcTJ4aANGQI\nVFeHNN3AgWGY7aqqsG3QoDzcKRHpTRR0uqnYQaczGxs38vq613lz/Zss/mAxb214i0XrF7F8y/Id\nHhdPDoLGPUhvHgm1o8KydW/YMga2jIW6vSA5gBHDKph48AD22yfBqFEwemSGfUpXs1fjMvasfZfB\nG5dRsX4FtmJ5mPl02bKuT8Ww997h+dHYsSEQDRwYAtFee4Xm2MiRMHx4CFzq7i2yW1LQ6abeGnQ6\ns7V5K6+ve53X173O0g1LeXfzu7yz8R1Wbl3JxsaNOB/yWhqroXYk1I2EuhHQsAc07AmN1Viqksqy\nCgZXDGC/wXGOLVnHlPRKRiVrGdjYQnn9VkqbtlCeqqOspZbyxs0MWPMOsZbmrn13LBaaYFVVoSWV\nSIQU3h57hKC0554hQI0eHQJZdXV4HmUW9h06NCwdjqgqIvmkoNNNu1vQ2ZGMZ9jctJn19etZU7eG\nVbWrqKmtYeXWlby/9X1WbFnBuvp1NCYbaWhppDHVQIZM97+wuQqaBoNHrRWPEW8YzH5rB3Loujij\nm+IMajGqWoyhqRSjU7WMTG5leNNWhjTVUdW8g3TehzF4cAhIgwa1LQMHhqWyMgSykpKwVFaG7VVV\noaU1fHiY3XXPPUNOUiM9iHSJgk439aWg82FlPMOGhg2srlvN6trVrKtfxwcNH/BBwwdsbNxIfbKR\nLfUNbK5vYFN9PZsbaqlrqaWRzTTHNuCW3qXvT6RhjwaobIG4QyID5alQNrw+LCNrYXQt7L0VBjeB\nAeZQkoGhjWGJ9+B/vlQiQaokdNSIpdPE0iEoJyvKaRlQTqqygtTwoWT23gsfNQpKy4g1NGD1DVgy\nRaxyIImBVSQGDiZeVk68vIJ4eQVeUkq6NEGyJEamJEFZopzSRBkWi4Vu7RUVYUkk2gJfLBaCZDZ4\nuocJAJPJsE/2GI3JJ0WgoNMJM5sO/ByIA7939+vabe+3QWdXuDtbmrewpWlLa0ovnUmzqWkT6+vX\ns3LzOtZu2Uxtcy1bm2upbaqnvjFNY2OG+sY0zckUTakWWlJJWjKNJKknaXW0UE8qkyTpLaRpJk0S\nt7AQywlyloZYBstAdVMISIObYVC0VLbAwGgpyUBJGkrTMCAZtle1QHUjjKiHEXWwZwOU7VoMLZqW\nmNGciNGcgKaE0VRiNJUkaClJkEzEtkm4ZsxJG6RjELM4FislbqV4IkFjaYymEqO5BGKWIPs/j8XI\nGKTiRjphpMpLSFeUki5PYPGw3WIxPJHAS+NkyhLEEglKMnFK0jESaSOZaiGZaqIl1UwmESM2sJJY\nZSWJygFYIo4lYhCPkW5oJl3XSKauCSuJERtcScmQgcQGlJGqayRV20iqoREvjZOuLMEHlhCvqGBg\nWRUDywZRUTqAlkySplQTTalmLBEnUVpOvKQMSktpLomRjIc/uEriJZTGSigjQcKNmEMiGnQ/ZU46\nZqRjkIiXUBIvoSRWguOkM2nSHv7PkoglKImVEI/FyXiGVCZFKpMiZjESscQ2S9zixCxGU6qJ+mQ9\n9S31recBKImVMKR8CEPKhzC4fDAt6RYakg00JhtJe5q4xYnH4rg7DcmGsC3VSMxilMRKKI2Xttaz\nJF5CIpbA3Vv/fcYs1noOw1rr6jil8VLK4mWUJcpw99ZtLekWmlJNNKebaUo1MX7YeEYPGg0o6HTI\nzOLAW8CJQA3wAnCGu7+Zs4+CTmT+/PlMnTq12NXoVCYT/tgPf/A7m+qaWL+llvVbaqlrTFLXkKKx\nKc3mhkY2NWxmU+NmtjRvwT2D42TIkEylaE4maUq20NDSTG1jE/UtjTS2NOEZpzSVoTTjNK1dTdne\nI0knANKUNjZT1tRMRWMTezZsZlTLZkY3byXhTn0iTl0iQTLuVGSaqfRmBmSaKfU0pZl0+JkyylIx\nypMxStNglsYsQ8xDMKxIQkUqtPay4hmoTLYFz4xBMgbJeGjtVaTCcfnuijEfmJrn7yikZAya4+GP\nka78oZGycM+TMXgmE+5F3MGBljg0J8I2gJiHxQ1SMUhHv5LLU2EpS4fylnhYIPw3T2TCcRkL581Y\n27mz+7ZEdchYVPdU+P8OhO/LWDhHPDqfEa6zKTpHIrreslTbdot+9WVs2yUda6tLblh5/gv/m/Nu\n+C2wa0Gno6Em+4qjgLfd/T0AM7sTmAG8uaOD+qveHnRyO7olEkZFRQWj9qwAhvf4d82ePZvZs2d3\nuM0damvbhsAbMCAspaUh89XcHJZMBtLpsKRSbeXZ7FiGNE3pBppb0qxvdhobnYx76FMRB8dpbDDq\n6436eiCWIlbSQqykhSSNNCQbaWiuo6WpnkR9CVaXIFZneFMj1rwJmrdgqUZisRiGRb9kYsQzEMs4\nKW+kxbeSYiuebqS8BcqbndJkGrckKUuSiSd5YvmbPD/6IGIZSKTSlCZbKG1poSSZJEaGmKeJeYa4\npylJpyhNp4l7mmQc0vHwSzZjCfCwJDxNeaqZ8nQz5ekWYu7EyWDutMQStMRLaEqUEvMMA5LNDEi1\nUJZJ0RQroSWeoCVWQkkmTUUqxYBUklJPYZbBSBMjAxh4uOKYOwnPkPAM5Zk0FakMJRmnJCe4pw1S\nZmTMskcTdw8pX3cSDokUVAALgJNy/r9QngY6eOVtR8rS4Y+J3dErG3omlduXg85o4P2czyuBo4tU\nF+kjzNr6K7RXVhYewXRNHKjqwZrlR3r2bC7pJAAXinvHfTzcQ1DPtoDdwyOuWCws6XR4D7qlBRqa\nYWOj01SbJFnfQqyshHh5CfGSWGvnyewjtey5Mmkn2Zgi1Zgk2ZDkg99fy9KvX068NA7upBtbSDW0\nkG5sATMsHsNihmccT6XxZAqARFUFpYPKSQwoJZNMk2poIVXfDGYkyhMkKhLEYkYm42SSGdLJDMn6\nFlpqm0nWNZNuSpJpTuLNLWRSGWLlpcQrSomVl4Y/xtwxzxBLxIiVJkiUxbGYkawLx6fqm0N5ZRkl\nA8uIlYYLtZiFa804ns6QSaYhkyGTyuDpDGQyrfcC4LgjR/bIf8++HHSUNxPpAzrrVJjtPd+ZeLyt\nv0V0BFAaLV36ZqAkWmCPh8sYd3hlzvaKjg7qgsqd79KH9eVnOscAs919evT5+0AmtzOBmfXNixcR\nyTN1JGjHzBKEjgQnAKsIKdltOhKIiEhh9dn0mrunzOwbwCOEBPqNCjgiIsXVZ1s6IiLS+/TLERfN\nbLqZLTazpWb2vWLXp5DMbIyZPWFmb5jZ62b2zah8qJnNM7MlZvaomQ0pdl0LxcziZvaymT0Yfe6X\n98LMhpjZ/5jZm2a2yMyO7sf34vvRv5HXzOwOMyvrL/fCzG4ys7Vm9lpOWafXHt2rpdHv1Gk7O3+/\nCzrRS6O/BqYDE4AzzOwjxa1VQSWBi939EOAY4OvR9V8CzHP3g4DHos/9xUXAItp6PPbXe/EL4CF3\n/whwGLCYfngvzGxf4HxgsrtPJKTnZ9F/7sXNhN+PuTq8djObAMwk/C6dDlxvZjuMK/0u6JDz0qi7\nJ4HsS6P9gruvcfd/ROt1hJdlRwOnArdGu90KnFacGhaWme1NeOfv94Q+stAP74WZDQY+6e43QXgm\n6u5b6If3AthK+ONsQNQhaQChM1K/uBfu/jSwqV1xZ9c+A5jj7snoRfy3Cb9jO9Ufg05HL42OLlJd\niir6i24S8HdghLuvjTatBUYUqVqF9jPgu7DNkNv98V7sB6w3s5vN7CUz+52ZVdIP74W7bwR+Cqwg\nBJvN7j6PfngvcnR27aMIv0Ozdvr7tD8GHfWcAMxsIHAPcJG71+Zuiwak6/P3ycxOBta5+8u0tXK2\n0V/uBaEn62TgenefDNTTLn3UX+6FmR0AfAvYl/BLdaCZnZW7T3+5Fx3pwrXv8L70x6BTA4zJ+TyG\nbSN1n2dmJYSAc7u73xcVrzWzvaLtI4F1xapfAX0MONXMlgFzgOPN7Hb6571YCax09xeiz/9DCEJr\n+uG9OAL4m7tvcPcU8CfgWPrnvcjq7N9E+9+ne0dlneqPQedFYJyZ7WtmpYSHYA8UuU4FY2YG3Ags\ncvef52x6ADg7Wj8buK/9sX2Nu1/q7mPcfT/Cg+LH3f0r9M97sQZ438wOiopOBN4AHqSf3QtCB4pj\nzKwi+vdyIqGjSX+8F1md/Zt4AJhlZqVmth8wjvAifqf65Xs6ZvZZ2ubZudHdrylylQrGzD4BPAW8\nSlsz+PuE/6P8ERgLvAec7u6bi1HHYjCzTwH/7u6nmtlQ+uG9MLOPEjpUlALvAOcQ/o30x3vxH4Rf\nrhngJeBfCSO09vl7YWZzgE8BexCe31wB3E8n125mlwLnAilCuv6RHZ6/PwYdEREpjv6YXhMRkSJR\n0BERkYJR0BERkYJR0BERkYJR0BERkYJR0BERkYJR0BGJmNkphZ7qwszmm9nkaP296B2hHe2/030K\nycxm9LNR2mUXKeiIRNz9QXe/blfPE41M3OWv7WS9K/v3Bp8jDGsv0iUKOtLnRUMeLY5GUH7LzP7b\nzKaZ2bPRpFRHRvv9i5n9Klq/xcx+Ee3zjpl9ISo3M/u/0eRer5rZ6VH5VDN72szuJwwf074O15vZ\nCxYmzpvdhTqfZWZ/tzC53G86mqOki/scGV3DP6J9K82sPLoXr0YjSk9tf/3R5z+b2T9F63VmdnV0\nnufMbLiZfQw4Bfi/UR3278J/DunnFHSkvzgA+AlwMDAemOnuHwe+A1zayTF7RfucDFwblX0e+Chh\nkrMTCb9w94q2TQK+6e7jOzjXZe5+ZHTsp8xsYmcVjdJVpwMfc/dJhKFYzuzGPqWE+aK+6e6HAycA\nTcDXgbS7HwacAdxqZmVs34rK/TwAeC46z1PA+e7+N8LYW99x90nu/m5n1ySS9WHSACK7s2Xu/gaA\nmb0B/DUqf50whH17TjSoobu/aWbZ+UM+AdwRDe++zsyeBI4kTPy1wN2Xd/L9M83sfMK/uZHAR4DX\nOtjPCMFhCvBiGG+SCmDNh9wHQnBd7e4Lo+uoi67/48Avo7K3zGw5cBA71uLuc6P1hcCn29VHpEsU\ndKS/aM5ZzwAtOeud/TtoyVnP/mJ1tv8lm20R1Hd0kmj03X8HjnD3LWZ2M1C+k/re6u6dtcA+zD6d\n6egaUmyb/citYzJnvf09623PmaQXU3pN5MN5mtBqiZnZnsA/EUbo3tFf+4MIAWlr1GL67A72dcIc\n9F+Mzo+ZDTWzsR9yH4C3gJFmdkS0T5WZxaNrODMqO4gwcvBbhNGDD4+eW41hJ9MOR2qj6xPpErV0\npL/Y0fMKz/nZWW8yB3D3e83sWOCVqOy77r4uesbS4V/87v6Kmb1MmKflfeCZHVY0pPMuBx6NOgck\ngQsI0yd/mH1azGwm8CszqwAaCM+hrgduMLNXCa2bs909CTxrYUK7RcCbhDRaZ/ci+/lO4HdmdiHw\nJT3XkZ3R1AYiIlIwSq+JiEjBKOiIiEjBKOiIiEjBKOiIiEjBKOiIiEjBKOiIiEjBKOiIiEjBKOiI\niEjB/P+pD5KVVuG+3wAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "pop1_sfs_folded_scaled = anhima.sf.site_frequency_spectrum_folded_scaled(pop1_biallelic_ac)\n", "pop2_sfs_folded_scaled = anhima.sf.site_frequency_spectrum_folded_scaled(pop2_biallelic_ac)\n", "pop3_sfs_folded_scaled = anhima.sf.site_frequency_spectrum_folded_scaled(pop3_biallelic_ac)\n", "\n", "fig, ax = plt.subplots()\n", "anhima.sf.plot_site_frequency_spectrum(pop1_sfs_folded_scaled, ax=ax, label='pop1', plot_kwargs=dict(lw=2), m=m)\n", "anhima.sf.plot_site_frequency_spectrum(pop2_sfs_folded_scaled, ax=ax, label='pop2', plot_kwargs=dict(lw=2), m=m)\n", "anhima.sf.plot_site_frequency_spectrum(pop3_sfs_folded_scaled, ax=ax, label='pop3', plot_kwargs=dict(lw=2), m=m)\n", "ax.set_ylim(bottom=0)\n", "ax.set_ylabel('scaled site frequency')\n", "ax.set_xlabel('minor allele frequency')\n", "ax.legend(bbox_to_anchor=(1, 1), loc='upper left');" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAEKCAYAAADO/lZMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd4FFUXxt+7KbQgSYAAoQUQRRA/ihUEEQQBpVtQUFQs\niIIFKSJIaAKKKB2kCYhAEETpRYogvYr0DiGhBQIhJJDsvt8fd2uySTawmw3J+T3PPJm9c2fmzGR3\n3jn3nnuuIglBEARBELyDwdsGCIIgCEJuRoRYEARBELyICLEgCIIgeBERYkEQBEHwIiLEgiAIguBF\nRIgFQRAEwYv4etsAb6KUkrFbgiAIdwBJ5W0bcgq53iMmCc6YAbZtC0ZH68+5cOnXr5/Xbcgui9wL\nuRdyL9JfBPeS64UYAPDLL8CcOcDatd62RBAEQchliBADQIMG+u+aNd61QxAEQch1iBADQP36+m8u\nFuJ69ep524Rsg9wLG3IvbMi9EDyFys3t/UopkgSMRqBIESA2Fjh5EggL87ZpgiAI2RalFCjBWm5D\nPGIA8PEBLG+7udgrFgRBuBuUUpQl7SWt+yZCbEGapwVBEO4ab0d0Z9clPUSILViE+K+/gFzcXC8I\ngiBkLSLEFipXBooVA86fBw4d8rY1giAIQi5BhNiCUtI8LQiCIGQ5IsT22DdPC4IgCEIWIEJsjyWx\nx7p1ekiTIAiCIDjh/PnzaN68OUqWLAmDwYAzZ87c8bFEiO0pV06PIb56Fdi719vWCIIgCNkUg8GA\npk2bYv78+Xd/LDfYk7OwNE+PHQucO+ddWwRBEIS7JiwsDEOHDkWVKlUQHByMd955B7du3QIATJo0\nCRUrVkThwoXRokULREdHW/czGAwYPXo0KlSogKJFi6JHjx7WoUghISHo1KkTHn300bu2T4Q4JY0b\n679TpwKlSgFPPAFMmACYTN61SxAE4R5GKfctd8Kvv/6KlStX4vjx4zhy5AgGDRqENWvWoHfv3pg3\nbx6io6NRtmxZtG3b1mG/hQsXYufOndi1axf++OMPTJ061Q13wxFJcZny+k0mYO5cICICWLECSEjQ\n5XXqANOmARUqZL2hgiAI2Yi0Ulw6faZat7nv/JmVrXLlyuHLL7/E+++/DwBYtmwZunTpgnr16qFI\nkSIYOnQoACA+Ph5BQUE4duwYypQpA4PBgOXLl6NRo0YAgPHjx2P+/PlYvXq19djJycnw9/fHqVOn\nUKZMmTRtSC8tqHjEKTEYgNdeA37/Hbh0SU+RWKwYsGED8L//AT/+CPz7r85LLQiCILgE6b7lTihd\nurR1vUyZMoiKikJUVJSDeBYoUACFCxfGObtuSWf7uRsR4vQoUABo1w7Yvx9o2xaIjwc++0wLclAQ\nUKgQUKsW0KULMH06cPSoty0WBEEQnGAf1XzmzBmEhoYiNDQUp0+ftpbHx8cjJiYGJUuWTHM/+23u\nQpqmM3P98+cDU6YAp04Bp08DN2+mrlOpEtCyJdCmDeCGTnxBEITsxp00TXuTsLAwFCpUCEuXLkW+\nfPnQvHlz1KtXD88++yxee+01rFq1CpUqVUKPHj2we/du/P333wB0sFaDBg0wb948xMXFoVGjRujW\nrRveffddAEBiYiKSk5Nx33334dChQyhTpgzy5s3r1Ib0mqZFiO/0+kng8mVgzx5g505gxw6dkevq\nVVudjz4CfvgB8PNzj8GCIAjZgHtNiMuVK4dOnTphxowZiIqKQsuWLTF+/HjkzZsXEydOxHfffYer\nV6+idu3amDBhAkJDQwFoIR41ahR++OEHXLt2DW+//TaGDRsGg8Fg3Q5Y7weUUjCmkYNChDgN3P6l\nSUrSfcm//w789BNw+7ZOEhIRAQQHu+88giAIXuReFOIpU6agvmV4qosYDAYcO3YM5cuXv2sb0hNi\n37s+umDDz0+PQ65fH3j9daBVK50u84kngPfeAxIT9VKmDNChA5Avn7ctFgRBELyMCLGneOopYNs2\noEUL3Xzds6fj9gEDgD59gHffBfz9vWOjIAiCkCbKnWOu0jtPdmxGyCqypBklPh4YPVr3J+fLp0V3\n4UJg1y69vVQpoG5doEoVPRVjgwZAwYKetUkQBOEuuNeaprMD0kecBl770pDAggVA377AwYOO28qU\n0X3MNWpkvV2CIAguIEKceUSI08DrXxqjEdi+XY9TPnAAWLkS+O8/IG9eYNIkoH1779kmCIKQBiLE\nmUeEOA2y3Zfm1i095GnKFP25c2fg6691Zi9BEIRsgghx5hEhToNs+aUh9dCnLl30cKg8eYA339QZ\nvSpVcm/CVkEQhDtAhDjziBCnQbb+0uzcqSOr//zTVpYvHxAaqpf33gPeeMNxH1Jn/QoLE8EWBMFj\niBBnHpn04V6kZk3gjz+AQ4eADz7Qea0TEoDjx3XSkDffBGbPttU3GoGOHYHy5XWKzfh479kuCIIg\nuIwIcXbnwQf1fMixscD161qYBwzQ2zp0ANau1SL89tt6mkZAe9F16wIemCVEEARBAJYsWYKnn34a\nQUFBKFGiBN577z3cuHHjjo7lMSFWSpVWSq1VSu1XSv2nlOpqLg9WSq1SSh1RSq1USgXa7fOlUuqo\nUuqQUqqRXXlNpdQ+87aRduV5lFJzzeVblFJl7bZ1MJ/jiFLqTU9dZ5ZSsKAW5r59gU8/1X3IrVoB\nrVsDM2fq2aKmTdNe8a5dOqPXDz8AvXrpZuxu3YC4OG9fhSAIwj3P9evX8fXXXyM6OhoHDx7EuXPn\n0L179zs7GEmPLACKA6hmXg8AcBjAQwC+BdDDXN4TwFDzemUAewD4AQgDcAy2PuxtAB43ry8F0Ni8\n3hnAOPP6qwDmmNeDARwHEGhejgMIdGIj71mMRrJNG9sUnQEB5IYNetvFi2StWs6n8qxfn0xI8K7t\ngiDc05ifnc6e+940K03Kli3LIUOGsHLlygwKCuLbb7/NxMREkuRPP/3E+++/n8HBwWzevDmjoqKs\n+ymlOGrUKJYvX55FihRh9+7daTKZnJ5jwYIFrFq1apo2pHXPSHouxSXJ8wDOm9dvKKUOAigJoDmA\nZ8zVpgNYB6AXgBYAZpNMAnBKKXUMwBNKqdMACpLcZt5nBoCWAJabj9XPXD4fwBjz+vMAVpKMBQCl\n1CoAjQHM8czVegGDQXvBcXHA7t06W1etWnpb0aI6x/W33wLR0UDJkkBICNCvn54h6vXX9UQUvpLh\nVBCErEH1d18AKftlPiDs119/xcqVK5E/f340a9YMgwYNwrPPPovevXtj1apVqFy5Mr744gu0bdsW\n69evt+63cOFC7Ny5E3FxcXjuuefw4IMPomPHjqmOv379ejz88MN3dD1Z8iRWSoUBqA5gK4BiJC+Y\nN10AYBkkGwpgi91ukdDCnWRet3DOXA7z37MAQDJZKXVNKVXYfKxIJ8fKWeTLByxfrvuIU4pq3rx6\nDLI9Tz4JPPOMztzVqZNOGiLR1YIg5HCUUvj4449RsqSWga+++gpdunRBdHQ0OnbsiGrVqgEAhgwZ\ngqCgIJw5cwZlypQBAPTs2ROBgYEIDAzEp59+itmzZ6cS4lWrVmHGjBnYtm0b7gSPC7FSKgDaW/2E\nZJx9Em2L6+9pG3I0Srnu2T7yCLB4MdCwoU4acuuWDgQrUMCzNgqCkOu5Ey/WnZQuXdq6XqZMGURF\nRSEqKgo17NIJFyhQAIULF8a5c+esQuxsP3u2bNmCdu3aYf78+bj//vvvyDaPCrFSyg9ahGeSXGgu\nvqCUKk7yvFKqBICL5vJzAErb7V4K2pM9Z15PWW7ZpwyAKKWUL4BCJGOUUucA1LPbpzSANc5sDA8P\nt67Xq1cP9erVc1Yt51C7ts5z3aYN8MsvOqjrt9+Ahx7ytmWCIGRT1q1bh3Xr1nnbjLvizJkzDuuh\noaEIDQ3F6dOnreXx8fGIiYmxes6Wug+Zn49nzpxx2LZ79260aNECP//8M5599tk7Ny6tzuO7XQAo\n6P7cH1KUfwugp3m9F1IHa/kDKAcdYGUJ1toK4AnzMVMGa403r7eFY7DWCehArSDLuhMb0+xYz/H8\n9x/50EM6gCt/frJfP/L338lDh8ikJG9bJwhCNgb3YLDWI488wsjISMbExLB27dr86quvuHr1ahYt\nWpR79uxhYmIiu3btyjp16lj3U0rxueee49WrV3nmzBlWqlSJkyZNIknu27ePISEhjIiIcMmGtO4Z\n9VPYY0L8NACTWVx3m5fGZpFcDeAIgJX2AgmgN3S09CEAz9uV1wSwz7xtlF15HgARAI5C9y+H2W17\n21x+FECHNGx0/T+ZE4mLI9u3Tx1Z7etL3n8/2bgx+emn5Nmz3rZUEIRsxL0mxGFhYRw6dCgrV67M\nwMBAvvXWW0wwjx6ZMGECK1SowODgYDZr1oznzp2z7qeU4ujRo1m+fHkWLlyYX3zxhTVq+u2336aP\njw8DAgKsy8MPP5ymDekJsaS4zMXXD0BL76JFwLp1ekrGgwcBu6YaAED16sC2bWn3Rc+bB/Turccr\nO4kmFAQhZ3GvpbgsV64cpkyZgvr162dqP4PBgGPHjqF8+fJ3bYOkuBTSRimgeXNgxAhg2TKdq/rm\nTT0d4x9/6PmRd+8GRo1Kva/JpJOLvPIKcOwY8MUXOvuXIAiC4DIixEJq8uUDqlTRAj1unC7r29fR\nU46L0xm9Bg3SY5pLl9ZpOC31BUEQ7nFUFg3vlKbpXHz9LvPKK7r5+cUXdR7rhQt1is0zZ4DAQGDu\nXF3v+ed1MpFTp4D8+b1qsiAInuNea5rODsg0iGkgXxoXiYrSw5uuX9ezQu3cqcurV9ciXLGi7mt+\n/HFgxw5g5Eiga1fv2iwIgscQIc480kcs3B2hocDQoXp9507tBY8dC2zfrkUY0H3Nffro9W+/1clC\nrl/X2b3q19cCLQiCIKRCPOJcfP2ZwmQCPv9cp9Ps21fnrnZW53//04Fer7yi81pfvqy35csHzJql\nZ4tyxtixwPjxus+5ZUvPXYcgCHeNeMSZR5qm00C+NB5g9mw9qYSF2rV1INecOdprHjZMR1fbB0FM\nmAB8+KFeVwoYMgTo0UPyYAtCNkWEOPNI07SQdbzyig7aqlFDD3/asAH49VctrqQW2CZNdDmg02x2\n7qzXW7fWdXr1At5+WzdvC4Ig5HDEI87F15/lzJ8PdOgAxMfrz48+qscoG43aU+7RQ88M1b69Hsv8\nyiu2iGxBELIN4hFnHvGIhexBmzZ6aFO/fkBQkA7gsvQ59+ih67Rqpb3lgAA9Z/LSpY7HuHJFe9gW\nMRcEQfACa9euxSOPPIKgoCAEBwejUaNGOHDgwB0dSzziXHz9XuXGDWD6dMDfH3j33dT9wSNGAN26\nAeXK6eCv/Pl1FHadOsC//+p5lZculfHKguAFxCMGLl68iKSkJJQsWRJJSUno06cP1q9fjy1btjit\nLx6xkP0ICAA++gh47z3nQVldu+r5k0+eBAYPBm7f1h71v//q7evX6+jqxMSstVsQhHuOsLAwDB06\nFFWqVEFwcDDeeecd3DLHoEyaNAkVK1ZE4cKF0aJFC0RHR1v3MxgMGD16NCpUqICiRYuiR48elskt\nEBISYp0S0WQywWAwoESJEndmYFqzQeSGBdl0phDBzD//6Nmg/PzIF1/U6yEh5LJlZLFi+nPTpuSt\nW962VBByFbiT2ZdSzvJ2N0smKVu2LKtWrcrIyEheuXKFtWvXZp8+ffjXX3+xSJEi3L17N2/dusUu\nXbqwbt261v2UUqxfv751GsQHHniAkydPtm4/ffo0AwMDaTAYWLVqVcbExGT6npnvjPcF0VuLCPE9\nwHvv2X58BQqQ27fr8n37yMKFdflTT5GbNnnXTkHIRdxrQhwWFsaJEydaPy9dupQVKlRgx44d2bNn\nT2v5jRs36Ofnx9OnT5PUQrxixQrr9nHjxrFBgwapjn/lyhW2b9+ezZs3T+fy0xZiaZoWsjdDhwLF\nigE+Pjp469FHdfnDDwOrVunEIps3A7Vq6UCv//5LfYzz53W2ryVLstZ2QRBsuFOK74DSpUtb18uU\nKYOoqChERUWhTJky1vICBQqgcOHCOHfuXLr7pSQoKAjDhw/HokWLcP0OZqATIRayN8HBwN69wKFD\nQNOmjtuqVweOHNGpNfPn15NRVK2qxXr4cB193bEjULYs0LOnnrTi1VeBixe9cy2CIHiNM2fOOKyH\nhoYiNDQUp+1mlYuPj0dMTIy179fZfvbb7ElKSoLBYECePHkyb1xarnJuWCBN0zmHqCiyc2eyYMHU\n789KkY0a6aZtQDdpz5njbYsF4Z4Fd9I07UXKli3LRx55hJGRkYyJiWHt2rX51VdfcfXq1SxatCj3\n7NnDxMREdu3alXXq1LHup5Tic889Z+0jrlSpEidNmkSSnD9/Pg8fPkyj0ciLFy/y5ZdfZrNmzdK0\nIa17RmmaFnIMJUrofNUXLwILFuhkIOXK6cjsw4eBFSuAffuABg2AmBigbVtg40ZvWy0IQhaglMLr\nr7+ORo0aoUKFCqhYsSL69OmDBg0aYODAgWjTpg1CQ0Nx8uRJzJkzx2HfFi1aoGbNmqhevTpefPFF\ndOzYEQAQFRWFxo0b47777kONGjUQFBSE6dOn35l9zKC9XSnVFcBMklfv6AzZmNw05k0wQ+rJK378\nEahbF1i3Lu2c1hERuh96wAAt9IIgALj3xhGXK1cOU6ZMQf369TO1n8FgwLFjx1C+fPm7tuFuxxEX\nA7BdKRWhlGqslGTiF+5hlALCw3Vmr7//1kLrjN27gXbtgMmTgcce01M+CkI2IjsKnnBnZCjEJL8C\n8ACAqQDeAnBUKfWNUqqCh20TBM9QqJAtpeZXX6WOwkxM1Pmuk5N1sNi5czqj1y+/ZL2tguCEzWc3\n476h96H7yu4iyB4kq/xOl/qISZoAnAdwAYARQBCA35RS33nQNkHwHF266GFRO3boWaLs+eor4MAB\n4IEHgOPHgfff1zNBvfEGMGOGd+wVciwmmmA0GTO1T//1/XHj9g0M3zwcnyz/xKNifDPpJjac3uCx\n42cFJ0+ezHSzNAAYjUa3NEtnRIZCrJT6RCm1E8C3AP4B8DDJDwHUBNDaw/YJgmcoUADo3Vuv9+2r\nJ58AdJ/xDz/occszZwKBgXq+5G++0dtHjfKKuULOwUQT+qzpg3o/10P5keWRd1BeFBteDOeun8t4\nZwD7LuzDiuMrkNc3L/x9/DF622h0XdbVY2Lcb20/1P25LoZsGOKR4wuuecTBAFqTbEQygmQSYPWS\nm3nUOkHwJB98AJQurZOAhIRoD7lJE91U3acP8Pjjup5SwKefAgULAjt3AkePOh7n22+BSpWAefOy\n/hpyKTeTbmLm3plYdTyNPv5szPJjyzF4w2CsP70eJ2NPIsmUhJiEGETsj3Bp/xFbRgAA3q3+Ln5/\n9Xf4+/hjzPYx6LKsi9vFeFf0LozYMgIGZUDDCg3demzBhitCvByANWJaKXWfUuoJACB5Z3M+CUJ2\nIE8ePQ+yUnp6xYsXdf9wrVq6edqefPn0JBOA4xzJMTF6WsfDh/WQqfbtgdhYz9l85oxOYpJLOXPt\nDHqt7oXSP5TGmwvfRJNZTVz2JC14u091xGYtpN2e6obDHx/GtBbTAAALDy/McN/ouGjM+ncWFBQ+\nffJTNK3YFAtfXYg8PnkwdvtYt3rGyaZkvPvnuzDRhE+e+ASPhj7qluMKTkhrgDFtA7T3ADDYffYB\nsDuj/e6FBdl08LmQxcTFkRcvkufP68QgycnO6y1ZohOCVKliKxsyRJdVrEjmy6fXS5Ui332X/OIL\nctAgctEix2MajeTs2WS9enp/+0krbt4ke/fW+bP37XM8/8mTZKFCpMFAjh7ttsu/V1hxbAX9B/oT\n4SDCwQKDCxDh4MD1A1PVPR17molJiQ5lJpOJvVb1YvCwYE7bPc0jNkZei+T1xOtpbt97fq/V9qsJ\nV0mSsQmx9BvgR0N/Ay/FX3Kofz7uPC/HX7Z+/nL1l0Q42Hpua4d6y44us96brku70mQykSSTjEk8\nfPkwbyWnnhhl2dFl7Lq0a6pzWvh247dEOFj2h7KMuxXnsA3pJPSQJe3F2T2jfmpkLMROyv7NaL97\nYREhFjLFrVtkcLD+2ezbR96+rUUX0DNCHT5MPvEEnWbHLVOG/OYbcsECslo1x20PP6xnmvr7b/KB\nB2zlFSqQV/XDmklJZO3ajvt99JEuz2Fsi9zGPdF7HMqirkex6LdFiXCwxewW3Hx2M1ccW2EVCqPJ\naK277OgyGvobWHFURe6K2kWSNJqM7Ly4s1XEEQ5+v+n7NG3Yf3E/+6/rzzOxZxzKjSYjR28dza5L\nuzoIJElO2D6BPv19GDwsmKO3jmaSMfX/5q2FbxHhYJelXRzKG85oSISDP+/+2eGaA4cGMs/APPxo\nyUc8eOkgg4YGEeHgP2f+SXXspUeWWsX41XmvssXsFrxvyH1EOFhhZAUuP7qcJHk7+Ta/WPGF9T40\nndXUKtwWjl85znyD8hHh4LKjy1KdKz1RkeUOtCjDCsDvALoC8APgD+ATAAu9bbhbLl6EWMgsltmg\n+vQh587V6w8+qL1cUgvjn3+SEyaQw4aR3btrQU0pzKVKaW/Zsk0p27bKlcmqVfX6iy/qYw8cqD+X\nKEGOGUP6++vPjRqRV65k6hIux1/m7eTbHrg5aXMr+RYbzWzEdvPbORUoC6uOr6Khv4E+/X04fvt4\nkmSyMZn1p9cnwsEG0xsw2ahbF4wmI8N+DCPCYRWZZGMyq4ytYhUZ/4H+HLllJN/9410iHMwzMA87\nLepk3d57dW8HETKZTBy/fTzzDspLhIPBw4L5x6E/SJJXbl5hs1+bWfctMbwElx9dzmRjMj9f/rmD\nyCMcrDy2MlceW2k9dnRcNP0H+lOFKx6LOeZw3WO3jSXCwVZzWlnL7MXSfnly8pNp3r8lR5Y4tBog\nHFYxRjj4csTLfHLyk0Q46NPfhwHfBBDh4JitY6zHiLsVxzpT6xDhYLv57ZyeR4TYzVqUYQWd0GMu\ngIvmZTaAEG8b7paLFyEWMstff+mfzf33k7Vq6fWxY9Pfx2jUHvOLL5IPPUQOH66boElbU7Svr553\nuV8/MjGRPHGCDArSx2/XjvTx0esrzQ/2f/4hixbVZWXLklu2uGT+n4f+pP9Af4cHflYwfc90qxh8\nsuwTp3VOXj3J4GHBDiLyxYovOGDdACIcLPptUUZdj3LYZ9D6QUQ42GZuG5LklF1TiHAw7McwB8FF\nOJhvUD6uOr6KJDlz70z69PchwsEaE2vws+Wfcd7+eWw5p6W1/gOjH7Cuv/vHuyw/sjwRDgYNDeLj\nkx63brMIv98AP07bPY0LDiyw1kU4+MaCN3g5/jL7rulLhIMt57RMde1nr50lwsH8g/Pz5u2bjLkZ\nY216/2XvL2z7W1uqcEWEg/MPzE/3Xq8+vpod/+jIKbum8HTsad5Ovs1hG4cx/+D8VptKjSjFjac3\n8rf9vxHhYN5Bebn/4n6eu36O1SdUJ8LBkO9CePHGRafnECF2sxZ52wCvXrwIsZBZkpPJYsVo9V4L\nFdJ9zHfLyZN6sWfZMkdP+bPPUu/z2GN6m6+vFnijkWmx9uRa5hmYhwgHVbji6djTd2+3C5hMJlab\nUM1BFCftnORQJ/52vLVOk1+acMquKfQd4Ouwj7Mm0nPXz9Gnvw99B/jy5NWTLPl9Sat4keRv+39j\noSGFGPBNANedXOew76LDi6weYUoPcs6+OTSajPx+0/f0G+Bn3VZzYk2evHqSycZkDtkwxLotaGgQ\n155caz12YlIih24YavWsQ74LsTYrrz+13ul9qjmxJhEO/nnoT4avDSfCwednPm/dfujSIa44tsLp\nvq5w6uoptpvfju0XtHdoVn974dvWF4rSI0pbm7KPXD6S5rFEiN2sRRlWAEIAfAVgEoBp5mWqtw13\ny8WLEAt3wscf0yqO3bp59lwDBujzPPIImZCQevutW+Snn1rtOdiwutNgs+3ntltFp9CQQkQ4OGDd\nADI+nryug4uSjcnccnaLQ3+rO1h7cq1VjMZtG2f1Hv8+9TeTjEk8d/0c281vR4SD94+63xrEtPr4\namuzao+VPdI8fovZLYhw8MHRDxLhYPUJ1R2uITYhNs2ApOuJ17ny2Ep+veZrPjfjOb4U8RJPXj3p\nUGdb5DY+NfkpdlnahQlJjv+DXVG7+Pnyz3n48mGnxz9y+QjrTqvrIOQp+2MtWDz/V+e9mqFou5Pr\niddZYWQFq41PTX4qTU/Yggixm7UowwrAZgDDALwC4CXz0sbbhrvl4kWIhTth40b90zEYUnux7sZk\nItevJ69dS7fa7MGv8Zo/SIB7X3jUwTPeEbmdD38ZyKfeAX/8vDYP9HqX0/4HHi7uR5PBoL36vXvZ\nfWV3Ihz8bPln6ZwpbYwmI09ePcnouGiHcku/avjacJLkZ8s/I8JB3wG+NPQ3OERB77vgGCl+LOYY\nZ/07K91+5cWHFzt4tPb9stkBo8nIMVvH8LGfHuPfp/5Os54lotqy1J5SO03RdjdbI7ey2HfF2G5+\nO968fTPD+iLE7l1cmX1pD8lq6Va6R7HMFDJjhp5TfuJEoGhRb1slZHtIPc64ZEk9zaKHMdEEg0p7\nyP9t422E/RiGCvujsWImkD8Z2P1KXVT+ZQXmDn8LFSbMRe0zae4OAEgqXRJlX7+A6HzJAIAV7Veg\nUd4q+lpLlcK+C/sw78A8HLp8CIdjDuN07GkUzFMQRfIXQeF8hRGTEIPDlw8jITkBCgrfNfwOnz/1\nOY5dOYYHxzwIfx9/nPnsDEIKhCDZlIyWc1piydElAICQAiEoW6gsBj47EM/f/3ym74/RZETYyDBE\nXo/Ec+Wfw6o37r0kH4B2isqPKo9TsacAAEtfX4omFZtk6fldza2c3kxCQuZxRYgHAdhMcknWmJR1\nWIT4+eeBlSt1RsP27b1tlSBoSGLc9nHoubonqharim5PdUOrSq3gY/BxqDfnvzl4bf5rqFK0Cr4z\nPYf6n45EHiMQGeSDUld16s7EfH7wr/IIDKVKAyVL4k/DUQxOWIlH672GsSOPAdu34+8yQIePQnE5\nNgoj/861bLXtAAAgAElEQVSPtzclQJGIL1EEi4pewe8PmBBRBUA6j99iBYrhQvwFAMAHNT8AAEzc\nOREdq3fE5OaTrfWMJiMuxF9A0fxF4efj53iQa9d0itGAgNQnuH0b8PNLNXXl9D3T8e2mbxHxUgSq\nhFRx8Q5nPz5b/hl+3PojqhWvhl3v78qySQcyiwixm8nIZQZwA4AJQCKAOPNy3duuvDsWmJumf/yR\nBMi2bSkI2YJL8ZcchspYlvIjy6dKRlFrSi0iHNbhPou/fZdJCiTAiwV9eOLLD639wBaOXD5ibQ5e\ns/EXRhbU9eMbP8fo4DwkwCQfxYSAvLQfdrW5z1vcGrmVl+Iv8XTsae6M2snlR5dz89nN1r7dOfvm\nWIPCLMu+8//q4LM9juODUxETQxYvrqPSb6ZoIt24kcyTh+zV667ubXbmTOwZtp7bmtsit3nblHSB\nNE27dcnQI87JWDziY8eAihX1FLUXLwK+vt62TMiJbDyzEb8d+A296/RGSIEQp3VI4s/Df+LDJR8i\n+kY0AvMGYmzTsYhNjMWIzSNw/OpxAMDwhsPRrVY37IrehZo/1UShPIUQ+XkkAvy1F7lm2tdIPH4Y\nz/YYj3z3BTs919NTn8Y/Z/9BAb8CeOhUPDZN94Xfbd00vaukAW83M2FfCFDlEjD+dkM8PWUV4O8P\nbNoE1KyZ7rVuOrsJLea0wOWbl/G+75OYuMIP2LBBT7bx339AWJjzHXv21Lm7AT35xqefWm4M8NRT\nwNatgMEA7NkDVK2arg2C5xCP2M1kpNTQ+ajfAPC1+XMZAI97+w3CHQvsgrUqViQB/dItCO7mwo0L\n1jGyD415KNV4WFJH59pH2D499WmHIUbJxmRr4geEgz/t+MmaqelOAqwm75xsPVbx4cV5c+4vepzz\nt99yxo6p1gQYc/bN0Tt8+KH+kZQr51ISkZOR+/lPq0dpsoyBtgzFathQB6GlJDKSzGvngYeE2IaG\n/fEH7T1z1q/v/BhClgDxiN2rRRlWACYAGAfgkPlzMIAd3jbcLRdvJ8SffKLvRu/eFAQHTCYTt0Zu\nZb+1/VJlRDKZTBy4fiCfnPwkd5zbkeYxXpn3ChEOa5RwxVEVefbaWZpMJq47uY4vRbxkFcXgYcH8\ncfOPaUYKj9462joW2G+An9NMTa5wPfG6NcnD5J2TU13XkiNL+N+F/2yFiYlkzZr6h9K8OXn2rE7r\nuWdP6mbkpCSyaVNao8u7dNF1CxfWZVOmpDbogw/0tjZtbKlChwzREeCWTGNff21LdLJgQaavWXAP\nIsRu1qIMK5gneIDdRA8A9rp0cGAqgAsA9tmVhQOIBLDbvDSx2/YlgKMADgFoZFdeE8A+87aRduV5\noLN+HQWwBUBZu20dABwxL2+mYZ/1i7Vihb4b1apREEjqJBOTdk5ijYk1rCIZ8l2IQx7kr9d87ZC5\nad7+eamO8/vB3639sTvO7bAmrgj7MYyVxlSy7u8/0J/dV3a39rWmhyWjFMLBF2a9cMfXOGffHH69\n5mtr2sgMOXGCDAxkqpSdoaF6mBWpPdWOHXV54cLkzp22/WfNojURSmSkrfzIEZ09zGAgDx4kV63S\n9YKCyPHj9Xrp0vplYMwYWj1zZ2OrybTLBbcgQpz1QrwVdjMuASgKF2dfAlAHQPUUQtwPwOdO6laG\nnunJD0AYgGOwRXVvszSHA1gKoLF5vTOAceb1VwHMMa8HAzgOINC8HAcQ6OSc1i9WYiKZP7++I/bP\nByH3kZCUwJFbRrLYd8UcvFSLgAYODeTms5s5dMNQq5f7/MznrXUHrBtgTShxNeEqSwwvQYSDo7aM\nIknG3Izhoz89aq1fYngJfr3ma569dtZlG00mE/uu6ctCQwpx05lNHrkPabJ0KVmypBbf++/XAmnx\nfAcPJvv21Z/z5SM3b05puE71CZBNmpD//ae96bZtddk779jq1a1LhyZtixedlKQnygC0p/3XX+TW\nreTq1WTPnuT//qe3PfaYy6k/SepAsbOu/w8yRVwcOWkSeeGCZ46fxYgQZ70QtwfwJ4BzAL4xe5iv\nuHwCLaophbibk3pfAuhp93k5gCcBlABw0K68LYAJdnWeMK/7ArhkXn8NwHi7fSYAaOvknLSneXN9\nRyY5Zt8Tcgkmk4lTd01lqRGlHDIhzdgzgwlJCUxMSrTmIrZPFTlz70yaTCZ+v+l7az7gAoML8LGf\nHrOmLaw1pZaD1xmbEMv+6/pzwYEFdzUBQ1YlfEiXpCTyyy/p4CEbDLpf1xmRkeR99znWB3RE9Gm7\ntJvr19u2PfCA40xTq1en3j+t5a239BSX6dk/cqS2yWDQLwUpp6C8W9q107Y89BB5yS7L161buvWg\nZEny1VfJn37SrQ5ZgcmkPZA7QIQ4i4VY33M8BOBj8/JQpk7gXIhPAdgLYIrFUwUwGkA7u3qTAbQx\nN0uvsiuvA2CReX0fgFC7bccAFAbQDcBXduV90hB/2jNhgr4jrbI2H75AMjZWP7fv8Llw15hMJofZ\nbqqOq8qFBxemErokYxLbL2hvrTdxx0SH7YsOL2LZH8o6DN3xH+jPAxcPuMXOW7fIGzfccij3s2yZ\nrQ944sT0665frz3iihV1nuy0AjQs/czzUjf5c/hw8oUX9LzOjz6qPeBu3XSz9qVL+uXAMktV8eLO\nvdFNmxynpTQYbOutWjm+GJBavH75RXv+rk5BuWEDHV4MHn1UDye7fl3PnuXs5aF/f9eO7YwbN3RL\ngzMOHdIzeb3wAlmkiL7eF1/UM4ZlYkpNEeIsFmLoKOkyAMqalzIAyrh8gtRCHAKdEkABGARgirnc\nK0Lcr18/6zJ37loCZECA41ztOQGTSacV9gRHjuhJg5w5ESaTbUrdtLh+3fYsrFUrfeflaMxR9l3T\nN9V8tdFx0ey6tCvbL2jPa4mp00EmJOgUzCaTiTE3Yxxy6RpNRutMPb4DfDl552QaTUbevEkeP57a\nBkvKwgUH0g4Wuhx/metPref47ePTTWvojIsXdRxSyudifLyeIdHHh3zqKd0CnJmW1yzh8uWMxwqn\nJCmJjI4mTSYajTo4u4cltXRcHLl9O0kd6/XII5lssTp6lHz8cf3lSpkoYPFim/CWLavfBM+c0bnE\n8+ix1AwO1i8YpP4HvPmmTSzffz/jyO3kZNuXu1Mnsnx5vW55eQD0LFqLFulZvFq1ss20NWtWJi7U\nzN69tnNMneq47cgRW/+bs6VUKT1ntpMf7Nq1ax2elSLEWS/E/5kFzxIslQxgv8snSCHEaW0D0AtA\nL7ttywE8AaB4iqZpa7Ozuc6T5nX7pmlr87X580QArzo5P1NiCc5cvTrVpnsCZy+1587pINSAAHLt\nWved67//yNdesz3LgoLIf/+1bY+P1y/bSukH620nLbBJSWTjxqmfBzucBCDvu7CPId+GWD3NZr82\n4/pT6xm+Ntw6ZZwlR2/cLduMSL8uiqLfyx1o+Phh+va1zbbz0JiH+NGSj6wRy3kG5uHiw4tJ6m7L\nRx7R9kyenNoWT3Hlim6FBVJPttStG50+P7/7Luvs8zS//Wa7rjVrHLe1akVrC/ahQ5k46IkTNgFa\ntEiXHT2qA8YALbwp31KjomzeuFLkF1/Y+p7z57cNsxo0yLZPUpI2ev9+W5mlma10aRrj4vWbXYkS\ntossV46x249w2zY7TR81ynahW7emfV1xcY4TfMye7Si0AQG2Zu7kZP32BpANGui6J07ot77vvrON\n37Ts162bfqE6eVLXuXnT4aVDhDiLhTjVDkANixfrYv2UHnEJu/XPAPxqXrcEa/kDKAcdYGUJ1tpq\nFmWF1MFaFlFuC8dgrRPQgVpBlnUntjElPXvqu1KkiP4d9u5NLl/udEKbbMewYbqVr2VLW8vUrl26\n+8nyG7vvPv3S7I5zWWJo/Py0pwboGQKPHNEv1bVr284LkE8+SZ46ZTuGyUS+957tfm/caJviN29e\nLYCW3/7e83sZPLSIFtEPHya+ypcq61TtMc1Y6ns9jVu9n+sx/nY8P5uwkOhR2LHulwFE7/wOZQUG\nF+CaE7Ynf+fONrsNBjIi4u7vmT1RUfq79tdftrLbt/Uz0v6eWbpZt27VdhgM+lm/aJHNRoNBR/27\nSmys610ARqPWlalTyRkz7m7o7o4d5LF0RlkZjWSVKrZrr1XLdr6dOx3vS+3aac/4aDKR69alGOo8\nYoTesWRJffMtwV6tWqV9UUajFlr75uqKFXXTz4IFth/A+PHk6NE6ittSr2lT7XGbm+r/+TSChQrp\nybRM/+7TdjzxBM/viWaFCnqXhx4ix40j466bbEO5ihd3HkD222/ac/bz029uderYzv3mm7a3lrp1\n9XUMG6Y/h4Y6HwNuMukvUcovoP0yapS1ugixl4VY/w/wn4v1ZgOIAnAbwFkA7wCYAeBf6D7ihQCK\n2dXvbW5ePgTgebtyy/ClYwBG2ZXnARAB2/ClMLttb5vLjwLokIZ9TMn+/bb51u2X0FDd5XToUPoP\no6Qk3UUTF5d+PZNJi9LKlfo3PWuWFp41a+6sCXnQIEd7lSJbt7a9ID/9tO23WaKEoyDak5ysH5Yb\nN+rf+qRJjrElpH4O+frqc3TurLvREhLI556j5eXf6k2WKkXOnKn/AnrkS7t2Oj6ldWtaRXeTOfA3\nMdE28gXQuR/+2L6DhQabJ4xv14TBIQkMLHmeeK4n8UUx4p1aRNn1BEifokeZ58tQHd0cHmYV2rJ9\nGvL3bVv5408xrFPXRPjcYrXmGzlg3UC+HPEyt0baPI+FC2l9wXjnHdv6stTT4ZLULQ4REfq+2f/P\n4+J0V+jevY7lK1fqXBWWa3znHf1s7NRJfw4JsXm/QUHaebPoxhdfOJ77669t9Zw1o5P63Fu26LqW\n6Yt9fMhKlciXXtKtkTt32oTt2jVy+nTdhZsypmr8eOfnSI+EBB3cDJAFCqQOpLYwZ47tO1OkiF5f\nskRvswRad+xocyjtdMFKZKSt27V0abuXzuRkWxO1pR/7wQcznNWKpG4eq1BBB1PFxtrKR49O/aAo\nWzZV829i7WcZWMhkLerblzTdTuLlSybri4elNRrQjvqiBbd18zWgPfGYGNt5Dx8mCxZMfW5fX3L0\naP46y8QP2lxiclHzvNkffGDtK9/1zTJ26qT1tnRpfa727fUlWl9stm/XzfiVKulKwcF6f7ux3yLE\nWSzE5v5Wy9LdLK4rvG24Wy7eiRCT+gt59Kh+uPbsSesbq2UpXlx7nQMH6ubDRo20UNv/mAD9W6lS\nRTe9tm6t8xS89JJ+ebW0ijlb/Px0K9Inn+hgzj//1L+NxYv1G3Pv3voZsHevtnXgQFrFd+RILY6W\n+BfLC3Jion4gWn7bDz6or2/NGt0CNWuWFkjLM8p+qVbNluDIZLK9fHfq5HjfbtywebSWc1hiXS5f\nJps1S31spbTg22Myae+rUPlDRKs3iK99tKC+1oxPPZ3ICxd0nehoLWojRujA2Bo1zP+DIge1QIeD\n6JOHzb/5wWF+2pgY23XOmeN47shI/dwByO+/1+f5/HP9OV++1H2yJpOtq8/yktO8ue7isHekKlfW\ngterl82RqlHD1hVp+T7kyaOFymjU8TT22ypUSP2SZjTaRKpq1dSBXIcP2/7nlsXf39E2y1KsmP4u\n53VML81SpfQLkeUeHDzo9GfjlEOHbC26liUwkNy927FecrJ+7gM6zmv4cL1evbq+5xYRv3CB/P13\n22f7WSjnzLHl+rDc44AA/bshqd8g/fxsGw7YAugSE8mPPtLOZeXK2ub69Z23DJ8+TW6zpIK2RIs/\n8YT+Iicn6zdXS+KRAgX4WcN9BPSzwPKM6N7d9lL00EP6uzx3rq31ODCQjNx72dZkXL26/uLevGnr\nP3v5Zf0P37dPN50cPmx1/AGy/xOLHW58dIv3nf7fLUvp0vr7mWbAuDRNe2zJuIJOwNHPvHwFoB2A\nvN423C0Xn4YQp8RkIv/+m+zQwfZDT2sxGPTDKuXDzNlStCj5zDNkixa6r/XNN/XDOb0fS8rF4rEo\npcXLwrFjWpBHj3b0xmJjbd5qWkuJEvq50rKlfsEHtLgkJ5M//6w/h4Q4b+G6elU/zBs21F1LKe/j\n6tXaQ540Sdu2yckQ2F1Ru9j2t7bWoUDo60u80Ikd3rmVYZNqTIw+/vPtDrLQK59xyNR/ndabOJFW\nkbGI19WrtqGrzz9v8xBMJptnXKOGY5Po8uW0vnRZvDjL4uur66csNxh0E2VyshY1+yZ8+/icS5cc\nuxVS9plaiI219SuXLatf4Fav1q0kFqEvXFh/HxYv1mKekKDFcMYM3T1gabGwLHXq6O5N+zH17dvb\n7oErwYyLFmmxtLxEbN6sv1OW7769oFvyfISF6WPfvGnzfC33wH6uh1desX0PS5VydEKbNNEC/frr\ntvs9cqR5x5Ej9ckXLrQe69o1LbrOfgsFC+qgZwuLF2sNB/RLenIy9Vums+avmze5ePplq+6fPq3F\n1v6FvVw5x3tsP8z62WdJ45lIPVbbIsaWi6pY0cGbN5l0wKTluPny6b/H6uu+H2PZMFYNu05A/x+X\nLNHOxtGjZHi4Y6s6oJ8RX36pvwNLlujYjzhb2IUIsZsXrxvg1Yt3UYjtMZm0lzFjhvaGBw/WL6PH\njzsGI5lMWhR279YPpHnztAcaEaG7YsxBok6JjdU5E4YN017n88/rN/Tnn9cPzX799I/JPo/CzJmu\nX8P589rLa91aC0/lyrpZecSI1E3vhw/bXj7ef9/WbG8v+nfK9nPbOWPPDC49spTbIrdxyZElbDC9\ngbU52W+AH9//833+suQE5893b2rh5GQtKIBuYYiI0C0dlod7ysjtmzd1qwegH6YWLMI9bJi27+BB\n8tdfyX/+sWV9vH1bP8DbttUZIlMGzBmNekTM/Pmp7dywQYtBysCtlBw4QJYpQ6di8tZbWivSw2TS\nntAvv6TdbREbq4XSIoonT2pde/FF/X2yD7adN8/WKvPaazbNSEy0NR0XLapfPkeMsDl+9pkvx461\nXUPBgo7XcP68Y/O+xYscP972PTGZ9Cggy/alS807271JRUXZPPbixXWf/b59+ndryTFSoIB+CRox\nIvVLctOm+r4kJGjB+uQT/X1askTfH8vLxOjRNtsjIrQYlyjhvDvh/Hnb7+z776mV2iLGgH7L37uX\nsbH6ezZpks1Wg0G/LE+apD+XLhzPG18OYr+X91sF1tnLrNGou1Hefz9tZ+P77231RYjdrEUZVgAW\nQSf0WORs3dsXcFcXfwdCnN04eVIHR3mStWsdm7qfeebuRPF83HmHsbgpl4BvAvjZ8s94JvaMuy7B\nKZs2MdXDpnbttJteLV60JbeEZXhoYKBrXY13iqv3OjlZP5h79dJ9yjVruj/6f+PGtFtsihfXIjNz\npq1O9+6p7Y+P19+hlPtXqOD4Mnvrlq1Fpm/f1LZcuqSDEU+fTj8mY8gQfYxixRyHEh85YnuxqFgx\ndR6N5GTbaCV7L7Z/fy3Ylu6N0FCb5+9sqVUrdWDZ6dOO3c0p+fNPva+/v3kkgr0YT5nCWbNsrR2W\nxc/P1s1jMtniriyjp6zHyoDERO1c9OunW4IaNtTdBnaNCCLE7taiDCsAo6DzOTcD0NzcR/wjgGcA\nPOPtC7iri88BQpxVTJ2qvy2+vo6jMzKDyWTiuG3jWGhIIetwoTZz27DRzEasPqE6q46ryiEbhriU\na9lddOigr6tQId0Ml1YkLqlFwvIsnDTJNrrFmUjkZCxNoAEBupvyp58cYwMsS79+aYtjcrJ+ERo/\nXrf6NGzoGEFuYcsWPfTtbpKYJCfbhL9ZM23Tzp02r/Oxx1J3o1gwGsl336XVEbWPKThxwrGbp3p1\n/V3o0UO/0Pn7a0/+Tn8v779P6wvKgQPUbxv79nHRItuLQbVqunXsm29SD98+ccLxBcGdw9xEiN2s\nRRlWAHa6UnYvLiLEmWPpUsf+Mgu3k29z7n9zeerqqTT3TTIm8Z2F71i93sa/NL6jGYPcTXy8bsqL\nSj0roVNmz9a/GktAV/78qaPKczomkxYG+yZOo1EHElpiFoYM8Z59zjh92hb09uGHtqDjRo0c+z6d\nYTTq7ghnyapu3NCjHlIm4CJ1c3VGx06PGzds3mxAgG5tWL/eFn9i32eeFpb5MZ55xr1DMEWI3axF\nGVYADgKoYPe5POwSbNzLiwixe/hoyUfWsbijt452iFAm9SQKlhzN+Qbl45x9c7JHjuQ7wGh0zIiY\nUd9tbuPiRd1cnB2xvERZlrZts38GvRs3dB+7ffcw4FpSL1LX2bjx7l4InCFC7N7FkjAjTZRSjQH8\nBOCkuSgMwPskV6S74z2AUooZXb+QPrP+nYX2v7eHQRlgogkA8HSZp9GnTh/4+/jDRBMGbxiMtafW\nIjBvIJa8vgS1StfystV3x7JlQNOmgJ8fcPIkULKkty0SXOWdd4Bp04CPPwZGjgQMBm9blDEkMHo0\n0K0bkJwMvPwyMHs24OPjPZuUUiCpvGdBziJDIQYApVReAA+aPx4iecujVmURIsR3x78X/sWTk59E\nQnICJrwwASEFQvDhkg9xIf5CqrolAkpgRfsVqFqsqhcsdS8kMHYsUKIE0KaNt60RMgMJnD0LlCnj\nbUsyz/btwNatwHvvAXnyeNcWEWL34opHXADA59ATPbynlKoI4EGSi7PCQE8iQuwasYmxmLhjIsZu\nH4vE5ETUL1cf9cvVx3ebvsOxK8fQ4X8dMK3FNCilcCXhCvqt7Yc9F/bAR/nAoAwIKRCCbxp8g/JB\n5b19KYIguAERYvfiihBHANgJ4E2SVczCvInk/7LCQE8iQpw+1xKvYeDfAzFx50TcuH3DaZ3/Ffsf\nNnXchPx++bPYOkEQvIUIsXvxdaFOBZKvKKXaAgDJeKXk/ud0dkfvxsvzXsbxq8cBAA3KNUD3Wt0R\nFhiGv07+hdUnVuNi/EXMaDVDRFgQBOEucEWIbyml8lk+KKUqAMgRfcS5mdOxpwEAoQVD4efjZy0n\niZ92/oRPln+CW8ZbqFa8GiY3m4yaoTWtdR4s8iA6P9Y5y20WBEHIibgixOHQ8/6WUkr9CqA2gLc8\naJPgYfqs6YPBGwYDABQUigUUQ17fvEhISkBCcgKu37oOAPig5gf4sfGPyOub15vmCoIg5GjS7SNW\nShkAvAzgLwBPmou3kryUBbZ5nNzYRzx221h8vOxj+CgfFA8ojugb0dZhRxYC8wZiTJMxaPdIOy9Z\nKQhCdkb6iN2LK8FaO0nWTLfSPUpuE+LfD/6ONhFtQBDTWkzDW9XeQrIpGedvnEeSMQl5ffMin18+\nBPgHwNfgSmOJIAi5ERFi9+KKEA8FcBk633S8pZzkFc+a5nlyqhDH3YrDn4f/xIJDC3Az6SaKBxRH\nUN4gjNs+DreMtzDo2UH4qu5X3jZTEIR7FBFi9+KKEJ8CkLISSd7zg0JzmhCfuHoCvVb3wqIji5CY\nnOi0TqeanTDuhXGQyHdBEO4UEWL3kqYQK6VeJjlPKVWe5IkstitLyElC/N/F/9BoZiNE34gGANQp\nUwevVnkVZQPL4sKNCzh/4zyKBxTHW9Xego/Bi7nxBEG45xEhdi/pCfFuktWVUrtI1shiu7KEnCLE\n289tR+NZjXEl4QrqhdXDjJYzULpQaW+bJQhCDkWE2L2kJ8SroZukHwOwIcVmkmzuYds8zr0uxLGJ\nsVh8ZDE6L+mMuNtxePGBFxHxUgTy+eXLeGdBEIQ7RITYvaQXGtsUQA0AMwEMB2B/0+9d9boHiY6L\nxtz9c3H91nUkJici/nY8tp7biu1R261Dj9o+3BYzWs5wSM4hCIIgZH9cCdYKIXkxi+zJUrK7R5xk\nTMKYbWPQb10/xN2OS7Xdz+CHp0o/hdaVWuPjxz+Wvl9BELIE8Yjdi0vTIOZUsrMQbzyzEZ0Wd8L+\nS/sBAE3ub4KaJWoir29e5PXNi0pFKuGZsGcQ4B/gZUsFQchtiBC7F8nakM0wmoz4ZsM3CF8fDhNN\nKB9UHqMaj8ILD7zgbdMEQRAEDyBCnI24cOMC2v/eHqtPrIaCQq/avdCvXj/J9SwIgpCDyVCIlVIP\nAhgHoLh5PuJHADQnOcjj1uUiNp/djNYRrXH+xnkUzV8Uv7T+BY0qNPK2WYIgCIKHMbhQZxKA3gBu\nmz/vA/CaxyzKhczYOwP1ptfD+RvnUbdsXezptEdEWBAEIZfgihDnJ7nV8sEc3ZTkOZNyD0aTET1X\n9USHhR1w23gbHz32EVa/sRqhBUO9bZogCIKQRbjSR3xJKXW/5YNS6iUA0Z4zKedzJeEKpu2ehnE7\nxuHE1RPwNfhidJPR6PRoJ2+bJgiCIGQxrowjrgDgJwBPAYgFcBJAO5KnPG6dh8nq4Usk0WdNH/yw\n5QckJCcAAMICwzC1+VQ8W+7ZLLNDEAThbpDhS+7FFY/YRLKBUioAgIHkdaVUOU8bltMgia7LumLM\n9jEAgEYVGuHjxz5G04pNJRGHIAhCLsYVj3g3yeopynaSrOlRy7KArPKISaLbym74YcsP8Pfxxx9t\n/0Dj+xt7/LyCIAieQDxi95KmR6yUeghAZQCFlFKtoXNNE8B9AGRgq4vcuH0D/df1xw9bfoCfwQ/z\nX5kvIiwIgiBYSa9p+kEAzQAUMv+1EAfgPU8ada8TdysOP+/5GYuPLsa6U+tw23gbvgZfRLwcgRcf\neNHb5gmCIAjZCFeapp8iuTmL7MlSPNE0ff3WdTSY0QA7onboc0DhyVJPok/dPmhasalbzyUIguAN\npGnavaQ3H3FPksOUUqOdbCbJrp41zfO4W4gTkxPRZFYTrDu1DuUCy6F/vf5ofH9jFC1Q1G3nEARB\n8DYixO4lvabpA+a/O2Gbf9hy47PnlEVeJMmYhFd/exXrTq1DiYAS+OvNv1AuSILLBUEQhPRJM7MW\nyUXmvz+TnE5yOoCZAH43r2eIUmqqUuqCUmqfXVmwUmqVUuqIUmqlUirQbtuXSqmjSqlDSqlGduU1\nlY1zD1UAABMESURBVFL7zNtG2pXnUUrNNZdvUUqVtdvWwXyOI0qpN129IXcCSby/+H38efhPBOcL\nxqo3VokIC4IgCC6RYYpLpdSvSqn7lFIFoPNMH1BK9XDx+NMApAwR7gVgFckHAPxl/gylVGUAr0JH\najcGME4pZfHAxwPoSLIigIpKKcsxOwKIMZf/AGCY+VjBAL4G8Lh56Wcv+O4mYn8Eft7zMwr4FcCy\ndstQJaSKp04lCIIg5DBcyTVdheR1AC0BLAMQBuANVw5OcgOAqymKmwOweNTTzccFgBYAZpNMMmft\nOgbgCaVUCQAFSW4z15tht4/9seYDaGBefx7ASpKxJGMBrELqFwK3cOHGBXy09CMAwIjnR+Dxko97\n4jSCIAhCDsUVIfZVSvlBi98ikkm4uz7iYiQvmNcvAChmXg8FEGlXLxJASSfl58zlMP89CwAkkwFc\nU0oVTudYboUkOi/tjJiEGDxX/jm8V0NGdQmCIAiZwxUhngjgFIAAAH8rpcIAXHPHyc0hy/ds4FfE\n/ggsOLgABf0LYnKzybC1pAuCIAiCa2SYa5rkKACjLJ+VUqcB3M0MBReUUsVJnjc3O180l58DUNqu\nXiloT/aceT1luWWfMgCilFK+AAqRjFFKnQNQz26f0gDWODMmPDzcul6vXj3Uq1fPWbVURF6PtDZJ\nD280HGUDy2awhyAIwr3JunXrsG7dOm+bkWPJMKHHXZ9Ae9CLSFY1f/4WOsBqmFKqF4BAkr3MwVq/\nQgdXlQSwGsD9JKmU2gqgK4BtAJYAGEVyuVKqM4CqJD9USrUF0JJkW3Ow1g4ANaCHXO0EUMPcX2xv\n2x2NIz5/4zye+fkZHIk5goblG2JF+xXiDQuCkGuQccTuxZXZl+4YpdRsAM8AKKKUOgsdyTwUQIRS\nqiN0k/crAEDygFIqAnr8cjKAznYq2RnAzwDyAVhKcrm5fAqAmUqpowBiALQ1H+uKUmoggO3mev1T\nivCdcin+EhrMaIAjMUdQrXg1zHlpjoiwIAiCcMd43CPOzmTWI76ScAX1p9fH3gt7UaVoFaztsFay\nZgmCkOsQj9i9pDf7UhvoQCrLrEsOkFzgQbuyJW8tfAt7L+zFA4UfwOo3V4sIC4IgCHdNek3TzaAF\nOARALdiCnZ4FsAlArhLixUcWY9GRRbgvz31Y/cZqFA8o7m2TBEEQhBxAmkJM8i0AUEqtAlCZZLT5\ncwnYkmjkChKTE/Hp8k8BAOHPhKN0odIZ7CEIgiAIruHKOOLSAM7bfb4APWQo1/D9pu9x/OpxVC5a\nGR8//rG3zREEQRByEK5ETa8GsEIp9St0f/Gr0CkjcwVnrp3B4A2DAQCjm4yGn4+fly0SBEEQchKu\nCHEXAK0A1DF/nkjyd8+ZlL34YuUXSEhOwMuVX0b9cvW9bY4gCIKQw3AlsxaVUrsAxJFcpZTKr5Qq\nSDIuC+zzKgsPLcS8A/OQ3y8/hjca7m1zBEEQhByIK9Mgvg9gHoAJ5qJSABZ60qjswOWbl/HB4g8A\nAEMaDEGZQrmqW1wQBEHIIlwJ1voIwNMArgMAySPQQ5pyNF2WdcHF+IuoW7auBGgJgiAIHsMVIb5F\n8pblg3lyhRydjmv+gfmY898cFPArgGktpsGgXLlNgiAIgpB5XFGY9UqprwDkV0o1hG6mXuRZs7zH\n5ZuX8eGSDwEA3zb8FuWDynvZIkEQBCEn44oQ9wJwCcA+AB8AWAqgjyeN8iYD1g/ApZuX8GzYs+j0\naCdvmyMIgiDkcGTSB7vrPxV7Cg+MfgDJpmT8++G/eDjkYS9aJwiCkD2RSR/cS3qTPuxLZz+SfMQD\n9niVfuv6IcmUhPaPtBcRFgRBELKEND1ipVRYejuSPOV+c7IWe494/8X9qDq+KnwMPjj88WHpGxYE\nQUgD8YjdS3qTPpzKQju8Tp+1fUAQH9T8QERYEARByDJcSejxlFJqu1IqXimVpJQyKaWuZ4VxWcWW\nyC1YeGgh8vvlR5+6OTYOTRAEQciGuBI1PQbA6wCOAMgLoCOAcZ40Kishid5/9QYAfPLEJzLPsCAI\ngpCluJSpguRRAD4kjSSnAWjsWbOyltervo6qIVXRvVZ3b5siCIIg5DIyHL6klPobQEMAkwFEQ89N\n3IHk/zxvnmexD9YiCaUk9kAQBCEjJFjLvbjiEb9prvcxgJvQkz608aRR3kBEWBAEQfAGrnjEAQAS\nSBrNn30A5CF5Mwvs8ygpE3oIgiAIGSMesXtxxSP+C0A+u8/5Aaz2jDmCIAiCkLtwRYjzkLxh+UAy\nDlqMBUEQBEG4S1wR4nilVE3LB6XUowASPGeSIAiCIOQe0sysZcenACKUUtHmzyUAvOo5kwRBEAQh\n9+DS7EtKKX8ADwIggMMkkzxtWFYgwVqCIAiZR4K13IsrKS5fAZCX5D4ArQDMVUrV8LhlgiAIgpAL\ncKWPuC/J60qppwE0ADAVwATPmiUIgiAIuQNXhNho/vsigEkkFwPw85xJgiAIgpB7cEWIzymlfoIO\n0FqilMrr4n6CIAiCIGSAK5m1CkBP8vAvyaNKqRIAqpJcmRUGehIJ1hIEQcg8EqzlXlyKms6piBAL\ngiBkHhFi9yJNzIIgCILgRUSIBUEQBMGLiBALgiAIghcRIRYEQRAELyJCLAiCIAhexGtCrJQ6pZT6\nVym1Wym1zVwWrJRapZQ6opRaqZQKtKv/pVLqqFLqkFKqkV15TaXUPvO2kXbleZRSc83lW5RSZbP2\nCgVBEAQhY7zpERNAPZLVST5uLusFYBXJBwD8Zf4MpVRl6IQilaHHNI9TSllC58cD6EiyIoCKSqnG\n5vKOAGLM5T8AGJYVFyUIgiAImcHbTdMpx6E1BzDdvD4dQEvzegsAs0kmkTwF4BiAJ8zJRQqS3Gau\nN8NuH/tjzYfOky0IgiAI2Qpve8SrlVI7lFLvmcuKkbxgXr+A/7d37zF2lGUcx78/KCCgiNWmIJe0\nEpBiAMFShFpYwdRCQkFagQqCgqg0iH+gEYFIDcZIjAkgAS+USwBDiJRSsNxh5RJoKeXaUqRKgSJ3\noeUSoYXHP+bd7rCe7c5Cz7y7s79P0vQ9c94z88yTOec5857Zd2Bkan8WWF567XJgqxbLn0vLSf8/\nCxARq4EVkoav870wMzP7CIZl3Pb4iHhe0gjgFklLyk9GREhq+7RXM2bMWNPu6Oigo6Oj3Zs0MxtU\nOjs76ezszB1GYw2IKS4lnQG8CRxP8bvxC2nY+Y6I2FHSKQAR8ZvU/0bgDODp1GdMWj4N2CciTkh9\nZkTEfZKGAc9HxIge2/UUl2Zm/eQpLtetLEPTkjaR9InU3hSYCDwKzAGOSd2OAWan9hzgCEkbShoN\nbA/Mj4gXgJWS9kwXb30buLb0mq51TaW4+MvMzGxAyTU0PRK4Jl34PAy4IiJulrQAuErSccAy4DCA\niFgs6SpgMbAamF46lZ0OXAJsDMyNiBvT8pnAZZKeBF4Fjqhjx8zMzPpjQAxN5+KhaTOz/vPQ9LqV\n+8+XzMzMhjQXYjMzs4xciM3MzDJyITYzM8vIhdjMzCwjF2IzM7OMXIjNzMwyciE2MzPLyIXYzMws\nIxdiMzOzjFyIzczMMnIhNjMzy8iF2MzMLCMXYjMzs4xciM3MzDJyITYzM8vIhdjMzCwjF2IzM7OM\nXIjNzMwyciE2MzPLyIXYzMwsIxdiMzOzjFyIzczMMnIhNjMzy8iF2MzMLCMXYjMzs4xciM3MzDJy\nITYzM8vIhdjMzCwjF2IzM7OMXIjNzMwyciE2MzPLyIXYzMwsIxdiMzOzjFyIzczMMnIhNjMzy8iF\n2MzMLKNGF2JJkyQtkfSkpJ/ljsfMzKynxhZiSesD5wGTgJ2AaZLG5I1q4Ors7MwdwoDhXHRzLro5\nF9YujS3EwDhgaUQsi4hVwJXAwZljGrD8IdPNuejmXHRzLqxdmlyItwKeLT1enpaZmZkNGE0uxJE7\nADMzs74oopn1StKXgRkRMSk9/jnwfkScVerTzJ03M2uziFDuGJqiyYV4GPAEsD/wb2A+MC0iHs8a\nmJmZWcmw3AG0S0SslnQicBOwPjDTRdjMzAaaxp4Rm5mZDQZNvlhrjSoTe0g6Nz3/sKTd6o6xLn3l\nQtKOku6V9F9JJ+eIsS4VcnFkOh4ekXSPpF1yxFmHCrk4OOXiQUkPSNovR5ztVnUSIEl7SFot6dA6\n46tThWOiQ9KKdEw8KOn0HHE2QkQ0+h/FsPRSYBSwAfAQMKZHnwOBuam9J3Bf7rgz5mIEMBb4FXBy\n7pgz52Iv4JOpPWmIHxeblto7U/yNfvbY685Dqd/twPXAlNxxZzwmOoA5uWNtwr+hcEZcZWKPycCl\nABExD9hc0sh6w6xFn7mIiJcjYgGwKkeANaqSi3sjYkV6OA/YuuYY61IlF2+VHn4ceKXG+OpSdRKg\nHwF/BV6uM7iaVc2Fr5xeB4ZCIa4ysUerPk380PUkJ936m4vjgLltjSifSrmQdIikx4EbgJNqiq1O\nfeZB0lYUBemCtKipF9lUOSYC2Dv9ZDFX0k61Rdcwjb1quqTqG6XnN7smvsGauE8fVuVcSPoqcCww\nvn3hZFUpFxExG5gtaQJwGfD5tkZVvyp5OBs4JSJCkmjuGWGVXCwEtomItyUdAMwGdmhvWM00FM6I\nnwO2KT3ehuLb3dr6bJ2WNU2VXAwVlXKRLtD6MzA5Il6rKba69eu4iIi7gGGSPt3uwGpWJQ9fAq6U\n9BQwBThf0uSa4qtTn7mIiDci4u3UvgHYQNLw+kJsjqFQiBcA20saJWlD4HBgTo8+c4CjYc2MXK9H\nxIv1hlmLKrno0tRv+l36zIWkbYFZwFERsTRDjHWpkovt0hkgknYHiIhXa4+0vfrMQ0R8LiJGR8Ro\nit+JT4iI3t5Dg1mVY2Jk6ZgYR/HnsP+pP9TBr/FD09HLxB6SfpCe/2NEzJV0oKSlwFvAdzOG3DZV\nciFpC+B+YDPgfUk/BnaKiDezBd4GVXIB/AL4FHBB+rxZFRHjcsXcLhVzMQU4WtIq4E3giGwBt0nF\nPAwJFXMxFThB0mrgbRp4TNTFE3qYmZllNBSGps3MzAYsF2IzM7OMXIjNzMwyciE2MzPLyIXYzMws\nIxdiMzOzjFyIrdEkHbS229m1aZudXZNeSFrW12xDVfr06D9B0iJJCyVt9FHjNbO8XIit0SLiuog4\n66OuR1J/Jr+JXtpV+ldxJPDriNg9It7pWtjPGM1sgHAhtkEpTb23RNLFkp6QdIWkiZLukfQPSXuk\nft+R9PvUvkTSOanPPyVNScsl6beSHpX0iKTD0vIOSXdJuhZY1CKG8yXdL+kxSTMqxHyUpHnpJup/\nkPR/77+++kj6HvBN4ExJl0vatxTjY5LWS/syP90V5/ulfTwv5ewWSX8r7f+aM3JJYyXdkdqbSroo\nxbOwa07llNNZkm5IuT6rFN8kSQ9IeihtR6nPZ9Lz66m40XzT5qk2+9D8DdoGs+0opl5cTDEt5+ER\nMT4VjFOBb7R4zRapzxiKuXOvBg4FdgV2AUYA90u6M/XfDfhCRDzdYl2nRcRrktYHbpW0c0Q82irQ\ntL3DgL0j4j1J51Oc2V7Wnz4RcaGk8cB1ETFLUkc5xlR4X4+IcWnY+m5JNwO7U9wZZwywRcrZzK7V\ntswunAbcFhHHStocmCfp1vTcrsAXgXeBJySdm9p/AiakWDZPdym6PO3HOcDXgIcaOE+12YfmQmyD\n2VMRsQhA0iKgq0g8Boxq0T8obtVGmjd3ZFr+FeAvUcz3+pKkvwN7ACuB+b0UYYDDJR1P8T7akqLI\ntSrEAvanuHPPAhXzVm8MvNDPPj3X2aUc40RgZ0lT0+PNgO2BCaV9fF7S7b2st2wicJCkn6THGwHb\nUuTxtoh4A0DSYop8Dwfu7IolIl5Pr7sIuJaiEB8LXFxh22ZDhguxDWbvlNrvU5yRdbV7O7bfLbW7\nillAr/ejfqvVSiSNBk4GxkbECkkXAx/rI95LI+LUddCnHF+rGE+MiFt6xHsgH9zHcns13T9T9dyH\nQyPiyR7r2pMP5v49iny3PLOOiOWSXpS0H8UXnGmt+pkNVf6N2Azuoji7XU/SCGAfYD5rvxXkZhQF\ncGU6sz5gLX0DuA2YmtaPpOEqbrPYnz5lvcV2EzBd6cItSTtI2gS4s7SPWwIdpdcsA8am9pQe6zpp\nzQal3day7QDuA/aRNKor/tLzFwKXA1eF7zRj9gE+I7bBrOcHequrlaOX5WvaEXGNpL2Ah9Oyn0bE\nS+k3297O8h6W9CCwBHgWuHutgRZD4acDN6cLsFYB04Fn+tOnwv5dSDFMvFDF+PZLwCFpH/ej+G34\nGeBeugvqL4GZklYCnaX1nQmcLekRii/t/wImt9hmV/yvpN+oZ6X4XwS+np6+jmJI2sPSZj34Nohm\nQ1AaSr8+Iq6uaXtjgd9FxL51bM9sMPEZsdnQVcu3cEmnAD8EvlXH9swGG58Rm5mZZeSLtczMzDJy\nITYzM8vIhdjMzCwjF2IzM7OMXIjNzMwyciE2MzPL6H/5K2xKj8F+7AAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 } ], "metadata": {} } ] }