{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from collections import defaultdict\n", "\n", "import numpy as np\n", "import seaborn as sns\n", "import matplotlib.pyplot as plt\n", "\n", "from Bio.PopGen.GenePop import Controller as gpc\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "ctrl = gpc.GenePopController()\n", "my_pops = [l.rstrip() for l in open('hapmap10_auto_noofs_2.pops')]\n", "num_pops = len(my_pops)\n", "pop_iter, loci_iter = ctrl.calc_allele_genotype_freqs('hapmap10_auto_noofs_2.gp')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "num_pops = len(my_pops)\n", "fix_pops = [0 for i in range(num_pops)]\n", "num_reads = [defaultdict(int) for i in range(num_pops)]\n", "num_buckets = 20\n", "MAFs = []\n", "for i in range(num_pops):\n", " MAFs.append([0] * num_buckets)\n", "for locus_data in loci_iter:\n", " locus_name = locus_data[0]\n", " allele_list = locus_data[1]\n", " pop_of_loci = locus_data[2]\n", " for i in range(num_pops):\n", " locus_num_reads = pop_of_loci[i][2]\n", " num_reads[i][locus_num_reads] += 1\n", " maf = min(pop_of_loci[i][1])\n", " if maf == 0:\n", " fix_pops[i] += 1\n", " else:\n", " bucket = min([num_buckets - 1, int(maf * 2 * num_buckets)])\n", " MAFs[i][bucket] += 1" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7MAAAIhCAYAAACcx+2TAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XucXHV5+PHPACYQSVZUQpSiaat5pCAKSBBQCYoXbC0t\ntlpFVH6Vi0WUVqU2KBiLJVWRikVqiYJ3LYiiUhDQKhirEaVQxD4BMbVWWpBCuJOEzO+P79kwrHtJ\n9szuzNn5vF+v88rMmXN55snZnX3m+z3fb6vdbiNJkiRJUpNs1esAJEmSJEnaUhazkiRJkqTGsZiV\nJEmSJDWOxawkSZIkqXEsZiVJkiRJjWMxK0mSJElqHItZSVLfiIiFEbExIr49ymvnVq89thexdUv1\nHu8e47VjIuIvpzum6tx7R8T51eN9IuLs6vGzhtdLktRPtul1AJIkjfAA8NSIeFJm/hwgIh4NPAeY\n0ZOjZ+ZHe3juHwJ/XD3dDfiNav3VHeslSeobFrOSpH7zEPAF4HDgtGrdYcCXgbcObxQRRwPHV9v/\nL/CmzLwxIs4D1gJPB3YB/gP4k8y8NyKeC7wPmAOsA96ZmV+PiNcDLwe2BRYCPwfOAt4ELAI+mJkf\nrM77LuBPgA3A6uq8/xsR3wJuB54GnA38EXADsDfweOBTmfnuKvytq5bPxcBjgLdn5oUR8W7gcZl5\nfEQsAj4K7AhsBE7NzH/qTFS1/e7VNguAa4E/zcy7I2I34O+Bx1K+BDg9Mz8VEdsD5wJPqY77Q+AY\n4EDgw8AhwHuAeRHxMeBTwIcz8+kRMVTl5RnVMS8BlmbmQxHxQPX/9ULgicCHMvNDEbEA+CTwuCrs\nizPzZCRJqsluxpKkfvQp4DUdz18LnDf8JCKeD7wdWJKZzwQ+Syl2h+0FvBjYlVJY/XFEPA44H3hz\nZj4DeB3w6YhYWO3zHOD1lOJ1J+CVmfl84KXAqdV5jwReAjyrOsb1HXG1gf/LzN0y8++r5wuBA4A9\ngVdGxO9W224LXJaZe1MK9Pd1xD7c+vx54AuZuXsVw99UhehI+1FaTp8GrAdOjoitga9QCspnUArU\nv4mIZwN/CGyfmXsC+1TH+K3hg2XmL4B3AVdl5p+OONeZwG2Z+XTgWZSi9m3Va7Oq155DKeSXR8Rs\n4Cjgp9V7fS6l1X3uKO9DkqQtYjErSeo7mfkjYGNE7BURuwBzM/PH1cstSkH5+cy8vdr+E8DOVWHa\nBi7NzPWZuQH4d0rr5L7ATZn5g2qfG4CVwJJqnx9k5n9nZhv4GXBZdb6bgW2rrs6HAB/PzPur184E\nXhARj6qeXzXirfxjZm7IzLsohfSLq3Oty8wvVdtcC8zv3CkidgD2AFZUsf4iM5+SmfeMOH4bOD8z\nb63i/lh1jkXA7Mz8crX/LcAXq7xdBewWEf8CvAP4u8z86YjjthjdSyitvWTmOuAfqpwMu6j69xpg\nNqUF/BLg5RFxMaUF+B2ZOeo9w5IkbQmLWUlSvxpunX0NpZtqpxa/XnC1gOGi8oGO9e0xtgfYmodv\nuXlwxGsbRtl+5HG2qvYfXjey2HxoxLmGj7l+lPhGO/eme4Qj4qkRse0oMY08x0OM/vm+NbBNZq6h\ndDE+DZgHXBERLx9l+9FsNSLWzvwB3A9QFdYAreqe298E/pHSUr0qIvbbzPNJkjQmi1lJUr/6NPAK\n4JWUbsTD2sDXKd12Hw+buv/+CriJ0YvWNvC9smnsU+2zG6Xb67fG2Ge0Y3wdODIi5lTr3gx8u2ql\nZMRxWsDhEdGqWlr/GPjqZpyrVbVc/pDS7Zmqdfq7lOLzEdsCvx8RQxGxFaVL71eABNZFxB9W+z+R\nct/x5RFxLHBuZl6Wme+o3tNuPHJwrQ08/MVAp68Dx1XHnA0cDVw+3nuJiOXAuzLzIuAE4MfAUyfI\ngSRJE7KYlST1mzZAZv6SMoDS6sy8c8RrVwBnAN+MiOuBI4Dfq1oEh5dHqLok/zHw4Yi4DvgM8PrM\nvGmMfdqjPP4YcAWldfEG4JmUgarG2mc2sIpSSH8kM/9llO06n3fG8WrgFRHxb5QC9U8z89ZR9vsf\n4GLgJ5SBr/6m6l79B8BbIuJaSsG5LDO/TWnl3joiboiIHwBzgQ9RCuPhc38XeFpEfHFETG8G5kfE\nvwPXVed87zjvqU35f3pmtc8PKN22P4ckSTW12u0ZPcuBJEk9Ud2TevbIEYi7fI53Aztl5hun6hyS\nJPWrcafmqQa0+DjwZMq3y6cCvwC+RpmOAMo3zedHxFGU7kYbKNMHXBwR21G6ie0I3A28LjN/NSXv\nRJKkwTNqK7QkSYNg3JbZat69PTLzL6r7fa4FlgFDw/PtVdstoIz6uDewHfAdypD9b6IM//+eiHgl\nsF9mnjBVb0aSJEmSNBjGbZmlTCNwQfV4K8roi3tTBtA4FLiRMpjDYmBlZq4H1kfETZQpBQ4A/rba\n/1LKvHWSJEmSJNUy7gBQmXlvZt5TTW5+PnASZSCLt2XmgZRBHE6hDB6xtmPXu4EhyqiLd41YJ0mS\nJElSLRO1zA5PB3AhcFZmfj4ihjJzuHD9EvBh4EpKQTtsLnAnpZCdO2LdRK6nTBEgSZIkSZp5NmdK\nvAlNNADUTpR7Yf+sYzqBSyPizZn5A+Bg4GpKa+17qznntgV2pRSlK4GXUobiP4RS9E5k98m8kWk0\n2uT22nLmsT5zWJ85rM8c1mcO6zOH3WEe6zOH9ZnD+gYmhxMNAPUhypx82bH6HcDplPtnbwGOrroi\nv4EymvFWwHsz80vVaMafAJ4APAi8epQ58ppmYC6OKWYe6zOH9ZnD+sxhfeawPnPYHeaxPnNYnzms\nb2By6DyzW25gLo4pZh7rM4f1mcP6zGF95rA+c9gd5rE+c1ifOaxvYHI47gBQkiRJkiT1I4tZSZIk\nSVLjWMxKkiRJkhrHYlaSJEmS1DgTzjMrSVJTtVqtWcDCXscxnswkIhZN8WnWtNvtdVN8DkmSppWj\nGW+5gRkdbIqZx/rMYX3msL6+zmGr1Vq0+LBTcs7Q/F6H0jP3rb2VVRcui3a7vbrXsUyhvr4OG8Q8\n1mcO6zOH9Q1MDm2ZlSTNaHOG5rP9Djv3OgxJktRl3jMrSZIkSWoci1lJkiRJUuNYzEqSJEmSGsdi\nVpIkSZLUOBazkiRJkqTGsZiVJEmSJDWOxawkSZIkqXEsZiVJkiRJjWMxK0mSJElqHItZSZIkSVLj\nbNPrAKR+02q1ZgELex3HRDKTiFg0hadY0263101mR3O4yaRzKEmSpPG12u12r2NomjbQ6nUQM0Df\n5rHVai1afNgpOWdofq9D6Zn71t7KqguXRbvdXj2Z/c1h/Rw2SN/+LEO5FpcceVZuv8POvQ6lZ+65\n47/51rnHzfRrsa+vwwYxj/WZw/rMYX0Dk0NbZqVRzBmazyD/8dsN5lCSJElTyXtmJUmSJEmNM27L\nbEQ8Cvg48GRgNnAq8BPgPGAjcD1wXGa2I+Io4GhgA3BqZl4cEdsBnwZ2BO4GXpeZv5qi9yJJkiRJ\nGhATtcweDtyWmc8DXgKcBZwOLK3WtYBDI2IBcDywP/Bi4LSImAW8Ebi22vaTwDun5m1IkiRJkgbJ\nRMXs+cDJHduuB/bKzCurdZcABwP7ACszc31m3gXcBOwBHABcWm17abWtJEmSJEm1jNvNODPvBYiI\nuZTC9p3ABzo2uRsYAuYBa8dYf9eIdZIkSZIk1TLhAFARsQvwTeCTmfk5yr2yw+YBd1IK1rkd6+eO\nsn543USupwwn3a8LfRDDTFj6No+ZmYgqD+awhjo5bNBCH8TgtTiBAbgW6YMYZsJiHs1hPyzmcDBy\n2BXjFrMRsRNwGXBiZp5Xrb4mIg6sHh8CXAmsAp4bEbMjYgjYlVKUrgReOmLbiexOuRe3Xxf6IIaZ\nsPRtHiMiEFUezGENdXLYoIU+iMFrcQIDcC3SBzHMhMU8msN+WMzhYOSwKyaaZ3YppWvwyRExfO/s\nW4AzqwGebgAuqEYzPhO4ilIgL83MByPibOATEXEV8CDw6m4FLkmSJEkaXBPdM/sWSvE60pJRtl0B\nrBix7n7gFTXikyRJkiTp10x4z6wkSZIkSf3GYlaSJEmS1DgWs5IkSZKkxrGYlSRJkiQ1jsWsJEmS\nJKlxLGYlSZIkSY1jMStJkiRJahyLWUmSJElS41jMSpIkSZIax2JWkiRJktQ4FrOSJEmSpMaxmJUk\nSZIkNY7FrCRJkiSpcSxmJUmSJEmNYzErSZIkSWoci1lJkiRJUuNYzEqSJEmSGsdiVpIkSZLUOBaz\nkiRJkqTGsZiVJEmSJDWOxawkSZIkqXEsZiVJkiRJjbPN5mwUEfsCyzPzoIjYE/gqcGP18kcy8/yI\nOAo4GtgAnJqZF0fEdsCngR2Bu4HXZeavuv4uJEmSJEkDZcJiNiJOBF4D3FOt2hv4YGZ+sGObBcDx\n1WvbAd+JiMuBNwLXZuZ7IuKVwDuBE7r7FiRJkiRJg2ZzWmZvAg4DPlU93xtYFBGHUlpnTwAWAysz\ncz2wPiJuAvYADgD+ttrvUuBdXYxdkiRJkjSgJrxnNjMvpHQdHvZ94G2ZeSBwM3AKMBdY27HN3cAQ\nMA+4a8Q6SZIkSZJqmcwAUF/KzGuGHwN7UgrWuR3bzAXuHLF+eN1ErgfafbzQBzHMhKVv85iZiajy\nYA5rqJPDBi30QQxeixMYgGuRPohhJizm0Rz2w2IOByOHXTGZYvbSiNinenwwcDWwCnhuRMyOiCFg\nV0pRuhJ4abXtIcCVm3H83YFWHy/0QQwzYenbPEZEIKo8mMMa6uSwQQt9EIPX4gQG4FqkD2KYCYt5\nNIf9sJjDwchhV2zWaMaV4Qr6WOCsiFgP3AIcnZn3RMSZwFWUAnlpZj4YEWcDn4iIq4AHgVd3K3BJ\nmslardYsYGGv45hIZhIRi6bwFGva7fa6KTy+JElqqM0qZjNzDbB/9fha4DmjbLMCWDFi3f3AK2pH\nKUmDZ+Hiw07JOUPzex3HuI5ZfgVLjjxrSrry3rf2VlZduCyA1VNxfEmS1Gxb0jIrSZpGc4bms/0O\nO/c6DEmSpL5kMTvD2DVxE7smSpIkSTOYxezMY9dEuyZKkiQ9gg0em9jgMYNYzM5Adk2UJEnSCDZ4\n2OAx41jMSpIkSQPABg/NNJOZZ1aSJEmSpJ6ymJUkSZIkNY7FrCRJkiSpcSxmJUmSJEmNYzErSZIk\nSWoci1lJkiRJUuNYzEqSJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiVJEmSJDWOxawkSZIkqXEsZiVJ\nkiRJjWMxK0mSJElqHItZSZIkSVLjWMxKkiRJkhrHYlaSJEmS1DjbbM5GEbEvsDwzD4qIpwDnARuB\n64HjMrMdEUcBRwMbgFMz8+KI2A74NLAjcDfwusz81RS8D0mSJEnSAJmwZTYiTgTOAWZXqz4ILM3M\n5wEt4NCIWAAcD+wPvBg4LSJmAW8Erq22/STwzu6/BUmSJEnSoNmcbsY3AYdRCleAvTLzyurxJcDB\nwD7Aysxcn5l3VfvsARwAXFpte2m1rSRJkiRJtUxYzGbmhZSuw8NaHY/vBoaAecDaMdbfNWKdJEmS\nJEm1TGYAqI0dj+cBd1IK1rkd6+eOsn543USuB9p9vNAHMYy5ZGZuRo5nvCoP5rAGc1ifOayvTg4x\nj5vUzWMDFvoghpmwmMcZnEN/HxYD8Puwr6/DETHWNpli9pqIOLB6fAhwJbAKeG5EzI6IIWBXSlG6\nEnjpiG0nsjul9bdfF/oghjGXiIjNyPGMV+XBHNZgDuszh/XVySHmcZO6eWzAQh/EMBMW8ziDc+jv\nw2IAfh/29XU4IsbatqSYHa6g3wosi4jvUkZDviAz/xc4E7gK+AZlgKgHgbOB3SLiKuANwLJuBS5J\nkiRJGlybNTVPZq6hjFRMZt4ILBllmxXAihHr7gdeUTdISZIkSZI6TaabsSRJkiRJPbVZLbOSJEma\nnFarNQtY2Os4JpKZRMSiKTzFmna7vW4Kjy9pwFjMSpIkTa2Fiw87JecMze91HOM6ZvkVLDnyrCkZ\n8fa+tbey6sJlAayeiuNLGkwWs5IkSVNsztB8tt9h516HIUkzivfMSpIkSZIax2JWkiRJktQ4FrOS\nJEmSpMaxmJUkSZIkNY7FrCRJkiSpcSxmJUmSJEmNYzErSZIkSWoc55mVJEmSpAm0Wq1ZwMJexzGR\nzCQiFk3hKda02+11U3j8zWYxK0mSJEkTW7j4sFNyztD8XscxrmOWX8GSI8/KqTj2fWtvZdWFywJY\nPRXH31IWs5IkSZK0GeYMzWf7HXbudRiqeM+sJEmSJKlxLGYlSZIkSY1jMStJkiRJahyLWUmSJElS\n41jMSpIkSZIax2JWkiRJktQ4FrOSJEmSpMaxmJUkSZIkNc42k90xIn4ErK2e3gycBpwHbASuB47L\nzHZEHAUcDWwATs3Mi2tFLEmSJEkaeJMqZiNiW4DMPKhj3VeApZl5ZUScDRwaEd8Djgf2BrYDvhMR\nl2fmuvqhS5IkSZIG1WRbZp8BzImIr1fHOAnYKzOvrF6/BHgR8BCwMjPXA+sj4iZgD+DqemFLkqTp\n0Gq1ZgELex3HeDKTiFg0xadZ0263/TJekvrIZIvZe4H3Z+bHIuKpwKUjXr8bGALm8XBX5M71kiSp\nGRYuPuyUnDM0v9dxjOmY5Vew5MizcqqOf9/aW1l14bIAVk/VOSRJW26yxexq4CaAzLwxIm4H9ux4\nfR5wJ3AXMLdj/VzgjgmOfT2w2yTjmi7tXgcwlszkmOVX9DqMnsvMSf9RYw4Lc1ifOayvTg6r/c0j\n3bkWt99h526G1Dj+PNdX9+e5Qfry70Svw8Kf5fq69LPc6sIxJj2a8ZHA6QAR8URKkXpZRBxYvX4I\ncCWwCnhuRMyOiCFgV0qxOp7dKW+uXxf6IIYxl4iICfI7EKo8mMMazGF95rC+OjnEPG7itVifOayv\n7s9zQxb6IAavw3H4s1xfl36Wu2KyLbMfA86NiOF7ZI8EbgfOiYhZwA3ABdVoxmcCV1EK56UO/iRJ\nkiRJqmtSxWxmbgCOGOWlJaNsuwJYMZnzSJIkSZI0msl2M5YkSZIkqWcsZiVJkiRJjTPZe2YlSZKk\nadGE+Y5hWuY8dr5jqYPFrCRJkvpd3893DFM757HzHUu/zmJWkiRJfW/O0PyBn+9Y0iN5z6wkSZIk\nqXH6rmW21WpN5X0GtU3DvRDg/RCSJEmSNK6+K2b7/X6IqbwXArwfQpIkSZI2R98Vs94PIUmSJEma\niPfMSpIkSZIax2JWkiRJktQ4FrOSJEmSpMaxmJUkSZIkNY7FrCRJkiSpcSxmJUmSJEmNYzErSZIk\nSWoci1lJkiRJUuNYzEqSJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiVJEmSJDXONlN9gojYCvgIsAfw\nIPCGzPzpVJ9XkiRJkjRzTUfL7B8AszJzf+AdwOnTcE5JkiRJ0gw2HcXsAcClAJn5feBZ03BOSZIk\nSdIMNuXdjIF5wF0dzx+KiK0yc+NoG9+39tZpCKl/deP9m0NzWJc5rM8c1tet928evRbrMof1mcP6\nzGF95rC+fnv/rXa7PaUniIjTge9l5vnV8//KzF2m9KSSJEmSpBltOroZrwReChARzwaum4ZzSpIk\nSZJmsOnoZvwl4IURsbJ6fuQ0nFOSJEmSNINNeTdjSZIkSZK6bTq6GUuSJEmS1FUWs5IkSZKkxrGY\nlSRJkiQ1jsWsJEmSJKlxpmM040aJiBOBE4DfzMwHI2JH4B+A7avlBuB44O3AA5n5/mq/S4ANmfmy\n6vm7gTsy80PT/y56KyJ2A/4WmEPJ2T8DnwCuBX7YsWkbOBg4GbglMz/acYzvAa/IzJ9PV9z9ZJwc\nfjYz9+vY7lhgp8xcFhHrKFNhtap9zsjMz0x78H2iZg4BtgO+npmnTG/k/Scifgt4H7AzcB9wP3Ai\nsBiIzPyriFgDLMrMddU+TwPOzsyDehJ0n4qIbwGfA04Cfkr5PbgtcAbwE+DD1ab7Ad8HNgLvz8x/\nnvZg+0BELAGOycxXRcRsYA3wxMxsR8R+wHeAfTPz6ojYFrgJ2IXymbKgOsbTgIuAYzPzX3rxPnqp\nyuE3gVdl5hc61l9H+Uw+CPhPyrU27K3A44APAIurv4d2Bi4BXpyZt0xT+H1hM3K4hPK7cF1ELAAu\npfzO/CXV9Vtt/0fAKcAhmfmLaX0TfSYiPgDsDSygfE7fDNwGXAy8lvK3zCxgWWZeHhGvp/q86U3E\n/Scizgd+mJnLq+dzgauBHwO/DfwfJY+PA07PzPNmYh4tZn/dayh/aPwJ5Q/ftwOXDRdaEXEGcCxw\nGeWPkfdHxBxgHvDoiJidmQ9SfrEdPf3h91ZEPIaSvz/MzJ9GxFbA+cCLgB+P9odtRIw2pPbADrM9\nQQ5HavNwrm4fzm9EzANWAwNZzHYjh9Vx/iEi3pSZfz/lQfep6vfbRcAbMvP71bp9gLOA8zo2Hfkz\nO7A/wxMYzsunM3MpQETsAFyXmbtQCgsi4mfAC4e/HBhgm66jqqD6N2BP4EeUOez/qfr3asoXAFdW\nhW4bNn2pdQHw2uHrd0D9B+Xvmi8ARMTTKQUElByPeq1FxKXAGRHxFuDzwJ8PWiHbYbwcUq3bmfLF\n6bsy8ytVETz82qsoXxI8PzNvm66g+1Vmvg0gIl5HKa6WVn+7/BDYNTM3RMQTgFUR8ST8TBnNscAP\nI+KizPwJ8H5KA9wewNsz8zLY9BnzY8pn9ozLo92MO1S/dG4EPgocV63+H+CPIuIFEbEd8DbgQ5QP\nzt2qbV4AfAu4Cnh+9e3wTpm5evqi7xuHAt/IzJ8CZOZGyjds3+xpVM2yJTlsVctIQ8AdUxZh/+tG\nDgFOB145JRE2x8soudxUCGTmD6qif2TeWmM81q/rzM8OlBZv/bqR19HlwHOrx88HlgKHVM+XUFrE\nAIiIZ1Dmun/FgBeybUrPqCdVxQKUL+47v+wc6+f1JErr2UWUL/a/MWVR9rfNyeFC4ApKwf+Vjv2I\niCMovf5eYCE7quHrbx2lNfbPIuK3qy9OfjszZ1wB1g2ZeTvwJmBFVcMszMwzqpc7f6afQOlRNXL9\njGDL7CO9AfhYZq6OiAcjYjGl69cdlBbaxZQuTX+Wmb+IiH+LiN2BlwCfpRQQh1AumG/14g30gScA\nP+tckZn3RsR64HciorOL19WZ+fZxjjWov7y2JIdP5OEP08dVr20FPJ3ypcugmmwOR7oVePzUhNgY\nCyndYQGIiC9Tftc9AfhUx3Yt4LKIGO6qOAe4d5pibKJXR8SzKV077wOO6HE8TXE5cEpEfA64NzN/\nFhGtiHg88DzgI9V284BzKX8cP6Y3ofadLwKHUVpn9qHchvGk6rXOn92HMvNggKp17BxKXo+Z3nD7\n0ng5PB+4B9ipY/sW5cuXnSlfWj1qugJtosx8ICKeTyn8L4mIWcBySmujRpGZX4uIPwQ+DhxQrW4B\n74uIk4AnU26R/OPqtRn3t7XFbKVqgj8E2DEijqf8sfYmSlfjT2TmuRHxKOAvgb8D/ojyofo8Stem\nN1HuezoJuJ1yX8kg+k9gr84VEbGQ8sv+hjHun7sfmD1i3fY8/C3SoNnsHEbEMTz8wdnZzXgu8N2I\nuCIzVzJ4JpvDkZ4MDPR9TcB/Ac8afpKZfwAQEf9Kyc3Tqpce0VUxIgL/AAEgIranjLGwoVrVpty7\nPWPuWZoumfnvEfFUypfIw5+zX6d0z942M/+3WreR0kPj8cAXI2LfAW4RG26J+RxwdkTcTOlJ1vna\nWN2MF1J6pL0d+HREHFT1dBk04+Vw2JGUL0BXRsSPMjOr9bdQxgc5ipLDQ2xpHF3VrXi7zDy+ev5U\n4NKI+E5vI+t7nwTmdNwC0KbqZhwRh1C+dLm5Z9FNMbsZP+w1wIrMfHFmHgLsS7m/7q3A4QCZuZ7y\n7cYD1T6XA68CVmdmOzPvp7TiPh8Y1K44XwNeUg0YQ/UFwAd5uEv2aH4E/H5EbF3t89vA7Mz81VQH\n26e2JIdjdZG9B7iTwf0WuHYOq/ts30b542WQXQQcHBH7Dq+IiKcAv0Ep9scy47oy1XAe8JzqmpoP\nPLq34TTedZSeVMPF7CXAW4DOHhf3ZOZ/ZeY1lPu7Px0RA31NZubPKNfem3lkrwoY/XfgLMr9oSdU\ng1n+nDJ40cAaI4fDubu+GtTpL4Dzq1vOAG7KzHWZeRalp8BJ0xlzwzwB+FT1BSCUa+5XwIO9C6kR\nRvvd1gLIzEuALwP/OM62jWbL7MP+lFLQApCZ90fEFymtEodFxAmUIvZW4I3VNj+tbvY/p+M4lwO/\nn5kD2b0uM++ubuY/p/rDbS7wFcofG8tHdO8EeH1mXhERz6HcxH4X5QdtYLvcTZDDkXnpHLzosVV+\nh0dH/X5mfmt6ou4vXcjhRsoXAZdl5senKey+VHXPfhnl5/cJlM+NhyjdwLYHYpzdbX0oTgfOrB6f\nTxlhcscJ9jF3RRt4UUT8oGPdV4EDMvM/quc/oPQQWDpiPwAy8/SIeBHwTuCvpzjeftT5O+4LwGsy\n86bqi+Ph1zq7GUO5TeX5lAG1hu9D/jPK5/Q3MvPKaYq9X2xODgHIzC9GxEsoXbM/wSN/lv8fcE1E\nXJWZ356e0BuhDZCZP4qIvweujIj7ga2BczLzxojYH38vjuUR12DHumF/TbnuXjrGto3Wardn1PuR\nJEmSJA0AuxlLkiRJkhrHYlaSZpiIWBgRGyPi17qxRcS51WuPneSxl1XTTMxIEfHuiDh7Gs93UkT8\nZ0R8bAv2+VZEvLz6f75nc7evF2l9EfG1qvt/L2NYExF7TbzlqPvO6GtfkprIe2YlaWZ6AHhqRDwp\nM38OEBGPBp5DjftlMnOmDwAz3ffe/D/gVZn53S3Yp/Oep82Jt1/ukeqHOLz2JWkGsZiVpJnpIcpA\nJYcDp1XrDqOMavhW2DRi8xmU0dvnUgZfewPwr8AVlLmg/zIiDqbM2bk38D7g36tBdR6gjBL9e5R5\nPd9Omcvu6cAvgZdl5n3VwDKPz8z/q867kTJlyh5VbP9NGWn6PspoqW+mDCz1xcz8i5FvLCLWVPG8\ngDLd0hd8CmMKAAAgAElEQVSqOJcAH87Mp1fbbXoeEe8Gfhv4Lcrcwt8HLgNeB/wmcGJmfr7KQUTE\nt4DHAddQ5ha/pxrw78PVOR8FfD4zT6umL7mKMtr9QuB5HVPEEBG/AZxNGf25RZnu7QMR8QXKqNAf\nj4iTM/OfOvZ5dLXPU4HHAncDr87M1SPz0bHPSZT/462ANVXct4zYZn/KvI2Ppgx09u7MvHiUYz1I\nuVaeAbyaMlXa31U52Ro4s5qybtRrKDO/GxFPpAyA8wTKYIqP7zj+MuAPKKO73k4ZDPB/RsRwXvVw\nEWUU6MuAN1dzn+46RjxLKIMX3VO9x32qmQhGy9fRwPGUn5X/Bd5UDTSzPeX/eX9gA/DlzDypiuff\nM/P00Y4nSZp+djOWpJnrU3SM0g68ljJNzLB9gQWZ+ezM3I0yV907qjkQDwdeGxGHUiZjf1Vm3soj\nW9dmAb/MzD0oI3euoEyR8juUubp/fzNifBbw15m5K6Wg+CvgpZR5go+LiAWj7NMGHp2Zz6MUHMdH\nxHjTBA07gDI/6a7AC4FdM/NAyjzhyzq2+y3g5VVR3KKMggslnx/PzGdRcvfCiBieiH5n4D2ZGZ2F\nbOUzwDeqPB0AvCYiXpmZr6QU/a/uLGQrLwH+LzP3y8ygjNj7prHeWES8FtgdWJyZe1JG7l4xYpsd\nKF8CvCYz96bMw3p2ROwyyiEfBXwlM59GmQrnAsq18SxgCfC2arqmxYxyDVXHOAv4bmbuThkJN6o4\ndqFcJ8/KzH0oReriMd7aMyj/V79TLcdExDbjxAPli5E/ycxnjlPIPp/y5cuSzHwm8FlK8Q7wHsq1\n/TTgmcABEXEgvW9VliSNYMusJM1Q1TQHG6t7BG8D5mbmjyNi+PV/jYjbI+KNlAJuCXBX9dr/RMRR\nlCmN3pWZY01a/8Xq35sprVa3AETEzygtihP5WWZeWz3+KXBnZm4Abq+m6nos8D+j7HdRFecvI+LW\nzTzX5Zl5dxXfL4HhKUduHrH/FzPz9urxucD7I+I9wIHADhExPL3LoynF1g8oLXj/OvKEVQvr/sDB\nVbx3VS18h1BazkdVTe/xs4g4HngK5f9mvK7IvwfsA1xd/f9uDWzX8XoL2A9YAFw0fA1QWmefTmk5\nHemq6t9FlOvj4x37bQs8MzM/GhHvGu0aorSc/0X1fn4WEZdX638BXEuZKuIS4JLM/OYo529TWrHv\nBYiIT1Jac785VjxAAv+VmaO9n85cvITSsn57Fd8nIuJDVSv7C4A/r77UWV+9JyLi9eMcU5LUAxaz\nkjSzDbfO3kZpNdskIn6X0lXzA5RWqf/gkS25u1MKyX0ZW+dk9qO2glVa1TlnjbM/lKJwc9zf8bhd\nHX/432Ejz7VuxPOx4u2cb3Orar+tq+f7ZeYDABHx+CqOHYEHM3Mjv26rKqbOuLZmgs/fqjg8itLd\n9TOUrrgLx9llK2B5Zn602n8WpQvuyG1+kpnP7jjPzpQW8dEMDy61NeVLhj079lsA3DnGNXR4tVmb\nR/YA2wBQFYkHRsTelFbXMyLiXzLzhFFieKjj8dbVMcaMB3h2R9zjGfl/MrzuUYy4BqscPbAZx5Qk\nTTO7GUvSzPZp4BXAKyldKYe1KK2FX60KoB8Cf0hVtFVdNt9MuU/2MRHx5o79ttRtlFZDKPd0TpXb\ngCdFxI4R0aK04k3G70fEYyJia+BoSsvh3cD3ePh+4yFKy+W4Xak79juuY78jgMvH2w94EXBeZp4L\nrK7Os/U4238dOCoi5lbP3025X3VYu4rjqRHxvCqWPSjF5xMmiCWBByLi8Gq/XSgtq3sx+jU0XKhf\nSsnf8H3DLxg+b0RcD/xHZi6nFMN7jHLeFvCKiJgVEdtSusl/dZx49hzlGKNpU/L1yuoLCSLiSOBX\nwE2U+8VfFxGtiJhN6X3wvM08tiRpGlnMStLM1IbSDZcyMNHqzLyz47U28A+UFrJrgH+mFFgLI2Ie\npTXwTVW34dcDJ0fEM3nkfYMjH491T+GbgbMi4oeUrqC/HOMYoz3fbJl5A/BR4GpKl99f8shRf8c7\ndud2NwAXU+4V/T/KgElQBkJ6dkRcRxlA6nOZ+bnNiPtw4AUd+12QmZ8YZ3soLZ3HRMTVlO7IX6Z0\nNx4r7hXA14DvVYXiMyiDW22Smb8CXg68LyL+jfJFxxFjdMltd+y3jnJ/7Rsi4lpKIfiuagTmUa+h\natfjgN+JiBso911fWx3vOuCfKF2if0C5vv58jBjuoXxpcB3wHUqBv36MeP61Y79xZeYVlIGrvlnl\n6wjg96pW42WU1vhrgR8BF2fmlzb32JKk6dNqt/29LEmS+ktEnEvpFv2+XsciSepPk7pntup6dQ5l\nUIg2cCzlvqfzKPcaXQ8cl5ntagCRoyn3oJw62hQAkiRJkiRtiUm1zFZTNbwsM99QDVc/PA/g6Zl5\nZUScTen28z3KkPt7U0ZV/A5lKP6Rg3BIkiRJkrTZJnXPbGZeBBxTPV0I3AHsnZlXVusuoQwKsQ+w\nMjPXZ+ZdlIEVRhvkQZIkSZKkzTbpAaAy86FqrrwPUQYK6Rzh8m5gCJgHrB1lvSRJkiRJk1ZrntnM\nfH1E7ASsokxYPmweZb63u4C5HevnUlpxx3M9sFuduCRJkiRJfWsyU/39mskOAHUE8BuZeRplwviH\nKEPsH5iZ3wYOAb5BKXLfW83Tti2wK6VYHc/uk4lpGrXpUvIHnHmszxzWZw7rM4f1mcP6zGF3mMf6\nzGF95rC+gcnhZAeA2o4ycvEC4FHAaZSJ188BZlHm6DuqGs34DZTRjLcC3tsxV1tTDczFMcXMY33m\nsD5zWJ85rM8c1mcOu8M81mcO6zOH9Q1MDp1ndssNzMUxxcxjfeawPnNYnzmszxzWZw67wzzWZw7r\nM4f1DUwOJz0AlCRJkiRJvWIxK0mSJElqnFqjGUuSJElTrdVqzQIW9jqOiWQmEbFoCk+xpt1ur5vC\n40uN4j2zW25g+qBPMfNYnzmszxzWZw7rM4f1mcPu6Ns8tlqtRYsPOyXnDM3vdSg9c9/aW1l14bJo\nt9urex3LFOvb67BBBiaHtsxKkiSp780Zms/2O+zc6zAk9RHvmZUkSZIkNY7FrCRJkiSpcSxmJUmS\nJEmNYzErSZIkSWoci1lJkiRJUuNYzEqSJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiVJEmSJDWOxawk\nSZIkqXEsZiVJkiRJjWMxK0mSJElqnG0ms1NEPAr4OPBkYDZwKvAL4GvA6mqzj2Tm+RFxFHA0sAE4\nNTMvrh21JEmSJGmgTaqYBQ4HbsvMIyJiB+BaYBlwemZ+cHijiFgAHA/sDWwHfCciLs/MdTXjliRJ\nkiQNsMkWs+cDF1SPtwLWUwrWiIhDgRuBE4DFwMrMXA+sj4ibgD2Aq2tFLUmSJEkaaJMqZjPzXoCI\nmEspbE8CtgXOycxrImIpcArwb8Dajl3vBoYmOPz1wG6TiWsatXsdwAxhHuszh/WZw/rMYX3msD5z\n2B19mcfM5JjlV/Q6jJ7LzOx1DNOkL6/Dhun3HLa6cZDJtswSEbsAFwJnZebnI2IoM4cL1y8BHwau\nBOZ27DYXuGOCQ+8+2ZimSZsuJX/Amcf6zGF95rA+c1ifOazPHHZH3+YxIhYtOfKsQSnkxhQR0W63\nV0+8ZaP17XXYIAOTw0mNZhwROwGXASdm5nnV6ksjYp/q8cGUrsSrgOdGxOyIGAJ2pbS8SpIkSZI0\naZNtmV1K6S58ckScXK07ATgjItYDtwBHZ+Y9EXEmcBWlcF7q4E+SJEmSpLome8/sW4C3jPLSc0bZ\ndgWwYjLnkSRJkiRpNJPqZixJkiRJUi9ZzEqSJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiVJEmSJDWO\nxawkSZIkqXEmO8+sJEmSJA2MVqs1C1jY6zgmkplExKIpPMWadru9bgqPv9ksZmcYf8g26ZsfMkmS\nJM0ICxcfdkrOGZrf6zjGdczyK1hy5Fk5Fce+b+2trLpwWQCrp+L4W8pidubxh6zPfsgkSZI0M8wZ\nms/2O+zc6zBUsZidgfwhkyRJkjTTOQCUJEmSJKlxLGYlSZIkSY1jMStJkiRJahyLWUmSJElS41jM\nSpIkSZIax2JWkiRJktQ4k5qaJyIeBXwceDIwGzgV+AlwHrARuB44LjPbEXEUcDSwATg1My/uQtyS\nJEmSpAE22ZbZw4HbMvN5wEuAs4DTgaXVuhZwaEQsAI4H9gdeDJwWEbPqhy1JkiRJGmSTapkFzgcu\nqB5vBawH9srMK6t1lwAvAh4CVmbmemB9RNwE7AFcPfmQJUmSJEmDblLFbGbeCxARcymF7TuBD3Rs\ncjcwBMwD1o6yXpIkSZKkSZv0AFARsQvwTeCTmfk5yr2yw+YBdwJ3AXM71s8F7pjg0NcD7T5e6IMY\nxlwyMyfI70Co8tDz/48pXuiDGJq+mENz2A+LOTSH/bL0bR79+6bw7xuvw37QpeuwKyZVzEbETsBl\nwImZeV61+pqIOLB6fAhwJbAKeG5EzI6IIWBXSrE6nt0p99z260IfxDDmEhExQX4HQpWHnv9/TPFC\nH8TQ9MUcmsN+WMyhOeyXpW/z6N83hX/feB32gy5dh10x2Xtml1K6C58cESdX694CnFkN8HQDcEE1\nmvGZwFWUwnlpZq6rG7QkSZIkabBN9p7Zt1CK15GWjLLtCmDFZM4jSZIkSdJoJn3PrCRJkiRJvWIx\nK0mSJElqHItZSZIkSVLjWMxKkiRJkhrHYlaSJEmS1DiTnZpHkiQNgFarNQtY2Os4xpOZRMSiKT7N\nmna77fSCktRHLGYlSdJ4Fi4+7JScMzS/13GM6ZjlV7DkyLNyqo5/39pbWXXhsgBWT9U5JElbzmJW\nkvpQE1rDYFpaxGwN6wNzhuaz/Q479zoMSTX4ubKJnysziMWsJPWnvm8Ng6ltEbM1TJK6ys8VP1dm\nHItZSepTtobV14SWCO/3lDRd/FzRTGMxK6nrmlBAgF2ZBkTft0R4v6ckSZNjMSuNYCG2SZ1CrO8L\nCLAr06CwJUKSpJnJYlb6dRZiXSjELCAkSZI0lSxmpVFYiEmSJEn9bateByBJkiRJ0paymJUkSZIk\nNY7FrCRJkiSpcWrdMxsR+wLLM/OgiNgT+CpwY/XyRzLz/Ig4Cjga2ACcmpkX14pYkiRJkjTwJl3M\nRsSJwGuAe6pVewMfzMwPdmyzADi+em074DsRcXlmOu+iJEmSJGnS6rTM3gQcBnyqer43sCgiDqW0\nzp4ALAZWZuZ6YH1E3ATsAVxd47ySJEmSpAE36XtmM/NCStfhYd8H3paZBwI3A6cAc4G1HdvcDQxN\n9pySJEmSJEF3B4D6UmZeM/wY2BO4i1LQDpsL3DHBca4H2n280AcxjLlkZk6Q34FQ5cEc1mAO6zOH\n9dXJIeZxE6/F+upeiw1Z6IMYvA7H4c9yfeawvi79PuyKbhazl0bEPtXjgyldiVcBz42I2RExBOxK\nKVbHszvQ6uOFPohhzCUiYoL8DoQqD+awBnNYnzmsr04OMY+beC3WV/dabMhCH8TgdTgOf5brM4f1\nden3YVfUGs24MlxZHwucFRHrgVuAozPznog4E7iKUjgvdfAnSZIkSVJdtYrZzFwD7F89vhZ4zijb\nrABW1DmPJEmSJEmdutnNWJIkSZKkaWExK0mSJElqHItZSZIkSVLjWMxKkiRJkhrHYlaSJEmS1DgW\ns5IkSZKkxrGYlSRJkiQ1jsWsJEmSJKlxLGYlSZIkSY1jMStJkiRJahyLWUmSJElS41jMSpIkSZIa\nx2JWkiRJktQ4FrOSJEmSpMaxmJUkSZIkNY7FrCRJkiSpcbaps3NE7Assz8yDIuIpwHnARuB64LjM\nbEfEUcDRwAbg1My8uGbMkiRJkqQBN+mW2Yg4ETgHmF2t+iCwNDOfB7SAQyNiAXA8sD/wYuC0iJhV\nL2RJkiRJ0qCr0834JuAwSuEKsFdmXlk9vgQ4GNgHWJmZ6zPzrmqfPWqcU5IkSZKkyXczzswLI2Jh\nx6pWx+O7gSFgHrB2lPWSJEkDodVqzQIW9jqOiWQmEbFoCk+xpt1ur5vC40saMLXumR1hY8fjecCd\nwF3A3I71c4E7JjjO9cBuXYxrKrR7HcBYMpNjll/R6zB6LjOzxr7mEHPYDeawvjo5rPY3j3gtdkPd\nHB7xV59lztD8bobUdccsv4IlR55V62duLPetvZVPnfbqSe/vdVj4s1yfOayv7mdzpTXxJhPrZjF7\nTUQcmJnfBg4BvgGsAt4bEbOBbYFdKcXqmFqt1mFdjKnrMjMjIqb4NJP+5jIiFk3VB1GTRES02+3V\nk9zXHGIOu8Ec1lcnh9X+5hGvxW7oRg6332HnbofVKF6H9ZnD+sxhfXU/m7upG8XscCvlW4FzqgGe\nbgAuqEYzPhO4inJ/7tLMHLdIW3zYKdnP31xO5beWUL65XHXhsgD64gKRJEmSpH5Uq5jNzDWUkYrJ\nzBuBJaNsswJYsbnHnDM0n0H/5lKSJEmSNL46oxlLkiRJktQTFrOSJEmSpMaxmJUkSZIkNY7FrCRJ\nkiSpcSxmJUmSJEmNYzErSZIkSWoci1lJkiRJUuNYzEqSJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiV\nJEmSJDWOxawkSZIkqXEsZiVJkiRJjWMxK0mSJElqHItZSZIkSVLjWMxKkiRJkhrHYlaSJEmS1Djb\ndPuAEfEjYG319GbgNOA8YCNwPXBcZra7fV5JkiRJ0uDoajEbEdsCZOZBHeu+AizNzCsj4mzgUODL\n3TyvJEmSJGmwdLtl9hnAnIj4enXsk4C9MvPK6vVLgBdhMStJkiRJqqHb98zeC7w/M18MHAt8ZsTr\n9wBDXT6nJEmSJGnAdLuYXU1VwGbmjcDtwE4dr88F7uzyOWeczEygPZml2nfgmcP6zGF95rC+OjnE\nPG7itVifOazPHNZnDuszh/XV/Wyulq7odjF7JHA6QEQ8kVK8XhYRB1avHwJcOca+qkREAK3JLNW+\nA88c1mcO6zOH9dXJIeZxE6/F+sxhfeawPnNYnzmsr+5nc7V0Rbfvmf0YcG5EDBesR1JaZ8+JiFnA\nDcAFXT6nJEmSJGnAdLWYzcwNwBGjvLSkm+eRJEmSJA22bnczliRJkiRpylnMSpIkSZIax2JWkiRJ\nktQ4FrOSJEmSpMaxmJUkSZIkNY7FrCRJkiSpcSxmJUmSJEmNYzErSZIkSWoci1lJkiRJUuNYzEqS\nJEmSGsdiVpIkSZLUOBazkiRJkqTGsZiVJEmSJDWOxawkSZIkqXEsZiVJkiRJjWMxK0mSJElqHItZ\nSZIkSVLjbDPVJ4iIrYCPAHsADwJvyMyfTvV5JUmSJEkz13S0zP4BMCsz9wfeAZw+DeeUJEmSJM1g\n01HMHgBcCpCZ3weeNQ3nlCRJkiTNYFPezRiYB9zV8fyhiNgqMzeOtvF9a2+dhpD6Vzfevzk0h3WZ\nw/rMYX3dev/m0WuxLnNYnzmszxzWZw7r67f332q321N6gog4HfheZp5fPf+vzNxlSk8qSZIkSZrR\npqOb8UrgpQAR8Wzgumk4pyRJkiRpBpuObsZfAl4YESur50dOwzklSZIkSTPYlHczliRJkiSp26aj\nm7EkSZIkSV1lMStJkiRJahyLWUmSJElS40zHAFCNEhEnAicAv5mZD0bEjsA/ANtXyw3A8cDbgQcy\n8/3VfpcAGzLzZdXzdwN3ZOaHpv9d9FZE7Ab8LTCHkrN/Bj4BXAv8sGPTNnAwcDJwS2Z+tOMY3wNe\nkZk/n664+8k4OfxsZu7Xsd2xwE6ZuSwi1lFGD29V+5yRmZ+Z9uD7RM0cAmwHfD0zT5neyPtPRPwW\n8D5gZ+A+4H7gRGAxEJn5VxGxBliUmeuqfZ4GnJ2ZB/Uk6D4VEd8CPgecBPyU8ntwW+AM4CfAh6tN\n9wO+D2wE3p+Z/zztwfaBiFgCHJOZr4qI2cAa4ImZ2Y6I/YDvAPtm5tURsS1wE7AL5TNlQXWMpwEX\nAcdm5r/04n30UpXDbwKvyswvdKy/jvKZfBDwn5RrbdhbgccBHwAWV38P7QxcArw4M2+ZpvD7wmbk\ncAnld+G6iFgAXEr5nflLquu32v6PgFOAQzLzF9P6JvpMRHwA2BtYQPmcvhm4DbgYeC3lb5lZwP9n\n77zj3SavP/w4O0AIYe8NL3sHKGGGEQIoagnYQClQSgATCgXBjzZQCmUWKmapmWUPOwQqxN4rzBDK\n5oWw904CJGT698eRE+fmDtvXS7nn+WByrfHq69eS7aOzTrfWPmSMOYTo+6YxipsPY8xo4CVr7bnR\n837AOOANYA3ge2QelwB8a+11C+I8qjE7PwciPzT2Q374ngg8WDC0jDEXAkcCDyI/Rs43xiwELAos\nbIzpba2dhnywHV5/+Y3FGLMYMn+/sda+Z4zpBowGdgPeaO2HrTGmtSpkXbYyWQdz2JI8c+fqu8L8\nGmMWBd4BuqQxW405jMa53BhztLX2XzUX3aREn28BcJi19vlo2UDgMuC6ok1bXrNd9hrugMK83GSt\nHQVgjBkAvBr1YC9cwx8AuxZuDnRh5pxHkUH1P2BTYDzS9i8X/TsOuQHwZGTo5mHOTa3bgYMK528X\n5W3kd00WwBizIWJAgMxxq+eaMeZ+4EJjzLHAbcBxXc2QLaK9OSRatgJy4/Sv1tq7IiO4sG5/5CbB\nYGvtN/US3axYa08AMMYcjBhXo6LfLi8B61prZxpjlgNeMMasjH6ntMaRwEvGmMBa+xZwPuKA2wg4\n0Vr7IMz5jnkD+c5e4OZRw4yLiD503gWuAEZGi78E9jHG7GyM6QucAFyMfHGuH22zM/A48BQwOLo7\nvIy19p36qW8aXOARa+17ANba2cgdtkcbqipelDOHiejRkv7ADzVT2PxUYw4BfCBVE4XxwUHmco4h\nYK19MTL6W85boo2/lfkpnp8BiMdbmZ+W59FDwHbR34OBUcDQ6PmOiEcMAGPMxkh7wGQXN2TzSGTU\nypGxAHLjvvhmZ1vX68mI9yxAbuw/UjOVzU0pc7gq8DBi8N9VtB/GmN8hUX87qyHbKoXzbzrijT3K\nGLNGdONkDWvtAmeAVQNr7XfA0cDVkQ2zqrX2wmh18TW9HBJR1XL5AoF6ZuflMOAaa+07xphpxpgt\nkdCvHxAP7ZZISNNR1tpPjTH/M8ZsAOwO3IIYEEORE+bxRryAJmA54IPiBdban40xM4D1jDHFIV7j\nrLUntjNWV/3wKmcOl2ful+kS0bpuwIbITZeuSqVz2JKvgSVrIzE2rIqEwwJgjPkv8lm3HHBj0XYJ\n4EFjTCFUcSHg5zppjCMHGGO2RkI7pwC/a7CeuPAQ8DdjzK3Az9baD4wxCWPMksD2wL+j7RYFrkV+\nHC/WGKlNxxhgb8Q7MxBJw1g5Wld87c6y1u4CEHnHrkLm9Yj6ym1K2pvD0cBPwDJF2yeQmy8rIDet\netZLaByx1v5ijBmMGP73GWN6Aeci3kalFay1dxtjfgP8BxgULU4A5xljTgZWQVIk943WLXC/rdWY\njYhc8EOBpYwxf0R+rB2NhBpfb6291hjTEzgJuAjYB/lS3R4JbToayXs6GfgOySvpinwEbFa8wBiz\nKvJh/2Yb+XNTgd4tli3C3LtIXY2S59AYcwRzvziLw4z7Ac8YYx621o6l61HpHLZkFaBL5zUBnwBb\nFJ5Ya38NYIx5FpmbdaJV84QqGmMM+gMEAGPMIkiNhZnRojySu73A5CzVC2vta8aYtZCbyIXv2QeQ\n8Ow+1tqvomWzkQiNJYExxpiturBHrOCJuRXIGGPeRyLJite1FWa8KhKRdiJwkzFmpyjSpavR3hwW\n+D1yA3SsMWa8tdZGy79A6oOMQOZwqHoaWycKK+5rrf1j9Hwt4H5jzNONVdb03AAsVJQCkCcKMzbG\nDEVuurzfMHU1RsOM53IgcLW1doi1diiwFZJf5wG/BbDWzkDubvwS7fMQsD/wjrU2b62dinhxBwNd\nNRTnbmD3qGAM0Q2AC5gbkt0a44Fhxpju0T5rAL2ttd/WWmyTUs4cthUi+xMwka57F7jTcxjl2Z6A\n/HjpygTALsaYrQoLjDFrAisixn5bLHChTJ3gOmDb6JxaGli4sXJiz6tIJFXBmL0POBYojrj4yVr7\nibX2ZSS/+yZjTJc+J621HyDn3jHMG1UBrX8G9kLyQ/8UFbP8GCle1GVpYw4Lc/d6VNTpeGB0lHIG\nMMFaO91aexkSKXByPTXHjOWAG6MbgCDn3LfAtMZJigWtfbYlAKy19wH/Ba5sZ9tYo57ZufwBMWgB\nsNZONcaMQbwSextj/oQYsV8D6Wib96Jk/6uKxnkIGGat7ZLhddbaH6Nk/quiH279gLuQHxvntgjv\nBDjEWvuwMWZbJIl9MnKhddmQuw7msOW8FBcvWjya30J11OettY/XR3VzUYU5nI3cCHjQWvufOslu\nSqLwbAe5fpdDvjdmIWFgiwCmnd3V+yD4wCXR36ORCpNLdbCPzp2QB3YzxrxYtCwEBllr346ev4hE\nCIxqsR8A1lrfGLMbcApwRo31NiPFn3FZ4EBr7YToxnFhXXGYMUiaymCkoFYhD/ko5Hv6EWvtk3XS\n3iyUMocAWGvHGGN2R0Kzr2fea/lQ4GVjzFPW2ifqIz0W5AGsteONMf8CnjTGTAW6A1dZa981xmyD\nfi62xTznYNGyAmcg590ebWwbaxL5/AL1ehRFURRFURRFUZQugIYZK4qiKIqiKIqiKLFDjVlFURRF\nURRFURQldmjOrKIoilI1ouqn7wNPWWt3aLHuWuBgYElr7ffRsp5I9elXouJ7xdvPBl5HcnQLvGit\nPbzFdtcBb1lr/1HdVzNn/LuB0dba6yNNc/SXu3+L5acjxWFaFuNpb6zZwLLW2q9bLHeAXay1x5Y6\nViVErR4OBx621v6hxbpDkaJpPZB+m8cUVXDucDtjzFJIn/c1onX3ACdZa/NRl4FRSO93gMktzy9F\nUfFrY44AACAASURBVBSl66HGrKIoilJtfgHWMsasbK39GMAYszCwLfMXnvgN8AqwmTFmnaKiPgV2\nLMFwrHXxh84WzGh1f2ttpZVh56tGaa0NkcJIteZQYH9r7TPFC6Oe66cBmyIFrm4GjgPOL2O7C5GK\nsHsbY3oDDwKHIP1itwGOs9beVqPXpSiKosQQNWYVRVGUajMLqfj5W+CcaNneSHsAr8W2RwG3ABOQ\nCslHtlhfahuBX0W9bxdFjKATrLWzWnpSC8+BjYCzgPeADZBe1yOttY8bY5ZHqpAuh1S0X7K1Axpj\n/oBUt++G9Bc/2lpry9j/OuC1qNruL8hc7QosD1wctUNpjTOMMZtHxz3FWnuPMeYQYLi11jHGPA48\nAwxCeis/hXjEuwOXRsunIx7037esvm+MWRHIIK2XEkiv9X8aY7JIS6b/GGNOtdbminZzgcBa+100\nxhVIBed5jNkOtrsDeBrAWjvNGPNGpB/EmF3EGPN/SFeBE6y1r7cxP4qiKEoXQXNmFUVRlFpwI0Xt\nzoCDkH6rczDGrIf09M4hxt/vjDGLtxjnMWPMy0WP1gzDBGIADgY2ATYGRpSgcUvgn9bazYBrEI8h\nSF/SZ6y1GyDG9nzth4wxO0Svabto/4IxVtL+EcUe217AN9babYF9kFZIvdrY711r7ebI/F5fNCfF\n3t/VozDcDZF52QH4FbCDtXYja+0WiDG7YSvj3ww8Yq3dCDF8DzTGpKy1KeBz4IAWhiyIkftp0fPP\nomUtaXM7a+0dhfBpY8ymSB/3O40xCwFvAWcVvVf3Rd5+RVEUpQujxqyiKIpSday144HZxpjNjDEr\nAf2stW+02CwN3GOtnWitHQd8ABzRYpsdrbWbFj2+beVweeBGa+1Ua+0M4CbEw9kRH1lrX43+fhko\nGNI7Exne1toPkP7hxSSAPYE1gWeMMS8D/wAGGGMGlLB/WwRFWnoDbRlrl0djvwG8iRipeeZ6sfNE\nIcfW2p8Qr/cA4FVgljHmeWPM34Ex1trnigeODMRtEIMca+3k6LXMk8/cCq39npjVyrIOtzPGDAEe\nQDzdr1prp1hrdy9otdaOBn4ABnagSVEURVnAUWNWURRFqRUF7+yBwA3FKyKj6SBgkDHmA2PMB0hY\n7khjTPcKjjW76O9uSBhtgUR0zJaezqlFf7c0Bou/H9syym4sGNnAZsDWwMRW9p+vCFIbTAWw1hY8\nrG2FWBe/1gTzvtZ5xorIAwlr7STEa+0RhYIbY/7UYr9u0ZjFx+5Ox2lJHyPvX4EVmNcDW9J2xpjj\nkXNlP2vtzdGylaMCUMW09boVRVGULoQas4qiKEqtuAlIAikkL7aY3yK5j8tba1ez1q4GrA4sEm1f\nDglgP2NML2NMHyQ/9L5o3TfM9eDtXeJ49yMVewv5o4NbrM8jebn7G2OWjZaNAB6MDNGW++9c5uvp\niEOisTcD1gKeZ37Dd77nxpg9gUeAZ621pyNG40bFG1lrfwSeA0ZGx+gP/I6Ovct3AcOMMUsZYxLI\n67+znO0iQ/YoYCtr7aNF+0xB8oQHRtvtAfQFXuhAk6IoirKAo8asoiiKUm3yANbaz5Ew2HestROL\n1iWQQk8XFHkhiTyHlwDHFo9T4vHeR4oHjQeesNYWPMHHAJcZY15C8mk/b6mzlecjgfWMMW8C/0Gq\nLbd8bQ8iocUPGWNeQbzPvylh//ZeQ3vPi1ndGDMeuApIRXPbsmJya+PdB7wBvG6MeREJTz6tlfF/\nC+xsjHkVMZRvb9lWqCXW2teAvwOPIvmtM5D5wRgzzBhzT3vbRV7zvyPh1XcW5Uj/JQotTwJXGGNe\nB04GftNa2x9FURSla5HI5zv+rWCM2Qo411q7kzFmTSR/ptD/b2TUA24Ecod1JnBmVF2xL3Jnfing\nR+DgNvKdFEVRFEVRFEVRFKVkOvTMRmXwr0LulgJcAIyy1m6P3F13ozCrPyJFI4YA50R3WdPAK9G2\nNwCnVP8lKIqiKIqiKIqiKF2NUsKMJyB5RoX8m82stU9Gf98H7ILkI4211s6IKh9OQPJwBiG5Q0T/\n7lIt4YqiKIqiKIqiKErXpUNj1lp7B/NWYiwuKvEj0B9pUj+pjeWTWyxTFEVRFEVRFEVRlE7RUan9\n1ihuCbAo0oZgMtCvaHm/VpYXlnXEL8wNaVYURVEURVEURVEWLNpqP1cWlRizLxtjdrDWPoE0UX8E\nKY9/ljGmN9AHWBcpDjUW2AN4Mdr2ydaHnMt97zzWe+jaO1XlxdWJ4t6EcUE115646QXVXA/iphdU\ncz2Im16In+a46QXVXA/iphdUcz2Im16Ip+aqUE5rnkLZYw843RjzDGIM326t/Qppp/AUYtyOstZO\nAzLA+saYp4DDgNM7OsjNr95JMptepgxdiqIoiqIoiqIoShejpNY89SSZTeeBC3KpjNdoLSUSxzsh\nqrn2xE0vqOZ6EDe9oJrrQdz0Qvw0x00vqOZ6EDe9oJrrQdz0Qjw1V4VyPLN1YYmFBgAclcyml2u0\nFkVRFEVRFEVRFKU5aTpjdvh6Q0Hybk9qsBRFURRFURRFURSlSWk6Y3bHVX8F8BFwZDKbXqHBchRF\nURRFURRFUZQmpOmM2R7dewCcgbTn+XNj1SiKoiiKoiiKoijNSNMZsxE3AO8Dhyez6ZUaLUZRFEVR\nFEVRFEVpLprSmM2lMjMQ72wvYFSD5SiKoiiKoiiKoihNRtO15iEqLZ3MpnsAbwGrAGvlUpmPGiur\nTeJYCls115646QXVXA/iphdUcz2Im16In+a46QXVXA/iphdUcz2Im17opGZjzP8BfwJWs9ZOM8Ys\nBVwOLBI93gT+CJwI/GKtPT/a7z5gprXWiZ6fBvxgrb24E6+lLJrSMwuQS2VmAn8HegInN1iOoiiK\noiiKoijKgsiBwK3AftHzE4EHrbVDrLWDgJ+AI4EHge0AjDELAYsCKxljekf77QjcV0fd9KjnwSrg\nVsSQ/X0ymz4nl8p80GhBiqIoiqIoiqIo1cLxgvOBfSvdf+kBffn6h6kftlg8OvTdEzva1xizI/Au\ncAVwE3A98CWwjzFmAvAMcAIwG3GErh/tujPwOGLQDjbGPAYsY619p9LXUQlN65mFebyzPYBTGixH\nURRFURRFURRlQeIw4JrICJ1mjNkSuBC4BfHQfgbcCaxgrZ0F/M8YswGwO3Av4okdCmyNGLd1pWlz\nZgtPktl0d+A1YG3A5FKZ9xolrA26XFx9g4ib5rjpBdVcD+KmF1RzPYibXoif5rjpBdVcD+KmF1Rz\nPYibXqhQszFmADABGId4XlcA/od4Zx+z1s42xvQETgI2sdbuY4w5EnGIHgZsDvQBHkZCkF+21t5V\nhddTMk3tmQXIpTKzgNOB7sBfGyxHURRFURRFURRlQeBA4OooN3YosBWwG+ABvwWw1s5ACkD9Eu3z\nELA/8I61Nm+tnQr8AAwGHqmz/uY3ZiNGA28Av0tm02s3WoyiKIqiKIqiKErM+QNwY+FJZJiOAZ4E\n9jbGvGSMGQv8DsmbxVr7HuLBvbdonIeQqsY/10t4gYrCjI0xvYCrgTWBGcAxwM/AdYiL+nVgpLU2\nb4wZARwOzATOtNbe08HwrbrJk9n0PohRe1Mulfld2aJrR5cJRWgwcdMcN72gmutB3PSCaq4HcdML\n8dMcN72gmutB3PSCaq4HcdML8dRcFSr1zI4Aplhrt4n+vhbwgVHW2u2RyXSNMcsiPYm2AYYA50SG\ncCXcAbwKHJDMptepcAxFURRFURRFURRlAaBSY3Y94H6AqPLVCsBga+2T0fr7gF2AgcBYa+0Ma+1k\nJMF4o0oOmEtlZgN/izSfWqFuRVEURVEURVEUZQGgUmP2f8BeAMaYrYGlgIWK1v8I9Ef6Dk1qZXml\nBMDLwH7JbHr9jjZWFEVRFEVRFEVRFkx6VLjff4B1jTFPAWMBCyxZtH5RYCIwGehXtLwfUu2qI1pN\n5M2lMoz77FXOezrDr1ba/PWKlNeGputvVAKqufbETS+o5noQN72gmutB3PRC/DTHTS+o5noQN72g\nmutB3PRC/DRXJce30gJQvwKWsNbebYzZAjgfMVwvsNY+YYy5HCnN/CRS3Wog0oPoOWBja+30doZv\nN4E5mU0ngBeALYCNcqnMa2W/gOoSx4Rr1Vx7YqXX8YIhfzts6/u3WHeZ2GiOiNU8Ez+9oJrrQdz0\nQvw0x00vqOZ6EDe9oJrrQdz0Qjw1V4VKjdnFgSywMNJzaAQSsnwV0AvpRTQiqmZ8GFLNuBtwlrX2\nzg6G7/DNSGbTewD3AHfkUpnhZb+A6hLHk0c1157Y6HW8YCDwbPduie6zZufXD333zUZrKoPYzHNE\n3PSCaq4HcdML8dMcN72gmutB3PSCaq4HcdML8dRcFSoyZmtMKcZsAngWaey7aS6V+V89hLVBHE8e\n1Vx7YqHX8YK+wHigUCH8EWDX0Heb7oOhDWIxz0XETS+o5noQN70QP81x0wuquR7ETS+o5noQN71Q\ngWZjzI7Ao8D+1tps0fJXgZeAnYCPkLarBTxgCeCfwJbW2mnGmBWQ4r9DrLVfdOZFVEKlBaAaSi6V\nySOVjQFOa6AURYk7ZyKG7KWbr7M0wM7A3g1VpCiKoiiKotSDt4H9Ck+MMRsyt6hvHtjVWrtT0WO8\ntfYhpKvNhcaYnsBtwHGNMGQhpp5ZmOOdfRrpYbtFLpV5qdbC2qBL3L1pAuKmuen1Ol6wA/AY0jJr\nkyv+vPPPR5z7yHTgK2Cd0HenNFRgaTT9PLcgbnpBNdeDuOmF+GmOm15QzfUgbnpBNdeDuutNZtPn\nA/tWuv9SCy2+yjdTvv+oxeLRuVTmxLb2McbsABwJrA3sZK2dbIz5B5JCujKwI7COtXZaK/v2QIoA\nfwc8a609o1LtnSWWnllQ76yidAbHC/oB1yEf2AeFvjtl+aUWAfCBlYC/NE6doiiKoiiKUifGMDcq\nbyDwTNG6B40xj0WPhwsLrbUzkVpJuyC/JxtGbD2zMMc7+wSwHbBlLpV5sZbC2iBud5tANdeDptbr\neMEVSGG2c0LfHRUtzkdG7ttIq631Q999r1EaS6Sp57kV4qYXVHM9iJteiJ/muOkF1VwP4qYXVHM9\niJteqDxn9ghgFJABzkbSzZ4B9gd2AExrXWiMMasiocYZxBDeyVo7u+V29SC2nlmYzzt7eiO1KEpc\ncLxgKGLIvkqL6yb03Z+Q5P7ewIX1V6coiqIoiqLUC2vtB0iHmmOAG1usns9ANsb0Qrra/MlaezHw\nMXPtsboTa2MWIJfKPAY8DgxNZtNbN1iOojQ1jhcsDlwDzEDCi+fLgwByyDXlOF6wRx3lKYqiKIqi\nKPUhHz1AjNMVrbUTWqwrDjN+zBjza6SS8ZPW2vujbY8C9jfGbF9P8QViHWZcIJlNb4+EGz+YS2WG\n1ERV23SJUIQmIG6am1Kv4wW3IKEjJ4e+e3aL1XM0O16wAfA/4ANggzaM3magKee5HeKmF1RzPYib\nXoif5rjpBdVcD+KmF1RzPYibXoin5qoQe88sQC6VeRLpj7lbMpse1Gg9zYzjBck7HnsXxwsWiPde\nKR3HC/ZFDNnngPPa2zb03deBfwFrAsfXXp2iKIqiKIqilMeCZNBo7mwHOF6wJ3DbtXe/CXCj4wU9\nGyxJqROOFyyLJOlPBQ4OfXdmCbudBnwNnOJ4wUo1lKcoiqIoiqIoZbPAGLO5VGYs8CCwcxR2rBTh\neMF6wK3AtDVW7A9wAHC74wV9GipMqTmOFySAK4ElgD+HvvtOKfuFvjsROAlpnv3P2ilUFEVRFEVR\nlPJZYIzZCPXOtoLjBUsAIdAP+P25R20L8BAwDLg3aseiLLgcDDjAY0jocDncgIQlJx0v2KnawhRF\nURRFURSlUhYoYzaXyjwH3AfsmMym9Yc3EIUSjwZWB84Kffe2Pr17gBg3dwI7AQ9FVW6VBQzHC1YG\nLgZ+BH4f+m5ZPcCi7f+IFBa4VEPTFUVRFEVRlGZhgTJmI+Z4Z5PZdJes6tWCixGD9b/AqYWFUXXa\nJNJPaivg8SivUllAiIp8XQssCvwp9N2PKhkn9N1xwFXA+sDI6ilUFEVRFEVRlMqpqDWPMaYbcDWw\nNjAbGAHMAq6Lnr8OjLTW5o0xI4DDgZnAmdbaezoYvtOlpZPZdAjsBeyaS2Ue7sxYJdC0pbAdLzgK\nuAx4FRgU+u5P0ariFizdgEsQI2UCsEulRk+Nadp5boOG63W84GjgUuBuYFjoux1d7G1qdrxgSeAd\noDuwdui7X1VTaydo+DyXSdz0gmquB3HTC/HTHDe9oJrrQdz0gmquB3HTCxVoNsbsCDwK7G+tzRYt\nfxV4CdgRMNba6caYZYH7kW4YnwNHWGv3j7bfB3EmDrXWftr5l1IelXpmdwMWttZuC/wdOBvwgVHW\n2u2RyXSjF/5HYBtgCHCOMaZX52V3SME7+/eu6p11vGAwYqR+gxgyP7W2XVEY6dlIG5anHS8wdROq\n1ATHC9ZGPnC+A0aUYMi2S+i73wKnIF7eczuvUFEURVEURWkwbwP7FZ4YYzZECn9StGwF4AHgVGvt\nLS3W7Q/8GRjcCEMWoEeF+00F+htjEkB/YDqwlbX2yWj9fYjBOwsYa62dAcwwxkwANgLGdU52++RS\nmfHJbDoA3EjHA7U8XrPheMGawO2Il3zvjjytkaFzsuMFk4B/AE85XjAk9N2Xa69WqTaOF/QArgf6\nIm14vqzS0FcgURiHOF5wZei7z1ZpXEVRFEVRlC7LWHf4+cC+le7fe+mlmPb1Nx+2WDx6UDDmxHZ2\nywOvAGsbYxa11k4GDgRuBlaOtlkVCJCI20eL9sMY8zvgaGBna+2kSrV3lko9s2OBPog1fwXiASz2\ngP6IGLmLApNaWV4PTov+7VLeWccL+gN3AQOAI0PffbrUfUPfPQ84ElgSeMzxgkG1UanUmBOBrYHb\nQt8dXa1BQ9+dhXxogRSD6l6tsRVFURRFUZSGMAbYO/p7IPBM0brRwE/AMkXLEsB2iINjANDQ4qCV\nemb/D/G4nmyMWRFp+VH8QhYFJgKTkXYwBfoBP5QwfqdCIgFyqQz+2Ct5/tOXt/zzdiPLquBaAZ3W\nWw1mzc6zxbrLMO6tr3C3X4PD3A2uAa5pY/NWNYe+yxPjP+WCW8f379mj29Pj7ddsZpaunejyaIp5\nLoO66/3g80n06J5g0YV78a8TB+9HUehIibSrOfRd/Fte4vGXPt185D4bz6xcaVXR86L2qObaEze9\nED/NcdMLqrkexE0vqOZ6UFe9g4Ix1RhmlRbPT4gerXLDDTeQzWY57rjj9jvttNN48cUXrx05ciSb\nbrrpTvfeey/PP/88l156KUsssQT777//Le+///4tq6++OjfccAPnnXcet9566wqjR4/m0Ucf/Saf\nz5NIlO07rIqzsdICUGcBk621/zDGLIwUfHoHONta+4Qx5nLgEeBJpJ/pQMST+xywsbV2ejvDVy3p\nOplNb4gUPxoHbJlLZWpxYjZNkrjjBecjJ+0DwF6h77ZlbHSo2fECB7kbkwD2D333jmpqrYCmmecS\nqbtexwt6AS8AGwN7hr57b5lDlKTZ8YLlkOt9OlIM6rtytVYRPS9qj2quPXHTC/HTHDe9oJrrQdz0\ngmquB3HTC5UXgDrCWru/MeYp4Csk/3UNYH9gB+YWgBqO1CTaEon+Ky4AdRfwgrX2zGq9mHKoNMz4\nfGDr6IU/AvwFCT883RjzDOLxvd1a+xUSglzYblQHhmxVyaUyrwE5YAukuvECi+MFByOGrAX2a8eQ\nLYnQd0NgKGKwjI7GV5qbvyGG7NUVGLIlE/ruF0gY/+LAGbU6jqIoiqIoilIz8sz1QGeBFa21E1pZ\nh7V2DPA88O+W64BDgSOMMTvUXHErVOSZrTFVvRuSzKbXQzzHrwCb1cA72/C7N44XbIOEek8Btgx9\n990OdilZs+MFWyEFvQYAx4S+e2lntHaChs9zmdRVb/Q+PQN8DGwU+u6PFQxTznnRi6hoALBFA4uF\n6XlRe1Rz7YmbXoif5rjpBdVcD+KmF1RzPYibXoin5qpQqWc2NuRSmTeB24BNgF83WE7VcbxgZeBO\npP9nsgRDtixC330eCTP4CrjE8YJTHC/okhdLs+J4wULADcj1fEiFhmxZhL47HTgmOua/9JxQFEVR\nFEVR6s0Cb8xG/B1pU3NaMpteYF6z4wULI+WylwaOC333oVocJ/Td14BtgY+QsNLz1HhpKs5BPKQX\nhb77RL0OGp1vY5A+0gfW67iKoiiKoiiKAl3EmM2lMm8DtyA9bvfuYPNY4HhBN6SX6CbAlcC/anm8\n0HcnIAbt20hu7uXamqXxOF6wE+IhfRsY1QAJHtJ3+jzHCxZtwPEVRVEURVGULkqXMGYjzmDB8s6e\nCgwHngD+GPpuh7nAyWy6WzKb/ssN/xtDJb13Q9/9FNgeeBk4HLjJ8YKG9pbqykTG47XALOCg0Hen\n1ltD6LsfIZ7hZZFzUmkDxwtOO/LcR3C8wDRai6IoiqIoyoLAgmDUlUQulXkHuBFYH9i3wXI6heMF\n+yKVaz8A9onyF9slmU33Ql7/2Xfbh0GM0bIJffcbYDAwFulheqfjBX0rGUvpNBciPcXODn33xQbq\nOB94HzjW8YL1GqijaXG84ETgb5998xPAf9WLrSiKoiiK0nm6jDEbcQbixfpbMpuOZYis4wWbIeHF\nPwHDQt/9tqN9ktl0P+Ae4ADg+YV79gX4ZzKbbtlcuSRC350IDAEeBPYE7nW8oF8lYymVEfUBPhTx\nkjekr1eB0Hd/Af6EtOS6RPOp58Xxgt8D5wGf7jJwZYB1gOujVAFFURRFURSlQrrUj6lcKvMeYgiu\nC6QaLKdsHC9YFin41Ac4IPTd1zvaJ5lNLwM8DuwC3AUMPmTTJMAiwNWVhBsDhL77MzAMuAPYEXjE\n8YIlKhlLKY9onq9CegAfVIpnvg7cjbRw2pkFJC+9GjheMAx5r74Hdjt6341B2mj9GunPrSiKoiiK\nolRIlzJmI84EZiLe2R6NFlMqjhf0QVrwrAj8JfTdsKN9ktn0Wkjv0c2QH9TDc6nMlO1X3QrEU7sL\nMKJSTaHvTkNuClwHDASecLxguUrHU0rm38AywF9LuaFRD6Kc7WMRA/vCqF1Ql8bxgu2RJuTTgD1C\n332re/duINfMJ8AZjhcMbaBERVEURVGUWNPljNlcKvMBUjRnbWD/BsspiShs80pga+AmJGSxXZLZ\n9EDEkF0dOB04IpfKzARIJBIARwCTAL/ScGOA0HdnAn8ALkHykZ92vGC1SsdT2sfxgv2AJPLe+g2W\nMw9Rj2MfWIku7nV0vGBjIERCr/eO+jUDc/LOf4MY/rc4XrBmY1QqiqIoiqLEmy5nzEacBcwATo2J\nd/YE4HfAC8CIjioXJ7Pp3ZHQ4sWBI3OpzGm5VGaefXKpzGeIJ20R4KpKw40BQt+djeRMnoEYz085\nXrBupeMprRN5vf8NTAEODn13VoMltcZZwKfAiY4XrNFoMY3A8YLVgQeAfkgY+AMttwl99yXgSGAx\npIjaIvVVqSiKoiiKEn+6pDGbS2U+Aq4B1gQObLCcdnG8YC/gH8BnwK+jYjttksymD0I8Qt2QsOIr\n2tn8BuBeYFfgsM7oDH03H/ruqYjhvQLwZFSsSqkCkXf+amAAcGLU97fpiHKpTwB6I9WWuxRRXvtD\nSBj4saHv3trWtqHvXof0h94AuEYLZymKoiiKopRHlzRmI85GwvxOTWbTTdkr1fGC9YFbgF8AN/Td\nL9raNplNJ5LZ9EnMrXS8ay6V+W9740fe2sOZG268cmc1h77rR2MuATzmeMF2nR0zbiSz6d4Tvvuw\nol6+7fAHYA/gYeDyKo5bC3JIZIDjeMEeDdZSNxwv6A/cj0QnnBH67qUl7HY88DQSOn5CDeUpiqIo\niqIscHRZYzaXynyC5KGuBhzUYDnzEVWsvQsJVfx9FJbYKslsuhtwEXAuEuK5bS6VebqU40ThxsdF\nx6m4unExoe9ehbQBWgh4wPGC3Ts7ZlyI3ovbRz38D5Bc5U7jeMGqiJdzEnBoFNbdtERh8H9E2mBd\n7HhB7wZLqjlRgba7gI2Rmw1/K2W/0HdnIH2vPwfOdbxgl5qJVBRFURRFWcDossZsxDlIpdG/JrPp\nXo0WU8Dxgp7A7cz18GTb2jaZTfcGbgWOAd4AfpVLZd4o85DXIW1VdkU8gJ0m9N3bkPYjCeAuxwv2\nqca4MeBEYK/o778ms+mjOjNY1Iv0OiS3+ZjQdz/pnLz6EFVZ/hcSyn98g+XUFMcLegC3Adsj1+3R\nHeW1FxP67pfAcMT4z2oBNUVRFEVRlNKoyJg1xhxsjHksejxnjJlqjNncGPO0MeZJY8y/jTGJaNsR\nxpgXjTHPGmP2rK78zpFLZT5HvCirAIc0Vo0Q5c1divRuvQM4ra1tk9l0IawxCTwFbJdLZT4t95hF\n4caTgQuqEW4MEPruPcBQ5IZB1vGC31dj3GYlmU1vjxRA+vzMnU8E+Br4VzKb7owhfwywA9Jf+MZO\ni6wvpyFzcIrjBSs1WEtNiK7XKwAXeAQ4sJLCXKHvPgeMRIq23aGtjRRFURRFUTqmImPWWnu9tXYn\na+1OwDgkpPBUYJS1dnvEG+caY5aN1m0DDAHOMcY0jQc04h9ITuopkZez0RyFtM15BamE2mpIaTKb\nXh54krlG7265VOaHSg8aGcGFcOMrq5XvGfru48BgYCLwH8cLjq3GuM1GMpteBvHOAaTWXnJ1EEP+\nZ+DmZDa9U7ljOl6wDhI98C1weDnevmYg9N2JwElIuPk/GyynVpwDHIp8Dv4m6r1cEVF4/pXAJsCV\nWhBKURRFURSlfToVZmyM2QJYz1p7NbC5tfbJaNV9wC7AQGCstXaGtXYyMAHYqDPHrDa5VOYLpN3J\nSlQpxLZSHC/YGbgY8Wa5UWXY+Uhm0wbpM7oRoj2ZS2XarXJcItcint4hyA/0qhD67ouId/EL+rXN\n7AAAIABJREFU4CLHC05dkH6oJ7Pp7kihruWAvxTylXOpzHikn2gC+G8ym96k1DGj0NUbgD7AEaHv\nfl114fXhBuA5IOl4QdkGfTPjeIGHGOvvAHuEvvtjFYY9Bpmv30Z/K4qiKIqiKG3Q2R6ro5hb5KbY\nOPkR6A8sihStabm8I+rqgbrS/QdH330Ki/Ra+LLps2Zc1qt768WN8/k8s6ZOZcYPE5k+8Qemfz+R\nmZMnA/k8iW7RDCRIJBLyd6IbiQRE/5P/CttFyxLd5N8ffpzGJtM+Yfqs2ew/ZJ2lV1pm0Q+/e/b5\nOdvKmAk+m/wV636dYMrMX9hhtV8xaOUtjkp063bUxP+9ImMWbTv3eNFbE63r3rcvfZZdJt+t17xO\n8lwqw7dTvse7/wyAq7+d8v3VSy60eFXmOPRdvvj2Z0654hm+/n7K6e72a5yez+dFa+k0pWdy+Hp7\nMObNe9l8+Q35v23T5wHnRavyuVSGsR+/yMXP/qfnYn0Wffnrn75l6UWW7HDMA3dfh5vuf5sdN18R\n74DNx9RSfz6fZ9rXX5Po1p3eSy1Z1TkOfZcJn0zk+IufYKVl+j06c9ZsenSveqp+3c+LR8d9DMAS\n/ftw3tHbrb304guVc7OhTb2h7/LdpKn86cInmPzz9Itee+/bizZco+PzpQ405bXXAXHTHDe9ED/N\ncdMLqrkexE0vqOZ6EDe9ED/NVXFsJfL5yl63MWYx4Glr7QbR80+stStFf7uIZ/ZBYHdr7cho+R3A\nmdba8e0MnadKL65UxrrDe908dMBF3WaTXu2zaVdu/fqUccCyRY/liv7uW09tNSQPfIJ4ywuPd4EJ\nV/1mie2m9O3+b8RLu0eUU1sVHC9YAenDuS7S6/eIEnMM635elEIymx6CRCJ8BGxWFOo9j95kNn0M\n4nV/FxiUS2W+aWtMxws2BV4AvgI2DH234vDxlox1hyeANYDNgM2jx2ZI/1qi1/FE0eP9QcGYTr//\njhdcgeRlHxf67kWdHa+Iup8XUe/n/yI55tuFvltOwbWS9DpesC3wGPADsHmDC3815bXXAXHTHDe9\nED/NcdMLqrkexE0vqOZ6EDe9EE/NVaEzxuwwYGdr7bHR87sA31r7hDHmcqQYypOI4TIQCZd8DtjY\nWju9naGr8mZEP9oHML8x2pqBukQHw81EDIsvWzy+WPuE4//1zj8vSEaau0X/Fj/aXTaLRLdnFt/o\n6Cnde6+/+PTJjw6c9Padre03YcXeW365ZI99E3lmrfLF9JtX/HqGbWP8DjUsvcvO//f1w488jlSa\nXbG1Fzyld+KX7/r36NNrRv7JZX6YeS+RoQu8NygY02r4c6k4XrAUYihvhvQk/V3ou+2dE9CEF2ky\nm14JeBnJMx6US2XGFa2eT28ymz4H+DOSX7lTLpX5qeWYURubccAGwO6h7z5Qqb6x7vBuyHvc0nBt\nGR3xLjB+8a22TH3//AvfI0WICnyGXMcF49ZWYtw6XrAkEo7bHVg79N2vyh2jDep6XkRG5kPRcXcO\nfffZMocoWa/jBSORitDjEKO5GqkEldB0114JxE1z3PRCjDQ7XrDXSsv0Cz/56sfdQt99qNF6yiA2\nc1xE3DTHTS+o5noQN70QT81VoTPG7AnAdGvtJdHztYCrgF7Am8AIa23eGHMY4pHpBpxlrb2zg6Hb\nfTPGusP7AsvQsYG6LNB6vPBcfqDIOH1/+V5rfrZ0z4GLTJ1986Z26vVF674bFIxpq7dnp04exwv+\nCXiIcbdXSy9lVIjpZOAM4Dtgr1wq81ylx4uYo3msO3whpAXQmsBa0b9rzk5gEnlWaOOFfU4Lby5z\nDd2S8gYdL+gPhMB2iGdzn9B3p5SiuRlIZtM9EePuV8DIXCrz7xabtGbMJhBv9O+RqAUnl8rMY8Q7\nXnAukod5eei76VL1RIbrWsw1WjcHNkVC/Yt5B3ip6PHyoGBMIRUgP9Yd3h1YH8lxLjyWKtr/K+Y1\nbt9s59qYB8cLjgIuA64Lfbdala3rdl44XrAR8toXBoaFvntfBcOUY8wmgP8gldavBf7QoCJgTXXt\nlUjcNMdNL8REcxTp8jRSiG4iMDD03QmNVVUysZjjFsRNc9z0gmquB3HTC/HUXBUqNmZrxQ/jX86/\nefqZh9C2kdpRzu10pNDQfF7UFs+/GhSMmcfTkcymlwA+RCrQrp5LZdozrgpUfPI4XnAI8iP1bWDr\n0HeL84sLhYUuBdJI+OeQXCpjKzlWC0rSfPD1R6T7Tpv972W/nfHSLs//eGsiMnSjxyptjPElrRu6\nEwYFYyYXbxi1HxkD7I4YCU7ou5Npnaa6SJPZtI/0T70NOKCVUOxW9Saz6R5IiOqewM3AQblUZjaA\n4wXbIC2WPgQ2Dn13Ps8tQGRwrs38husiLY5vmddw/V/L96AjzVGEwzrMa9wuV7TJd8xr3L42KBjT\nati44wXdES/jJsA2FXg1S9JcC6Ler2OR135g6Ls3VzhUWXodL+iD/BDfHDgq9N1MhcftDE117ZVI\n3DTHTS/EQLPjBcshKRsr7LrlyomHXvgYon7sVSrYVmuafo5bIW6a46YXVHM9iJteiKfmqtB0xuxY\nd3hbgr6hYwP1S2BiZ3L8ktn0mYgn9IRcKuOXsEtFJ4/jBYOARxHDeavQd99toaMPYuzsjbTp2SPq\ni1sNStIceRLvB3YDDs2lMtcW1o11h/cGVqOFRzd6rErrlbK/poWh+3P33h9es9Kwk6b06PtrxODa\nPfTdbyvVXA+S2fRvkHZIFhiYS2Va+1HUpt5kNr0QEoa/NXBBLpXxHC9YGPgfks+6Q+i7T8Ecw3Ud\n5jVcN0G8gwVmIzdEWhqurRrD7dDhHEfG7ZrA9sw1bov7Ek9EjK+CcfvyoGDMzMLK6Lx/OtK4VSU9\nWcvV3FkcL1gG0bwm8KfQdy/uxHBl63W8YGVkvvoDO4W+O7YTx6+Eprn2yiBumuOmF5pcc3Sz9HEk\nzenPoe+e63jBJUiV8DuRaKCSokoaSFPPcRvETXPc9IJqrgdx0wvx1FwVms6Y/SR3e/7jm289jHkN\n1K8HBWNm1OP4yWx6ceADYBqwWi6V6Sg/tNIfpy8iubpDQt99pIWGAUCAhOA+Bvwml8pMmm+gyilZ\nczKbXhl4HTGYNoj60bbLWHd4L8Sgbc3QXQ3Jm5yHaYkev3zba7E+P/VYaNIKv3xz1SKzpr5MZPQO\nCsb8UI7mWpLMptcAxiMh7FvmUpnX29i0Xb1RFMDTiKF64tQXdl81kZ89ctUpX1yX+uKRJ5hruG6M\nhMcVmI2E8Rcbrq90No+5FM1tMdYdvirzem5XL1r9I+LRLBi3485d86D/AAcixb+ubITmUnG8YFHk\nB/GmwFmh757SySErvfm1E5Kr+w1SEKpaN7ZKoSmuvTKJm+a46YUm1ux4QTcgC+wDXAccGvrubMcL\neiEpHjsCfw1998yGiSyNpp3jdoib5rjpBdVcD+KmF+KpuSo0nTFLE7wZyWz6dOBU4KRcKnNeB5uX\nGza4CGLEbAwcHfruZS2OvSLiDV0fKY50UC6VmVaG/FIoS3Mymx4BXAnci+TsVnzSjHWH90RClOcx\ndPOwVp7EGt3It+bR/aHvCssPmPrZ5y8gFWQnlfnvT9WoyBt5y59BDJtDcqnM9e1s3lHud0+7cu/B\nny3T89YlJs4asNTnPVl6yk/5HvnZxfvMQkLiCkbreMRwLSX8vRKqVXxtReY1btcuWj1lRqL7uOcG\nbLD1J32W+aUbs9e5NJP+ohOHq9nnRRTiex/yw/dK4Mgq5Kx2Ji3hOOAC4FlgxxKKplWLhn8mV0Dc\nNMdNLzSxZscLzgBOQVIgdo2ulTyQiAoQjkN6y7uh74aNU9ohTTvH7RA3zXHTC6q5HsRNL8RTc1VQ\nY7YVktn0Ykje4kzEO9tebk05BV26AaOR0OHLkRy4OW9AMpteD3gAqTB8CXBcIZ+yypRrzCYiXbsC\nv8+lMtfVQBNjkiN6PLzkQL9HftYxS02fOGmLiW8F/Wf+vDiwZs/+/deZMWnSNKB3BUPPRjyEkyjf\nEJ7z78UHLH0xcARwTS6VOayDYxYX2eqF3JwoDhXeqPi1zErALz17fLrw9JkPMNdwfXVQMGZqBa+3\nUmpy7Y11hy+HRBkUjNv1C+tmkZjVnXxxWPKzZb7mmmiO8ntHA4WQ8mQVQqKhc8ZsArgJOIAyC4R1\nkoZ/JldA3DTHTS80qWbHCw4EbgTeQ2pRFNJW5uh1vGAzJGJkOpLu8HYjtJZAU85xB8RNc9z0gmqu\nB3HTC/HUXBXUmG2DZDZ9KnA6MCqXypzTzqblGLMFj+/jwG6h784JnU5m04OQ6r4DkGq251ezv2sL\nyp7jFuHG6+dSmc9qIQzm8UB9h+TQjiPSHOXqLho9+rfzb3vrFqX1nN52mdkNZvZIzOw9Pf9xQnJD\n2zJ+J61+5IjL37/8qiuYa7j2KhpqBvAa8NIzAzbYcMJSA7aetImdObsH04HBuVTm+XK1VYm6XHtj\n3eFLzUh03/HNfqtdvuwv3y2+9PQf8om5x52BFGwpGLfPdJD7W3XNkdF4JXAYEua/RxXb4nRKb5QH\n+AwS2TEi9N2rq6SrPZriM7lM4qY5bnqhCTUX1aKYihiyxUbqPHodL/gtcnPIIgZtNVN5qkXTzXEJ\nxE1z3PSCaq4HcdML8dRcFdSYbYNkNt0fyZ0FWDWXynSqyq7jBUkkh+d9YMvQd78rOtavgVuBHsAf\ncqnMDZ3RXgIVzXEymz4cuAK4B2krU7OTx/GCPyAGxc9Iy6InqNJ5ERUxWpiODeJFgf5TeyVW/GZA\nj517z8iz5MSZH3afzULRuoVbPcC8TAdeRTythXDh1wcFY6Y5XuAilY1f6rP5Q2cnus8ajbSLGlSl\nqtXlUtdrz/GCXYCHFpr5y/NHf5g7u9tcz+2mzL3ZMBOZs4Jx+3SLisy1MGbPAkYh79lO7VTYroRO\n640qK49DqldvH/purW9+NMVncpnETXPc9EKTaY6ui+eRPtlDW+knO5/eotZ4dyMhx81WEKqp5rhE\n4qY5bnpBNdeDuOmFeGquCmrMtkMymx4FnAX8NZfKtFUookO9jhdsjrRcmYm0BHij6BhHAP9G7iTv\nk0tl7q+G9g6o1JhNIMUzdqHjnNFO43jBvkhF51kH7r5On5vuf3sD4M169tpMZtOLIJ7CdYF9c6nM\n7YV1Y93hPYB+tGIAr3XsH29+9+JLNwPeGBSMmS+3Mcrbej3aZ7PQd99MZtOHIb2aPwK2qWL16lKp\n+7XneMHtwHDgoNB3bwQY6w7vDwxibsXkLZAbPSCRAS8juXBPrX/GaXe88dfTdovWd2/xb1t/t7ns\ng77LbfVN7wGDe82eMXHdnz68s/fsGdPLHKfd9QutsvLqUz76+BIkL/7xSsPIHS/YDcnn/QIpCPVV\nJeOUSNN8JpdB3DTHTS80keaob/kzwHpAOvTdy1vZrDVjtgdyHe0CnBH67qm11lomTTPHZRA3zXHT\nC6q5HsRNL8RTc1VQY7Ydktl0P8Q72x3JnZ3Yymbt6o363L0ILA8MC3337mjsBPC36PENsGculXmx\nuq+gTSqe42Q2vQpihM2ixuHGAI4XDAVuZ25F3w+QcOy7gCeLQ7WrTfQe3Qj8Frgkl8ocW8bubc5x\nFMZ6O5I7fULou3NaQCWz6VOAM5Aw5O3bOOdqRSOM2VWAt5DwbNOaF3SsO3wR4FfM9dxuybxh283E\nbOSm1UzkGpnn7269e684e9qcem6/ICkH9yM/qN8tp1CZ4wUnAeciN8p2ruG10DSfyWUQC82OF/QF\nTtl6g2VHPff6l8uEvvt1ozWVQVPMcWSQhki/8otD3/1TG5u2qtfxgiWQ7+jVgOGh795RK60V0BRz\nXCZx0xw3vaCa60Hc9EI8NVcFNWY7IJlNF34wnpZLZU5vZZP2jJY+SGjklsBJoe+eF43ZA/HGjkCM\nsyG5VObd1saoEZ2a48ibfDl1CDcGcLxg8eMP2Oy7C24ZPxr5wdIvWjUZMQJC4L7Qd7+v5nGLwqqf\nRwzLcqrHtndeFHK1nkLCWOcUFooM6EuBkYj3cUgulalWvmZHNOTac7zgr8DfAT/03RM62n6sO7wv\n0qN365X2S579yW25v9KG8djGsvn+fnbABltMWGjFs4Ep6//4/ojNJr9jS9mv5b8dGaOzZ8zIP7vP\nfoOBoci5vGHR6veZa9g+1lG7peimSA5pP3JJ6Lvl3Gwph6b6TC6RptfseMGOSCTGmtGit4Bd6tx2\nqTM0xRxHvWP/iFw3TjuF2tr7TN4IqRI+G8m1faO17RpAU8xxmcRNc9z0gmquB3HTC/HUXBXUmO2A\nKMz0A8QTtGoulfmhxSZt3e1NADcg/TRvBA4OfTefzKYXQvJjhyHhknvkUpkva/gSWqOzxmwC6Xm5\nM9I66MZqCWuHPNJWoRcSfjoMcIBVo/WzkJZHdwFh6LudujmQzKY3Q8LWfgY2zaUyH1eit+VCxwtW\nQDzbPYGNQt99v5VjdwduQ4yUO4BkLpWpRiXdjmiUMdsHaUG0MjInb5WxezVyULcBHo6e7hr67tjO\njNcB8+iN2hjtHj12RcLUQXKtn0R+oN8PvNWaoRy1+noOqRB9cOi7tci3b6rP5BJpWs2OFywGnI8U\nGJsNXDRsu9WPv+up90Eq8O4c+u5HDZRYKg2fY8cLjgIuQz4/tukgv72jKKp9kZtDE5C6Fi2/6xtB\nw+e4AuKmOW56QTXXg7jphXhqrgpqzJZAMps+AfnxcWYulflri9VtGS0Fj+5ziPftl2Q2vTjiRSz8\neN67g7Y/taLTc5zMpldFQmFnIuHGtfYmtJbvlAA2QIzaYYgHvLDN28wNR362nLYqUWuml4DVkZsN\n91VR733AEKRn6RXtaCjucXoFkK61B5wGXnuOFzjIe/UwUum71NfaKc2OF2yAeMj7IQVg7ql0rBJp\nU2/Ug/lXiGE7FNikaPXHzDVsHxkUjJnzueF4wVpImGRvYFDou+OrITSZTXcDjjtm60P/ue0qA5vq\nM7kEmu57BMDxgr0R42tZ5PPzsNB3X8jn8/lhJ9x1JtIb9RNgcOi7ExootRQaOsdR3vi9wPeI8flh\nB7uUUt/ibOAvyHW2V5XacXWGpjyPOyBumuOmF1RzPYibXoin5qqgxmwJJLPphZEQwL5I7ux3Ratb\nM1ocIAA+AwaGvvtllGt6P7AOcAvSr7WcsNVqUpU5TmbTRwIZpBLksBobW6X8EFkW2BMxbndlbp7t\nd0hIdAg8EPpumzcQIq/zHcCvgbNzqczJ1dLreEEhPPsBpNpmu/MVVdR+AmnD0laYezVppDGbQM6j\nPYB9Qt8dU+KuFWt2vGBVpNfk8hQVoKoxJeuN+vMOQYzb3ZC2XSA3kJ5GjNv7gNfPXfOgPZDz+xOk\nINS3849YOslsuidwLZIvDnBALpW5tTNj1pmm+h6Jaif8C8mTn46E1Z8f+m7hO6AQefIX4GyksNcu\noe++2Qi9JdLIz4v1kLDg3ojh/0wJu5XyHdIduY6GAueGvvuXzmrtJE11HpdI3DTHTS+o5noQN70Q\nT81VQY3ZEklm04Xep+fkUplRRavm0Rt5ep5FikZtG/ru+GQ2vSFiyC4P+MD/5VKZRrYAqJYx2w0J\nNx4M/C6XytzU2THboSzNUWGVwYhh6yBzD9LH9DHkB0vYMpwvmU0fj7xHTwC75FKZmdXQ63jBGsAr\n0fE3CH23pMJZyWx6OSTceVXgyFwq06Y3two02tOyFhKC/RWwTui7U9raNplNbw14J2131D6bL79h\n2ZodL1gaMQjXAo4PfffCCmWXS0VzHFXO3pK5XtstilZ/Btz/9ICNFn5xsfX2m9a916PAkNB3Kzp3\no1SIHHJjaNxCPftuMWXG1OlI/vbjlYzZAJrieyS6SfMH4J9I5fKnkf7Ab7fYdI5exwuOBS4CvkWi\nFF6un+KyaFRawpJIHYPVgQND3725xF1L0huFgb+AfDakQt/NVaq1CjTFeVwmcdMcN72gmutB3PRC\nPDVXhYqNWWPMXxAjoSdyx3kscB2SA/Q6MNJamzfGjAAOR7wJZ1prOwrja8o3I5lN90W8s/2Q3NmC\n56P4R8iSyJfgakAy9N3RyWx6B8RL2x/wcqnMBXUXPz9Vm+NkNr0aEi43HQk3/qIa47ZCZzxwCWAz\n5ubZblq0+lWiPNs+WzzQM9Et/zjiyd20k6+l+LzojlSt3ZbyfnwBkMym10aur8WR9k13dkJXezT8\n2isK8zsz9N2WIf0ks+mlkPD9QwESiQT5fH5kLpX5dxnH6Ifc0Nic+ntfqjLHY93hSyPe2qGI93YJ\ngNkk8p/1WSrxQ89+T2/043vHAK8MCsaUfOMsCrEPkXP1fmCfU3f8009/f/yiGcAUpAdysxTHaY9m\nOJfXRHpl7wT8CJwEXNFGL9OWN79GIOkFk4Dd69BLuBLqPseOF/RGbqBuR/mtdErWG3l+n0f6XW8T\n+u4r5WqtEg0/jysgbprjphdUcz2Im16Ip+aqUJExa4zZETjeWjvMGLMw8H9IfpdvrX3SGJNBQimf\nQ/qSbo6E6D4NbGGtbS+8tmnfjGQ2fQxwMXBeLpU5KVpcXJjoQaR1yOmh756WzKb3QfqkJoCDmyhM\nr6pznMym00h15hBwaxRuXDXNjhesBOyFGLeDgV70mE6fDcbOpue0xMwvVj955qdrX9yeZ7AcvY4X\neIhnZgywbyV9cpPZ9EDEAOsB7JZLZZ7shLa2aPi153jBwki+81LA+qHvvgdzimIdgfR9Xgy5gXLx\nor0XuXrytJ8A/gGM6ijiIfoxfC/yvl+DeMnqGZ5S9Tke6w7vjnzGDp1NYs8E+eIE16+YWyH5oUHB\nmDYrfiez6WWibTcBskhxt+lAPplN/xb5LPsE2LoBPZDLpZEhsD2A44HTgT7I5+JRoe9+2s5uraUl\nHAhcj9xE2DP03Vpc852hrnMc3ZS8FjgYiRzYv40bA21RbnTPr4E7gQ+BLULf/a79PWpCwz+TKyBu\nmuOmF1RzPYibXoin5qpQqTF7NjJp6yPVN08EAmvtitH6YYjX4AFgD2ttOlp+B3C2tXZcO8M37ZsR\nFeV5D/kxvVoulfkayDte0A3JhTwcMViSfbe8/yjgEqQa7m9yqczDbQzbCKptzHZDCvfsBByYS2XK\n8jyWSE3OC8cLFiExe7fe6z13SbeFJ68w45O1mPnFGiA9QB9CfojeHfpuuV7awk2O9ZFiUpOQ8OJv\nKtWazKZ3Q3J/fwa2y6Uyr1U6Vhs0xbXneEEKqeYchr47LAopvgzxrk8GTgUuy6UyM7/86Zv8Mfec\n+i4SEngzcGhbueiRhzwLDAf+i9xYqDSMvFJqPscnHn7x1t2Y/ejqUz7vte6PH0zqTn7xaNVsxNtU\nMG5fKnhto4JuDyFtYi4Hji6qoJ0HEsls+s/AOUi4/Pa5VKa9yrGNplEhsJsiN0k2RfqH/xHIlXDD\npK1Cgvsg1e9nAL8OfffB6iruFPU2Zgvn3wvAjqHvTi1ziLL1Ol5wGtIL/hHEQ77AfV7UgLhpjpte\nUM31IG56IZ6aq0KlxuxVQMG7tTryg38Ra+0K0fqdkDDA+4ENrbV/jpZfD9xgrX2kneGbLom3mPvf\nfZz/jM+yl9mFgzYZDsDdT7/PFXe+xurL9+fckYP47zv3cudb99O/z6KM2v5oVhuwUoNV156vf/oW\n74Ez6dGtOxfsfioD+vZvtKSSuf2Ne8m9HrLJsuvjrrQf4976mhfe/JKPv5xbJ2qtlRZjy/WXZav1\nl+X/2TvvOCfK7Y1/k+303hFEZEUpKkhx7YqAGkZdScSrXFQUI/ZYr73XWK5eY+9XTXTVMerFclXQ\ntWFX1KBiAem9LVvz++O8C2FJ7/O783w++cAmM5MnbyYz73nPc57Tv2c7LJbY14uGxiYu/Occfl20\njstPHsWYIT1T5vrB759x76eP07GsPTccehFdW3dO+Zj5hmAwyOWej/juj8WMOnQN362RksED+o3m\nxOHH0KHFubW+diO3fnA/P6/6jSHdyrmwYgatist2OOb9Vd8y6+PfGbJLZ649bSzFRQXZ+khZx8ff\nLeGmJz6jR6cybji6H3XzvmPtV1+z/qcANEkyq7BdOzrutSdNg3fm/vVzWMomjt19Ao4hk8Ke38Fg\nkEe/eJ63fp3D0O67cdn+MyksKMz2R8tL1NY38tybP/Hy7F9pagpyyMi+nDppCO1aF6d87Lk/LOXm\nJ+cSDMKlU0cyOg3XEaPho28Xc/OTc+nSvhT3eQfSqV1pVt63qSnITU98xqfzlnL0gbtw6qQhWXlf\nEyZMmPgfQVqC72SD2ZuBFYFA4E719zfALoFAoI36WwMOQ2S3EwKBwEz1/EtI3Wy01hF5vbJg9zpL\nkD50nYEB/9jr2iVXPvhxI7DKUrJpbOnwD64EpgE/AxN8Ds8OfUTzABkZY7vX2dzz71Xg6DTLjTPF\n+VAkI7UIqZPdKiWzufQBbGv7cwAi8QVpk+JXj/f9bq02HF+1qn8N8KTfrU1LI+dmk6oAsF9I/Xaq\nyIvfnt3rLGja1O5qS8nmKy2FDQSDfG+xcKbP4fkgzObNWcNWSGb2aESCfITP4dkq67S59OuAK4Gv\nkazOuix8lHDI2hjbXPr1SKuXNxGZamO1VtkRuTY397bt1UxqU5l1YZuapseRRcjPKvSq7TKzsFXu\n/RLym3gKmJaFllHJIJvjfBDwMJLV/h043e/W3k7wMFH52lz6och1tRipu/cmRTa9yMoY21z6CKR9\nVhNiqvh1kodKiq/NpbdDFA27ASf53VomjQ5bIi+uyQnCaJyNxhdMztmA0fiCMTmnBckGs0cC5wYC\ngcPLy8t7Ic6vPwB3BgKB2eXl5Q8gspw5SKCwD1I79Akw3Kg1szaX3h7YqXjXL2YUdFwxs3Ft17nF\nC8fss7Gmvs5SunFi6bAPL0QMWeYCR/ocnqQlpRlGpgJDK/K9HwT8zefwPJvGw6eds93XShqUAAAg\nAElEQVTr7AV8hbQ8OcDn8HwSaVvlcDkBCW6PQKTmABuRYOFV4I3mlii/LFwbPP/u2Q1Ii42h6Q6e\n7F7nbYi8/zPgEJ/DsykNh835b8/udY5G6q/3DjZZ6+oXDipuXN73Cv8dx94YYZeWgdY/gTORxYmJ\nPofne5tLb651/xXpw7os4x8kMrIZZBUg5+URwM1+txbqwo7juTMO67yu4dV+i+tKh/9cE2i7uWkg\n2xZsViOLkbOG33X7E9+cf9FgpHZz8+eDW/HRnq1fD1oso4DrfQ5PIiY82ULGx1ldE24HpiOB1t3A\nVX63lsxvMSZfm0vfD6n3bg2c4ndrTybxPulENsa4N3KN64nIrF9N4XCpmAiWKx7FpLGXcxzI+TU5\nCRiNs9H4gsk5GzAaXzAm57QgFTfjW5EaSSviPvo7sjpdjAS2pyk34+lILakVuDEQCMRyYs1VrVMB\nkqXYKeTRr8Xfom+0NFEybA6WojrqvjuQYMGWs0qHfDwVaZ3xH8Duc3g2ZvszJICMjbHd6xyAOATX\nIu7GS9N06HTX+RYC7yKumOf5HJ574t3X5tKLgAq2ZW0HqpeakDY6/r7d2966cNkGgHF+t5b2emm1\ncPAEcBJyzmk+h6c+xcPm7EKoXIpvRtqYADzdsLLX9fULhn0ItEFa9SwMs+t2nFWf4IsQQ6h19YsG\n3teweODlwFJkEpprpUS26ww7IotruxDSv9fudVYi/a4BHD6H55VqrbIdcCjb2v9ErY9osBJsKLRY\nrE2sLW4IrgBqUAFvEv+P+FqFXpXMeZ3Rcba59GMQFUpPRAkw3e/WPkvhkHHxtbn0fZDFs46A0+/W\nHkjhPVNFpse4NbIgvjdwod+tuVM8ZEp8bS79SESNswgxhFqeIp94YMTJqdE4G40vmJyzAaPxBWNy\nTgv+Z/rM2lx6G3YMTkP/7oP0hg2HDcAf6vEn8GfRgG8GFHZZctqInnvyxZKvmw1ongROS0NQkWlk\n9IS3e50zkXZNOmJ+lY6TLN3B7C1Iq4wqYHKyHJXDZjnb2v7siyzcAPzL79bOSgPdsLB7nUVI5m0C\n6ZF8Zv1CqLKppyMuxR2RwOCsZrdmm0ufhjiY+vxuzRHmEGE5273OE4JBniRoKaz/bcjmxlW9x/rd\n2reZ+hwJIOtjrHpfN6sORpeNmjUGaRmzGZjkc3jea7lPtVZpAQYDE3rajnIv8b/2INBKPcqAVvUF\ndFzfpmBIQSMFbTY3ri5sokC9lnqh6PZoJLHAeGP5JRf+I3DrHd0q9Kq0qmNsLr0ncm07FmlHdh1w\nu9+tRVMbRYXd69zrjH1O/PKBuc8UhphuReMwDDHc6wq4/G4tV+3eMnYuK1PFF4FjSJ/reMp8bS79\ncuAGRI02zu/WMn2vN+Lk1GicjcYXTM7ZgNH4gjE5pwX/L4JZdePrQfRgtWOE3ZuAxUiQujVYDf07\nnERUBRIBpKcsSFbp8jytH2uJTAezViTreSBwQppaEqWNs93rPApZYf8FGOlzeNImAVa9ho84flz5\nk8+/HWidYnufmLB7nW0Qafcotm8ZlQyyeiEMlRTTwqW4eRv1264GxgCH+N1ay8ArkgvsGGu7le8V\nD/y61FLYANI+7I48+H3mSnkyGfAV9vplZVGfX7ogvZQn+ByeaM7yzYjIWTlNv6u2OdDn8HxerVUW\nIkFtGS0C4DD/j/ZaPNvFwg9I4DEHmF2hVyXVO1otWp2KtNhqj7SZO83v1n5K5njNsHudNqTNTCmy\n+Pe3eEoGbC59N+R33wupBb8xyy2mILPB7M3ApUh/7vGpLBaEIB3BrAV4AXFEv9fv1s5JA69oMOLk\n1GicjcYXTM7ZgNH4gjE5pwWGCGaV3KgvkYPVPkBRhONtYvsgtWXA+leyq6t2r/O4ssLSF2oatpzj\nc3juTeYYOULGT3glN/4OaXGTDrlxWjirFiRfIhPh0T6H55tUjxkBWbuo2L3OLkjANwi4wOfw3JXk\nobLCWfG9Gak1BHgauDjSOaIMYOYigcleLX6v4a4XuyOGMe0Luv15XnH/Hy4FeiMZtfPiyX5lEDm5\n2di9Tkvj6u4fF3RaNjpYX7QFa+OIF0741w9x7h6Vs93r1BBTqJXA2GyZ3qnscSk7Brmdd/rblP/+\n+e/n3kbKAUKD3p+R4HY2EtyGk65vB5tLH4hksg9GVDqXAA8m2ON0B9i9zpOR0pzaXTvv3OrnVb+B\ntPGaFE8fX5tL3wUJaPuhFlON3jMZwObS/46UUPwMjPG7tYi9kRNEWvgqldfHwBCkdvnxVI8ZBUac\nnBqNs9H4gsk5GzAaXzAm57Qg74LZl9//JfiYf94/2T5YjdR/JIgY7OyQTQ15rMnkDb4p2BS0WqxG\nO3myFbScBdyL9PQ8NtcSWOVE/SEwEpGDP5LK8WIg21nO/ki9bk+SN9/KdMa+ADgNuIkwkuJosLn0\nBxE58vl+t3Z3yEvbcba59J2QcegNTPO7tSftXmcfpK54CPAyMj6J9qhMF3Il5fYApzXVlm6u+3F0\nq2Bd2bV+t3ZNnIeIyTmktGA+sG+oK3iOEAQs1VplETACUYkcCOwHtA3Z7jdCglvg9wq9Kghgc+mF\nwAXAtUjQ7AfO9Lu1RaQAVdd9CRKArgaOfHbyfR+f8MJZjyLZ30XAUfEstNlcel8koN0VMT87L4sB\nbdrPZZtL3x/5PBuRQHZ+Gg+fNr5qIeFzZKHkAL9b+zQdxw0DI05OjcbZaHzB5JwNGI0vGJNzWpB3\nwazNpYcSqiGC9JdtWdVwbVGyCSOePNkKZq3Ae0hbmyk+h+f5FA6XjmD2XuAsstNSJBdByzBETtkK\nmQy/leAhMmkMNhoxzBmBZLeupIWkOBqUfHs+Utc+KMSReCtnm0vvimRky4GL/G7tjpD374BkDw9G\nsiqT0tjSKBFke5GjBHgGOA74smFpvxPq/xw8C+gPaHG6w8bF2e513orIuT8CDsvhggFE4Kwk0Huy\nLbjdn23O5AALgdmLSrsumNV1zHErizvsjsWyAjgbqdtO6Zqhrolu4Dz1XuN9Ds+PQFC91mxethEx\n5noj1jFtLr0HUkO7B5LpdfrdWjbUB2k9l1WA+Cki4x7vd2vvpuvYCunmOx5xl14KjPC7tXSZHYbC\nnF9kHkbjCybnbMBofMGYnNOCvAtmP/5ucfCmJ+aOQILVVTmoA0oURjx5simB3QXJwG1G5MbJtkVJ\nibPd63QAzwPzEHlxOlrZREOu5KQHIm6nDcBBcdZENiPtnMNIip9BJMUJ1y7aXLoTqbF9wu/WTlZP\nBwGLzaW3RWo3RwK3+d3aDrXDKrB7HJiCBMYTc9AHOpu/vTZIAD8OyTpO8jk8620ufU8k4KwHRvnd\nWiDGoeINZq3I9ztFva89h5LuuDhXa5VWYCgquA3CARbo0vz6FmtxTUGwcVZRsPFtZAx/bM7cJgq7\n11mMyGenIJL58SG9kLfyVU7TzyBGWuf4HJ5/xTq2Wux5C9gL6bc8ze/W4looSgHpzHR2QBaZdkP6\n9D6cjuO2QCYyyRcBtyG/p4PTVNsbCnN+kXkYjS8YiLPq/uB47IrDn+7ascwQnBUMM8YhMCLntCDv\nglmM92UYjS9kPzvU3OfzJeC4JDOiSXO2e53liCQMYB+fw5OSaUucyNl5Yfc6j0VMSlYBFT6H5+c4\nd02nyVYBEsDejEiKvwdmxiMpjgTVPutzJLO2r9+tfQwEbS69FHgNOAwJVk+NtAimAq6bEJnncqQf\ndCIBf6rIliqiE5I1Go04Xh8fmim1ufQTkKDnR2C0361tiHK4uDmrBYNZSK/pe3wOz3lJfYDUkfA4\n21z6QZZg8OHOdesG7rz5r1Uj1/30ffuGTbsB3UM2W4Eyk1KP7yv0qpi1s3avsy3inD4OCXxsPocn\ntBZ0O75KyfAq0A25drpiLQyogPA/iFlaFXBCOoKraq2yGPGl6Bvy6N1v6okz/3jqmaHAvGQDfMW7\nEDlXxwF3+t2aK1XOEZCJYNaCtLg6HqmjPiOdx8dA8wsln594/1E3vt6ldSdDcLa59IJHrxjX0K1j\nK0PwDYEhzgvVJ9oLVHRoW8LaDbV7+93aV7nmFScMMcYtYETOaYEZzKYOo/GF7AezVsSVcn9kUu1N\n4jBJcbZ7na0Q6doQ0uesHA9yel7Yvc4zkDrJ35AaxngkcOky2RqFZFCbJcXNLsUpt7GwufQKpO75\nC2D0K7dPajj6old9gB2Z/FfGk5FSNZ73IooBezxyzjQhG+ZrvZHs/B6IpP7UcHJum0u/C5G7voz0\noI0UlCXEWUm6P1Tvn4ohWSqIm7MKAm9HFl+agLuBq/xubZMymRrENlnygUg9djNWI9L25uD2mwq9\narugU/VQfgNRDbyGyIdbupzvwFfVwb+GjKMfuX5F7V+uFAqvIosJbyC/hy2Rtq/WKguQOvu+YR47\nqX+7R9pf4U/gdfV4r0KvitvBXQWD/wKcyGc9OoMS6UwZVrVCDPj2BM7wu7UH03h4Q8wv7F5nN6SF\n0lFtiluzsW7TwT6H5/0c04oKm0sfh/zuhwOPAOdmuvtAOmD3OkuerrxnS0lhcV6fFzaXfijwHNJC\nrNpioSIYZAPyG093CUEmYIjfXgsYkXNaYAazqcNofCE39ZwDgW+R4GF3n8OTaMP5ZIPZx4FpgMfn\n8JyZ6P4pIOfnhd3rvAa4GvgaaZuyPsYuqUq5uyBZz+nqOElLiqPB5tKfAk4Czjhi3/4PvPHR7yAB\nxXi/W4u7TtPudR6N3GyLgBk+h+fRdPKMgEybbA0E3kZqYu9BgsmwQaqSf72FBD6X+93aTREOmzBn\nu9e5EyIb7YUEb75E9k8D4uJsc+nHIMFUT6QcYrrfrX0WaXsV3O7M9sFt/5BN1iGB/Gxg9qNa51Ub\nWxf8BzFoegw5z8IttoTla/c62yOtew4HvkIyun/F+EytgJcIBse3a9j0wfGL3/5Hp/oNXQkfsPYi\ncn/1WsSMamGYx+Jdzzv7y5/vvtcLjGdb3fEWRO7/OvB6hV71RwyuZyPGVd8C+8VQCKSKjP32bC69\nP6IaaYfIjavTdOic30diwe51TkDk892BuQUW6z6NwaYGwJlhk8WkYHPpQxFp+AQg2KVDmWXl2hqQ\nEiSH363NyyW/aLB7ncOBV1oXlfXfVF/jBB7OsTv/DlDt9P6B9OBuQEz0/nXxSSObbnv682alyEl+\nt5bte0KiyPvfXhgYkXNaYAazqcNofCF39ZznIlmPKmBygnLjZCbUJyMTyC8QuW02zcJyfl4o2dcD\niAvwe0iNaLQxSHbBoFlSfBPQiTRIiqPB5tJ7Ij2ey4BC4BvgIL9bW5vosexe51gk69UZufleY1Rj\nMDXReROZVF4F3BDrs9hcejfk99EbOMLv1maF2SzZ82I4sshQghhCfZDoMVJAVM7qHLoPOBaoQ777\n25OR5VZrlTuxfXA7sPm1ukJLcHHXIktNiXV2+R9bLrcGmVuhV4V7j4h8VU/ze4EZwF/AUec+u/wb\nJIAMzaBufQShbxOW/gUErRFoN6pjhQtUmx8rYsiHg4hjdCGwL3CkeuwRss08tmVtP6rQq0J7SE9E\nsrErEKl71MA3Fdi9zs7OfU5a6Zn7dCefw7MmE+9hc+mHIItDKxFDqKiLDnEi5/eRSLB7naXALcC5\nSP39ZcBd1xx8fuM17921GrkX3AVclA8Bl82l90J+5ycDVsQ1+6KqW476svLS1+5DzCFrkM/zSL75\ntSjfj8eBspLCEmobagE+A870OTxf5JScgs2ld0Za7U1EriF2v1v7RL0cVNnaV4A2SCY8n1ta5u1v\nLwqMyDktMIPZ1GE0vpC7YDZUbpxotiZRqeMwRF68Bdjb5/D8lsB7pQN5cV6oQPMF4Bj175QoE4tk\nFgxGIZmtkaRZUhwNNpfuAu7o0bkVS1dt7pmKk6jd6xyE1HnujEwWZmSQf0bOC7vXuR8SGLRD2h3d\nH+++Npe+DxJ01gAj/W7t1xabpFKvPg6Ru25AFpR+TOY4SSAsZyVrPQW4AwkGPwRO87u1tNXRV2uV\nvX7tUzyjpsR6ac8V9cWd12/3c6tBMtbNsuRPK/SqLS35VmuVbWgRnC7tXDixrsgyqu2mpmCHjY21\n1iClESgEgWVBWLSotFuPpSWd+9RaixYO3fDr1e0bNv2ITDKXtpRDJ4GwY1ytVfZjW2B7CGzluRZZ\nbHn95R4H/BZo0/8NxOTqwAy2tsHudR6HXKO6Ic7DZ/ocnpcz8V42l968YPsZ8rkiSrzjRF7cR1rC\n7nXugdQKDwN+Aqac++zyLcDxg1znXX3xr0+Xb2pVoCOGXm8g951YyqCMQPUFvgi4EHH6n6f+nqUC\n1mYTwWOQxe8OSJ3n6X63lhPOoVD38OuRxYKNwEkPTrrl5RmvXvoscAJSGnE/cKXP4Ul4QTddsLn0\n0cgcoy/yOz/R79ZCOwY0j/NeSG1/d2QB/Ip8WzhQyMvfXgwYkXNaYAazqcNofCG35kS7Ipm0TYi7\ncbxy40RMaNohkq9dAc3n8MTTfiTdyJvzQq2gv4UsIvwLODtCxi6RMc6KpDgSlGnMCY9fefiTXTqk\n7pBo9zq7I8HgSORGPNnn8GRC8pj288LudR4BvIjIpacmUxduc+nTkED+O2Cs362Fun2nKj+fCjyJ\ntFUbm6VzZAfONpc+EHgIac+0ATEBezBKrXBSUPL155Hsz9/PfXb5f5H2ZM2Z26Ehm9cCn3Yfd9gB\ny95+5w22BbAdiICaYgsbWhdQ2Bj8vtP6xjnsmFH9qzn7q0zTHkayUd8B40LaWqWKmOdFtVbZChnv\n5uB2p+YdF5d0od5a6O1fs/RW4OtUTKTCwe519kCud8cCtQf0G10y549PaxGlwIvIdTCt7XTUYskT\nwFQkMJqe4kQ9b+4jsFXtMxOpNS21NgYfnP7yytlldcHpyMJFM37ZXGK569Gju2hNBZbDkQDSls1F\nZXWPOBkJBLsjCxlXIm74oVL/rWNsc+n9kCB9X2ABIjvOpkHgdlD+A/8GjgB+AY72OTzzUJztXuch\nSCBbDiwDXMCzGVYXbQd1zs8E7kSUUlcDN4a5roaO8wDkPjsQVX6RBff1RJFXv704YUTOaYEZzKYO\no/GF3JsTnYfIj170OTyT49wtLs7qZusFJgO3+xyei5Mmmhry6rxQN8U5yET6Cp/Dc2OYzWJyzrak\nOA6kbZxVKxsvMnH4CnE6Tnfgldbzwu51TkFMnuoRp/CkjaxsLv1fwJnIGEwJmYSnzNnudV4O3ICM\n64EZWigIRejEqRCp27oWyRK+Bpzpd2sL0/2mdq/zNETaXwMcG67Xc7VW2RlZWGoObvdk2/iuJ7r0\nd9E9J3TbAzF56oFIpc+P1q9Z1bD9E5lwBoBDcyGBrdYqLWsL2+z1Y5v+L/XZsrxfny3LmywS8AMs\nRjJ4rwPvVOhVUY2uokHdA6Yi95iOSPb9VJ/DE7B7nYMRs58KYA1yXjyZzsm/zaWXIdfakcBZfrcW\ns7VSFOTNfUSZPD0GHNl6c+PqidXr3+y9ov4AthmivQc81n38uKeXvfl2PVAUhN+/HNxqwcfDWh/S\nWGBZCRzjc3g+zCRPFVxNRALu3RGfjtuBO/xuLdx5td0YKy+Ba4FLkZrPS4C7s509VOfqK4gB3VuI\neWazRH4rZ+UgfwESqJch6rczs6GCUYZzDwMORF4/xe/W3omwectx7ob83rca4+WZAVfe/PYSgBE5\npwVmMJs6jMYXch/MWhGJ3X6Ik+wLcewWbzDbbCjyIXBIpuWuUZB354Xd6+yFSBx3AqaHMTyKytnu\nde6DrAJnVVIcA+kODgsRF+jpSCZxYponBekMvs9Egpn1wFGpThJtLr0YMe+pAC7yu7U71EvpCGYt\nwIPAaYike1KGz5tQSdujSP/VFcDZgC/dE1P1+S5HskArgSN8Ds/cePat1io77HXfPWu+Ouvc9hV6\nVVyyRmWw9RqyOPUGMtGNuECgJvi3IVLLBUhA+3s87xUFCZ0XisPTwN+A587/9dmZJcGG8UjGdiJS\ntw5SwzybbSZSv8T7HmpcHkTMfTYhgYhHmaA1Z7OswBnArUjt3tvA6T6H5/d43ycWbC69L6IO6gQc\n5ndrs5M8VF7cR+xe5wSCwSd6rajvPur7Tct2Wlrf2SJZuA2I6sJToVf9oDYPVmuVfYGLkd97aV2h\nZe3Hw1q3+35gWUNDoeV0n8PzZCZ4qt/7HUiWuAkJvq/2u7XFUXaLVJIwDlEcNQdd01rIZjMGu9dp\nQzKybZFA/LIW5UE7cFbu5/8EbMjiphvxTQhV2aQNNpe+B+J/Uo60G3P43dqiKLuEU8u0Ucc4HJmb\n2PxubVUm+CaBvPjtJQgjck4Lkg5my8vLv0ScG0Fujjcj8pomVLYmEAgEy8vLT0MMaBqAGwKBwOsx\nDm20L8NofCEPOCu58bdIDcjuPodnRYxd4skajkKC2LXAXrFcPzOMnI9xONi9zt2QMeqIrJKHSrDD\ncg4jKf43YuqRFUlxDGRCtmsBrkDMQtYggVe6sgnpCgyb+S0Dxvscnm/SwA2bS+8BfInI8sarVfa0\njLFaKNCRzPdjyIJKRlZTa+sbg8dd+tqtSPBWgEy4XZmYKCm1wj1I5vMP4HCfwzM/wcMkU6/eDsmi\nT0BKN47yOTwRJ5MqmLxaPRYhAW2iPEORaDB7BRLsf4I4/m6tJ1UtgkaxTY68Z8iu89lmIvVBOAOt\nMAHqW0iAGmoqtR3fMIHvP5DFubSYFdlc+gGIydAapBb9zyQOk+uF59KyLU13Dly4xTlsfk2wy7rG\nZi7fIxLuZ8Jk0bdyrtYqeyDSVyfQenOJJfjF4FaWHweU3VlTar0oktN6olCLBzcCJ6r3/g9wsd+t\nfR/H7hHHWF0Pn0Z6mP+F9G7OmAopZFHsOqQE4VSfw/NsmE2jGcZNQoLafkjbrHN8Do+eTp42l34i\n8ttphciLL/W7tViLk5EWDYqR+8HfkL7nE5L8raQbeTmHiwEjck4Lkgpmy8vLS4GPAoHA3iHPvQrc\nEQgE5pSXl3sQPfwnyE1lBCJ/+BAYGQgEojlGGu3LMBpfyBPOdq/zfORC+ILP4bHH2DxW1rATMgnf\nCZlMRpK6ZAt5McbhYPc6RyMZOCviMtvcRqLlZK8AOBVZqGqWFJ/lc3iSzTJkAhkbZ7vXOQ2RUDUC\nf/M5PFVpOGyq9adW5DdzLvA7MM7n8MSduYoHNpc+FsmMrUcm4b+RXin3e0h2/xqfw3NtqsdUdaEl\nIY9hPbu0fnPJyk0gYzTD79Z2kPumA0ri9xTS5/g7YILP4YmWBYqEpM4LtUDwTyRYWILUJUZ1NrW5\n9IuRoG8ZkjWMZ8IfDnFztrn0yUiLoT+BUbHqdqu1yt7IoseRwDhk0gySCXwbyUr/p0KvWqoWRh9B\n6pLXAucTXjocLptlQSbR9yDXuE+QAOIH0gCbSz8TCfq+RFoPxd06TCFn95F7LjtpUquapsd2WVjb\nuaQhSBAaLFJrfD/wYZQa5x04V2uVXYDzg3CuBVrXFFuY37/0xyAcNvP2Z5L5vQBgc+ntETnweUgJ\nwTfAhVGkrnHxbfEeViTDf73a7lqkJjStDs3q2vgEUIn8To7xOTxfJsPZ7nW2RoLiCxEvBT9wbqo1\nyzaXXooYnM1A7g8n+93aS3HuHm3RwIooR1zIosGEFK5LKUFdEw6ZMlR757nv9AOAj6OVceQZ8nbe\nmWkkG8yOZpuhRyHyo3kxEAj0Ua9PQmQDbwJHBAIBp3r+JeCmQCAQraDeaF+G0fhCnnBWwdJsRNYY\nS24cbSXSitSQHUmaJshpQF6McSTYvc6JyA1uPbB/qKmEer2lpPhq4L4cS4rDIaPjbPc6xyMTuNZI\nbeI9KR4yab6qRcsjSD3gPCQjmxH1gc2ln46svH/14i1H7VVSVNBSHlaITB5Dg8i4/raUre9cMviz\n8yyFDZ3qF+/8esOi8vmJ7B/m/4Ut+Vst0BTkTuCqFmZWaYPKjL6MSBo/QDL4ybqJpnJeWJDFjTuR\nWt0TYmVibC79LKTdzyrgcL9bizRpjoa4OCu37DmIOmtfv1v7LpE3qdYqS5Ha4uas7YDm1zaWWRfN\n26W0x2+9SwqXdyx8OWi1zIyiGIl2D+mGBLTHIzLnG4BbfQ5Pwq2aQqGy4Q8ji4LPAFMTlLhn9T5S\nrVUWBeHoDa2tN7bb1LQrwJZiy8aCxuCdRY14KvSqeAyzInKu1io7bimyXGQNclFxQ7CwrtDSVFdk\nua9NTdN1FXpV3KoJVdd6OnAN0AUJgC4HnkkiyIz3PK5A+pL3RepS/xZDvhw37F7nAKQ+dijyW5kc\nwxwz3rKrwchiysHIteEGwJ1Mm0KbS98ZuRfujajqjvO7tZ8TOERMzs1dCpBFqUl+t5a1dm7KKPME\nZGEk1KRvDcqFHZjlc3iyIjVPEnk978wkkg1mhwCjA4HAo+Xl5bsiNVAlIcHswUj7g1nA0EAgcKl6\n/kngqUAg8N8ohzfal2E0vpBHnFVblG+QgGmPKHLjaBORS5Hs4dtIjWPOe9qRR2McCSEus4uAfX0O\nz58hkuLTyD9JcThkfJztXufeyI2sB1KHdHEK8rhkM3BliKTUhrScOsLn8KxOkkNcsLn0h4Hp7dsU\ns25j3VK2DyAj9S+NC5bSjZTs/ikUNFA3f2+a1nWNtnkd0mKrNuQR7e8N7nMPOHXQTh0zucjRHalV\n3RuZhJ7gc3gSzbqFIh3y80nIZLsMyXDcHU3GbXPppyCLI+uBiX639nGCbxnP5LQv0qamGzI5jVVm\nFBXVWqUFKF/ZvuDUhkLLjK6rG9oWBLeSWWYReenrwNsVetW6FrvHU6oyCamZ74Vk2k+Nt/Y5Emwu\nvQQJfsYAF/jd2l0J7J6V+0i1VtkLOL0JzrBKiQGLuhXVr2lbcMfQX7dcFdofOA7E5PzcjOM7rmtT\n8MYuC2vHtKoN0mRhszXIfYC7Qq+KGMSpxQENURYMQuYNtyAGTcmaByWiMOiE1MII4hkAACAASURB\nVOAfjdTGT/W7tf8k+b4A2L3Ow5Breyck8Dw/jkXjuDmrxa4pyGJXd8QEbqbP4Yk2D98ONpduQxQo\nHRBJ8FmZUhkoCfPjyOLXFL9beyXB90kIyvXcqR5dESXWC+eMOfn4f37yuIcQF3bkM3zCtrKHb7Lp\nHB0H8n7emSkkG8wWA9ZAILBF/f0ZsFcgEChSf2tIjcFbwIRAIDBTPf8SUjcbbRU4n04ME1nAa4F3\neOrrKsb2HcH5+05PaN8fls/n2vfvpmNpe247/B+0K22bIZb/P/HqT2/xzDcv07tdD8YPPBDf96+x\nsW4Tfdv15NQRx7N7t0G5ppgXWL5pFTfPvo+/NixlbN8RzBz9d4oLirLy3pvra7jtAw8/rPiZYd0H\nc2HF6ZQWRWoxmj7UNzRy57NfEvhzDUUFVoqLCigstFJUYKWoUP4uUn8Xtvi7qNBKUZGVogL1XKGV\nYvV3ofr/0i2L+Pf8J7BarJy55xn079h36/bFRbJPYYEViyW/7s3LNq7ghtn3smzjCg4dsB/TRxxP\ngbUg17QAWLD6D279wMOaLes4fJcDOHlve1Rus79cxJ3PfUlxoZWrTh3D0IFd0salpraBS+77gN8W\nr2e6NgTtgF1SPmZDYwMv/TiLl3+cRWNTIwd234tJll3Z8vU81nzxJfXrJH61FBTQbvfBdBw5go4j\n96asd++4z6PNdTU8881LvLPgQywWC0cOOhTHEBslhcVJ8161robz75rNuo21XHv6WPYc1C3pY6UL\nwWCQ9d/PY8kbs1j1yafQ1ERdkZV5O5dQv+8Qpk1w0rGsfUbf/8157/DVS88w4ofNtK5pwlpcTPfx\n4+h9jEZJ587bbT//zzU85p/HvAWrsFotTBjTjymH70aHtiUZ4xiJ9+vVv/Hoq/NoaGzimIMGctLE\nwRQVJra+FwwGeX3+uzz9TRVWi5Xpex/PobvslyHWsKluM97v/bz5y2yCwSAVO41k6p7HRf2OGxub\nePo/P1L13i8UF1o549hhjBvdL2Mcm/HFT8u45cm51NU34qwczoSx/dP+Hr+vWcjr89+l+s/PaWhq\noHVxKw4bsB/jdz2QLq06bd0uGAyycN1ivlzyPV8t+Z7AygU0BWUtu1NZB/bqOYS9ew1haLfyrNyX\n/x8iLTf4ZIPZGcCwQCAws7y8vBdicrAAuC0QCMwuLy9/QD03B8mW7YOs6n8CDDdrZnOOvOKs5MZz\nkN5uk30Oz4thNgtX79Qd+BpZTTso05b/CSKvxjgS1KqtG6kzg/yWFIdD1sZZ1WXriAv3HKTn35ro\ne+2AhPgq+eN/kAzgi8CJyUjEUkQm65KPRT7XMmBMC8OeVJARznavc09EcdQdqaG7Ok0r8+l0ue6L\n1JQOQ7g6fA5PRJdkm0s/BskMNQLH+N3arDjfKloNXAHwEjAJkas7U3WQVqUPjwFDEDXJGT6HZ2um\nt1qrtCL+HM1y5JEhuy/oeeTEAUte/8/xSKb491h9be1e50FI5noXZH4z3efwvJcs/5Ba9A1sq0WP\nhbSfx9VaZTvgJKQN1+4AG1pZV342pHWXn/qX1jUUWi4B/pkt9Ynd6xxf0Bj07b6gpt2+X29aV1of\nbI+oMR4Dbrll4FQrohY6Xu2iI4ZDPyXJLyW+zVDOyV6kn/1nSBZxQTz7KqXNg8j3sBRp45WIMiKV\nsoS9EfXBKESVcSVwf8u6UJtL74koPQ5Eetwe53drqRgNJmoYtw+ifumCzEmuT8M1pAC5NpwPHKSe\nno/UAT/Vwvk5kiFmJ6DZhX0C27uwv4/K2vocnl9T4ZoobC590GNXHB7o2rEs7+edmUCywWwhIgNo\nXqK5GKm9eRgoBn4ATlNuxtOR2gYrcGMgEHg5xuENEQSEwGh8IQ85273OciQw3YC4G7esSwhnTvQW\nUq92kc/huYP8Qt6NcSSomuNbDx2w34X/XfBhrzyWFIdDVsdZ1dU8g5h0/IDI2hNxXkxEHrYTshg4\nCJlUn5EjCX2m65LPQeoVfwT2S5N8Ou2c7V7nwchEug1wts/hSaV/aEukla/d62wLPI+YKH2HOB1H\nPE9tLn0iEnwWAPY4pX3RgtnbgIuQRe2JcTidRoSa+F+H9NK0IkHAxdECdNjqpDsRmXQejrQ5acYK\nJAD5VP37WYVetcPClN3rbIXUZbrUez+s3jup2mibS5+ujvEtUj8cq547bedFtVY5BAlgT0LO4frN\nJZa3ZlW0H7ywe9EALJYfELl8qs7oCXNWtZ2vWRuDA0b8uPmLsd9u6mSBnZuwNH3fdgAfdRxqXVvc\nbi5i7pRuN+Gkx1j1Wf0XMqbrgel+txa11aDd6+yD1NqPRM69Y5PwPkiHkeBpSHlWR2Tu5fQ5PJ8A\n2Fz6Qcj1oztyXTjF79ZayvYTRcKcbS59EFKv2h/p331WMsZb6no4DfEWaJaIvIP0op4VYeEmnrKE\nAmA02xbQhoe8HGCbHPnDVOvvI0H1670ROHXILp0tN5+5nyHmnemG2Wc2dRiNL+QpZ7vX2Vz87/U5\nPMe3eLllMHs90p5ER1z/8u1EzssxjgGTcxxQEwE3YhSxBKlf/TrO3eM17tgNCWT7IPVhl+XwHM9G\nXfIdSMDwAeJGviXGLrGQ7uDwOKR+HCQ7Hk9v7ESQieC7EJmsnYVkf2w+hyei+aLNpR+MmMKVAif6\n3drzMd4iLGebSz8VWXyZD4zxu7VE1QtbYfc6D0BqFAcCvwKnJZMdrdYqi4fcdH3t9/+48iIkIzWK\nbYvxzfiZ7QPcryv0qlrFY6TiMQxYDJyZbLsTm0u/H6nP8yLZvGi/65TOi2qtshg4Bmkbtb96emEQ\nHnh2YkfLyo5FVyI18Pcji8LJ1p2GIllfgM5Iz9ED2VK8aOAnPTuNXjG/VZf6dQTlmM9a4MYKvSqd\nfb+T5hsKm0v/OzKGrZDFlvPD1ZTavc4K5DN2R5yLnUle69JyvVDKn1uRQI9gkIdrv91/SbC29RVI\nm82LkXrkqPeeaq2yEPlMvZF689B/ewPdO43aZ+jqz+ZehvSTnVuhV8V1rqkM8X+QQPElxHgrrjFT\nvXfPQtoKtkc8FZ5B/ARiuSUnsyjTh20u7IexvQv7W0hg+4bP4Ynq5h4PlAnamYi7dntg3h3n7L9H\neb9ORpvDpQVmMJs6jMYX8pSzWuX6ABgLHNeiFUqo0+4E5OL2G7B3Cg6imURejnEMmJwTQEhrqQ1A\npc/heTuO3eJZ7R2JyEM7A5f4HJ7bUuWaIrIRzFoRSZsdaeEyJcUelOmU7Z6BTFQ3ItLyd9Nx3BbI\npJT7HCSorUVaTEVUR9lc+r7ItbUtcKrfrT0e5dA7cFYZnbeRLNWYBN1OQzm3RSbZTmRSfTdwZYrB\n1nZ8VeZ2HySz0hzghhYQ1iMZq0+BzzaWWb94TOt8bNBquRJRoPmQHp4JTUxVX83/IuUKl/jdWrTf\nd1LnRbVW2QdRxJ2GGNeBTKbvf/KoTnPXtit8BMlarwRO9jk8ryX6HlGQrGzXQkH98UX9fniwsMuS\ntsG64mDjH4PunfnFB5+WNtVfgiwkBJFg8IYKvSotvbWT5dsSNpe+G7JAMQxRQzj8bm1r4G33Omcg\nDuJWRGVwbwoLlOlerNs/GLQ8YLEEdw/WF1G/aNe1jSt72y79+d/ViOlTaGAaNlgluingRkQN0IxG\nxPjz45DHb5Gk/6oF0yuINHgOoPndWth5nyqd2hdZbD5W8VqGZNAfiGIy2hKpZr8jurADn7Mta/tF\novc6m0s/FGnJtjvi/HwV4FEKGKPN4dICM5hNHUbjC3nMWcmNvwHWIe7GzXLjIGBR9WBfIZOtfWP1\nVMwh8naMo8DknCDsXqcdeBq5YZ7ic3iejrFLVL5KyvoqsqI7w+fwPJIurikgK2Osbv5vIRkkt8/h\nuTCFw6XMWU2KrlaP5YikPJkWNvEg01LuoxDZYCsk2+KONJG2ufQRyPfQCZH1RZJTb8fZ5tJ3RQK/\nNsA4v1ubnSTXCcBDSAuUH5Df1afJHKsFoo6xqrndFQlqmwPcPZE+nc1YW1tkmffjzqX9/uxZ3GdZ\np8I1m8sKzgWeSSQwsbn07sAXSEBwRJQ65bjPC+X0fAiSrdEQyfg6pCTMU6FXzbd7nUeqv7si8s1p\nPocnnnY7iSAZOWkFosoaA8H6on4/fFzQbeH+Fgu1wNRzn11ehTi5X8G2WuhXkaA2JbfpZPhGgs2l\nlyGqHSewGTi7ZNicZ6ylm+8BzkDK8expWBBLmXO1VlmCCkgXtOpVsbCs28WtSld0alewKth6S6Ol\n/cbGLW02N2ERpUYk1CJtkf5CFAuh/zb/f3GFXlVTu2pV8PNTTp+MJCvGIvXtoa5qyxFfnebg9vMK\nvWqrDF+5gj8DHIfq6R3aGkm1rpuMBLH7qKe/RhbyvEl4TaRzQdQClLMtsN2fbS3llrHNhf2tGP4G\n/ZHz61jF72HgCr9baw7QjTiHSwvMYDZ1GI0v5Dlnu9d5IXA78JzP4TlBPR20e50liIHGGETm5ckV\nxziQ12McASbnJKCkkDqygn05cHOUiW1Evnav82gk4LAgtWtV4bbLAbI2xspc40NgMHCuz+H5Z5KH\nSnVVvQBZyZ+BmP+M9zk8vyR7vDiQjez3XogxVC8kWDwrksmbzaUPRTKs3YGL/W7t9jCbbeVsc+kd\nkYnoICSj+1gS/DohE8+pSFuOm4Eb02h4lvAYqx63w9kW3I5GJM9bsb61lbVtC5a239D4aPtNTbOA\nL+ORUCqDmw+Q/p/7+N1auPMrJudqrbI9MmZnArupp79Czt/nK/SqTaru+DZEclkHpGryFA1xj7Na\nALkFmZwDvABc5ndrv6oFmOeQxZGrgBvOfXY5iPnOlUj2DSQov6FCr0rWADLtvz2bS68EHqWwtn3J\nbnNXWFtt7IrUSR/tc3jiMf6KhYic1aJMV2JnUzuH27/54JtLrWwqswL82ml9w/uFjfzBjoHrmlhG\napE4q2B6L2Q+1xzg9g3Zfofs7bO9Dv/jz1Y97kFk838C48tGzVqOqBDOUp8riNyP7wbm5Ev2OxR2\nr7M9MA4JbI9AWpeBXPc+YFvWNuBzeII2l968CHkJsrjwEXB2mP7gOZ8P5QpmMJs6jMYX8pyzmkh+\niFzkjlWyuKDd67wLcaF7DpHL5d3JG4K8HuMIMDknCbvXuTuyuroTYlRxdkt3SIWwfO1e5zSkJq8G\nmfC8kzm2CSPbJlv9kMCoO1Ju8FISh0mas8oQ/xuZYH+NZGTTnb1qiWxlv3sjAe2eSLA62efwhDV2\nsbn0ckQO2xsxQbquRe1cELCo2q1ZSFbwdr9buzgJXpVI8NUd+BLJxqZLRtqMtIxxtVbZCSVPriu0\nHBi0sH9JfXBr9jYIjRbJHIXW3/5YoVftYFxjc+nNvb7nAWP9bm1DvJyrtcphSAB7ItAaCVJ9yDh+\n2hxk2L3Oocg9cw8k0z3F5/B8m+znjwMxx9nm0rsgAaoTyVB9hJg7befoq7j7kfrm55HzokZloQ9G\nMrUHq83fRxzG30sgwIqLbzKo9Fxjs7Zd86KlqK64cU23jfV/DRz/6vUnfZTKMVVtar89rr/ml3lX\nXnMSOwarvYGebMv6hcN6YHETliV/lPXot7Sk84BNhaWbu9eucQ/d8OvrSJC67J4Tuh0K3IcYJi1G\nMp4vZjI4rNYqe7MtsA2bvQ3CJz+02bn0m659Dl/Zb/2Wpp5LLBYLJYiU+TFkkSYdTsLZuiZHdGEP\nBlkQ3NwuUP/XwBFN6zp3I1iwBDHWezZCHXNezIdyATOYTR1G4wsG4KwMcL5GZFK7uypOX+mufgjg\nJ2Afn8OzMZf84kDej3EYmJxTgN3r7IW0EhiOTMCmtLD6hzB87V7nBYh0aDViJpUOSWU6kfUxVlnE\nOcik7DCfw1Od4CGSNaFpj6zqHwi8hywsRHXNTROymf1uAzyLyDbnIU7Hv4fb1ubSd0YC2p2RzN6l\nIZOooM2lW5HFm9ORcatMxG3U7nX2QCbMlYhk8WpEAh1uIShVZGSM3ZefZFnapej8strgdV1X17fu\nu6x+Q5e1DcUWMVZqxkakTq45uP20Qq/6C8Dm0u9GXFZfRtqfhGZLW2azipGxmglUqKf/QL6DRyv0\nqq31gEraeA5Se1yCBLkX+RyeHYyJ0oyI42xz6aWK0+VAO8TU6xLgpUgmQ8qk6GUkE/sZ8pvc6rhf\nrVVWIEHtBPXUR8ANwKw4g9q0nxd2r/NE4OFgkJLGlb3m1P829ECw1COByD+jGSpVa5VFiHvvQPXY\nNeT/OxM5UG0gstR3678VetVGVd/7IrLA8Rkw2e/WdnA7Vwt7l6pHCVJ+cJbP4UmmFj4ZZUQJsvA2\nFhgbhLGWkOxtkwVWtC9iY3HJDz3Wb76r9Zamd4lSe5tpvumAuiZODNYXTcHaeIiloKkAINhkqQfe\ntliDOtL6J5wDdt7Mh7INM5hNHUbjCwbhbPc6L0ImUG+WFZWOr6nfshkY5XN45uWYWjwwxBi3gMk5\nRdi9znbIJGEcMkmw+Rye5SGbbOWrJpvXIxO7xYiLbz6e27m6qY9HpFbrkPr4QAK7J8zZ7nWGuma+\nCJyUBlfleJHt7HcBsoByLlKrNinSIorNpfdBAtpBSOB5rgq4gjaXfgFigvY1sL/frcW1yKjO/ZMQ\nKWBHoBo4NcHvOFFkui65O2LwM9naFKwbsKj24fEfrf+hsImRiDx5cIv3Xwx82oj1c73HAZN/b9Vj\nzzpr8ZV+t3ZDS87VWuVObDN0apYkzkKMyd5omfVVE+LHkQBvBZLRTKfJUzTsMM5q0WMK0i92J2Th\n7lrgAb9bi9myRJUYPYTIqRch5+tXodtUa5X7INdSTT31BRLUvlqhV0WTU6ezNrIQWTy4AMmAnuBz\neF63ufTxiLdCV+DVvjXLZvztrzfbsy1IDQ1a+yP1zi2xAun3+nNfx+SpC70vzGD7YHVljM8JgM2l\nH4+4jbdGztcLY30Hdq9zV7XteGTR6RbglgSvj6moZcoQBcJ5rTc37t5zZT39F9f91XtJE21qGnoX\nst3Hjlp7mw2+qUCVbFwLnImlqcDaYfnHxf1+/NFSXFuB1N0242tk8fz1IT/XfHbo3A1t9n35hTUW\nqzVv5kPZhBnMpg6j8QWDcFYTrmpkIgAwNQ6DnXyBIca4BUzOaYDd6yxGjBmmIpmHCSH1lkHEyKwA\nCQzOQCYo4yJlx/IAORtju9d5MiId+w0Ym4BzbEKc1WTtTSTz4UFk4tns6ZurBYOZiCtmHRK8vxhu\nO2VY9A4wBPk+Tr/y1NEN1z/6aRAxMNnH79YWxfmezVL8icAmJOtzf4bqN0ORLdng0UiQ2ROpkzzF\n5/B8Ua1VtkMkhM3OyaMReehWciuL22MNBt/pXL++Cvh0j2uv+nLe1de9AkxCTObWIOP/QIVeFbaG\nu4XJ0yzErTjTMvlQbDfOyuHaDeyNBEL3ADdHcqONBLUAcjFSS12DlBrt0A+5WqscjgS1xyke3yFB\nbVU4qXdLvslC1Xw/jyxk/tRxXcPkqa+vbkQFq1usxcNWFbWb1LpxS6d2DRsjWf8uQwWs6t+tjwq9\nKrQcIBmTrRLEZOssRClwqt+t+RL4fBZEFXA3Imf+FcnSRjIva4lkFxhnIvfJzkj22Yu01vkc4IRz\nnz2oXf0mvWftynZD1i/4tnvd6o6WyLW3zUHugjiyt1m9JttcekFBU+P04mD9TSWNdZ3aNNYsHLzx\n94dHrAssQMz4Om0qte68qcy6e0OBZWdrMNi5pC5oKa1roqQ2iBUo3msI+1xzbV7Nh7IFM5hNHUbj\nCwbirJqpv3/Ergd3m7a33RCcFQwzxiEwOacJ6sZ/HSJ/W4lIOT9lm5HZU4ADWV2dkI6+cxlETsfY\n7nVejdRsfgEcFGeJQdyc7V7nCCQj2xWRuV6fg3r8XC4YHIFMENsggeVt4T6/zaV3RgL+EcBrZSUF\nR9XUNm4BDvC7tZiOsqo2bAaitmmD1OyensVFnGxKuTsgJobTkdZCbuCalq2FVI3gaGDUFmvRodZg\ncGRxMKzC+gtEJuyNZC6lsle3I5P/OiTwuzcLiwQt0VxLPRjJUtrU8/8GLve7tT9SObhaLPg34sx9\nGXBruPO1WqscDPwDOAFZCPgJyQw/V6FXhQ5y0ueFMggb8HvP4sMWdy26qs3mps7dV9ev7La6YbMK\nqHY4bq2lcMPykk5t1xS1DRYGG98u3/jnYwU0zUcC1pY105GQEGebS++H1FOPQkoLjvO7tZ/i3T8U\nqm3W1UgNbQHSKul8n8OzMF2c7V7n3og/igNxE1+N9O/9Vzh5rc2l745cm/oA95y34Lk7Spvqx7DN\nXGoE28v+48nepnJeFCKGkJ1CHp1b/L31UW8p6N1ksfYobqoviPcNg1DfaGVTTYk1uLGVtfXmUmvx\n5mE7c8YZN+XdfCgbMIPZ1GE0vmAwznavs1C5bhqGMwYbYwWTc5qhegvej2Qkjn+q8m59atV5sxD5\n34eIDDkf+ySHItfBrAWRxZ2CyKq0OGoq4+Js9zoPQ+rxWiEO6Q+mSDdZ5HqMhyPGUH0QIzJnOKdj\n1e/xDbY5ydr9bu2FOI6/K6JWOBDpi3gB8ESWFw1yUft9KCKPHYBk2E7zOTzvR9r+6POr/taxfuMz\nO9UsXXrwqi9m9Tn80GnL3nxrNDA3WibJ7nUOQ0yedkeClRMybPIUEWs2bAlOvebNBxBJdAHSgeBC\nv1v7PF3vYfc690R8Cfog8t3TIrleV2uVA5GgdypSb7oAye4+VaFX1RHjvKjWKlsh319o7Wrz//tE\n2PcvtmVVQ7Osv6qa1f2Q76tZwn+S360tCXOcSIj7XLa59IlIS5tOyFg5/W4tGdntdlDmXB6kdnsT\nsuB4TySH9FiclVppEhIkH6Ce/hHJBD8Tq8e0zaX3RZQIuyMZ8ml+t1YLO9beEt05+RPg431fefGX\nj44+rhgpg2gZgEYMTNVrof2qo6IRS9PmglJrTUEJFljWvn7jV8XBhsVIAN/8WNXi79XAphCTNwsw\n5KFJt3zboax93s6HMgkzmE0dRuMLJudswGh8weScEagWE16gtHfbHta/NiwFqQW1x7pB5wlyPsaq\nh6Afqdl6GOnBG+3mFZOz3et0IJO7IGLWlYxrcrqQD2PcCxnjvZEJ9nHhFlpsLr0N4D5r8vDTx4/p\nH2uMC5AMy/VIS4lXkEWDRCbu6UKupNytEZXGeUiG8EHgkigu0rcDFwKv6bdPOspqtUQLAKxsM3kq\nRkoXLo7H5EnVsRap/YqRzFVxio+uZSUF59XUNoJkQi8GXotmeJQsVF3wK0hm+yPgmBb+BNuhWqvs\nh5hNnaq4/gncOvaF5/718eQpbRDX3nABa+9wxwvCwrVtCxoXdSvqv75NQV2bzU33DP+55mlEwhoz\nWFRKh8eQAG45MNXv1t6M8+PHPJdtLr0ACTCvQDL1ZwMPp/O7UOffVEQR0AVZSHH6HJ4P4uWsfCZO\nQc7jndXTbyJB7FuJKAtsLr0T0n+4ArmGHet3a2EN/JQqIrQt0HbZW2txMU11Mcu5Q1FD9OBz6/Pr\nC1pteKnnQZUrizuc02ApaI3F8gVwjt+tpeR2TR7cR3IFM5hNHUbjCybnbMBofMHknDHYvc5RSOar\nK+IkOy3KCna+IS/GWMnbZiO9Ca/wOTw3Rtk8VhbgbKR2bwOS6X0/jVSTQb6McWtEwqkhWZEjo/TF\njDXGzfW1+yDmNWcBL+SwpVqus9+jkKz3ECRz5/Q5PP6W29lceiGS/R530N59eP/LRVeyfcBYAhRb\nSja1Kx74zf7W1ut7BhsK6+r/GPx946re64k/8IzWviVpdGhTwtqNtU7gEb9by4Qr9VYoafWjiLnU\nH4jS5bto+1Rrlb0QV+EZQFmUoCUILGT7zOrPwC//HdV2xfcDyx5BJNS/IdeQqO8bDjaX3uw4fTuy\nsHArcKXfrcW6N0Q9l20uvRtynzlU8Zvsd2tfJMovXqh64ZuRbLwFaTd1kc/hWRGy2Xac7V7nzshn\nPxVoC2xBym/u8Tk8PyTLxebSy5DM7CSk3/JEv1uLWcbTMnvbepcB9k2/LniX8EFpy4B1TYVeFZdL\nuM2lH4kE6gOREqTLgMcTcYGPgry4j+QCZjCbOozGF0zO2YDR+ILJOaOwe507nTd2+h93f/xIQQ7q\n2FJB3oyxMgT5GOk7+Xefw/NUhE0jZQFCHaSXIvXK6e5nmgzyaYwLkLrWC5AgVPM5PB+H2TTSGBcj\nE7TLkQn6v4HzfA7PyoyRjg85H2M1Npci2bIiZNJ9bsuMosowzUXkrTvA2mE5xTt/h6Wonsa1Xahb\nMBQatiaVGpBMXMtHbYTnU3m0PGat76YjZ5eVFGbTmduCnGvXI8ZGU+Jxbq7WKrsBF5T17XNJzcJF\nb7O94dLPSIuXHdx67V7nICQjPBjJ/jl8Ds+qVD6DzaWPQM6FgYjMdYrfrf0eZZeI57KSMHsRczE/\n8He/W1uTCr94Yfc6xyBlNXshZmWXAQ+r+11QZXL3Q9QaGqJUWIIoCh5K1zVCLQh5kJr1BcB4v1sL\na5gWBWm9Xthc+iAkiJ2IyJrvA65J1AgtBnJ+jcsVzGA2dRiNL5icswGj8QWTczZgNL6QZ5yVKVw1\nspp/hM/heTvMZjtwVm0zHkAyAb8grZAiZR2zjbwaYwC71+lE2nE0IE7yLZ1Pw43xSCQbOxTJPs7w\nOTyvZ4FuPMibMbZ7nXsgdeBjkAzPucCzoVlrm0tvf/X0MWuvfeSTw1DBoqV0o6Vkt7nnWIprpwSD\n1AbrSm+o/W4/D02FzcFkfYs+tblAruTcxyGZvVIk83pnnCqARIyJJiK1ru2BuxBJd1qyzzaX3ha5\nPp2A1JWf6ndrkUofduCssrwXINldEPOrO7J9PqjrrBNxkG6HtKg756zR0z6579MnvkDkvCCmZnch\nao2E9LzxQI3HtcCViIx7ot+tfZnAIdJyHttcejtk8eo8ZAHrXURSnIk2Er8IxAAAE1VJREFUfHlz\njcs2Ugpmy8vLuyEn5KGIY98T6t/vgZmBQCBYXl5+GtIfrQG4IRAIxLqxGe3LMBpfMDlnA0bjCybn\nbMBofCEPOdu9zv2RVjG1wP5hsqstJW1lyCRUQ+5ZR0Srr8sB8m6MAexe5wTEBbUtkv26OSRA2MpZ\nje+1gAvJtjyETPTD1oXmCHk1xioDfhbisNsKcdQ+w+fw/BmyWegYD0eko80mT1OSkbZmAbl05h4J\n6EhW8jFEyh0rUIqnvt6C1NvehCwanB5FFZI0VAA2DcnatUKynC6/W2uZId6OszJmexw4BlGcHO93\na7PTzS8RKBXNHUhw3owmJKt9N/BhNkoObC79TGQ8NwHH+N3aO3HumtJ5rGrS/4aoXHogMvgLgJcz\nUUOukFfXuGwi6WC2vLy8CLnJDUYmCLcDdwQCgTnl5eUepID7E+AtZCWmDHHvHBkIBKJdXIz2ZRiN\nL5icswGj8QWTczZgNL6Qp5ztXudk5B60GBjTojVEaBDQETEF2Q8JgI/1OTzxtsDIFvJyjGGra+lr\nwE7IhPkMFSAEAYtaWHgUMcxZAEz3OTzv5YpvFOTlGKvawYeAwxCZ7CXAAyHSzAIkc3sLUut6L2Ig\nFVeNXg6Q69rk3khAOwKYA1TGkK/Gqv1ujZzfDkRtcIzP4YnZiioVqLZGPqS++hvA4XdrgZBNtnK2\nufQ9gRcRA6v3EYlyNvsKR4Xd6zwEcB0x6JAj3pj/7oBcqGFsLv04pNzBghhtPR/HbkmfxzaXPhLp\n3z0WqQW+Gbjd79Yy/ZvNy2tcNpBKMHs3YlJwGdLQ+L+BQKCPem0ScDgS0B4RCASc6vmXgJsCgUA0\ni3ajfRlG4wsm52zAaHzB5JwNGI0v5DFnu9d5AdLDcx6wX4j7bnOg1Qu5Dw1BatL+nglJWxqQt2MM\nW51j/cBI4D2g8slj71r995fOvx84E+F/N3Clz+FJue1HhpC3Y6wyf9OAO5H+lB8C0x+cdMtPM169\n9C1kPrUcONnn8LyRM6LxIefjbPc6WyEmRMchCyy2KKZCEfnavc7+SCZxOFLacJzP4clKoKiMjO5C\njKo2ATP9bu1J9XIQsNhc+ilI/+FSJGC6KtOmWykgp+eFzaUfhCxytAPO97u1u2PskjBfZbx1I1LK\nYkEWGS5MtbdyAsj5by9XSCqYLS8vnwb0DgQCN5aXl7+H6OP/GwgEeqvXD0astmcBQwOBwKXq+SeB\npwKBwH+jHN5oX4bR+ILJORswGl8wOWcDRuMLecxZBQF3IZmr9xFDp1oko7UbEsj2Q7JZ5+Wx8Vbe\njnEzVIDwDCJlnN+lVadBKzevBnE9PsXn8HySS35xwAhj3AORRFYCtW2KW5dsrNsEIkE+2efwxHRl\nzQPkxTgrs6GrgauA9YhR06wwm0YyMjsIeAFpOfMQcHYuFsJsLn0yUl/dDmklduYLNx+5YfJlrz8O\nnIzU157kd2sxTa9yjJyfFzaXPhyJS3og8t9Lo0h+4+Zrc+lFyKLetUg99TykLvbdlEknhpyPca6Q\nbDA7Gxm0IGJlPR/YKxAIFKvXNUQy8xYwIRAIzFTPv4TUzUYrws47RyoTJkyYMJGfaGpq4q6PH+HT\nRV9RsdNIzh5zMgtW/8nNc+5jQ90mpgzVOHrweCyR23WaiBNNwSae/fYVXv3pbQosVrTB46ncfSJF\nBUW5pvb/Cp8s/JJHv/SyuW4zJw4/lgm7HmSev0niwz/m4vnsKRqCjUzbc3LMsQwGg8z6+X2e/PpF\nLMApex/PuIH7Z49wGCxdtYnbn/mc+X+upVeX1hQXFfD7kvUM7NOeS6buQ4/OrXPKz0hYtnozVz/0\nEX+t2MQhI/tytn1PCgusSR/v6/nLeeiV71m4bAOty4o4ccJuTBzbn4IUjvk/hrRc2FJ2M1aZ2TOQ\nmll3IBCYXV5e/gBiWT4HeBvpM1eK1NAON2tmcw6Tc+ZhNL5gcs4GjMYXDMBZGRC9A+wLvFhSWHJc\nbUNtE2LU8mhu2cWFvB/jUNi9znF3TrjqrT7texqGM8Yb4zYPabdu6FDazjCcFfJunO1e52hEYtod\ncQw+J6TPd2h9fQnS0uVkRNZd6XN4Psw+4x1hc+nFiIT1QvWUB7ggjDlUviJvzgubS++K+ACMQlQP\nk/1urWV5RFS+NpfeHylxOVZt+xDSI3hFpH2ygLwZ42wjXcHsDGQQH0YMCn4ATlNuxtMRN2MrcGMg\nEHg5xiGN9mUYjS+YnLMBo/EFk3M2YDS+YBDOdq+zM/ARMKiooIj6xnrN5/C8mmteccIQY9wCRuNs\nNL5gck4b7F5nX6TuezjSHmWyz+FZzfb19S8BoxHH82NamMrlBWwu/YCrp4+ZPXJw97wb4xjIq/PC\n5tLbIDLyCUj7oCP9bi3UKCwsX5tLb4WYtF2MJOk+As5OsO1PppBXY5xNmH1mU4fR+ILJORswGl8w\nOWcDRuMLBuKsnGGvve6QC0/aresuhuCsYJgxDoHROBuNL5ic0wq719kGqTs9GvgZOMrn8ATsXudY\nJJDtidSFn57HbtGQx2McBXnHWdW6PgqcBASACX639rt6eTu+qm3ScUg2ti+wBOln/GwGW+0kirwb\n42zBDGZTh9H4gsk5GzAaXzA5ZwNG4wsm52zAaHzBeJyNxhdMzmmHMoa6AenEsfboweM7vPLjm3VA\nIRKc3JWN/qcpIq/HOALykrPqB3sL8t0vQQLab9m+/dFQpNXOQUif4TuBm/xuzWzxlicwg9nUYTS+\nYHLOBozGF0zO2YDR+ILJORswGl8wHmej8QWTc8Zg9zpPQlyCi4E1iNvx27llFTcMMcYtkNecbS79\nfCRIXQdofrf2vs2ld0Ycis9ESiVfQ+qUf84d06jI6zHOJMxgNnUYjS+YnLMBo/EFk3M2YDS+YHLO\nBozGF4zH2Wh8weScUdi9zrGHDKj46N0F1QN9Ds+vueaTAAwzxiHIe842l34C8ATQdMxBA0tefv+X\nVUBnRJJ+nt+tmT2e8xRmMJs6jMYXTM7ZgNH4gsk5GzAaXzA5ZwNG4wvG42w0vmByzgaMxhdMzhmD\nzaUfjtRPtwY2AtcB9/jdWtZ7DCcBQ4xxJmAGs6nDaHzB5JwNGI0vmJyzAaPxBZNzNmA0vmA8zkbj\nCybnbMBofMHknFHYXPqex48r/+r5twO9/G5tSa75JADDjHG6YQazqcNofMHknA0YjS+YnLMBo/EF\nk3M2YDS+YDzORuMLJudswGh8weScDRiNLxiTc1pgzTUBEyZMmDBhwoQJEyZMmDBhIlGYwawJEyZM\nmDBhwoQJEyZMmDAczGDWhAkTJkyYMGHChAkTJkwYDmYwa8KECRMmTJgwYcKECRMmDAczmDVhwoQJ\nEyZMmDBhwoQJE4aDGcyaMGHChAkTJkyYMGHChAnDwQxmTZgwYcKECRMmTJgwYcKE4VCYzE7l5eUF\nwMPAIKSv0RlALfAE0AR8D8wMBALB8vLy04DTgQbghkAg8HoaeJswYcKECRMmTJgwYcKEif9hJJuZ\nPQpoCgQC+wFXADcBbuAfgUDgAKRpr1ZeXt4DOBvYFxgP3FxeXl6cOm0TJkyYMGHChAkTJkyYMPG/\njKSC2UAgoAMz1J/9gTXAiEAgMEc99x/gMGAfoDoQCNQHAoH1wC/AsJQYmzBhwoQJEyZMmDBhwoSJ\n/3kkXTMbCAQay8vLnwDuAf6NZGObsQFoD7QD1oV53oQJEyZMmDBhwoQJEyZMmEgaKRlABQKBaUA5\n8AhQGvJSO2AtsB5oG/J8WySLGw2WGK/nG4zGF0zO2YDR+ILJORswGl8wOWcDRuMLxuNsNL5gcs4G\njMYXTM7ZgNH4gjE5pwVJBbPl5eUnlZeXX6b+rOH/2jv3WLmqKg5/vS2tPEp9AaIxxgguCUYCrQ8Q\ne1E0WsDIQ0nw1UJpa9BqSIhi4zNqNFEr9dXwlJJUDRpQiTwEJKCikMZGSQgLkL9M+KNBbWkVqLfX\nP/aeMPfm7HM7pzPnrIm/L2ly5uyz7nxZXbPm7Jlz9sAUsM3MJvO+FcB9wIPA28xskZktAY4jLQ4l\nhBBCCCGEEEI0Zt709PTAQWZ2MGnl4pcBBwFfBx4hrXC8EHgYWJNXM76YtJrxBPA1d795OOpCCCGE\nEEIIIf5faTSZFUIIIYQQQgghuuSA7pkVQgghhBBCCCG6QJNZIYQQQgghhBBjhyazQgghhBBCCCHG\njgWjfgIzmwB+CLwBeBa42N3/1jf+XuDzwH+B69z9mlKMmR1DWnhqH2lV5I+7+1Bv+m3oexBwHfAq\nYBHwVXe/xcxOBG4BHsvhm939xmH6NnXO+//M878D/IS7r24jx02dzWwVsDIfcjBwAnAU8BpGnOe5\nfPMxhwB3Ahe5u3dZxwfgHLqWq5zzvk5qeUDfC939UTNbCazKw63W8f44m9kFwKdIr72HgEtIS/5H\n7clVvgsIXMdVznnBxMg9uSrPKwnak83sPOAzwDSw1d2/G70nF5xD9+Qq57w/ZE8u5HgVHdXx/jj3\nHXcV8JS7fzbyeXLBN3QdVznnx2F7cpVz5PNkM7sUWA3syLvWAo8Dm2fHDJrjNr6ZPRtY6O6nAJcD\n3+4N5OLeCLwLmATWmtmROWZRRcxGYIO7LyedXL0viO+HgB3Z6z3A93PIUmCju789/xv6C7Sh8xFm\n9gKAPrfVOaSNHDdxPtLdr+/5AtuA9e6+i3byXPTNzstIP0f1atKbZF1M5zmucQ5byyXnjmt5EF+y\n55YO67jWOa9U/xXgNHc/FVgCnEXcnlzyDVvHJefgPbnSOWpPNrP5pF9ZOB04GbjEzF5Ct3Xc1PnD\nxK3lKucXR+3JJd+O67jWuc99HfB6YpxfNPENW8cl58g9ueQctSdnTgI+0ufwGHBOIWagHLcxmX0r\ncDuAuz8ALOsbOw543N13uvte4PfA8hxzW0XMSe5+X96+DXhnEN+fAV/Ix0wAe/P2UuBMM7vXzK4x\ns8NG4NvEeZL0ac0hZnaHmd1tZm/Ox7eR4ybOy3uDeYJwvOdvmEkvkFHnuc4X0k9SnQ34fsREyHHJ\nOXItl5y7rOUmvkBndTyX8zPAye7+TH68IO+L2pOrfP9D7DouOUfuySVnIF5Pdvcp4HXu/jRwBDAf\neI5u67iJ87PAjQSt5Zo8h+zJNb5A2J6MmZ0CvAm4knRiXxfTdU8u+Yat4xrnyD255NwbC9WTM0uB\nDWb2OzO7fI6YgXLcxmT2cGBX3+Op/FV0b2xn39jTpE98q2LmM/M/a3c+dtgM7Ovue9x9t5ktBn4O\nfC6PPwBc5u6TwBPAF0fg28gZ2AN8093fDXwM2Npijps699gAfKnv8YOMPs91vrj7/e7+9/2IiZLj\nSufgtVzKc5e13MS3Rxd1DDXO7j7t7jsAzGw9cKi731mI6TzHBd+7ItdxyZnAPbnGuUfEnrzPzM4F\ntgP3kPIbvSfPdv535FouORO7J1f59gjXk83saNIk8BPMzF/InlzyjVzHNTkO25NrnHuE68nAT4B1\nwDuAU83szELMwDluYzK7C1jc/5zuvi9v75w1thj4VyFminTt9Oxjh82gvv8EMLNXAr8Ftrj7T/P4\nze6+PW//AjhxBL5NnR8FtgLkr/qfAo6mnRw3dcbMXgi81t3v7RtvI891voPEtFXHpeefyzlyLZfo\nspab5rirOoY5nM1swsy+RboU77yamAg9ueQbuo4LzpF7cl2ew/Zkd78JeAXpHr2PFmJC9eQK59C1\nXHAO3ZMLOY7ak98PvBS4lXSv7wctrbsQtSdX+Uav45Jz5J5cl+eoPXmTu//D0xWXv84OQ6njNiaz\nfwDOADCztwB/7Rt7BDjWzF5kZgtJl5LeXxOz3cwm8/YK0r1oXfv+0cyOAn4DfNrdr+87/nYze2Pe\nPp10/fooGNgZuJB8bbqZvZxULE/STo6bOpO37571t9rIc53voDERclxJ8Fou0WUtN/GF7uoY5na+\nknSSd44/f1lp1J5c6TsGdVyV48g9ueQMAXuymR2eL6db6GnRkD3AVE1M5zkuOUeu5Zo8h+zJNb4Q\ntCe7+/fcfZmneyC/QVq0aktNTKc5rvD9sbvfELmOCzm+AbiIoD25lOc8HLEnLwEeMrNDzWwe6dvZ\nbTUxA+V43vT00BfgmkGW7q1uBanJLQUOc/erzews0lflE8C17r65KsbTKqDHAleT7kN7GFjjw1+l\nrYnvJuADzLwvbgVgwA9I9wY8Cax1993D9D0A5wXAj0gry0FqMH9qI8dNnXPcZcBznldMzPtOYMR5\nnsu377h7gHW5Xjur4wNwDl3LBefOarmJb37cSR3P5Ux6c9nGzDeOK4BfzY6J0JMLvpuA04DzCVjH\nBecrSJ9Uh+zJJWd3/2XUnmxma0grZ+4F/gKsz8eF7ckVzp8EvkPgnlzI83yC9uQqX08riYfsybPe\nR1YC5u4bIp8nF3zH5dyi3zn0eXKVc34ctSdfAFxKWgvgLnf/8rDqeOSTWSGEEEIIIYQQYti0cZmx\nEEIIIYQQQggxVDSZFUIIIYQQQggxdmgyK4QQQgghhBBi7NBkVgghhBBCCCHE2KHJrBBCCCGEEEKI\nsUOTWSGEEEIIIYQQY4cms0IIIYQQQgghxg5NZoUQQgghhBBCjB3/AzKlEdZrZKssAAAAAElFTkSu\nQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, axs = plt.subplots(3, figsize=(16, 9), squeeze=False)\n", "axs[0, 0].bar(range(num_pops), fix_pops)\n", "axs[0, 0].set_xlim(0, 11)\n", "axs[0, 0].set_xticks(0.5 + np.arange(num_pops))\n", "axs[0, 0].set_xticklabels(my_pops)\n", "axs[0, 0].set_title('Monomorphic positions')\n", "\n", "axs[1, 0].bar(range(num_pops), [np.max(vals.keys()) for vals in num_reads])\n", "axs[1, 0].set_xlim(0, 11)\n", "axs[1, 0].set_xticks(0.5 + np.arange(num_pops))\n", "axs[1, 0].set_xticklabels(my_pops)\n", "axs[1, 0].set_title('Maximum number of allele reads per loci')\n", "\n", "for pop in [0, 7, 8]:\n", " axs[2, 0].plot(MAFs[pop], label=my_pops[pop])\n", "axs[2, 0].legend()\n", "axs[2, 0].set_xticks(range(num_buckets + 1))\n", "axs[2, 0].set_xticklabels(['%.3f' % (x / (num_buckets * 2.)) for x in range(num_buckets + 1)])\n", "axs[2, 0].set_title('MAF bundled in bins of 0.025')\n", "None" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [ "#for pop_data in pop_iter:\n", "# pop_name, allele = pop_data\n", "# print(pop_name)\n", "# for locus_name, vals in allele.items():\n", "# geno_list, heterozygosity, allele_cnts, summary = vals\n", "# print(locus_name, vals)\n", "# break\n", "# break" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2446/NA20127\n", "1375/NA12264\n", "NA18594/NA18594" ] } ], "source": [ "exp_hes = []\n", "for pop_data in pop_iter:\n", " pop_name, allele = pop_data\n", " print(pop_name)\n", " exp_vals = []\n", " for locus_name, vals in allele.items():\n", " geno_list, heterozygosity, allele_cnts, summary = vals\n", " cexp_ho, cobs_ho, cexp_he, cobs_he = heterozygosity\n", " exp_vals.append(cexp_he / (cexp_he + cexp_ho))\n", " exp_hes.append(exp_vals)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [ "fig = plt.figure(figsize=(16, 9))\n", "ax = fig.add_subplot(111)\n", "sns.boxplot(exp_hes, ax=ax)\n", "ax.set_title('Distribution of expected Heterozygosity')\n", "ax.set_xticks(1 + np.arange(num_pops))\n", "ax.set_xticklabels(my_pops)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.9" } }, "nbformat": 4, "nbformat_minor": 0 }