{ "metadata": { "name": "abundance-hists" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "def read_hist(filename):\n", " return numpy.loadtxt(filename)\n", "\n", "transcripts_all = read_hist('transcript-reads.fa.hist')\n", "transcripts_dn = read_hist('transcript-reads.fa.keep.hist')\n", "transcripts_trin = read_hist('transcript-reads.fa.normalized_K25_C20_pctSD100.fa.hist')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 17 }, { "cell_type": "code", "collapsed": false, "input": [ "print 'Unnormalized data:', transcripts_all[:,2][-1], 'total k-mers'\n", "print 'diginorm normalization:', transcripts_dn[:,2][-1], 'total k-mers remaining'\n", "print 'Trinity normalization:', transcripts_trin[:,2][-1], 'total k-mers remaining'" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Unnormalized data: 3163778.0 total k-mers\n", "diginorm normalization: 1976319.0 total k-mers remaining\n", "Trinity normalization: 190815.0 total k-mers remaining\n" ] } ], "prompt_number": 29 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(transcripts_all[:,0], transcripts_all[:,1], label='all data')\n", "plot(transcripts_dn[:,0], transcripts_dn[:,1], label='diginorm')\n", "plot(transcripts_trin[:,0], transcripts_trin[:,1], label='Trinity norm')\n", "title('k-mer abundance plot of normalized data')\n", "xlabel('k-mer abundance')\n", "ylabel('N of k-mers with that abundance')\n", "legend()\n", "axis(xmax=40)\n", "savefig('normhist.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEVCAYAAADQC4MUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XtYVNX6B/DvgCIIJGAyoNzk4sCAKCAiKg6aJzFDDTUw\n8QIqRpaZ/lLLjqCWlnbyUonp8VJKoWiklUppjehJBwULGZBSSVK5iHgBBLm9vz+IHaPAMIjI5f08\nD4/M2nut/c7Fedlrr72WiIgIjDHGWAvRetIBMMYY61g48TDGGGtRnHgYY4y1KE48jDHGWhQnHsYY\nYy2KEw9jjLEWxYmnjbKxscGxY8eedBhNoqWlhcuXL7fIsdra69SSr01ERAR69+6NQYMGtcjxmktk\nZCSmTp0KAMjKyoKhoSGa+64QTT43O3fuhI+PT7Mev73jxNNGiUQiiESiJx1Gq9deX6dH/bK7evUq\nPv30UyQlJeH06dPNGNnjV/v9tLKyQmFhYbO/x4/rc1M7aXZknHjYI6msrHzSIbAmOHnyJCQSCUxM\nTB77sZr7M8L3vLd9nHjagfT0dNja2mLPnj11bo+MjMTkyZPxyiuvwNTUFMOGDUN2djY+/fRT2Nvb\nY8SIEUhMTBT2LykpwbZt2zBw4EAMHToUsbGxwn/2nTt3YujQoYiIiIC1tTWWL1/+0PESExPh7e0N\nY2NjeHt745NPPkFFRYXKPgkJCejXrx+cnZ0RExMjtP/gX4R//vkntLS0UFVVBQDw9fXF+++/j1Gj\nRsHc3BwLFizA7du3hf1PnTqFZ599FjY2Nvj44481iktLSwvR0dHo378/7O3tsW7dOpX6P//8M4KD\ng2FiYgJnZ2ecO3cOAHDr1i2sX78ezs7OGD16NH744Yd63ilgxowZmD9/PgICAmBmZobFixfj5s2b\nde577949REVFoW/fvnj22Wfx7bffAqh+v8PDw3Hq1CkYGhrWmzwKCgrwwQcfwMHBARMnTsTx48cB\nANu2bUNISAgUCgUMDQ3rfA9r3ucVK1agV69e8PPzw6lTp9S2DfzzeQsPD4e5uTl27twJX19frFq1\nCv/617/Qo0cPTJ8+HeXl5Xj55ZfRs2dPhISEIDs7W2jj9ddfh5WVFcRiMebMmYPffvutzudY+/NR\n83rU/Ojq6qJ3797CvocPH8bYsWMhkUiwbt06FBUVCdsa+tw8qKioCKtWrYKFhQWeeeYZ5OTkqGyv\nL/YjR45g9erV2LNnDwwNDeHm5gYA2LFjB6RSKYyMjDB+/HjhfW7XiLVJNjY2dOzYMUpKSiIrKyv6\n/vvv6903IiKCdHR0aMeOHZSfn09jx44lJycnmj9/PuXn59PKlStpxIgRwv5vvPEGBQUFUWZmJv36\n66/k4uJCP/zwAxER7dixgzp37kxvvfUW3b59m0pKSh46XlJSEikUCqqoqKD//e9/ZG1tTT/++KOw\nXSQSkbe3N6WkpNDx48fJxsaGjhw5QkREkZGRFBwcLOybmZlJIpGIKisriYhIJpORpaUlHT16lK5e\nvUqenp703//+l4iIbt26RV27dqXPP/+crl27RlOmTKHOnTvTsWPHGh3XiBEjKD09nZKSksjQ0JAu\nXrxIRETJyclkampKX375JZWVldHFixfpypUrRET0wgsv0Lx58ygnJ4cSEhKoZ8+e9Mcff9T5Xkyf\nPp309fWFGIOCgigoKEglhkuXLhER0bJly2j48OGUkZFBx44dIxsbG/r555+JiGjnzp00dOjQet9z\nIqJp06bRiy++SH/99Rft37+fTExMKDMzs1H1d+zYQTo6OrR8+XIqKCigiIgIlf0bajsiIoI6d+5M\nGzdupJKSEiopKSGZTEbW1tYkl8vp0qVL1Lt3b3J2dqbPP/9c+EyuWLFCaH/37t1UUFBAt27dojff\nfJOGDBkibIuIiBA+Iw9+PmqUl5eTTCajt99+m4iIDhw4QK6urnTq1Cm6fv06vfjii8I2dZ+bBy1c\nuJBGjRpFly9fpgMHDpBYLCYfH59GxR4ZGUlTp05Vae/777+ny5cvU1lZGX311Vekp6dHRUVF9b43\n7QEnnjbKxsaGli1bRhYWFnT8+PEG942IiCBXV1fh8e7du6lTp05069YtIiK6evUq6erqUmlpKVVV\nVZGNjQ1lZWUJ+69bt45eeeUVIqr+QtLT06P79+83OtalS5fSq6++KjwWiUS0bds24fFbb70lbK/9\npUL08BeLr68vvfbaa8L21atXU2BgIBER7dmzR+UL4NKlSyQSier9Aqkrrv379wuPR40aRVFRUURE\ntGjRInr99dcfauPu3btkbm5O9+7dE8rmz59Pa9asqfOY06dPV4kxIyODjI2NhedXO/H069eP4uPj\nVeKdN28eEVW/Dw0ljoqKCurevTtlZGQIZVOmTKGPPvqoUfV37NihEtf169epc+fOVFRUpLbtiIgI\nsrW1VWnP19dXiJ2IaPbs2eTu7i48jo6OJplMVmcsd+/eJX19fbpx44bQvrrE8/LLL5O/v7/w+KWX\nXqLo6Gjh8blz50gqlRKR5p8bqVSqsm3q1Kn1vpYNxV6foUOH0r59+xrcp63jrrY2iojw2WefYciQ\nIRg2bJhQHh0dLXQ1jBkzRih3dXUVfheLxbCwsICRkZHw+P79+7h58yYuXLiArKwsuLq6wtjYGMbG\nxoiIiMD//vc/oX6/fv2go6NTb2zXrl3Dyy+/DFdXVzz11FNYt24dUlJSVPbp37+/8Lubm5tKN446\nteuamZnh2rVrAACFQoF+/foJ22xtbdGtW7cmx2Vubo7r168DAORyOYYMGfJQLCdPnsSNGzfQs2dP\n4fXatm0bTp48WWfsIpFIJcY+ffqgvLwc6enpKvsVFhYiJSUFHh4eQpmHhwdOnDhR/wtTS3p6Ou7f\nv48+ffo0qT4AODs7Q0ur+ivC3NwcFRUVyM3NbVTbXl5eD7VX+3mLxWKVz6SpqanwPgLVXX1jxoxB\njx49YGVlhZKSEpw/f75RcX/22WdISEjAl19+KZQdPXoU4eHhwns0fPhw/Pnnn8jNzVX7uant7t27\nSE9Pf+jzW5umsZ88eRKTJ0+GtbU1jIyMkJiY+NDnsr3hxNNGiUQifPbZZ7hy5QoWLFgglE+ZMgWF\nhYUoLCzE999/L+zbWBKJBBYWFkhLS8OtW7dw69Yt3LlzB7/++quwT6dOnRps491330V5eTkOHTqE\nO3fu4I033hCu0dSouT4CAMnJyRg8eDAAwMLCArm5uXXup46Xl5dKnJcuXcKdO3c0iqs+w4cPrzOZ\neHt7o0ePHsjNzRVer7t37+LAgQN1tkNEKjFmZGSgc+fOcHJyUtnP0NAQrq6uOHv2rFB29uxZ4Y8M\nbW3tBi+yOzo6okuXLsjIyKiz/qNoTNva2tpNbv+vv/7CggUL8Pbbb+PKlSvIysqCnp5eowYVnDhx\nAsuWLcOBAwdgYGAglI8YMQJbt24V3qNbt26huLgYYrFY7eemtqeeegqOjo4PfX5r/o+pi71Tp04q\nz4OIMGfOHMhkMiQnJ+P27dsYOHBgux9AwYmnDTM0NMSRI0eQkJCAt956q979NPkQa2lpITAwEIsX\nL0Z6ejqqqqpw6dIlJCQkNLqN69evw8TEBN27d4dcLscXX3zx0D7bt29HamoqTpw4gT179uD5558H\nUP0Fcfr0aSQnJyMjIwOffvppo5/Ps88+i+TkZERHR+P69etYvny5SpJsTFwPHqfmWEFBQdi7dy/2\n7t2LsrIyXLx4EVlZWTAyMsLQoUOFL5rKykqkpqaqJIwHnTt3TohxxYoV8PPzE84sahs3bhzWrl2L\n33//HXK5HF999RXGjx8PoPoM448//lC5QF5bp06dMGbMGERERODatWv45ptvcOTIEaH+o2hq2w9+\n4dbnxo0bICKYmZmhsLAQb7/9Nu7fv682rr/++gsvvvgidu3aBXt7e5VtU6dOxZo1a3Dy5ElUVlbi\nxo0bOHjwIAD1n5sHPffcc/jwww+RmZmJ7777TuV+H3Wxe3h4IC0tTSgrKyvDjRs3IBaLoaurix07\ndkChUKh9rm0dJ542rlu3bvjxxx9x+PBhRERE1LlPXfckNPQ4MjISw4cPR3h4OExMTDBp0iRh5E5j\n7m+IjIzEr7/+CgsLC6xduxavvvqqSh2RSISwsDBMmTIFc+bMwbvvvot//etfAKq7OSIjI/Hiiy9i\n8uTJmDVrVoOx1o7HyMgI8fHx2LFjB7y9vTFw4EBYWFhoFFd9r1v//v0RHR2Nr7/+GqampggICMCt\nW7cAAJs3b4a1tTUmTpyIHj16ICwsDHfv3q3ztRGJRJg9ezb2798Pd3d39OrVCxs3bqwzhkWLFmH8\n+PEICAjAe++9h48++ggymQwAIJVKMX78eDg7O8PU1LTOY3300Ufo168fZDIZvvjiC8TGxsLGxuah\n51ZfnA297k1pu7737cHt7u7ueOWVVzBixAgMGzYMLi4usLS0rLduze/Hjh1DXl4eJkyYIHQ39+3b\nFwAwevRorFixAp988gl69OgBb29vYSSnus/NgyIiIjBkyBAMHToU69evx7x584Rt6mKXyWTo06cP\nevfujQEDBqBLly7YsGEDVqxYAXt7e5w9exZBQUH1Hru9EFF7P6djrBUJCQmBhYUFVq5c+aRDYeyJ\nafYzntDQUIjFYuEvDQBIS0vD888/j/79+8Pf31/lQurGjRvh4OAAqVSq0oeenp4Od3d32NraYunS\npUJ5eXk5Zs6cCWtra/j6+qqMoY+NjYVEIoFEIsG+ffuE8uzsbMhkMlhbW2PWrFl80yN7YvjvPMYe\nQ+IJCQnBkSNHVMpWrFiBadOm4ddff8VLL72EFStWAADy8vKwadMmHDt2DFFRUSqnrAsXLsTixYtx\n5swZHD9+XOgzj4uLw507d5Ceng4/Pz+8++67AICqqiosWbIE+/fvR2xsLJYsWSK0tXLlSowZMwZp\naWkoKChAXFxccz9txhqlvU7hw5gmmj3x+Pj4wNjYWKWsW7duuHnzJqqqqnDz5k1hu0KhgJ+fH6ys\nrCCTyUBEwsXSjIwMBAYGonv37ggICBAuuCkUCgQHB6Nr164ICwsTypVKJVxcXODi4gJXV1dIpVIo\nlUoA1Xesh4WFQV9fH8HBwR3i4h1rnXbs2CH84cVYR9UigwvWrl2LDRs2wNjYGJ988gnWrFkDoDoh\n1B5GKpFIoFAocPHiRZULplKpVJjIMDExEVKpFABgYmKC3NxclJaWQqFQCOW165SUlCAvL0+4Z8XJ\nyanNTYrIGGPtScM3ZDST0NBQvPbaa5gzZw4+/fRThIaGYu/evXX2d9fVDUFEQnntIa41j+tT063R\n2H517gJhjLGm0ei2jccYh+DkyZMIDQ1Fp06dMHPmTOGeEC8vL6SlpQn7XbhwAZ6enrC3t1e5iTAt\nLU24E7p2nYKCAmH8+4Nt1dTR1dWFqampMPS1dlt1qUlsrfknIiLiicfAcXKMHCfHWfOjqRZJPMOH\nDxdu1jpw4IBwz8bAgQMRHx+PrKwsyOVyaGlpwdDQEED13dExMTHIz89HXFycSuLZvXs3iouLsWXL\nFmERK6lUitTUVJw/fx4pKSlQKpVwdnYW6mzZsgXFxcWIjo5ucwtfMcZYe9LsiWfy5MkYPHgwMjIy\nYGlpiR07duCdd97BN998g379+uHQoUPC8GixWIzw8HCMGDECr7zyCjZs2CC08+GHH2LNmjXw9PSE\nj48PBgwYAAB44YUX0K1bNzg5OeHIkSN45513AFRP0bF69WpMmDABEydOxPvvvy+09c477+DQoUNw\ncnKCsbExXnjhheZ+2owxxhqJbyCtRZPrQU+SXC6Hr6/vkw5DLY6z+bSFGAGOs7m1lTg1/e7kxFNL\nW0k8jDHWmmj63dkio9oYYx2XiYmJMLiHtW3GxsYoKCh45Hb4jKcWPuNhrPnx/6v2o773UtP3mGen\nZowx1qI48TDGGGtRnHgYY4y1KE48jDFWh507d8LHx0d4rKWlhcuXLzeqbmRkJKZOnfq4QmvzOPEw\nxlgz02TexxkzZuDf//73Y4ym9eHEwxhjrEVx4mGMdVjvv/8+7O3t0b17d0yZMgUnTpxoUjs3btzA\nokWLYGZmhgkTJgjritWYNGkSzM3NYWlpiQULFghddlu2bMGXX36JNWvWwNDQEOPGjWvWuForTjyM\nsQ7L3t4eJ0+exNWrVzFgwAC89NJLTWonPDwc2dnZ+PXXXzF27FhERUWpdLeNGTMGFy9eRGJiIvLz\n87Fs2TIAQFhYGKZMmYLFixejsLAQBw4caNa4WitOPIyxJ04kap4fTU2cOBFmZmbQ09PD/PnzIRKJ\nkJSUpFEbFRUVOHr0KJYvXw4zMzNMnz4d7u7uKvvMmDED+vr6MDc3x7Jly3Do0CFUVVUJ2x+8+bI5\n4mrNOPEwxp44oub50dTBgwcREBCAnj17wsTEBNnZ2UhJSdGojfT0dFRVVcHW1lYoc3d3V0kmH374\nIUaOHAkTExN4enri9u3buHLlirD9wcEIzRFXa8aJhzHWIRUXF2P27NmYPn06Lly4gIKCAvTq1Uvj\n6X0cHR2hpaWFS5cuCWVJSUlCMlEoFPjoo4+wbt06ZGdn48yZMwD+OcvR1tZWOftprrhaM048jZSZ\nCbz++pOOgjHWXAoLC1FUVARzc3NUVVVh9erVuH79usbtdO7cGSNHjsTy5cuRk5OD3bt349dffxW2\nX7t2Dfr6+jA1NUV2drZwfaeGh4cHUlJSUFFR0axxtWaceBqpvBw4fPhJR8EYay5mZmZYvXo1pk6d\nin79+qGsrAxDhw4VtotEIpUusIbuzdm0aRNMTU3Rv39/xMXFITw8XNg2fvx4jBgxAv3794e/vz8C\nAwNV2ho7diy0tLTQq1cvBAQEqI2rXaBmFhISQqampuTi4qJSvn37dnJ0dCSpVEqLFi0Syjds2ED2\n9vbk5OREJ06cEMrT0tLIzc2NevfuTW+//bZQXlZWRqGhoWRlZUUymYyys7OFbXv37qU+ffpQnz59\nKDY2Vii/fv06DRs2jKysrGjmzJlUUVFRZ+wNvRzXrhGZmTX+dWCMVXsMXzPsCanvvdT0PW72T0RC\nQgIlJyerJJ7z58/ToEGD6Pfffyciory8PCIiys3NJYlEQleuXCG5XE5ubm5CndGjR1NMTAzl5+fT\nkCFD6MyZM0REtGfPHpowYQIVFxfT6tWrae7cuUREVFlZSba2tnT+/Hn67bffyM7OTmgrPDycPvjg\nAyoqKqIXXnhBJSnV1tCLd+cOkYFBE18UxjowTjztR3MlnmbvavPx8YGxsbFK2eHDhzFz5kw4ODgA\nAHr06AGg+qKbn58frKysIJPJQETCjVcZGRkIDAxE9+7dERAQAIVCIdQJDg5G165dERYWJpQrlUq4\nuLjAxcUFrq6ukEqlUCqVAIDExESEhYVBX18fwcHBQh1N6OsDxcVArWuAjDHGmqBFViD94Ycf4Ozs\njAEDBqB///5YsGABpFIpEhMT4eTkJOwnkUigUChgbW0NU1NToVwqlSI6Ohpz585FYmIi5syZA6B6\nZcPc3FyUlpZCoVBAKpWq1Dl9+jRsbW2Rl5cHIyMjAICTkxPWrVtXb6yRkZHC776+vsJ659ragK4u\nUFJSnYQYY6yjksvlkMvlTa7fIomntLQUBQUFOHHiBI4ePYpXX30VP/30U70r2T2IiIRyqu4eVNlW\nn5qLgw3t86DaiedBhoZAUREnHsZYx1b7j3IAWL58uUb1W2RU26BBgxAYGAg9PT34+/vjwoULKC0t\nhZeXF9LS0oT9Lly4AE9PT9jb2yM3N1coT0tLg5eXFwCo1CkoKIBYLIauru5DbdXU0dXVhampqbDm\ne+22NGVgABQWNqkqY4yxv7VI4vH29sbhw4dBRFAoFLCzs4Ouri4GDhyI+Ph4ZGVlQS6XQ0tLC4aG\nhgCqb8qKiYlBfn4+4uLiVBLP7t27UVxcjC1btmDQoEEAqrvWUlNTcf78eaSkpECpVMLZ2Vmos2XL\nFhQXFyM6OlqooykDg+ozHsYYY4/gUUY41CUoKIjMzc1JR0eHLCwsaPv27VRRUUFz5swhR0dHGj9+\nPCUmJgr7r1+/nuzs7MjJyYkSEhKEcqVSSW5ubmRjY0NLliwRysvKyigkJIQsLS0fGk69Z88ecnBw\nIAcHB9q7d69Qfu3aNRo2bBhZWlpSaGhok4ZTExENHkxUa8Q3Y6wRHsPXDHtC6nsvNX2PRX9XYoDa\n60GjRgFvvAH4+bVgUIy1cZpeZ2WtV33vpabvsdqutps3b+KDDz7A2LFjAVRfI9m2bZsGobYfNYML\nGGPtU81qoCdPnoSjo2Oj6pw4caLR+7JqahPPsmXLYGhoiD///BMA4ODg0OBw5PaMBxcw1r7VjIQd\nOnQoLly40Kg6Pj4+jd6XVVObeM6ePYtXXnkF2traAIBOnToJv3c0PLiAsfavtXQL1kwa2h6pTTzu\n7u7466+/hMdff/01fHx8HmtQrRUnHsbalytXriAsLAxmZmaYPXu28GUvl8thaWlZ536zZs1CcHAw\n/v3vf9e5r42NDTZv3gxvb29YWVkhMjIS5eXlwna5XI4JEybAwcEBa9asEW71AAAtLS188cUXcHNz\ng6OjI44fPw4LCwts3rwZvXv3hlQqxU8//YSEhAR4enrCyckJX3755eN+mZqd2sQzf/58zJ07F1eu\nXIG9vT02b96M1zvo+gCceBhrXyZMmICnnnoKqampkEgk2Lt370OzUtfsZ2RkhNTUVDg7OyM2Nrbe\n2apFIhE2b96MjRs34tixY/j888+RkJAAAMjMzMQLL7yAyZMnQy6XIzk5GW+88YZK/f/+97/YuXMn\nlEoliAh5eXlIS0vD2bNn8dJLL2HKlCnYsGEDYmJisGnTJoSFhakktrZA7cwFEokEBw8eRF5eHior\nK2Fubt4ScbVKhobAjRtPOgrG2h/R8iasW10Himh8N1lubi6USiVOnjwJXV1d/N///R82btyodr83\n3ngD69evb7DtadOmwdPTEwAwatQo/Pjjj3jmmWfwzTffYPTo0Zg4cSIA4N1338WgQYNQVVUFLa3q\n84DZs2ejX79+QltVVVWIiIhA9+7dERISgmXLlmHGjBmws7ODnZ0devbsif/9738qMwm0dmoTzyef\nfIIpU6YIc6fdunULX331FV555ZXHHlxrw2c8jD0emiSM5pKYmAh7e3vo6uoKZe7u7nXu5+DgoHa/\n2vr37y/8bmZmJqxO+ssvv6jcwG5vb4+KigoolUr07dsXAB6aWcXc3Bzdu3cHAIjFYgBQSUxisbjN\nLRSntqtt69atKrNNGxsbY8uWLY81qNaKR7Ux1n54enri4sWLKCkpEcqSk5Pr3O+PP/5AaWlpg/s1\nxpAhQ3D27Fnh8R9//AFtbW1hlhWgegBXe6c28Tz11FMqF78KCgqgp6f3WINqrfiMh7H2w8zMDM7O\nzoiIiMCNGzfw0UcfCXNE1h7ZVrNfZGQk8vPzsWHDBuTk5DTpmGPHjkV8fDy+/vprXLt2DREREfD3\n9xe62ToKtc82ODgYgYGBiImJwVdffYWgoCDMmDGjBUJrfTjxMNa+xMbGoqCgAC4uLrhw4QICAwMB\nPLzsdWxsLPLz8+Hs7Izz589jzJgx6Natm7C9oWWxa7dla2uL2NhY7Nq1CzKZDK6urvjoo48abOfB\nsoaO1VaonTKHiHD8+HHExsYCACZOnAiZTNYuM7S6aR+Sk4FZs6r/ZYw1TnubMoeIYG5ujkOHDqm9\n1tPeNNeUOWo7E0Ui0UNrL3RUfMbDWMeUkJCAPn36oHPnzvjkk09QVVXV4ZJOc1J72nLu3DmEhobC\nyckJvXv3Ru/evWFra9sSsbU6PLiAsY4pIyMD/fv3R58+fXD9+nXEx8c/6ZDaNLVdbTKZDGFhYRg+\nfDh0dHSE8qeffvqxB9fS1J0u3r0L9OrFyYcxTbS3rraOrMW62oqLixEUFNRh52erTV8fuHcPqKoC\n2uElLsYYaxFqvz79/f0xe/ZsHDt2DMnJycJPfUJDQyEWi4WboWr7z3/+Ay0tLRQUFAhlGzduhIOD\nA6RSKU6ePCmUp6enw93dHba2tli6dKlQXl5ejpkzZ8La2hq+vr4qwxpjY2MhkUggkUiwb98+oTw7\nOxsymQzW1taYNWsWKisr1T3tOmlrA7q6QK1h/4wxxjSlbqU4mUxGvr6+D/3UJyEhgZKTk8nFxUWl\nPCsri0aNGkU2NjZ08+ZNIiLKzc0liURCV65cIblcTm5ubsL+o0ePppiYGMrPz6chQ4bQmTNniKh6\nldEJEyZQcXExrV69mubOnUtERJWVlWRra0vnz5+n3377jezs7IS2wsPD6YMPPqCioiJ64YUXKDY2\nts7YG/FykKkpUU6O2t0YY39rzP8r1jbU915q+h6r7WqTy+UaJTIfHx9h7Z7aFixYgDVr1mDcuHFC\nmUKhgJ+fH6ysrGBlZQUiQlFREQwMDJCRkSGMqQ8ICIBCocCAAQOgUCgQHByMrl27IiwsDKNGjQIA\nKJVKuLi4wMXFBQAglUqhVCrh7OyMxMRErFq1Cvr6+ggODsapU6eEuZI0VTPA4O+ZKxhjjGmoUVcq\nMjIy8PHHH2PFihXCjyYOHDgACwsLuLq6qpQnJibCyclJeCyRSKBQKHDx4kVhbjigOomcPn1aqCOV\nSgEAJiYmyM3NRWlpKRQKhVBeu05JSQny8vJgZGQEAHBychLaagoeUs0YY49G7RnPqlWrcPr0aSQn\nJ2PSpEk4cOAAnnvuuUYf4N69e1i1ahV+/PFHoYz+Hv1A9YyOeBARCeVEpFKvrjZqt6XpaIvIyEjh\n97ruX+LEw1jHc+LECcyePbtRK41qsm9bJZfLNe4NU6GuL27AgAFUUVFBUqmUiIiuXr1KQ4YMabBO\nZmamcI0nJSWFTE1NycbGhmxsbKhTp05kbW1NOTk5dPDgQZo3b55Qr1+/fnT37l0iIurdu7dQ/uGH\nH9Inn3xCREQLFiygr7/+moiIbt68SR4eHsJxxo4dK9Tx9/en1NRUIiJyd3engoICIiLat28fLVy4\nsM64G/GtThTzAAAgAElEQVRy0KhRRIcPq92NMfa3xvy/ehL09fXJwMCADAwMSCQSkZ6envD4yy+/\nbLbjWFtb07Fjx5qtvSepvvdS0/dYbVebSCSCtrY2HB0dkZqaim7duqmMSlOnb9++yM3NRWZmJjIz\nM2FhYYHk5GSIxWIMHDgQ8fHxyMrKglwuh5aWFgwNDQEAjo6OiImJQX5+PuLi4oSpwr28vLB7924U\nFxdjy5YtwhTjUqkUqampOH/+PFJSUoTrOzV1tmzZguLiYkRHR6tMS64pPuNhrH0oKipCYWEhCgsL\nYW1tje+++054PHnyZGG/R12CurXcx9SqltJWl5mWL19OBQUF9MMPP5BEIqGePXvSxo0b690/KCiI\nzM3NSUdHhywsLGj79u0q23v37i2MaiMiWr9+PdnZ2ZGTkxMlJCQI5Uqlktzc3MjGxoaWLFkilJeV\nlVFISAhZWlqSTCaj7OxsYduePXvIwcGBHBwcaO/evUL5tWvXaNiwYWRpaUmhoaFUUVFRZ+yNeDlo\n+nSiB54SY6wBjfl/9aTZ2NgIZyU///wz9erVi6KiosjBwYGmTZtGP//8M1lYWAj7W1tbU1RUFA0a\nNIgsLS0pIiKCysrKhPo1+wYHB5OWlpZwNrVmzRoaM2YMffzxxyrH79u3L33zzTcPxZWZmUkikYj2\n799Pjo6O1LdvX9q1a5ewvby8nHbt2kVeXl40aNAg2r17N5WXl9f5PKZOnUqRkZEUFBRE4eHh1KNH\nD/Lx8aHr16/TJ598QnZ2djR8+HBSKBT1vk71vZeavsca7V1VVUUlJSUaHaAtacyLN3cuUQN5lzH2\ngLaYeDp16kShoaGUnZ1NJSUlDyUeGxsb6tevHyUmJtLvv/9ONjY2dPToUaH+g/vW7mrbu3cveXl5\nCY9//fVX6t69u5AwaqtJPIGBgZSVlUXx8fHUpUsX4Xt4+/bt5OrqSmfOnKGkpCTq378/7dixo97n\nERERQTo6OrRjxw7Kz8+nsWPHkpOTE82fP5/y8/Np5cqVNGLEiHpfp+ZKPPUOLti/f79wiljXBf+A\ngIDHcgbW2nFXG2OPQXNN9d9MXVqVlZWIjIyEmZlZvfvUt7y1Ov7+/pgzZw4uXboEOzs77Nq1C0FB\nQQ0uALdo0SJYWlrC0tISNjY2OH78OEaNGoVvvvkGc+fOxYABAwAAc+fORVxcnLB0TV3Pw9HRUdj+\n4osvYsaMGfjll19gZGSEkJAQvPfee7h//z66dOmi9rk0Vb3P9Ntvv4VIJMLt27dx5MgReHl5QSQS\n4fTp0xg9enSHTTyGhpx4GGt2reAaSG1isRiWlpYN7lN7eWtzc3NcvHixUW3r6urixRdfxK5duxAR\nEYGYmBjs379fo2Ndu3YNQPVS2suWLRO2eXh4YNGiRQ0+j9q3tYjFYlhYWAi3m4jFYty/fx83b95E\nz549G/V8mqLewQU7d+7Ejh07cO/ePSQlJUEul+Pnn39GcnIyiouLH1tArR2f8TDW/mm6/DQ1kDi1\ntbUf2j59+nRER0fj6NGj6Nq1qzB4SlMPLqV99uxZDBs2THj84PNoLYvIqR3Vlp2dDQsLC+Fxr169\nkJ2d/ViDas14aQTGmCY8PDyQlJSkUubt7Q2RSIT/+7//w7Rp05rc9rhx47B582YkJSXh3Llz2Lx5\nM8aPH1/v/g0lyJakNq3Pnj0bfn5+mDhxIogIcXFxCAsLa4nYWiU+42Gs/WvMEtQPbqu9vfbvL7/8\nMiIiIvD+++/jnXfewYIFCwBUXyNatmwZDhw4oHEsNYKDg6GtrY1XXnkFIpEI8+fPR1BQUL11H4yz\nvn0eN7Xr8QBAcnIyDh8+DJFIhNGjR8PNze2xB/YkNGa8/fffA59+Chw61EJBMdbGtZb7WFqbL774\nAv/973+RkJDwpENptBZbjwcA3N3deZnXv/HgAsbYo7p37x7WrVun8byX7YXaazxHjx7FiBEjYGRk\nBENDQxgaGuKpp55qidhaJe5qY4w9ivj4ePTs2RPu7u7w9/d/0uE8EWq72gYMGIANGzbA29sbWu18\n2c3GnC7+/jswZgzwxx8tFBRjbRx3tbUfzdXVpjaT6OjowMPDo90nncbiMx7GGHs0aq/x+Pj4YPz4\n8Zg0aZJwk5FIJOqwN5By4mGMsUejNvHk5ubCzMwMJ0+eVCnvqIlHXx+4dw+oqgL4JJAxxjTXqOHU\nHUVj+yn19YG8vOp/GWMNMzExwa1bt550GKwZGBsb17ksTrMPpw4JCXnoAACwffv2Rh+kvamZvYAT\nD2PqabJ+F+sY1CaeMWPGCMnm5s2b2LNnDzw8PB57YK0ZX+dhjLGmU5t4Jk6cqPJ4ypQpGDFixGML\nqC3gxMMYY02n8eXxlJQUVFVV1bs9NDQUYrEYffv2FcrefPNNODk5wd3dHfPnz0dJSYmwbePGjXBw\ncIBUKlUZwJCeng53d3fY2tpi6dKlQnl5eTlmzpwJa2tr+Pr6IicnR9gWGxsLiUQCiUSCffv2CeXZ\n2dmQyWSwtrbGrFmzUFlZqenTVsGJhzHGmk5t4jEwMBBmLDAyMsLSpUuxevXqevcPCQnBkSNHVMqe\nffZZKJVKnD17FsXFxfjyyy8BAHl5edi0aROOHTuGqKgozJs3T6izcOFCLF68GGfOnMHx48eFqb/j\n4uJw584dpKenw8/PD++++y4AoKqqCkuWLMH+/fsRGxuLJUuWCG2tXLkSY8aMQVpaGgoKChAXF6fB\nS/QwnjaHMcaaTm3iKSoqQmFhIQoLC3H79m389NNPGDlyZL37+/j4wNjYWKXsX//6F7S0tKClpYVR\no0bh+PHjAACFQgE/Pz9YWVlBJpOBiFD09zd6RkYGAgMD0b17dwQEBEChUAh1goOD0bVrV4SFhQnl\nSqUSLi4ucHFxgaurK6RSKZRKJQAgMTERYWFh0NfXR3BwsFCnqfiMhzHGmq5RXW1JSUl47733sGrV\nKiQnJz/SAbdu3SrMT5SYmAgnJydhm0QigUKhwMWLF2FqaiqUS6VSnD59WqgjlUoBVA/TzM3NRWlp\nKRQKhVBeu05JSQny8vKEm1+dnJyEtpqK1+RhjLGmUzu4YMOGDdizZ49ww+i8efMwadIkvP766xof\nbMWKFTA0NMSkSZMA1L0oUV1rQRCRUE5EKvUaGjtes/aEJuPLIyMjhd99fX3h6+v70D58xsMY68jk\ncjnkcnmT66tNPFu3bsUvv/wizEgdFhaGwYMHa5x4du7cifj4eBw7dkwo8/LywtGjR4XHFy5cgKen\nJwwNDZGbmyuUp6WlCUvDenl5IS0tDRKJBAUFBRCLxdDV1YWXlxe+/fZblTpTp06Frq4uTE1NcevW\nLRgbG6u0VZfaiac+nHgYYx3Zg3+UL1++XKP6arva7O3t8UetqZgvXboEe3t7jQ5y5MgRrF27FgcP\nHoSurq5QPnDgQMTHxyMrKwtyuRxaWlowNDQEADg6OiImJgb5+fmIi4tTSTy7d+9GcXExtmzZgkGD\nBgGo7lpLTU3F+fPnkZKSAqVSCWdnZ6HOli1bUFxcjOjoaKFOU/HgAsYYewRUj+eff56ef/55GjFi\nBGlpaVHfvn2pb9++pKWlRSNGjKivGgUFBZG5uTl17tyZLCwsaNu2bWRvb09WVlbUv39/6t+/P4WH\nhwv7r1+/nuzs7MjJyYkSEhKEcqVSSW5ubmRjY0NLliwRysvKyigkJIQsLS1JJpNRdna2sG3Pnj3k\n4OBADg4OtHfvXqH82rVrNGzYMLK0tKTQ0FCqqKioM/YGXg4V69cTzZvXqF0ZY6zda+x3Z41652pr\nqP9OJBJBJpM9nkz4BDX2etC2bcD//gd04FmDGGNM0GxztdV1UZ1V42s8jDHWdGqv8fz222946aWX\n8PTTT6NTp07Q0tLq0EtfA5x4GGPsUagd1bZy5Uq8+eabSE9PR0ZGBrZu3YqysrKWiK3V4sEFjDHW\ndGrPeC5fvgwvLy9oa2tDX18fixcvxt69e1sitlaLz3gYY6zp1J7xGBgY4P79+3jmmWcwd+5cWFtb\no2fPni0RW6vFMxcwxljTqV2B9M8//4RYLIZIJEJMTAyuXbuG6dOnw8LCoqVibDGNHZlx/Trg4QFk\nZ7dAUIwx1sppOqqNl76upbEv3t27QK9efNbDGGOA5olH4/V4WPWS1/fuAQ0sS8QYY6wenHiaQFsb\n0NUFaq1nxxhjrJHUJp7Y2NhGlXU0PMCAMcaaRm3iWbVqVaPK2ruLBRcx45sZwmMeUs0YY01T73Dq\nw4cP49ChQ7h27RrmzZsnXDi6ceNGhxxO3UW7C45e/mcJB048jDHWNPUmnp49e8LDwwMHDhyAh4eH\nkHhsbGzg7e3dYgG2Fqb6psgrzkMVVUFLpMWzFzDGWBPVm3j69euHfv364aWXXoKOjk5LxtQqdenU\nBQY6BigoKcDTXZ/mMx7GGGsitTMXXL9+HZ999hni4+Nx69YtANVjti9fvvzYg2ttzAzMkFuUKyQe\nHlzAGGOaUzu4ICIiAm5ubqioqEBcXByee+45hIWFtURsrY7YQIycohwAfI2HMcaaSm3iSUlJwYsv\nvgiRSARnZ2esX78eX331Vb37h4aGQiwWo2/fvkJZYWEhxo0bBysrK4wfPx5Ftb6xN27cCAcHB0il\nUpw8eVIoT09Ph7u7O2xtbbF06VKhvLy8HDNnzoS1tTV8fX2Rk5MjbIuNjYVEIoFEIsG+ffuE8uzs\nbMhkMlhbW2PWrFmorKxsxEvzMDMDM+QW5wLgxMMYY02lNvHo6emhsrISMpkMq1atwldffQUDA4N6\n9w8JCcGRI0dUyqKiomBlZYU//vgDFhYW2Lx5MwAgLy8PmzZtwrFjxxAVFYV58+YJdRYuXIjFixfj\nzJkzOH78OM6ePQsAiIuLw507d5Ceng4/Pz+8++67AICqqiosWbIE+/fvR2xsLJYsWSK0tXLlSowZ\nMwZpaWkoKChAXFycBi/RP8T6/5zx8OACxhhrGrWJZ/369bh37x7eeecdEBFOnDiBqKioevf38fGB\nsbGxSlliYiJmzpyJLl26IDQ0FAqFAgCgUCjg5+cHKysryGQyEJFwNpSRkYHAwEB0794dAQEBKnWC\ng4PRtWtXhIWFCeVKpRIuLi5wcXGBq6srpFIplEqlcPywsDDo6+sjODhYqKMpMwMz7mpjjLFHpDbx\nDBw4EIaGhjA1NUVkZCS2bt0KV1dXjQ5y5swZODo6AgAcHR2RmJgIoDqJODk5CftJJBIoFApcvHgR\npqamQrlUKsXp06cBVCcRqVQKADAxMUFubi5KS0uhUCiE8tp1SkpKkJeXByMjIwCAk5OT0JamHuxq\n48EFjDGmObWj2jIzM7Fr1y6cOnUKpaWlAKpHtf3000+NPohG02WLRHXWryknIpX2GmpbJBJpPGtq\nZGSk8Luvry98fX2Fx7W72viMhzHWUcnlcsjl8ibXV5t4XnvtNXh7e2PZsmXo3LkzgLqTQ0M8PT2R\nnp4ONzc3pKenw9PTEwDg5eWFo0f/mQ3gwoUL8PT0hKGhIXJzc4XytLQ0eHl5CXXS0tIgkUhQUFAA\nsVgMXV1deHl54dtvv1WpM3XqVOjq6sLU1BS3bt2CsbGxSlt1qZ14HlQznBrgxMMY67ge/KN8+fLl\nGtVX29WWlZWFpUuXwtvbGwMGDMCAAQPg4eGh0UG8vLywfft2lJSUYPv27Rg0aBCA6m68+Ph4ZGVl\nQS6XQ0tLC4aGhgCqu+RiYmKQn5+PuLg4lcSze/duFBcXY8uWLUJbUqkUqampOH/+PFJSUqBUKuHs\n7CzU2bJlC4qLixEdHS3U0VTt4dQ8uIAxxpqI6nHz5k3Kz8+n1atX07Jly+jixYt08+ZN4ac+QUFB\nZG5uTjo6OmRhYUHbt2+nu3fv0tixY8nS0pLGjRtHhYWFwv7r168nOzs7cnJyooSEBKFcqVSSm5sb\n2djY0JIlS4TysrIyCgkJIUtLS5LJZJSdnS1s27NnDzk4OJCDgwPt3btXKL927RoNGzaMLC0tKTQ0\nlCoqKuqMvYGXo/rYFWXUaUUnqqisoOPHiXx8GtydMcY6BHXfnQ+qdwVSGxubBrvUMjMzH1MqfHIa\ncz2ox9oeSA1PxbXfxZg5Ezh3roWCY4yxVkrTa+n1XuP5888/AQClpaXQ1dVV2VYzyKAjqhlSbWAg\n5q42xhhrArXXeAYPHtyoso6iZkg1Dy5gjLGmqfeMJzs7G9evX8e9e/eQnJwsDGnOy8tDly5dWjLG\nVqVmSLW3KScexhhrinoTzw8//ICdO3fi2rVrWLhwoVBubW2NlStXtkhwrVHNkOquXYF794CqKkBL\n7XkjY4yxGvUmnunTp2P69OnYt28fJk6c2JIxtWpifTFyinOgrQ3o6gIlJYC+/pOOijHG2g61f6tz\n0lH14HxtPG0OY4xphjuJNCQ2EPPsBYwx9gg48Wio9hkPz17AGGOaUztXGwD89ddf+OWXX1QmCZ02\nbdpjDay14sXgGGPs0ahNPEuXLsXBgwcxePBg6OjoCOUdNfF01+uO26W3UV5ZDgODzpx4GGNMQ2oT\nT1xcHM6dO9eh792pTVtLG093fRo37t2AgUFPHlzAGGMaUnuNx9XVVZg+h1WruYmUu9oYY0xz9Z7x\n+Pv7AwBKSkrQt29fDBw4UFjSWiQS4eDBgy0TYStUM8CABxcwxpjm6k08NbMV1DXrqKYLwbU3NbMX\n8BkPY4xprt7EU7O63KJFi7BmzRqVbYsWLYJMJnusgbVmNQvCceJhjDHNqb3G8+OPPz5UVnu56o7I\nTP+fGap5cAFjjGmm3sQTFRWFvn37IiMjA3379hV+rK2tMXLkyCYdbOvWrRg8eDA8PDwwf/58AEBh\nYSHGjRsHKysrjB8/HkW1TiE2btwIBwcHSKVSnDx5UihPT0+Hu7s7bG1tsXTpUqG8vLwcM2fOhLW1\nNXx9fZGTkyNsi42NhUQigUQiwb59+5oUfw0+42GMsUdQ39Kkt2/fpszMTAoMDKQ///yTMjMzKTMz\nk+7du9ekpVFv3rxJNjY2VFRURJWVlTR69Gg6cuQIffDBB/Tqq69SaWkpzZ07l9auXUtERLm5uSSR\nSOjKlSskl8vJzc1NaGv06NEUExND+fn5NGTIEDpz5gwRVS99PWHCBCouLqbVq1fT3LlziYiosrKS\nbG1t6fz58/Tbb7+RnZ1dnTE28HKoOHb5GMl2yGjPHqJJk5r0cjDGWLvR2O/OGvWe8YhEItjY2ODT\nTz+FoaEhnnrqKTz11FMoKSlBQUGBxglOT08PRIQ7d+6gpKQE9+7dg5GRERITEzFz5kx06dIFoaGh\nUCgUAACFQgE/Pz9YWVlBJpOBiISzoYyMDAQGBqJ79+4ICAhQqRMcHIyuXbsiLCxMKFcqlXBxcYGL\niwtcXV0hlUqhVCo1fg41xPpiXgyOMcaaqN7EM3nyZACAh4fHQz8DBgzQ+EB6enqIioqCjY0NzMzM\nMGTIEHh5eeHMmTNwdHQEADg6OiIxMRFAdRJxcnIS6kskEigUCly8eBGmpqZCuVQqxenTpwEAiYmJ\nkEqlAAATExPk5uaitLQUCoVCKH+wTlP8s/w1Jx7GGNNUvaPavv/+ewBotptHb9y4gfDwcKSlpcHY\n2BiTJk3Cd99999BQ7YbUNYyb/l4Zteb32u011HZ9Q8IjIyOF3319fYXRfbUZ6xmjuKwYOnr3UVjI\nMzowxjoWuVwOuVze5Ppqp8wJDg6GTCaDj4+PcGbSFImJiRg0aBDs7e0BAJMmTcKJEyfg6emJ9PR0\nuLm5IT09HZ6engAALy8vldFzFy5cgKenJwwNDZGbmyuUp6WlwcvLS6iTlpYGiUSCgoICiMVi6Orq\nwsvLC99++61KnalTp9YZZ+3EUx8tkRZM9U1R2ikXRUVWGr8WjDHWlj34R/ny5cs1qq92OHVoaCiu\nX7+O1157Db1798aECROwfv16jQP18fHB2bNnUVBQgPv37+Pw4cN49tln4eXlhe3bt6OkpATbt2/H\noEGDAAADBw5EfHw8srKyIJfLoaWlBUNDQwDVXXIxMTHIz89HXFycSuLZvXs3iouLsWXLFqEtqVSK\n1NRUnD9/HikpKVAqlXB2dtb4OdRmZmCGEq1c7mpjjDFNNWYEQnl5OZ06dYree+89srS0pD59+mg2\n5OFvO3bsoGHDhtGAAQPonXfeocrKSrp79y6NHTuWLC0tady4cVRYWCjsv379erKzsyMnJydKSEgQ\nypVKJbm5uZGNjQ0tWbJEKC8rK6OQkBCytLQkmUxG2dnZwrY9e/aQg4MDOTg40N69e+uMr5EvBxER\nPRf9HO05d5AMDDR5BRhjrP3R5LuTiEj0d6V6PfPMMyguLoa3tzeGDh0KHx8flYv77Uld0wPVZ+bB\nmRjYcxBe8ZqN8nJAi5fUY4x1UJp8dwKNnJ26c+fOSE1NRUpKClJTU1FSUvJIQbYHZgZmuHEvF7q6\nwL17TzoaxhhrO9QOLli3bh2A6hkGdu7ciZCQEOTk5OD+/fuPPbjWTKwvxu83fxeGVBsYPOmIGGOs\nbVCbeD7++GOcOHECSUlJ6N27N0JDQ+Hj49MSsbVqZgZmOJF1gpdGYIwxDalNPKWlpVi4cCHc3d3R\nuXPnloipTeDF4BhjrGnUJp4333yzJeJoc2pmLxBz4mGMMY3wWKwmEhuIhcXgeGkExhhrPE48TdSt\nSzeUVZZB1/Aen/EwxpgG1CaeoqIiVFZWAgByc3Pxyy+/PPag2gKRSAQzAzN0MuLZCxhjTBNqE8+w\nYcNw//593L17F15eXlixYoWwiFtHJzYQAwY5nHgYY0wDahNPVVUVunbtii+++AKhoaE4cuQIn/X8\nzczADNSVz3gYY0wTahNP9+7dcezYMXz++ecIDg4GAJ654G9ifTEqdHN4cAFjjGlAbeL5z3/+g127\ndmHWrFmwtbXFpUuXMHz48JaIrdUzMzBDmQ53tTHGmCYavI+nsrIS0dHR2Llzp1BmZ2eHjRs3Pu64\n2gQzAzOUaJ/nxMMYYxpo8IxHW1sbCQkJKOS+pDqJ9cW4J+IzHsYY04TamQuGDBkCf39/TJw4ET17\n9hSWmg4ICGiJ+Fo1MwMzFIEHFzDGmCbUJp6CggLY2NggKSkJSUlJQjknnurh1Lcr+IyHMcY00vxr\n0dWvqKiIpk2bRg4ODuTk5ESnT59ucAXSDRs2kL29PTk5OdGJEyeE8rS0NHJzc6PevXvT22+/LZSX\nlZVRaGgoWVlZPbQC6d69e6lPnz7Up08fio2NrTM+TV+OwvuFpLtSj/r1r9KoHmOMtSeafneqHdWW\nmZmJ8PBwuLm5AQBSUlLw7rvvNinJRUREwMrKCikpKUhJSYGjoyOioqJgZWWFP/74AxYWFti8eTMA\nIC8vD5s2bcKxY8cQFRWFefPmCe0sXLgQixcvxpkzZ3D8+HGcPXsWABAXF4c7d+4gPT0dfn5+QpxV\nVVVYsmQJ9u/fj9jYWCxZsqRJ8T/IQMcAIohQWMqnPIwx1lhqE09kZCT8/f2Fx3379sVXX33VpIMd\nPXoUb7/9NnR1ddGpUyd069YNiYmJmDlzJrp06YLQ0FAoFAoAgEKhgJ+fH6ysrCCTyUBEKPq7Tysj\nIwOBgYHo3r07AgICVOoEBweja9euCAsLE8qVSiVcXFzg4uICV1dXSKVSKJXKJj2HB5l2NcNdymmW\nthhjrCNQe43n999/x3PPPYelS5cCqD570NHR0fhAV69eRWlpKcLDw5Geno6AgADMmzcPZ86cgaOj\nIwDA0dERiYmJAKqTiJOTk1BfIpFAoVDA2toapqamQrlUKkV0dDTmzp2LxMREzJkzBwBgYmKC3Nxc\nlJaWQqFQQCqVqtQ5ffo0nJ2dH4ozMjJS+N3X1xe+vr4NPi8zQzOkIBeAg6YvCWOMtUlyuRxyubzJ\n9dUmnqFDhwqDCu7fv4+oqCiMGjVK4wOVlpbi999/x9q1azFy5EjMmTMHe/fuRXX3YOOIRKKHyujv\nUXY1v9dur6G262oLUE08jWFuKEZipxxUVQFaPNc3Y6wDePCP8uXLl2tUX+1X5fz587Fp0ybk5OTA\n1tYWSqVS5XpLY9nb20MikcDf3x96enqYPHkyjhw5Ak9PT6SnpwMA0tPT4enpCQDw8vJCWlqaUP/C\nhQvw9PSEvb09cnNzhfK0tDR4eXk9VKegoABisRi6uroPtVW7zqMyN6yeofrevWZpjjHG2j21iadX\nr17Ytm0b/vrrL2RmZmLr1q3o2bNnkw7m4OAAhUKBqqoqfP/99xg5ciS8vLywfft2lJSUYPv27Rg0\naBAAYODAgYiPj0dWVhbkcjm0tLRgaGgIoLpLLiYmBvn5+YiLi1NJPLt370ZxcTG2bNkitCWVSpGa\nmorz588jJSUFSqWyzm62phDri9HZmIdUM8ZYY6ntaissLMQPP/yAU6dO4f79+wCqu6maMm3Ohx9+\niGnTpqG0tBQjR45EUFAQqqqqEBwcDIlEAnd3d3zwwQcAALFYjPDwcIwYMQI6Ojr47LPPVNoJDg7G\nW2+9haCgIAwYMAAA8MILL+DIkSNwcnKCra0tYmJiAFTPwLB69WpMmDABAPD+++9rHHt9zAzMoN3t\nLCcexhhrJBGpucgyefJk6OnpwdvbGzo6OsI1lenTp7dUjC1GJBJpdM0JAA5cOIBp67fh+MsH0b//\nYwqMMcZaMU2/O9We8SiVSqSkpDxSUO2Z2ECMqq7c1cYYY42l9hpPeHg43nvvPVy+fBkFBQXCD6tm\nZmCGCt1cXpOHMcYaSe0ZT9euXbFw4UJs2rRJuH9HJBLh8uXLjz24tkCsL0ZZ5xwUFhKAuodoM8YY\n+4faxLN8+XKkp6fD2tq6JeJpc/Q660EbusgrvA3A+EmHwxhjrZ7arjZ7e3vo6em1RCxtlj6JkVOY\nq35Hxhhj6s94jI2N0a9fP4wcORJGRkYAmj6cur0yFJkhtzgHgOOTDoUxxlo9tYln9OjRGD16NIB/\nhs4PE2cAACAASURBVMzVN91MR9VN2ww3SviMhzHGGkNt4pkxY4bwe1JSEjw8PB5nPG2ScWcxCop4\nhmrGGGsMjaa1nD179uOKo017Ws8Mdyr4jIcxxhpDo8Sj6V39HUUPPTHuVvEZD2OMNYZGiSciIgIA\nJ6AHmRuaoVjEiYcxxhpDbeL597//Lfw+fvx4VFZWYsqUKY81qLam51NmuKfFXW2MMdYYahPPX3/9\nhdWrVwOoXgguICAAffr0eeyBtSUWxtWzFzDGGFNP7ezUVVVVmDJlClxdXfHTTz/hueeewxtvvNFS\n8bWopsxODQDKC/fh8qUhKpeXQkvEy5AyxjoWTb876008SUlJwv065eXlmDNnDgYPHoxZs2YBANzd\n3Zsh3NalqYknOxvotdEEef/+HU93ffoxRMYYY62Xpt+d9f55vnDhQuFnyZIlMDIyQnp6ulDWFJWV\nlXBzc4O/vz+A6kXmxo0bBysrK4wfPx5FtdYW2LhxIxwcHCCVSnHy5EmhPD09He7u7rC1tcXSpUuF\n8vLycsycORPW1tbw9fVFTs4/XV+xsbGQSCSQSCTYt29fk2JviIEBgCIz5PC9PIwxph61oP/85z/0\n0ksvkb+/PxERffDBB/Tqq69SaWkpzZ07l9auXUtERLm5uSSRSOjKlSskl8vJzc1NaGP06NEUExND\n+fn5NGTIEDpz5gwREe3Zs4cmTJhAxcXFtHr1apo7dy4REVVWVpKtrS2dP3+efvvtN7Kzs6s3vqa+\nHBUVRJjuSz9cPNqk+owx1pZp+t3ZYhckrl69ikOHDmHWrFnCKVliYiJmzpyJLl26IDQ0FAqFAgCg\nUCjg5+cHKysryGQyEJFwNpSRkYHAwEB0794dAQEBKnWCg4PRtWtXhIWFCeVKpRIuLi5wcXGBq6sr\npFIplEplsz43bW1Au8QMWQV8xsMYY+q0WOJ54403sHbtWmhp/XPIM2fOwNGxemJNR0dHJCYmAqhO\nIk5OTsJ+EokECoUCFy9ehKmpqVAulUpx+vRpANVJTCqVAgBMTEyQm5uL0tJSKBQKofzBOs2pS7kZ\nrt7iIdWMMaZOvXO1xcbGYtKkSbh8+TJsbW0f6SDfffcdTE1N4ebmBrlcLpSTJqMg6piYlGpNWEpE\nKu011HZDk5xGRkYKv/v6+sLX17dR8elWinH1Np/xMMbaP7lcrvJdrql6E8/777+PSZMmYcKECTh3\n7lyTDwAAv/zyCw4ePIhDhw6htLQUd+/exdSpU+Hp6Yn09HS4ubkhPT0dnp6eAAAvLy8cPXpUqH/h\nwgV4enrC0NAQubn/nFWkpaXBy8tLqJOWlgaJRIKCggKIxWLo6urCy8sL3377rUqdqVOn1htr7cSj\nCX0yQ05RepPqMsZYW/LgH+XLly/XqH69XW2Ojo7w9fVFZmYm/P39VX7Gjh2r0UH+v707j4uq3B84\n/plhlcUl11IU8ecFBmKRWNRrIpWZiaCWW/LKoMUuXsur9euW3otZesm8mnUr9SdX04wys8w0yhQV\nFyDADbFCNEkRwQ2QnXl+fxwYQAfZmUGet695zZwzZ/nOF5wvzznPec6SJUvIzMzk7NmzREdHExAQ\nwMaNG/H19SUqKoqioiKioqLw8/MDwMfHh5iYGM6fP09sbCxqtRpbW1tdXNHR0eTm5rJt27ZahWfT\npk3cvHmTNWvW6Lal0Wg4efIkJ06c4Pjx46SmpuLi4tKo+BvCVtW78p48kiRJ0p3U2eL59NNPOX78\nOBMnTmT+/Pm1Dl019348Veu/+OKLzJgxA0dHR4YMGUJkZCQAvXv35sUXXyQgIABzc3NWr16tW/fd\nd99lxowZ/P3vf2fq1Kk88MADAEyYMIHvv/8eZ2dnHBwciI6OBsDExISlS5cyadIkQGnJtYYuJn24\nVCQLjyRJUn3qHbkgJyeHnj17Ul5eDoCpab238Gm3mnoBKcC4aVnEuXhyfYEsPpIkdSwtdgFplby8\nPKZMmYKDgwMODg5MnTqVjIyMZgV5N+pu2ZOCiitUaCsMHYokSZJRq7fwLFmyhPHjx5ORkUFGRgZB\nQUG8/fbbbRFbu9LZxpROqm7kFuYaOhRJkiSjVm/hSUpKYtq0aZiammJqasrkyZNJSkpqi9jaFRsb\n6IwdP1/82dChSJIkGbV6C09gYCAvv/wyycnJJCUlMW/ePN1Ya1I1Gxv4c8lbhO8M50bxDUOHI0mS\nZLTq7VyQl5fH+vXr+e677wAYN24cTz/9NJ07d26TANtSczoXrFoF6elQ/uhfKCgt4JMJn7RwdJIk\nScapxW6L0BE1p/BERUFcHLz/8U2GrBnCW6Pe4kmXJ1s4QkmSJOPT4r3apIaxsYGCArA2t2bjhI3M\n3jWbi/kXDR2WJEmS0ZGFp4VUFR4An74+/OWBvxD6TWiTW1CSJEl3K1l4WoiNDeTnV0+/PuJ1rhVf\n46OfPzJcUJIkSUaozmEI6hr0rWq4m3/84x+tE1E7ZWtb3eIBMDMxY+OEjQyPGs5DAx/CsYej4YKT\nJEkyInW2eKytrbGxsan1UKlUrFu3TjemmlSt5qG2Kn/q/ife9H+TGdtmUFZRZpjAJEmSjEyDerXl\n5eWxatUq1q1bx+TJk5k3b16tG7LdLZrTqy0rC4YMUZ5rEkLw+ObH8e7rzSL/xg0dLkmS1B60aK+2\nK1eusGDBAtzd3SkrKyM5OZnIyMi7sug0l74WDyg/kHXj17H659Uc+aPl73wqSZLU3tRZeObPn4+P\njw+2trYcP36cRYsW0a1bt7aMrV2xtobCQtBqb3/vXtt7+fDxDwnZFsLN0pttH5wkSZIRqfNQm1qt\nxtzcHDMzs9tXUqnIy8tr9eDaWnMOtYHS6rl0SXnWZ+bXM7EwtWD1uNX6F5AkSWqH5MgFzdDcwtOn\nDxw9qjzrk1eSh8fHHvz70X8T7BTc5P1IkiQZE6MduSAzM5NRo0bh4uKCv78/mzdvBiA/P5+goCD6\n9+9PcHAwBTVOlKxatYrBgwej0WiIi4vTzU9LS2PIkCE4ODjwxhtv6OaXlZURFhbGgAED8Pf359Kl\n6puybdmyBUdHRxwdHfnyyy9b5TPWdZ6nSmeLznw68VNm7ZglRzWQJKnjEm0kKytLpKSkCCGEyMnJ\nEQMHDhR5eXkiMjJSzJ49WxQXF4vw8HCxbNkyIYQQ2dnZwtHRUfz+++8iNjZWeHp66rb12GOPiejo\naJGbmyuGDx8uEhMThRBCfP7552LSpEni5s2bYunSpSI8PFwIIURFRYVwcHAQJ06cEMeOHRODBg3S\nG2Nz0zF0qBDbt9e/3KLYReKhDQ+JCm1Fs/YnSZJkDBr73dlmLZ4+ffrg4eEBQI8ePXBxcSExMZGE\nhATCwsKwsLAgNDSU+Ph4AOLj4xkzZgz9+/dn5MiRCCF0raFffvmFKVOm0L17dyZOnFhrnRkzZmBl\nZcXzzz+vm5+amoqrqyuurq64ubmh0WhITU1t8c/46qvw979D5V3C6/T6iNcpLi9mxeEVLR6DJEmS\nsTPIkDnp6emkpqbi4+NDYmIiTk5OADg5OZGQkAAoRcTZ2Vm3jqOjI/Hx8aSnp9fqzq3RaDhyROmm\nnJCQgEajAeCee+4hOzub4uJi4uPjdfNvXaclBQVBr16wdu2dlzNVm7Jp4iYiD0aSkpXS4nFIkiQZ\nszqHzGkt+fn5TJkyhRUrVmBjY9O4nhCVw/XUJITQzRdC1Nrenbatb1sAERERutf+/v74+/s3Ij5Y\nuRIeeQSmTYOuXete1r6rPSvHrGT6V9NJej4JKzOrBu9HkiTJkGJjY4mNjW3y+m1aeMrKypg0aRIh\nISEEBQUB4O3tTVpaGp6enqSlpeHt7Q2Ar68vu3fv1q17+vRpvL29sbW1JTs7Wzf/1KlT+Pr66tY5\ndeoUjo6OXL16ld69e2NpaYmvry/ffvttrXVCQkL0xliz8DSFmxsEB8Obb8K//33nZaffP51d6buY\n98M8PnpcDiYqSVL7cOsf5XWN7VmXNjvUJoQgLCwMV1dXXn75Zd18X19foqKiKCoqIioqCj8/PwB8\nfHyIiYnh/PnzxMbGolarsbW1BZRDctHR0eTm5rJt27ZahWfTpk3cvHmTNWvW6Lal0Wg4efIkJ06c\n4Pjx46SmpuLi4tJqn3XxYvjkE/j11/qX/eCxD4hJj+Gb09+0WjySJElGpaV7N9TlwIEDQqVSCXd3\nd+Hh4SE8PDzErl27RF5enhg/fryws7MTQUFBIj8/X7fOypUrxaBBg4Szs7PYv3+/bn5qaqrw9PQU\n9vb24rXXXtPNLy0tFc8884yws7MTI0eOFFlZWbr3Pv/8czF48GAxePBg8cUXX+iNsSXTsWyZEIGB\nDVv24PmDovey3uJC3oUW278kSVJbaex3p7yAtIbmXkBaU0kJuLjARx8p53zqExEbwcHMg8TMiEGt\nkrdJkiSp/TDaC0g7GgsLePddmDu3/u7VAAseXEBhWSErj6xs/eAkSZIMSBaehiorg5TGdX2u6l69\nZk39y5qqTdk0Qeli/X36900MUpIkyfjJQ2013LG5eOYM+Poqz126NHibx48rh9pOn4aGDO59KPMQ\nwdHBbJ+2Hb9+fg3eT0sqK1MGO718ufqRk1N7uqgI7rsP7OygX7/azz16gFr+SSNJHYYcJLQZ6k1e\nSAg4OUGN8eEa4oUXlNsm1Ne9usrO33YS+k0oe57eg6anpv4VmqGoSCmOyclKgy45GU6dgu7doWdP\npcXWq1ft1716gaUlXLgAf/wBmZm1n/PzlSLk5wejRikPBwflOidJku4+svA0Q73JO30aHnwQMjLq\nvveBHpcvg0YDBw+Co2PD1tl4bCNv7HmDuNA4+nfp3+B91efsWdixA37+WSkyZ84otXTIkOqHmxtY\nNeN61sJC+P135fPu2QN794K5eXURCghQWkaSJN0dZOFphgYlb8oU8PaG+fMbte1ly2D/fqhxHWu9\nVhxeweqk1cSFxtHDqkej9ldFCKUF89VXsG2b0iIZNw6GDVOKjIuL0hGiNQkBv/yiFKA9eyA2Vjla\n+eijMHEijBwJpm0+hoYkSS1FFp5maFDyTpyA0aOVVk+nTg3edlX36qVL4YknGn7Y6fWfXmd3xm72\nPL0HG/OGtbKEUFo0X32lPIqKYMIE5Uv+z38GE5MGh90qtFpITYXvvlPiy8iA8eOV+B5+WDmMJ0lS\n+yELTzM0OHkTJijHjObMadT29+2D0FDl/MmcOfDkk/W3NoQQPL/jeX6//jvfTvsWC1P9K1y/rrQm\nfvgBdu5UauKkScqXuZeXcZ9fOX8evv5aKUJHj8KYMUrcjz0GlYNVSJJkxGThaYYGJy8pSekrfeZM\no49TVVQohWHVKjh5EmbNUjof1HXXUoBybTmTt0zG3MScTyd+ionahLIyiI9XCs2PPyrbGj5caYyN\nGQPOzsZdbOpy+TJs364UoYMHlRbQ1Knw+OPNO+8kSVLrkYWnGRqVvLFjleNDs2Y1eX+pqfDBBxAd\nrZx3mTNHOX1Upbi4uhvzH5eKee3kGLqWudAz8QP271Ph4KAUmkceUYrO3XaI6to15bxUdDQkJCgp\nnzpVOTfU2uelJElqOFl4mqFRyTt8WLn3wW+/gZlZs/Z77RpERSlFyNJSGeng8mWl8NTsyty1zw32\n9g9gYBcH/m/S+7ja36GZdJe5fBm2blWK0IkTSoNz6lSlh1wz0y9JUjPJwtMMjR6r7eGHYfp05cRN\nC6ioUM5xdO6sFJrOnW8/XFZUVsSb+99kXfI63nnkHZ52f7rOewvdrS5cgC1b4PPPld5yY8Yoheix\nx5ScSZLUtmThaYZGF559++DZZyEtrc37A6dkpRC2PYweVj1YE7gG+672bbp/Y3HxotJFfft2OHAA\nhg5VilBgoLxWSJLaiiw8zdCk0akffFDpHfDUU60T1B2Ua8tZfmg5yw4tY+GDC5ntMxsTtYH7ShtQ\nQQHExMA33yhdte3tlfNCI0YoBUn2kJOk1iELTzM0qfD8+CO89JLSrcxAA5T9euVXnvv2OUrKS1g3\nfh0uvVrvJnftRXk5xMUpP564OKUjoqOjch1T1ePeew0dpSTdHWThaYYmFR4hlD+n581TLswxEK3Q\nsjZpLQv2LmD6/dOZ4DSB4XbDMTORZ95BuYA3KUk5HBcXp3TVvucepTfgAw8oozi4uzdqJCRJkirJ\n+/HUYf/+/Tg7OzN48GDef//9ltuwSgULF8JbbymX5LeB2NjY2+apVWpeeOAFjr5wlG6W3Xjlx1fo\n9W4vntzyJOuPrie7ILtNYqsvTkOxsFCGCfrf/1XOCeXmKhetDhsGP/wQy9y5SocOJyelv8i778JP\nP8HVq4aOXGFMubwTGWfLai9xNlaHKTwvvfQSq1evZvfu3fznP/8hNze35TY+dqzSuaAxA7E1w51+\nGft27kuEfwSJzyWSFp7G2P8Zy3e/fYfjB454r/UmIjaCI38c4WrR1Ra722pT4jQ0tRpcXZXTc15e\nsSQkwI0b8MUXynVCmZkQEaGcJ+rZU2kRjR8P4eHKsEebNiljzp05owxJ1NqMOZc1yThbVnuJs7E6\nxNCMN27cAODBBx8EYPTo0cTHx/P444+3zA5UKliwABYvVi4ssbIy/IBoQB+bPjzj+QzPeD5DaUUp\nB88f5LvfvmPWjlmcvX6Wcm05dp3tsOtiR/8u/ZXXldPdO3XH2twaKzMrrM2ssTa3xsLE4q7uum1m\nBm73C9xctDw9XQtaLdqyCnIva7lwvoKLmRVc/EPLpQsVHIqvIDtLS/bFCnKytZibarmnm6B7N+X5\nnq5aunVVprt2EVhZKa0uy07QyRIsLFVYWlZPm5kr+zczFZiaVHajr/mHQXa2cv8Klar6oVbXntY3\nX99rtbr265rzqh4mJrWnJakFdYjCk5iYiJOTk25ao9Fw5MiRlis8oPTh/eAD5Yx1UZHSArKyUgZN\ns7KqfpiZ3X5xTmO+zIVQ7jmwZ0/1tL5lbnk2B0YJwSgAOgEaykUFpeWllGgvUFqRQUlFaeWjhHJt\nOYXaCgqEFq22Aq3QIoTARKXGRGWCWqWuDF2FClChfLkp/5TPlHOtlFOfLFfeF1QthUoIZVkhKucr\n7yuvRfW0br6ontaK2+YjBOoa02o981SVedC3nYoKLeVvvYlaq6xToQKhVqFVgVCpsFWr+JMK/ket\nQls5X1vjtehcOa9MhfYyygOoQIUW5bXuR1H1+YFSAaUI8gQI3a+AqMyAEqgyX5BdWMGJL/+v+rNT\n+fmo8Zkqt62CyvfQ5UM3r8a6VdNV75toq5+r3jetjLtCpTy0uoeqxuvq6YIyQda7ixEqFaLqPape\nqxBUP4vK/Nacp61cr+q3uuq1qNwGKNsEVY33qpatsd0a63PLMgDnbpSw97/La71Xc/lbXys/l9r/\nT3XbFvqWr71ffVN1raNQ3vjtWik7Ny6va6HKTTbmj0FVjVcqPSHp35bQ9x3VjAMmHaLwNEaL/UVf\nWqo8rl9vme3dYtH5862y3fpVfZU2zIfXSlsvlBb0Vs3/RQKoEDUmjMPHhW1zDlGvmt/wuhn6p5eX\nGU/O7mTDjfbxu7m5nfwfaowOUXi8vb155ZVXdNOpqamMGTPmtuVkBz9JkqTW1yEO3nbp0gVQerad\nO3eOH3/8EV9fXwNHJUmS1DF1iBYPwMqVK3nhhRcoKytjzpw59OjRtDt6SpIkSc3TIVo8ACNHjiQt\nLY309HTm3HIDt1a7xqeF2dvb4+bmhqenJz4+PoYORyc0NJTevXtz//336+bl5+cTFBRE//79CQ4O\npqCgwIARKvTFGRERQb9+/fD09MTT05Pvv//egBFCZmYmo0aNwsXFBX9/fzZv3gwYXz7ritPY8llc\nXIyvry8eHh74+fmxYsUKwPjyWVecxpZPgIqKCjw9PQkMDASamEshCQ8PD7Fv3z5x7tw54ejoKHJy\ncgwdkl729vbiypUrhg7jNvv37xfJycnC1dVVNy8yMlLMnj1bFBcXi/DwcLFs2TIDRqjQF2dERIRY\nvny5AaOqLSsrS6SkpAghhMjJyREDBw4UeXl5RpfPuuI0tnwKIcTNmzeFEEIUFxcLFxcX8euvvxpd\nPoXQH6cx5nP58uVi+vTpIjAwUAjRtP/rHabFU5ea1/gMGDBAd42PsRJG2AFixIgRdOvWrda8hIQE\nwsLCsLCwIDQ01Chyqi9OMK6c9unTBw8PDwB69OiBi4sLiYmJRpfPuuIE48ongFXlrWsLCgooLy/H\nwsLC6PIJ+uME48rnH3/8wc6dO3n22Wd1cTUllx2+8NR1jY8xUqlUBAQEEBwczPbt2w0dzh3VzKuT\nkxMJCQkGjqhu77//Pn5+fkRGRpKfn2/ocHTS09NJTU3Fx8fHqPNZFWdVhx1jy6dWq8Xd3Z3evXsz\ne/Zs+vfvb5T51BcnGFc+586dy7Jly1DXuKi4Kbns8IWnPTl48CDHjh1j6dKl/O1vf+PSpUuGDqlO\nxvRX2p28+OKLnD17lpiYGM6cOcPq1asNHRKgHDefMmUKK1aswMbGxmjzWTNOa2tro8ynWq3m2LFj\npKen8+GHH5KSkmKU+dQXpzHlc8eOHfTq1QtPT89a+WtKLjt84fH29ub06dO66dTUVPz8/AwYUd3u\nrRzH39nZmfHjx/NtG40N1xTe3t6kpaUBkJaWhre3t4Ej0q9Xr16oVCq6dOlCeHg427ZtM3RIlJWV\nMWnSJEJCQggKCgKMM5/64jTGfFaxt7dn7NixxMfHG2U+q9SM05jyeejQIbZv387AgQOZNm0ae/bs\nISQkpEm57PCFp71c41NYWKhrZufk5BATE6P3Ilhj4evrS1RUFEVFRURFRRltMc/KygKgvLyczZs3\nM3bsWIPGI4QgLCwMV1dXXn75Zd18Y8tnXXEaWz5zc3O5Xjl6yJUrV/jhhx8ICgoyunzWFacx5XPJ\nkiVkZmZy9uxZoqOjCQgIYOPGjU3LZWv0emhvYmNjhZOTkxg0aJB47733DB2OXhkZGcLd3V24u7uL\ngIAAsW7dOkOHpDN16lRx7733CnNzc9GvXz8RFRUl8vLyxPjx44WdnZ0ICgoS+fn5hg5TF6eZmZno\n16+fWLdunQgJCRH333+/8PLyEnPnzjV4r8EDBw4IlUol3N3dhYeHh/Dw8BC7du0yunzqi3Pnzp1G\nl8/jx48LT09P4ebmJkaPHi02bNgghBBGl8+64jS2fFaJjY3V9WprSi7ljeAkSZKkNtXhD7VJkiRJ\nbUsWHkmSJKlNycIjSZIktSlZeCRJkqQ2JQuP1CGdO3eu1mChxmTmzJls3bq1VbZtY2PTKtuVpMaQ\nhUeS2ohW27A7iLbYXXDbeNuS1FCy8EgdXkZGBkOGDCEpKem29/z9/VmwYAEeHh54enqSnp7OE088\ngaurKx9//LFuud27d/Pkk08ydOhQlixZoptvY2PDwoUL8fDwuG0MwLVr1+Lj44OXlxevvvoqpaXV\ntzg+dOgQXl5e+Pv76wZdXL9+PX/96191y4wbN479+/fr9rN48WJcXFyYPn06V69eBeDixYuEhYXh\n5ORUK66CggIefvhhhgwZwtixY9m3bx+gtAQ1Gg3h4eFoNBpmzZpFWVkZoNwKYd68eXh6euLl5cXZ\ns2cB2LJlC+PGjWPEiBGsWbOmCT8BqaORhUfq0H755ReeeOIJNmzYgJeX123vq1QqsrOzSU5OJjg4\nGB8fHyIjIzly5AhLlixBCEFhYSGRkZFs3LiRgwcPcuLECV2xKCwspGfPnhw9epRhw4bV2vakSZNI\nSEjg559/prCwkL179wLKqACJiYnExsbyzjvv8Nxzz+mNvWbrpbCwkPvuu4/U1FSsra3ZsWMHAP/6\n179wcnIiLS2N4uJi3fKdOnVi27ZtJCcn8/HHHxMREaF77/Tp00ycOJGTJ09y7tw5Dh8+DMDChQsZ\nNGgQKSkpHD58mD59+nDu3Dm+/PJLvv76a3766Sc2b96su9pekuoiC4/UYV2+fJng4GA2b958x/M9\n06ZNQ61WM3ToUFxcXBg0aBA2NjbY2dmRmprKrl27OHXqFEOHDsXLy4ujR4/qioharWbmzJl6t5uR\nkcGMGTNwdXVl586dxMTEAEpBmTBhAra2tvj4+CCE4MKFC3c8TGZqaspTTz0FQEBAgK5YxMTEEBoa\nikqlIjQ0VLe8iYkJ7733HsOGDSMwMJDExETdLUL69u3LQw89hFqtZuTIkRw+fJiysjL27t2rK4Lm\n5uZ06tSJrVu3kpCQgLe3N76+vly8eJE9e/Y08CcgdVQd5tbXknSrrl27MmDAAA4cOKAb1j00NJSU\nlBT69u2razV07doVUL5sq15XTZeUlKDVahk9ejT//e9/b9tHp06d6Ny5s979z58/nzfeeINNmzbx\n3nvvcfToUd17tw4oolKpsLS0pKSkRDev6nAagIWFBZaWlgCYmZnVat3oG5wkNjaWAwcOEBMTg7W1\nNb169dIVnls/482bN3XbuXVbWq2WmTNn8s9//lPvZ5QkfWSLR+qwzM3N+eqrr/jkk0/47LPPAIiK\niiIlJUVXdOqjUqkIDAzkwIEDuhF6r169yvnz5+td9+LFiwwePJhr167x2Wef6Vo0Qgi++eYbCgoK\nSExMRK1Wc9999+Hn58eRI0coLS3l5MmTDbrvyZgxY9iwYQNarZb169fr5l+4cIG+fftia2tLdHR0\nrSJ2KyEEZmZmjBo1irVr1yKEoKSkhKKiIqZOncrWrVt1n/fChQvk5OTUG5fUscnCI3VYKpUKKysr\nduzYwYoVK+otNiqVSu/hLktLS9auXcvChQtxc3Nj9OjRunsl3enw2OLFixk3bhyPPvooo0aNqrUf\nLy8vRo4cySuvvMLatWsBGDBgAIGBgXh4eLBo0SL8/f1rraMvztdee41Tp06h0WiwsLDQzQ8ODub6\n9es4OzsTFxeHRqPRu62a02+//Tbp6em4u7szfPhwsrOzsbOzIyIiglmzZuHm5sbkyZMpKCi4lnWM\noQAAAEBJREFUYx4lSQ4SKkmSJLUp2eKRJEmS2pQsPJIkSVKbkoVHkiRJalOy8EiSJEltShYeSZIk\nqU3JwiNJkiS1qf8Hsyxmkpe8Af0AAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 42 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(transcripts_all[:,0], transcripts_all[:,1], label='all data')\n", "plot(transcripts_dn[:,0], transcripts_dn[:,1], label='diginorm')\n", "plot(transcripts_trin[:,0], transcripts_trin[:,1], label='Trinity norm')\n", "title('k-mer abundance plot of normalized data (zoom)')\n", "xlabel('k-mer abundance')\n", "ylabel('N of k-mers with that abundance')\n", "legend()\n", "axis(xmax=40, ymax=20000)\n", "savefig('normhistzoom.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEVCAYAAADKN2OaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXdYFOf2x79LUxCUojSRJkhTEAwgUQSJiS0gUWMDY4sa\n9BfFYElMImCMRpNri7FeSyyxXfXqTdHEgj0goDECGhEVC6AUBRSQ8v7+mDBhlWVhWZid5XyeZx/Y\nKe98Z2d2z5xz3ve8EsYYA0EQBEEoiIbQAgiCIAhxQ4aEIAiCaBRkSAiCIIhGQYaEIAiCaBRkSAiC\nIIhGQYaEIAiCaBRkSBqAra0tTpw4IbQMhdDQ0EBGRkazHEtsn1NzfjbR0dGws7NDz549m+V4yiIm\nJgZjx44FAGRmZsLAwADKHjnQkPtm27Zt8Pf3V+rx66KoqAgODg4oLS1ttmPWRUJCAvr27Su0DB4y\nJA1AIpFAIpEILUPlUdfPqbE/Xvfv38d3332HpKQk/P7770pU1vTUvJ7W1tYoKipS+jVuqvumphFU\nlLVr12LYsGFo3bq1klQ1Dh8fHzx79gwJCQlCSwFAhkSUVFZWCi2BUIBz587ByckJxsbGTX4sZd8j\nLXnc8osXL7Bq1SpMmzZNaClSTJs2DUuXLhVaBgAyJAqTlpYGe3t77N27t9b1MTExGD16NKZNmwZT\nU1P06dMHWVlZ+O677+Dg4ICgoCCpp4mSkhJs3rwZPj4+6N27N/bv389/ebdt24bevXsjOjoaNjY2\niI2NfeV4CQkJ8PPzg5GREfz8/LBmzRpUVFRIbXPmzBl4eHjAzc0Ne/bs4dt/+Yntzp070NDQQFVV\nFQAgMDAQX331Ffr37w8LCwt89NFHePLkCb/9xYsX8dZbb8HW1hbffvttg3RpaGhg165d6N69Oxwc\nHLBixQqp/U+dOoXw8HAYGxvDzc0Nly9fBgAUFBRg5cqVcHNzw8CBA/Hrr7/KuFLA+PHjERkZiaFD\nh8Lc3Bzz5s1DXl5erds+f/4c69atQ7du3fDWW2/hf//7HwDuekdERODixYswMDCQaQzy8/OxdOlS\nODo6Yvjw4Th9+jQAYPPmzZgwYQLi4+NhYGBQ6zWsvs4LFy5Ex44dMWDAAFy8eFFu28A/91tERAQs\nLCywbds2BAYGYvHixXjzzTfRoUMHjBs3DuXl5fjggw9gaWmJCRMmICsri29j5syZsLa2hpmZGaZO\nnYo//vij1nOseX9Ufx7Vr9atW8POzo7f9pdffkFISAicnJywYsUKFBcX8+vqum9epri4GIsXL4aV\nlRXeeOMNZGdnS62Xpf3o0aNYsmQJ9u7dCwMDA3h6egIAtm7dCldXVxgaGiI0NJS/zrWRkpKCqqoq\n2NjY8LplnXNFRQV27tyJnj17ws/PD7t27ZK63+Pi4jBs2DA4Ojpi2bJlKCgo4NdpaGhgx44d6N69\nOzp27IgVK1YgOzsb/fv3h5WVFaKjo6Xa8vX1xcmTJ+v83JoNRtQbW1tbduLECZaUlMSsra3ZTz/9\nJHPb6OhopqOjw7Zu3cpyc3NZSEgIc3FxYZGRkSw3N5d98cUXLCgoiN9+1qxZbNSoUez27dvsypUr\nrGvXruzXX39ljDG2detWpq2tzT755BP25MkTVlJS8srxkpKSWHx8PKuoqGDnz59nNjY27LfffuPX\nSyQS5ufnx65evcpOnz7NbG1t2dGjRxljjMXExLDw8HB+29u3bzOJRMIqKysZY4wFBASwTp06sePH\nj7P79+8zb29v9u9//5sxxlhBQQHT09Nj33//PXvw4AELCwtj2tra7MSJE/XWFRQUxNLS0lhSUhIz\nMDBg6enpjDHGkpOTmampKfvhhx/YixcvWHp6Ort79y5jjLF33nmHzZgxg2VnZ7MzZ84wS0tLdvPm\nzVqvxbhx41ibNm14jaNGjWKjRo2S0nDr1i3GGGMLFixgffv2ZTdu3GAnTpxgtra27NSpU4wxxrZt\n28Z69+4t85ozxth7773HRowYwe7du8cOHDjAjI2N2e3bt+u1/9atW5mOjg6LjY1l+fn5LDo6Wmr7\nutqOjo5m2trabPXq1aykpISVlJSwgIAAZmNjw+Li4titW7eYnZ0dc3NzY99//z1/Ty5cuJBvf+fO\nnSw/P58VFBSwOXPmsF69evHroqOj+Xvk5fujmvLychYQEMDmz5/PGGPs8OHDzN3dnV28eJE9fPiQ\njRgxgl8n7755maioKNa/f3+WkZHBDh8+zMzMzJi/v3+9tMfExLCxY8dKtffTTz+xjIwM9uLFC7Z7\n926mq6vLiouLaz32Dz/8wF5//fVa1718zlu2bGHu7u7s0qVLLCkpiXXv3p1t3bqVMcZYRkYGMzQ0\nZPv372f3799nI0eOZOPGjePbkkgk7M0332Q3b95kJ0+eZJqamiwoKIidPn2av35nzpzhty8tLWWa\nmposKyurVm3NCRmSBmBra8sWLFjArKys2OnTp+vcNjo6mrm7u/Pvd+7cybS0tFhBQQFjjLH79++z\n1q1bs9LSUlZVVcVsbW1ZZmYmv/2KFSvYtGnTGGPcD4yuri4rKyurt9ZPP/2U/d///R//XiKRsM2b\nN/PvP/nkE359zR8Jxl79oQgMDGQffvghv37JkiVs5MiRjDHG9u7dK/WFvnXrFpNIJDJ/EGrTdeDA\nAf59//792bp16xhjjM2dO5fNnDnzlTYKCwuZhYUFe/78Ob8sMjKSLVu2rNZjjhs3TkrjjRs3mJGR\nEX9+NQ2Jh4cHO3bsmJTeGTNmMMa461CXIaioqGAmJibsxo0b/LKwsDC2fPnyeu2/detWKV0PHz5k\n2trarLi4WG7b0dHRzN7eXqq9wMBAXjtjjE2ePJl5eXnx73ft2sUCAgJq1VJYWMjatGnDHj9+zLcv\nz5B88MEHLDg4mH8/ZswYtmvXLv795cuXmaurK2Os4feNq6ur1LqxY8fK/Czr0i6L3r17s//85z+1\nrlu6dCkbPnx4retePueQkBC2YcMG/v2mTZtYSEgIY4yx5cuXs9GjR/Prbt68yUxMTKTuw4MHD/Lr\nHR0d2UcffcS/nzx5MouOjpY6vpmZGfv999/rPLfmgEJbDYAxhg0bNqBXr17o06cPv3zXrl28mzt4\n8GB+ubu7O/+/mZkZrKysYGhoyL8vKytDXl4erl+/jszMTLi7u8PIyAhGRkaIjo7G+fPn+f09PDyg\no6MjU9uDBw/wwQcfwN3dHW3btsWKFStw9epVqW26d+/O/+/p6SkVNpFHzX3Nzc3x4MEDAEB8fDw8\nPDz4dfb29mjXrp3CuiwsLPDw4UMAXBigV69er2g5d+4cHj9+DEtLS/7z2rx5M86dO1erdolEIqWx\nS5cuKC8vR1pamtR2RUVFuHr1Knr06MEv69GjB86ePSv7g6lBWloaysrK0KVLF4X2BwA3NzdoaHBf\nSwsLC1RUVCAnJ6debfv6+r7SXs3zNjMzk7onTU1N+esIcKG1wYMHo0OHDrC2tkZJSQn+/PPPeune\nsGEDzpw5gx9++IFfdvz4cURERPDXqG/fvrhz5w5ycnLk3jc1KSwsRFpa2iv3b00aqv3cuXMYPXo0\nbGxsYGhoiISEhFfuy2psbW35e1LeOV+4cEHm/fPyOgcHB1RUVCAlJYVf9vL1evl9zetVVlaG3Nxc\nPuQmJGRIGoBEIsGGDRtw9+5dfPTRR/zysLAwFBUVoaioCD/99BO/bX1xcnKClZUVUlNTUVBQgIKC\nAjx9+hRXrlzht9HS0qqzjUWLFqG8vBw///wznj59ilmzZvE5jmqq8wsAkJycjNdffx0AYGVlhZyc\nnFq3k4evr6+Uzlu3buHp06cN0iWLvn371moc/Pz80KFDB+Tk5PCfV2FhIQ4fPlxrO4wxKY03btyA\ntrY2XFxcpLYzMDCAu7s7EhMT+WWJiYn8Q4OmpmadSWdnZ2e0atUKN27cqHX/xlCftjU1NRVu/969\ne/joo48wf/583L17F5mZmdDV1a1Xkv3s2bNYsGABDh8+DH19fX55UFAQNm3axF+jgoICPHv2DGZm\nZnLvm5q0bdsWzs7Or9y/1d8xedq1tLSkzoMxhqlTpyIgIADJycl48uQJfHx8ZJ6rs7Mzbt26Va9z\n7tWrl8z75+V1N2/ehKamJtzc3GR8snWTkZEBfX19mJubK7S/MiFD0kAMDAxw9OhRnDlzBp988onM\n7erzBaxGQ0MDI0eOxLx585CWloaqqircunULZ86cqXcbDx8+hLGxMUxMTBAXF4ft27e/ss2WLVtw\n7do1nD17Fnv37sXbb78NgPvC//7770hOTsaNGzfw3Xff1ft83nrrLSQnJ2PXrl14+PAhYmNjpYxe\nfXS9fJzqY40aNQr79u3Dvn378OLFC6SnpyMzMxOGhobo3bs3/8NRWVmJa9euSX1JX+by5cu8xoUL\nF2LAgAH8k39NhgwZgq+//hp//fUX4uLisHv3boSGhgLgni5v3rwplTCuiZaWFgYPHozo6Gg8ePAA\n//3vf3H06FF+/8agaNsv/4DK4vHjx2CMwdzcHEVFRZg/fz7Kysrk6rp37x5GjBiBHTt2wMHBQWrd\n2LFjsWzZMpw7dw6VlZV4/Pgxjhw5AkD+ffMygwYNwjfffIPbt2/jxx9/lBpvIk97jx49kJqayi97\n8eIFHj9+DDMzM7Ru3Rpbt25FfHy8zGO7uLhAU1MTd+/elXvOQ4YMwfr165GUlITLly9j/fr1/DUa\nMmQIjh07hoMHD+LBgweIjo5GcHBwrfdhNXVdv/j4eLzxxhsy921OyJAoQLt27fDbb7/hl19+QXR0\ndK3b1NYnvq73MTEx6Nu3LyIiImBsbIx3332X75lSn/71MTExuHLlCqysrPD111/j//7v/6T2kUgk\nmDJlCsLCwjB16lQsWrQIb775JgAurBATE4MRI0Zg9OjReP/99+vUWlOPoaEhjh07hq1bt8LPzw8+\nPj6wsrJqkC5Zn1v37t2xa9cuHDx4EKamphg6dCjfy2X9+vWwsbHB8OHD0aFDB0yZMgWFhYW1fjYS\niQSTJ0/GgQMH4OXlhY4dO2L16tW1apg7dy5CQ0MxdOhQfPnll1i+fDkCAgIAAK6urggNDYWbmxtM\nTU1rPdby5cvh4eGBgIAAbN++Hfv374etre0r5yZLZ12fuyJty7puL6/38vLCtGnTEBQUhD59+qBr\n167o1KmTzH2r/z9x4gQePXqEYcOG8eHdbt26AQAGDhyIhQsXYs2aNejQoQP8/Pz4nory7puXiY6O\nRq9evdC7d2+sXLkSM2bM4NfJ0x4QEIAuXbrAzs4Or732Glq1aoVVq1Zh4cKFcHBwQGJiIkaNGiXz\n2Nra2pg1axb/gFXXOYeHh2PWrFmYNm0aIiIiEBkZibCwMACAnZ0d9u/fjx07diAgIADu7u5Yvnx5\nrdeqtmUvX4P169dj3rx5MnU3JxLWkEdnghAhEyZMgJWVFb744guhpRAipbi4GN27d8e1a9dUYlBi\nQkICPv74Y5Xp/qt0j+TevXvo27cv3NzcEBgYyCeiioqKMGTIEFhbWyM0NFQqPLB69Wo4OjrC1dVV\nKiaelpYGLy8v2Nvb49NPP+WXl5eXY9KkSbCxsUFgYOArfcoJoib0rEQ0Fn19faSnp6uEEQG4ke2q\nYkSAJjAk2traWLFiBVJSUvCf//wHn332GYqKirBu3TpYW1vj5s2bsLKywvr16wEAjx49wtq1a3Hi\nxAmsW7dOymWNiorCvHnzcOnSJZw+fZqPgR86dAhPnz5FWloaBgwYgEWLFin7NAg1Ql1LthCEqlB3\nVyAFMDc353sRtG/fHm5ubrh06RISEhLw2WefoVWrVpg4cSKWLFkCgEsYDRgwANbW1rC2tgZjDMXF\nxdDX18eNGzcwcuRIAMDQoUMRHx+P1157DfHx8QgPD4eenh6mTJmC/v37K/s0CDVi69atQksgCLWm\nSZPt6enpSElJgY+PDy5dugRnZ2cAXHe66qRbfHy8VDdMJycnxMfHIz09XSqh6erqyhe6S0hIgKur\nKwDA2NgYOTk59ephQhAEQSgfpXsk1RQVFWHkyJFYsWIF9PX1GxSnri0MwRjjl9fsIlr9vr7tEARB\nEHXT0Lxik3gk5eXlGDZsGMaOHYshQ4YAALy9vfmRxGlpafD29gbADWhLTU3l971+/Tq8vb3h4OAg\nNUguNTWVH7lbc5/8/HyYmZmhVatWtWqpNjoNel26BOblpdi+DXxFR0c3eJ/+O/rjl5u/NIu+xugU\n4lWXTlNThqws4TWqy+epSi/SqbyXIijdkDDGMGnSJHTt2hWRkZH8cl9fX2zZsgUlJSXYsmULP7GP\nj48Pjh07hszMTMTFxUFDQwMGBgYAuBDYnj17kJubi0OHDkkZkp07d+LZs2fYuHGj8icJ0tAAFPxA\nmwM9bT08L38utAzR0a4dIGOoCUEQjUDphuT8+fPYuXMnTp48CU9PT3h6euLo0aOIiIhAZmYmnJyc\n+PpLAFc/JiIiAkFBQZg2bRpWrVrFt/XNN99g2bJl8Pb2hr+/P1577TUAwDvvvIN27drBxcUFR48e\nxWeffabck5BIgHqW8RACXW1dMiQK0LYtIKMKB0EQjUDpOZLevXvLrKUkqxbSzJkzMXPmzFeWu7q6\nIjk5+ZXl2tra2LJlS+OE1oVE0mweSWBgYIP30dPWQ0l5ifLF1IEiOoWgLp3t2qmOIVGHz1OVIJ3C\nQiVSaqMZQ1uKGpLm9kjE8gWQZ0hUJbSlDp+nKkE6haXJem2JGhUPbVGORDEotKUcjI2NpWb2I8SJ\nkZER8vPzldIWGZLaUPFku66WLp5XkCFpKKoU2hIzBQUFCvfuIVQHZQ6PoNBWbYjAI2nuHIk6oEqh\nLYJQJ8iQ1EYzJtsVgUJbikGhLYJoGsiQ1IaKh7bIkCgGhbYIomkgQ1IbKh7a0tWicSSK0LYthbZa\nMtu2bYO/vz//XkNDAxkZGfXaNyYmBmPHjm0qaaKHDEltiMAjKamgHElDIY+EUJSGJKbHjx+Pzz//\nvAnVqB5kSGpDxT0SCm0pBiXbCaJpIENSG5RsV0so2a7+fPXVV3BwcICJiQnCwsJw9uxZhdp5/Pgx\n5s6dC3NzcwwbNkxqRlcAePfdd2FhYYFOnTrho48+4kNkGzduxA8//IBly5bBwMCAL1qrLF2qChmS\n2lDx0BbV2lIMCm2pPw4ODjh37hzu37+P1157DWPGjFGonYiICGRlZeHKlSsICQnBunXrpMJbgwcP\nRnp6OhISEpCbm4sFCxYAAKZMmYKwsDDMmzcPRUVFfFkoZelSVciQ1IYIQls0jqThUGireZBIlPNS\nhOHDh8Pc3By6urqIjIyERCJBUlJSg9qoqKjA8ePHERsbC3Nzc4wbNw5eXl5S24wfPx5t2rSBhYUF\nFixYgJ9//lmqxuDLAzaVoUuVIUNSGxTaUkvatAHKyrgX0XQwppyXIhw5cgRDhw6FpaUljI2NkZWV\nhatXrzaojbS0NFRVVcHe3p5f5vX3/ETVfPPNN+jXrx+MjY3h7e2NJ0+e4O7du/z6l5PzytClypAh\nqQ0VD22RIVEMDQ3AzAzIzhZaCdEUPHv2DJMnT8a4ceNw/fp15Ofno2PHjg0u5+Ls7AwNDQ3cunWL\nX5aUlMQbh/j4eCxfvhwrVqxAVlYWLl26BOAfL0RTU1PKO1GWLlWGDEltqHhoi8aRKI6FBZCVJbQK\noikoKipCcXExLCwsUFVVhSVLluDhw4cNbkdbWxv9+vVDbGwssrOzsXPnTly5coVf/+DBA7Rp0wam\npqbIysri8yPV9OjRA1evXkVFRYVSdakyZEhqQ8U9ktZarfGi8gWqmOoaO1XF0pIMibpibm6OJUuW\nYOzYsfDw8MCLFy/Qu3dvfr1EIpEKOdU1NmTt2rUwNTVF9+7dcejQIURERPDrQkNDERQUhO7duyM4\nOBgjR46UaiskJAQaGhro2LEjhg4dKleXOiBh6uRfvYREIlHMfczOBjw8gBpzxqsabRa3waPZj9BG\np43QUkRFRATQtSswfbrQSsSLwt8rQqWQdR0Vub7kkdSGiifbAcqTKAp5JAShfMiQ1IaKh7YAypMo\nioUFoGbhaYIQHDIktaHiyXaA6m0pCnkkBKF8yJDUhgg8EgptKQZ5JAShfMiQ1IZIPBIyJA2HPBKC\nUD5kSGpDBMl2qrelGB06AE+eAC9eCK2EINQHMiS1IZLQFtXbajgaGoCpqUr37CYI0UGGpDYotKXW\nUJ6EIJSLXEOSl5eHpUuXIiQkBACQmpqKzZs3N7kwQRGBR0LdfxWH8iQth+rZCs+dOwdnZ+d67XP2\n7Nl6b0twyDUkCxYsgIGBAe7cuQMAcHR0xIoVK5pal7CQR6LWkEfScqgui9K7d29cv369Xvv4+/vX\ne1uCQ64hSUxMxLRp06CpqQkA0NLS4v9XW0SQbKdxJIpDhRtbFqpSzqW6iKM6IteQeHl54d69e/z7\ngwcPwt/fv0lFCY4IQlvkkSiOpSV5JOrK3bt3MWXKFJibm2Py5Mn8j3dcXBw6depU63bvv/8+wsPD\n8fnnn9e6ra2tLdavXw8/Pz9YW1sjJiYG5eXl/Pq4uDgMGzYMjo6OWLZsGQoKCvh1Ghoa2L59Ozw9\nPeHs7IzTp0/DysoK69evh52dHVxdXXHy5EmcOXMG3t7ecHFxwQ8//NDUH5PSkWtIIiMjMX36dNy9\nexcODg5Yv349Zs6c2RzahEMEoS3KkSgOeSTqy7Bhw9C2bVtcu3YNTk5O2Ldv3ytVf6u3MzQ0xLVr\n1+Dm5ob9+/fLrAYskUiwfv16rF69GidOnMD333+PM2fOAABu376Nd955B6NHj0ZcXBySk5Mxa9Ys\nqf3//e9/Y9u2bUhJSQFjDI8ePUJqaioSExMxZswYhIWFYdWqVdizZw/Wrl2LKVOmSBkqMaAlbwMn\nJyccOXIEjx49QmVlJSwsLJpDl7CIJLR1r/Ce/A2JVyCPpGmRxCo4T+5LsOiGfQdzcnKQkpKCc+fO\noXXr1pg9ezZWr14td7tZs2Zh5cqVdbb93nvvwdvbGwDQv39//Pbbb3jjjTfw3//+FwMHDsTw4cMB\nAIsWLULPnj1RVVUFDQ3uOX3y5Mnw8PDg26qqqkJ0dDRMTEwwYcIELFiwAOPHj0fnzp3RuXNnWFpa\n4vz58wgMDGzQ+QuJXEOyZs0ahIWFwdTUFABQUFCA3bt3Y9q0aU0uTjBEEtqiHIlikEfStDTUACiL\nhIQEODg4oHXr1vyyl+dar97O0dFR7nY16d69O/+/ubk5P3vihQsX0LNnT36dg4MDKioqkJKSgm7d\nugEAfH19pdqysLCAiYkJAMDMzAwApAyNmZmZ6Ca+khva2rRpE4yMjPj3RkZG2LhxY5OKEhwRhLYo\nR6I4pqZAfj4gsugBIQdvb2+kp6ejpOSfB6zk5ORat7t58yZKS0vr3K4+9OrVC4mJifz7mzdvQlNT\nE25ubvwyLS25z+uiR64hadu2rVTyKD8/H7q6uk0qSnCqY6Uq7JVQiRTF0dTkSqXQ6Hb1wtzcHG5u\nboiOjsbjx4+xfPly5Px9kWv23KreLiYmBrm5uVi1ahWys7MVOmZISAiOHTuGgwcP4sGDB4iOjkZw\ncDAf1mopyD3b8PBwjBw5Env27MHu3bsxatQojB8/vhmkCYgIDAl5JI2D8iTqyf79+5Gfn4+uXbvi\n+vXrGDlyJIBXp9ndv38/cnNz4ebmhj///BODBw9Gu3bt+PV1TcNbsy17e3vs378fO3bsQEBAANzd\n3bF8+fI623l5WV3HEgtyp9pljOH06dPYv38/AGD48OEICAgQhcVt1JSgGhpc7ENFx8ycyzyHj49/\njHMTzwktRZQEBwPvvw8MGSK0EvGhblPtMsZgYWGBn3/+WW6uRJ1Q5lS7coN3EokEgYGBoupBoBRU\nPOFOHknjoDIpLZszZ86gS5cu0NbWxpo1a1BVVdWijIiyketWXL58GRMnToSLiwvs7OxgZ2cHe3v7\n5tAmLCqecKdxJI2DyqS0bG7cuIHu3bujS5cuePjwIY4dOya0JFEjN7QVEBCAKVOmoG/fvtDR0eGX\nt2/fvsnFNZZGueCtWgGFhdxfFeTuk7vw3+qPzFmZQksRJRs3ApcuAZs2Ca1EfKhbaKul0qyhrWfP\nnmHUqFHqX1/rZVTcI6FxJI2DPBKCUB5yDUlwcDAmT56MsLAwqfEkah9PVPHR7ZQjaRyUIyEI5SHX\nkJw6dQoSiQSLFi16Zblao+LJdl1tXZSUl4AxphbdB5sb8kgIQnnINSRxcXHNIEMFUfHQloZEAzqa\nOiitKIWutpoPEG0CTE2BvDygogJoAQOPCaJJqddX6MaNG/j111+lRrgvWLCgyUSpBCrukQD/5EnI\nkDQcLS2gfXtudHvHjkKrIQhxI7f77+LFizFnzhwsXboUBQUF2LZtGx49etQc2oRFxT0SgPIkjYWK\nN7ZcGjKdLk29Kx+5huTQoUM4dOgQ2rVrhxUrVuDs2bO4cuVKc2gTFhVPtgNUb6uxUJkU9UJfXx8G\nBgYwMDCAhoYG9PT0+Pe7d++W2rYh0+m+vK2trS1OnjypVO1ip14j2zU1NeHs7Ixr167B1tYW+fn5\nzaFNWEQS2iJDojjkkagXxcXF/P92dnbYvHkzgoKCXtmuoqKiURV5VWUcTWPPQ5nI9UjefvttFBQU\n4IMPPsDw4cPh5OSEiIgImdtPnDgRZmZmfC1+AIiJiYGVlRU8PT3h6emJX375hV+3evVqODo6wtXV\nFefO/VM3Ki0tDV5eXrC3t8enn37KLy8vL8ekSZNgY2ODwMBAhat2ykUkoa2SchpLoijkkbQM4uLi\n+Oltu3TpgkmTJjVoOt2a244dOxaZmZkIDg6GgYEBvv76a7z99ttYs2aN1DHd3d1x+PDhV7TcuXMH\nGhoaOHjwIFxcXODu7o6dO3fy6ysqKrBz50707NkTfn5+2LVrl9R0wTXPY+LEiYiNjcXo0aMxbdo0\nmJqaok+fPsjKysJ3330HBwcHBAUFISEhQemf6SuwBlBVVcVKSkrq3ObMmTMsOTmZde3alV8WExPD\n/vWvf70A4LO0AAAgAElEQVSybU5ODnNycmJ3795lcXFxzNPTk183cOBAtmfPHpabm8t69erFLl26\nxBhjbO/evWzYsGHs2bNnbMmSJWz69OkytTTw9KQxMWHs0SPF928Ggr4PYsdvHRdahmhZt46xyZOF\nViE+GvW9aiZsbW3ZiRMnGGOMnTp1imlpabGJEyeyrKwsVlJSwk6dOsWsrKyktvfw8GAJCQnsr7/+\nYra2tuz48eP8/i9vW902Y4zt27eP+fr68u+vXLnCTExMWHl5+Su6bt++zSQSCRs5ciTLzMxkx44d\nY61ateJ/V7ds2cLc3d3ZpUuXWFJSEuvevTvbunWrzPOIjo5mOjo6bOvWrSw3N5eFhIQwFxcXFhkZ\nyXJzc9kXX3zBgoKCav2MZF1HRa6vTL/owIEDvAtX2ziFoUOH1rqfv78/7ty5U5vBemVZfHw8BgwY\nAGtra1hbW4MxhuLiYujr6+PGjRt8CeihQ4ciPj4er732GuLj4xEeHg49PT1MmTIF/fv3r4+9bDgi\nCG1Rva3GYWkJ/Pij0CrUEGWNa1Li96+yshIxMTEwNzeXuY2s6XTlERwcjKlTp+LWrVvo3LkzduzY\ngVGjRtUZdpo7dy46deqETp06wdbWFqdPn0b//v3x3//+F9OnT8drr70GAJg+fToOHTrET91R23k4\nOzvz60eMGIHx48fjwoULMDQ0xIQJE/Dll1+irKwMrZqw3JPMM/3f//4HiUSCJ0+e4OjRo/D19YVE\nIsHvv/+OgQMHyjQksvj222+xf/9+vPPOO5g2bRoMDAyQkJAAFxcXfhsnJyfEx8fDxsaGn9oXAFxd\nXbFr1y5Mnz4dCQkJmDp1KgDA2NgYOTk5dX5IMTEx/P8NqmIsktAWGRLFoRxJE6GCD2BmZmZSoaza\nqDmdroWFBdLT0+vVduvWrTFixAjs2LED0dHR2LNnDw4cONCgYz148AAAN3VvzaEVPXr0wNy5c+s8\nD3d3d6n1VlZWMDQ05N+XlZUhLy8PlpaWtWqJi4tr9HhBmYZk27ZtADjLnJSUBFdXVwBc7iIyMrJB\nB4mIiMCCBQtQWFiIOXPmYMOGDZg9e7bMgmEvU9MrYoxJ7VdbGzWpaUgahAg8Eqq31TgoR9JyaGhS\nuq7fFU1NzVfWjxs3Du+99x569eoFPT29V+Zpry/VU/f26NEDAJCYmIg+ffrw618+D2VUtXj5ATs2\nNrbBbchNtmdlZcHKyop/37FjR2Q18DHO1NQUEokE7dq14101APD19UVqaiq/3fXr1+Ht7Q0HBwd+\nikwASE1N5S9MzX3y8/NhZmbWNC4beSRqj5kZkJvLjW4niPrSo0cPJCUlSS3z8/ODRCLB7Nmz8d57\n7ync9pAhQ7B+/XokJSXh8uXLWL9+PUJDQ2VuL+9BurmQa0gmT56MAQMGYPny5fjXv/6FQYMGYcqU\nKQ06SLXhqaiowA8//IBBgwYBAHx8fHDs2DFkZmYiLi4OGhoaMDAwAMDF/fbs2YPc3FwcOnRIypDs\n3LkTz549w8aNG9GzZ88Gaak3NI5E7dHSAkxMgJYwvralU58pb19eV3N9zf8/+OAD/PjjjzA2Npaa\nVve9997Dn3/+ifDw8AZrqSY8PByzZs3CtGnTEBERgcjISISFhcnc92WdsrZpcuqTkU9KSmKLFi1i\nX375JUtOTq5z21GjRjELCwumra3NrKys2ObNm9nYsWNZt27dWI8ePdisWbNYXl4ev/3KlStZ586d\nmYuLCztz5gy/PCUlhXl6ejJbW1v28ccf88tfvHjBJkyYwDp16sQCAgJYVlaWTC31PL3asbJiLDNT\n8f2bgc9OfsZi42KFliFqundnLDFRaBXiolHfKzXm+++/Z/7+/kLLqDeyrqMi11fuxFZiplEDh6yt\ngbNnARsb5YpSIkvOLkHhi0IseWOJ0FJEy6BBQEQEN4c7UT9UZUCeKvH8+XP06tULCxcuRLBIbiZl\nTmwlN7R1/PhxBAUFwdDQkC830LZt2wYdRJSIINlOoa3GQ/OSEI3l2LFjsLS0hJeXl2iMiLKR25Xh\n448/xqpVq+Dn5wcNDbl2R32gZHuLgLoAE42lf//+ePLkidAyBEWuZdDR0UGPHj1alhEBRJFsJ0PS\neKgLMEE0Hrkeib+/P0JDQ/Huu+/yg1wkEkmDBySKDhGEtqjWVuOxsABqlH4jCEIB5BqSnJwcmJub\nSxVUBGSXSFEbRBDaohIpjYc8EoJoPHINSfUI9xaHSDwSMiSNg3IkDcfIyKh5xiYQTYqRkZHS2pJr\nSCZMmCD1vvoG2rJli9JEqCQi8EjIkDQec3NuQGJlJaCpKbQacdAi5iMiGoRcQzJ48GDeeOTl5WHv\n3r18HRi1RiTJdqq11Ti0tQEjI+DxY86oEATRcOQakuHDh0u9DwsLq3XWMbVDBKEtGkeiHKrzJGRI\nCEIxGtyn9+rVq6hS8ZCPUqDQVouB8iQE0TjkeiT6+vp8aEtTUxNeXl5YsqQFlOQQSWiLDEnjoZ5b\nBNE45BqS4uLi5tCheoghtKWli5LyEpmzWBL1gzwSgmgc9ZrtJSkpCUePHoVEIsGAAQPg5eXV1LqE\nRwShLW1NbUgkEpRXlUNHU0doOaLF0hK4elVoFQQhXuTmSFatWoUPP/wQrVq1go6ODmbMmIFVq1Y1\nhzZhEYFHAlB4SxmQR0IQjUOuR7Jp0yZcuHCBr/g7ZcoUvP7665g5c2aTixMUEXgkwD+GxLC1odBS\nRIuFBeVICKIxyPVIHBwccPPmTf79rVu34ODg0KSiVAIRJNsBqrelDKiUPEE0DpkeSXVd/efPn8PH\nxwdubm4AgJSUFKmJ4tUWkYS2qN5W4zE3B3JyaHQ7QSiKTEMSFRXF///5559LrWsRPYREFtoiFEdH\nBzA0BHJzATMzodUQhPiQaUhahNdRFyLxSMiQKIfqPAkZEoJoOHJzJH/88QfGjBmD9u3bQ0tLCxoa\nGi1jql0ReSRUb6vxUJ6EIBRHbq+tL774AnPmzEFaWhpu3LiBTZs24cWLF82hTVhEkmynelvKgXpu\nEYTiyPVIMjIy4OvrC01NTbRp0wbz5s3Dvn37mkObsFBoq0VBHglBKE69am2VlZXhjTfewPTp02Fj\nYwNLS8vm0CYsIgptkSFpPBYWQGqq0CoIQpzI9Ui2b9+OqqoqxMbGwt/fH5qamuo/qRUgmtAWjSNR\nDlS4kSAUR65HYmtry/8/fvz4JpSiYmhoiMIjoXEkyoHKpBCE4jR4PpIWg4g8kucVZEgaC3kkBKE4\nZEhkQcn2FkX16HYROKEEoXLINST79++v1zK1Q0TJdsqRNJ5WrQADA250O0EQDUOuIVm8eHG9lqkd\nIgltUY5EeVAXYIJQDJnJ9l9++QU///wzHjx4gBkzZoD9/aP6+PHjltH9VyShLYNWBnha9lRoGWpB\n9aBEDw+hlRCEuJBpSCwtLdGjRw8cPnwYPXr04A2Jra0t/Pz8mk2gYIgktGWub47s4myhZagF5JEQ\nhGLINCQeHh7w8PDAmDFjoKPTAqdxFYlHYqFvgawi+vVTBlQmhSAUQ+44kocPH2LDhg04duwYCgoK\nAHBl5DMyMppcnKCIxCOxMLBAdnE2GGMto7x/E2JpCVy/LrQKghAfcpPt0dHR8PT0REVFBQ4dOoRB\ngwZhypQpzaFNWESSbG+t1RptdNogryRPaCmihzwSglAMuYbk6tWrGDFiBCQSCdzc3LBy5Urs3r27\nObQJi0hCWwCFt5QF5UgIQjHkhrZ0dXVRWVmJgIAALF68GHZ2dtDX128ObcIiktAWAFgaWCKrOAvd\nzLoJLUXUkEdCEIoh1yNZuXIlnj9/js8++wyMMZw9exbr1q1rDm3CIiaPxIA8EmVgYQFkZ4vm+YEg\nVAa5HomPjw8AwMDAADExMU2tR3UQkUdioW+BrGIyJI2ldWtAXx/IywM6dBBaDUGIB7mG5Pbt29ix\nYwcuXryI0tJSAFyvrZMnTza5OEERSbId4AzJnad3hJahFlTnSciQEET9kWtIPvzwQ/j5+WHBggXQ\n1tYGgJbRzVRkoa2L9y8KLUMtqM6TuLsLrYQgxINcQ5KZmYkff/yxObSoFhTaapF07Ajcvy+0CoIQ\nFzKT7fn5+cjLy8OYMWMQHR2NW7duIT8/n3+pPWIKbVGyXWl07gzcuiW0CoIQFzI9Ei8vL6kQ1vbt\n26XW3759u+lUqQIimSER+McjodHtjcfREThwQGgVBCEuZBqSO3fuAABKS0vRunVrqXXVSXe1RkQe\niUErA0ggQdGLIrRt1VZoOaLGwQG4eVNoFQQhLuSOI3n99dfrtUztEFGyHaDwlrJwcADS00V16QlC\ncGR6JFlZWXj48CGeP3+O5ORkPmzy6NEjtGrVqjk1CoOIku3AP+Etp/ZOQksRNYaGgK4uN+2uubnQ\naghCHMj0SH799VfMnj0bDx48QFRUFGbPno2oqCjs2bMHX3zxhcwGJ06cCDMzM3Tr9k+5jqKiIgwZ\nMgTW1tYIDQ1FcXExv2716tVwdHSEq6srzp07xy9PS0uDl5cX7O3t8emnn/LLy8vLMWnSJNjY2CAw\nMBDZ2U00F4eIQlsA55E8LKL6Hsqg2ishCKJ+yDQk48aNw6lTp7B161acOnWKf23btg39+vWT2eCE\nCRNw9OhRqWXr1q2DtbU1bt68CSsrK6xfvx4A8OjRI6xduxYnTpzAunXrMGPGDH6fqKgozJs3D5cu\nXcLp06eRmJgIADh06BCePn2KtLQ0DBgwAIsWLWrUByATsYW2qHCj0qA8CUE0DLk5kuHDhzeoQX9/\nfxgZGUktS0hIwKRJk9CqVStMnDgR8fHxAID4+HgMGDAA1tbWCAgIAGOM91Zu3LiBkSNHwsTEBEOH\nDpXaJzw8HHp6epgyZQq/XOmINLRFNB7ySAiiYcg1JMrg0qVLcHZ2BgA4OzsjISEBAGcUXFxc+O2c\nnJwQHx+P9PR0mJqa8stdXV3x+++/A+CMkqurKwDA2NgYOTk5KCsrU75osXkkBmRIlIWjIxkSgmgI\ncke2KwPWgB/k2sZB1BwfwRiTak9e2zULTQYGBiIwMLC+QsTnkVBoSymQR0K0JOLi4hAXF9eoNupl\nSO7du4cLFy5IFW1877336n0Qb29vpKWlwdPTE2lpafD29gYA+Pr64vjx4/x2169fh7e3NwwMDJCT\nk8MvT01Nha+vL79PamoqnJyckJ+fDzMzszp7kSlcsViEyXbySJRDdY6EMe42IAh15uUH7NjY2Aa3\nITe09emnn2LQoEE4efIkEhMTkZiYiEuXLjXoIL6+vtiyZQtKSkqwZcsW9OzZEwBXov7YsWPIzMxE\nXFwcNDQ0YGBgAIALge3Zswe5ubk4dOiQlCHZuXMnnj17ho0bN/JtKR2xhbbII1EaxsaAlhbw+LHQ\nSghCJDA5uLi4sNLSUnmb8YwaNYpZWFgwHR0dZmVlxbZs2cIKCwtZSEgI69SpExsyZAgrKirit1+5\nciXr3Lkzc3FxYWfOnOGXp6SkME9PT2Zra8s+/vhjfvmLFy/YhAkTWKdOnVhAQADLysqSqaUepyeb\nDz5g7LvvFN+/mamqqmI6X+iw5y+eCy1FLfDxYez8eaFVEETzo8jvpuTvHWUyatQoxMbGwslJfAPd\nJBJJg/IzUkyfDri6cn9Fgs1KG5wadwr2RvZCSxE9YWFA//5AAyK4BKEWKPK7KTNHEhwcDAAoKSlB\nt27d4OPjw3frlUgkOHLkSCOkigCRJduBv+duL8oiQ6IEaCwJQdQfmYYkKioKQO3WqUVUmBVZsh2g\nsSTKxMEB+PlnoVUQhDiQaUiqs/hz587FsmXLpNbNnTsXAQEBTSpMcESWbAeocKMyobEkBFF/5Pba\n+u23315ZVrPLrtoiwtAWeSTKo2YXYIIg6kamIVm3bh26deuGGzduoFu3bvzLxsamzlpbagOFtlo0\nJibc37w8YXUQhBiQGdoaM2YMBg4ciI8//hhLly7l8yRmZmbQ1dVtNoGCIaIZEquh0JbykEj+GeHe\nvr3QaghCtZHpkUgkEtja2uK7776DgYEB2rZti7Zt26KkpITmbFdRyCNRLpQnIYj6IdMjGT16NH76\n6Sf06NHjlV5aEokEGRkZTS5OUCjZ3uKhmlsEUT9kGpKffvoJwD9zt7c4RJhs76DXAQWlBSivLIe2\nprbQckSPgwPw669CqyAI1Udur63w8HBs2rQJ169fbw49qoMIQ1uaGprooNcBOc9y5G9MyIU8EoKo\nH3INycSJE/Hw4UN8+OGHsLOzw7Bhw7By5crm0CYsIgxtARTeUiaUIyGI+iG3jHxQUBD69OmDxMRE\nnDx5EuvXr8e1a9cQGRnZHPqEQ4ShLYAS7sqkQwegvBzIz+cqAhMEUTtyDckbb7yBZ8+ewc/PD717\n90ZiYqLU7IVqC3kkLZ6aXYB9fIRWQxCqi9zQlru7O7S1tXHt2jVcvXoV165dQ0lJSXNoExbySAhQ\nnoQg6oNcj2TFihUAgKKiImzbtg0TJkxAdnZ208yTrkqIMNkOcIbkSs4VoWWoDZQnIQj5yDUk3377\nLc6ePYukpCTY2dlh4sSJ8Pf3bw5twiLCke0AF9r6Jf0XoWWoDQ4OwMmTQqsgCNVGriEpLS1FVFQU\nvLy8oK3dgsYmUGiLAGdINm4UWgVBqDZyDcmcOXOaQ4fqQcl2AhTaIoj6IDfZ3mIRqUdirm+OR88e\noYqJT7sqYmYGlJQAT54IrYQgVBcyJLIQabJdR1MHbVu1Re7zXKGlqAU1uwATBFE7cg1JcXExKisr\nAQA5OTm4cOFCk4tSCUQa2gK4udsfFj0UWobaQIaEIOpGriHp06cPysrKUFhYCF9fXyxcuFD9R7UD\nog1tAZQnUTaUJyGIupFrSKqqqqCnp4ft27dj4sSJOHr0aMvwSkQa2gKo55ayIY+EIOpGriExMTHB\niRMn8P333yM8PBwAWsbIdpGOIwHII1E21fO3EwRRO3INyb/+9S/s2LED77//Puzt7XHr1i307du3\nObQJC3kkxN+QR0IQdVPnOJLKykrs2rUL27Zt45d17twZq1evbmpdwiPiZLuFvgVO3z0ttAy1wdIS\nKC4GCguBtm2FVkMQqkedHommpibOnDmDoqKi5tKjOlCynfgbiQTo3Bm4dUtoJQShmsgd2d6rVy8E\nBwdj+PDhsLS0BGMMEokEQ4cObQ59wkGhLaIG1XkST0+hlRCE6iHXkOTn58PW1hZJSUlISkril6u9\nIRFzaMvAAtnF2bzRJxoP5UkIQjZyDUnN/EiLQsShLT1tPeho6uBJ6RMY6RoJLUctcHQELl4UWgVB\nqCZye23dvn0bERER8Pzbp7969SoWLVrU5MIER8QeCUDhLWVDHglByEauIYmJiUFwcDD/vlu3bti9\ne3eTilIJROyRAJRwVzY0loQgZCPXkPz1118YNGgQ/76qqgo6OjpNKkolEHGyHSCPRNl07MhVAC4u\nFloJQagecg1J7969+SR7WVkZvv32W/Tv37/JhQmO2ENb5JEoFQ0N6gJMELKQa0giIyOxdu1aZGdn\nw97eHikpKZgxY0ZzaBMWsYe2yCNROhTeIojakdtrq2PHjti8eTMqKipaTlgLUIvQVlJWkvwNiXpD\nCXeCqB25hqSoqAi//vorLl68iLKyMgCARCJR/zIpIi7aCFBoqylwcAASE4VWQRCqh1xDMmXKFOjq\n6sLPzw86OjotZ5CbGngkFNpSLo6OwJ49QqsgCNVDriFJSUnB1atXm0OLakHJduIlKEdCELUjN9ke\nERGBL7/8EhkZGcjPz+dfao/Ik+3tWrVDeVU5nr14JrQUtcHKCsjLA57RR0oQUsj1SPT09BAVFYW1\na9fyiXaJRIKMjIwmFycoIg9tSSQSWBpYIqs4Cw7GDkLLUQs0NQE7OyAjA+jWTWg1BKE6yDUksbGx\nSEtLg42NTXPoUR1EnmwH/s6TFJEhUSbV87eTISGIf5Ab2nJwcICurm5zaFEtRO6RAH/nSSjhrlQo\nT0IQryLXIzEyMoKHhwf69esHQ0NDAC2o+6/YDYk+JdyVjYMDcOWK0CoIQrWQa0gGDhyIgQMHAuAM\nSIvq/qsOoS3ySJSKszOwY4fQKghCtZBrSMaPH8//n5SUhB49ejSlHtVBTUJb1/OuCy1DrfDzA1JS\nuN5bJiZCqyEI1UBujqQmkydPbiodqgeFtoha0NUFAgKAY8eEVkIQqkODDAlr5A+rra0t3N3d4enp\nCR8fHwBcCZYhQ4bA2toaoaGhKK5Rp3v16tVwdHSEq6srzp07xy9PS0uDl5cX7O3t8emnnzZKk0zU\nIbRFyfYmYfBg4KefhFZBEKpDgwxJdHQ0AMUNikQiQVxcHC5fvoyEhAQAwLp162BtbY2bN2/CysoK\n69evBwA8evQIa9euxYkTJ7Bu3TqpisNRUVGYN28eLl26hNOnTyOxKQogkUdCyGDwYM4jqawUWglB\nqAZyDcnnn3/O/x8aGorKykqEhYUpfMCXjVBCQgImTZqEVq1aYeLEiYiPjwcAxMfHY8CAAbC2tkZA\nQAAYY7y3cuPGDYwcORImJiYYOnQov49SUQOPxETPBIVlhSirKBNailrRqRM30dXvvwuthCBUA7mG\n5N69e1iyZAkAbmKroUOHokuXLgodTCKRICgoCKGhoThy5AgA4NKlS3B2dgYAODs7855KfHw8XFxc\n+H2dnJwQHx+P9PR0mJqa8stdXV3xe1N8o9Ug2a4h0YCZvhmyi7OFlqJ2DB4M/Pij0CoIQjWQ22tr\ny5YtCAsLw5IlS3Dy5EkMGjQIs2bNUuhg58+fh4WFBdLS0hAcHAwfH58Ghclq63Ysb/+YmBj+/8DA\nQAQGBtbvYGoQ2gL+6QJsY9jCKhM0MYMHAxERwN/PWAQhWuLi4hAXF9eoNmQakqSkJP6HOzIyElOn\nTsXrr7+OgIAAJCcnw8vLq8EHs7CwAAC4uLggJCQE//vf/+Dt7Y20tDR4enoiLS0N3t7eAABfX18c\nP36c3/f69evw9vaGgYEBcnJy+OWpqano2bOnzGPWNCQNQg1CWwBVAW4qevYEHj4EMjMBa2uh1RCE\n4rz8gB0bG9vgNmQakqioKCkPwNDQEGlpaYiKigIAnDp1qkEHev78OSorK2FgYIDHjx/j2LFjmDVr\nFgoKCrBlyxYsW7YMW7Zs4Y2Cj48P5syZg8zMTGRkZEBDQwMGBgYAuBDYnj170K9fPxw6dAgrV65s\n8InLRQ1CWwANSmwqNDWBAQOAn38GPvhAaDUEISwyDUljXZ2XycnJwTvvvAMAMDExQVRUFDp16oSI\niAiEh4fDyckJXl5eWLp0KQDAzMwMERERCAoKgo6ODjZs2MC39c033yA8PByffPIJRo0ahddee02p\nWgGoRdFGgAxJUzJ4MLBrFxkSgpCwxg4OUWGqS7ooxNGjwIoVoh95tjFpIxIeJODfIf8WWorakZ8P\n2NgAjx5xAxUJQh1Q5HezQeNIWhRqlmwnlI+xMeDpCTQwyksQaodMQ7J//34AUP8JrGRByXaiHtAo\nd4Kow5B89dVXAIBhw4Y1mxiVgpLtRD2oHk+iBrcKQSiMzGS7s7MzAgMDcfv2bQQHB0utk0gk/IBC\ntUVNku1m+mbIe56HiqoKaGnIHTZENBA3N+5vSgrQtauwWghCKGT+suzatQtXr17F0KFDMXv2bKnk\nS4uZj0QNHjO1NLTg3N4ZSQ+T4GvlK7QctUMiAd5+mwtvkSEhWip1PqK6u7vj4sWL6NChAyoqKrgd\ntFrIU62aJNsBoJ99PxzPOE6GpIkYPBj46itg3jyhlRCEMMjttVVYWIiRI0fC3t4e9vb2GDVqVMtI\nwKtJsh0A3rB7AydunxBahtrSty83/W5+vtBKCEIY5LoXixcvRkhICHbt2gWA68315ZdfYvPmzU0u\nTlDUJLQFAH1s+mDkf0bieflz6GnrCS1H3Dx4AJw8CbRtCxgZAUZG0DUywpuvG+HYUT2MHlMj7FtS\nwm3/8uvhQ66EcHg40L07d68RhIiROyCxe/fuSE5OhoYG57xUVlaiR48euHLlSrMIbAyNGpB4/jww\ndy73Vw3w3+qPz/t8jrc6vyW0FHGSlwcsXQps3gwEBQGlpUBBAf+qyC0AKiuh1d6IMzJ5ecDz54Cl\nJVdzvubL0pLLzu/cCbRpA4wdC4SFAVZWQp8lQSj0uynXIwkODkZkZCTGjx8Pxhh27NjxSi8utUSN\nQlvAP+EtMiQNpKgIWLkSWLUKePdd4M8/OUPwEg8zAT/PUmQmFUDzWSE3obuJSd3eRmws96CyYwfg\n4cF5J+HhwLBhnDEiCJEgN0cyZ84cODg44JNPPsH8+fPRuXNnzJkzpzm0CYsaJduBfxLuRD0pLeUM\niKMjcP06N4vVunW1GhGAqwBs0rE14jMtACcnoH17+SErDQ3A3x/YuJELeU2bBhw+zDUWHg4kJzfB\niRGE8qFaW7KIjwc+/BD4e6ItsVNeWY72X7dHxowMmOiZCC1HdamoALZv57wFDw9g0SLA3b1eu37y\nCWc7Fi9upIa8PGDbNs6QOTtzIdZ+/SiXQjQLVGtLmahRsh0AtDW10du6N07docJQMsnN5TyE778H\ndu8GjhyptxEB/hlP0mhMTICoKODWLc4zmTUL8PLiNP3dDZ8gVAkyJLJQs9AWAPSzo/CWTDIzOSMS\nEADExQGvv97gJnr25CJU9+4pSZOODjBuHHD1KucZbdjAhdq+/RZ49kxJByGIxkOGRBZqlmwHKE8i\nk5QUoHdvYMoUbmShgiEkTU2gf39usiuloqHBjXqMiwP27OH+2tlxvcieP1fywQii4cjstSVrusXq\n8igLFixoGkWqgpqFtgCgq2lXFL0owp0nd2BraCu0HNXgwgXgnXeA5cu5LriNZPBg7rd+6lQlaKsN\nX1/gwAGuA8CCBZyH8vnnwKRJgLZ2Ex2UIOpGpkfSpk0b6OvrS70kEgk2b97Mz2Ko1mhqAi9eCK1C\nqehOyyYAAB8CSURBVEgkEq4bcAaNcgfAJTSGDOFyIkowIgA3/W5cXDM4Cs7OwL59XC+vgwcBFxcu\nh6JmXjQhDurVa6uwsBCrV6/G5s2bMWLECERFRcHU1LQ59DWKRvXaKi0FTE2B27e55KeasDl5M47f\nPo7dw3YLLUVYtm/nekMdPsw95SuRkBAuxDV9ulKbrZuTJ7luY2VlXLexgQOplxehEErvtZWXl4fP\nPvsMHh4eKC8vR3JyMpYuXSoKI9JoWrcGAgOBX38VWolS6WffDydvn0QVa8FPrt98w4WDTp1SuhEB\ngPnzgWXLgPJypTctm6AgbqxLTAwwZw7Qpw+QlNSMAoiWjExDMnv2bPj4+MDAwABXr15FbGwsjIyM\nmlOb8AwapHbT39kY2qBtq7a49uia0FKan6oq7kd2yxbg3DkuHNQE9OwJODhwFVCaFYkECA3lenlN\nnMjF2fbta2YRREtEZmhLQ0MDOjo60K4lgSeRSFBYWNjk4hpLo0JbANcl1MsLyMnhciZqwgc/foAu\nJl3wkd9HQktpPp49A957j7uWhw83ebjy5EkgIgJITRXw1rlyhYuzTZvG1binUBdRD5Qa2qqqqkJp\naSmKiopeeYnBiCgFa2vAwkJtRrdX08++X8sqK3//Phfq0dcHTpxolpxX376AsTHXwUowuncHLl4E\n9u4FJk9u5lgb0ZKgcSTyGDy4CQYGCEtf2744e/csXlSqV6+0Wrl0iYs1jRzJlR1p1apZDiuRcLmS\nxYsF7kXesSNw9iyQnc0l4J88EVAMoa6QIZGHGuZJTPRM4GjiiIQH6uVpvcK+fdz1W7OG66HVzKGd\nt9/mjIjgzyH6+lw4z80N6NULuHNHYEGEukGGRB6vv8598R4+FFqJUlHrUe6MAQsXcon1337jEtAC\nUO2VfPmlCoxt1dTkSuFPncrd02oWriWEhQyJPLS0gDffBH75RWglSkVt626VlABjxnBuQHw8lycQ\nkOHDuWK+cXGCyviHGTO4svVvvy1wAodQJ8iQ1Ac1zJP0su6FP3L+QFFZkdBSlMetW9zYH4mEGyNi\nbi60ImhqAh9/zHklKsPbbwPHjgEffcQZlpISoRURIocMSX0YMIDr7aNGJVP0tPXgbemNM3fPCC2l\ncTDGPe6HhnJJ9eHDgV27AF1doZXxhIUBN29yDlJjyS/Jx49//YjPTn6GpeeW4kDqAfyR/QeKXxQ3\nrCFPT6578KNHgLc38McfjRdHtFhoYqv64usLLFnCjSBWE7488yVyS3Kxov8KoaU0nNJSrrbUypVc\nt9aZM7m5O9q0EVpZraxZw6VrDh+u/z6MMdx9ehfnMs/xr8ynmfC18kVPq54oKS9Ben460vPTkVGQ\ngbat2sLB2IF/eZp7IsguCLradRhVxriRk1FRnOsUGclVGyZaLIr8bpIhqS+xsdz83d98o5z2VID4\n+/GY/L/JuBpxVWgp9Scri5vydsMGoEcP7ofvzTdVfrBdSQlgb89V3OnWre5tr2RfwdLzS3H27llU\nVFXA38YfvTv1Rm/r3vAw94CWxqtFu6tYFbKKsnjDkl6Qjgv3LuBK9hUE2QUhpEsIBncZDNM2Msob\n3b7NGWI9Pa6bdMeOjT9pQpSQIXkJpRqSxERg7FggLU057akAFVUV6PB1B6RNT4O5vvD5hDq5d48r\nm374MDB6NDcNsrOz0KoaxLJlXDTphx9qX19UVoTouGjsvLoTn/p/isFdBqOzUWd+6gZFyHueh59v\n/owjfx3Bb7d+Q1fTrghxCkGIUwicTJyk266o4LzuNWuAtWuBYcMUPi4hXsiQvIRSDUlVFWBpyc1f\nYW+vnDZVgNA9oRjhNgJjuo0RWkrtFBZyEzitX8/VHImKAkRa862wEOjcmRts7uDwz3LGGP57/b+Y\neXQmguyC8PWbX6NDmw5KP35ZRRlO3TmFIzeO4MiNI2ij0wbDXIbhXdd30d28+z9GJT6eS+wEBHDz\ntLRrp3QthOpCc7Y3JRoa3MhgNeu9pbLlUioquBBWly7c/LV//MFNNytSIwIAbdtyZa+++uqfZXef\n3EXInhDMPzkf29/Zjm2h25rEiABAK61WGOAwAGsHr8W9Wfewa+guVLJKDNs3DF3WdMEnJz5BclYy\nmI8PcPky1/W9SxdgxQquPD1ByIA8koawfz+wdataGZPrudfx1o63cDfybqNCKEqDMa6SwJw5nAf4\nzTdcDyM1IS+Pm9Qw8XI5/nN/BZadX4ZZPWdhTq850NHUEUQTYwzJWcnYn7of+1P3QwIJhrsOx7uu\n78IrVxuSTz8F/vwT+OILboyOGhUwJV6FQlsvoXRD8vQp0KkTV7dIT0957QoIYwxWK6wQNy4OjiaO\nwoq5fBmYPZtLqH/9NVfeRBWMm5KZ9Hki9pROgJ9bR2wI+Q6djTsLLYmHMYbL2ZexP3U/9qXsQ3u9\n9ogJiMGArDaQfPwxV0V5yRKaOEuNIUPyEko3JAA34G3OHG6Qopow6cgktG3VVrhuwJcucT9OFy9y\nEzNNmsSFVdSQvdf24v9++T94Zq3Go1OjcOyoBGZmQquqnSpWhQOpBxB7OhZtdNogpk80BqS+gGT+\nfG720K++4sbuEGqFQr+bTI1pktNbtoyxiAjltysgj589Zl2+7cJW/766+Q5aVcXYqVOMvfkmY506\nMbZ6NWPPnjXf8ZuZqqoqtjBuIbNeYc3+yP6DVVUxFhPDmIMDYxkZQqurm8qqSrbv2j7m9p0b897o\nzX5OPcKqNm1izMqKsdBQxv74Q2iJhBJR5HeTPJKGkpLCeSO3b6uVa3+74DZ6b+2Nlf1X4l23d5vu\nQIwBP/7I1VfPz+cGwYWFATrC5Aeag9KKUrx/5H38lfcXDo86DAsDC37dd99xztjRo0DXrgKKrAdV\nrAoH0w4i9nQsWmu1xkLfTzDgt9uQfP0118MrOhpwdRVaZrNQXlmO9Px0XM+9zr3yrqOorAidjTvD\n0diRHxRq1dYKGhJx9Wmi0NZLNIkhYQywteUS7m5uym1bYK5kX8FbO97Cvnf3IdA2ULmNV1RwnRWW\nLOGStfPnA0OHqn3i9vGzx3hn7zuwMLDA96HfQ0/71dzanj3cwPxDh7jCvKpOTYPCGMNouxBMOP8c\nFht/gOTNNzmD0qWL0DLrRRWrQlFZEYpeFOF5+XOUlJdwfytKUFJewv99Vv4Md5/eRdrjNFzPvY47\nT+6gU7tOcG7vDJf2LnBu7wx9HX1kFGQgPT8dN/NvIj0/Hfkl+bAztIODsQM6G3eGeRtzdGjTAR30\nOvB/TduYQl9HXzU6u4AMySs0iSEBuD6cdnZcrkTNOHn7JEb9ZxSOv3cc7mbujWvs4UOuOODRo8Dx\n49zT6vz5XO0yFfnSNCWpj1MRvDsYo7uOxsK+C+t8Mj12jBvvum0b18dADFSxKly4dwGHrh/CwbSD\naFsGfJXWEf2OpEArOASSzxdwA2cE4tGzR0h4kID4B/HIKMjA09KneFr2lP/7pPQJil/8f3vnHtTk\nlf7xb8Il4RoUBIRAYBHlJiSyXKxrReyitSK01V7XHYvt6K6t06vbna2z7Dh111pHGXesXSurtT91\nq1bbOlpaLxRQuVTwhoAiRhAohnsgEAg5vz+OeZNAKArFvOD5zDxzzvsmeXny8OZ83+e85z2nA452\njnAVucLRzhEOtg5wsHOAg60D3b5Xd7BzgEwi44RjysQpENkOvUhaZ08nqlqqcKP5BqpaqtDQ2QBV\npwoqjcqs1Ol1nLB4OHpw1n87wC0AgRMCRzVuTEj6MWpCcuwYHZbKm7nBf10OXD2Ad79/F2fTzkLm\nJrv/D2q1wNmzVDi++44ucfv73wPz51N7hKbd+OHmD3j5q5fxcdLH+GPUH+/rM/n5dO7JzZtpb99Y\nghCCSw2XcKT8CE5eOISns25jVb4OLfMfh81rq+CR+BTs76PhHS6aXg1K6ktQUFvAiUdrdytifGIQ\nJ43D1IlT4SZ2g0QsgUQk4UpXkStshNbPirt6u3C38y4aNY3m1tUIVaeK2y5TlcPDyR2pIalInZaK\n3/r89lfPZJiQ9GPUhESjAby8aEM5Tp/63Zq/FZ9e+BR5r+TB3XGQNc4JofeMTp2iGUdODhAaSjOO\n+fPprLLjdPTVYOj0OnxS9Ak+zP0QB5cexGzZ7Af6/LVrNHx//jPt7uLRJMYPRGVzJU4U/B8cP9uD\nOXk1sOnR4dgMJ+TPDoQ2dCr83WTwc/WDv8QfzvbO0BP9oNar7zVmEyYZhaFs7mqGslWJsElhiPON\nQ6xvLOJ84xDsHjzm7k/o9TSRv3HD3CorgaoqAAI9eiYVQhR1FLopRwH7Dkg7UxCKVEx3nYNJE+1B\nCJ2ovKeHXtsZ6gYjhN6SHMzWrmVCYsaoCQlAx9GnpQFLR/HGtJVZ+8Na5FXn4eQfTxr79qurjcJx\n+jR9nmbePGpPPAG4DyI645xuXTd2X9yNTec2wdfFF5kpmZgyccrQH7TA7dt0NpiiIuCVV4BVq8b4\nrDyEoK+kGJ17d0F08Ct0i+1wJTEcZ+K9UOzUDk2vBkKBcIDZCGwgFAhhK7Q1yyBMMwqJWAI3sRum\nuk+F2FZs7W96XxACqFTA9etUJEzLmzfpDAhTptAHV00tKIiumqzV0nEqjY1AcXU5vq/+GudbjuLn\n3gr49cyHl34GPIXT4GU7FZ52v4GT2N5MKAA6YXZ/gTHYxx8zITFjVIXk3/8GLlygT7o/bLq66KXr\nlSvA5cv0rJLJ6CAAg0mlgJ2d+ecIoWdfeTlQUWEsKyroMV1c6Jnq7Ay4uIA4O+FM409ot9NjsW8i\nhGeygZYWo3DMmzfGW7iR09bdhk9++gQZBRmI8YnBX2b9BbP8Z/0qx755k04xtns3EBtLs5QFC8b4\n+AS9nj4vdOAA8OWXgL8/HXSxcCEQGcnLe2c9PfT6Sa02NraWGmKtlq5u0NVluVSraVZx/bpx9png\nYGNpMBeX4flZp67DiRsnUKoqRUVTBa43XUdNWw38JH6Y5j4N0zymYerEqZgycQpk9zJCS/d5WNdW\nP0ZVSKqq6BCburrRWb+hr4822o2N9My7fNkoHEolPfumT6c/Pnd3eqYrlXRYslJJn7738aGi4u1N\nXy8vp2ISEkJt2jRj6egIdHQYTa0GOjqga29FZt422NqLIX/+TciTlkFo82h1V1miXl2PrQVb8Vnx\nZ3gq+CmsnbUWEZ6jM363qwv43//oUOGmJpqhpKUBHh6j8uceHjodXcny22/pUtYaDVXKJ5+k2a2b\n20Nxo7sbaGgw/nwMPyFD/e5d+lOSSH65S8jennZFisXG0rTu5ESvu4KDH17i3tPXg6qWKlQ0UmGp\naKpAVUsVlK1K1Kpr4eHoAZlEhgC3AMjcZAiQBGBVzComJKaMqpAA9KnemzfptClSKTXTulRKG+j2\ndmptbQPrbW20dehvbW00x3V3p2eeQTQiI2nDP9RzFz099B6OUkmnHPHzo6IxadIDX/V19HRga/5W\nHLx2EKpOFZ4JfQZLw5bid/6/48WNytFGT/Ro6GjA7bbbuN16G6duncKha4ewLGoZ3o5/+8EGJIyQ\noiIqKEePAlFRdL21+HhajvmxDDduUEE5cQLIy6NzrD35JB2wERFBW+P7hBCaqCuV1GprqSCoVLQ0\nta4uesszIIAOxjSUhrpUOj5v9fXp+1CnroOyVYnbbbe58rPFnzEhMWXUhUSvp2dmTQ1ttO/cGVjv\n6qKXMq6u1CzV3d0H2oQJvDx7rzddx8HSgzhUdgj16npOVB6XPT5mREWr00Ldo0a7th3t2naotcZ6\nu7Yd9R31qG6r5oTjTvsdSMQSyCQyyNxkkHvJsfK3K+HhaL2UoK0NKCykI70KCmjp4GAUlfh4QC6n\nvZR8hRCaCTQ2Gq+vDKZp6oL7lWz4l57Ab6rPYFJbJZonBOHu5Cg0+UahVRaF9sAo6D29YWdnzChM\nzdbW2NPr60vFwtPT3CZNoj9DHvaoWY1HqmsrJycHK1euhE6nw5o1a/DGG28MeM+oC8mvQHZ2NhIS\nEqztxpBY8vNG0w0cunYIh8oOobqtGoFugZjoMBETHCZgosNEWhcb6w62DtATPQgILQkx29bpdWjX\ntlscoWPYr+3TQqfXmVlvXy9X197UwibQBgQEhJABZR/pAwBIRBK4iFzgKnKFq8gVLva07iJygbez\nNxWNe8Lh5+r3y8vV/krxHAmE0OTYICr5+cDVq/R6xLQP3lAPCrq/C/zh+qnT0R7StjZjYnz7tnlZ\nXU2FzsvLeG3l4mKsG0wsBvo6u+FcfQ0Tay7Bo+4SvH6+BB/VJegEdrjtFoUcBxGiQ2IhmiKFJMwX\nk+RSuIZJ6QF4xFj4vT9SQqJQKJCRkQGZTIb58+cjLy8PHv06jceCkKSnpyM9Pd3abgzJUH7eab+D\nOnUdmruaOWvpakFzt3G7W9cNAQQQCoQQCO6VEHB1G4ENHZVjGJHTb3SOq8gVYlsx7IR2sBXampmd\nDd330Ycf4a8f/BUCgYA7tmkpFAhhb2Nv9aeIH8b/Xa+njbhhCKlhdNCNG7Qh9/CgPa/29oBIZLks\nK0tHaGg6CKHHIwRm9d7eAbfV0NFB998bswGplGYFMplxTIih7uQ0gi9ICP2Cly8jfft2pMvldLu2\n1tgrIBBQB3x8gMmTB5q3Ny1dXR9KWjIWfu/DaTf513dyH7S1tQEAHn/8cQBAUlISCgoK8NQ4mpF3\nrCF1lULqKrW2GxDbiiERj89nex4UoZAOivL3pwPsTNHp6DiRrq6BzxyYll9+CTz3HD2WQEDNtG5r\nS8Wi34A/iMUPoV0WCOi9Pz8/evOofwNNCO0nq6mhX7a+nlp1NU3dDNv19XRwi5ubscvZ0O1sqEsk\nxi/a3wxf2tGRBsRghkA9AoxJISkqKkKIyXrdYWFhyM/PZ0LCYNwntrZUYIbi6lVgyZLR92dUEAiM\nIjDUjJidnbQfzmCGgTCm9Tt3jKmXJdNoqCDpdLTU642iYmNDy54eYPt2um0ww+s2NtRnQ9oHGOv9\nMwSDQBkU3XTb9Nj9TSg0FzhToRuB6I3Jrq2TJ09i165d2L9/PwBgx44dqK2txfr1683eZ+3uCwaD\nwRiLPBJdWzExMXjPZMLE0tJSLFiwYMD7xqBGMhgMxphjbE1Ecw/JvfmtcnJyoFQq8cMPPyAuLs7K\nXjEYDMajyZjMSABg69atWLlyJXp7e7FmzZoBI7YYDAaD8XAYkxkJAMyZMwdlZWWorKzEmjVrzF7L\nyclBaGgogoODsW3bNit5ODQBAQGIjIyEQqFAbGystd0BAKSlpcHLywvTp0/n9qnVaqSkpMDf3x+p\nqano6OiwoocUS36mp6dDKpVCoVBAoVDgu+++s6KHlJqaGsydOxfh4eFISEjAvn37APAvpoP5ybeY\ndnd3Iy4uDnK5HPHx8diyZQsAfsVzMB/5FksDfX19UCgUSE5OBjDMWD7w4rxjALlcTn788UeiVCrJ\ntGnTiEqlsrZLFgkICCBNTU3WdsOMnJwcUlxcTCIiIrh9GzduJK+//jrp7u4mq1evJps2bbKihxRL\nfqanp5PNmzdb0auB1NfXk5KSEkIIISqVigQGBpL29nbexXQwP/kY087OTkIIId3d3SQ8PJxcv36d\nd/G05CMfY0kIIZs3byYvvfQSSU5OJoQM7/c+ZjOSwTB9xkQmk3HPmPAVwrMBAbNnz8aECRPM9hUW\nFmLFihUQiURIS0vjRTwt+QnwL57e3t6Qy+UAAA8PD4SHh6OoqIh3MR3MT4B/MXV0pEsadHR0QKfT\nQSQS8S6elnwE+BfLO3fu4Pjx43j11Vc534YTy3EnJIM9Y8JHBAIBEhMTkZqaim+++cba7gyKaUxD\nQkJQWFhoZY8GZ9u2bYiPj8fGjRuhVqut7Y4ZlZWVKC0tRWxsLK9javDTMICFbzHV6/WIioqCl5cX\nXn/9dfj7+/MunpZ8BPgXy7feegubNm2C0GQG8+HEctwJyVji7NmzuHTpEv75z3/i7bffxs8//2xt\nlyzCt6uowfjTn/6EW7duISsrCzdv3sSnn35qbZc41Go1nn/+eWzZsgXOzs68jampn05OTryMqVAo\nxKVLl1BZWYnt27ejpKSEd/G05CPfYnns2DF4enpCoVCYxW84sRx3QhITE4Py8nJuu7S0FPHx8Vb0\naHAmT54MAAgNDcXixYvx7bffWtkjy8TExKCsrAwAUFZWhpiYGCt7ZBlPT08IBAJIJBKsXr0aR44c\nsbZLAIDe3l48++yzWLZsGVJSUgDwM6aW/ORrTAE6WGXhwoUoKCjgZTwBcx/5Fstz587hm2++QWBg\nIF588UWcPn0ay5YtG1Ysx52QjJVnTDQaDZfaqlQqZGVlWXyokg/ExcUhMzMTXV1dyMzM5K0w19fX\nAwB0Oh327duHhQsXWtkjenW3YsUKRERE4M033+T28y2mg/nJt5g2NjaitbUVANDU1ITvv/8eKSkp\nvIrnYD7yLZYbNmxATU0Nbt26hQMHDiAxMRF79+4dXixHYxSAtcnOziYhISEkKCiIZGRkWNsdi1RV\nVZGoqCgSFRVFEhMTya5du6ztEiGEkBdeeIFMnjyZ2NvbE6lUSjIzM0l7eztZvHgx8fPzIykpKUSt\nVlvbTc5POzs7IpVKya5du8iyZcvI9OnTSXR0NHnrrbd4MSIuNzeXCAQCEhUVReRyOZHL5eTEiRO8\ni6klP48fP867mF6+fJkoFAoSGRlJkpKSyJ49ewghhFfxHMxHvsXSlOzsbG7U1nBiOSbn2mIwGAwG\nfxh3XVsMBoPBeLgwIWEwGAzGiGBCwmAwGIwRwYSEwWAwGCOCCQljXKBUKs0mcOQTy5cvx+HDh0fl\n2M7OzqNyXAbjQWBCwmAME71ef1/vG82VOtkqoAw+wISEMe6oqqrCjBkzcOHChQGvJSQk4IMPPoBc\nLodCoUBlZSWWLFmCiIgI7Nixg3vfyZMnsXTpUsycORMbNmzg9js7O2PdunWQy+UD5nDbuXMnYmNj\nER0djbVr16Knp4d77dy5c4iOjkZCQgI3Cd7u3bvxxhtvcO9ZtGgRcnJyuL+zfv16hIeH46WXXkJz\nczMAoK6uDitWrEBISIiZXx0dHXjiiScwY8YMLFy4ED/++CMAmqmFhYVh9erVCAsLw6pVq9Db2wuA\nTh3/zjvvQKFQIDo6Grdu3QIAHDx4EIsWLcLs2bPxn//8Zxj/AcajBhMSxriioqICS5YswZ49exAd\nHT3gdYFAgIaGBhQXFyM1NRWxsbHYuHEj8vPzsWHDBhBCoNFosHHjRuzduxdnz57FlStXuMZfo9Fg\n0qRJuHjxIh577DGzYz/77LMoLCzETz/9BI1GgzNnzgCgT40XFRUhOzsbH330EV577TWLvptmFxqN\nBj4+PigtLYWTkxOOHTsGAPjXv/6FkJAQlJWVobu7m3u/g4MDjhw5guLiYuzYsQPp6enca+Xl5Xjm\nmWdw9epVKJVKnD9/HgCwbt06BAUFoaSkBOfPn4e3tzeUSiUOHTqEo0eP4tSpU9i3bx/3RDaDMRhM\nSBjjhrt37yI1NRX79u37xfslL774IoRCIWbOnInw8HAEBQXB2dkZfn5+KC0txYkTJ3Dt2jXMnDkT\n0dHRuHjxIicKQqEQy5cvt3jcqqoq/OEPf0BERASOHz+OrKwsAFQgnn76abi4uCA2NhaEENTW1v5i\nt5StrS1efvllAEBiYiLX+GdlZSEtLQ0CgQBpaWnc+21sbJCRkYHHHnsMycnJKCoq4pZU8PX1xbx5\n8yAUCjFnzhycP38evb29OHPmDCdq9vb2cHBwwOHDh1FYWIiYmBjExcWhrq4Op0+fvs//AONRZcwu\ntctg9MfNzQ0ymQy5ubncNNhpaWkoKSmBr68vd1Xv5uYGgDaehrphW6vVQq/XIykpCf/9738H/A0H\nBwe4urpa/Pvvvvsu/va3v+GLL75ARkYGLl68yL3WfwIJgUAAsVgMrVbL7TN0XwGASCSCWCwGANjZ\n2ZllH5Ymo8jOzkZubi6ysrLg5OQET09PTkj6f8fOzk7uOP2PpdfrsXz5cvz973+3+B0ZDEuwjIQx\nbrC3t8dXX32Fzz//HPv37wcAZGZmoqSkhBORoRAIBEhOTkZubi43A2pzczOqq6uH/GxdXR2Cg4PR\n0tKC/fv3cxkHIQRff/01Ojo6UFRUBKFQCB8fH8THxyM/Px89PT24evXqfa37sGDBAuzZswd6vR67\nd+/m9tfW1sLX1xcuLi44cOCAmSj1hxACOzs7zJ07Fzt37gQhBFqtFl1dXXjhhRdw+PBh7vvW1tZC\npVIN6Rfj0YYJCWPcIBAI4OjoiGPHjmHLli1DiodAILDYvSQWi7Fz506sW7cOkZGRSEpK4taK+aXu\nqPXr12PRokWYP38+5s6da/Z3oqOjMWfOHLz33nvYuXMnAEAmkyE5ORlyuRz/+Mc/kJCQYPYZS36+\n//77uHbtGsLCwiASibj9qampaG1tRWhoKPLy8hAWFmbxWKbbH374ISorKxEVFYVZs2ahoaEBfn5+\nSE9Px6pVqxAZGYnnnnvO6uvJM/gPm7SRwWAwGCOCZSQMBoPBGBFMSBgMBoMxIpiQMBgMBmNEMCFh\nMBgMxohgQsJgMBiMEcGEhMFgMBgj4v8BE6+rBA5PeH4AAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 41 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(transcripts_all[:,0], transcripts_all[:,1] / transcripts_all[:,2][-1], label='all data')\n", "plot(transcripts_dn[:,0], transcripts_dn[:,1] / transcripts_dn[:,2][-1], label='diginorm')\n", "plot(transcripts_trin[:,0], transcripts_trin[:,1] / transcripts_trin[:,2][-1], label='Trinity norm')\n", "title('k-mer abundance plot of normalized data (zoom)')\n", "xlabel('k-mer abundance')\n", "ylabel('Fraction of total k-mers with that abundance')\n", "legend()\n", "axis(xmax=30, ymax=0.3)\n", "savefig('normhistzoom.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEVCAYAAAD6u3K7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdYU2f7B/BvAickIoqgDEVBRQUHoChKHSD6otaCs26r\nbVXUt+KuXa+gdllbZ23BOqu0tf7qbpWKFhAHKOJWHIi0TiyIYWc8vz9STgkkhJEBeH+uK5ckZ92H\nxNw8W8AYYyCEEELKEJo6AEIIIbUTJQhCCCEaUYIghBCiESUIQgghGlGCIIQQohElCEIIIRpRggDg\n4uKC48ePmzqMahEKhUhLSzPKtera78mYv5uwsDC0bt0avXr1Msr19CU8PByTJ08GAGRkZMDKygr6\n7vlelc/N9u3b0bdvX71evyJSqRSurq4oLCw02jUrkpSUhP79+5s6DB4lCAACgQACgcDUYdR69fX3\nVNMvpb/++gsbN25EcnIyzp49q8fIDK/0+9mqVStIpVK9v8eG+tyUTm7V9c0332DUqFEQi8V6iqpm\nfHx8kJeXh6SkJFOHAoASRK2iUChMHQKphoSEBHTo0AE2NjYGv5a+PyMv8zjZ4uJirFu3DrNnzzZ1\nKGpmz56NlStXmjoMAJQgyrlx4wbatGmD3bt3a9weHh6O8ePHY/bs2bCzs0O/fv3w6NEjbNy4Ea6u\nrggICFDL/gUFBdiyZQt8fHzQp08f7Nmzh/9PuX37dvTp0wdhYWFwdnbGsmXLyl0vKSkJvr6+aNKk\nCXx9ffH1119DLper7RMfHw9PT0906tQJP/30E3/+sn9hpaenQygUQqlUAgD8/f3x+eefY9CgQXB0\ndMSCBQvw/Plzfv8zZ84gMDAQLi4u2LBhQ5XiEgqFiIqKgpeXF1xdXbFmzRq14//44w9MmjQJNjY2\n6NSpE1JSUgAA2dnZWLt2LTp16oQhQ4bg999/1/JOAVOnTsW8efMwcuRIODg4YMmSJfj777817puf\nn49vv/0WXbp0QWBgIA4dOgRA9X7PmjULZ86cgZWVldYv+aysLKxcuRLt2rXD6NGjERcXBwDYsmUL\n3nzzTSQmJsLKykrje1jyPi9fvhwtWrTA4MGDcebMGZ3nBv79vM2aNQuOjo7Yvn07/P398emnn+I/\n//kPmjVrhilTpkAmk2HmzJlo3rw53nzzTTx69Ig/x9y5c9GqVSvY29sjJCQEly5d0niPpT8fJb+P\nkodYLEbr1q35fY8cOYLg4GB06NABa9asQW5uLr+tos9NWbm5ufj000/h5OSEAQMG4PHjx2rbtcV+\n9OhRfPbZZ9i9ezesrKzQtWtXAMC2bdvQsWNHWFtbY/jw4fz7rMm1a9egVCrh7OzMx63tnuVyOXbt\n2oVevXrB19cXUVFRap/32NhYjBo1Cu3atcMXX3yB7OxsfptQKMTOnTvh5eWFFi1aYM2aNXj8+DEG\nDRoEJycnhIWFqZ2rZ8+eOHHiRIW/N6NhhLm4uLDjx4+z5ORk1qpVK/brr79q3TcsLIyJRCK2bds2\n9uzZMxYcHMzc3d3ZvHnz2LNnz9iKFStYQEAAv//8+fPZuHHj2L1799jFixdZ586d2e+//84YY2zb\ntm2M4zj2/vvvs+fPn7OCgoJy10tOTmaJiYlMLpezU6dOMWdnZ3bs2DF+u0AgYL6+vuzy5cssLi6O\nubi4sKNHjzLGGAsPD2eTJk3i97137x4TCARMoVAwxhjz8/NjLVu2ZDExMeyvv/5iPXr0YJs3b2aM\nMZadnc0aNGjAduzYwR48eMAmTpzIOI5jx48fr3RcAQEB7MaNGyw5OZlZWVmxO3fuMMYYu3DhArOz\ns2M//PADKy4uZnfu3GH3799njDE2YsQIFhoayh4/fszi4+NZ8+bN2e3btzW+F1OmTGGWlpZ8jOPG\njWPjxo1Ti+Hu3buMMcaWLl3K+vfvz1JTU9nx48eZi4sL++OPPxhjjG3fvp316dNH63vOGGNvvPEG\nGzNmDPvzzz/ZL7/8wmxsbNi9e/cqdfy2bduYSCRiy5YtY1lZWSwsLExt/4rOHRYWxjiOY+vXr2cF\nBQWsoKCA+fn5MWdnZxYbG8vu3r3LWrduzTp16sR27NjBfyaXL1/On3/Xrl0sKyuLZWdns8WLF7Pe\nvXvz28LCwvjPSNnPRwmZTMb8/PzYBx98wBhj7MCBA8zDw4OdOXOGPXz4kI0ZM4bfputzU9bChQvZ\noEGDWFpaGjtw4ACzt7dnffv2rVTs4eHhbPLkyWrn+/XXX1laWhorLi5mP/74I5NIJCw3N1fjtX/4\n4Qf2yiuvaNxW9p63bt3KPDw82Llz51hycjLz8vJi27ZtY4wxlpaWxqytrdmePXvYX3/9xcaOHcum\nTJnCn0sgELD//Oc/7Pbt2+zEiRPMzMyMBQQEsLi4OP79i4+P5/cvLCxkZmZm7NGjRxpjMyZKEEyV\nIJYuXcqcnJxYXFxchfuGhYUxDw8P/vmuXbuYubk5y87OZowx9tdffzGxWMwKCwuZUqlkLi4uLCMj\ng99/zZo1bPbs2Ywx1ReHRCJhRUVFlY71ww8/ZO+88w7/XCAQsC1btvDP33//fX576f/8jJX/AvD3\n92dz5szht3/22Wds7NixjDHGdu/erfYf9e7du0wgEGj9j64prl9++YV/PmjQIPbtt98yxhh79913\n2dy5c8ud48WLF8zR0ZHl5+fzr82bN4998cUXGq85ZcoUtRhTU1NZkyZN+PsrnSA8PT1ZdHS0Wryh\noaGMMdX7UNEXvFwuZ7a2tiw1NZV/beLEiWz16tWVOn7btm1qcT18+JBxHMdyc3N1njssLIy1adNG\n7Xz+/v587IwxNn36dNatWzf+eVRUFPPz89MYy4sXL5ilpSXLzMzkz68rQcycOZMFBQXxzydMmMCi\noqL45ykpKaxjx46Msap/bjp27Ki2bfLkyVp/lxXFrk2fPn3Y//3f/2nctnLlSjZ69GiN28rec3Bw\nMIuMjOSff/fddyw4OJgxxtjq1avZ+PHj+W23b99mtra2ap/DvXv38tvbtWvHFixYwD+fPn06CwsL\nU7u+vb09O3v2bIX3ZgxUxQRVPWxkZCR69+6Nfv368a9HRUXxxc2hQ4fyr3t4ePA/29vbw8nJCdbW\n1vzzoqIi/P3337h58yYyMjLg4eGBJk2aoEmTJggLC8OpU6f44z09PSESibTG9uDBA8ycORMeHh5o\n1KgR1qxZg8uXL6vt4+Xlxf/ctWtXteoLXUof6+DggAcPHgAAEhMT4enpyW9r06YNGjduXO24HB0d\n8fDhQwCq4njv3r3LxZKQkIDMzEw0b96c/31t2bIFCQkJGmMXCARqMbZv3x4ymQw3btxQ208qleLy\n5cvw9vbmX/P29sbJkye1/2JKuXHjBoqKitC+fftqHQ8AnTp1glCo+u/m6OgIuVyOJ0+eVOrcPXv2\nLHe+0vdtb2+v9pm0s7Pj30dAVcU1dOhQNGvWDK1atUJBQQGuXLlSqbgjIyMRHx+PH374gX8tJiYG\ns2bN4t+j/v37Iz09HU+ePNH5uSntxYsXuHHjRrnPb2lVjT0hIQHjx4+Hs7MzrK2tkZSUVO5zWcLF\nxYX/TOq659OnT2v9/JTd5urqCrlcjmvXrvGvlX2/yj4v/X4VFRXh2bNnfNWXKVGCgOqLJjIyEvfv\n38eCBQv41ydOnAipVAqpVIpff/2V37eyOnToACcnJ1y/fh3Z2dnIzs5GTk4OLl68yO9jbm5e4Tk+\n/vhjyGQy/Pbbb8jJycH8+fP5NoQSJfX3AHDhwgW88sorAAAnJyc8efJE43669OzZUy3Ou3fvIicn\np0pxadO/f3+NX/q+vr5o1qwZnjx5wv++Xrx4gQMHDmg8D2NMLcbU1FRwHAd3d3e1/aysrODh4YHz\n58/zr50/f57/Y8DMzKzCxlo3NzdYWFggNTVV4/E1UZlzm5mZVfv8f/75JxYsWIAPPvgA9+/fR0ZG\nBiQSSaUap0+ePImlS5fiwIEDaNiwIf96QEAAvvvuO/49ys7ORl5eHuzt7XV+bkpr1KgR3Nzcyn1+\nS/6P6Yrd3Nxc7T4YYwgJCYGfnx8uXLiA58+fw8fHR+u9urm54e7du5W65969e2v9/JTddvv2bZiZ\nmaFTp05afrMVS0tLQ8OGDeHg4FCt4/WJEsQ/rKyscPToUcTHx+P999/Xul9l/mOVEAqFGDt2LJYs\nWYIbN25AqVTi7t27iI+Pr/Q5Hj58CBsbG9ja2iI2Nhbff/99uX22bt2Kq1ev4uTJk9i9ezdee+01\nAKr/yGfPnsWFCxeQmpqKjRs3Vvp+AgMDceHCBURFReHhw4dYtmyZWjKrTFxlr1NyrXHjxuHnn3/G\nzz//jOLiYty5cwcZGRmwtrZGnz59+C8EhUKBq1evqv3nKyslJYWPcfny5Rg8eDD/l3ppw4YNw6pV\nq3Dr1i3Exsbixx9/xPDhwwGo/hq8ffu2WkNraebm5hg6dCjCwsLw4MED7N+/H0ePHuWPr4nqnrvs\nF6M2mZmZYIzBwcEBUqkUH3zwAYqKinTG9eeff2LMmDHYuXMnXF1d1bZNnjwZX3zxBRISEqBQKJCZ\nmYmDBw8C0P25KevVV1/Fl19+iXv37uHw4cNq4yV0xe7t7Y3r16/zrxUXFyMzMxP29vYQi8XYtm0b\nEhMTtV7b3d0dZmZmuH//vs57HjZsGCIiIpCcnIyUlBRERETw79GwYcMQHR2NvXv34sGDBwgLC0NQ\nUJDGz2GJit6/xMREDBgwQOuxxkQJopTGjRvj2LFjOHLkCMLCwjTuo6lPd0XPw8PD0b9/f8yaNQs2\nNjZ4/fXX+Z4alekfHh4ejosXL8LJyQmrVq3CO++8o3aMQCDAjBkzMHHiRISEhODjjz/Gf/7zHwCq\n4n14eDjGjBmD8ePHY9q0aRXGWjoea2trREdHY9u2bfD19YWPjw+cnJyqFJe235uXlxeioqKwd+9e\n2NnZYeTIkXyvj4iICDg7O2P06NFo1qwZZsyYgRcvXmj83QgEAkyfPh2//PILunXrhhYtWmD9+vUa\nY3j33XcxfPhwjBw5Ep988glWr14NPz8/AEDHjh0xfPhwdOrUCXZ2dhqvtXr1anh6esLPzw/ff/89\n9uzZAxcXl3L3pi3Oin7v1Tm3tvet7PZu3bph9uzZCAgIQL9+/dC5c2e0bNlS67ElPx8/fhxPnz7F\nqFGj+GrWLl26AACGDBmC5cuX4+uvv0azZs3g6+vL99zT9bkpKywsDL1790afPn2wdu1ahIaG8tt0\nxe7n54f27dujdevW6N69OywsLLBu3TosX74crq6uOH/+PMaNG6f12hzHYf78+fwfThXd86RJkzB/\n/nzMnj0bs2bNwrx58zBx4kQAQOvWrbFnzx7s3LkTfn5+8PDwwOrVqzW+V5peK/seREREYMmSJVrj\nNiYBq8qfxITUIm+++SacnJywYsUKU4dC6qjc3Fx4eXnh6tWrtWKwXFJSEt57771a08214gpwQmox\n+tuG1FTDhg1x584dU4fB8/HxqTXJAaAqJlKH1depPwipLaiKiRBCiEZUgiCEEKJRpdsgCgsLa0Uj\nDlC1sQiEEEL+VaWu+rp2uHjxIoYOHcoPPrp48WKtmP2wpF99dR5jxjDEL9iLk11t8ce9P2p0LkM8\nwsLCTB4D3R/d38t2by/D/VWVzgTxySefYOXKlfxUEl5eXmqzTdZFEglQpOAgYkLkFeeZOhxCCKmV\ndCaIhw8fonPnzvzzoqIiNGjQwKBBGZpYDBQoRLBQCpAvyzd1OIQQUivpbIMIDAzk58LJyMjAhg0b\nMGzYMIMHZkgSCVBYxEGkqJ0Jwt/f39QhGBTdX91Vn+8NqP/3V1U6SxChoaFISUmBQqHAkCFDYG1t\njTlz5hgjNoMRi4ECOQeREsiT1b4qpvr+IaX7q7vq870B9f/+qkpnCaJJkyYIDw9HeHi4EcIxDokE\nyJeLwFEVEyFqbGxs1FZDI3VTkyZNkJWVVePz6CxBvPHGG2rLUGZnZ+Ott96q8YVNqaQEwckZJQhC\nSsnOzjZ5Txt61PyhrySvM0FcvnyZ78EEqDJTcnKyXi5uKhIJkC/jYK5gtbKKiRBCagOdCcLZ2Rm3\nb9/mn9+6davC6XvrgtIJgkoQhBCimc42iNmzZ2PIkCEYOHAgGGOIiYnBt99+a4zYDEYsBvJkIkoQ\nhBBSAZ0liEGDBuHy5csYMGAABg4ciCtXriAwMNAYsRmMRALkyTiYyRU0UI6QemD79u3o27cv/1wo\nFCItLa1Sx4aHh2Py5MmGCq1Oq9RcTA0aNEC/fv1QWFiIZ8+eAQBatWpl0MAMSSwG8oo4mMmVVIIg\n5CVXlbndpk6dipYtW740i1TpTBA//fQTPvroI5iZmUEkEvGvX7lyxaCBGZJEAkiLRBDKFZQgCCFE\nC51VTJ9++in++OMPpKam4sqVK/yjLhOLgbxiDkKZnHoxEVJHfP7553B1dYWtrS0mTpyIkydPVus8\nmZmZePfdd+Hg4IBRo0YhNzdXbfvrr78OR0dHtGzZEgsWLOCrqjZt2oQffvgBX3zxBaysrPgZJfQV\nV22kM0HY2trCysrKGLEYjUQCSAs5CORyKkEQUke4uroiISEBf/31F7p3744JEyZU6zyzZs3Co0eP\ncPHiRQQHB+Pbb79Vq2YaOnQo7ty5g6SkJDx79gxLly4FAMyYMQMTJ07EkiVLIJVK+SmI9BVXbaSz\nisnNzQ39+vXDsGHD+PEQAoEACxYsMHhwhiKRALlF/5QginJ1H0AIAQDoaymWasw8jdGjR/M/z5s3\nD2vWrEFycjK8vb0rfQ65XI6YmBhcuHABDg4OmDJlCrZs2aK2z9SpUwEAlpaWWLp0KXx8fKBUKiEU\nCv+JXT14fcRVW+ksQdjb22PkyJEwNzdHbm4upFIppFKpMWIzGLEYyC8UgpmZobiIShCEVBZj+nlU\nx8GDBzFy5Eg0b94cNjY2ePToES5fvlylc9y4cQNKpRJt2rThX+vWrZval/6XX36JgQMHwsbGBj16\n9MDz589x//59fnvZRm19xFVb6SxB1Kc5mEpIJEBhIQCOg6yQEgQhtV1eXh6mT5+OTZs2Yfv27bCy\nskLr1q3L/TWvi5ubG4RCIe7evYu2bdsCAJKTk9G6dWsAQGJiIlavXo3o6Gi0b98ef/75J9q3b89f\nx8zMDEqlUu9x1VY6E0RWVhb27NmD6Ohofn4PgUCAEydOGDw4QxGLgYICAByH4iJqpCaktpNKpcjN\nzYWjoyOUSiU+++wzPHz4sMrn4TgOAwcOxLJly/DFF18gJiYGFy9e5BPEgwcPYGlpCTs7Ozx69Ihv\nfyjh7e2NI0eOQC6Xw9zcXG9x1VY6q5g++ugj5OTk4Pr165g7dy6sra3h5+dnjNgMRiJRJQiBSARO\nAcgUMlOHRAipgIODAz777DNMnjwZnp6eKC4uRp8+ffjtAoFAreqnorEN33zzDezs7ODl5YV9+/Zh\n1qxZ/Lbhw4cjICAAXl5eCAoKwtixY9XOFRwcDKFQiBYtWmDkyJE646rrBExHWahr165ISUlBly5d\ncOnSJRQWFqJv374mnbBPIBDUqAjHGGBmBijsHeH2Zi6Slv6FxuLGeoyQkLqppv+3SO2g7X2s6vur\ns4rJwsICANCrVy9s374drq6udf4DJBCoShHMnIOVUII8WR4lCEIIKUNnFdOHH36I58+f491330V8\nfDxWrFiBr776yhixGZREAjBOhEYCCY2FIIQQDXSWIIKCggAA1tbW2L59u6HjMRqxGFCacWgkEFKC\nIIQQDbQmiNLrTpfUW5VurFm/fr1hIzMwiUSVIBoKzGlGV0II0UBrFZO3tze6d+8Oc3NznDx5Era2\ntrCxsUFCQgLMzMyMGaNBiMWAQsjBSiimEgQhhGigtQRRMtx848aNiImJga2tLQBVyWLIkCFGCc6Q\nJBJAUSBCQ1AVEyGEaKKzkdrMzAzPnz/nn+fk5NSrEoQlRDSjKyGEaKCzkXrZsmUYOHAgOnfuDAC4\ndu0aIiMjDR6YoUkkgOIFB0sBRyUIQgjRQGeCCAwMxJ07d3D27FkIBAL06tWLn9WwLpNIALlQBEtQ\nIzUhdVHJ6m6DBg3CtGnTcPPmTZ3HnDx5EtOnT6/UvqQSVUwAoFQqYW5uDplMhoSEBMTHxxs6LoMT\niwGZgEMDZk4lCELqoJLpNfr06VPpL/y+fftScqgCnSWI9evXY9WqVejYsaPakqP9+vUzaGCGJpEA\ncgGHBjBHvpwSBCF1UW2Z1aFk8r76RmcJYtOmTbh+/Tqio6Nx6NAh/lGR+Ph4uLu7o127dtiwYUO5\n7VFRUfD09ISnpycmTJiAW7du8dtcXFzg4eGBrl27wsfHpxq3VDliMSCDCBJGVUyE1AX379/HjBkz\n4ODggOnTp0MulwMAYmNj0bJlS437TZs2DZMmTcL//vc/jfu6uLggIiICvr6+aNWqFcLDwyGT/Tt5\nZ2xsLEaNGoV27drhiy++4Ge0BgChUIjvv/8eXbt2hZubG+Li4uDk5ISIiAi0bt0aHTt2xIkTJxAf\nH48ePXrA3d0dP/zwg6F/TXqlM0G0atWq3JqtusydOxeRkZGIiYnBxo0b8ezZM7Xtbdq0QXx8PC5d\nuoRBgwZhxYoV/DaBQIDY2FikpKQgKSmpStetCokEKAaHBjCjKiZC6oBRo0ahUaNGuHr1Kjp06ICf\nf/653CyuJftZW1vj6tWr6NSpE/bs2aN1dleBQICIiAisX78ex48fx44dO/gq9Hv37mHEiBEYP348\nYmNjceHCBcyfP1/t+M2bN2P79u24du0aGGN4+vQprl+/jvPnz2PChAmYOHEi1q1bh59++gnffPMN\nZsyYoZaAajudZaJGjRrBy8sLgYGBakuOahtJnZOTA+DfKqjAwEAkJiZi6NCh/D6+vr78z0OHDuWz\newljFBtVJQgOYqakBEFIJQmW6WfNURZWtf/jT548wbVr15CQkACxWIxFixZp/A4qu9/8+fOxdu3a\nCs/9xhtvoEePHgCAQYMG4dixYxgwYAD279+PIUOG8EuKfvzxx+jVq5fa8qPTp0+Hp6cnfy6lUomw\nsDDY2trizTffxNKlSzF16lS0bdsWbdu2RfPmzXHq1Cn4+/tX6f5NRWeCGDx4MAYPHqz2WkVzrZ87\ndw5ubm78844dO+Ls2bNqCaK0TZs28fM9lZw7ICAArVu3xltvvYXg4GCdN1EdEglQrOQgZnLkyXIM\ncg1C6puqfrHrS1JSElxdXSEWi/nXunXrpnG/du3a6dyvNC8vL/5nBwcH3L17FwBw+vRp9OrVi9/m\n6uoKuVyOa9euoUuXLgCAnj17qp3L0dGRH1Rsb28PAGoJxN7evk4tKKQzQZSMqDaEmJgY7Nq1C6dP\nn+ZfO3XqFBwdHXHjxg0EBQXBx8cHDg4O5Y4tvRSqv79/lTOyRAIUQQQbBZUgCKntevTogTt37qCg\noAASiQQAcOHCBXh4eJTb7/bt2ygsLOSTxIULF/gv9Kro3bs3EhMT+ee3b9+GmZkZOnXqxL9W2xum\nY2NjERsbW+3jdd5dyVJ8pQkEAqSlpWncv0ePHli8eDH//Nq1a+VKIABw+fJlzJw5E0ePHuWrrgBV\nBgYAd3d3BAcH49ChQ5g+fXq542u6VrZYrCpBWLBiShCE1HIODg7o1KkTwsLCsHjxYuzcuRNPnjwB\noF4lXbJfeHg4Fi1ahKioKDx+/Lha1wwODsby5cuxd+9e9OzZE2FhYQgKCqpT48DK/vG8bNmyKh2v\n807PnTvHP6KjozF+/HiNX9glGjdWLbwTHx+P9PR0HDt2rFwxLCMjA6NGjUJUVBRcXV351/Pz8yGV\nSgEAmZmZiI6O1phc9EEiAYqUHERKIfViIqQO2LNnD7KystC5c2fcvHkTY8eOBVB+udE9e/bg2bNn\n6NSpE65cuYKhQ4fy30sl+2tT+lxt2rTBnj17sHPnTvj5+cHDwwOrV6+u8DxlX6voWnWBziVHy1Iq\nlfD09MSVK1e07hMXF4eZM2dCJpMhNDQUoaGh/PQcISEhmDZtGvbt24dWrVoBUC0knpSUhLS0NIwc\nORIAYGtri4kTJ+Ktt94qH7QelkWMigIkny1F76FZGOAci6uzr9bofITUB/VtyVHGGBwdHfHbb7/p\nbIuoT/S15KjOBJGcnMxnwcLCQsTFxeHUqVM4fPhwFUPWH318iPfuBQo/XIEhQ5/Au+VvSJurucqM\nkJdJfUgQ8fHxaN++PTiOw9dff42NGzfi6dOnpg7LqIy2JvXChQv5BCEWi+Hr66tx8FtdIxYDz5Uc\nOAVoNldC6pHU1FSMGTMGMpkMo0ePRnR0tKlDqrN0JoiatIDXZhIJUCAXQaQENVITUo9Mnz69wnZS\nUnk6E8Tz58+xZcsWHD16FAAwZMgQvP3222qNPnWRRAIUKjiYy1XdXMsuqUoIIS87nb2YwsPDcf/+\nfXz22Wf47LPPcP/+fYSFhRkjNoMSi4F8OQehXAFzoTmKFEWmDokQQmoVnSWIo0eP4tq1a/wqcl27\ndlUbKFJXqaqYOEAmgyVniXxZPsTmYt0HEkLIS0JnCWLUqFFYv349srKykJWVha+//hqjRo0yRmwG\nJRYD+TIRIJOhAdeAxkIQQkgZWksQDRs25Ovk8/LysGjRIgCqfsWWlpb45JNPjBOhgUgkQJ6MA4qL\n0YBrQA3VhBBShtYSRG5uLqRSKaRSKZRKJRQKBRQKBZRKJT/auS5TlSD+qWISWVKCIKSeOHnypNqE\nofra92VU6UlFrly5gvj4eP5R10kkQG5xqSomGgtBSK3VsGFDWFlZwcrKCkKhEA0aNOCf//jjj2r7\nVmVZ0bL7uri44MSJE3qNvS7T2Ui9b98+rFixAmlpaWjdujUuXbqEgQMH4vfffzdGfAbDcar1IJRF\nxXwjNSGkdiq9aFnr1q2xZcsWBAQElNuvpkt/1paR5LVlCVOdJYgNGzbwy/SlpKTg5MmTdX4MRAmB\niIOySEZtEITUUbGxsfwyn+3bt8fbb79dpWVFS+87efJkZGRkICgoCFZWVli1ahVee+01fP3112rX\n9PDwwIGX+7EsAAAgAElEQVQDB8rFkp6eDqFQiL1798Ld3R0eHh7YtWsXv10ul2PXrl3o1asXfH19\nERUVpbZsaun7eOutt7Bs2TKMHz8es2fPhp2dHfr164dHjx5h48aNcHV1RUBAgEFX3QQqkSBycnLQ\nqFEj2NnZISsrC71798bVq/VjYjuhxb8JgnoxEVI3PXnyBOfOnUN8fDw/KWhpFS0rWtrOnTvRqlUr\nHD58GFKpFIsXL8aUKVPUvuQvXbqEhw8fal0ADQB+/vln/P777/jyyy8xbdo0FBYW8udftWoVPz/U\nl19+qXbu0vexadMmMMawd+9e+Pj44MaNG2jSpAkGDBiAO3fuIDExEQEBAXj//fdr8qvTqVJrUmdn\nZ2P06NHw9/fHgAED1JYMrcsEIhFYUTE1UhNSWQKBfh56pFAoEB4eDgcHB7WV5EorWVa0Xbt2/LKi\nlREUFIRbt27xq8zt3LkT48aNq7D6591330XLli0RGBgIFxcXxMXFAQD279+P//73v+jevTu6deuG\n//73v9i3b1+F9+Hm5oapU6fC1tYWY8aMwe3bt9WWND19+jSKigw3yFdngti3bx+aNGmCWbNmYf/+\n/Vi2bBm2bt1qsICMyUzMgRVTFRMhlcaYfh56ZG9vr1alpEnpZUUdHR3x4MGDSp1bLBZjzJgx2Llz\nJxhj+OmnnzB58uRqXev06dPw9vbmt3l7e+PkyZMV3kfpFfPs7e3h5OTEL7Bmb2+PoqIi/P3335W6\nl+qo0tJIbdq0QZ8+fQwVi9EJLf5NENSLiZC6qaqNuRU1QpuZmZXbPmXKFERFRSEmJgYNGjQotwBa\nZfXu3Rvnz5/nn58/fx79+vXjn5e9j9owN1zdWTvPAMwkIrWpNgghLzdvb28kJyervebr6wuBQIBF\nixbhjTfeqPa5hw0bhoiICCQnJyMlJQUREREYPny41v1rQ2+qlzpBmEtoJDUhdV1llv4su6309tI/\nz5w5E4cPH4aNjY3a8qJvvPEGrly5gkmTJlU5lhKTJk3C/PnzMXv2bMyaNQvz5s3DxIkTtR5bNk5t\n+xiSzhXllixZgpUrV+p8zZj01Vd5sl8GNl3vjR1HPkTK4xREvla+BwQhL5PaMg6gtvn++++xefPm\nOjNIWF8ryuksQWgaEFfZHgC1nXkDEQRyqmIihGiXn5+PNWvWYPHixaYOxei0Johvv/0WXbp0QWpq\nKrp06cI/nJ2dMXDgQGPGaDDmEg4CWTGNgyCEaBQdHY3mzZujW7duCAoKMnU4Rqe1+X/ChAkYMmQI\n3nvvPaxcuZIvltjb20MikRgtQEPiGnAQyKmbKyFEs0GDBuH58+emDsNktCaIxo0bo3Hjxvjpp58A\nqEZU5+TkIDMzE4BqAF1dxzXgIFTQbK6EEKKJzjaIuLg49O/fHy1atEC3bt3g4uKCV1991RixGRxn\nKYKZopjGQRBCiAY6E8TKlSuxY8cOtG3bFk+fPsWuXbvUBnfUZWJL1TKqlkIxlSAIIaQMnUMQHz9+\njFatWsHS0hJ5eXmYMGECVqxYYYzYDE4sBhRCDg0gogRBCIAmTZrUihG8pGaaNGmil/PoTBA2NjaQ\nSqV49dVXMXr0aLRo0QLu7u56ubipSSSAXCiCJTjqxUQIgKysLFOHQGoRnQli//79kEgk+OijjxAb\nG4sHDx5UODy8Lvm3BGFOJQhCCClDZ4Jo2LAh/7O/v78hYzE6iQSQCzhYMTMUyguhZEoIBS/17COE\nEMLT+W0YExODgIAAWFtb82vANmrUqMJj4uPj4e7ujnbt2mHDhg3ltkdFRcHT0xOenp6YMGECbt26\nVelj9UksBuQCEYRyBcTmYhTICgx6PUIIqVOYDt7e3iwhIYEpFApdu/K8vLxYXFwcS09PZx06dGCZ\nmZlq20+fPs2eP3/OGGNs+/btbNKkSZU+9p+5oyodS0X27mXsUYPWjN25w5p+0ZQ9zX2ql/MSQkht\nVNXvTp0lCJFIBG9vbwiFlat6ycnJAQD069cPzs7OCAwMRGJioto+vr6+/LrWQ4cO5Vdcqsyx+iSR\nADJwgIxGUxNCSFla2yB++eUXAKov6+HDh+P111/nVzISCAQYOXKkxuPOnTsHNzc3/nnHjh1x9uxZ\nrWu4btq0iZ/jpKrH1pRYrJ4gaLAcIYT8S2uCOHToEN8f2sHBAQkJCWrbtSWIqoiJicGuXbtw+vTp\nKh8bHh7O/+zv71+tBnSJBCiGCCguphldCSH1TmxsLGJjY6t/Al11UCdPnqzUayWeP3/OvLy8+Ofv\nvPMOO3z4cLn9Ll26xNq2bctu375d5WMrEXalXLzI2CWJD2NnzrC+W/uyuPQ4vZyXEEJqo6p+d+ps\nWAgNDa3UayVK2hbi4+ORnp6OY8eOlVvDNSMjA6NGjUJUVBRcXV2rdKw+SSRAMStVxUSD5QghhKe1\niunMmTM4ffo0nj59itWrV/PTfWdmZsLGxqbCk65duxYhISGQyWQIDQ1F06ZNERmpWq0tJCQEy5cv\nR1ZWFmbOnAkA4DgOSUlJWo81FLEYyGT/rEtNM7oSQogarQmiuLgYUqkUCoUCUqmUf93Nza3CEgQA\n+Pn54caNG2qvhYSE8D9v3rwZmzdvrvSxhiKRAEXKf9elpkZqQgj5l9YE4efnBz8/P0ydOhUuLi5G\nDMl4+ARB3VwJIaQcnW0Q9TU5AKoqpkIFB1ZM61ITQkhZL/XEQ+bmgEwggqKQ1qUmhJCyXuoEAQBK\nMw6yvH9KEHIqQRBCSAmds7lmZWXh8OHDOHPmDAoLCwGoRlJv3brV4MEZAzPjUJyvaoN4nPfY1OEQ\nQkitoTNBvPPOO7C0tERAQAA4jgOAerXilNJcBHk+jYMghJCydCaIS5cu4dq1a8aIxSQYx0GeXwxL\nkS01UhNCSCk62yDGjRuHLVu28NVL9Y45B3kBdXMlhJCytJYgGjZsyFcl5eXlYdasWbCwsACgqmJ6\n8eKFcSI0MMaJ+ARBA+UIIeRfWhNEbm6uMeMwHREHRSHN5koIIWXprGIaMGBApV6rszgOikKqYiKE\nkLK0liAKCgqQn5+PzMxMZGVl8a8/ffpUbW6muk4g4qAoLKReTIQQUobWBBEZGYl169bh4cOH8Pb2\n5l93dnbGvHnzjBKcMQhEIrDCYjSk2VwJIUSN1gQxb948zJs3Dxs2bMCcOXOMGZNRCS04KIuokZoQ\nQsrSmiBOnDiBgIAANG/eHHv37i23XR9LjtYGQgsOymJqgyCEkLK0Joi4uDgEBASorU1dWn1JEAKx\nCCxbBgszCzDGkFecB0uRpanDIoQQk9OaIJYtWwYA2L59u7FiMQkzC9WCQQKBAH2d+yImLQbD3IaZ\nOixCCDE5nd1c27Zti4kTJyIiIqJeTrlhJubAZDIAQHD7YBy8ddDEERFCSO2gM0Fcu3YNM2bMwN9/\n/41Fixahbdu2GD58uDFiMwoziQgo/idBdAjG4VuHoVAqTBwVIYSYns4EYW5uDo7jYGZmBqFQiGbN\nmsHe3t4YsRmFmYQDZMUAgNZNWsPe0h5JD5JMHBUhhJieztlcGzVqhC5dumDBggWYNm0amjZtaoy4\njIaTcBDIZfzz4A7BOJB6AL4tfU0YFSGEmJ7OEsSPP/6Ivn374ptvvsG4ceOwdOlSxMTEGCM2ozDX\nkCAOplI7BCGECBhjrDI73rx5E7/99hvWrl2Lp0+fmnT6b4FAgEqGrdPJj+Ngs+Z/6PR3PABAyZRw\nWu2EuKlxaGfbTi/XIISQ2qCq3506SxCjRo1C27ZtERoaivz8fOzcuRPZ2dk1CrI2EVlyECj+LUEI\nBUIEdQjCoVuHTBgVIYSYns42iPfeew9du3aFubnOXeskrgEHYakEAai6u35x+gss8F1goqgIIcT0\ndJYgevToUW+TAwBwDUXlEkRA6wBcfHwRf+f/baKoCCHE9HQmiPpO1ICDubJY7TUJJ8GA1gPw2+3f\nTBQVIYSY3kufICwacjArU4IA/unNRKOqCSEvMZ0JIiEhgV9+9PDhw/j000/VFhDSJD4+Hu7u7mjX\nrh02bNhQbvvNmzfh6+sLsViMr776Sm2bi4sLPDw80LVrV/j4+FTlXqrFoiEHM1Y+QQxtNxTH7h5D\nkbzI4DEQQkhtpDNBzJo1C5aWlrh37x7ef/99CIVCTJ8+vcJj5s6di8jISMTExGDjxo149uyZ2nZb\nW1ts2LABixYtKnesQCBAbGwsUlJSkJRk+BHNFlaiclVMANDMshk623XGH+l/GDwGQgipjSo11YZA\nIMC2bdswe/ZsvPfee0hPT9e6f05ODgCgX79+cHZ2RmBgIBITE9X2adasGbp37w6O4zSeQ19jHCrD\noiEHcw0lCAAY1mEYDZojhLy0dCYIFxcX/O9//8PPP/+MCRMmQKFQoLi4/F/cJc6dOwc3Nzf+eceO\nHXH27NlKByQQCBAQEIDhw4fj4EHDfzmLrThwkEFTTioZVW3MhEUIIbWFzv6rUVFR2L17N3788Uc0\nbtwYGRkZWLx4scECOnXqFBwdHXHjxg0EBQXBx8cHDg4O5fYLDw/nf/b394e/v3+1rmcmEUGEYhQX\nAxYW6ts6NO0AS5ElLjy6AO/m3ppPQAghtVRsbCxiY2OrfXyFU23I5XIMGjQIx48fr/QJc3Jy4O/v\nj5SUFADAnDlzMHjwYAwdOrTcvsuWLUPDhg2xcOFCjedasGAB3N3dy7V56HOqDchkkIskyHsuR+PG\n5TcvPrYYDbgGWOa/TD/XI4QQE9HrVBsl7Q8VtTmU1fifb9n4+Hikp6fj2LFj6Nmzp8Z9ywaan58P\nqVQKAMjMzER0dDQGDx5c6WtXi7k5zKFAQb7mXxq1QxBCXlY6q5iaNGmCbt26ISAgAI6OjgBUWWj9\n+vVaj1m7di1CQkIgk8kQGhqKpk2bIjIyEgAQEhKCx48fo0ePHnjx4gWEQiHWrVuH69ev4+nTp/xa\n17a2tli4cCFatmypj/vUTiBAMTgUSmWAo6jcZl8nX/z14i9k5GSgVeNWho2FEEJqEZ2zuZZek7qk\neCIQCDBlyhRDx6aVXquYAOQLLZGR9ARu3Rtq3D51/1R0b94d7/i8o7drEkKIsVX1u7PS032npaWh\nTZs21Q5Mn/SdIKTm1rh3/B48/Jpo3L73xl5EnI/A75N/19s1CSHE2PQ+3XdsbCx69uyJgIAAAEBK\nSgqCg4OrH2EtJBdwKMzVPBYCAALbBuLsX2eRU5hjxKgIIcS0dCaIVatW4eDBg2jSRPXXddeuXZGW\nlmbwwIxJLhShWKp9bEdDUUP0de6L6LvRRoyKEEJMS2eCyM3Nhb29Pf9cKpWiUaNGBg3K2JRmHIrz\ntJcgANUaEQdSDxgpIkIIMT2dCWLYsGFYv3495HI54uPjMXPmTIwdO9YYsRmNQqg7QQR1CMKR20cg\n0zDzKyGE1Ec6E8Ts2bPRqFEjuLi4YOXKlXj11Vcxc+ZMY8RmNEpzEWT5FX/xN7dqDlcbVyRkJBgp\nKkIIMS2d4yDEYjGmTp2KqVOnGiEc01CacZDlaW+DKFGyRkT/1v2NEBUhhJiWzhLE1atXMWfOHHTv\n3h1dunRBly5d4OHhYYzYjIaZc5AV6K46Cu4QjAM3D9DkfYSQl4LOEkRISAhmzJiBOXPmQCQqP9K4\nPlCaiyDXUcUEAF3sukDJlLiWeQ2d7TobITJCCDEdnQkCAMaPH19vkwMAgOMgz9ddxSQQCDDMTTU3\nEyUIQkh9pzNBrFq1CpMmTcLgwYP5ifgEAgE/Z1K9YM5BXli53knB7YPx4YkP8UHfDwwcFCGEmJbO\nBPHjjz/i0qVL4DhOrRRRrxKEiIOiEm0QANDPuR9S/07FI+kjOFo5GjgwQggxHZ0J4ujRo7h27Vr9\nrmISiaAs1F3FBACcGYfBroNx+NZhTPeueG1uQgipy3T2Yurfvz/OnDljjFhMRsBxUFSyign4Z42I\nW7RGBCGkftOZIBISEtC/f3+0bNmy3nZzFVhwUBZVPkEMdh2MuPQ45BXnGTAqQggxLZ1VTL/99hv/\ns76n2a4tBBYiKHRMtVGatdgaPZ164ve7v2OE+wgDRkYIIaajswTh4uLCP6Kjo/mf6xOhiAOrZBtE\nieEdhmN/6n4DRUQIIaanM0GU9u233xoqDpMSWnBgxVWbhC+4QzAO3zoMuVJuoKgIIcS0qpQg6isz\ncdUTRMvGLdGmSRucvH/SQFERQohp6UwQhYWF/M8HD6p67mRlZRkuIhMQSkSArGpVTICqmmnfzX0G\niIgQQkxPZ4IYOXIkZDLVX9ctW7bEo0ePMHDgQIMHZkxmYg6oYgkCAEa4j8D+m/vrZcM9IYToTBAj\nRozAmDFjoFAokJ6ejkGDBuHzzz83RmxGYy7hAHnVE4R7U3eIzcVIeZxigKgIIcS0dHZznT59OoqK\nijBs2DDcv38fERER6N27tzFiMxpziQgCWdXHNAgEAgx3G479N/ejm2M3A0RGCCGmozVBfPXVVwD+\nHfvw559/wtPTE2fPnkViYiIWLFhgtCANzbwBB4GsekuJDncbjpDDIVjef7meoyKEENPSmiCkUikE\nAgH/fMSIERAIBMjNzTVKYMZkLuEgUMjAGFDqliull1MvZOZl4k7WHbjauBomQEIIMQGtCSI8PNyI\nYZiW0EIEC6EMRUWAWFzFYwVCDHMbhgM3D2DhKwsNEyAhhJgAjYMAAI5DA7NilOrRWyU0qpoQUh9R\nggAAjoPYTIaCguodHtA6AFeeXMGT3Cf6jYsQQkyIEgQAcBwszGTVLkFYmFtgkOsgHLp1SL9xEUKI\nCWltg5gzZ47WgwQCAdavX691e3x8PEJCQiCXyxEaGlruXDdv3sSbb76JlJQUfPLJJ1i4cGGljzUI\nkQhiYXG1SxAAMMJtBHZd3oVp3abpLy5CCDEhrQnC29ub78VUdqSwQEdXn7lz5yIyMhLOzs4YNGgQ\nxo8fj6ZNm/LbbW1tsWHDBuzfX77eXtexBsFxsBBWv4oJAIa4DsGMQzMgLZLCysJKf7ERQoiJaE0Q\nU6dOrdYJc3JyAAD9+vUDAAQGBiIxMRFDhw7l92nWrBmaNWuGX3/9tcrHGgTHQSSofhUTADQWN8Yr\nLV/B0TtH8Xqn1/UXGyGEmIjOkdRZWVnYs2cPoqOjkZ2dDUBVgjhx4oTG/c+dOwc3Nzf+eceOHXH2\n7NlKfclX5djS3XD9/f3h7++v8/xaiUQQC4qRXYMSBKAaNLc/dT8lCEJIrRAbG4vY2NhqH68zQXz0\n0UdwcXHB9evX8fnnn2PHjh3w8vKq9gX1Ra/jNPRQggBUa0R8cPwDFCuKITIT6Sc2QgipprJ/PC9b\ntqxKx+vsxXTmzBm8++674DgOwcHBiIqK4qf91qRHjx64efMm//zatWvo1atXpYKpybE1wnHgBDVr\ngwCA5lbN0aFpB8Slx+knLkIIMSGdCcLCwgIA0KtXL2zfvh3nz5+vcHrrxo0bA1D1RkpPT8exY8fQ\ns2dPjfuWPU9VjtUrkQgcal6CAGjQHCGkHmE6HDp0iGVnZ7Nbt26xKVOmsIEDB7ITJ05UeExsbCxz\nc3Njbdu2ZevWrWOMMRYREcEiIiIYY4w9evSIOTk5sUaNGjFra2vWsmVLJpVKtR5bViXCrpoLF9h9\nG08WGVnzU93MvMmaf9WcKZSKmp+MEEL0qKrfnYJ/DtIqISEBffr00fmaMZXMMKs3V6/ikd9Y/Lz0\nGubOrfnp3De6Y8fwHfBp4VPzkxFCiJ5U9btTZxWTpoFqRhm8ZkwcB3OmnyomAPwaEYQQUpdp7cV0\n5swZnD59GpmZmVi9ejWfdTIzM2Fra2u0AI1CJIK5smYjqUsb4TYCU/dPxacDPtXPCQkhxAS0liCK\ni4shlUqhUCgglUohlUqRm5sLNzc3fP/998aM0fD0XILo3rw7copykPosVT8nJIQQE9DZBpGeng4X\nFxfI/llxjeM4owRWEb23QTx5gnzXLnj/radYt04/p5z962w4N3bGkj5L9HNCQgipIb23QSgUCowd\nOxZt27ZF27ZtMW7cOKSlpdUoyFpHJIKZUn8lCODfUdWEEFJX6UwQn376KYKDg5GWloa0tDQMGzYM\nn3zyiTFiMx6Og1CuvzYIAPB38cfNZzfxSPpIfyclhBAj0pkgkpOTMX78eJibm8Pc3BxjxoxBcnKy\nMWIzHo6DUFHzkdSlicxEeLXdqziQekB/JyWEECPSmSCCgoIwb948XLhwAcnJyVi4cCGCgoKMEZvx\ncBzMFDIUFuixXQP/jKqm7q6EkDpKZyP1ixcvsH37dn5q7tdeew1TpkxBo0aNjBKgJnpvpAagNDPH\nEL8CRJ/QXyO8tEiKFqtb4M/5f6KxuLHezksIIdWh90bq6OhohIaGIjo6GtHR0ZgzZw6io6NrFGRt\nxMw5yAtkej2nlYUV+jr3xZE7R/R6XkIIMYZKNVJX5rU6j9N/ggBU1Uz7bu7T+3kJIcTQtI6kPnLk\nCH777Tc8ePAAoaGhaiOpmzdvbrQAjYVxIigKivV+3uAOwVh8bDGK5EWwMLfQ+/kJIcRQtCaI5s2b\nw9vbGwcOHIC3tzcYYxAIBHB2doavr68xYzQOA5Ug7Bvao7NdZxxLO4bX2r+m9/MTQoih6GykLi4u\nhkhUu1ZHM0QjtbxFK/SSncT5p856PS8ARJ6PxNG7R7FvLFU1EUJMR++N1LUtORiKwEIEZaH+SxAA\nMKHLBMSmx+Kh9KFBzk8IIYagM0G8LAQiDspC/bdBAKreTGM6jcHWlK0GOT8hhBiC1gQxefJkAMDa\ntWuNFowpCUQcmEwGpdIw5w/xDsF3F76DQqkwzAUIIUTPtCaI1NRU3L9/H1u3bkVWVla5R30j4Di4\ntZUhPt4w5+/m2A12lnaIvlv/xpAQQuonrb2YQkNDMXLkSKSmpsLb21ttm0AgqJczuo4KKkZkJODv\nb5hLhHiHIDI5Eq+2e9UwFyCEED3S2Ytp5syZiIiIMFY8lWKIXkzo1w/Sd1eg5SQ/3L4NNGum39MD\nQG5xLlqtaYXLsy7DqZGT/i9ACCEV0HsvpoiICEilUuzevRs///wzpFJpjQKstTgOVmIZhg8Htm83\nzCUaihpiXOdx2HJhi2EuQAgheqQzQezbtw/du3dHfHw8YmNj0aNHD+zbVw/784tEQHExQkKATZtg\n0MbqzSmbIVfKDXMBQgjRE61tECU2bNiAEydOoEWLFgCAhw8fYtKkSRgxYoTBgzMqjgNkMvTqBUgk\nwB9/AAMG6P8yng6eaGHVAkduH0FQh3o2bTohpF6p1DgIoVCo9rPe6/9rg38ShEAAhIQAkZGGu1RJ\nYzUhhNRmOhPE7Nmz0b9/f4SGhmLOnDno378//vvf/xojNuMSiQCZaiT1pEnAsWPAkyeGudTYzmNx\n5q8zyMjJMMwFCCFED3QmiNGjR+PMmTPo1asXXnnlFZw5cwajR482RmzGxXFAsWokdePGwMiRhmus\nbsA1wIQuE7D5wmbDXIAQQvRAZzfX2sgg3Vzffhvw9QWmTQMAJCUB48cDt28DQgNMSHL16VUM2jUI\n9+fdh7lQZ1MQIYTUmN67ub40/mmDKNGjB2BlBRw/bpjLdbbrDBdrFxy+ddgwFyCEkBoySIKIj4+H\nu7s72rVrhw0bNmjc5/3330ebNm3g7e2Nmzdv8q+7uLjAw8MDXbt2hY+PjyHC0+yfbq4ljNFYPaPb\nDGqsJoTUWgZJEHPnzkVkZCRiYmKwceNGPHv2TG17UlISTp48ifPnz2PRokVYtGgRv00gECA2NhYp\nKSlISkoyRHialSlBAMDEiaoSxOPHhrnkmE5jkPQgCenP0w1zAUIIqQGdCSImJgYBAQGwtraGlZUV\nrKys0KhRI6375+TkAAD69esHZ2dnBAYGIjExUW2fxMREjB49GjY2Nhg/fjxu3Lihtt0kzSIaEkSj\nRsDo0cBWA83SLeEkmOQxCd9d+M4wFyCEkBrQmSDee+89rFixAllZWZBKpZBKpXjx4oXW/c+dOwc3\nNzf+eceOHXH27Fm1fZKSktCxY0f+ebNmzfjJ/wQCAQICAjB8+HAcPHiwyjdUbWWqmEqEhADffWfY\nkdVbU7ZCpjDMYkWEEFJdOrvPiEQieHt7qw2WqynGmNZSwqlTp+Do6IgbN24gKCgIPj4+cHBwKLdf\neHg4/7O/vz/8azoFK8cBRUXlXu7eHbCxAX7/HRg8uGaX0KRjs45wtXHFwdSDGNVxlP4vQAh5acXG\nxiI2Nrbax+vs5rpkyRJcunQJr7/+OqytrVUHCQQYOXKkxv1zcnLg7++PlJQUAMCcOXMwePBgDB06\nlN9nw4YNkMvlmD9/PgCgbdu2uHv3brlzLViwAO7u7pg+fbp60Ibo5vr550B2NrByZblNmzYBR44A\nhpqCatflXfj+0vf4ffLvhrkAIYTAAN1cnzx5AgcHByQkJODw4cM4fPgwDh06pHX/xo0bA1D1ZEpP\nT8exY8fQs2dPtX169uyJX375BX///Td++OEHuLu7AwDy8/P52WIzMzMRHR2NwYb4s10TDW0QJcaP\nB2JjgYcGWlJ6dMfRSHmcgrtZ5ZMkIYSYis4qpu3VGE68du1ahISEQCaTITQ0FE2bNkXkP/1FQ0JC\n4OPjgz59+qB79+6wsbHBrl27AACPHz/mSya2trZYuHAhWrZsWeXrV4uWNghANR5izBhVY/VHH+n/\n0mJzMSZ7TMZ3F77D5wM/1/8FDIgxIDcXaNhQ1TWYEFJ/6KxievLkCVavXs2XGoKDg7FgwQLY2dkZ\nJUBNDFLFFBEBpKRoHfhw4QIwYgSQlgaYmen30gBw89lN+G33w5/z/4TITKT/C1STUqnq5nv/vvaH\nQqH6nbi5Ae7u6v+2bavKvYQQ06vqd6fOEsTnn38OOzs7vqFj69at+Oyzz7BmzZpqB1krVVDFBADd\nugF2dkB0NPCqAVYMdWvqBvem7th/cz/GdBqj/wtUglKpSoTR0UBcnCoZ/vWXam4qZ+d/H+7uqgb7\nkrQQVpoAABywSURBVOfW1qrmm5s3/31s26b6NyMDcHFRTxoDBwL/zB5PCKnFdJYgPD09cenSJf65\nUqlE165d1V4zNoOUIHbuVH0z/lPdpcnmzcChQ8CBA/q9dIkfr/yIzSmbcfwNA83vocGjR6oeWtHR\nqhlsmzUDBg1SrYXRvj3QsqVqfYzqKioC7twBbtxQJYyrV1XX69YNmDxZNSmilZX+7ocQol1Vvzt1\nJoi5c+fCyckJb731Fhhj2LFjBzIyMrBu3boaB1tdBkkQP/2k6qa0e7fWXXJzgVatgMuXAScDLCld\nJC9CyzUtceqtU2hn207/F4DqCzshQZUQoqOBP/9UJYNBg4DAQNX9GVpBgSrR7twJnDwJDB2qShYD\nBwLmNG8hIQaj915MS5YswaNHj9CnTx/07dsXDx8+xHvvvVejIGslHVVMgKohdtw4YIuBlpS2MLfA\nW13fwqcJn+r1vAoFsH8/8NprqhLCRx8BDRqoml2ePgX27FFNYmuM5ACoSiRjxqiSxO3bQK9ewNKl\nqtLKggWqpqC6N8cwIfUPTfdd4tAh1YCHCrrwAsClS6ov2nv3DPPX7ouiF/D41gMbX92Ioe2H6j6g\nAgUFwI4dwOrVqnaCuXOBIUNUA/9qo9RUVQ3frl2ApaWqVPHGG4Cjo6kjI6R+0FsV08qVK7FkyRLM\nmTNH40XWr19f/ShryCAJ4uhRYM0aVb2LDr16AR9+CAQZaEnpP+79gUn7JuHyzMuwbWBb5eMzM4Fv\nvlE9evYEFi8G+vSpO91QlUrg1Cng+++B//s/VTvFokWqRm5CSPXprYqpZK4kb29vdO/enX94e3vD\n29u75pHWNpWoYiph6GnA+7fuj9EdR2POkfLJuSK3bwOzZwMdOgAPHqh6Ih08CPTtW3eSA6BaoKlv\nX9UcWHfuqHpB9e+vSshxcVT9RIjRMB12795dqdeMqRJhV11cHGN9+lRq19xcxmxsGNu/X/9hlMgr\nzmPt1rdjP1/9Wee+p04xNmIEY02bMvbRR4w9fmy4uEwlP5+xiAjG2rVjrEcPxnbvZkwmM3VUhNQt\nVf3u1Lm3l5dXpV4zJoMkiDNnGPPxqfTusbGMtW/PWHAwY+np+g+HMcZOZ5xm9qvs2WNp+W/8ggLG\nfviBsVdeYax1a8Y2bFAlrvpOoWBs376X774J0YeqfndqbYM4cuQIfvvtN+zevRvjxo3j660yMzMh\nlUrx66+/Gq2UU5ZB2iCSk4Hp01UjxSqpqAhYtQpYu1ZVzz9/vv5HDb9//H1cz7yO/WP3QyAQ4OJF\nVS+qH39UjSWYMUM1wtsQo7tru9OngS+/VHWVnTkTmDNHNZiREKKZ3togmjdvDm9vb4jFYr7doXv3\n7pg5cyb27t2rl2BrlSq0QZSwsFB1GU1KUtWNd+2q+lefwv3CcedZGqau2Qlvb2DYMMDWFjh/XjXg\nbPTolzM5AMArrwB796oatDMzVaO0581Ttb8QQmpOZzfXFy9ewNLSEmb/fAspFAoUFRWhQYMGRglQ\nE4OUIG7eVH37pqZW63DGVOPs5s1TNaiuWlWzv2aVSlWy2bIFOJCYAtm4QGzxuYBxr7Z8aROCLg8f\nAl99pZrmY+xY4L33VFOBEEJU9D5QLjAwEAUFBfzz/Px8DBw4sHrR1WbVKEGUJhCoumNev64ajNa5\ns2ogWmVXomMMyMpSDRL75BOgXTvVuIUePYB7Z7riw4Gh2PH8bQiF1IVHm+bNVQkiNVU17qNbN+Dt\nt1U9oQghVVfluZgAoEuXLrhy5YpBA6uIQUoQGRlA796quSf04MoVYNYsVc759ltV9dOTJ6rZT9PT\n1WdDLXkuFKq6dPr6qr7Yunf/t3uqTCHDK1tfwdtd38bM7jP1EmN9l5UFrF8PbNyomkrkww9pLAV5\nuel9LqYZM2YgODgYr732GgDg0KFDOHDgADZv3lyzSGvAIAni0SPVt/jjx3o7pVKpGsn83nvAixeq\nSelKZkB1cVGfIbVkVtSKXM+8jn7b+iFxWiLa2rTVW5z13YsXqiSxdi3Qr5+q3cjT09RREWJ8ek8Q\n169fx+zZs/H06VMwxmBnZ4eIiAh+FThTMEiC+PtvVb1OVpZ+zwsgP1/1rz6abb46/RUOpB7AH1P+\ngJmQGiOqIi9PNcDxyy8Bb2/VvE/+/nVrECEhNaH3BFHi8ePHEAgEsLe3r3Zw+mKQBPHihWqRgn+W\nPK2tFEoF/Hf4Y4TbCCzwXWDqcOqkggLVNB7r1qmanubNUy0rKxabOrJqkstVVaN376oaXO7eVS3Q\nwXGqCcPK/lv2NYFA1Wdb26Ow8N+flUpV7wsHB9UkWQ4O6j83akQZtxYzSILIy8vD8ePH8fz5c/61\nN954o3oR/n975x4VZbX+8e87MAM44AVQQeSimCHXQRxG9JhIHeioKB3No5UrwzrW0spKXbXKky3X\n8WStMqt1SvmB2FHs1CGzXLmwTJQUgQI0ATXAEZKLBCiDMwxz2b8/NvPODAxyGxxm3B/WXu/9nWfP\nHvb33ft99rOtwLAIhEoFjBtH/xlGOFUtVZD9nwx5T+VhxnjWqT5YCKFzYHzwAR3+sm4dfW/k42Nr\nyyzQ2Wms/E2FoKqKvsCaOJFO3xccDEybRn2htVqaNJo7L/V66rNtKbm6mm8D1Ke4oYF2yzY0mK/r\ndOai4e/fsy/Vy4uJiI2wukCkpaUhPT0d1dXVmDt3Lk6cOIHk5GQcPHhwyMYOlmERCJ2OPk311+3I\nxnxS9AkySjOQvzYfzgI2icJQqaigL7Q//5x6O7/4In0lZTMUCiA/n44CzMujA18mTaKVv6kQBAcD\nU6aMnOZPe7tRNOrqqPNH9zlqtVoaW767cBjyxQRk2LC6QMyZMwe5ubmIjo5GWVkZrly5gg0bNuD4\n8eNDNnawDItAAHTEWVPTyI2HbQIhBEkHkhA2IQzv/vldJhJWoqWFBgn8+GNaV23cSIMEDvvYkxs3\n6ExOBkG4dIm+KJk3j6a4ONp94wjcutVTOORyOsetwSfZVABNl76+1N2PMSisLhBSqRRFRUVYtGgR\n9u7dCx8fH4SHh6OiomLIxg6WYROIzZuBY8do6O/hmDLOyjS0N+Dxrx6HQq1AZkomQseH2tokh0Gj\nAbKzaQT4GzfoRFErVgASiRUebgmhoXfz842i0NBAh4YbBEEqNXbp3EsYBgSZdqOZLtvagKlTqRta\nXBz9ziIj2VSE/WSgdWef36pUKkVrayuefPJJzJs3D0KhEMuWLRuSkSOWd96ho9zmzqVCETqyK1wf\ndx98v/p77P1lLx7Y9wA2zdmETXM2sdaEFRAKqSisXEkHL37xBQ1rIhBQoXj0UVpH9Ussbt2i8Vjy\n84Fz54CCAurzLJPRiTrWrwciIu7dmCmmcBztYvLyAmJjex5vb6dCUVxMv89PPqGtkVmzqFjExdHk\nNfB5VBg9uWMLghCC2tpaBHTNRalQKNDa2spv24pha0EY+Owz2pr46isqFnaA/KYcT3/zNG6pbyFz\naSbCJoTZ2iSHgxBaL33xBU1CIRWLFSto/c5xoO+yysupEBjStWu0u2j2bJpkMvo+gWEdWlup6J49\nS0WjoIB2RRkEIzaWPuxZO5KmHWLVLiZCCCIjI206atoSwy4QAO1mWr2aBkNasmR4P8tKEEKQVpyG\n1398HS/Nfglb5m5hrQlr094O1NaC1NRCnleLypO1aD1fAz9Si+lutfBS1oKb7AfOIAZxcTTuilBo\na8vvHXQ6oKyMCsbZszRSc3U1MH067SOMijIu77GWxrCMpF60aBGWLl06ZOOsxV0RCIB2CyxdCmzf\nDjz99PB/npW4dvManvn2GbSoWpCZkonwCeG2Nsl+0OvpmIIrV+h7gt9+o10aNTV0v0pFXTcNKSAA\nZLI/rqj88U2JPzJ/DIACHkhMpOE9HnqIek8zbIxKRUWjtJROLF9aCly4QF/8mwpGRAT1CnPQ9z9W\nF4gZM2bg8uXL8PLygk+XgzjHcbhw4cLQLB0Cd00gAFpBJCUBTz1FYzTYifsdIQTpJel47cRr2Cjb\niC1zt0DoxJ5iAVARqK83CsBvvxkFobqaPlXedx9N06dT75mAAJr6cMEkhAYLPH6cTm+el0d7N5KS\naIqNZe9TRwx6PfWeMgjG+fPAxYv0QcDgUtw9TZ0KuLnZ2vJBYzWBqKmpQUBAAORyucWbBgUFDcnQ\noXBXBQKgHiZ/+QvtMvj4Y7t6mVhzqwZ///bvaFI2Ye/ivZjpOxOcnYjcoNHraWREudxyqqmhT46m\nImBYnzYNEIutZopaTeeryMmh6do1ICGBikViIo3JxRhhaDT0N2JoPZomuZw6spiOR5k6labgYNpc\nHMH/X1YTiOjoaJSUlAAAli1bhuzsbOtYaAXuukAA1L3ukUeAMWOArKyRMzCpHxBCsK90H/5x8h8Q\ncAIkTEng0+TRI9+dF4TQLoLmZsuppYUuGxrMBSAoyHIKCLCqCAyEhgY6evv4cbrUammE2ZAQujSs\nBwYyd/8RiU5HWxgGwaiuNqaqKioOBsEwiMbUqbTbytfXZr87A8MiEKbr/eH06dNYt24dtFotXnjh\nBTz//PM9znnttdfw3//+F+PGjcPBgwcREhLS72ttIhAAfRx88kk6QvTIkWHrXM7NzUV8fLx1bqbX\n0yh17e0gHR2oVl7HqRuF+KHuJxz//RS8xN5ULIISsGDKAniP8r7z/bRaGn1QqaT3VSpp5Q3QlpWT\nE+1DMaxb2M49cQLx06fTir21lS57WzeIgKn7o5cXHcxouu3lRWMETZlCa1cbTmjV3/IjhI6xqKgw\npkuX6LKlhTZsDMIREkLrmsmT6QOsrcTDqr/NEciQ8mcYw2EqGgbhuHqVPh04OdGwKBMn0nAkhvXu\n+7y9hyWuldXHQQyGF198EXv27EFgYCCSkpKwatUqeHsbK57CwkLk5eXh559/Rk5ODjZt2oSjR4/2\n61qb4uJCWw8vv0z915OTe8ar6W2d42jlqtMZY+T0sp57+DDif/2VNnX7Smo19azpLalU9KnF3R2c\nSIRgtRrBHR1I7egA6ewEcbmNTuF+3HbahzauE7dcRRCJx0AsHgtXHQdRhwZOKjU4gxhotbTyFYvp\nctQoY5+sTmdMhjxZ2Jer1SJ+yhRayXt6UqH19KQV/P33m+83iIENK/yB0t9KhuOM9UL30xUK+i7D\nIByHDtHGUW0tLVY/PyoW/v6Wl+PHD09PBxOIO2D6ECOV9jxOCC3YxkYqFo2NxvWSEvPt5mYaF87T\nk4pF94ch031jxgDu7nRsjbu7MVmhK7xXgbhw4QI8PDwAACqVil+n3wOHtrY2i9fdunULAPDAAw8A\noDPSFRQUYNGiRfw5BQUFWL58OTw9PbFq1Sq88cYb/b7W5ggEdHjtl1/SJwO1mg6EMo16aRr90rAO\nGJ+kDRE1e9uur6ePkiIRdY80JLHYfFsopOJj+qPo/kNxc+v1cZPT68Gp1XBVqeDa0YHR7W0ory3G\nL1fPoPz3EtRq/kCt5ibq9LfgNsYLo8cFwGvsJPh4+MLXwxc+Yh/4evhiongixCIxXJxc4OrsChdn\nF7N1oUBo/t5j2zaaGL3i4UHHfs2a1fOYSgX8/jtNtbV0efEi9cw2bLe20uIfPdqYPDwsr4vFBCI3\nDZxdNBC6dMLZRQMnUSecRRoIhJ1wEmkgENL1qj9qkHM5F4TTgnBa6IgWWn3vSU/00BM9CCF0CdLr\nNkDrFg6c2bqlpYATDCg5cU69HxMYj8lvypF3LQ/OAmc4CZzgxDnxS0v7TK+1lDgYbeU4DpzYDZga\nBC54So889UCt7r1btbGRjrdpbqbd3wqF8aFQoQC5fdtYN5jWBwOkV4HQ6XQDvhkAFBUV8d1FABAa\nGopz586ZVfKFhYVYvXo1vz1+/HhUVVXh6tWrfV47IuA4OjpquLhbFahAQAWkqwUgxCRETQ9BFB4z\nO02j06BJ2YR6RT3q2+vR0N6AekU9LjVfwkn5STTeboRSo0SHtgNqrRpqnZpfdmg7oNPrIHIS8YKh\nylchfVd6v/55BZxR3Az/RKYViOn2QCAwNrMNlZOlfd3394eGXxpwZM+Rvm3o5TMs7Te1kRDCLwGA\njCYgoV2JEIwFMJroodHpoNXpcEOvQ71eB41OCx3RQW9I0EF/Wwei1IMjzhAQITi9CBwRgtOJAL0Q\n0IkAnRBEKwLRCaEraMBBZTWIzhnQ08QRmgSgybjuBA5OtAJEVwUJAdC1LTDZFhgqSM4gFAQAodsm\n6+b79QCnB4He4rpxn44uOT2AbktOx19DoIeqoBpf764Ega5LBHUgnK7rPEv7zD+XQG9xHzi9ofRM\n8kn4/Jr/MPrxe3bvSkG9/DYJ4KZVwUOtgntnE9zVgEcngLN939oUmzjcEUJ69IMN1LPG0T1x3nrr\nLVubYFXUXX8GFMdH9rwbQ6XhqPVmJrwbEGihgxaAqs9zUWw+LW9XVQf7iIPcN+pz12xsgXXer6q6\n0o0h3MPqAiGVSrF582Z+u6ysDA8//LDZOTKZDOXl5UhKSgIANDU1YerUqfD09OzzWgC2eUHNYDAY\n9xhW94UYM2YMAOqNJJfL8f3330Mmk5mdI5PJkJ2djebmZmRlZfHTl47tmpT5TtcyGAwG4+4wLF1M\nH3zwAdatWweNRoMXXngB3t7e2LNnDwBg3bp1iI2NxZ/+9CfMmjULnp6eOHDgwB2vZTAYDIYNIHbE\nqVOnSEhICJk2bRr58MMPbW2O1QkMDCQRERFEIpEQqVRqa3OGzFNPPUUmTJhAwsPD+X1tbW1kyZIl\nxN/fnyxdupQoFAobWjg0LOXvzTffJH5+fkQikRCJREKOHTtmQwsHT01NDYmPjyehoaFk/vz55ODB\ng4QQxym/3vLnKOWnUqlIbGwsiYqKIjKZjLz//vuEkIGXn10JhEQiIadOnSJyuZzcf//9pKmpydYm\nWZWgoCDS3NxsazOsxunTp0lxcbFZBbpz506yYcMG0tHRQdavX0/effddG1o4NCzlb9u2beS9996z\noVXWob6+npSUlBBCCGlqaiJTpkwhbW1tDlN+veXPUcqPEEJu375NCCGko6ODhIWFkStXrgy4/Oxm\nML/pGInAwEB+jISjQRzoBfy8efMwrtto88LCQqxduxYuLi5ITU216zK0lD/AMcrQx8cHEokEAODt\n7Y2wsDAUFRU5TPn1lj/AMcoPAEZ1DS5tb2+HVquFi4vLgMvPbgSit/EVjgTHcUhISEBKSgq++eYb\nW5szLJiWY0hICAoLC21skfX56KOPMHv2bOzcuRMKhf2781ZWVqKsrAyxsbEOWX6G/BkcYhyl/PR6\nPaKiojBx4kRs2LABAQEBAy4/uxGIe4EzZ87g/Pnz+Ne//oWXX34ZDQ325UvfHxzl6aw3nnvuOVy9\nehU5OTmoqqrinTPsFYVCgb/97W/YtWsX3N3dHa78TPMnFosdqvwEAgHOnz+PyspK/Pvf/0ZJScmA\ny89uBEIqleLSpUv8dllZGWbPnm1Di6yPr68vADoHx5IlS/Dtt9/a2CLrI5VKUVFRAQCoqKiA1FLM\nGjtmwoQJ4DgOY8aMwfr163H48GFbmzRoNBoNli1bhtWrV/MThjlS+VnKnyOVn4GgoCAsXLgQBQUF\nAy4/uxGI/oyvsGeUSiXfnG1qakJOTo7FQYL2jkwmQ0ZGBlQqFTIyMhxO5Ovr6wEAWq0WWVlZWLhw\noY0tGhyEEKxduxbh4eHYuHEjv99Ryq+3/DlK+f3xxx+4efMmAKC5uRnHjx/H0qVLB15+w/kW3drk\n5uaSkJAQEhwcTHbv3m1rc6xKdXU1iYqKIlFRUSQhIYGkp6fb2qQhs3LlSuLr60tEIhGZPHkyycjI\ncBg3SUKM+RMKhWTy5MkkPT2drF69mkRERJCYmBjy0ksv2a1XWl5eHuE4jkRFRZm5fDpK+VnK33ff\nfecw5XfhwgUSHR1NIiMjSWJiItm/fz8hZOBurn1OOcpgMBiMexO76WJiMBgMxt2FCQSDwWAwLMIE\ngsFgMBgWYQLBYDAYDIswgWCMaORyOSIiImxthkXWrFmD7OzsYbm3+yCmh2QwrA0TCAajG3p9/+ZG\nG85ZDR19xkSGfcAEgmE3VFdXY+bMmfjll196HIuPj8cbb7wBiUSC6OhoVFZWYvny5QgPD8enn37K\nn/fDDz/g0UcfRVxcHHbs2MHvd3d3x9atWyGRSHrE+EpLS0NsbCxiYmKwZcsWdHZ28sfOnj2LmJgY\nxMfH84HPMjMz8fzzz/PnLF68GKdPn+Y/Z/v27QgLC8Njjz2GlpYWAEBdXR3Wrl2LkJAQM7va29vx\n0EMPYebMmVi4cCFOnToFgLasQkNDsX79eoSGhuLZZ5+FRqMBANTW1uKVV15BdHQ0YmJicPXqVQDA\nl19+icWLF2PevHnYu3fvIEqAca/BBIJhF1y+fBnLly/H/v37ERMT0+M4x3FobGxEcXExUlJSEBsb\ni507d+LcuXPYsWMHCCFQKpXYuXMn/vOf/+DMmTP49ddf+UpdqVRi/PjxKC0txZw5c8zuvWzZMhQW\nFuLnn3+GUqnEyZMnAdDRuEVFRcjNzcU777yDZ555xqLtpq0BpVKJSZMmoaysDGKxGEePHgUAvP32\n2wgJCUFFRQU6Ojr4893c3HD48GEUFxfj008/xbZt2/hjly5dwl//+ldcvHgRcrkc+fn5AICtW7ci\nODgYJSUlyM/Ph4+PD+RyOf73v//h66+/xokTJ5CVlcWPGmYweoMJBGPEc+PGDaSkpCArK+uO7yNW\nrVoFgUCAuLg4hIWFITg4GO7u7vD390dZWRmOHTuG8vJyxMXFISYmBqWlpXxlLxAIsGbNGov3ra6u\nxhNPPIHw8HB89913yMnJAUAr/kceeQQeHh6IjY0FIQTXr1+/Y/eQs7MzHn/8cQBAQkICX6nn5OQg\nNTUVHMchNTWVP9/JyQm7d+/GnDlzkJycjKKiIj70vZ+fHx588EEIBALMnz8f+fn50Gg0OHnyJC9W\nIpEIbm5uyM7ORmFhIaRSKWQyGerq6vDjjz/2swQY9yrDMuUog2FNxo4di8DAQOTl5fGhilNTU1FS\nUgI/Pz/+Kdwwp7lIJOLXDdtqtRp6vR6JiYnYt29fj89wc3PD6NGjLX7+pk2b8Prrr+PAgQPYvXs3\nSktL+WPdAxFwHAdXV1eo1Wp+n6EbCQBcXFzg6uoKABAKhWatBUtBDXJzc5GXl4ecnByIxWJMmDCB\nF4juebx9+zZ/n+730uv1WLNmDd58802LeWQwLMFaEIwRj0gkwldffYXPPvsMhw4dAgBkZGSgpKSE\nF4e+4DgOycnJyMvL46NZtrS0oKamps9r6+rqcN9996G1tRWHDh3iWwiEEBw5cgTt7e0oKiqCQCDA\npEmTMHv2bJw7dw6dnZ24ePFiv+ZMePjhh7F//37o9XpkZmby+69fvw4/Pz94eHjg888/NxOb7hBC\nIBQKsWDBAqSlpYEQArVaDZVKhZUrVyI7O5vP7/Xr19HU1NSnXYx7GyYQjBEPx3EYNWoUjh49il27\ndvUpChzHWezmcXV1RVpaGrZu3YrIyEgkJibyc27cqVto+/btWLx4MZKSkrBgwQKzz4mJicH8+fOx\nefNmpKWlAQACAwORnJwMiUSCt956C/Hx8WbXWLLz1VdfRXl5OUJDQ+Hi4sLvT0lJwc2bNzFjxgz8\n9NNPCA0NtXgv0+1//vOfqKysRFRUFObOnYvGxkb4+/tj27ZtePbZZxEZGYkVK1agvb39jt8jg8GC\n9TEYDAbDIqwFwWAwGAyLMIFgMBgMhkWYQDAYDAbDIkwgGAwGg2ERJhAMBoPBsAgTCAaDwWBY5P8B\nuXYyWJfTlGoAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 43 }, { "cell_type": "code", "collapsed": false, "input": [ "missing_all = numpy.loadtxt('transcripts.fa.hist')[0][1]\n", "missing_dn = numpy.loadtxt('transcripts.fa.diginorm.hist')[0][1]\n", "missing_trin = numpy.loadtxt('transcripts.fa.trinnorm.hist')[0][1]\n", "\n", "print 'Missing', missing_all, 'true k-mers in the sequence reads'\n", "print 'Missing', missing_dn, 'true k-mers in the diginorm reads'\n", "print 'Missing', missing_trin, 'true k-mers in the Trinity norm reads'\n", "print 'Ratio between dn and trin', (missing_trin - missing_all) / (missing_dn - missing_all)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Missing 96.0 true k-mers in the sequence reads\n", "Missing 103.0 true k-mers in the diginorm reads\n", "Missing 363.0 true k-mers in the Trinity norm reads\n", "Ratio between dn and trin 38.1428571429\n" ] } ], "prompt_number": 50 }, { "cell_type": "code", "collapsed": false, "input": [ "1976319 / 3163778.0" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 52, "text": [ "0.624670567909632" ] } ], "prompt_number": 52 }, { "cell_type": "code", "collapsed": false, "input": [ "190815.0 / 3163778.0" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 53, "text": [ "0.06031238601444223" ] } ], "prompt_number": 53 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }