{ "metadata": { "name": "dens" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "ls" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "01-skiplist.ipynb ipythonblocks_demo.ipynb\r\n", "02-coinflips.ipynb \u001b[34mpyskiplist\u001b[m\u001b[m/\r\n", "03-hyper-log-log-counter.ipynb pyskiplist.tar.gz\r\n", "04-bloom-filters.ipynb \u001b[34mpython-skiplist\u001b[m\u001b[m/\r\n", "\u001b[34mLog-Log-Sketch\u001b[m\u001b[m/ reads-bias-nonrandom.dens\r\n", "README.rst reads-bias-random.dens\r\n", "dens.ipynb reads-nobias.dens\r\n", "france-pyskiplist-51aebed83570.zip stupid example.ipynb\r\n", "\u001b[34mhyperloglog\u001b[m\u001b[m/\r\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "ls *.dens" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "reads-bias-nonrandom.dens reads-bias-random.r3.dens\r\n", "reads-bias-random.dens reads-bias-random.r5.dens\r\n", "reads-bias-random.r1.dens reads-nobias.dens\r\n" ] } ], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "nobias = numpy.loadtxt('reads-nobias.dens')\n", "bias_rand = numpy.loadtxt('reads-bias-random.dens')\n", "bias_nonrand = numpy.loadtxt('reads-bias-nonrandom.dens')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(nobias[:,0], nobias[:,3], label='no bias')\n", "plot(bias_rand[:,0], bias_rand[:,3], label='random errors, with bias')\n", "plot(bias_nonrand[:,0], bias_nonrand[:,3], label='random errors, bias')\n", "axis(ymin=0)\n", "xlabel(\"Position in read\")\n", "ylabel(\"Average graph density within 10 steps\")\n", "title(\"Effects of errors on assembly graph density\")\n", "legend(loc='lower right')\n", "savefig('/tmp/bias.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEVCAYAAAABwEUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdcFGf+B/DPbKEvsKAIgrhgwVhBVIIVvZgYlRSxxG6M\nGnM5U9QU7y6COU9TjCb+vNSzRPRM5KKemghqlNg1il0sIIgCFqQssH33+f0x2YGVsiywFPm+X695\nLbszO893htnvPPNMeTjGGAMhhJDHmqixAyCEEGJ/lOwJIaQFoGRPCCEtACV7QghpASjZE0JIC0DJ\nnhBCWgBK9jb617/+hS5duqBt27YAgAMHDmDw4MGQyWQ4f/58I0dXUUFBAWbMmIHWrVvjnXfeaexw\nSCUyMzMhEolgMpkqHR8XF4epU6c2cFS227BhAwYNGlTr78+YMQMffPBBPUZUUffu3XHo0CG7ltFU\nSRo7gKZGoVDg/v37EIvFwmcvv/wyVq9eDZ1Oh0WLFuHcuXMIDg4GACxbtgwzZ87EjBkz6lSuSCRC\nWlqaMN/68vPPPyM3NxfZ2dlwcHCo13mThsFxXGOH0CA4jrP7sl66dEn4Oy4uDunp6YiPj7drmU0F\nJftHcByH3bt3Y9iwYRXGpaSkQCaTCQmZMYajR4/iq6++qpey7XF/25EjR9CnT596SfQGgwESieUm\nYzQaLXaM1tg6PbHPdlGZpvC/oXs87YeacWpo//79GDZsGO7evQuZTIZJkybB3d0dWq0WYWFh6NSp\nEwC+2eTzzz9Ht27d8Oyzz2Lv3r3CPBhj2LVrF1544QV4enqib9++uHPnDgYPHgwA6NWrF2QyGRIS\nEqBWqzFr1iwoFAp4e3tj8ODBVf4Qrl+/jjfeeAOBgYF48803cePGDQD8YfG6deuwYsUKyGQyHDhw\noMJ3DQYDtm7dimHDhiE0NBRr166FTqcDACQnJyMgIABff/01OnfujJkzZ2LJkiWYOHEiXnvtNfj5\n+eH7779HQUEBPv74Y3Tq1Aljx47Fb7/9Jsw/Li6uwvRXrlzBmDFj4OPjA19fXyxYsKDK9Z6cnIyY\nmBh06tQJn3zyCQoKCoRxIpEImzdvRmhoKDp27IhVq1ZVOZ9Tp04hMjIScrkckZGRWLNmDQwGgzD+\no48+QmhoKDw8PNCzZ09cvnwZAL+zfPbZZ+Hl5YWAgAB89tlnwncuXLiAuXPnIjAwEAsWLEBWVpYw\nTqFQ4KuvvkJkZCTatGmDd999F6WlpRg/fjz8/Pzw1ltvobi42CLGhIQEhISEICIiAomJiRbjzDXe\nUaNGYc2aNRbjevbsif/973+VLvfx48fx9NNPo3379li9ejUUCoWwHVT2v/n999+rXU8ikQgbNmxA\nr1690K1bN/zwww8VtssPP/wQ/v7+GDFiBI4fP17l/+TWrVuYM2cOfH19MXv2bItyarJ+v/76a0RG\nRiIwMBBxcXHQ6/UAAJVKVeG3U/57v/76KxITE7F8+XL8+OOPkMlkCAsLw3//+1/06dPHIoaVK1fi\nhRdeqHIZmhVGLCgUCrZ///5KxyUnJ7OAgACLzziOY+np6cL7F198kb3xxhvs7t277NChQ6xt27bs\nxo0bjDHGtm3bxjp16sR27drFjEYjO3/+PHv48GGl81mzZg2bNGkSKyoqYgaDgR05cqTKmNu3b8+W\nLl3K8vLy2PLly5lCoRDGzZgxg33wwQdVfveLL75gw4YNY5cuXWJpaWksKiqKffvtt4wxxg4ePMgk\nEgmbOXMmy83NZWq1msXGxjKpVMpWr17N1Go1U6vVbNq0aWz8+PHs9u3b7KeffmJeXl4sIyODMcYq\nnX7s2LFs9erVTKfTsdLSUnbixIlKY7t58ybz9PRkCQkJ7M6dO2zChAls+vTpFut+2LBhLDU1lZ05\nc4bJZDKWlpZW6bzOnDnDTp48yQwGAzt69Chr374927dvH2OMsUuXLrGQkBCWnZ3NGGPs6tWrLDc3\nlzHGWJ8+fdj27duZ0WhkhYWFLCUlhTHGWF5eHpPL5WzHjh2sqKiILVu2jPXv318oT6FQsLCwMHb2\n7Fl2/vx55u7uzvr27ct27tzJcnJyWEREBNu4cSNjjLGMjAzGcRwbPXo0u3nzJtu+fTuTy+Xs6tWr\nwjqcMmUKY4yxrVu3soiICKGcc+fOMW9vb6bX6yssc0FBAXNxcWGbNm1iOTk5bNq0aUwqlbJff/21\nyv9NdevJvM4jIyPZhQsX2G+//cYUCgVLTExkjDG2fv165uDgwJYsWcLy8/NZbGwsGzhwYKX/D8YY\nCw8PZwsWLGAPHjxgn376KXNwcBC21Zqs3169erFTp06x69evW/xuq/vtKBQKYfnj4uLY1KlThXFa\nrZZ5eXmx1NRU4bPQ0FC2bdu2KpehOaFk/4j27dszNzc35unpKQz//ve/GWN88qsu2SuVSubn58dU\nKpUw/s0332SffPIJY4yx8ePHs1WrVlVa7qPJfvXq1ezpp59mV65cqTbelJQU5ufnZ/GZv7+/kJRm\nzJjB/v73v1f5/QEDBrCjR48K77dv385GjhwpLC/HcSwrK0sYHxsby4KDg4X3BoOBeXt7s2vXrgmf\nTZ48ma1cubLS6RljbMyYMWzhwoXs7t271S7bypUr2cSJE4X3N27cYN7e3sxoNDLG+HX2008/CeOf\neeYZ9tVXX1U7T7O//e1v7C9/+QtjjLHz58+z9u3bs+TkZGHeZr1792YrVqxgBQUFFp9/++23bPbs\n2cJ7g8HAfHx82L179xhjfFIxrwPGGBs+fDgbM2aM8P6f//ynsOMyJ3tzEmKMsYkTJ7IVK1YwxiyT\nvVqtZnK5XNipLViwgL3++uuVLuOPP/7IBg0aJLxPT0+3KKey/01164kxfp2vXbtWeL9o0SJh/Pr1\n65lcLhfWYU5ODpNKpaykpKTCfO/evcucnJyYWq0WPmvXrp2Q7Guyfj/77DNh/Kuvvsree+89xlj1\nv53yyb78ejWbO3cu+9vf/sYY4ysBcrmc6XS6atdRc0HNOI/gOA7/+9//UFBQIAyvvPJKjb575MgR\nPHjwAG3btoVcLodcLse6detw5MgRAHyTxIABA2o0r1deeQVRUVEYPXo0evTogbVr11Y63dGjR9G7\nd2+Lz/r06SOUWZ3S0lIcO3YMo0aNEuKdMWMGjh07JkzTpk0btGvXzuJ7ERERwt+pqanQarXo3Lmz\n8Fl4eDgOHz5c6fQAsGrVKqhUKnTv3h0jRoywaPYp79ixYwgPDxfed+zYEQaDQWhiAYDQ0FDhbz8/\nP2RnZ1c6r+zsbMydOxc9e/aEu7s7Vq1ahQsXLgDgm0H++c9/4v3334e/vz8WL14MlUoFANi4cSPO\nnz+P4OBgjBs3Trjiav/+/di8ebOw3lq1aoXS0lKLKz169eplsR7Lv/fx8akQa/llCQsLq7QJxMnJ\nCePHj0d8fDwYY/jhhx+qvFLn1KlTFvMMDg6Gh4eHxTSP/m+qW081ibNbt24Qifi04ufnB4PBgHv3\n7lUaW8eOHeHk5CR8Vn47rsn6rep/X9PfTmWmT5+O//znPwCA+Ph4TJgwAVKptMbfb8oo2dejyMhI\ntG7dGvfu3RN2FEqlUmhPHTp0aJVJWCQSWbR9uri4YNGiRUhPT8e6deswf/58XLlypcL3Bg4ciJSU\nFIvPzpw5U6NL4FxdXREREYGkpCQh3sLCQot28UdPyHIcZ3ESr0uXLnB0dMS1a9eEz06fPm3RTvro\nSb/AwED861//wt27dzF+/HhMnDix0ssOBwwYgNOnTwvvb9y4AbFYjG7dulW5TFVdzbF06VLo9Xr8\n8ssvKCoqwttvv21R5uTJk3H8+HGcOHECe/fuxfr16wHwyWvjxo3Izc1Fjx49MGvWLADAsGHDMG3a\nNItKQUlJCcaOHVtlbMzKycezZ88Kf6ekpKB///6VTjd9+nRs3rwZ+/fvh4uLS4WEbdavXz+cO3dO\neH/z5k0UFRVZTPPo/8baerIlzur07dsXaWlpUKvVFvMys3X91ua3I5FIKvxPnnzySTg4OODQoUPY\nsmVLs7jktaYo2VfC2o+yKp6enhg4cCD++te/4tatWzAajbh06ZKQsF566SV888032LNnDwwGAy5c\nuID8/HwAfG24fGL7+eefkZaWBpPJBFdXVzg4OFjUgsxCQ0Ph4OCA5cuXIy8vD5988gkkEolQ67G2\nLFOnTsXixYuRkpICk8mE7Oxsi5PKj3p0fhKJBKNGjUJsbCyys7OxY8cOJCYmVntSa9OmTXjw4AEY\nY3B1dYWbm1ul0z3//PNISkrCtm3bkJ2djdjYWERHRws1x8piq2p5c3Jy4OXlBW9vbyQnJ2Pjxo3C\njuH06dM4efIk9Ho9nJ2dIZFIIJPJoNfrsXnzZiFBuri4QCaTAQDGjx+Pbdu2YceOHSgtLUVpaSl+\n/vlnlJSUVLnc1qxevRoZGRnYtWsX9u7di9GjR1c6XWRkJDiOw8KFCzFt2rQq5/f0008jJSUFW7Zs\nQW5uLj788MMKO+9HVbaeHrVu3TpcunQJhw8fxo8//lhlnNXx9fVFt27dEBsbiwcPHmDlypUWRwB1\nWb+7d++u0W8nPDwcV65cgVartfh86tSp+Mtf/gIHB4da7ciaKkr2lYiOjoZMJhOGmJgYYdyjNcdH\n33/99ddo3749xo4di9atW2POnDlQKpUAgOeeew4ff/wx1qxZA29vb8yePRsajQYAsHDhQqxYsQJy\nuRwJCQm4ceMGhg8fDg8PD8yePRtLly6t8hr8xMREZGdnIywsDLdv37a4ksPatcuzZ8/GzJkzsXjx\nYnh5eWH48OG4fv16tcv76GcrV65Er169MGTIEGzcuBEJCQlQKBRVTp+UlITu3bujTZs22LRpE779\n9ttKE3hQUBASEhIQHx+PIUOGoGfPnli5cqVNsZnFxcXh3LlzCAgIwKeffoq//OUvwjilUok5c+bA\ny8sLQ4cORb9+/TBlyhQA/I4pKCgIgYGBOHnyJD7//HMAgFwuR1JSEg4ePIjOnTujU6dOFjuQypQf\n92isHMdh0qRJGDFiBJYuXYr4+Hihaayy5Zo2bRouXrwoxFkZT09PJCYmYu3atXjyySfRu3dveHp6\nCk05lc23svX06DRz5szB5MmT8eqrr2Lp0qUYPnx4lfOrbn0kJCQgPz8f3bt3x9WrVzFhwgRhnK3r\nt3zZaWlpNfrtDBkyBJ07d0ZQUJDFVThTp07F5cuXq123zRHHaluNJYQ0mo0bN+Lf//63TXeDXr58\nGQMHDkR+fn6tb16y181/TYlKpYKvry/Onj2LDh06NHY49cZuNftr164hLCxMGDw8PLB69Wp7FUdI\ni6FSqbBq1aoaPf5i165dUKlUuH79OmJjY/GnP/2pxdyRW1uff/45hg8f/lglesCOd9CGhIQIJ3JM\nJhP8/f3x4osv2qs4QlqEpKQkTJgwATExMYiOjrY6/c6dOzFt2jR4eHhgxowZmDNnTp3Kf9x3FAqF\nAl5eXvjhhx8aO5R61yDNOHv37sWHH35Yo8sBCSGE1L8GOUH7ww8/YNKkSQ1RFCGEkErYvWav0+ng\n7++PK1euoHXr1paFP+aHhIQQYi+2pm671+z37NmD8PDwConezHxtdHMcYmNjGz0Gir/x46D4m9/Q\nnGNnrHb1c7sn+y1btmDixIn2LoYQQkg17JrsS0tLsX//fowZM8aexRBCCLHCrp2XuLq6Ii8vz55F\nNKqoqKjGDqFOKP7GRfE3nuYce2016h20HMfVuv2JEEJaqtrkTno2DiGEtACU7AkhpAWgZE8IIS0A\nJXtCCGkBKNkTQkgLQMmeEEJaAEr2hBDSAtj1pipCCHmsvPUW8MsvgKcn4OFR9urkBDg4WA5iMSCR\nlL06OQFeXvzg7c2/Ojvz8+W4skEm46evZ3RTFSGE1ARjgK8vsHUrn7iLioDCQv5VqwV0urJBqwWM\nRn4wGPhXlQooKADy84GHD/lBq+Xnax4AYM8e4Mknqw2lNrnT5t2HVquFo6OjrV8jhJDm7fp1PskP\nGdLYkdSK1Tb7iRMnQqlUwmg0IiIiAp06dcK6desaIjZCCGk6Dh8GBg1q7ChqzWqyv3LlCtzd3bF9\n+3aEh4fj+vXrWLt2bUPERgghTcfjnuxdXFygUqkQHx+PKVOmwMnJCcXFxQ0RGyGENB2Pe7KfN28e\nevfuDZlMhv79+yMzMxMeHh4NERshhDQN2dmAUgk88URjR1JrNl+NwxiD0WiEpB4uDaKrcQghzcIP\nP/DDjh2NHQkAO12NU1hYiLVr1yIxMREA8Oyzz+KVV16h2j0hpOU4dAgYPLixo6gTqzX7t956CyaT\nCdOmTQMAxMfHg+M4fP7553UvnGr2hJDmoEcPYN06oG/fxo4EQO1yp9Vk36VLF1y+fBlisRgAYDQa\n0a1bN1y9erX2kZoLp2RPCGnq8vMBhYJ/tcOdrbVhl56qYmJisHr1auTn5yM/Px9r1qxBTExMrYMk\nhJBm5ehRICKiyST62rJas3dzc4NKpQLHcQD4E7Surq78lzkOSqWy9oVTzZ4Q0tS9+y7g5gYsXtzY\nkQjscoK2pKSk1gERQkizd+gQ8NFHjR1FndXoEccnT57ER38sbFZWFk6dOmXXoAghpEkoLQUuXuSb\ncZo5q8l+2bJl+Pzzz/H9998D4Jt1/vznP9do5qWlpZg+fTo6d+6Mrl274sSJE3WLlhBCGtLJk0Cv\nXmWPIm7GrDbj7Nq1C0ePHkV4eDgAwMvLCzqdrkYzj42NRWBgIL755htIJBKUlpbWLVpCCGlIzfwR\nCeVZTfYBAQEWyT01NRWdO3eu0cz379+P48ePw8nJCQDoRixCSPNy+DDfYcljwGqyf/XVVxEdHY37\n9+/j5ZdfxuHDh/Hdd99ZnfGdO3eg0Wjw2muvITU1FWPGjMGbb74pJH6zuLg44e+oqChERUXZvBCE\nEFLv9Hq+GWfAgMaOBMnJyUhOTq7TPGr0bByVSoU9e/bAZDIhOjq6QsKuTFpaGjp37oz//e9/eOqp\np/Dqq6/iqaeeEu7EBejSS0JIE3byJDBnDnD+fGNHUoFdbqqaOnUqXFxcEBMTg3HjxsHJyQlTp061\nOuOOHTsiJCQE0dHRcHZ2xsSJE7Fnzx6bgiOEkEZz6NBj014P1CDZX7p0yeK9SqXClStXajTzTp06\n4eTJkzCZTPj555/x1FNP1S5KQghpaHv3Ao9Rzqoy2S9btgwymQwXL16ETCYThu7du2PKlCk1mvmK\nFSvw5ptvonfv3nBycsJLL71Ub4ETQojdqFTAiRPAsGGNHUm9sdpm//777ws3VNV74dRmTwhpivbs\n4e+a/e23xo6kUnZpsx89erTwyITdu3dj2bJlyM/Pr12EhBDSHCQlAc8809hR1Curyf61116Dq6sr\nMjIysGjRIohEIsyePbshYiOEkMaRmNjykr1EIgHHcVi/fj3+/Oc/4/3330dmZmYDhEYIIY3g1i3+\n2fVhYY0dSb2yelOVQqHABx98gISEBJw8eRJGo7HGj0sghJBmJykJGD4cENXoOZHNhtWl2bRpE4KD\ng7FlyxZ4eHggOzsb77zzTkPERgghDe8xbK8HangHrd0Kp6txCCFNicEAtG4NpKYCvr6NHU2V7HI1\nDiGEtBgnTwLt2zfpRF9blOwJIcTsMW3CASjZE0JImaQkYMSIxo7CLqpssy8pKcHXX3+NkydPCt0Q\n9u3bF08++STmzp0LNze3uhdObfaEkKbi4UMgOBh48ABwcGjsaKpVr232U6ZMQWZmJhYtWoTU1FSk\npqZi0aJFyMjIqPGzcQghpNnYtw8YPLjJJ/raqrJmHxgYiBs3bsDR0dHic7VajZCQEGRlZdW9cKrZ\nE0KaipdfBsLDgb/8pbEjsapea/ZhYWFYuHAhzp07B41GA41Gg7Nnz+Ldd99F2GN2ZxkhpIVjjH+k\n8WN6chaopmZfXFyMr776CqdOncLvv/8Oxhj69u2LiIgIvPbaa5DJZHUvnGr2hJCmICUFGDsWSE8H\nOK6xo7GqNrmTbqoihLRsjPGPRxg9utl0Lt5gN1WtX7++Nl8jhJCmZ+tW/gqcZtBWXxe1qtm3a9cO\nt2/frnvhVLMnhDQmpRLo2hX48UdgwIDGjqbGapM7q3zqZY8ePar80v37920qhBBCmqS4OL4Jpxkl\n+tqqMtnfv38fiYmJkMvlFcb179/frkGRisx7cSMzwmAyCIPeqIfOqIPWqIXWoIXWqIXGoEGxthhF\n2iIotUootUqU6krBcRzEnBhikRhiTgwXqQvauLVBG9c2wqurg2sjLykhDeTiRWDTJuDSpcaOpEFU\nmexHjRqFkpKSSi+zHDJkiF2DsgcTM6FEVwKNQSN8xoE/6y4WieEgdoCD2AFSkRQcx8FgMqBIU4Qi\nbREKNYUo1hZDa9QKCVZv0sNoMpbN648z+AaTAWq9GmqDGiq9Cmq9GirDH696FdQGNdR6tfB9g8kA\nIzNCZ9ShRFeCYm0xinXFKNYWQ21QV1gOESeCVCSFRCSBRCSBWCSGo9gRjhJHOEmchL/dHd3h7ugO\nD0cPuDu6w0XqAgYGo8kIIzPCaDKiVF+K+6X3ca/kHu6V3sO9kntwlDhC4akoGzwU6NmmJ3r79YaH\nk4ed/0uENBDGgD//GViyBPDxaexoGkSjX42z8dxGPFQ/RL46Hw/VD1GoKRSSpflVZ9TBxEwwMROM\nJiNMzASpWAoniROcJc5wljrDUewInVEHjUFj8V1z8lTr1XB1cIWTxAkcODDwi80YE5KtzqiDwWSA\nRCSBiZng7ugOTydPIWE6ShwhEUkgFUkhFUsh5sQV2s7EIr7G7Cxx5l+lznCR8K/mWJ0lzpCK+YQt\n5sRC4nZzcIPMUQaZgwwyRxmcJc4QcSJhR2JvjDHkq/ORWZgpDDcLb+L83fM4d/cc2sraIrxtOHr7\n9kawPBiBHoEI9AiEj6tPg8VISL34/ntgzRrgxAlALG7saGzWLC+9nPzTZHi7eMPb2Rtezl7wdPIU\nkmX5xGhufhBxInDga97mpK4xaKAxaOAgduB3AOUSqzlxukhdIOKsX3xkYibojXo4iB0ogZVjNBlx\nNe8qTuecxtm7Z3Gr6BayirKQVZSFYm0xAj0CESwPRgevDugo78i/enVEB3kHOEocrRdASEM5dQp4\n7jlg1y6gb9/GjqZWmlyyVygUcHd3h1gshlQqFR6oJhROV+M8FlR6FW4V3kJ6QTrS89ORXpCOtPw0\npOWn4VbRLQS4ByDEOwSdvTujh08PhPmFoVvrbrQTIA3LYACWL+dr9F9+CcTENHZEtdbkkn1QUBDO\nnDkDLy+vygunZP/Y0xl1yCjIwLWH13At7xou3L+As7lnkV6QjhDvEIT5hSG0TShCfUPRy7cXPJ08\nGztk8jjKyACmTgUcHfkmnICAxo6oTppksj99+jS8vb0rL5ySfYul1qtx8f5FnM09i/P3+HMCF+5d\nQCuXVujl2wvdWndD19Zd0bV1V3Rp1QUuUpfGDpk0RxoNEB8P/PWvwPvvA2+//Vh0JG6XZL9//34s\nW7YMKSkpMBqNQkFKpdLqzIODgyGTyRAUFISZM2fiueeeq3PAdqPXA/fuAXfv8sODB0BeHv+M67w8\nID8f0On4Q8Hyg8nEn9k3mfgB4E/4iMX8RiUWA05OgJtb2eDiwn9Xo+EHtRrQavnpzYN5Ho6OloNY\nzJdjNJa9mgeDoexzmQyQy8sGT0++XBcXwNmZf330Ua6MlcWlVpfFxxj/vBDzYF5fOh0/aLX898rH\nLRbz8y+/3G5u/PeKi/mbWYqL+eGP9WrS6ZCnzEW28g5uGQuQbriPVF0Oruqy4e7hgx7yLujm2Qld\nPDogRBYED7GL5XKb/yd6veWrmTl2jgMkEn4Qi/lXqZT/P5UfJBJ+mczLXf5/U35ZHRws/0eSKi9y\nIw3l0iXg3/8GNm8GevUCVqwAQkMbO6p6Y5dk36dPH3zxxReIjIyEyMY9Ym5uLvz8/JCamoro6Ggc\nOXIEvuX6duQ4DrGxscL7qKgoREVF2VQGAD4x3bvHv5qTkF7PJ6r8fD5hP3zI/52fzyeaoiJ+UCr5\nZF5UxHc07OsLtGnDX47VqhU/eHsDXl5lP2SJBEwsgd4kht4ogtYggt4ggk7PwWgEmMEIZjKBGYyA\nwQipUQNnYwkcDaVw1JdAoiuFARJoOGdomBNUJidomKOw0+BMRoCZYNIaoC3WQqfUQleihaFUC5Pe\nCJFEDE4sAicRg5OIYOLEMEICI8QwQAIT4+CkV8JVWwAXbQGctQVw1BRBqlNBrFdDqlNBolNBZNRV\neOiTUSSFQewEvcQJOhE/MPbHhmViwgbGOTpA7OwAiYsDpK4OkDpLIOZMEDMjRDBCBBPEeg0k2lKI\nVSXgSkuAkhIwqRR6Jxk0Du4oFbujmLnBIC5br5BIIBYBLqYSuBiUcNIpIdUUwahRQi0yoJhpUMTU\nKDTw83Jz8YSHixwebt7wdG0FqaMzmFgCk1gKk0gCk0gCiQMHMVduMzfvJM07B6OR32FptWU7OI2G\nH2fekZffoZff2RoMZTs88+DgALRrxw+Bgfzg7c3vUKRSfryDA/9ZQADg7w+4u9u+3RNLGRnAnj18\nTT4ri39k8cyZfIckzVxycjKSk5OF90uWLKn/ZN+/f38cOHAATk5OtQrSbP78+XjiiScwe/bsssJr\nuncqLQVu3ACuXeOH69eB27fLauFaLZ+cXVzKfkxSKZ+cvbwAb28wL2+UOHqjUOSFErEHlHBHAfNA\ngdEDeSYv3DO2grJULFQ4dbqKlebi4rJ9RFER//t3cuKLMb+aK4PlK4R6Pb8fUqn4V7W6rLIvk/GD\ni0vFo0uJhM8BHh784O7OL9qjlVeOK6ugmiudRiM/zjyY85Z5AMpezRjj5/HoMjk4WORiAEBBQdnB\nT14e/95c0Tfvb83LqlLxMTk78zm0VSs+DwYEAH5+/DzLHyxptUBhIb9fLijgX4uK+M/NByXOLgxG\njzSoPH+HxvsU9K1/h7H1OaDUB8jrAlF+F0gKn4CkKATa3I7gSvzgLhNZHGTIZGV/m9ezp2fZq5NT\nxcp8+e1ejyLtAAAgAElEQVShfLxabVnOF2tK0Up9G96q25AXZ8FDmQUnVT5EBj04ox4igx4igw5O\nJQ/gUpANxwd3AI4Da+vPH4W5uwPu7mDuHoBMBpOTC0wOzjA4OMPo6AKTgxNEDhKIHSTgHCQQO0oh\ncnECJ3MD5y6D2MMNnMyND5oxMNMflxmbGIw6I4w6A0w6I4x6I6DVwVFfAqn2jx2yeeP/Y4EZx1cm\nOIkYIgdJ2VGQeWMrPzg68r9DHx+bj24Y44tVqfihtBS4f5+vw5mH/Pw/DrA4BmdjCVx0hfDNv4KQ\nm3vQJXMPnDWFSG3/DC6GjMWtriPh5CaBszO/zZgPusz7WScnPjWY63Ours3iYZcCu9Ts33vvPZw/\nfx7jxo2Dp6enUNCYMWOqnbFKpYLRaIRMJsODBw8QFRWFxMREtGvXzjLgQ4fKatz5+XzmyM3lk7j5\nVakE69AR+uDO0ASGoMQ/BEWe7VHg6IuHUl88NHhAWcyhpITfSMyvSiWQnc0Pubn8j7hNG/6HbP5x\nu7vzP3oP/nclvHdwsNyOJZKy6cyDI11MUiOM8TsAlYr/UUmltZuPyVS28ygt5edr/vFKpQAnNiBH\nlYm0wqu49vAqrubxQ3pBOoo0RWgnUyDANRhtHILgIW4LGdrCxegHR50fJKp20Bd7oqiI39EUFvKJ\nu3wl3mTiE375HV/5ljZza45IVNYSZh7MO+Xy60St5rdRZREDK1LCTZkDN2MhZEwJd6aEjCkhgxJu\nIhVcOTVcRSq4cio4QQPOZITYZICY6SEx6eHINHBFCVxZCWQohhtKIALfrMjAgf1xA6HhjyNAI/fH\nkSAnhZLJUMzcoBK5QS12g4FzAPtjgTmTERLOCI6Z4CAywFGkh6PIAAeRHiJmhMhkBMeM4ExGODIN\nfNg9eCEfD+GNXPjhIbyhFznCIHKAUeQAg9gBWs4ZSsigZDIoTTIUmdwAvQHu4hJ4SkrgISmFp6gY\nXg7FkEuK4ckp4caK4WxQwkFdBEdNEQwSJ2idPVEoD8KtLiOQ+cSzyAsIBScWwWAoq1iZX807Y/MO\nWa0uO+jPy+P/x56eZduSeTD/n8q3akql/HZsrii4upZ9r3yLoLd32c6kdWs+95grOXXNHXZJ9jNm\nzBBmXp61J19mZGTgxRdfBAB4e3tj8uTJmDlzZoWAi0MHQO3shVJHLxRLvVAAL9w2+iFD7YerRX64\n9NAP1/JbQ28SQybjV6yrq2WyNifo8uPN07Rty6/ctm0pObdkpbpSZBRmID0/HRmFGcgtyUVucS5y\ninOQW5KLrKIsOIod0dGrozB0kHcQ/m7l0qpZ3XdhMlmeoqjJ9OZWLKOxbOcllfI7L8Ysd2AaTdnp\nCvMgHFmaDOAe3IfoXi5Y3kMYVDoY1DoYVToYSrXgNGqINSWQqIohVhVDpC6B2EkKkdsjGdT8wy7/\nIzcfetXzeRGVij96NB+ZmgeO448MnJzKXvX6skqleSh/BG1uQTbvSMxHwLm5fINEbi5/VGFu6fP3\nL2vNCwgAevfmF7U6Te5qHKuFcxx69WLC+UO5nF8JAQGWQ6tW/MbXjH5rpJlhjOGB6oFwf8CN/BsW\n9wzojXp08OqAAPcA+Lr58oMr/+rj6gMfVx+0dm0NuZO8We0USMMzGvkGi6wsvtXhzh3L1//7P/6c\ncnXqNdl//PHHeO+99zBv3rxKC1q9erVNBVVaeFO6GoeQahSoC5BekI6c4hzcLbkrDLkluXhQ+gD3\nS+/jful9qPQqtHZtjUCPQLRzbye8+rv7o41rG/i6+aKNWxvIHGS0UyC1Vq+POO7atSsAIDw8vMKM\naSMlLY3cWY4+zn2sTqc1aHGv9B5uF93GbeVtZBVlIa0gDYeyDuFuyV3hoXNGkxEKTwU6e3dGJ+9O\n6OzFvwa4B8DPzY+ePkrqXaM341DNnrREJboSZBRk4PrD67iRf0N4zVZmI7ckFw5iB/i5+cHf3R9B\nnkEIlgdbDN7O3lTpasHs0mafkZGB+Ph4HD9+HBqNRijowIEDtY/UXDgle0IqYIyhUFOInOIc3FHe\nEZ4+erOAH9Lz08HALE4id/bujBDvEHRp1QVy54p9UJDHi12S/ejRoxEZGYlhw4ZB+sc1cxzHITw8\nvPaRmgunZE+IzcyPojY/eC4tPw3X86/jWt41XM27CmepM0K8QyweQx3oEYj2Hu0RJA+Cg9jBeiGk\nSbNLsu/ZsycuXLhQp8CqLJySPSH1ijGG3JJcXMu7hszCTItHUWcWZuKO8g4C3AP4I4FWIego72ix\nQ/B08qTmoWagXpN9fn4+GGP47rvvoFarMW3aNIsuCqt6kqVNhVOyJ6RBlX8K6fWH15GWnyacSM4q\nyhJOHAfLgxEkD0KwJ3+OwNPJk+/JTSyFVCSFg9gB7o7ukDvL4SSp2931xHb1muwVCkW1e/iMjAzb\noquscEr2hDQpRZoiZBZmIqMwQzhHcLPgJpRaJXRGHfQmvdDvcZG2CAXqAohFYsid5PBw8oCD2EHo\nec3cq5uTxEnoVMj8t7n7TEcxPzhL+Z7dzB0XOUocodQqka/OF4ZCTSHUBr6zInOnRQaTAVKxVOhS\nVCqWwsvZCwHuAQiQBSDAPQD+7v5wFFe8o9JB7ABHiSP/KnaEi9Sl2fSxYJdmHI1GU+G5OJV9VhuU\n7Alp3hhjUOlVKNAUoEhTBIPJYDGYuwo1D+ZkrTVooTVqhVdzr3MqvQoqvQpaoxbuju7wcvaCl7MX\nvJ294eHoARepi7DDcJI4QSwSQ2/U831Dm/idUL46H3eUd3BbeRt3lHeQrcyG3qSvELfOqIPWqOVf\nDVqo9CpIxVKh1zxvF2/IneSQO8vh6eQJT0dP/tXJE3JnuTBO7iQXuiCViBrmiad2Sfa9e/dGSkqK\n1c9qg5I9IaSpYIyhVF+Kh6qyPrEL1AUo1BSiUFOIAk0BCjT8+wI1/7d5vEqvgtqghogTwVniDDcH\nNwS4B6CdRzu0c+cHhacCXVp1QUevjnU+gqjXm6pyc3ORk5MDlUqFlJQUMMbAcRzu378PR3rIDCHk\nMcNxHNwc3ODm4Ib2nu1t/j5jDHqTHhqDBkqtkj+6+OPmutvK2/jt1m+49vAabhXeQjuPdujSqguC\nPIPgL/MXmpv8Zf5o79neLldMVVmz37BhA77//nucPn0affqU3TnYvn17TJkyBU899VTdC6eaPSGk\nhdEZdbhZcBOpD1KRWZiJ7OJsflBm447yDn4Y+wP6+ferdh52acb56aefEGOnjnkp2RNCiO3qNdnH\nx8dj6tSp+OyzzyyuyjE358yfP79u0YKSPSGE1Ea9ttmrVCoAQHFxMd1kQQghzZzVZhy1Wg1nZ2f7\nFE41e0IIsZld2uw7duwIHx8fDB48GIMGDcLAgQPh4eFRp0CFwinZE0KIzezWU9WtW7dw5MgRHDly\nBL/88gvkcjnOnTtX60CFwinZE0KIzeq1zd7szp07OHr0KA4fPoxz586hW7duGDRoUK2DJIQQ0vCs\n1uxFIhH69u2LRYsW4fnnn6/Xk7VUsyeEENvZpRnn/PnzOHz4MA4fPoysrCx06tQJgwcPxqxZs+oU\nLEDJnhBCasNubfbFxcU4evQoDh06hE2bNgEAsrKyahdl+cIp2RNCiM3s0mbfp08faDQa9O/fH4MH\nD8bhw4fRvr3tz40ghBDSeKzW7O/fvw8fH59aF2A0GtGnTx8EBARg165dloVTzZ4QQmxWm9wpsjZB\nXRI9AHzxxRfo2rUr3YVLCCGNyGqyr4s7d+7gl19+waxZs6gGTwghjciuyf7tt9/Gp59+CpHIrsUQ\nQgixwuoJ2vDwcMycOROTJk2y6HDcmt27d8PHxwdhYWFITk6ucrq4uDjh76ioKERFRdW4DEIIaQmS\nk5OrzaM1YfUE7Y0bN7B+/Xps3boVffr0wcsvv4ynn37aahv8X//6V8THx0MikUCj0UCpVCImJgYb\nN24sK5xO0BJCiM3sdp09AJhMJuzevRuvvfYaRCIR5s6di/nz59foiZi//fYbVqxYQVfjEEJIPbDL\n1TgAfxft/Pnz8c477yAmJgYJCQngOA5PP/20TcERQghpHFZr9uHh4fDw8MCsWbMQExNj0dn4iy++\niO3bt9e+cKrZE0KIzezSjHPz5k0EBwdbfJaRkYGgoCDbI3y0cEr2hBBiM7s044wdO7bCZ2PGjLGp\nEEIIIY2ryksvU1NTceXKFRQWFmLbtm3CXuTBgweQyWQNFiAhhJC6qzLZX7t2Dbt27UJRUZHFVTTt\n27fHmjVrGiQ4Qggh9cNqm/3x48cRGRlpn8KpzZ4QQmxWrydoP/74Y7z33nuYN29epQWtXr26dlE+\nMh9K9oQQYpt6fZ59165dAfCXXpqvkTfPnK6ZJ4SQ5qXGd9CaaTQaODk51U/hVLMnhBCb2eXSy0mT\nJkGpVMJoNCIiIgKdO3fGunXrah0kIYSQhmc12V++fBnu7u7Yvn07wsPDcf36daxdu7YhYiOEEFJP\nrCZ7FxcXqFQqxMfHY8qUKXByckJxcXFDxEYIIaSeWE328+bNQ+/evSGTydC/f39kZmbCw8OjIWIj\nhBBST2w+QcsYg9FohERitd8T64XTCVpCCLFZvV56aabX63H8+HEcP34cGo1GKGjx4sW1i5IQQkiD\ns5rs582bh8zMTAwZMgRubm4NERMhhJB6ZrUZp2vXrrh06ZJdOg2nZhxCCLGdXa6zHzp0KA4ePFjr\noAghhDS+GtXsr169Cn9/f3h6evJf4jhcuHCh7oVTzZ4QQmxml56qMjMzK/1coVDYVFClhVOyJ4QQ\nm9mlGUehUMDR0RFHjx6FQqGAq6srJWhCCGlmrCb7b7/9FhMnTsSSJUsAADqdDlOmTLF7YIQQQuqP\n1WQfHx+PvXv3wtXVFQDg7+9Pj0sghJBmxmqy9/DwsLjsMisrCwEBAXYNihBCSP2ymuynT5+OyZMn\no7CwEEuWLMHo0aMxa9YsqzPWaDSIiIhAaGgonnzySaxatapeAiaEEGK7Gj0bJzMzEz/99BNMJhNe\neukltGvXrkYzV6lUcHFxgVarRXh4OHbs2IGOHTuWFU5X4xBCiM3s8mwcgL8iZ8GCBTYH5OLiAgAo\nKSmBwWCAo6OjzfMghBBSd1Umezc3tyr7muU4Dkql0urMTSYTwsLCcPnyZXz++ec1PiIghBBSv6pM\n9iUlJQCAZcuWQaPRYObMmQCADRs21LiGLhKJcP78eWRmZmLkyJEYMGAAwsLCLKaJi4sT/o6KikJU\nVJSNi0AIIY+35ORkJCcn12keVtvsu3TpgtTUVKGWbzKZhEco2GLhwoXo2LEj5s6dW1Y4tdkTQojN\n7HIH7YABA7BixQo8fPgQeXl5WLVqFQYMGGB1xnl5eSgsLAQAPHz4EHv37sXzzz9vU3CEEELqh9UT\ntP/4xz/w0UcfoX///gCAESNGYOnSpVZnnJubi+nTp8NoNMLX1xcLFy6En59f3SMmhBBiM5u7JazX\nwqkZhxBCbGaXZhxCCCHNHyV7QghpASjZE0JIC2D1BK1er8fx48dx/PhxaDQaAHx70eLFi+0eHCGE\nkPphNdnPmzcPmZmZGDJkCNzc3BoiJkIIIfWsRn3QXrp0yeIxx/VWOF2NQwghNrPL1ThDhw7FwYMH\nax0UIYSQxldlzb5Hjx4A+McjpKamwt/fH56envyXOA4XLlyoe+FUsyeEEJvVJndWmewzMzOrLUCh\nUNgcYHXzI4QQUjP1muzLy8vLQ1JSEjiOwzPPPANvb+9aB2lROCV7QgixmV3a7Ddv3ozIyEgcP34c\nx44dQ2RkJDZv3lzrIAkhhDQ8qzX70NBQJCYmwtfXFwBw7949PPPMMzh37lzdC6eaPSGE2MwuNXsv\nLy+o1WrhvVqthpeXl+3REUIIaTRWb6pq1aoVwsPDMWjQIDDGcOTIEQwfPhzz5s0Dx3FYvXp1Q8RJ\nCCGkDqw242zYsKHil/44hOA4DtOnT6994dSMQwghNrPb1Tj2QsmeEEJsV5vcabUZJzMzE9988w2S\nkpJQUFAgFHTz5s3aRUkIIaTBWT1BGxsbi7CwMBgMBmzfvh0jR47EnDlzGiI2Qggh9cRqM05YWBjO\nnj2LXr164fTp0wCAPn364Pz583UvnJpxCCHEZnZpxnF2dobRaMSQIUOwbNkyBAUF0aOOCSGkmbFa\ns//999/RpUsXqNVqfPnll8jOzsa8efPQs2fPuhdONXtCCLFZvV+NYzQa8f777+PTTz+tc3CVFk7J\nnhBCbFbvd9CKxWIcOnQIxcXFdQqMEEJI47LaZj9gwABER0dj7Nix8PPzA8DvVcaMGVPt927fvo1p\n06bh/v37aN26NebMmYNJkybVT9SEEEJsYrXNfsaMGfyEHGfx+fr166ud8d27d3H37l2EhoYiLy8P\n/fr1w/nz5yGTycoKp2YcQgixWZO+gzY6Ohrz58/H0KFDywqnZE8IITazy6WX5geeARCehxMUFISR\nI0ciJCSkRoWkpaXh8uXL6Nevn03BEWLm5eUl3MFNSEshl8uRn59fL/OymuzFYjEOHz6M5557Dowx\n7N69G8XFxdi2bRvGjRuHN954o9rvFxcXY8KECVi1ahVcXV0rjI+LixP+joqKQlRUlM0LQR5/BQUF\ndBRIWhxzRTs5ORnJycl1m5e1Zpx+/fohMTFReIZ9fn4+RowYgX379mH48OE4depUld/V6/UYNWoU\nRo4cibfeeqti4dSMQ2qIthXSElW13dul8xKxWAylUim8VyqV4DgOHh4e0Ov1VX6PMYZXXnkF3bt3\nrzTRE0IIaThWm3FiY2MxdOhQ9OjRAxzH4eLFi1izZg1KS0vx1FNPVfm9o0ePYtOmTejZsyfCwsIA\nAMuXL8eIESPqL3pCCCE1UqOrcQwGA06cOAGO4xAREQGJxOo+omaF06E5qaHHfVvJzMxEcHAwDAYD\nRKKKB9zLly/HzZs38d133zVCdKSx1GczDnVeQpqFx31bsZbsScvUoG32hBBCmj9K9oTUgUKhwNdf\nf43IyEgEBgYiLi7O4sKF5ORkxMTEoFOnTvjkk0+s3iuQkJCAkJAQREREIDExUfg8Li4OU6dOFd6P\nGzcOfn5+aNeuHebPn2/Rc9yRI0fw7LPPwsvLCwEBAfjss8/qcYlJc1WjZJ+bm4vNmzcDAB48eICM\njAy7BkVIc8FxHL7++musXr0av/76K77//nscOnQIAJCRkYEXX3wREydORHJyMlJSUvD2229XO79N\nmzYhMTERixYtwqRJk3Dt2jWhnPJGjRqFtLQ0nDp1Cnl5eVi8eLEw7u2338arr76KvLw8XL58GcOG\nDavnpSbNkdUzrd9++y3+85//ICcnB5MnT4ZOp8OUKVNw9OjRhoiPkCZv2rRp6Nu3LwDgmWeewb59\n+/CnP/0JO3bswLPPPouxY8cCAJYuXYonn3wSJpOpynb5t99+G0FBQQgKCsLWrVuxe/duhISEVGif\nNT+zytXVFYsXL0a/fv2E+ZpMJqSnp0OpVMLT01O4Go60bFZr9vHx8di7d69w96u/vz898pg0ORxX\nP0NthIaGCn/7+voiOzsbAHDs2DGEh4cL4zp27AiDwYDLly/XaF5hYWE4fvx4pdOtWLECTz31FLy8\nvNC3b18UFhbi1q1bAICNGzfi/PnzCA4Oxrhx4+qlC1HS/FlN9h4eHha1kKysLAQEBNg1KEJsxVj9\nDPVpwIABQr/NAHDjxg2IxWJ069atyu+cPXtW+DslJQX9+/evMM3JkyexcuVKrFq1Crm5ufj9998B\nQKj9d+vWDRs3bkRubi569OiBWbNm1dcikWbMarKfPn06Jk+ejMLCQixZsgSjR4+mjYeQGnjuueeQ\nlJSEbdu2ITs7G7GxsYiOjq720srVq1cjIyMDu3btwt69ezF69OgK02RnZ8PV1RU+Pj7Izc21aK/X\n6/XYvHkzioqKAAAuLi4WjxUnLZfVNvtx48ahb9+++Omnn2AymfDzzz+jXbt2DREbIc0Ox3HCydTg\n4GAkJCRgzZo1ePfddzFr1izMmTOn2u9OmjQJI0aMgKenJ+Lj49G5c+cK833hhReQlJSE0NBQtGrV\nCkuXLsWPP/4ozGfTpk2YN28epFIpBg8ejM8//9yOS0yaC7qpijQLtK2Qlqg+b6qyWrOXyWQVZhwc\nHIyRI0fi9ddfp/Z7QghpBqzW7JctWwadTidc6hUfH4+SkhK0adMG6enp+Ne//lX7wqm2RmqIthXS\nEjXos3GeeOIJpKamCu9NJhO6deuGK1euoEePHrh06ZJNBVoUTj9gUkO0rZCWqEGfjRMZGYkvvvgC\n+fn5yM/Px5o1axAREQGO4+rt6ZeEEELsy2rN/s6dO1i+fDn27dsHAHj66afx/vvvw9vbG2fOnMHA\ngQNrXzjV1kgN0bZCWiJ6xDFpcWhbIS1Rg16No9PpcPDgQezdu9fiiX3r1q2zqSBCCCGNx2qb/d//\n/nfs2rUL27dvR69evXDlyhW0adOmIWIjhBBST6wm+19//RX/93//B2dnZ7z55pvYs2cPfv3114aI\njZAW79Hn2JOaO3z4MLp06VLl+MzMTOEpoTUxY8YMfPDBB1WOl8lkyMzMtDXMBmM12YvFYnAch7Cw\nMOzbtw9FRUVQqVQNERshLd6jz7EnNTdo0CBcvXpVeK9QKHDgwIFaz6/8IysqU1xcDIVCUev525vV\nNvs5c+YgPz8fb731FhYuXIicnBz84x//aIjYCGlWDAZDi7gcubLltHXZG2NdtfST/NXW7E0mE2Qy\nGby8vNCnTx8kJyfj+vXrmDBhQkPFR0iTplAo8OWXX6J///7w9PSE0WjERx99hI4dO8Lb2xuTJ0/G\n4cOHhek3bNiAgQMH4sMPP4S/vz9GjBhh8cz6Bw8e4N1334Wvry9iYmJQUlJiUV5KSgpmzJgBhUKB\nDz74ADk5ORaxfPXVV4iMjESbNm3w7rvvorS0FOPHj4efnx/eeuutavuiuHDhAubOnYvAwEAsWLAA\nWVlZVS5neno6RCIREhIS0L17dwwfPhwAsHPnTgwfPhw9evTA119/LbQCmJtMHp1+/vz56NKlC+Ry\nOfr164f79+9bXefTp0/HypUrAfBPABWJRPjyyy8BAOnp6fD29gbAdwlpfmjj1KlTkZWVhejoaMhk\nMqxYsUKY344dO/DEE0+gZ8+e2LRpU7VlFxcXY8yYMfD19cV7772Hhw8fCuNEIpHQPeTPP/+MsLAw\neHh4YPjw4di4caPFfGqz3HXGrAgLC2Mmk8naZLVSg+IJYYw13W1FoVCwrl27skOHDjGNRsMYYywh\nIYHl5uYylUrFVq5cyQICAoTp169fzxwcHNiSJUtYfn4+i42NZQMHDhTGx8TEsClTprDc3Fy2YcMG\n5urqyqZOncoYY6y0tJS5ubmx7777jt2/f5+98cYbbMiQIRaxhIWFsbNnz7Lz588zd3d31rdvX7Zz\n506Wk5PDIiIi2MaNGytdjry8PCaXy9mOHTtYUVERW7ZsGevfv3+Vy5mRkcE4jmMvvPACS09PZxqN\nhh04cIAFBgayffv2sevXr7M//elPLDY2ljHGKkyvVqvZ7t27WVRUFMvLy2Mmk4mlpKQwpVJpdZ2v\nW7eORUdHM8YY27x5M+vQoQObMGECY4yxtWvXshdeeIExxtjBgwct1r1CoWC//vqr8N4c04QJE1hW\nVhZLSkpijo6OTK1WV1ru9OnTmaurK/v+++9ZdnY2e+mll9hLL70kjOc4jqWnpzPGGEtOTmaXLl1i\nBoOBJSYmMplMxm7cuMEYYzYtd1XbfW1+D1a/sWTJErZw4UJ28eJF9vDhQ2GoiZdffpn5+Piw7t27\nV154E/0Bk6anqW4rCoWCffjhh1WON5lMrF27duz06dOMMT7Zy+VyZjQaGWOM5eTkMKlUykpKSphe\nr2ceHh5CwmCMsUGDBgnJftu2bSwyMlIYV1paylxcXFheXp4Qy8qVK4Xxw4cPZ2PGjBHe//Of/2TT\np0+vNM5vv/2WzZ49W3hvMBiYj48Pu3fvXqXLaU6Uhw4dEj5744032KJFi4T3+/btYz179qxy+p07\nd7LevXuz33//vcr1V5m0tDQml8uZyWRic+fOZd98842Q1KdNm8ZWrVrFGKt5sj9z5ozwWUhICNuz\nZ0+l5U6fPp0NGjRIeH/t2jWL/2X5ZP+oKVOmsBUrVti83PWZ7K02mq1btw4cx+G///2vxec16XT8\n5Zdfxrx58zBt2jTbDzkIsQG3pH5OZLJY29t0IyIiLN7v3LkTGzZswIkTJ6BWq1FSUoILFy4IXRR2\n69ZN6MDEz88PBoMB9+7dQ2lpKUwmE4KDg4V5hYWFCfe3HD16FL179xbGubi4oFOnTjh27Biio6MB\nAL169RLGt2nTBp06dRLe+/j44ODBg5Uuw/79+7F7924kJCQIn+n1ehw6dEjoQ/fR5Xz0s2PHjuH9\n998X3oeHh+PixYsWTUflpx85ciSysrLw8ssvo7S0FK+99hoWLFhQbecuANChQwe4urri3LlzOHz4\nMD744AOsXbsW169fx6FDh/DWW29V+/1Hle8K0s/Pz6JprDyO4yzWb+fOnaHX65Gamlqh97HLly9j\nxYoVOHbsGO7evQudTicsV22Xu66sJvu6XEo0aNCgJn0pEnl81CZJ15fyJxpLS0sxe/ZsfPvtt9iw\nYQNkMhmCgoJqdGKwS5cuEIlESE9PR4cOHQDwbfRBQUEAgIEDB+Ljjz+2KOvGjRuVdl1oVpNyAWDY\nsGHw8vLCV199VeU0lZ1QLf+ZuRvGmJgYAMDp06fRo0cPyGQyoW27/PRisRivv/46Xn/9daSmpmLU\nqFHo2rUrRo0aZTXeIUOGICEhAXq9Hm3btsWQIUOwYcMGFBQUWCTv8sRicZ1O0DLGcO7cOeH9tWvX\nIJVK8cQTT1SYduHChejTpw9+++03+Pr6YsqUKULZdVnuurC6K9Fqtfjxxx/x+uuvA+D70dy9e7dd\ng+oPQL4AABHdSURBVCKkuSouLkZJSQn8/PxgMpmwfPnyKmuKj5JKpXjqqaewZMkS3L17F5s2bbJI\nLsOHD8fly5exbt063L9/H3//+9/Rt29f4YRkXUyYMAHbtm3Djh07UFpaitLSUvz8888VThBX5/nn\nn8eWLVtw4MABpKWl4dNPP8WLL75Y5fTJycm4ePEijEYj3NzcIBKJhC4U4+LiMHTo0Cq/O2TIEKxZ\nswaDBw8GAERFRWHNmjUYNGhQlZdHhoeH48yZM1aXo7odwtmzZ7F582bk5OTgww8/xIgRIyqtkefk\n5KBVq1bw8PDAzp07sXPnzhottz1ZTfaxsbFISUlBcnIyAKBt27b429/+Vm8BxMXFCYO5DEKaK19f\nXyxfvhxTp05Fr169oNPpLB4WWNm12uXff/nll/Dx8UFoaCi2b9+O1157TRjn6uqKAwcO4LfffkPf\nvn3h7OyMzZs3VxtP+XlXd524p6cnkpKScPDgQXTu3BmdOnXCxo0bq72u/NFxUVFRWLVqFZYtW4YX\nXngBzz//PN55550qp7979y7GjRsHT09PPP/885gxY4aQvG/fvl3tQxYHDx6MkpISYfoBAwZArVYL\n7ysrc+7cudi9eze8vLyEq3kqW76qlpnjOMyZMwc//fQTevfuDX9/f6xevbrS73322WfYunUrAgMD\nsWXLFrz66qs1Wu6qJCcnW+TK2rD6ILSIiAicPHkSYWFhOHv2LACgZ8+euHDhQo0KyMzMRHR0NC5e\nvFix8BZ+3SupOdpWWpawsDAcOHAAcrm8sUNpVA36ILSQkBChp3oAOHHiBMLCwmwqhBBCbGGuWJL6\nY7UZZ968eXjhhRdw584dDB06FDNnzsQbb7xRo5lPnDgR/fv3x/Xr19GuXTusX7++zgETQgixXY2f\nZ3/mzBkYjUb069ev/gqnQ3NSQ7StkJaoQZtxoqOjMXHiRDz//PNwdXW1aeaEEEKaBqvNOAsWLMDh\nw4fRtWtXxMTE4L///S80Gk1DxEYIIaSe1LgZx2Aw4ODBg/juu++QmJgIpVJZ98Lp0JzUEG0rpCVq\n0GYcAFCr1di5cye2bt2KlJQUTJ8+3aZCCCGENC6rNfvx48fj5MmTGDFiBF566SUMGTKk3p7hQLU1\nUlO0rZCWqD5r9laz9syZM3Hz5k188803GDp0KI4ePSo8OoEQYl/ULWFF1roTXL58OWbPnt3AUTV9\nVptxRowYgZSUFGzZsgVbt25FUFCQ8KAjQoh9UbeEtlu0aFFjh9AkVZnsr127hi1btuDHH39E69at\nMW7cODDG6Pk1hFSBuiVs2t0StnRVNuM88cQTSElJQVJSEg4dOoR58+ZBLBY3ZGyENHnULWHDd0to\nlpCQgJCQEERERCAxMVH4/NGmr3HjxsHPzw/t2rXD/Pnzha4DAeDIkSN49tln4eXlhYCAAHz22Wc1\nLr/ZqapXk+3bt7Px48ez9u3bs1dffZXt37+ftW/f3ubeUapTTfGEWGiq2wp1S9jw3RKa5zV69Gh2\n8+ZNtn37diaXy9nVq1cZY4zFxcWxKVOmWKzzkpISlpOTw6ZOncomT54sjOvTpw/bvn07MxqNrLCw\nkKWkpFgtvyFVtd3X5vdg9RvFxcVs06ZNbNSoUczFxYXNnTuXJSUl2VxQpYU30R8waXqa6rZC3RI2\nfLeE5nmV72Jw4sSJQrd/sbGxFsm+vBs3blis/969e7MVK1awgoICm2JoKPWZ7K1ejePm5obJkydj\n9+7duH37NsLCwvDRRx/Z82CDENtxXP0MtVBZt4RjxoxB27Zt4eXlhdzcXItHglfVLWHq/7d3/0FN\n138cwJ8T+SHpEUQUhUtUZPJznxAQ+qEip5w3GB0R4GXI8PK8Lk3T8sDoOjqzyxPBBLqCQ5Eoqz8S\ncYo7WIrYIBmjQyRULCIyNYmNH+fE9/cPjs8XFBSWY3za63G3Oz+ffbbP8z3ni/nZm/eruXnUtoRD\n7teWcMjdbQmHb7u5uaGjo2PUMahUKpSUlMDZ2RnOzs5wdXVFT08PTp06NeY4795XU1PDt14ExteW\nUKFQICUlBXPnzsUnn3wy5gyb0QzvSMVx3IjLYcPt3r0bkZGRcHFxQXBwMLq6uvDrr78CAA4ePAid\nToe5c+ciPj4eOp1u3OcXmglNmHdxccHrr7+OyspKc+UhxDSMPZybCUZrS5icnIwLFy7g77//xtNP\nPz3htoRD6uvr+T8///zzIzotPey2hK+99hpu3rzJ3wwGA99/9u5xjrZvqC3hkOFtCUc7fqg9388/\n/4zy8nLk5eVBqVSOKy8wchnk+vr6UV8HjUaDPXv2ICsrC52dnairqwPw/9fF19cXBw8eRGdnJ/z9\n/bFu3bpxn19ozNvhlhArQ20JJ6ctIQDk5OSgra0NZWVlqKiogEwmu+eYjo4OPPLII3Bzc0NnZycy\nMjL4+4xGI0pKSvh+HY6OjpPSHtBSqNgT8hBRW8LJaUsoEomwevVqREVF4cMPP0RxcTEWLFhwzzhj\nY2MREREBqVSK6OhoJCQkjMhx6NAheHp6QiwWQ6PRYO/evWOeU+jGvRCaWU5OvwJPxoneK9aF2hIO\nepjLJVCxJ4JA7xVijSZ1bRxCCCHCR8WeEEKsABV7QgixAlTsCSHECtCyc0QQnJ2dablfYnUe5mwk\nmo1DCCECM+Vm45w6dQoLFy6El5cX9u3bZ85TWYTQ1/an/JZF+S1HyNlNZdZiv2nTJnz22WdQqVTY\nv38/rl+/bs7TTTqhv2Eov2VRfssRcnZTma3YD6038eKLL+KZZ57BihUroNFozHU6Qggh92G2Yl9X\nVweJRMJv+/j44McffzTX6QghhNyH2b6gValUKCgoQGlpKQAgPz8fHR0dyMzM/P/JaXYFIYSYZKKl\n22xTL4ODg0eseNfU1ISoqKgRx9BMHEIImRxmu4zj5OQEYHBGzpUrV3Dy5MlRO90QQggxP7P+UtXe\nvXuxfv16GI1GbNy4Ea6uruY8HSGEkDGYderlkiVL0NzcjIsXL2Ljxo38fqHNv1coFHjiiSfg7+/P\n79Pr9ZDL5RCLxYiNjZ1QR5/J1t7ejmXLlsHX1xdLly7Fl19+CUA4Y+jv70doaCikUikWL16MrKws\nAMLJP2RgYAAcxyE6OhqAsPLPmTMHAQEB4DgOISEhAISTv6enB8nJyViwYAF8fHyg0WgEk72lpQUc\nx/E3Jycn5OTkwGAwTDi/RdbGEdr8+5SUFBw/fnzEvry8PIjFYrS2tsLDwwP5+fkWSvdgtra2yMrK\nQlNTE7799lvs2LEDer1eMGNwcHBAVVUVGhoa8MMPP6CgoACtra2CyT8kOzsbPj4+/MQEIeUXiURQ\nq9XQarWora0FIJz877//PsRiMRobG9HY2AiJRCKY7N7e3tBqtdBqtTh37hwcHR3x0ksvITc3d8L5\nJ73YC3H+/QsvvHDPGhW1tbVITU2Fvb09FArFlB7Dk08+CalUCgBwdXWFr68v6urqBDUGR0dHAIDB\nYMDt27dhb28vqPy///47jh07hnXr1vETE4SUH7h3QoVQ8qtUKqSlpcHBwQHTp0+Hk5OTYLIPp1Kp\nMH/+fMyePdu0/GySnTx5kiUmJvLbeXl5bMeOHZMdY8La2tqYn58fvy0Wi1lfXx9jjLGenh4mFost\nFW1CWltbmaenJ9Pr9YIaw8DAAAsICGA2NjZs3759jDFh/R28/PLLrL6+nqnVaiaTyRhjwsrv6enJ\nAgICmFwuZ99//z1jTBj529vbmbe3N0tOTmYhISFs165drLe3VxDZ75aSksL279/PGDPttacljk3E\nBDhtVK/XIyEhAVlZWZg5c6agxjBt2jTodDpcvHgRubm50Gq1gsl/9OhRuLm5geO4EZmFkh8Azpw5\nA51Oh48++ghbtmzBn3/+KYj8/f39+OWXXxAXFwe1Wo2mpiYcPnxYENmHu3XrFsrKyhAfHw/AtPfO\npBf74OBgXLhwgd9uamrC4sWLJzvGvxYcHIzm5mYAQHNzM4KDgy2c6P6MRiPi4uKwZs0ayOVyAMIb\nAzD4ReGqVaug0WgEk7+mpgZHjhyBp6cnkpKSUFlZiTVr1ggmPwC4u7sDABYuXIiYmBiUlZUJIv/8\n+fPh7e2N6OhozJgxA0lJSTh+/Lggsg+nVCoRFBSExx9/HIBp/3Ynvdj/V+bfh4aGorCwEH19fSgs\nLJzSP7AYY0hNTYWfnx/eeustfr9QxnD9+nV0dXUBAG7cuIGKigrI5XLB5N+5cyfa29vR1taGr776\nChERESguLhZM/t7eXuj1egDAtWvXcOLECURFRQkmv5eXFzQaDe7cuYPy8nJERkYKJvuQ0tJSJCUl\n8dsm5TfHtaUHUavVTCKRsHnz5rHs7GxLRJiQxMRE5u7uzuzs7JiHhwcrLCxk3d3dLCYmhs2ePZvJ\n5XKm1+stHXNMp0+fZiKRiAUGBjKpVMqkUilTKpWCGUNjYyPjOI4FBASwFStWsAMHDjDGmGDyD6dW\nq1l0dDRjTDj5L1++zAIDA1lgYCCLiIhgBQUFjDHh5G9paWGhoaEsMDCQvf3228xgMAgmO2OMGQwG\n9thjj7Hu7m5+nyn5Ldq8hBBCyOSgL2gJIcQKULEnhBArQMWeEEKsABV7QgixAlTsicXZ2NiA4zgE\nBQVh27ZtMBqNE3r8H3/8wf+yiU6ng1Kp5O8rKyvDxx9//FByPvfccw/leUw1c+ZMi56fCBvNxiEW\nN2vWLOj1ehiNRsTGxmLDhg2QyWQmPVdRURHOnTs35VZTHRgYgI2Nzb96jqHXiRBT0Cd7MmXY2toi\nIiIC1dXV6O7uxnvvvQepVIrU1FRcunQJANDQ0IDly5dDKpXi2WefRU9PD65cuQJ/f38YjUZkZGTg\n66+/BsdxOHz4MIqKivDmm28CGPwfwKZNmxAYGIjNmzfj6tWrAIC1a9fi3XffRXh4OBYtWgSVSjVq\nvqFP1mq1GsuXL0diYiJ8fHyQnp4+6vFLly5Feno6Fi1ahJycHLS0tGDDhg0IDQ3FG2+8gRs3bgAA\nPv/8c4SEhCAoKAjvvPMObt26xedNTU2FRCLBzp07H94LTawSFXsyZXR1deHo0aNYuXIlioqK8M8/\n/6C+vh4RERFIS0sDMNgQJy0tDQ0NDaiuroaDgwP/eFtbW2RmZiIxMRFarRavvPLKiD7Hu3fvhoeH\nB3Q6Hdzc3LBnzx7+Pp1Oh8rKSmRnZ49ZWIc/1+nTp/HBBx9Aq9XiyJEjaG9vH/X4trY21NTUYPPm\nzdi2bRvS0tKg0Wjg6+uLL774AgAQFxeH2tpa/PTTT+jt7UVVVRUAYNeuXZBIJGhubkZ/f/+/eGUJ\noWJPpoC+vj5wHAe5XI7o6GgsWbIE5eXlWLt2LaZNm4aEhAScPXsWRqMRYWFh2L59Oz799FPcvn37\nnksjjLExF4lSKpVQKBQAgNTUVJSVlQEYLMrx8fFwcHBAWFgY6uvrH5g5JCQE3t7esLe3R3h4OGpq\nakY9bvXq1bCzs8Nff/2F6upqxMTEgOM45Ofn48yZMwCAy5cv49VXX4Wfnx+OHTuGiooKAMCJEyeg\nUCggEon43ISYyqxtCQkZjxkzZkCr1d6z/+6iLRKJsH79ekRGRqKkpAQBAQETXod8rB8Ejz76KIDB\n1TUHBgYe+DzD+xvY2dmN+cn7qaeeAjB4zd7FxWXUcW7duhXp6ek4dOgQsrOzodPpHpiXkImiT/Zk\nSpLJZCguLsbAwAC++eYbhIeHY/r06bh06RLmzZuHjIwMSCQS/lr+kDlz5uDatWv89vBiuWrVKhw4\ncAB37txBYWEhYmJizD6OofO7u7vD09MT3333HRhjMBqNOH/+PIDBa/NeXl64efMmSktL+cdGRUXx\neYuKisyelfy3UbEnFjf8WviQ5ORkzJo1C0FBQVCpVPx19OzsbPj7+yMkJAQSiQTh4eEjniMsLAx6\nvZ7/glYkEvH3bd26Fb/99hs4jsPVq1exZcuWUTOMludBx4znMbm5uaiqqoJUKgXHcTh79iwAIDMz\nEzKZDCtXrsSyZcv447dv347z58/Dx8cH9vb2Y56DkPGgqZeEEGIF6JM9IYRYASr2hBBiBajYE0KI\nFaBiTwghVoCKPSGEWAEq9oQQYgX+BzWwZezpSF7cAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "bias_rand_r1 = numpy.loadtxt('reads-bias-random.r1.dens')\n", "bias_rand_r3 = numpy.loadtxt('reads-bias-random.r3.dens')\n", "bias_rand_r5 = numpy.loadtxt('reads-bias-random.r5.dens')\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(bias_rand_r1[:,0], bias_rand_r1[:,3], label='r1')\n", "plot(bias_rand_r3[:,0], bias_rand_r3[:,3], label='r3')\n", "plot(bias_rand_r5[:,0], bias_rand_r5[:,3], label='r5')\n", "plot(bias_rand[:,0], bias_rand[:,3], label='r10')\n", "axis(ymin=0)\n", "legend(loc='lower left')\n", "xlabel(\"Position in read\")\n", "ylabel(\"Average graph density within given radius\")\n", "title(\"Effects of errors density for different radii\")\n", "savefig('/tmp/foo.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEVCAYAAAABwEUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8VFX6P/DPnV5TJoVACEkQSAg1EKqUgAp8ERYVRFAE\nDOi6vxV3Rd1dXRVQV5e1oAjWRZBiYe2gEmqkiKAEQq9JgISWnsn0uXN+f9zMTYYkTCaZSSHP+/U6\nr2EmM/c8Nxmec++5557DMcYYCCGE3NQkzR0AIYSQwKNkTwghbQAle0IIaQMo2RNCSBtAyZ4QQtoA\nSvaEENIGULL3o+XLlyMxMREdOnQAAGzfvh0jRoyAXq9HVlZWM0dXU0lJCWbPno2IiAg8/fTTzR1O\no0kkEmRnZwds+7t27UJiYqLfthfo78fChQvx4IMPAgAuXLgAvV4P90jr2v72CxYsQHx8PAYPHuz3\nWFqC1NRUrFixAgCwbt06jB07VvyZXq9Hbm5uM0XWRBipt9jYWKZWq5lOpxPLvHnzGGOM2Ww2ptfr\n2blz58T333bbbWzlypWNrpfjOI/t+suaNWvYmDFjmM1m8/u2m0Ogfk91iY2NZdu2bWvw5/31/ajL\nwoUL2YwZM2r92fV/+4sXL7KwsDBWVFQUsHjqkpOTwziOYzzPB7Se1NRUtmLFioDW0ZLRkb0POI7D\nxo0bYTQaxbJ06VIAQGZmJvR6PTp37gwAYIxhz549uPXWW/1SNwvAvW+7d+9GSkoKFApFo7fldDpr\nvMbzvE/b8PX9zY3juAb/XRr7/XC5XA36nNv1f/vdu3cjISEBBoPB523V9rdvCF9+l63tu9IiNGtT\n08rExcXVeiS3ZcsWplarmUQiYTqdjk2fPp3pdDrGcRzTarWsS5cujDHGiouL2ZIlS1hSUhIbN24c\nS09PF7fhcrnY999/zyZNmsSCg4NZSkoKu3jxIhs+fLi4HZ1Ox9avX8/MZjObM2cOi42NZQaDgQ0f\nPpy5XK5aYz516hSbN28ei4mJYY8//jg7ffo0Y4yxWbNmMblczhQKBdPpdLXul8PhYF988QUbNWoU\n69OnD/vvf/8rHgnu2LGDRUdHs/fee4917dqVPfjgg2zhwoVs2rRp7NFHH2VRUVFsxYoVrLi4mP37\n3/9mXbp0YZMnT2YZGRni9hcsWFDj/ceOHWN33303i4iIYO3atWPz58+v8+/xww8/sAEDBrCEhAS2\nfv16jyP7+sT+wQcfsPj4eDZ06FD2448/itvdtWsXGzduHAsNDWXR0dHs9ddfFz/XsWNHxhhjM2bM\nYBKJRDzT+89//sPuvPNO9s4773jE2KtXL/btt996vGa1WplWq63x/cjLy2P//Oc/WWxsLJs9ezbL\nzMwUPzNr1iz217/+ld17773MYDDU+ve6du0ae/rpp1m7du3YPffcw5566inxyN599Ox0Omv87T/4\n4AOmUqmYVCplOp2OLVy4kDHGWFZWFvvjH//IYmJi2Pz589n58+fFumJjY9ny5cvZkCFDmFarZTzP\ns7Nnz7Knn36aderUic2dO5cdO3ZMfP/IkSPZq6++ysaMGcOioqLYE088wUpKShhjjMXExDCO48Sz\n5V9//bXGvtX2Xdm/fz8bPHgwCwkJYYMHD2bvvPMOczgc4meOHj3K7r33Xta+fXv27LPPehzZr1y5\nkg0bNkx8b1OfFTYHSvY+iIuLY1u3bq31ZxkZGWIicLv+C3T33Xezxx9/nF25coXt3LmTdejQgZ05\nc4YxxtjXX3/NunbtyjZs2MB4nmdZWVniKfX121m2bBm7//77WVlZGXM6nWz37t11xhwbG8tefvll\nVlhYyF599VUWFxcn/mz27Nns+eefr/Ozb7/9Nhs9ejQ7evQoO3v2LEtNTWUffvghY0xIfDKZjKWl\npbHLly8zi8XCFixYwORyOVu6dCmzWCzMYrGwmTNnsqlTp7KLFy+yr776ihkMBpaTk8MYY7W+f8qU\nKWzp0qXMbrczk8lU6398xhg7cuQICwsLY99//z07d+4cGz9+vMfvyVvscrmc/elPf2LXrl1jH330\nkcffLiUlhX3zzTeM53lWWloqJt3qyZ6xmo3/+vXr2aBBg8Tnhw4dYmFhYR4JqLrr/64jRoxgjz32\nGLt27RpbsWIFCwoKYhaLhTEmJHutVss+/fRT5nA4mNVqrbG9yZMnsxkzZrDLly+zVatWMa1Wyx58\n8EHGWM2ukuv/9qtWrfJIfoWFhSw0NJR9++23rKysjL3yyits6NChHvuelJTEdu7cyaxWK3M6nSwy\nMpKtXLmSlZeXs08++cTjdzVy5EgWExPDtm7dyvLy8tiAAQPYf//7X8YYY7m5uV67cWr7rhw4cIDt\n27ePOZ1OtmfPHhYbG8u2bNnCGBMOniIjI9nrr7/OCgoK2Pz585lCoaBkT+onNjaW6XQ6FhISIhb3\nF/b6RMCY5xeovLyctW/fnpnNZvHnf/nLX9h//vMfxhhjU6dOZUuWLKm13uu/iEuXLmVjxoxhx48f\nv2G8mZmZrH379h6vRUdHi8lr9uzZ7Lnnnqvz87feeivbs2eP+Pybb75h48ePF/eX4zh24cIF8ecL\nFixgnTt3Fp87nU4WFhbGTp06Jb72wAMPsDfffLPW9zPGxCPSK1eu3HDfFi9eLCYyxhjbtm2bx+9p\n6NChN4xdKpWygoICxphwFqDT6djJkycZY4z169ePvf766+KRp5u3ZG+xWFhoaCg7e/YsY4yxJ598\nkv35z3+ucx+qx1tQUMDUajWrqKgQf37rrbeyr7/+mjEmJPvRo0fXuS2Hw8GCg4M9vifDhw+/YbKv\n/re/Pvl9+OGH7OGHHxafu5P51atXxX1/8cUXxZ9v3ryZ3XHHHR4x9e3bl+3fv58xJvSXu69vMcbY\nq6++yu67775aY6tNbd+V6/3zn/9kjz32GGOMsX379rGYmBjxZ2azmSmVyjad7KnP3gccx+G7775D\nSUmJWObMmVOvz+7evRsFBQXo0KEDQkNDERoaio8//hi7d+8GAGRkZNS7/3bOnDlITU3FhAkT0KtX\nL3GEwfX27NmDfv36ebyWkpIi1nkjJpMJv/zyC+68804x3tmzZ+OXX34R39OuXTvExMR4fG7QoEHi\nv0+cOAGbzYZu3bqJr/Xv3x+7du2q9f0AsGTJEpjNZvTs2RPjxo3Dzz//XGt8+/fvR9++fcXnycnJ\nHrHv3bv3hrG3b98e4eHhAACZTIbw8HDk5+cDAFavXo2srCx07twZ9957b71HyqhUKkydOhVr1qwB\nYwyff/65OBrGm19//RWdO3eGVqsVX6v+t+I4rsbvqroTJ07A5XKJ14wAoF+/fg2+prB161asW7dO\n/P2Fh4fDZDJh586d4nuqx7N161bs2rVLfH9oaCjOnj3r8f7qf6+oqCjx911f1+9/fn4+Hn30UfTu\n3RtBQUFYsmQJDh8+DADYt28f+vTpI75XrVb7dSRVa0TJvokMGTIEERERuHr1qthQlJeX47vvvgMA\njBo1qs4kLJFIPP7TajQaPPPMMzh37hw+/vhjzJ8/H8ePH6/xuWHDhiEzM9PjtQMHDmD48OFe49Vq\ntRg0aBDS09PFeEtLS1FSUiK+RyaTeXyG4zhIpVLxeWJiIpRKJU6dOiW+9vvvv2PEiBHi8+rvB4BO\nnTph+fLluHLlCqZOnYrp06fXejFy4MCBOHjwoPi8+n7WJ/Yb6dGjB1avXo3Lly+jV69emDt3bq3v\nk0qlNZLprFmzsG7dOmzduhUajeaGCbq6wYMHIzs7GyaTSXztt99+8/hbXf+7qi4xMRESiQTnzp0T\nXztw4AA4jqtX/dcbPXo0Zs6c6XFgU1FRgSlTpojvqf73Hz16NFJTUz3ebzQa8eSTT3qty71fN2qY\nrv9uAcDLL78Mh8OBH3/8EWVlZXjiiSfE78rAgQM9GmmLxYKTJ0/Wb+dvUpTsfdTQI6WQkBAMGzYM\nzz77LM6fPw+e53H06FH8/vvvAIBp06bhgw8+wE8//QSn04nDhw+juLgYgHA07H4fAPzwww84e/Ys\nXC4XtFotFAoFVCpVjTr79u0LhUKBV199FYWFhfjPf/4DmUwmHmF525cHH3wQL7zwAjIzM+FyuZCf\nn4/NmzfX+f7rtyeTyXDnnXdiwYIFyM/Px7fffotNmzbhrrvuqnMba9euRUFBARhj0Gq10Ol0tb5v\n/Pjx2LRpE3744QdkZ2fjrbfealTsbg6HA+vWrUNZWRkAoWHV6/W1vrd///44cOCAx2tDhgwBx3F4\n6qmnMHPmTK/1uYWHh2PAgAF49tlnce3aNaxatQrHjh0Tx4J7+1vJ5XLcfvvtWLRoEa5cuYK1a9fi\n0KFDdb7f2/buu+8+fP311/j2229hMplgMpnwww8/oKKiotb333777Thy5AhWr16NkpISWK1WZGRk\neBy911Vnx44dERkZ6fEdr0+8ly5dgsFgQFhYGDIyMrB69WrxZwMGDIDNZsOSJUtQUFCA559/PiAj\n2loTSvY+mjhxIvR6vVgmT54s/uz6o6jrn7///vuIjY3FlClTEBERgUceeQTl5eUAgD/84Q9YvHgx\nli1bhrCwMDz88MOwWq0AgKeeegqvv/46QkND8b///Q9nzpzBHXfcgeDgYDz88MN4+eWXPU7fq9u0\naRPy8/ORnJyMixcvYtOmTR7x3ejI7+GHH0ZaWhpeeOEFGAwG3HHHHTh9+vQN9/f6195880306dMH\nI0eOxOrVq/G///0PcXFxdb4/PT0dPXv2RLt27bB27Vp8+OGHkEhqfk179uyJlStXYtGiRRg/fjxm\nzZrlsS1fY69u7dq1iI+PR6dOnbBv3z6PhqT65x599FFs3LgRBoMBb775pvj6zJkzceTIEcyYMaPO\nOmqLYd26ddBoNBgwYAAyMjKwbds2qNXqOn9X13v33XcRGRmJvn374ptvvsGf/vSnOuu7fnvXPw8J\nCUF6ejp27NiBbt26oWvXrli9enWdMUilUmRkZODUqVPo378/OnXqhDfeeMMjwdZVH8dxeP755zFn\nzhyEhoZi//79tf6urq974cKFOHToEDp27IjXXnsNjz32mPgeiUSCrVu3Ys+ePejTpw+USqVHN2lt\n+3+z41hbb+4I8bPVq1fjv//9r0d/NSHNjY7sCfEjs9mMJUuW3BTTT5CbCyV7QvwkPT0dHTp0QL9+\n/TBx4sTmDocQD9SNQwghbQAd2RNCSBsg8/6WwGkLV8AJISQQfO2UafYjeyZM2dAqy4IFC5o9Boq/\n+eOg+Ftfac2xM9awnvdmT/aEEEICj5I9IYS0AZTsGyE1NbW5Q2gUir95UfzNpzXH3lDNOvSyMSv9\nEEJIW9WQ3On1yL6iokJcAuzq1ase08QSQghpHbwe2ffr1w+7d++G0+lE7969kZiYiMTExBqzDDao\ncjqyJ4QQnwXkyN7lckGj0WD16tVIS0vDpk2b6OieEEJaGa83VYWFhWHbtm345JNP8MUXXwAQFgIg\nhBDSenhN9m+88QbeeustzJ07F507d8a5c+cwatSopoitTWCMwcEY7Iyh1OlEqdOJEocDJU4nKnge\n7hM19yPPGOwul/gZh8sFq8sFS7Vi5nmUOJ0odjhQVPlYzvOQApByHGSVRSeVIkapRIxKJTwqlbhF\nrUZ3jQYxSiXd4UzITaRNjMaxuVxCEnU6Ya682MxV1s8BkHEcVBIJlBKJ8MhxqOB5FFYmy0KHA8UO\nB2wuF+yMiY/OytjdKZED4GAMFTzvtTgqE7Wrsn45xyFEJkOoTIZQuRwhMhl0Uikk1bfPcZAAUEgk\nUFR+Rl4Zs0YigVoigVoqhVoiQahMBoNcDoNMhjC5HHqpFC4IjYWTMfCMoZzncdFqxQWbDRdtNlyw\nWnHWYsEJsxkVPI8EjQZJGg2GBAVhdGgouqnV1AAQ0gI0JHd6Tfa1HcVzHIft27f7Fl1tlXMc0ouK\ncMVux1W7HVfsdhQ6HGJCNLlcqOB52Fwu8IyJyYpnDAqJBFqJBFqpFNrKBGdzuWByuWDieZgqt1Hq\ndMLBGEJlMoTIZNC617tE5VQNAJyVCdzqcsHGGKwuF3RSKcJkMoTL5QiTy2GQy6GuTLLuZCurlvjc\nv0Qpx0EvlUJ3g6KVSqGUSCDjOEjRMucIKnU6cdJsxjGTCbvKyrC9pAQ8YxgVGoqRwcHoqdWim0aD\nMLm8uUMlpM0JSLKvvi5kcXEx1qxZg86dO2PRokUNi7J65RyH2w8dQjuFAlEKBdrJ5QiXy6GvPKrV\nSiTQVSZGKcdBWnlkK+U42K9L7BaXC8rrGgCtVIpQmQwaiaRFJtTWhDGGbKsV20tKsKusDKfMZpyy\nWCAF0E2jQaJGg55aLXpVliiFgn7nhARIQJL99RwOBwYOHIiDBw/6VFGtldPQy1aNMYYChwOnLRac\nMJlw1GTCkcrCAPTWapGs0yFZr0eyTocEjcbjbIgQ0jANyZ1eL9AWFxeL/7bZbMjIyEBQUJDv0ZGb\nDsdxiFQoEKlQYFhwsPg6YwxXHQ5kVVTgYEUFNhQW4sXcXOTbbOiv12NkSAhGBAdjSHAwdJXdaoSQ\nwPJ6ZB8XFyeejqtUKgwZMgTz5s1DcnKy143HxcUhKCgIUqkUcrm8xqrxdGTftpQ7nfi1vBw7y8rw\nc2kpDlZUoKdWi+HBwRhWWegaACHeNUk3ji/i4+Nx4MABGAyG2iunZN+mWXge+4xG7Cotxe6yMuwt\nL0eMUolhwcEYGBSEgXo9krRaSKnrhxAPfu3G2bZtG2677TZ89dVXtV5ou+eee+pVASVzUhe1VIrU\nkBCkhoQAEEZFZVVUYE9ZGXaUluI/Fy7gkt2OZJ0Og4OCMDIkBMOCgxEsa9YF1ghpler8X7Nz507c\ndttt2LBhQ4OTPcdxGD16NOLj45GWloY//OEPjYuW3NRkHIf+ej366/XiayUOBw5UVOCXsjK8efEi\n7jt+HIkaDVJDQjCqsu9fT8mfEK8C2o1z+fJltG/fHidOnMDEiROxe/duREVFVVXOcViwYIH4PDU1\ntU3OM03qz+ZyYX95OX6uHPu/32hEsk6H20NDcVtoKFL0eqgktEwDublkZGQgIyNDfL5o0SL/9dm/\n8cYbVW+q5ch+/vz5PlU0f/58dO/eHQ8//LDHdqmbhzSGmeexu6wMW0tKsK2kBCfMZiRoNMIZgk6H\n/no9emi14s10hNwM/NpnbzQawXEcLly4gPT0dNx+++0AhL78sWPHet2w2WwGz/PQ6/UoKChAeno6\nnnjiCZ+CI8QbjVSKMQYDxlQOArC6XDhcUYEDRiMOVFTgo8uXcdJsRnuFAj21WvTUatFdq0VnlQpx\nKhWiFApI6AIwaQO8duMMGzYMn376KTp16gQAuHjxIqZPn47du3ffcMM5OTm4++67AQgzZz7wwANI\nS0vzrJyO7EkTcDKGbIsFRytv/DphNiPHakWu1YpSpxOdlEp0VavRS6dDb60WvXU6JKjVkFN3EGmh\nAnJTldlshkKhEJ8rFAqYzWavG46Pj8ehQ4d8CoaQQJBxHLppNOim0eCeiAiPn5l5HhdsNpwym3HE\nZMK3hYV48fx5nLda0a1aA9BLq0UvnQ7RNA0EaaW8Jvunn34aI0aMwLhx48AYw+bNm/0yLw4hLYFG\nKkVi5dw+k8LDxdfNPI8TZjMOV1TgiMmEzSUlOGoyoczpRIxSiU4qFWIrH2OUSkQrlehYWYKkUmoQ\nSItTr9E4hYWFSE9PB8dxGDt2LMLCwvxTOXXjkFbGxPO4aLPhvNUqlny7HXk2G/JtNuTZbJAA6KJW\no5tGg65qNbqq1bhFrUaMUon2SiXND0QaLaB30JaWlqKsrEw8YnH34TcGJXtyMyp2OHDGYhGK2YzT\nFgtyrFZctFpR4HCgvUKBTioVums06KvToa9Oh946Hc0TROotIMn+559/xsKFC/H7779DqVSiuLgY\nSUlJOHr0aKOCBSjZk7bH7nIh32bDeZsNx00mHKyowKGKChwzmRCtVKKLWo34ypFCcSoVOleuHEZD\nR0l1AUn248ePx/vvv4+JEyfi4MGD+Pzzz7F79268++67jQoWoGRPiJuTMZw2m5FdOUoo12pFjsWC\nc1YrTpvNiFIo0EOrRY/K4aM9tVp012igpBFDbVJAkn2/fv2QmZmJoUOHIj09HXq9Ht27d8eJEyca\nFSxAyZ6Q+nAPHT1mMuGY2SysG1BRgWyrFfEqFXppteir06Ff5VQT4TRz6E0vIEMvQ0NDYTQaMX78\neEyZMgXR0dHo3r17g4MkhPim+tDRu6u9bnO5cLJyyOhBoxH/vnABmUYjQmQy9Ku8c7hb5YXiBLUa\nodQItGlej+xNJhNUKhWkUikyMjKQn5+Pu+66C1qttvGV05E9IX7lqlw+MtNoxAmzGacrl488bTZD\nJZGgt04nrB5WWbppNDSFdCvk924cp9OJsWPHYtu2bY0OrtbKKdkT0iQYY7hkt4urh7kvDOfbbLhF\nrRbPALqp1UjQaNBDq0UIzSbaYvm9G0cmk4HjOOTm5iIuLq4xsRFCmhHHcYiuvPlrfLX7ZCp4Hmer\nDRHNKC3F+5cu4bjZjGCpVLwoHCmXQymRQCmRQMFxUEkkMMjliJTLEaFQIFIuh4ZGDLVo9eqz79ev\nH0aPHo327dsDEL44S5cuDXhwhJDA0kml4lj/6lyM4YLNJlwUNplQ4nTC5nDA5nLBzhgsLheKHQ5c\nczhQYLfjqsMBDoCC4yCXSCDnOMg5DkqJBFqpFJrKR61UCrVEUlWq/UxXWfRSKVQSCYqdTly128VS\n5HSigudh4nmYeR4mlwtOxqCorMf9GKlQoItaLZbOKhXUlaOW3MfCHABlZQwqiaRNzIPktc9+1apV\nVW+uPHXgOA6zZs1qfOXUjUPITYExBrPLBbvLBQdjYrG5XLC4XB4J2szzsFS+buF5mCt/XlGtmF0u\nGGQytFMohCKXI1wuh66ywXA3IDKOg50x2F0u2Crru2K346zFIpZsiwX2anmGA+CCcIHbWhkHIDR8\n7eRyRCkUYolUKBAulyOisv6IyjMZg0zWrLOltrg1aL1WTsmeENICOFwuGHkeV+12XKk8U7lss6HA\n4UCBw4HCyscCux2FDgfKeR4GmQyRCgUi5HLoKxsgd2MUIpOhs0qFrpVTZoT5eSQUJXtCCGkCDpcL\nhe5uLIdD7F5yP5Y4nThnseB05fUQKcchSatFv8oFdfrpdOiu1TZ4niRK9oQQ0sIwxlDgcOCYyYQD\n7oV1jEbk22zoptHgFrUat1ROjXGLWo1Ber3XdZUp2RNCSCtR7nTiVOUUGecsFvHx7S5d0Pu6C+bX\nC0iyz8nJwZo1a7B3715YrVaxou3bt/tUUa2VU7InhBCfBSTZT5gwAUOGDMHo0aMhr7zIwHEc+vfv\n3/BI3ZVTsieEEJ8FJNn37t0bhw8fblRgdVZOyZ4QQnwWkGT/73//GxaLBTNnzkRoaKj4usFgaFiU\n1SunZE8IIT4LSLKPi4urdT3NnJwc36KrrXJK9oQQ4jMajUMIIW1AQ3Kn1wkhbDYbvvjiC/z5z38G\nAJw5cwYbN25sWISEEEKahdcj+3/84x9gjGHjxo04duwYTCYThg4diqysrMZXTkf2hBDis4Ac2e/Y\nsQOLFy+GQqEAAGi1WkrQhBDSynhN9gkJCSgrKxOf//rrr0hOTg5oUIQQQvzL63z28+bNw1133YW8\nvDyMGjUKV69exZo1a5oiNkIIIX5S79E4Bw4cAM/zGDhwoP8qpz57Qgjxmd+XJQSAiRMnYvr06Zg0\naZJfFhknhBDS9Lz22T/55JPYtWsXkpKSMHnyZHz55ZfihGiEEEJah3p34zidTuzYsQMfffQRNm3a\nhPLy8sZXTt04hBDis4B04wCAxWLB999/j/Xr1yMzM9Mv688SQghpOl6P7KdOnYp9+/Zh3LhxmDZt\nGkaMGAGpVFrvCnieR0pKCjp27IgNGzZ4Vk5H9oQQ4rOAHNnPmTMHn332mU8Jvrq3334bSUlJMBqN\nDfo8IYSQxqsz2W/btg233XYbKioq8N1334mvM8bAcRzuuecerxvPy8vDjz/+iH/+85948803/RMx\nIYQQn9WZ7Hfu3InbbrsNGzZsqHWK4/ok+yeeeAKvvfbaDS/mLly4UPx3amoqUlNTvW6XEELakoyM\nDGRkZDRqGwGb4njjxo346aefsHz5cmRkZOCNN96gPntCCPGDgPTZv/HGGzWO7OPj43HHHXdAd4MV\n0H/55Rd8//33+PHHH2G1WlFeXo6ZM2di9erVPgVICCGk8bwe2aelpWHLli24/fbbwRjD9u3bMWDA\nABw/fhzPPfccHnjgAa+V/Pzzz3j99dfpyJ4QQvwgIEf2Z86cwd69e9GxY0cAQH5+PqZNm4aff/4Z\nU6ZMqVeydwdHCCGkeXhN9iaTSZzLHgAUCgUqKioQGRlZ77toR44ciZEjRzY8SkIIIY3iNdk/+eST\nGDlyJMaMGQMA2LJlC5555hmYTCb06NEj4AESQghpvHqNxrl27Ro2b94MjuMwZswYRERE+Kdy6rMn\nhBCfNSR3BmzoZb0qp2RPCCE+C8gatIQQQlo/SvaEENIGULInhJA2wOtonK1bt+KVV15BZmYmeJ4H\nIPQX+WPxEnIdxgC7HTCbAYtF+DfHARKJUDgOkMkAudyzMAY4HDWL3V5VHA7A5RIKY8KjTAbo9UBQ\nkPCo0wn1EEJuOl4v0KakpODtt9/GkCFDIPFzIqhxkcFmA65eBQoLgaIi4bGwECgrE5Kf1Vr1yHGA\nUgkoFMKjXC68z/2ZggKgtFRIbEJlwqNUCqhUgFotPNZWlEohOVosVYnXahWSo0JRVeRyz0QskQA8\nX/U5d7FahYRrs1UlX6ezKvm6XMJzq1WIT6MR4lMohPjdydn9vuuTOsfVbACqx6lQCLFLpVVxcpyw\nLaNRKOXlQqx6PWAwAGFhVY+33AIkJQE9egDdugm/I0JIswnIaJyhQ4di+/btUAXgPzjHcWDjxgGX\nLgmlrAyIjAQiIoQkEx4ulKAgIfm5i1IpbMBm80ygwcHCZ92fCw0VElv1XXQnVXdxJ/LqxWYTkqNG\nU5V4VSpr11USAAAgAElEQVQhkVc/WrbbPRMxY0J97s9V/+z1DZM7+bobC6lUeK+v6wYwVtWQNZbL\nJfwNiouFxra4WGg0z54Fjh8Hjh0DcnKAjh2B9u2Fv1W7dsJjVBQQG1tVNBr/xEQIqSEgyf7vf/87\nsrKycO+99yIkJESsqD5THHutnOPANm4EoqOBDh2EBE3dCC2b3Q5kZwtnYFevAteuCeXSJeD8eaFc\nuCA00PHxwlnBLbcAXboIjz16CI0wIaTBApLsZ8+eLW68upUrV/oWXW2V0zj7m5PLJTQE2dnAuXNV\n5exZ4ewgPBxITgb69RMe+/QRGnyaP4mQeqGbqkjL53IJST8zEzh4UHjMyhJe791bSPy9ewOJiUBC\ngnDdgBDiwa/JfvHixfj73/+OefPm1VrR0qVLGxbldduhZE/AGHDlipD0s7KAI0eAU6eEolAIST8x\nseoicY8ewnUDOhMgbZRfpzhOSkoCAPTv37/Ghmm6YuJXHCdc8G3fHhg3rup1xoTuoFOngJMnhYvE\nP/0kPFZUVCX/nj2rGoEOHagRIKQW1I1DWqfi4qoRQkePCuXYMWE0Ubt2VY1H+/ZVI4bco4fatRNG\nD+n11DCQVikgffY5OTlYs2YN9u7dC6vVKla0ffv2hkfqrpySPfE3q1XoErp8uaq4Rwy5Rw9duSIU\nxjwbhY4dhQvF0dHCv2NihDMFuby594oQDwFJ9hMmTMCQIUMwevRoyCu/9BzHoX///g2P1F05JXvS\nnIzGqobh0iUgPx/Iy/MsV68KZwKdOlUl/+oNRFSUcG+HwSDcO0FIEwhIsu/duzcOHz7cqMDqrJyS\nPWnpHA6hIbh4Ubh/4NIlz7OGK1eEO7ZLSoRuIffNgBERVTcIuruPoqKqGomQEOpCIg3m12RfXFwM\nxhg++ugjWCwWzJw5E6HVboYx+GFIHCV7ctNwuYTpOYqKhLuO3aV6F1L17iWbTUj67hsK3Y+hoUJD\nEBwsFINBeJ3uSCbV+DXZx8XF3XDUTU5Ojm/R1VY5JXvSVpnNQtLPz6/qQrp0SWgwSkuFC83uqSvy\n8wGtVriO4C6dOgnTUrgfo6OpG6kNCUg3jtVqrTEvTm2vNQQle0LqgTHhjMF9HcHdpeSenuL8eeEs\non17IC6uqnTqVHWhuWNHoZuJ3BQCkuz79euHzMxMr681BCV7QvzEbhcagtxcoeTkVDUM7ke5HOjc\n2XOuos6dqxoFms201fDrTVWXL1/GpUuXYDabkZmZCcYYOI7DtWvXoHTPOkkIaRkUCiFxd+5c+88Z\nE7qE3PMVnT0L7NkDrFkjNAT5+cI1gpgY4azA3Si4t9munbDeAV1UbrXqPLJftWoVPvnkE/z+++9I\nSUkRX4+NjcWMGTNw++23N75yOrInpGVwT1538aJwZpCdXdUwZGcL3URWq9AV5L54fH1x36TmXoMB\nqH3BHUCYLpznhSnHeb7q/e5Hl0u4S9q93oLRKDx3T2vuLjxfNV24VCr8W632XJRHrxfqsVg8y/VT\nm1dfLIjjqkr1BX8YExq+Ll2Arl2Fcsstwkir8HDh4noTzNwbkG6cr776CpMnT25UYHVWTsmekNbD\n4RAWuXFfPK5eysuFhOxWPVFev+AOUJWY3Um6eoJ0f1ar9UzYWm3V2hDuIpV6NhrupO5ekMfdULgb\ngetL9UWL3Ku+VU/s7jUq3A0AIIyqOnNGODs6c0ZoEK9dE66rGI3CiKrISKEhSEioKt27+216b78m\n+zVr1uDBBx/EG2+84TEqx92dM3/+/MZFC0r2hJCbjNMpdJddvgycPi0U96R+J04IjVavXsJ8Tr16\nCQ1AQoJwZuQDv/bZm81mAIDRaKSJzwghpD5ksqp5mPr08fwZY8IoqiNHhJKeDrz9ttAg6HRC0u/W\nDfjHP4SuIT/z2o1jsVigVqv9XjFAR/aEEALGhHss3GcAd90lDKO9gYD02Xfp0gWRkZEYMWIEhg8f\njmHDhiHYx1OOOiu/LmDGGEwOEywOC2y8DTanDTbeBqfLCRdzwcVcYIzBxVzQK/WI0EQgVB0KCed5\nQYQxBqPdiHJbOSwOC6xOKyxO4dHBOyCVSCHlpJBKpJBJZFBKldDINdAqtNDKtdDINWBgsPN2sTh4\nBziOg4STiEUmkUEtU0MlU93w7IcxBovTApPdBLPDDJPDBJPdJD5W2Ctgdpg96rPzdjhdTsgkMrFI\nJVKoZCroFDroFXroFDroFDrIJDIwMDDGxEee8eBdPHjGw8Vc4F08nC6n+LrT5YSEk0Aj14j77v63\nWqaGRq7xul+EkOYRsJWqzp8/j927d2P37t348ccfERoaikOHDjU4ULFyjsPru97GzhPHkXXpOC45\njsPJmQGHBuCVgEMF5lSCYzLIZRKhyDnIZRyYwggLVwArM8KgMiBcGwaTzYpicynMfDkkvAawBYHj\n1ZC6VJAyNaRMBSkng1TGQyLjIZHy4GROuDgbHDDDzkzCI0zgIIEUCkihgIQpwDGhx4vBBQYXXODB\n4AQvscDF2SGDCkqJGlJODp5VJlXmhIvxcMIKKRSQMy1kTAuZSwupSwsJr4XEqQXn1AJ2DVwOFXi7\nArxNAaddAZdTComUh0TmhETmBCdzQqayQqGtgFRdASiNYPIKIXk7OTidgNPJgXdykHBSyCRSyKVS\nyKQSyKRSSDkZJJzwKOWkAOeCkzPDyZnF/XYwi/BvlwUOZoNCokJMUCf0aJeIhPAEJIQloKuhK6J0\nUQjXhCNEFUINAiFNLCDJPi8vDzt37sTOnTtx6NAhGAwGDB8+HM8880yjggWEgGWT/oQO8iQMiEvC\nnQN6YGRKJLRaDgoFxGKxCKPC3PNOXb4sdH3l5ADZuQ5kXy2ETVIErUKD5O4hGNA7CCn9ZEhKEs6Q\n7HZhlJbdLtyl7r5AX14uFIej6oK7+6J79UEC1R+rDyBgrHI6lBIeRaVWFJZZYHM6oNfIoNPKoNdK\nodNIoVWqIJdJPbbpHlTgHgigVAqDDdRqYRoUjUYYHOB0CnE7HMKj0Vi11rd7xl6l0nPeLYNB2F/3\n/rkHJbgHLFQf9VZ9kIT792SxCL8ns8UFo9WM7OJcsLCTMHQ9BXn7U7Drz8CMazCxQjhggYYLQ5A0\nEu21HRFviEFC+464JaIjovXR6KDvgA76DtQoEOJHAUn2EokEAwYMwDPPPINJkybV+z+s1WrFyJEj\nYbPZoFKpcN999+GJJ56oEbDRyKDT+RRzrSoqhGRJ+cT/GBOGWZ88KZRz54TGgTHAwWywoBDF9qu4\nUJKPS6Y8FNnzgOA8yAx5YNrLcKgugUltUDs7IBidEKWMQ6fgOHSLiEfPjnHoGhmDWyI7wBCshFTa\n3HtLSMsXkGSflZWFXbt2YdeuXbhw4QK6du2KESNGYO7cuV43bjabodFoYLPZ0L9/f3z77bfo0qVL\nowImLR9jwoy/164BJpPQEBeVm3Ch9BLOFlzAuaJc5FXkosCRg3JJLuyqPPDqy4AtGJyxI+TmTgh2\nJiBK2h3xuiR0j0hEXPsgcZEpd/HHQQIhrVHA+uyNRiP27NmDnTt3Yu3atQCACxcu1LuSoqIi3Hrr\nrdiyZQtiYmIaFTC5ObmYC1crriG7MB8nL5/HwYsncezacWQbT+Cy/SQUrlBoTN0hKUmE83IiKnK7\nQ23pit5x0ejVUyIuQatW17wxsvpNmEaj0Bh17y4Mc05KotmDSesTkGSfkpICq9WKoUOHiiNyYmNj\n67Vxl8uF5ORkHDt2DG+99RYee+yxRgdM2h4Xc+F86XmcKjqFEwUncLLoJE4UnMDpwrMothQjhIuF\n0hwPx7V4cDYDFJwaSokaSqkaKqkaKqUUapUEGrUEapUEHJPj2vlQXDgVhgsnDegYFobOMZoaN2cG\nBQlrkVQvHTrQBJKk+QUk2V+7dg2RkZGNCiw3Nxfjx4/HunXrkJycXFU5x2HBggXi89TUVKSmpjaq\nLtK2WBwW5JbmIqc0B9kl2SizlsHitMDsMMPitMDisIjDdt3FzttRYi1BkbkIxZZiFJmLIYEMkfJ4\nRMg7I0wSj1B0htbeGdKyLuCL4lBapEBhoTAcOi9PuFDfsaMwHFqp9Jz6RSKpeeFbo/Gcjj46Wrig\nbjAIRaFo7t8kackyMjKQkZEhPl+0aFFgunH84amnnkKXLl3w6KOPVlVOR/akhSi1liKnRGgwskuy\nkV2ajXPF53Cu5BzyyvPQQd8Bt4Tegna6dghTh0EnCQdnDQMzGSBz6SBjWkhdGkh5LeTQQqvQQKtU\nQ6fUQK2Qw2TyXOL24kVhNcPiYqEolULS1+uFaxHux6CgqgbBfXZhMAh31wcFVc1BplY3yfxbpIUI\nWJ99QxQWFkImkyEkJARFRUUYNWoU0tPT0b7anWGU7Elr4OAdOF92HueKz+Ga6RqKLEVCqTwzuP4m\nOfeNge4b6CScBDqFDiGqEISoQhCsDEaIKgTtdO0qh6dGwyCLhtYVDTUfBaktDGaTBBUVwhxjJSXC\nHFtFRVWNw/XzkFmtVUOGZTKhhIcLZxDuFQ87dKhqLEJDqx6DgoSGhRqL1qNFJfsjR45g1qxZ4Hke\nUVFReOCBBzBz5kzPyinZk5scYwwOlwNGmxFltjKUWctQai1FibUEVyuuIt+YL5TyfOSV5+Gq6SrK\nbeUI14SjnbYd2unaIVIbiUhtJCI0EeKjQW1AmCYMBrUBBrUBMolMnGCS54Wuo8LCmqseuhuLkpKq\nfxuNwn0V7kkmQ0KqzibcJTKyar30qChhNJREUnWvhvvRfY+G+9HlqjpTcReVynMm4erbcRens6rx\nct//4nIJDZt75cbSUmG0l8tVVXhe2L77rCcoSCgajecElxKJcJ9Kfn5VKS72bCylUiGO6jMsV1QI\nDWT1qf5jYoRY8vKqtlVSIjSs1UePBQd73uficgnbacj1n4Ak+/79+yMtLQ3333+/x4Lj/kDJnpCa\n7LwdBaYCXDVdxdWKqygwF+Ca6ZpYCswFKLYUi6XEUgKNXIMQVQhC1aHCoyoUYZowRGgiEK4JR4Qm\nAhHaCETroxETHINQVajHPTPu6ePLy4VE5W4M3GcUBQVVNzVeuSIkSsaqEqN7ynr3TYHuR46rOS29\nzVZzJmGp1PO6h0xW9XN3Euc4oSGqXjQazxsgJRLhLMd9M2H1M5/qxekUGjD3mU90tJCcqzeYTqcQ\nh7uR0umEUlRUNd1/drZwg2doqLAN9/WY0FDhfdVvgCwvr4pVKhX2Jz8feOQR4PHHvU6H4yEgyf7M\nmTNYuXIl1q9fj5SUFDz00EMYM2aMX+6GpGRPSOO5mAtGm1E8Yyi1lqLEUoIiSxEKTAUoMBeg0FyI\na6ZryCvPQ155HhwuB2KCYhATHIO4kDjEBccJj5Wlvb59jTmniP/l5ABvvAGsWwdMmQI89ZQw+aU3\nAe3Gcblc2LhxI/70pz9BIpHg0Ucfxfz58xs1IyYle0KaR7mtHHnlebhQdgHnS88jtywXuaVCyS7J\nRrmtHLHBsegc2hnxofEI14SLkwRq5VpoFVoopArIJDLIJXJxoj73JH42p02czM894aCUk4rvk3AS\ncTJCCSeBQqqAUqqEUqaESqaCXCKHyWFCmbUM5bZyj2K0G4ViM8LG2zy2I+Wk0Cq0CFFWneWEqEIg\n4STCRICVkwM6XU6xa829XZPD5PEe94SB109QqJarEaoKRag6FKGqUAQrg8EzHjbeJu67VCJF33Z9\nMTB6ILoYung9OC4oAJYvB959F/jxR6Da4oC1Cliyz8rKwsqVK/HTTz9h7NixuP/++7F9+3b89NNP\n2LVrl08VNjZgQkjgmewmMfHnlObUuBDtnqXV6XLCwTvgcDnAu3ghacuUUEqVYmPgTq7u5OliLjGZ\nuv/tTpLu2W7tvB06hQ5ByiCx6JV64VFR9aiUKT22w7t4mBwmj7OcUmspGGM1Zrutvu0gRRC0Cq3Y\nKLnfK5PIoJQJ+6KQKiCXyGF2mMXtl1hKUGYrg1QiFRqryv2283ZkXsnE/vz9MNqMGBA9ACkdUpAc\nlYzkqGTEh8bXeuZkMlV1f91IwPrsg4ODMXfuXEyePNljsfG7774b33zzjU8VelROyZ4QcpO7UnEF\nv+X/ht8v/45DVw7h4OWDKLWWok9UH/SN6is2AEkRSVDKlN43iAAl++zsbHS+bsX6nJwcxMfH+1RR\nrZVTsieEtEFF5iIcunJIKFeFBiC7JBvdwrphzd1r0Ktdrxt+PiDJvl+/fsjMzPR4LTk5GQcPHvSp\nolorp2RPCCEAhLvBjxUcQ0JYAvTKG4/H9OsatCdOnMDx48dRWlqKr7/+WtxwQUEB9DQxCCGE+JVa\nrkZKBy9XZhuhzmR/6tQpbNiwAWVlZdiwYYP4emxsLJYtWxawgAghhPif126cvXv3YsiQIYGpnLpx\nCCHEZ37ts1+8eDH+/ve/Y968ebVWtHTp0oZFed12KNkTQohv/Npnn5SUBEAYeum+IcC9cVpLlBBC\nWhefJ0KzWq1QqVT+qZyO7AkhxGcNyZ1eJ7+4//77UV5eDp7nMWjQIHTr1g0ff/xxg4MkhBDS9Lwm\n+2PHjiEoKAjffPMN+vfvj9OnT2PFihVNERshhBA/8ZrsNRoNzGYz1qxZgxkzZkClUsFoNDZFbIQQ\nQvzEa7KfN28e+vXrB71ej6FDhyI3NxfBwcFNERshhBA/8fkCLWMMPM9DJqtzIE/9K6cLtIQQ4jO/\nDr10czgc2Lt3L/bu3Qur1SpW9MILLzQsSkIIIU3Oa7KfN28ecnNzMXLkSOh0uqaIiRBCiJ957cZJ\nSkrC0aNHIQnA0vPUjUMIIb4LyDj7UaNGYceOHQ0OihBCSPOr15H9yZMnER0djZCQEOFDHIfDhw83\nvnI6sieEEJ8FZPGS3NzcWl+Pi4vzqaJaK6dkTwghPgtIN05cXByUSiX27NmDuLg4aLVaStCEENLK\neE32H374IaZPn45FixYBAOx2O2bMmBHwwAghhPiP12S/Zs0abN68GVqtFgAQHR1N0yUQQkgr4zXZ\nBwcHewy7vHDhAjp27BjQoAghhPiX12Q/a9YsPPDAAygtLcWiRYswYcIEzJ07tyliI4QQ4if1mhsn\nNzcXX331FVwuF6ZNm4aYmBj/VE6jcQghxGcBGXoZSJTsCSHEd36dCE2n09W51izHcSgvL/ctOkII\nIc2mzmRfUVEBAHjllVdgtVqRlpYGAFi1ahWUSmXTREcIIcQvvHbjJCYm4sSJE+JRvsvlEqdQuJGL\nFy9i5syZuHbtGiIiIvDII4/g/vvv96ycunEIIcRnAemznzNnDhITE5GWlgbGGD755BMcP37c6zq0\nV65cwZUrV9C3b18UFhZi4MCByMrKgl6vb1TAhBDS1gVkuoSXXnoJFy9exNChQ3HrrbfiwoULePnl\nl71uOCoqCn379gUAhIeHo0ePHvj999/rFZTBYADHcW2uGAyGev1+CCHEV00yGufs2bMYM2YMjhw5\nIt6JC9TdOrXVI/62ut+EEN8EZFnCxjIajbjvvvuwZMkSj0TvtnDhQvHfqampSE1NDXRIhBDSqmRk\nZCAjI6NR2wjokb3D4cCdd96J8ePH469//WvNyunI3kNb3W9CiG9a1E1VjDHMmjUL4eHhePPNN2uv\nnJK9h7a634QQ3wQk2TscDuzduxd79+6F1WoVK3rhhRduuOHdu3djxIgR6N27tzhs89VXX8W4ceO8\nBtxWk15b3W9CiG8C0mc/b9485ObmYuTIkdDpdPXe8LBhw+ByuXwKprVZtmwZVq1ahaNHj2L69OlY\nuXJlc4dECCG18prsd+7ciaNHj3pMc0wAp9OJ6OhoPP/880hPT4fFYmnukAghpE5eM/ioUaOwY8eO\npoilxYuLi8O7776LoUOHIiQkBJMmTcKkSZMQFhbW3KERQsgN1Xlk36tXLwDC9AjvvfceoqOjERIS\nAkDoLzp8+HDTRNiCcByH5cuX4/3338fAgQPFsx3qZyeEtHR1JvsNGzbUeI0uIALTpk3D8OHDPV6r\na3ZQQghpKepM9nFxceK/CwsLkZ6eDo7jMHbs2BbRbeGv/Opr2zVo0KBattG2G0BCSMvntc9+3bp1\nGDJkCPbu3YtffvkFQ4YMwbp165oithtizD/FVzJZzfaRjuwJIS2d19E4r732Gnbt2oWoqCgAwNWr\nVzF27Fg88MADAQ+upeN5Hg6HA06nEzzPw2azQSaTQSqVNndohBDiweuRvcFg8BhWaLFYaHbGSi+9\n9BI0Gg0WL16MtWvXQq1W41//+ldzh0UIITV4vYN26tSp2Lp1K4YPHw7GGHbv3o077rgDkZGR4DgO\nS5cubXjldAeth7a634QQ3wRkuoRVq1bVWRHHcZg1a5ZPFda2nfq+frNrq/tNCPFNi5oIrV6VU7L3\n0Fb3mxDim4DMjZObm4sPPvgA6enpKCkpESvKzs5uWJSEEEKanNcLtAsWLEBycjKcTie++eYbjB8/\nHo888khTxEYIIcRPvHbjJCcn4+DBg+jTp4+4hmxKSgqysrIaXzl143hoq/tNCPFNQLpx1Go1eJ7H\nyJEj8corryA+Pt6nqY4JIYQ0P69H9r/99hsSExNhsVjw7rvvIj8/H/PmzUPv3r0bXzkd2Xtoq/tN\nCPGN30fj8DyPf/zjH3jttdcaHVytlVOy99BW95sQ4puG5IobXqCVSqXYuXMnjEZjowIjhBDSvLz2\n2d96662YOHEipkyZgvbt2wMQWpV77rkn4MG1dDNmzMC2bdtgtVrRt29fTJs2DX/84x+bOyxCCKnB\na5/97NmzhTdeN7OjP9Zbbc3dOE6nE6dOncItt9wCpVKJ/fv3Y+zYsdizZw969OjRoG22hv0mhDS/\ngIzGqW26hLYqLi4Of/vb37B27VocPnwY5eXlkEgk4HkeHMdBJpNBpVI1d5iEEFKD1yP7efPmiUf1\n7vlw4uPjMX78eCQkJDSu8lZ2ZB8fHw+NRiMuS6hUKvH//t//w8qVK8HzPH744QfccccdDd5+S91v\nQkjL4vcLtIBwkXbPnj0ICwtDWFgYfvnlFxw+fBhz585t1IyXrZV7WUKlUgkAePfdd1FcXIxPPvkE\nU6dOxfHjx5s5QkIIqcnrkf3AgQOxadMmcQ774uJijBs3Dlu2bMEdd9yB/fv3N7zyRhzZc4v8szoU\nW1D/1jE+Ph4ffPABxowZU+vP77nnHiQnJ+P5559vUCx0ZE8IqY+A9NlLpVKUl5eLyb68vBwcxyE4\nOBgOh6NhkfqBL0nan2pbltDNZDKJI5YIIaQl8ZrsFyxYgFGjRqFXr17gOA5HjhzBsmXLYDKZcPvt\ntzdFjC1SQUEBtm3bhgkTJsDhcOD777/Hvn378PXXXzd3aIQQUkO95rN3Op349ddfwXEcBg0adMOj\nW58qb4UXaFesWIHRo0ejsLAQU6ZMQVZWFvR6PSZMmID7778fw4YNa/D2W+p+E0JaFlq8pJVrq/tN\nCPFNQEbjEEIIaf0o2RNCSBtQr2R/+fJlrFu3DoBwYTInJyegQRFCCPEvr8n+ww8/xPTp07Fo0SIA\ngN1ux4wZM+q18bS0NLRr1w69evVqXJSEEEIaxWuyX7NmDTZv3gytVgsAiI6OrveUxw899BA2bdrU\nuAgJIYQ0mtdkHxwcDImk6m0XLlxAx44d67Xx4cOHIzQ0tOHREUII8QuvyX7WrFl44IEHUFpaikWL\nFmHChAmYO3duU8RGCCHET7zeHXXvvfdiwIAB+Oqrr+ByufDDDz8gJibGbwEsXLhQ/HdqaipSU1P9\ntm1CCLkZZGRkICMjo1HbCPhNVbm5uZg4cSKOHDlSs3K6qcpDW91vQohvAnJTlV6vR1BQEPR6vVj6\n9OmDZ555Bnl5eQ0O9maQmpoKtVot/l66d+/e3CERQkitvCb7Z555BvPnz8eRI0dw5MgR/O1vf8O4\nceMQERGBV1999YafnT59OoYOHYrTp08jJibGL0sZthROpxMcx2H58uUwGo0wGo04ceJEc4dFCCG1\n8tqN0717d48k5nK50KNHDxw/fhy9evXC0aNHG155K+vGqb4sYVZWFvr164fZs2djzpw5ftl+S91v\nQkjLEpBunCFDhuDtt99GcXExiouLsWzZMgwaNEhcc7UtcR/JL168GEVFRZDJZHjmmWcQExODxx9/\nHFlZWc0dIiGE1Mprsn/xxRdx+vRpDB48GIMHD8bp06fx8ssvw2KxYNmyZU0RY4viXpZQpVJh8eLF\nyMnJwYEDB9ChQwf83//9H3ieb+4QCSGkhtY7xTHnn2UJ4cPue1uWsHfv3njppZcwadKkBoVC3TiE\nkPoIyLKEdrsdO3bswObNm1FSUiK+/vHHH/seoT81U1K8UdcVJWtCSEvltRvnueeew4YNG/DNN9+g\nT58+OH78ONq1a9cUsbVoZWVlSE9Ph9VqRWFhIV5//XUUFhbiD3/4Q3OHRgghNXhN9tu2bcM777wD\ntVqNv/zlL/jpp5+wbdu2poitRXM4HHj++ecRGRmJlJQUXLhwAd99953HPEKEENJSeO3GkUql4DgO\nycnJ2LJlC7p27Qqz2dwUsbU41efxDw8Px/79+5sxGkIIqT+vyf6RRx5BcXEx/vrXv+Kpp57CpUuX\n8NJLLzVFbIQQQvzkhsne5XJBr9fDYDDAYDA0eiIeQgghzcPr0Mt+/frhwIED4Pw11LF65a3sDtpA\na6v7TQjxTUNyhddk/+KLL8JoNGLWrFno0KGD+LrBYGhYlNUrp2Tvoa3uNyHENwFJ9nFxcbUe1ftj\n0XFK9p7a6n4TQnwTkGQfSJTsPbXV/SaE+CYgE6HZbDZ88cUX+POf/wwAOHPmDDZu3NiwCAkhhDQL\nr0f2//jHP8AYw8aNG3Hs2DGYTCYMHTrULzM80pG9p7a634QQ3wTkyH7Hjh1YvHgxFAoFAECr1VJC\nIoSQVsZrsk9ISEBZWZn4/Ndff0VycnJAg2otli1bhpSUFKhUKjz00EM1fr5v3z6kpKTAYDBgwoQJ\nKB3wMCEAAA2fSURBVCgoaIYoCSGkHsl+3rx5uOuuu5CXl4dRo0YhLS0Njz/+eFPE1qI5nU5ER0fj\n+eefR1paWo2fV1RUYNy4cRg/fjwOHToEpVKJadOmNUOkhBBSj+kSBgwYgB07duDAgQPgeR4DBw5s\nirhapOrLEh4+fBjl5eWQSCT4/fffayy+/uWXXyI8PBwvvvgiAOEsIDo6Gjk5OYiPj2+O8AkhbZjX\nI/uJEyfi008/RWJiYptO9EDNZQndM1zWdg3j1KlT6NWrl/i8ffv2MBgMOHXqVJPFSwghbl6T/ZNP\nPoldu3YhKSkJkydPxpdffgmr1doUsbVI7mUJlUql+FptN50VFxcjLi7O47XOnTujqKgo0CESQkgN\nXrtxUlNTkZqaCqfTiR07duCjjz5CWloaysvLmyK+OnF+mpSNpab69P5BgwbV3EYtR/ZhYWE4ceKE\nx2vZ2dkICwvzqT5CCPEHr8keACwWC77//nusX78emZmZmDVrVqDj8srXJO0vtS1LWNuRfUJCAtav\nXy8+v3TpEoqLi5GQkBDQ+AghpDZeu3GmTp2KxMREbN++HY899hjOnTuHd955pylia/F4nofVaoXT\n6QTP87DZbOB5HgAwefJkFBcXY9GiRTh//jwee+wxjB49mi7OEkKahddkn5aWhuzsbHzwwQcYNWoU\n9uzZI06d0Na99NJL0Gg0WLx4MdauXQu1Wo1//etfAACdToeffvoJGzZsQN++fWG32/H55583c8SE\nkLaqXhOhZWZm4rPPPsP69esRHx+PyZMnY968eY2vnKZL8NBW95sQ4puG5Io6++xPnTqFzz77DF98\n8QUiIiJw7733gjFGq1URQkgrVOeRvUQiwYQJE7Bs2TJ06tQJABAfH++XeezFyunI3kNb3W9CiG/8\nOhHa119/DbVajREjRuDRRx/Ftm3bKBERQkgr5bXPvqKiAt999x0+++wz7NixAzNnzsTdd9+NMWPG\nNL5yOrL30Fb3mxDim4CvVFVcXIwvv/wSn3/+ObZv3+5zgDUqryNgg8GAkpKSRm+/tQkNDUVxcXFz\nh0EIaeFummUJCSGE1C0gi5c0xs6dO9G9e3d07dr1prwRq7WPTKL4mxfF33xac+wNFdBk/5e//AUf\nfPABtm7diuXLl6OwsDCQ1TW51v6FofibF8XffFpz7A0VsGTvXt1qxIgRiI2NxZgxY7Bv375AVUcI\nIeQGApbsf/vtNyQmJorPk5KS8OuvvwaqOkIIITcQsAu0W7duxYoVK/DZZ58BAN5//33k5+fjpZde\nqqq8ltkiCSGEeOe36RIaa8CAAXj66afF58eOHcO4ceM83kMjcQghpGkErBsnODgYgDAiJzc3F1u2\nbKl14Q9CCCGBF7AjewB466238Mc//hEOhwOPP/44wsPDA1kdIYSQOgR06OXIkSNx4sQJnD17Fo8/\n/rj4emsbf5+WloZ27dp5LCBuNBoxadIkdOrUCXfddRcqKiqaMcIbu3jxIkaNGoUePXogNTUVn376\nKYDWsw9WqxWDBg1C3759MXjwYCxZsgRA64nfjed5JCcnY+LEiQBaV/xxcXHo3bs3kpOTMXDgQACt\nJ36TyYRZs2ahW7duSEpKwr59+1pN7KdOnUJycrJYgoODsXTpUlRUVPgcf0CTfV1a2/j7hx56CJs2\nbfJ47b333kOnTp1w5swZdOzYEe+//34zReedXC7HkiVLcOzYMXz55Zd47rnnYDQaW80+qFQq7Nix\nA4cOHcLPP/+MFStW4MyZM60mfre3334bSUlJ4sCE1hQ/x3HIyMjAwYMHsX//fgCtJ/4FCxagU6dO\nOHz4MA4fPozExMRWE3tCQgIOHjyIgwcP4sCBA9BoNLj77rvx7rvv+hx/kyf71jj+fvjw4QgNDfV4\nbf/+/ZgzZw6USiXS0tJa9D5ERUWhb9++AIDw8HD06NEDv/32W6vaB41GA0CYmM/pdEKpVLaq+PPy\n8vDjjz9i7ty54sCE1hQ/UHNARWuJf+vWrXj22WehUqkgk8kQHBzcamKvbuvWrejSpQtiYmIaFj9r\nYlu2bGHTpk0Tn7/33nvsueeea+owfJaTk8N69uwpPu/UqROzWCz/v717DWnq/+MA/p7mDRSiHwmG\nRcPEZbpL05nrQTUlReYWmKmRaO5BSERpWqIkhGAGUayLBNWYWFhGD3LaUkQH3jDKuUDNSo2EyEwy\n5w3X/PwfiMN+P/1XdtHTvq9n5/Y97+9RPzv7nuM5REQ0MTFBmzZtWqloP+TVq1fE5/PJarVyqg92\nu52EQiG5urrSlStXiIhbP4P9+/dTR0cHmUwmUiqVRMSt/Hw+n4RCIanVanr48CERcSP/4OAgBQUF\nUVpaGslkMiopKaHJyUlOZP+3w4cP07Vr14hoecd+RYZx/gbEwdtGrVYrkpKScOnSJXh7e3OqDy4u\nLrBYLHj9+jVKS0thNps5k7+6uhq+vr6QSCRfZeZKfgBoaWmBxWLBuXPnkJ2djffv33Mi//T0NF6+\nfImEhASYTCZ0dXWhsrKSE9kXmpmZgcFgQGJiIoDl/e788WIfHh6OFy9eOKa7urqwY8eOPx3jp4WH\nh6OnpwcA0NPTg/Dw8BVO9P/ZbDYkJCQgNTUVarUaAPf6AMxdKIyLi0N7eztn8re2tqKqqgp8Ph8p\nKSloaGhAamoqZ/IDgJ+fHwBg69atUKlUMBgMnMi/ZcsWBAUFIT4+Hl5eXkhJScHjx485kX0ho9EI\nqVSK9evXA1je3+4fL/Z/y/33ERER0Ol0mJqagk6nW9UfWEQEjUaDkJAQnDhxwjGfK334+PEjRkdH\nAQAjIyOoq6uDWq3mTP7i4mIMDg5iYGAAd+/ehUKhQHl5OWfyT05Owmq1AgCGh4dRW1uL2NhYzuQP\nDAxEe3s7ZmdnUVNTg+joaM5kn1dRUYGUlBTH9LLy/46xpW8xmUwkEAgoICCAtFrtSkT4IcnJyeTn\n50fu7u7k7+9POp2OxsbGSKVS0caNG0mtVpPVal3pmEtqamoiHo9HIpGIxGIxicViMhqNnOnD8+fP\nSSKRkFAopL1791JZWRkREWfyL2QymSg+Pp6IuJO/v7+fRCIRiUQiUigUdOvWLSLiTv7e3l6KiIgg\nkUhEJ0+epPHxcc5kJyIaHx+nf/75h8bGxhzzlpN/RV9ewjAMw/wZ7AItwzCME2DFnmEYxgmwYs8w\nDOMEWLFnGIZxAqzYMyvO1dUVEokEUqkUubm5sNlsP7T9u3fvHP9sYrFYYDQaHcsMBgPOnz//S3Lu\n3Lnzl7SzXN7e3iu6f4bb2N04zIrz8fGB1WqFzWbDvn37kJmZCaVSuay29Ho9nj17tuqepmq32+Hq\n6vpTbcwfJ4ZZDnZmz6wabm5uUCgUaG5uxtjYGM6cOQOxWAyNRoO+vj4AQGdnJ6KioiAWi7F9+3ZM\nTEzgzZs3CA0Nhc1mQ2FhIe7duweJRILKykro9XocO3YMwNw3gOPHj0MkEiErKwtDQ0MAgPT0dJw+\nfRpyuRxhYWGor69fNN/8mbXJZEJUVBSSk5MRHByMgoKCRdffvXs3CgoKEBYWhsuXL6O3txeZmZmI\niIjA0aNHMTIyAgC4ceMGZDIZpFIpTp06hZmZGUdejUYDgUCA4uLiX3egGafEij2zaoyOjqK6uhox\nMTHQ6/X4/PkzOjo6oFAokJ+fD2DuhTj5+fno7OxEc3MzPD09Hdu7ubmhqKgIycnJMJvNOHDgwFfv\nOb5w4QL8/f1hsVjg6+uLixcvOpZZLBY0NDRAq9UuWVgXttXU1ISzZ8/CbDajqqoKg4ODi64/MDCA\n1tZWZGVlITc3F/n5+Whvb8e2bdtw8+ZNAEBCQgKePHmCp0+fYnJyEo2NjQCAkpISCAQC9PT0YHp6\n+ieOLMOwYs+sAlNTU5BIJFCr1YiPj8euXbtQU1OD9PR0uLi4ICkpCW1tbbDZbIiMjEReXh6uXr2K\nL1++/GdohIiWfEiU0WhERkYGAECj0cBgMACYK8qJiYnw9PREZGQkOjo6vplZJpMhKCgIHh4ekMvl\naG1tXXS9gwcPwt3dHR8+fEBzczNUKhUkEgmuX7+OlpYWAEB/fz8OHTqEkJAQPHr0CHV1dQCA2tpa\nZGRkgMfjOXIzzHL91tcSMsz38PLygtls/s/8fxdtHo+HI0eOIDo6Gnfu3IFQKPzh55Av9UGwdu1a\nAHNP17Tb7d9sZ+H7Ddzd3Zc8896wYQOAuTH7devWLdrPnJwcFBQU4Pbt29BqtbBYLN/MyzA/ip3Z\nM6uSUqlEeXk57HY77t+/D7lcjjVr1qCvrw8BAQEoLCyEQCBwjOXP27x5M4aHhx3TC4tlXFwcysrK\nMDs7C51OB5VK9dv7Mb9/Pz8/8Pl8PHjwAEQEm82G7u5uAHNj84GBgfj06RMqKioc28bGxjry6vX6\n356V+buxYs+suIVj4fPS0tLg4+MDqVSK+vp6xzi6VqtFaGgoZDIZBAIB5HL5V21ERkbCarU6LtDy\neDzHspycHLx9+xYSiQRDQ0PIzs5eNMNieb61zvdsU1paisbGRojFYkgkErS1tQEAioqKoFQqERMT\ngz179jjWz8vLQ3d3N4KDg+Hh4bHkPhjme7BbLxmGYZwAO7NnGIZxAqzYMwzDOAFW7BmGYZwAK/YM\nwzBOgBV7hmEYJ8CKPcMwjBP4Hxd4aIMPurYGAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 13 }, { "cell_type": "code", "collapsed": false, "input": [ "posns = [ int(line) for line in open('biased-random-positions.out') ]\n", "hist(posns)\n", "xlabel('Error position within read')\n", "ylabel('Number of errors')\n", "savefig('/tmp/foo2.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEMCAYAAAAMMiuwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1YVHWiB/DvzIooSpq0IC2MCBIzowgzCoNwjUn37uIr\nphVU+Kiw6jXLzaTNm63h1lNr1vWlXfGl6M31Qdvirtkam9qIdGUGDfU+MJJ419BdRJGUAZEX/d0/\n2DlBgCAcZgb4fp5nngd+nJfvOeJ8OefMnFEIIQSIiIhkoHR2ACIi6jtYKkREJBuWChERyYalQkRE\nsmGpEBGRbFgqREQkG9lL5ebNmzAYDAgPD0dUVBQ2btwIALDZbIiPj4dKpcKcOXNQXV0tzbNlyxYE\nBwdDq9UiNzdXGrdardDr9QgMDMSaNWuk8YaGBqSkpGDUqFEwGo24dOmS3JtBRERdIHupDBo0CF99\n9RVOnjyJI0eO4N1338XZs2eRnp4OlUqFs2fPws/PD9u2bQMAXL58GVu3bsWhQ4eQnp6OFStWSMta\ntWoVXnjhBeTn5+PIkSM4fvw4ACArKwvXr1+H1WpFXFwcXn31Vbk3g4iIuqBHTn95eHgAAKqrq9HY\n2Ah3d3dYLBakpKTA3d0dycnJMJvNAACz2Yy4uDioVCrExsZCCCEdxRQXFyMhIQFeXl6YO3dui3mS\nkpLg4eGBJUuWSONERORcPVIqt2/fRlhYGHx8fPD0009DpVIhPz8farUaAKBWq2GxWAA0FYRGo5Hm\nDQkJgdlsRklJCby9vaVxrVaLvLw8AIDFYoFWqwUAjBgxAuXl5airq+uJTSEiorswoCcWqlQqcerU\nKZw/fx7Tp09HTEwM7uZuMAqFotWYEEIaF0K0WF57y25rOUREdGfduXtXj776KyAgANOnT4fZbEZE\nRASsViuApgvwERERAACDwYCioiJpnjNnziAiIgJjxoxBeXm5NF5UVASDwdBqnsrKSvj4+MDd3b3N\nDPYC6u+Pl19+2ekZXOXBfcH9wH3R/qO7ZC+ViooKXLt2DQBw9epV/O1vf0N8fDwMBgMyMjJQW1uL\njIwMREVFAQAiIyORnZ2N0tJSmEwmKJVKeHp6Amg6TZaZmYmKigpkZWW1KJVdu3ahpqYGO3bskJZF\nRETOJfvpr7KyMixYsAC3bt3CyJEjkZqaCl9fXyxbtgxJSUkICQmBXq/H+vXrAQA+Pj5YtmwZpkyZ\ngoEDB2L79u3Sst58800kJSXhP//zP5GYmIiJEycCAB5++GF88cUX0Gg0CAwMRGZmptybQUREXaAQ\nchzvuCiFQiHL4VxfYDKZYDQanR3DJXBfNOF++AH3xQ+6+7zJUiEiIkl3nzd5mxYiIpINS4WIiGTD\nUiEiItmwVIiISDYsFSIikg1LhYiIZMNSISIi2bBUiIhINiwVIiKSDUuFiIhkw1IhIiLZsFSIiEg2\nLBUiIpINS4WIiGTDUiEiItmwVIiISDayf5wwEbm2e+4ZAZvte2fHAOAGoKHfZ/D0vBdVVZVOzSAn\nfvIjUT+jUCgAuML/C1fI4QoZ3AA0OjVB82LjxwnfAUuFXI3rHCW4wv8LV3hCZwZ7BvtzJUvlDlgq\n5Gpc4yjBFTIArpGDGewZ5CoVXqgnIiLZsFSIiEg2LBUiIpINS4WIiGTDUiEiItmwVIiISDYsFSIi\nko3spXLhwgU89NBDGDt2LIxGI3bv3g0ASEtLg5+fH3Q6HXQ6HQ4cOCDNs2XLFgQHB0Or1SI3N1ca\nt1qt0Ov1CAwMxJo1a6TxhoYGpKSkYNSoUTAajbh06ZLcm0FERF0hZFZWViYKCgqEEEJcuXJFjB49\nWlRVVYm0tDTx1ltvtZq+vLxchISEiO+++06YTCah0+mkn02bNk1kZmaKiooKERMTI/Lz84UQQuzZ\ns0fMmzdP1NTUiNdff10sX768zSw9sHlE3QJAAMLJD1fI4Co5mMGeofnvaHfIfqQycuRIhIeHAwDu\nu+8+jB07Fvn5+fYCazW92WxGXFwcVCoVYmNjIYRAdXU1AKC4uBgJCQnw8vLC3LlzYTabpXmSkpLg\n4eGBJUuWSONERORcPXpNpaSkBIWFhTAYDACAt99+G1FRUVi/fj1sNhsAwGKxQKPRSPOEhITAbDaj\npKQE3t7e0rhWq0VeXp40j1arBQCMGDEC5eXlqKur68lNISKiTuixW9/bbDYkJCRg48aNGDJkCJYt\nW4a1a9eiqqoKzz//PLZv347U1NQ2j16a7o/UkhBCGhdCtJivrWXYpaWlSV8bjUYYjcaubxQRUZ9j\nAtDyubJbunXyrB319fXi3//938XGjRvb/PnJkydFdHS0EEKIffv2iRUrVkg/CwsLE1VVVUIIIUaP\nHi2Nv/nmm+IPf/iDEEKI5557Tnz66adCCCGuXr0qJkyY0OZ6emjziLoMLnL+3PkZXCUHM9gzNP8d\n7Q7ZT38JIZCSkoJx48bh2WeflcbLysoAAI2Njdi9ezemT58OAIiMjER2djZKS0thMpmgVCrh6ekJ\nAFCr1cjMzERFRQWysrKk02gGgwG7du1CTU0NduzYgaioKLk3g4iIuqKbf3i1cvToUaFQKERYWJgI\nDw8X4eHh4q9//auYP3++CA0NFRMmTBArV64UV69elebZtGmTCAoKEhqNRuTk5EjjhYWFQqfTiYCA\nALF69WppvL6+XixatEj4+/uL2NhYUVZW1maWHtg8om6Bi/xV6vwMrpKDGewZmv+Odgc/T4XIgfh5\nKs25Qg5msGewP1fy81SIiMhlsFSIiEg2LBUiIpINS4WIiGTDUiEiItmwVIiISDYsFSIikg1LhYiI\nZMNSISIi2bBUiIhINiwVIiKSDUuFiIhkw1IhIiLZsFSIiEg2LBUiIpINS4WIiGTDUiEiItmwVIiI\nSDYsFSIikg1LhYiIZMNSISIi2bBUiIhINiwVIiKSDUuFiIhkw1IhIiLZsFSIiEg2LBUiIpINS4WI\niGTTYans3bsXVVVVAICtW7di8eLFKCkpaXf6Cxcu4KGHHsLYsWNhNBqxe/duAIDNZkN8fDxUKhXm\nzJmD6upqaZ4tW7YgODgYWq0Wubm50rjVaoVer0dgYCDWrFkjjTc0NCAlJQWjRo2C0WjEpUuX7n7L\niYhIfqID48aNE0IIcfr0aWEwGMTu3bvFjBkz2p2+rKxMFBQUCCGEuHLlihg9erSoqqoS69evF08/\n/bS4efOmWL58udiwYYMQQojy8nIREhIivvvuO2EymYROp5OWNW3aNJGZmSkqKipETEyMyM/PF0II\nsWfPHjFv3jxRU1MjXn/9dbF8+fI2s3Ri84gcCoAAhJMfrpDBVXIwgz1D89/R7ujwSMXNzQ0A8P77\n7+Opp57C448/jn/+85/tTj9y5EiEh4cDAO677z6MHTsW+fn5sFgsSElJgbu7O5KTk2E2mwEAZrMZ\ncXFxUKlUiI2NhRBCOoopLi5GQkICvLy8MHfu3BbzJCUlwcPDA0uWLJHGiYjIuToslbCwMMyfPx/7\n9+/HY489hps3b+LWrVudWnhJSQkKCwsRGRmJ/Px8qNVqAIBarYbFYgHQVBAajUaaJyQkBGazGSUl\nJfD29pbGtVot8vLyAAAWiwVarRYAMGLECJSXl6Ourq6Tm0xERD1lQEcTvPfee/jqq6+wYcMGDBo0\nCGVlZdiwYUOHC7bZbEhISMDGjRsxdOhQNB1VdY5CoWg1JoSQxoUQLZZ3p2WnpaVJXxuNRhiNxk7n\nICLq+0wAWj5XdscdS6WxsRHjx49HUVGRNObr6wtfX987LrShoQHz5s3D/PnzER8fDwCIiIiA1WqF\nTqeD1WpFREQEAMBgMODgwYPSvGfOnEFERAQ8PT1RXl4ujRcVFcFgMEjzFBUVISQkBJWVlfDx8YG7\nu3ubWeTaUUREfZMRwA/PlevWrevW0u54+mvAgAHQaDQoKCjo9AKFEEhJScG4cePw7LPPSuMGgwEZ\nGRmora1FRkYGoqKiAACRkZHIzs5GaWkpTCYTlEolPD09ATSdJsvMzERFRQWysrJalMquXbtQU1OD\nHTt2SMsiIiIn6+hKvtFoFEqlUuj1ejFz5kwxc+ZMMWvWrHanP3r0qFAoFCIsLEyEh4eL8PBwceDA\nAVFVVSVmz54t/P39RXx8vLDZbNI8mzZtEkFBQUKj0YicnBxpvLCwUOh0OhEQECBWr14tjdfX14tF\nixYJf39/ERsbK8rKytrM0onNI3IouMgrfZyfwVVyMIM9Q/Pf0e5Q/Gsh7TKZTNLXCoUCQjRd24iN\nje2pnpONPS+Rq2i6Lujs30lXyAC4Rg5msGewP1d293mzw1KxM5vNUCgUiIyM7PLKHI2lQq6GpdKc\nK+RgBnsGuUqlw1d/mUwmLF68GA888AAA4OzZs9i5c2evOFIhIiLH6rBUNmzYgP379yMkJAQA8O23\n3+LZZ59lqRARUSsdvvnx+++/x8iRI6XvfXx8cO3atR4NRUREvVOHRyoLFizAtGnT8Mgjj0AIgays\nLCxcuNAB0YiIqLe544V6IQQuXLiAa9euYf/+/VAoFJg5cyZCQ0MdmbHLeKGeXA0v1DfnCjmYwZ7B\nIa/+EkJg/Pjx+N///d8ur8CZWCrkalgqzblCDmawZ5CrVO54TUWhUGDSpEn4y1/+0uUVEBFR/9Hh\n+1Q0Gg2Ki4vh5eUlXbBXKBQ4ffq0QwJ2B49UyNXwSKU5V8jBDPYMDnmfihAC27dvh0ql6vIKiIio\n/+A1FSIH4pFKc66QgxnsGXhNhYiIXA6vqRA5EI9UmnOFHMxgz+Cwe38dOHCgywsnIqL+pcPbtAQE\nBMDd3R1ff/01AgICMGTIEP71T0REbeqwVHbs2IHHH39c+ojJ+vp6JCUl9XgwIiLqfToslY8++gh/\n+9vfMGTIEADAz372M9hsth4PRkREvU+HpTJs2DAolT9MVlpaCj8/vx4NRUREvVOHpbJgwQI8+eST\nuHbtGtatW4eZM2fiV7/6lSOyERFRL9OpjxM+f/48PvnkE9y+fRuJiYnw9/d3RLZu40uKydXwJcXN\nuUIOZrBncPhn1PdGLBVyNSyV5lwhBzPYMzjkHfVERER3g6VCRESyabdUpk6dCgD4zW9+47AwRETU\nu7V7m5ba2lqYTCbs27cPiYmJEEL863xwE71e75CARETUe7R7of7LL79Eeno6vvzyS0ycOLHVz7/6\n6qseD9ddvFBProYX6ptzhRzMYM/gsFd//e53v8PatWu7vAJnYqmQq2GpNOcKOZjBnsGhLym2Wq3Y\nt28fFAoFZs+eDbVa3eUVOhJLhVwNS6U5V8jBDPYMDntJ8TvvvIOFCxdKt2pZtGgR3nnnnXanT05O\nho+PD0JDQ6WxtLQ0+Pn5QafTQafTtbid/pYtWxAcHAytVovc3Fxp3Gq1Qq/XIzAwEGvWrJHGGxoa\nkJKSglGjRsFoNOLSpUt3t8VERNRzRAeio6NFZWWl9H1lZaWYNGlSu9Pn5OSIb775RowbN04aS0tL\nE2+99VaracvLy0VISIj47rvvhMlkEjqdTvrZtGnTRGZmpqioqBAxMTEiPz9fCCHEnj17xLx580RN\nTY14/fXXxfLly9vN0onNI3IoAAIQTn64QgZXycEM9gzNf0e7o8MjleHDh+Pq1avS95WVlRg+fHi7\n00+ePBn33ntvW+XVasxsNiMuLg4qlQqxsbEQQqC6uhoAUFxcjISEBHh5eWHu3Lkwm83SPElJSfDw\n8MCSJUukcSIicr4OS+W5555DXFwcZs2ahVmzZiEuLg6pqal3vaK3334bUVFRWL9+vXTrfIvFAo1G\nI00TEhICs9mMkpISeHt7S+NarRZ5eXnSPFqtFgAwYsQIlJeXo66u7q7zEBGR/Dr8OOGpU6fi22+/\nRV5eHhQKBQwGQ4tb4XfGsmXLsHbtWlRVVeH555/H9u3bkZqa2ubRS/P3wtgJ8cN7ZIQQLeZraxnN\npaWlSV8bjUYYjca7yk5E1LeZALR8ruyODksFAJRKJaKjo7u8EvtRx7Bhw7B8+XI89dRTSE1NhcFg\nwMGDB6Xpzpw5g4iICHh6eqK8vFwaLyoqgsFgAAAYDAYUFRUhJCQElZWV8PHxgbu7e7vrlmtHERH1\nTUYAPzxX2j/lt6sccu+vsrIyAEBjYyN2796N6dOnAwAiIyORnZ2N0tJSmEwmKJVKeHp6AgDUajUy\nMzNRUVGBrKysFqWya9cu1NTUYMeOHYiKinLEJhARUWd06zJ/GxITE4Wvr69wc3MTfn5+4t133xXz\n588XoaGhYsKECWLlypXi6tWr0vSbNm0SQUFBQqPRiJycHGm8sLBQ6HQ6ERAQIFavXi2N19fXi0WL\nFgl/f38RGxsrysrK2s3SA5tH1C1wkVf6OD+Dq+RgBnuG5r+j3XHHNz82NjYiNDQUVqvVEf0mO775\nkVwN3/zYnCvkYAZ7BvtzZY+++XHAgAHQarUoKCjo8gqIiKj/6PBCfWVlJSZOnIjw8HDcf//9AJqa\nbN++fT0ejoiIepcOS+Xll19uNdbWy36JiIg6dUPJ+vp65OXl4cEHH8SNGzfQ2NiIe+65xxH5uoXX\nVMjV8JpKc66QgxnsGRxyTQUAPv30U0RFRWHRokUAgIsXL+Lhhx/u8gqJiKjv6rBUtm7diqNHj0pH\nJg888AAuX77c48GIiKj36fCaikKhgIeHh/T9lStX4OXl1aOh5OT86z9uABqcmsDT815UVVU6NcM9\n94yAzfa9UzO4wn4g6us6LJXHHnsMqampuHHjBj744AN8+OGHmD9/viOyycT55yqdncFmc3OBcgWc\nvx9cYR8Q9W0dXqgXQuDIkSP45JNPcPv2bTzxxBOIiYlxVL5u4UVRZmjJDUCjkzMAzt8PrvBvAbhG\nDmawZ5DrQn2nXv0FAP/4xz8AAD/72c+6vDJHY6kwAzO4agbANXIwgz2DXKXS4ekvs9mMxYsXSytR\nKpXYuXMnIiMju7xSIiLqmzo8UjEajXjjjTekEsnPz8fzzz8Pk8nkiHzdwiMVZmAGV80AuEYOZrBn\ncNj7VGw2W4tPZ9RoNNInNxIRETXX7umvTz75BEDTkcqMGTPw8MMPQwiBv/zlL4iNjXVYQCIi6j3a\nLZXPPvtMehlqYGAgTp8+DQAYPXo0vv/eue83ICIi19TpV3/1RrymwgzM4KoZANfIwQz2DA579dfF\nixexZ88eHDt2DHV1ddJKeet7IiL6sQ5LZfHixYiKisLSpUvh5uYGwBVufUJERK6ow9NfEydOhMVi\ngVLZ4QvFXA5PfzEDM7hqBsA1cjCDPYPD3lGflZUFk8mE+Ph4DB8+XBrX6/VdXqmjsFSYgRlcNQPg\nGjmYwZ7BYddUiouL8eGHH+L48eMYOHCgNP7VV191eaVERNQ3dXikMmbMGJw8eRJDhw51VCbZ8EiF\nGZjBVTMArpGDGewZHPaO+rCwMJSXl3d5BURE1H90ePrr2rVr0Gq1iIyMlK6p8CXFRETUlg5L5be/\n/a0jchARUR/Ad9T3fApmYAZmaJMr5GAGewaHvfpr6NCh0psd6+rq0NjYiKFDh6KqqqrLKyUior6p\nw1Kprq6Wvr5x4wY+/PBDXLp0qUdDERFR73RXb5P38PDAf/zHf2Dv3r3tTpOcnAwfHx+EhoZKYzab\nDfHx8VCpVJgzZ06LotqyZQuCg4Oh1WqRm5srjVutVuj1egQGBmLNmjXSeENDA1JSUjBq1CgYjUYW\nHBGRC+mwVD755BPpsXv3bixduhTh4eHtTr9o0SJ88cUXLcbS09OhUqlw9uxZ+Pn5Ydu2bQCAy5cv\nY+vWrTh06BDS09OxYsUKaZ5Vq1bhhRdeQH5+Po4cOYLjx48DaHqH//Xr12G1WhEXF4dXX321SxtO\nRETy6/D0V/PPVRk0aBBiYmIwc+bMdqefPHkyzp8/32LMYrHgpZdegru7O5KTk/H6668DAMxmM+Li\n4qBSqaBSqSCEQHV1NYYOHYri4mIkJCQAAObOnQuz2YyJEyfCbDYjKSkJHh4eWLJkCX75y192dduJ\niEhmHZbK+++/3+2V5OfnQ61WAwDUajUsFguAplJp/lHFISEhMJvNGDVqFLy9vaVxrVaLP/3pT1i+\nfDksFguWLl0KABgxYgTKy8tRV1cHd3f3buckIqLuabdU1q1b1+a4/ahl7dq1nV7J3bw8ra3b6gsh\npHEhRIvldbzstGZfG//1ICKiJiYAQFpamixLa7dUhgwZ0uoJvqamBu+++y4qKiruqlQiIiJgtVqh\n0+lgtVoREREBADAYDDh48KA03ZkzZxAREQFPT88Wt4YpKiqCwWCQ5ikqKkJISAgqKyvh4+PTwVFK\nWqdzEhH1P0YAP5RKewcUndXuhfrU1FSsWrUKq1atwuLFi1FbW4v33nsPiYmJ+Pvf/35XKzEYDMjI\nyEBtbS0yMjIQFRUFAIiMjER2djZKS0thMpmgVCrh6ekJoOk0WWZmJioqKpCVldWiVHbt2oWamhrs\n2LFDWhYREbkAcQcVFRVizZo1IiAgQKxdu1ZUVlbeaXIhhBCJiYnC19dXDBw4UPj5+YmMjAxRVVUl\nZs+eLfz9/UV8fLyw2WzS9Js2bRJBQUFCo9GInJwcabywsFDodDoREBAgVq9eLY3X19eLRYsWCX9/\nfxEbGyvKysrazQJAAMLJD2ZgBmZw3RzMYM/Q/HmzO9q9TUtqaiqysrKwZMkSPPXUU9IRRG/C27Qw\nAzO4agbANXIwgz2DvQp67JMflUolBg4cKH0ufYuZFIpecZsWlgozMIOrZgBcIwcz2DPIVSrtXqi/\nfft2lxdKRET9013dpoWIiOhOWCpERCQblgoREcmGpUJERLJhqRARkWxYKkREJBuWChERyYalQkRE\nsmGpEBGRbFgqREQkG5YKERHJhqVCRESyYakQEZFsWCpERCQblgoREcmGpUJERLJhqRARkWxYKkRE\nJBuWChERyYalQkREsmGpEBGRbFgqREQkG5YKERHJhqVCRESyYakQEZFsWCpERCQblgoREcnGoaUS\nEBCA8ePHQ6fTITIyEgBgs9kQHx8PlUqFOXPmoLq6Wpp+y5YtCA4OhlarRW5urjRutVqh1+sRGBiI\nNWvWOHITiIjoDhxaKgqFAiaTCQUFBbBYLACA9PR0qFQqnD17Fn5+fti2bRsA4PLly9i6dSsOHTqE\n9PR0rFixQlrOqlWr8MILLyA/Px9HjhzB8ePHHbkZRETUDoef/hJCtPjeYrEgJSUF7u7uSE5Ohtls\nBgCYzWbExcVBpVIhNjYWQgjpKKa4uBgJCQnw8vLC3LlzpXmIiMi5HH6kMmXKFMyZMwf79u0DAOTn\n50OtVgMA1Gq1dARjNpuh0WikeUNCQmA2m1FSUgJvb29pXKvVIi8vz4FbQURE7RngyJV9/fXX8PX1\nhdVqxaxZsxAZGdnqyOVOFApFq7GO509r9rXxXw8iImpiAgCkpaXJsjSHloqvry8AQKPRYPbs2fjs\ns88QEREBq9UKnU4Hq9WKiIgIAIDBYMDBgwelec+cOYOIiAh4enqivLxcGi8qKkJUVNQd1prWE5tC\nRNRHGAH8UCrr1q3r1tIcdvrrxo0bsNlsAIArV64gOzsbcXFxMBgMyMjIQG1tLTIyMqSCiIyMRHZ2\nNkpLS2EymaBUKuHp6Qmg6TRZZmYmKioqkJWVBYPB4KjNICKiOxEO8n//938iLCxMhIWFiSlTpoh3\n331XCCFEVVWVmD17tvD39xfx8fHCZrNJ82zatEkEBQUJjUYjcnJypPHCwkKh0+lEQECAWL16dbvr\nBCAA4eQHMzADM7huDmawZ2j+vNkdin8tpE9qugbj7M1jBmZghra5Qg5msGewV4FC8cPXXcF31BMR\nkWxYKkREJBuWChERyYalQkREsmGpEBGRbFgqREQkG5YKERHJhqVCRESyYakQEZFsWCpERCQblgoR\nEcmGpUJERLJhqRARkWxYKkREJBuWChERyYalQkREsmGpEBGRbFgqREQkG5YKERHJhqVCRESyYakQ\nEZFsWCpERCQblgoREcmGpUJERLJhqRARkWxYKkREJBuWChERyYalQkREsum1pZKTkwONRoPg4GC8\n/fbbzo7TC5icHcCFmJwdwEWYnB3AhZicHaDP6LWl8utf/xrbt2/HwYMH8cc//hEVFRXOjuTiTM4O\n4EJMzg7gIkzODuBCTM4O0Gf0ylK5fv06AODBBx/EqFGj8Itf/AJms9nJqYiIqFeWSn5+PtRqtfS9\nVqtFXl6eExMREREADHB2gJ6ncHYAuEaGdf96OJMr7Ad7BmfuC+6HH3BftMzgxAQKeTL0ylKJiIjA\n888/L31fWFiIuLi4VtMJIRwZi4io3+uVp7+GDRsGoOkVYOfPn8eXX34Jg8Hg5FRERNQrj1QAYNOm\nTVi6dCkaGhqwYsUK3Hfffc6ORETU7/XKIxUAiI2NhdVqRUlJCVasWNHiZ/31PSwXLlzAQw89hLFj\nx8JoNGL37t0AAJvNhvj4eKhUKsyZMwfV1dVOTuo4t27dgk6nw6xZswD0331RU1ODBQsW4IEHHoBW\nq4XZbO63+2Lnzp2Ijo7GhAkT8OyzzwLoP78XycnJ8PHxQWhoqDR2p23fsmULgoODodVqkZub26l1\n9NpSuZP++h4WNzc3bNy4EYWFhfjzn/+Ml156CTabDenp6VCpVDh79iz8/Pywbds2Z0d1mM2bN0Or\n1UoXIfvrvnj55ZehUqlw+vRpnD59Gmq1ul/ui8rKSrz22mv48ssvkZ+fj2+//RbZ2dn9Zl8sWrQI\nX3zxRYux9rb98uXL2Lp1Kw4dOoT09PRWf7y3p8+VSn9+D8vIkSMRHh4OALjvvvswduxY5Ofnw2Kx\nICUlBe7u7khOTu43++PixYv461//il/96lfSizb66744ePAgXnzxRQwaNAgDBgzAsGHD+uW+GDx4\nMIQQuH79Ompra3Hjxg0MHz683+yLyZMn4957720x1t62m81mxMXFQaVSITY2FkII2Gy2DtfR50qF\n72FpUlJSgsLCQkRGRrbYJ2q1GhaLxcnpHGPlypXYsGEDlMoffs374764ePEibt68iWXLlsFgMGD9\n+vWora2FuTZRAAAJqElEQVTtl/ti8ODBSE9PR0BAAEaOHImYmBgYDIZ+uS/s2tt2s9kMjUYjTRcS\nEtKp/dLnSoWazpEmJCRg48aNGDp0aL98afX+/fvh7e0NnU7XYvv74764efMmvv32W8ybNw8mkwmF\nhYXYu3dvv9wXV65cwbJly1BUVITz58/j2LFj2L9/f7/cF3Z3s+2deS9LnyuViIgInDlzRvq+sLAQ\nUVFRTkzkWA0NDZg3bx7mz5+P+Ph4AE37xGq1AgCsVisiIiKcGdEh/ud//gf79u3D6NGj8fjjj+Pw\n4cOYP39+v9wXY8aMQUhICGbNmoXBgwfj8ccfxxdffNEv94XFYkFUVBTGjBkDLy8vPProozh69Gi/\n3Bd27W27wWBAUVGRNN2ZM2c6tV/6XKn05/ewCCGQkpKCcePGSa9qAZp+OTIyMlBbW4uMjIx+UbKv\nvfYaLly4gL///e/IzMzElClT8NFHH/XLfQEAwcHBMJvNuH37Nj7//HP8/Oc/75f7YvLkyTh+/Dgq\nKytRV1eHAwcO4Be/+EW/3Bd27W17ZGQksrOzUVpaCpPJBKVSCU9Pz44XKPogk8kk1Gq1CAoKEps3\nb3Z2HIc5evSoUCgUIiwsTISHh4vw8HBx4MABUVVVJWbPni38/f1FfHy8sNlszo7qUCaTScyaNUsI\nIfrtviguLhYGg0GEhYWJVatWierq6n67L9577z3x4IMPiokTJ4qXXnpJ3Lp1q9/si8TEROHr6ysG\nDhwo/Pz8REZGxh23fdOmTSIoKEhoNBqRk5PTqXUohOjHJxOJiEhWfe70FxEROQ9LhYiIZMNSISIi\n2bBUiIhINiwVcqqf/OQn0Ol00uONN95wdqQ7+uc//4lHH30UAHDq1CkcOHBA+tlnn32G9evXOyTH\njBkzUFVVhevXryM9PV0aN5lM0s0zf2zx4sXS+xGcwWg04sSJE05bPzlGr731PfUNHh4eKCgouOM0\nt2/fbnGrlR9/39n55HD//ffj448/BgAUFBTgxIkTmDZtGgBg1qxZ7T6hy+3zzz8HAJw/fx5bt27F\nsmXLOpxn586dXV5fY2MjBgzo3tOFQqGQ7dMFyXXxSIVcUkBAAH73u99Br9fj448/bvH9n//8Zxw6\ndAgzZsxATEwM3nnnHWm+oUOH4re//S3Cw8Nb3fPNaDTixRdfRGhoKOLj46U7L1y/fl2aJyUlBefO\nnQMAnDx5ElOnTkV4eDj0ej1qampw/vx5hIaGoqGhAWvXrsWePXug0+mwd+9evP/++3jmmWcANB3R\n/PrXv0ZYWBhWrlyJ8vJyAMDChQvxwgsvIDo6GhMnTsTBgwdbbfuGDRukj2xYuXIlpk6dCgA4fPgw\nkpKSpP1z9epVrF69GufOnYNOp8NvfvMbKBQK1NbWIjExEVqtFmvWrGmx/d988420n1555RWMHTsW\nTzzxBCorK1vlSEtLw5IlSxATE4OFCxfi2rVrWLduHWJiYvDoo4/i5MmTAJrepR4dHQ2dTocFCxbg\n/PnzAID6+nq88sor0Gq1SExMRH19fb++HUp/wVIhp6qtrW1x+st+FKBQKGCz2XDixAkkJCS0+P6R\nRx7B0qVLsXnzZuzfvx87d+6UTuvcuHEDP/3pT3Hy5ElER0e3WJdCocC5c+dw4sQJPPnkk9JHUn/w\nwQe4fv06vvnmG0yZMgUvvvgigKYPgnvxxRdx8uRJ5ObmYtCgQdKy3Nzc8MorryAxMREFBQV47LHH\nWvwV/uabb8LPzw+nTp2Ct7c3/uu//kv62alTp3D48GFs3rwZr732Wqt98uCDD+Lo0aMAgOPHj6Om\npgaNjY04evQoYmNjpW1RKBRYv349goKCUFBQgDfeeANCCOTk5GDdunUoKCjAvn37cPHiRWkeuxs3\nbuD+++9HYWEhhgwZgv3797f575OTk4P//u//xq5du7Bp0yaEh4fj66+/xssvvywVlkajwdGjR1FQ\nUIAZM2Zg+/btAJruv3by5EkUFBRg8eLFOHbsGI9U+gGWCjnV4MGDUVBQID3s1ysAYMGCBS2ehOzf\n5+XlQaPRYMyYMbj33nvxyCOPYN++fQAApVKJhQsXtru+xMREDBw4EPPmzcM333yD+vp6fP7551i4\ncCGUSiUSEhJw7NgxNDQ0YNKkSVi9ejX+8Ic/oLGxET/5yU9aLEsI0e5f3gcOHEBycjIAICUlBZ99\n9hmApif2Rx99FIMGDcKkSZOkI4fm9Ho9Tpw4AZvNJk13/Phx5ObmYvLkya0y/FhkZCRCQkLg7u6O\n6OhofP31162mGTBgAJ588kkAwJQpU3Ds2LFW0ygUCsycORM//elPAQCffvop0tLSoNPpMH/+fBQV\nFeHmzZuora3FypUrERYWhldffRXZ2dnSPnjyySfh7u6OqVOnQqVStbmvqG/hNRVyWffff3+b3//4\nr10hhDQ2ePBg3HPPPe0us/mTcPPl/PjJWaFQYOnSpfj5z3+OP/3pTxg/fvxdf8ZGe4UzfPhwAE0F\neOvWrVY/d3Nzw+jRo/H+++8jOjoa48ePx+HDh1FSUtLiYx3a0/zzMgYOHIi6urpW07i7u0tHXm5u\nbrh582aby/L19ZW+vnXrFvbv39+qHH7/+9/Dy8sLx48fR2FhIR5++GEATfuwvf1NfRePVKjXMRgM\nOHPmDM6dO4fvv/8eWVlZmD17dofzCSGwd+9e1NfXIysrC3q9HgMHDsTMmTPx0Ucf4datW/j4448R\nHR2NAQMG4Ny5cwgKCsLatWuhVqulay12AQEBuHLlSovl202fPh0ffPABbt++jYyMjE7la27y5Ml4\n8803ERsbi8mTJ2Pbtm3Q6/WtpvPx8UFVVdVdLburnnjiCbz99ttSSZ06dQoA8I9//AOjR48G0PLF\nANOmTUNmZibq6upw+PBhfPfddw7JSc7FUiGn+vE1Ffv1jB9r/leuUqnE9u3b8cwzz2DGjBlISUmR\n/oK/01/DCoUCgYGBmDBhAj766CNs2LABQNNpNU9PT0yYMAEHDx6UrnNs3rwZoaGhiIyMhFqtlq7R\n2NcxadIk2Gw26UJ981c3paamorS0FDqdDuXl5Xjuuefa3Jb28k6ePBmXLl3CpEmT4O3tjcGDB7c4\n9dX8yCwhIQF6vV66UN+ZI4IfZ2hvnubjTz/9NIYNG4Z/+7d/w9ixY6VrJ8888wy2b9+OiRMnwt/f\nX5pnxowZGD9+PPR6PXbs2NHqGhf1TbyhJPUbDz30EN566602/+InInnwSIWIiGTDIxUiIpINj1SI\niEg2LBUiIpINS4WIiGTDUiEiItmwVIiISDYsFSIiks3/A6DueCNkVG8fAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 18 }, { "cell_type": "code", "collapsed": false, "input": [ "bias_rand_rep = numpy.loadtxt('repreads-bias-random.dens')\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 19 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(bias_rand_rep[:,0], bias_rand_rep[:,3], label='w/repeats')\n", "plot(bias_rand[:,0], bias_rand[:,3], label='no repeats')\n", "axis(ymin=0)\n", "xlabel(\"Position in read\")\n", "ylabel(\"Average graph density within 10 steps\")\n", "title(\"Effects of errors on assembly graph density\")\n", "legend(loc='lower right')\n", "savefig('/tmp/bias2.png')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEVCAYAAAABwEUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlcVOX+B/DPDPsywIALyDYQqImyCEq4gZZmJmlqppni\ndSn9ldXVsj20W5plWl69WV3NRNM0l9ISl3LEBVe2BFREEEVAkH0ZmOX5/XHuHBxZhoEZYOT7fr2e\n18w5M+c83zkw3/PMc55zjoAxxkAIIeShJuzoAAghhBgeJXtCCOkCKNkTQkgXQMmeEEK6AEr2hBDS\nBVCyJ4SQLoCSvY42bNiAvn37olevXgCAv/76CyNGjIBIJEJycnIHR9dQSUkJZs+eje7du+Ott97q\n6HBII7KzsyEUCqFSqRp9fdmyZZg5c2Y7R6W7LVu2YPjw4a1efvbs2fjwww/1GFFD/fv3R1xcnEHr\n6KxMOzqAzkYikeDu3bswMTHh5/3jH//AunXrUFdXh3fffRdJSUnw9vYGAKxYsQJz5szB7Nmz21Sv\nUCjE9evX+fXqy++//468vDzk5ubC3Nxcr+sm7UMgEHR0CO1CIBAY/LNevnyZf75s2TJkZmYiJibG\noHV2FpTsHyAQCHDw4EGMGjWqwWsJCQkQiUR8QmaM4fTp0/jmm2/0Urchzm87deoUQkJC9JLoFQoF\nTE01/2WUSqXGjlEbXd9PDPN/0ZjO8LehczwNh7pxWujYsWMYNWoU8vPzIRKJ8MILL8DOzg61tbUI\nCgqCr68vAK7b5KuvvoKfnx+eeuopHDlyhF8HYwwHDhzAxIkT4eDggEGDBuH27dsYMWIEACAgIAAi\nkQi7d+9GTU0N5s2bB4lEAicnJ4wYMaLJL8K1a9fw2muvwcPDA6+//joyMjIAcD+LN2/ejNWrV0Mk\nEuGvv/5qsKxCocCuXbswatQoBAYGYtOmTairqwMASKVSuLm5YePGjejduzfmzJmD5cuXY/r06Vi4\ncCFcXFzw448/oqSkBKtWrYKvry+mTJmCEydO8OtftmxZg/enpaVh0qRJ6NGjB5ydnbFkyZImt7tU\nKsXkyZPh6+uLzz//HCUlJfxrQqEQ27dvR2BgIHx8fLB27dom13P+/HmEhYVBLBYjLCwM69evh0Kh\n4F//7LPPEBgYCHt7e/j7+yM1NRUAt7N86qmn4OjoCDc3N3z55Zf8MikpKViwYAE8PDywZMkS5OTk\n8K9JJBJ88803CAsLQ8+ePbF06VJUVVVh6tSpcHFxwRtvvIGKigqNGHfv3o0+ffogNDQUsbGxGq+p\nW7xPP/001q9fr/Gav78/fv3110Y/d3x8PMaMGQNPT0+sW7cOEomE/z9o7G9z4cKFZreTUCjEli1b\nEBAQAD8/P+zcubPB/+XHH38MV1dXjB07FvHx8U3+TW7evImXXnoJzs7OmD9/vkY9Ldm+GzduRFhY\nGDw8PLBs2TLI5XIAQHV1dYPvzv3L/fnnn4iNjcXKlSvx888/QyQSISgoCL/88gtCQkI0YlizZg0m\nTpzY5GcwKoxokEgk7NixY42+JpVKmZubm8Y8gUDAMjMz+elnn32Wvfbaayw/P5/FxcWxXr16sYyM\nDMYYY3v37mW+vr7swIEDTKlUsuTkZHbv3r1G17N+/Xr2wgsvsLKyMqZQKNipU6eajNnT05N98skn\nrKioiK1cuZJJJBL+tdmzZ7MPP/ywyWW//vprNmrUKHb58mV2/fp1FhERwb777jvGGGPHjx9npqam\nbM6cOSwvL4/V1NSw6OhoZmZmxtatW8dqampYTU0NmzVrFps6dSq7desW27NnD3N0dGRZWVmMMdbo\n+6dMmcLWrVvH6urqWFVVFTt79myjsd24cYM5ODiw3bt3s9u3b7Pnn3+eRUVFaWz7UaNGsfT0dHbp\n0iUmEonY9evXG13XpUuX2Llz55hCoWCnT59mnp6e7OjRo4wxxi5fvsz69OnDcnNzGWOMXblyheXl\n5THGGAsJCWH79u1jSqWSlZaWsoSEBMYYY0VFRUwsFrP9+/ezsrIytmLFCjZkyBC+PolEwoKCglhi\nYiJLTk5mdnZ2bNCgQey3335jd+7cYaGhoWzr1q2MMcaysrKYQCBg48ePZzdu3GD79u1jYrGYXbly\nhd+GL774ImOMsV27drHQ0FC+nqSkJObk5MTkcnmDz1xSUsKsra3Ztm3b2J07d9isWbOYmZkZ+/PP\nP5v82zS3ndTbPCwsjKWkpLATJ04wiUTCYmNjGWOM/fDDD8zc3JwtX76cFRcXs+joaDZs2LBG/x6M\nMRYcHMyWLFnCCgsL2RdffMHMzc35/9WWbN+AgAB2/vx5du3aNY3vbXPfHYlEwn/+ZcuWsZkzZ/Kv\n1dbWMkdHR5aens7PCwwMZHv37m3yMxgTSvYP8PT0ZLa2tszBwYEv//3vfxljXPJrLtmXl5czFxcX\nVl1dzb/++uuvs88//5wxxtjUqVPZ2rVrG633wWS/bt06NmbMGJaWltZsvAkJCczFxUVjnqurK5+U\nZs+ezT744IMmlx86dCg7ffo0P71v3z42btw4/vMKBAKWk5PDvx4dHc28vb35aYVCwZycnNjVq1f5\neTNmzGBr1qxp9P2MMTZp0iT25ptvsvz8/GY/25o1a9j06dP56YyMDObk5MSUSiVjjNtme/bs4V9/\n8skn2TfffNPsOtXef/999uqrrzLGGEtOTmaenp5MKpXy61YbOHAgW716NSspKdGY/91337H58+fz\n0wqFgvXo0YMVFBQwxrikot4GjDE2evRoNmnSJH76008/5Xdc6mSvTkKMMTZ9+nS2evVqxphmsq+p\nqWFisZjfqS1ZsoS98sorjX7Gn3/+mQ0fPpyfzszM1Kinsb9Nc9uJMW6bb9q0iZ9+9913+dd/+OEH\nJhaL+W14584dZmZmxiorKxusNz8/n1laWrKamhp+nru7O5/sW7J9v/zyS/71l19+mb399tuMsea/\nO/cn+/u3q9qCBQvY+++/zxjjGgFisZjV1dU1u42MBXXjPEAgEODXX39FSUkJX+bOnduiZU+dOoXC\nwkL06tULYrEYYrEYmzdvxqlTpwBwXRJDhw5t0brmzp2LiIgIjB8/HgMGDMCmTZsafd/p06cxcOBA\njXkhISF8nc2pqqrCmTNn8PTTT/Pxzp49G2fOnOHf07NnT7i7u2ssFxoayj9PT09HbW0tevfuzc8L\nDg7GyZMnG30/AKxduxbV1dXo378/xo4dq9Htc78zZ84gODiYn/bx8YFCoeC7WAAgMDCQf+7i4oLc\n3NxG15Wbm4sFCxbA398fdnZ2WLt2LVJSUgBw3SCffvop3nnnHbi6uuKjjz5CdXU1AGDr1q1ITk6G\nt7c3nnvuOX7E1bFjx7B9+3Z+u3Xr1g1VVVUaIz0CAgI0tuP90z169GgQ6/2fJSgoqNEuEEtLS0yd\nOhUxMTFgjGHnzp1NjtQ5f/68xjq9vb1hb2+v8Z4H/zbNbaeWxOnn5wehkEsrLi4uUCgUKCgoaDQ2\nHx8fWFpa8vPu/z9uyfZt6m/f0u9OY6KiovDTTz8BAGJiYvD888/DzMysxct3ZpTs9SgsLAzdu3dH\nQUEBv6MoLy/n+1NHjhzZZBIWCoUafZ/W1tZ49913kZmZic2bN2Px4sVIS0trsNywYcOQkJCgMe/S\npUstGgJnY2OD0NBQHD58mI+3tLRUo1/8wQOyAoFA4yBe3759YWFhgatXr/LzLl68qNFP+uBBPw8P\nD2zYsAH5+fmYOnUqpk+f3uiww6FDh+LixYv8dEZGBkxMTODn59fkZ2pqNMcnn3wCuVyOP/74A2Vl\nZfjnP/+pUeeMGTMQHx+Ps2fP4siRI/jhhx8AcMlr69atyMvLw4ABAzBv3jwAwKhRozBr1iyNRkFl\nZSWmTJnSZGxMy8HHxMRE/nlCQgKGDBnS6PuioqKwfft2HDt2DNbW1g0SttrgwYORlJTET9+4cQNl\nZWUa73nwb6NtO+kSZ3MGDRqE69evo6amRmNdarpu39Z8d0xNTRv8TR577DGYm5sjLi4OO3bsMIoh\nry1Fyb4R2r6UTXFwcMCwYcPw3nvv4ebNm1Aqlbh8+TKfsKZNm4Zvv/0Whw4dgkKhQEpKCoqLiwFw\nreH7E9vvv/+O69evQ6VSwcbGBubm5hqtILXAwECYm5tj5cqVKCoqwueffw5TU1O+1aPts8ycORMf\nffQREhISoFKpkJubq3FQ+UEPrs/U1BRPP/00oqOjkZubi/379yM2NrbZg1rbtm1DYWEhGGOwsbGB\nra1to++bMGECDh8+jL179yI3NxfR0dGIjIzkW46NxdbU571z5w4cHR3h5OQEqVSKrVu38juGixcv\n4ty5c5DL5bCysoKpqSlEIhHkcjm2b9/OJ0hra2uIRCIAwNSpU7F3717s378fVVVVqKqqwu+//47K\nysomP7c269atQ1ZWFg4cOIAjR45g/Pjxjb4vLCwMAoEAb775JmbNmtXk+saMGYOEhATs2LEDeXl5\n+PjjjxvsvB/U2HZ60ObNm3H58mWcPHkSP//8c5NxNsfZ2Rl+fn6Ijo5GYWEh1qxZo/ELoC3b9+DB\ngy367gQHByMtLQ21tbUa82fOnIlXX30V5ubmrdqRdVaU7BsRGRkJkUjEl8mTJ/OvPdhyfHB648aN\n8PT0xJQpU9C9e3e89NJLKC8vBwA888wzWLVqFdavXw8nJyfMnz8fMpkMAPDmm29i9erVEIvF2L17\nNzIyMjB69GjY29tj/vz5+OSTT5ocgx8bG4vc3FwEBQXh1q1bGiM5tI1dnj9/PubMmYOPPvoIjo6O\nGD16NK5du9bs531w3po1axAQEIDw8HBs3boVu3fvhkQiafL9hw8fRv/+/dGzZ09s27YN3333XaMJ\n3MvLC7t370ZMTAzCw8Ph7++PNWvW6BSb2rJly5CUlAQ3Nzd88cUXePXVV/nXysvL8dJLL8HR0REj\nR47E4MGD8eKLLwLgdkxeXl7w8PDAuXPn8NVXXwEAxGIxDh8+jOPHj6N3797w9fXV2IE05v7XHoxV\nIBDghRdewNixY/HJJ58gJiaG7xpr7HPNmjULf//9Nx9nYxwcHBAbG4tNmzbhsccew8CBA+Hg4MB3\n5TS23sa204PveemllzBjxgy8/PLL+OSTTzB69Ogm19fc9ti9ezeKi4vRv39/XLlyBc8//zz/mq7b\n9/66r1+/3qLvTnh4OHr37g0vLy+NUTgzZ85Eampqs9vWGAlYa5uxhJAOs3XrVvz3v//V6WzQ1NRU\nDBs2DMXFxa0+eclQJ/91JtXV1XB2dkZiYiIeeeSRjg5Hb6hlT4iRqa6uxtq1a1t0+YsDBw6guroa\n165dQ3R0NB5//PEuc0Zua3311VcYPXr0Q5XoAUr2hBiVw4cPo1evXhg4cCAiIyO1vv+3336Dq6sr\nxowZg/79+2PdunVtqv9h31FIJBL88ssvWLlyZUeHonfUjUMIIV0AtewJIaQL6NALoT3sPwkJIcRQ\ndO2U6fCWvXpstDGW6OjoDo+B4u/4OCh+4yvGHDtjret57/BkTwghxPAo2RNCSBdAyb4NIiIiOjqE\nNqH4OxbF33GMOfbW6tChlwKBoNX9T4QQ0lW1JndSy54QQroAugftQ4gxoLwcKCriyr173LStLSAS\nAXZ2XHFyAhwcOjpaQkh7oGRvQIwBdXVcoi0r40ppKTcNAGZmgKkp9ygUAhUV3GvqUlEB1NTUF5kM\nqKoCKiu5UlHBPdbWAgoFIJfXP9rYAN261ReRiFv2/nXfvcvV6+lZX3x8gL59ueLhwb1OCDF+1GcP\nLilXVwMlJY2X4uL6JF1RUZ9kq6u5Ze//CHI595o6KQsEXCva3p5rRdvbc9OAZoJWKjVb3fb2XEvc\nyqq+WFpySVwkqm+l29hw89U7DVPT+tKSz11aCmRnAzdvco/XrwNXrnClpATo3RsYORKYMAEYOrRl\n6yWEGFZrcmeHJ/tffmF8wmQMMDHhkqJYXF+srLhkqE6KSiWgUnFF/Vyp5JJvVVV9olW3ZNUJuqKC\na12ruzbUpbiYq/f+OsViwNGx/tHBgUvCIlF9srW2rm/5qk8GNjXlXrO15RKxuXlHbd22Ky8H0tKA\nw4eBX38FcnKAp58GIiOB4GDulwC1/Alpf50q2V+9ehXTpk3jp2/cuIF//etfeO211+orFwjw7LMM\nAgH4olBwrc37W9YyGZdETUw0i1DIFfVza2suwdrY1CdbdXK+v9XcrRvXX60ujo5c65g0LycH+O03\n4NAhIDmZ2xn07w/4+wMDBnDP+/fntikhxHA6VbK/n0qlgqurK86fP69x8+rO0o1DWqe4GPj7b66k\npACpqcDly9xOtn9/rvX/xBNc9w/tTAnRn06b7I8cOYKPP/64wc22Kdk/fBgDbt/mdgDx8cDRo9xO\nICwMGD0aCA8HgoK44wuEkNZpTe5sl8NtO3fuxAsvvNDoa8uWLeOfR0REdMkz2x4mAgHg7s6VceOA\nf/2L65aTSrnEHxMD3LgBhIRwLf6hQ4FBg4Du3Ts6ckI6L6lUCqlU2qZ16Nyyr62thYWFRYvfX1dX\nB1dXV6SlpaH7A99oatl3TaWlwNmzwJkzXLl0iTvGMnAgV4KDuR1Az54dHSkhnZNBzqCdPn06ysvL\noVQqERoaCl9fX2zevLnFFRw6dAjBwcENEj3puhwcgLFjgY8/Bo4d4/r+T5wAZs7kDsavW8eN8/f0\nBJ57Dvj8c+5XQVYWdwCfEKI7rS37gIAAJCcn45dffsFff/2FNWvW4PHHH8fp06dbVMG0adPw1FNP\nISoqqmHl1LInTWCMG/N//jxXkpOBzEzuRDAPD+CRR4A+fQA/v/pCZwOTrsIgB2jDwsLw559/Yvr0\n6Xj77bcxZMgQ+Pv7IyUlRevKq6qq4OnpiaysLIhEIr0ETLo2mYxr4Wdmcid+paZyJT2dG1br6Qk4\nOwMuLlzp1Ys7MezRR2lIKHl4GOQA7aJFizBw4ECEhIRgyJAhyM7Ohr29fYtWbmNjg6KiIp0CIqQ5\nlpZc4n70UWD8+Pr5KhVw6xZX8vOBvDyunDgBfPcdtzOwsAD69eOW7d27vkgkNDqIPPx0PkDLGINS\nqYSpHs6bp5Y9aS+Mcck/LY1L/BkZwLVrXLlzh9sBPPkkdyxhyBDjPvOZPPwM0o1TWlqKTZs2ITY2\nFgDw1FNPYe7cuS1u3TdbOSV70gnIZMDFi9xlIQ4fBq5eBSIigOHDudFBQUHcZTMI6SwMkuzfeOMN\nqFQqzJo1CwAQExMDgUCAr776qvWRqiunZE86ocJCbvTPuXNAQgKQlMRdYmPgwPouIF9f7pGOA5CO\nYJBk37dvX6SmpsLExAQAoFQq4efnhytXrrQ+UnXllOyJEVCpuJFBCQlcq//atfpuIIEAcHMDXF25\ng8Gurlxxd+fmu7tzvwrUF8ojRB8Mkuzff/99dOvWjR86GRMTg7t37+LTTz9tfaTqyinZEyPGGHfV\n1Nzc+nLnDne5CHW5dYu7jLWHB3cgWCLhRgxJJPW/DhoZqEZIswyS7G1tbVFdXQ3B/5omjDHY2Njw\nFZar78TRCpTsSVdQUcEl/ezs+pKVxf06yMjgWv59+nDFxwfw8gK8vblH9b0PCLlfp70QWpOVU7In\nXZx6yOjVq9x5AzducDuCGze4IhJx1xEaPJgrgwbRcQJiwGR/7tw5HD9+HO+88w5ycnKQn5+PwYMH\ntzpQvnJK9oQ0iTGua+jChfoziS9e5C4ap078gwdzo4WsrTs6WtKeDJLsV6xYgb///htJSUlIT09H\ncXExxowZg4sXL7YpWICSPSG6Uqm4XwDqHcCFC9w9BDw9uesJ9elT/+jry/0KoIPDDx+DXS7h9OnT\nCA4ORmJiIgC0+HIJWiunZE9Im9XWcn3/6nsHq8v169zrjzxSX9zcuFFD6uLsTGcPGyODXC7Bzc0N\ndXV1/HR6ejp69+6te3SEEIOwsKi/JeT9GOOuKJqZWV9SU4EjR7hRQ3fucBeWE4u5pK8urq7c6CFP\nz/pia9sxn43oj9aW/bFjx7Bq1SqkpaVhzJgxOHnyJL7//nuMHDmy7ZVTy56QDqVUAkVF3PWE1OX2\nbeDmTc3CGJfwbW3r7+ncvXv9RefUj15eXDHEcFKZDKirqy9yOferpFu3rnd5C4MdoK2ursahQ4eg\nUqkQGRkJSz3dUJSSPSGdH2Ncoq2srC/l5dyZxnl59Reeu3OHG0mUlcXtFLy9uW4jBwfA3p4rdnbc\nwWQzM8DUtP5Rvf6Kivr137/evDwuwVtYcMuYm3OPdXXcuQ7W1tzOp3t3wNGRq0tdr4MDN+/+IhZz\nMVpbc+s0tuMaBkn2M2fORExMjNZ5rUHJnpCHD2Ncos7K4n4llJVplpoa7iY0cnn9o6Ul92vg/l8O\n6l8LvXpxjyJR40mZMe7uZ0VF3A6ouJirp7S0vpSUcPPvL9XVXJHLASsrwMZGMwZbW26epSW3Q7C0\n5Ir6fepiZ8ftZHr04O6u9r/TkAzKIH32ly9f1piurq5GWlqabpERQroMgaD+fgLtVZ9YzBVfX92X\nVyq5HdD9v1zUvzCqqrgD4DIZV2pquHkFBdxB8YqK+l85d+9y801MuFFQ6l8h6l8iIhF3LER9PET9\n3MOD24kYWpPJfsWKFVi5ciVqamo0bjzSvXt3LFq0yPCREUJIOzAxqT8e0VaMcTuJe/fqjyvI5dzz\nsjLuBLqbN4GTJ7nHnBzu14+9ff1O4NNPuaGz+qa1G+edd97BZ599pv+aQd04hBCiUnG/CHJyuDJy\nJHfQuTkG6bM/deoUAgMDYWtri4MHDyIlJQULFiyAo6OjThU1Wjkle0II0VlrcqdQ2xsWLlwIGxsb\nZGVl4d1334VQKMT8+fNbHSQhhJD2pzXZm5qaQiAQ4IcffsD//d//4Z133kF2dnY7hEYIIURftCZ7\niUSCDz/8ELt27cILL7wApVKpcUZtc6qqqhAVFYXevXujX79+OHv2bJsDJoQQojutyX7btm3w9vbG\njh07YG9vj9zcXLz11lstWnl0dDQ8PDyQkpKClJQUPProo20OmBBCiO4Mej37wMBAxMfHw8rKqvHK\n6QAtIYTozCAHaFvr9u3bkMlkWLhwIUJDQ7Fq1SrIZDJDVUcIIaQZWs+gbS2ZTIZr167hiy++wBNP\nPIGXX34Zu3btwqxZszTet2zZMv55REQEIiIiDBUSIYQYJalUCqlU2qZ1GLQb59FHH0V6ejoA4NCh\nQ9i6dSt27NhRXzl14xBCiM702o1TWVmJ1atX47nnnoOnpyc8PT0xZcoUrF69GpWVlS1aua+vL86d\nOweVSoXff/8dTzzxhE7BEUII0Y8mW/YTJ06Em5sb5syZg759+wLgblyyefNm5ObmYv/+/VpXfu3a\nNcyaNQsymQxPPPEEli9fDpv7LglHLXtCCNGdXi+X4OHhgYyMDFhYWGjMr6mpQZ8+fZCTk9P6SNWV\nU7InhBCd6bUbJygoCG+++SaSkpIgk8kgk8mQmJiIpUuXIigoqM3BEkIIaT9NtuwrKirwzTff4Pz5\n87hw4QIYYxg0aBBCQ0OxcOFCjcset7pyatkTQojODHZbQkOhZE8IIbprt5Oqfvjhh9YsRgghpIO0\nqmXv7u6OW7dutb1yatkTQojO9HoP2gEDBjS50N27d3WqhBBCSMdqMtnfvXsXsbGxEIvFDV4bMmSI\nQYMihBCiX00m+6effhqVlZWNDrMMDw83aFCEEEL0i0bjEEKIkelUlzgmhBDSeVCyJ4SQLoCSPSGE\ndAGU7AkhpAvQmuyPHTuGUaNGwcHBASKRCCKRCHZ2du0RGyGEED3ROhonJCQEX3/9NcLCwiAU6veH\nAI3GIYQQ3RlkNI65uTmCg4P1nugJIYS0H603HB8+fDgmTpyI5557Dg4ODgC4vcqkSZMMHhwhhBD9\n0JrsCwoK4OzsjFOnTmnMp2TfeShVStQqa1GrqOUfK+sqUV5bjoq6ClTUVqBaXg2BQAATgQlMhCYQ\nCoSwMrVCD5se6GHTA91tusPazLqjPwohxEC6zBm0jDHUKmshU8i4uiHgYzARmMDMxAxmQjMIBAL+\n/eqEqU6adco6KFQKyJVyKFQKKJmy/rP8b30KlQI1ihrUyGv4R5lChhqF5qNcKYeSKaFUKaFQKaBQ\nKVAlr0JFbQUq6ypRWVeJGkUNH4uaiqn49ytUCshVcjDGYGFqAQsTC1iYWsDcxBwicxFEFiL+0drM\nGowxqJgKSqaEiqlQLa9GYVUh7lbdRWF1IUwEJugl6gUPew+N0q97PwzoMQA25vX3DyaEdBy93rxk\n1apVePvtt7Fo0aJGK1q3bl3ronxgPd9f+h4lNSUokZWguKYYZbVlkClkXGL8X6KUq+RQMZVGMRWa\nwtLUki/mJuZQqBT8surl1Ymzsq4SpkJTWJjW31OXMQYGBqVKCbmKS+CmQlOYCc1Qq6yFlakV7Czs\nYGdhB5GFCOYm5jATmnHvMTGDUCCEAAIw1G9CdYvZysxK8/G+55amljAzMeNb2SYCE5gKTWFrbqtR\nLE0tIRQI+W2lXr86BlOhKUyE3LJtpd653am4g5yyHK6U5yC7NBuX715GemE6PB08EegcCP8e/vB0\n8ISbnRvc7NzgKnLV2K6EEMPS6yWO+/XrBwAIDg5usGJ14tGHM7fOwNHKEWJLMdzs3GBvYa+RFO9P\njEKBkEuwAkGDxF6rqIW5ibnGDsDS1BK25rYQWYhgY2YDMxOzZmNhjEGukkOulMPC1EIvSdRYCAQC\niCxE6GPRB3269Wnwulwpx5WiK0jKT0LK3RQkX03G7fLbuF1+G3cq7kBsJYbEQQIvBy94ib3g5eCF\nR8SP4BHHR+Bu5w4ToUkHfCpCiJpBu3EkEgns7OxgYmICMzMznD9/XrNyGnr5UFAxFfIr85FVkoWs\n0iz+8UbJDWSWZKKwqhASBwl8HH3Qr3s/+Pf0h39Pf/Tt1hfmJuYdHT4hRscg96DNyspCTEwM4uPj\nIZPJ+Ir++usvrSv38vLCpUuX4OjoqLeAifGpkdfgRskNZBRnIK0wDSkFKUgpSEFWaRZ8HH0woMcA\n9O/Rn38tskhsAAAgAElEQVT0dPDku68IIQ0ZJNmPHz8eYWFhGDVqFMzMzPiKgoODta7cy8sLFy9e\nhJOTk94CJg8PmUKGtMI0XL57GX/f/Zt7LPgbpbJS9OnWB32c/le69UG/7v3olwAh/2OQZO/v74+U\nlJRWBeTt7Q2RSAQvLy/MmTMHzzzzjGbllOxJI0plpbhadBVX73HlStEVpBWmIbs0G76OvghwDkBA\nzwAEOgciyDkITtaNNyYIeVjpNdkXFxeDMYbvv/8eNTU1mDVrlsYtCpvqmrlfXl4eXFxckJ6ejsjI\nSJw6dQrOzs4aAUdHR/PTERERiIiI0OkDkK6jRl6D1MJUpBSkILkgGUn5SUjMS4TYSoyBLgMx0Hkg\n+nXvB18nX/g4+tB5A+ShIZVKIZVK+enly5frL9lLJJJmR91kZWXpVNHixYvx6KOPYv78+fWVU8ue\ntJGKqZBZnImEvAQk5CfgatFVXLt3DTdKbqC7TXf4OvrCS+wFib2EGy0k9oKHvQecbZ2pS4gYLYN0\n48hkMlhaWmqd96Dq6moolUqIRCIUFhYiIiICsbGxcHd3b1PAhLSEUqXErfJbyLiXgazSLGSXZiO7\nNBtZpVnIKcvB3aq7EJmL4GzrDGdbZ/QS9YK32BveYm94OXjBW+yNXqJeNGSUdEoGSfYDBw5EQkKC\n1nkPysrKwrPPPgsAcHJywowZMzBnzpw2B0yIPqiYCsU1xcivzEdBZQFuld9qMGz0XvU9uNu78zsA\nLwcvuIhc0NOmJ3rY9EBP257obt1d6/kbhOibXk+qysvLw507d1BdXY2EhAQwxiAQCHD37l1YWGg/\nW9LLywtJSUk6BUNIexEKhOhm3Q3drLuhf4/+jb5HppBxvwbU5w+UZuFy4WUUVBbgbtVdFFQVoKi6\nCLbmtvy6ull3Qw+bHnAVufJnF7vauaKXqBecrJzolwLpME227Lds2YIff/wRFy9eREhICD/f09MT\nL774Ip544om2V04te2LkVEyFUlkpiqqL+FJQWYDcilzcLr/NP+ZX5qNUVgqxpZi/+Jy7vTt3lrH4\nEfg4+sBb7I1u1t30eoY6eTgZpBtnz549mDx5cpsCa7JySvakC1GoFLhXfQ+F1YUoqCzAzbKbyCzJ\nRGZxJq4XX8eNkhuoklehp01P9LTtCWdbZ7jYusDT3hMSB+4As6eDJ1xsXegXQhen12QfExODmTNn\n4ssvv9Roaai7cxYvXty2aEHJnpAHyRQyFFQWIL8yH/mV+bhTcQc3y27iZtlNZJdm42bpTRTXFMPd\n3r3BdYjUQ05tzW07+mMQA9Nrn311dTUAoKKign5WEtJOLE0t4engCU8HzybfI1PIcLP0Jn9A+Ubp\nDZzPPY+M4gxkFmfCwdIBvk6+8LSvvzKpm50b3O24g832lvbt+IlIZ6G1G6empgZWVlaGqZxa9oTo\nlYqpcLv8NjLuZeBW+S3+yqS3y28jpywHN0puwMrMCj6OPvxxAjeRm8ZOwcHSgRp4nZxB+ux9fHzQ\no0cPjBgxAsOHD8ewYcNgb6+flgEle0LaF2MMBVUFuF58nT9OwB9MLuceFSoFPB244wTq4wX2Fvb8\nDX7MTMxgbmIOOws72FvYw8HSAQ6WDrCzsIOp0JS/DDkxHIMkewC4efMmTp06hVOnTuGPP/6AWCzW\ny7BKSvaEdD5lsjLuOEHp/44TlN1EeW05f68HuUqOOmUdKmorUCor5Ut5bTkUKgUYGH9zHTOhGSxM\nLfj7S6jvpnb/o/o+FNZm1vy9LMxNzFFRV8Hf2KikpoS/sVGtopa/j4WSKTVuKGQqNIXYUoxeol78\nkFcXW5dGz5Y2NzHni4WJBazMrOBg6cDfX8PB0qHTHgg3SLK/ffs24uLiEBcXh6SkJDg6OmL48OF4\n99132xQsQMmekIeRiqn4223WKev4+yKrb8tZp6zj75esfv7grTxrlbWws7CDg6UDxJZiiK3E/I2N\n7r85kVAg1LhVaJ2yDiWyEtypuIPc8lzcqbiDvMo8KFQKjRgZGORKOR9fnbIO1fJqjbvmldeWw9bc\nFvYW9rC35H7BqH/JqKcdLLjn1mbWfExWplawMbfhdzSGOOnOIMleKBRi0KBBePfddzFhwgS9/jyj\nZE8I6ayUKiXKa8tRVluGMlmZxqP614z6ufoWquodWkVtBfIq81BQWQBHK0e42rnC094Tvk6+8HX8\nX3HyhYutS6tyqkGSfXJyMk6ePImTJ08iJycHvr6+GDFiBObNm6dzgA0qp2RPCHmIKVVKFFQVILc8\nF9ml2cgozuDKvQxcu3cNZbVl/LWZXGxd4CJywZKwJfAWeze7XoP12VdUVOD06dOIi4vDtm3bAAA5\nOTk6VdRo5ZTsCSFdWI28BvmV+cirzENeRR7uVNzB5H6T0UvUq9nlDJLsQ0JCIJPJMGTIEH5Ejqdn\n02OAdaqckj0hhOjMIMn+7t276NGjR5sCa7JySvaEEKIzg3XjGAole0II0V1rcqfQQLEQQgjpRCjZ\nE0JIF6A12QcHB2PDhg0oKSlpj3gIIYQYgNZkv3PnTuTm5mLQoEGYNm0aDh8+TP3shBBiZFp8gFal\nUuHgwYNYuHAhhEIhFixYgMWLF7fpiph0gJYQQnRnsAO0ycnJWLx4Md566y1MnjwZu3fvhkAgwJgx\nY1oVKCGEkPbV5M1L1IKDg2Fvb4958+Zh1apV/M3GH3vsMVy4cMHgARJCCGk7rd04N27cgLe35nUa\nsrKy4OXl1aIKlEolQkJC4ObmhgMHDmhWTt04hBCiM4N040yZMqXBvEmTJrW4gq+//hr9+vWjmxkQ\nQkgHarIbJz09HWlpaSgtLcXevXv5vUhhYSFEIlGLVn779m388ccfeP/997FmzRr9REwIIURnTSb7\nq1ev4sCBAygrK9PofvH09MT69etbtPJ//vOf+OKLL1BeXt72SAkhhLRak8l+4sSJmDhxIuLj4xEW\nFqbzig8ePIgePXogKCgIUqm0yfctW7aMfx4REYGIiAid6yKEkIeZVCptNo+2RJMHaFetWoW3334b\nixYtariQQIB169Y1u+L33nsPMTExMDU1hUwmQ3l5OSZPnoytW7dqrIcO0BJCiG70etXLAwcOIDIy\nElu2bOEPrqrfKhAIEBUV1eJKTpw4gdWrV9NoHEII0YPW5M4mu3EiIyMBALNnz9aYL5PJYGlp2arg\nCCGEdAytQy9feOEFlJeXQ6lUIjQ0FL1798bmzZt1qiQ8PBy//fZbq4MkhBDSNlqTfWpqKuzs7LBv\n3z4EBwfj2rVr2LRpU3vERgghRE+0Jntra2tUV1cjJiYGL774IiwtLVFRUdEesRFCCNETrcl+0aJF\nGDhwIEQiEYYMGYLs7GzY29u3R2yEEEL0ROd70DLGoFQqYWqq9Rpq2iun0TiEEKIzvY7GUZPL5YiP\nj0d8fDxkMhlf0UcffdS6KAkhhLQ7rcl+0aJFyM7ORnh4OGxtbdsjJkIIIXqmtRunX79+uHz5MoRC\n/d+bnLpxCCFEdwa5xPHIkSNx/PjxVgdFCCGk47WoZX/lyhW4urrCwcGBW0ggQEpKStsrp5Y9IYTo\nTK/XxlHLzs5udL5EItGpokYrp2RPCCE6M0g3jkQigYWFBU6fPg2JRAIbGxtK0IQQYmS0JvvvvvsO\n06dPx/LlywEAdXV1ePHFFw0eGCGEEP3RmuxjYmJw5MgR2NjYAABcXV3pcgmEEGJktCZ7e3t7jWGX\nOTk5cHNzM2hQhBBC9Etrso+KisKMGTNQWlqK5cuXY/z48Zg3b157xEYIIURPWnRtnOzsbOzZswcq\nlQrTpk2Du7u7fiqn0TiEEKIzgwy9NCRK9oQQoju9XgjN1ta2yVsJCgQClJeX6xYdIYSQDtNksq+s\nrAQArFixAjKZDHPmzAEAbNmyBRYWFu0THSGEEL3Q2o3Tt29fpKen8618lUrFX0KhzZVTNw4hhOjM\nIGfQDh06FKtXr8a9e/dQVFSEtWvXYujQoa0OkhBCSPvTmuz/9a9/4datWxgyZAiGDh2KnJwcfPLJ\nJ1pXLJPJEBoaisDAQDz22GNYu3atXgImhBCiO4OOxqmuroa1tTVqa2sRHByM/fv3w8fHp75y6sYh\nhBCdGaQbpy2sra0BcAd7FQoFHdglhJAOYtBkr1KpEBAQgJ49e+LVV1/V28lYhBBCdKP1HrRtIRQK\nkZycjOzsbIwbNw5Dhw5FUFCQxnuWLVvGP4+IiEBERIQhQyKEEKMjlUohlUrbtA6tffZyuRzx8fGI\nj4+HTCbjFhII8NFHH+lU0ZtvvgkfHx8sWLCgvnLqsyeEEJ3p9QxatUWLFiE7Oxvh4eGwtbVt8YqL\niopgamoKBwcH3Lt3D0eOHMGSJUt0Co4QQoh+aE32cXFxuHz5ssZljlsiLy8PUVFRUCqVcHZ2xptv\nvgkXF5dWB0oIIaT1tHbjvPLKK5g0aRIef/xx/VdO3TiEEKIzvV71csCAAQC4ETXp6elwdXWFg4MD\nX1FKSkobw6VkTwghraHXZJ+dnd1sBRKJROcAm1sfIYSQljHY9eyLiopw+PBhCAQCPPnkk3Bycmp1\nkBqVU7InhBCdGeQM2u3btyMsLAzx8fE4c+YMwsLCsH379lYHSQghpP1pbdkHBgYiNjYWzs7OAICC\nggI8+eSTSEpKanvl1LInhBCdGaRl7+joiJqaGn66pqYGjo6OukdHCCGkw2gdZ9+tWzcEBwdj+PDh\nYIzh1KlTGD16NBYtWgSBQIB169a1R5yEEELaQGs3zpYtWxou9L+fEAKBAFFRUa2vnLpxCCFEZwYb\njWMolOwJIUR3Brk2TnZ2Nr799lscPnwYJSUlfEU3btxoXZSEEELandYDtNHR0QgKCoJCocC+ffsw\nbtw4vPTSS+0RGyGEED3R2o0TFBSExMREBAQE4OLFiwCAkJAQJCcnt71y6sYhhBCdGaQbx8rKCkql\nEuHh4VixYgW8vLx0utQxIYSQjqe1ZX/hwgX07dsXNTU1+M9//oPc3FwsWrQI/v7+ba+cWvaEEKIz\nvY/GUSqVeOedd/DFF1+0ObhGK6dkTwghOtP7GbQmJiaIi4tDRUVFmwIjhBDSsbT22Q8dOhSRkZGY\nMmUKf6cpgUCASZMmGTw4Qggh+qE12RcXF0MikeDSpUsa8ynZE0KI8aAzaAkhxMgYZOil+oJnAPjr\n4Xh5eWHcuHHo06dP6yIlhBDSrrSeQWtiYoLTp0/DyckJTk5OOHPmDFJSUjBv3rxmr3h569YtjBw5\nEn5+foiIiMBPP/2k18AJIYS0nNZunMGDByM2Npa/hn1xcTHGjh2Lo0ePYvTo0Th//nyjy+Xn5yM/\nPx+BgYEoKirC4MGDkZycDJFIVF85deMQQojODHLzEhMTE5SXl/PT5eXlEAgEsLe3h1wub3I5Z2dn\nBAYGAuCuie/n58dfboEQQkj70tpnHx0djZEjR2LAgAEQCAT4+++/sX79elRVVeGJJ55oUSXXr19H\namoqBg8e3OaACSGE6K5Fo3EUCgXOnj0LgUCA0NBQmJpq3UfwKioqEBERgY8++ggTJkzQrFwgQHR0\nND8dERGBiIiIlkdPCCFdgFQqhVQq5aeXL1/euW5eIpfL8fTTT2PcuHF44403GlZOffaEEKKzTnWn\nKsYYoqKi0K1bN6xZs6bxyinZE0KIzjpVsj916hRGjBgBf39/fpz+ypUrMXbs2PrKKdkTQojODJbs\n8/Ly8Ndff2HGjBkoLCxEZWUlvLy8Wh0oXzkle0II0ZlBhl5+9913mD59OpYvXw4AqKurw4svvti6\nCAkhhHQIrck+JiYGR44cgY2NDQDA1dWVLnlMCCFGRmuyt7e3h1BY/7acnBy4ubkZNChCCCH6pTXZ\nR0VFYcaMGSgtLcXy5csxfvx4zJs3rz1iI4QQoictOkCbnZ2NPXv2QKVSYdq0aXB3d9dP5XSAlhBC\ndNaphl62qHJK9oQQojODXM9eJBI1WLG3tzfGjRuHV155hfrvCSHECGht2a9YsQJ1dXWYPXs2AG50\nTmVlJXr27InMzExs2LCh9ZVTy56Qdufo6IiSkpKODoO0gFgsRnFxcYP5BunGefTRR5Gens5Pq1Qq\n+Pn5IS0tDQMGDMDly5d1qlCjckr2hLQ7+t4Zj6b+VgY5qSosLAxff/01iouLUVxcjPXr1yM0NBQC\ngUCnq18SQgjpOFpb9rdv38bKlStx9OhRAMCYMWPwzjvvwMnJCZcuXcKwYcNaXzm1MAhpd/S9Mx76\nbNnTaBxCuhj63hkPfSZ7rf0wdXV1OH78OI4cOaJxUGfz5s06VUQIIaTjaO2z/+CDD3DgwAHs27cP\nAQEBSEtLQ8+ePdsjNkIIaVR8fDyGDh3a0WFo2LJlC4YPH97RYTRJa7L/888/8e9//xtWVlZ4/fXX\ncejQIfz555/tERshhDTq999/x9NPP91gvkKh6IBojIPWZG9iYgKBQICgoCAcPXoUZWVlqK6ubo/Y\nCCFdyA8//IBnnnmGn/b19cXUqVP5aXd3d6SkpAAADh06hHHjxgEAhEIhtm7diqCgIPTt2xcAkJKS\nggULFsDDwwNLlixBTk4Ovx6JRIL169dj8ODBeOSRR7Bx40bI5XL+9eaW/eyzz+Dj4wMnJyfMmDED\nJ0+eBACkp6dj4cKFiI+Ph0gkgqOjIwDuJk5PPfUUHB0d4ebmhi+//FLfm63lmBbff/89u3fvHrtw\n4QILDw9nvr6+bOfOndoWa5EWVE8I0bPO+r27ceMGc3BwYIwxlpubyzw9PZm7uztjjLHMzEwmFosZ\nY4zduXOHubq68ssJBAI2fPhwlpSUxGQyGSsqKmJisZjt37+flZWVsRUrVrAhQ4bw7/f09GR9+vRh\nJ0+eZElJSSwoKIht3LiRMca0Lrt7926Wl5fHqqur2Zo1a5ibmxv/2pYtW9iwYcM0PlNISAjbt28f\nUyqVrLS0lCUkJOi0TZr6W7Xmb9hsy16lUvF7qZCQEEilUly7dg3PP/98e+yHCCFdiJeXF0QiERIT\nExEXF4cnn3wSvXr1wtWrV3HixAmMGDECAPDHH3/gqaee0lh2/vz5CAgIgIWFBfbu3YspU6ZgwoQJ\nsLOzw9KlS3H9+nXcvXsXADeSZerUqRg2bBgCAgKwYMECHDx4EACaXLagoAAAMGXKFDg7O8PKygpv\nvPEGBAIBLl26BACNjo5RqVTIzMxEeXk57O3tERQUZLDtp02zyV4oFGLVqlU0TIuQLkYgaHtpjfDw\ncEilUpw8eRLh4eEIDw/HiRMnEBcXh/DwcABcsld34aiFhobyz48dO4bt27dDLBZDLBajW7duqKqq\nQlxcHP+ewMBA/nlQUBDi4+ObXVbdXfPbb79h0qRJ6NWrFxwdHZGXl8d3LTVm69atSE5Ohre3N557\n7jkkJye3bsPogdY++4kTJ2Lp0qW4fPkyfxZtY9dqIIQ8PBhre2mN8PBwHD9+HCdPnkRERASf/E+c\nOIHw8HDI5XLExcVh9OjRGsvdfzb/qFGjMGvWLJSUlPClsrISU6ZM4d+TmJjIP09ISMCQIUO0LltV\nVYX58+cjKioKV65cQXFxMVxdXfnGsImJSYOGsZ+fH7Zu3Yq8vDwMGDCgQ+8FojXZb968Gb/88gsi\nIyMRHBzMF0II0Td1spfJZOjVqxeGDRuG2NhYFBcXIygoCKdOnYK/vz9sbW2bXMfzzz+PvXv3Yv/+\n/aiqqkJVVRV+//13VFZWAuC6W/bs2YPTp08jJSUF3333HcaPHw8AmDp1apPLVlRUoLKyEi4uLlCp\nVFi5ciXu3LnD1xscHIyMjAy+Hrlcju3bt6OsrAwAYG1tDZFIZKhNp5XWZJ+dnY2srKwGpSXmzJmD\nnj17YsCAAW0OlBDy8PP19YVIJOLHq9vZ2eGRRx7B0KFDIRAIGh1yKXigz8jBwQGHDx/G8ePH0bt3\nb/j6+mLr1q38+wQCAV555RUsXrwYEydOxNy5c/mr+orF4kaXBQBnZ2esXLkSM2fOREBAAOrq6jQu\nF9OvXz9MnDgRfn5+6NGjBwBg27Zt8PLygoeHB86dO4evvvrKINutJbReLqG2thb79+9HXFwcNmzY\ngIyMDFy9epXfEzbn5MmTsLW1xaxZs/D33383rJxO2yak3Rnz987Pzw979uzhh1i2hpeXFzZt2oRR\no0bpMTLDaNerXkZHRyMhIQFSqRQA0KtXL7z//vstWvnw4cMhFot1CogQQhojl8sRFRXVpkTflWlN\n9sePH8eqVatgbm4OALCxsTHaVgEhxHiZmZlh6dKlHR2G0dJ6IbQ+ffrwBxgA4OzZs3odK7ps2TL+\neUREBCIiIvS2bkIIeVBLjzl2JlKplO9daS2tffYXLlzgh172798fBQUFiImJafGInOzsbERGRlKf\nPSGdBH3vjEe7XuJ40KBBOH78OC5dugSlUonBgwfrVAEhhJCOp7XPPjIyEj/99BP69u2rc6KfPn06\nhgwZgmvXrsHd3R0//PBDqwMlhBDSelq7caRSKX7++Wf88ccfCAkJwfTp0zF+/HhYWlq2vXL6OUlI\nu6PvnfHokNsSKhQKHD9+HN9//z1iY2NRXl6uU0WNVk7/dIS0O/reGY927bMHgJqaGvz222/YtWsX\nEhISEBUVpVMlhBBCOpbWlv3UqVNx7tw5jB07FtOmTUN4eDiEQq1d/S2rnFoYhLQ7+t61r2XLliEz\nMxMxMTE6L9uuLfs5c+Zgx44dMDExAcBdAmHnzp3YsGGDThURQkhHUSgUGlfG7Iq0NtHHjh2L5ORk\nvPXWW/D09MSHH35IpysTQgxCIpFg48aNCAsLg4eHB5YtW6Zxy0CpVIrJkyfD19cXn3/+OUpKShpd\nT3Z2NoRCIXbv3o3+/fvzl0Q+deoUZsyYAS8vLyxfvhxFRUX8MkKhEFu2bEFAQAD8/Pywc+dOjdZz\nc8u+/vrr8PDwQM+ePfHyyy/z162PjY3FypUr8fPPP0MkEvEnpP76668YOXIkHBwc4O3tjZ9++kl/\nG7EpTd3C6sqVKyw6Opr17duXDR8+nK1bt46/RZi+NFM9IcRAOvP3TiKRsICAAHb+/Hl27do1JpFI\n2LFjxxhj9bct3L17N7t9+zZ7/vnnWVRUVKPrycrKYgKBgE2cOJFlZmYymUzGkpOTmZubGzt69Cgr\nLi5mixYtYi+88AK/jEAgYGFhYSwlJYWdOHGCSSQSFhsbyxhjWpfdtm0bKy4uZiUlJeytt95iQ4cO\n5V9btmwZmzlzJj9dV1fHPD092dmzZxljjOXn57PU1NRGP0dTf6vW/A2bXEIgELDIyEh28+ZNfp5E\nItG5gmYr78T/dIQ8rDrz904ikbAvv/ySn3755ZfZ22+/zRhjbM2aNWz69On8axkZGczJyYkplcoG\n61En+7i4OH7ee++9xz799FN+uqioiHXr1o0pFArGGJfzNm3axL/+7rvvsldffbXZZeVyeYO6y8vL\nmY2NDSssLGSMMRYdHc1efPFF/nW5XM5cXFzYtm3bWFVVVbPbQ5/JvslOrL1792LHjh0YMWIExo4d\ni+eee44O6hDSRQiWt/K+gvdh0a3LF/ffMtDZ2RmZmZkAgDNnzuCxxx7jX/Px8YFCoUBqamqT98x4\n8HaFaWlp+OKLL/h5dXV1SEhIwKBBgxrUHRQUhFWrVjW7bGJiIgYNGoQtW7Zg9+7dOH/+PBQKBWpq\navD3339j5MiRDWIyNTXFnj178Pnnn2PRokV49tln8cEHH8DLy0un7aSrJpP9xIkTMXHiRFRWVuLX\nX3/F2rVrUVhYiIULF+LZZ5/FmDFjDBoYIaTjtDZRG9LQoUNx7tw5fjojIwMmJibw8/NrcpkHb1c4\nadIkvP32202+PzExEQMHDgTQ8HaFTS1769YtLF68GAcOHEBQUBCUSiVcXFz4xrGpqWmDhnJYWBj2\n7duH8vJyLF26FEuXLsXu3btbsBVaT+sBWltbW8yYMQMHDx7ErVu3EBQUhM8++8ygQRFCyIOeeeYZ\nHD58GHv37kVubi6io6MRGRnZ4qHgM2fOxLfffosjR46grq4OZWVlDRLs5s2bcfnyZZw8eRI///wz\nf5Om5pYtLCwEYwzOzs6oqKjAe++9h9raWn6dwcHBSEtL4+fdvXsXv/76K6qqqmBiYgJLS8t2uV2h\nTgPmHR0d8dJLL+Gvv/4yVDyEEMITCAT87QS9vb2xe/duxMTEIDw8HP7+/lizZk2zy96vX79++PHH\nH7Fr1y64ublhwIABOHz4sMZ7XnrpJcyYMQMvv/wyPvnkE34UT3PLDhw4EP/3f/+HUaNGYcSIEejf\nvz/c3d35dYaHh6N3797w8vJCSEgIVCoV1q5dC1dXV/Tt2xfFxcVYvny5XrZXc1p8uQSDVE4ndxDS\n7uh71zihUIjr16/D29u7o0PhtettCQkhhBg/SvaEEIKG3T4Pm659/jAhhPyPUqns6BAMilr2hBDS\nBVCyJ4SQLoCSPSGEdAGU7AkhpAugA7SEdDFisfihH3nysBCLxXpbF51URQghRqbTnVQVFxeHRx99\nFL6+vvj3v/9tyKo6hFQq7egQ2oTi71gUf8cx5thby6DJ/vXXX8e3336LY8eOYcOGDRp3dnkYGPs/\nDMXfsSj+jmPMsbeWwZJ9WVkZAGDEiBHw9PTEmDFjNC5PSgghpP0YLNlfuHBB4161/fr1w9mzZw1V\nHSGEkGYY7ADtsWPHsGnTJuzYsQMAsHHjRuTm5uJf//pXfeU0IoAQQlpF19RtsKGXgwYNwltvvcVP\np6amYuzYsRrvoZE4hBDSPgzWjWNvbw+AG5GTnZ2No0ePatwPkhBCSPsx6ElVX331FV5++WXI5XK8\n9tpr6NatmyGrI4QQ0gSDDr0MDw9Heno6rl+/jtdee42fb2zj7+fMmYOePXtq3MG+oqICEyZMgIeH\nB39j9s7q1q1bGDlyJPz8/BAREYGffvoJgPF8BplMhtDQUAQGBuKxxx7D2rVrARhP/GpKpRJBQUGI\njP4UrkIAAAkhSURBVIwEYFzxSyQS+Pv7IygoCIMHDwZgPPFXVVUhKioKvXv3Rr9+/XDu3Dmjif3q\n1asICgrii729PdatW4fKykqd4++Qa+MY2/j7f/zjH4iNjdWY980338DDwwMZGRlwc3PDxo0bOyg6\n7czMzLB27Vqkpqbil19+wQcffICKigqj+QyWlpY4fvw4kpKScOLECWzatAkZGRlGE7/a119/jX79\n+vEDE4wpfoFAAKlUisTERJw/fx6A8cQfHR0NDw8PpKSkICUlBX379jWa2Pv06YPExEQkJibi0qVL\nsLa2xrPPPov//Oc/Osff7sneGMffDx8+vME1Ks6fP4+5c+fCwsICc+bM6dSfwdnZGYGBgQCAbt26\nwc/PDxcuXDCqz2BtbQ0AqKyshEKhgIWFhVHFf/v2bfzxxx+YN28ePzDBmOIHGg6oMJb4jx07hvfe\new+WlpYwNTWFvb290cR+v2PHjsHHxwfu7u6ti5+1s6NHj7Jp06bx09988w374IMP2jsMnWVlZbH+\n/fvz0x4eHqympoYxxlhVVRXz8PDoqNB0kpGRwby8vFhFRYVRfQalUsn8/f2ZiYkJ+/e//80YM66/\nwZQpU1hCQgKTSqVs/PjxjDHjit/Ly4v5+/uzCRMmsF9//ZUxZhzx37p1i/Xp04dFRUWxwYMHs88+\n+4xVV1cbRewP+sc//sE2bNjAGGvdtqdLHLcSM8JhoxUVFXj++eexdu1a2NraGtVnEAqFSE5OxvXr\n1/Gf//wHiYmJRhP/wYMH0aNHDwQFBWnEbCzxA8Dp06eRnJyMlStXYvHixcjPzzeK+GUyGa5du4bJ\nkydDKpUiNTUVu3btMorY71dXV4cDBw7gueeeA9C6/512T/aDBg3ClStX+OnU1FQ89thj7R1Gmw0a\nNAjp6ekAgPT0dAwaNKiDI2qeXC7H5MmTMXPmTEyYMAGA8X0GgDtQOG7cOJw7d85o4j9z5gx+++03\neHl5Yfr06fjrr78wc+ZMo4kfAFxcXAAAjz76KJ555hkcOHDAKOL38fFBnz59EBkZCSsrK0yfPh2x\nsbFGEfv9Dh06hODgYHTv3h1A67677Z7sH5bx96Ghodi8eTNqamqwefPmTr3DYoxh7ty56N+/P954\n4w1+vrF8hqKiIpSWlgIA7t27hyNHjmDChAlGE/+KFStw69YtZGVlYefOnRg1ahRiYmKMJv7q6mpU\nVFQAAAoLC3H48GGMHTvWaOL39fXFuXPnoFKp8Pvvv+OJJ54wmtjVduzYgenTp/PTrYrfEH1L2kil\nUta3b1/2yCOPsK+//rojQtDJtGnTmIuLCzM3N2dubm5s8+bNrLy8nD3zzDPM3d2dTZgwgVVUVHR0\nmE06efIkEwgELCAggAUGBrLAwEB26NAho/kMKSkpLCgoiPn7+7MxY8awH3/8kTHGjCb++0mlUhYZ\nGckYM574b9y4wQICAlhAQAAbNWoU27RpE2PMeOK/evUqCw0NZQEBAWzJkiWssrLSaGJnjLHKykrm\n5OTEysvL+Xmtib9Db15CCCGkfdABWkII6QIo2RNCSBdAyZ4QQroASvaEENIFULInHc7ExARBQUEI\nDg7GW2+9BblcrtPyd+7c4U82SU5OxqFDh/jXDhw4gFWrVuklzqFDh+plPa1la2vbofUT40ajcUiH\nE4lEqKiogFwux8SJE7Fw4UKMHz++VevasmULLl261OmupqpUKmFiYtKmdai3EyGtQS170mmYmZlh\n1KhROHXqFMrLy/Hhhx8iMDAQc+fORWZmJgAgKSkJjz/+OAIDAzFw4EBUVVUhOzsbAwYMgFwux0cf\nfYSff/4ZQUFB2LVrF7Zs2YJFixYB4H4BvP766wgICMA///lPFBQUAABmz56Nt99+G0OGDEFISAiO\nHTvWaHzqlrVUKsXjjz+OadOmoV+/fnj//fcbfX9ERATef/99hISEYN26dbh69SoWLlyI0NBQvPLK\nK7h37x4A4Pvvv8fgwYMRHByMpUuXoq6ujo937ty56Nu3L1asWKG/DU26JEr2pNMoLS3FwYMH8eST\nT2LLli0oKytDwv+3dz+h7MdxHMefE9sclBwUOZBk+ZdvarVdRIo0c5A/SdEcnGlqTSnt4iC1khwk\nQg5yIUSyg9kibXaZE4oTDg4KNbPfYfnm9/tt8fv90u/Xb+/Hbft+Pp/v+7vDa58+W59PIEB9fT1O\npxOIH4jjdDo5PT3F6/Wi1+vV/hkZGbhcLrq6uggGg3R0dHx3zvHExAQFBQWEQiFyc3OZnJxUr4VC\nIfb393G73UmD9f1YBwcHjI2NEQwGWV9f5/r6OmH7y8tLfD4fg4ODDA8P43Q6OTo6ory8nNnZWQDa\n2to4Pj7m5OSEx8dHPB4PAOPj4xgMBs7Oznh+fv6DT1YICXvxD3h6ekJRFFpbW2lpaaG2tpbNzU36\n+vpIS0ujs7MTv99PJBLBZDLhcDiYmpri5eXlp6WRWCyWdJOo7e1tbDYbAP39/WxsbADxUG5vb0ev\n12MymQgEAh/WbDQaKS0tRafTYTab8fl8Cdt1d3ej1Wq5vb3F6/VitVpRFIWZmRkODw8BuLi4oKen\nh4qKCra2ttjd3QVgZ2cHm82GRqNR6xbid33psYRCfEZmZibBYPCn938MbY1Gw8DAAA0NDSwvL1NV\nVfXL+5An+yLIzs4G4rtrRqPRD8d5f76BVqtNOvPOz88H4mv2OTk5CZ/TbrczMjLC0tISbrebUCj0\nYb1C/CqZ2Yt/ksViYXFxkWg0yurqKmazmfT0dM7PzykuLmZ0dBSDwaCu5b8pLCzk7u5Off0+LJub\nm1lYWOD19ZW5uTmsVuuXP8fb/fPy8igqKmJtbY1YLEYkEiEcDgPxtfmSkhLu7+9ZWVlR+zY1Nan1\nzs/Pf3mt4v8mYS/+uvdr4W96e3vJysqipqaGvb09dR3d7XZTWVmJ0WjEYDBgNpu/G8NkMvHw8KD+\nQKvRaNRrdrudq6srFEXh5uaGoaGhhDUkquejNp/pMz09jcfjobq6GkVR8Pv9ALhcLiwWC42NjdTV\n1antHQ4H4XCYsrIydDpd0nsI8Rny10shhEgBMrMXQogUIGEvhBApQMJeCCFSgIS9EEKkAAl7IYRI\nARL2QgiRAr4B62L8f0tq6nMAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 21 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }