{ "metadata": { "name": "", "signature": "sha256:deb2fc17a02fc33140d49e35fa16fb21ebfca223c044ae7d79482e5412dff6ae" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "``anhima.loc`` - Locating variants, samples and genome regions" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import numpy as np\n", "np.random.seed(42)\n", "import random\n", "random.seed(1)\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "import sys\n", "# import anhima\n", "# dev imports\n", "sys.path.insert(0, '..')\n", "%reload_ext autoreload\n", "%autoreload 1\n", "%aimport anhima.loc" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "# simulate non-uniform variant positions\n", "n_variants = 1000\n", "p = 0\n", "pos = []\n", "for i in range(n_variants):\n", " gap = int(np.abs(np.cos(i/100))*100)\n", " p += gap\n", " pos.append(p)\n", "pos = np.array(pos)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant locations\n", "anhima.loc.plot_variant_locator(pos, step=10);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAABUCAYAAAClWsGgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfWmQXtV55nO0ICSBJQuzIxAYmUUGG2Njx2vseE0cJ3Gc\n2ElVKjVTNX9mqmaqpsoLzuY4To0zqcqMncROamKPHWeCsMHsIECA2EGgBQkhIXaEJIQEQkIS3ZK6\n7/w4/fo857nvud/9uj8h09ynSqU+3d9y7z3nvOvzvidUVYUOHTp06NBhophyuC+gQ4cOHTpMDnQK\npUOHDh06DASdQunQoUOHDgNBp1A6dOjQocNA0CmUDh06dOgwEHQKpUOHDh06DATTmv4YQug4xR06\ndOjQoYaqqoL+rtFD+Yu/qFBV8d/WrRW+9700fuKJCj/6URqvWVPhZz9L47vvrnD99Wm8ZEmFu+5K\n40suqfDww2n8D/9QYcuWNP7Wtyrs25fGf/qnFUZH87H9rOPR0Qp//udpvGVLhe9/P42XLq1wxx1p\n/Hd/V+Hll9OY77uq8s8aHq7wzW+m8fr1FX7ykzS+5poK99yT39fmzWn8l39ZYWgojb/+9fy+vv71\n/Lv/5E+a75PHIyP5te7fn1/rnj0Vvv3tNN6+vcJ3v5vGzzxT4V/+JY03bKjw7/+exitWVLjyyjS+\n/fb4LG187bUV7r8/jf/t3+Jn2Pjv/77Ctm1p/Fd/FZ+n95z7He/YUeE730njO++scNNNafz97+fr\nS+f4z/4s/XzgQIVvfCNf6z/4QRrfcEO+fr73vfjsbPzNb1bYuzeNL744zo2Nv/a1/LsvvjifU14D\nBw/ma2B4OL+2vXsr/PVfp/GLL8b1bOPNm/O1//jjFX784zRes6bCZZel8b335vt26dI4zza+8sq4\nDmz8r/8aP9PG3/lOnAsbf+Mb8R4GMce7d1f4279N4wceqHD11Wn8ox/FueL38t7iOR4ZyffOs89W\n+Md/zO+b1/Y//3N+n9/6VrwenlO+z698JZ/TL385/+6vfjWfU762PXvyfbtzZ4W/+Zs03ro137dP\nPZXv23Xronzl58T79o47Ktx8c76eWWb97GdxXdj4hz+M38HzUjnKpKdCYcyaBbz6ahofcwzw4otp\nfMopwObNafzWtwJPPJHG554LPPJIGr/zncBDD6XxhRcCK1ak8QUXAKtXp/HppwNPP53GRx6ZXw8j\nyK0efzzw/PNpPG8e8NJLaTx1KnDwoP9ZisceAxYuTOMHHgDe8540XrEi3oth61bgpJPSeP9+YMaM\n+HNV5dc7PAwccUR67e7dwNFHp/G2bcBxx6Xxpk3A/Plp/NRT8TkZnn4aWLAgjTdvjvNkePHFOI+G\nV18FZs707jpi5sz8mb/5zcDOnWl84on5cz7tNODZZ9NY18QZZ8Rr5s/jeQkBGB0tXw9j7tz8WvRa\nDxwApk/331uJH75hA3DOOWm8bBnwkY+k8b33Au97Xxpv2gScemr6rH374n4B4pxOnw5MGdtpL70U\n79PwzDPpvUDcI+eem8arVuXryVtv7353Gq9dC5x/fhrren3uuXwNvPRS3A8Gbw3wfpoxI96T4eij\n4zo1qFx4y1uAHTvSeNq0OBfjgc5pVeXX9sor+X7ha9c5fvhh4O1vT+MbbgA+85k01jl/8sm4fgFg\nZATYsyd9165dwFFHRTkCxH132mnpvatXA+96VxrfcQfw4Q+n8dKlwCc+kcY33gh86lNpfPPNwCc/\nmV/br/5qGt93X74e16xpXgMqB156KZcD+/YBs2en8chInLc2aK1QZs6MX2R405vyhaQb+rjjgBde\nSONTTokbz7BwIbBxYxq/4x25gnnPe+Lm4fHy5fnr16wpfz9jypR8Qc2bly/6adPiQ/Owa1e8V8Mj\njwCLFqWxTtaBA0kpjI7mC37Hjnzi1q3LF/V99wG/8itpfOedwIc+lMb33AN84ANprMLlUAgThm7o\nefPyZ37CCVGBGk49NQpMgyqUt70tXwNnnQU8+mgan3wysGVLGjcpmKlT87/ptR48WFYoL7wQjQ7D\ngw/mQvrJJ6Py8z5Llf66dcB556XxLbcAH/94Gl99NfBbv5WPP/e5NFbBdsstwMc+lsZ33QV88INp\nvHw5cNFFafzQQ/2tgRdfbF4DKohVobzpTVGQG1ShHH98Lgd03I+C8fapKhTeqwwVorq3nn46GWMj\nI/G+TYiq0tc5vfRS4ItfTOPFi4EvfSmNr7wyn/ObbsoVxF135fta19/q1VHeGdaty2XQ+vXA2Wen\n8aOPxr1l0PX73HNxbxl0DezdmwwiIK53U5a90FqhTJ+eW/HqBfQ79gTA0FAaH3sssH17Gp9zTnxw\nhgsvjA/ecP75uYJp8mD68VDY+gTqk1VVyfp89dX4vYaHHoqemOHWW3PhsHQp8Gu/lsZ33pkLC1UY\nKizUylq3Lrdu2wgT9VB4ISk8D4WfoxoRJ52UK4QFC3IvU40KVSiqcE48MVdYTfA8lJKVtXFjfQOa\ncFEvkv8GAEuWAJ/+dBpfey3wG7+RxnffDbz//Wn86KPxPg3qwe7cmW/u3buBOXPSWIXmyy9HY8qw\nY0fcOwadcxUmap3266GoYXnMMc1r4oQTci/2uOPyfV5VdSVWgr5ueDh5/wq12vm5a2Tg7rtzAX/Z\nZcAXvpDGt90GfPSjafzEE8CZZ8afR0fjnNgcDg9HBWXP9NVXoyy1tTg0FL/bZMiBA1Ee2diMUp4D\nljn2Ghb46lEMD+dyib0rIK459ppfEw/lcEOtEw2vqIfztrdFgerhqKPiQy19NuPZZ/OwEoesVEip\ntXjbbblrunZtbr1u356HsFSgq1WtE6thHF04Kiw8y0Rd3X48FB3rc9TxEUfk1uicObkwOu203KNR\nhfPWt0aBzp/Hwq3pWptCXqpQgLSBNWx44425Alm9OjcaWAEcOBA3vm129WbUY9Xx7t1xrfI9tN3Y\nBhU+Kix6eam9PBQNean3f9xxMVRrUIVy/PH53486KlrIbe9NDVWDWtVr1ybjS+/prrtyb+WWW3JD\nj/fprl3xnu2ZPvxw7i3oe6+9FvjN30xjzyP99V9PYw1n6fravDk3QFQZjgeqkF4TD2U80BCFLj5d\nXGqBzp0btb1hypSy4FdL5+yzYyzcwK615y219VAY6iHce28esmKBXlXxWdhC1IWgobWXX84tUxX2\n/Flt4SmcXiEv3nwqpEubebyYOrVuNPB6OOOMXKFoSIzhhbxMGOsGYaNh//5c8dx+e77B2eI3AW/P\n4fHHU5wdiAYFe6TXXgt89rNpfM01ubDRcJe+X0MhKlxGRvpfE7qu1KhRoT3RkJfu+V4eSxOaFMrW\nrXH9GPbuTcpZQ0AcGRgdjevD9uaaNbkR+LOfAb/3e2l8+eXA7/5uGqtCuf/+3MhcuTLPpyxfnkch\nNL9y8815fkW9I1U4mqMbD3p5PE0YuEJhAaQb/txz87DVO96RJ941Ma9hLX1/U1jrzDNzD0Wvha+z\nyUN5/vkUX1dhoyEpdiXVmmS3GKi71Sq4NMb7wAP5wlSrupd34WHXrlxpqULRzaphz0MB/k79/vnz\n8yT/KadEAe/BS+CasFXPjP+meS1WEmoEqGWr1ucdd+R/Vw/1scfyOdywIY+FqwLR/Mn99+fJWA1x\nag6vBH6Nt476CXnNmRPXlUFzm0qQUYVy7LG5AgKaQ2Cl+3vuuTyywNC9xVGH++8H3vve9DfNfzz+\neHrG+/fHfW4e30svxfu1tbR1a7w/u8YXXoj3Z+M9e+J7bTw6GuUQyxgNeeoa0evVPOrOnXlIdDw4\nbB6KWie9mF7nnx8fAI+V+bVyZRp7iXl+PXsaM2bECTdogpjB79OHx9paNyyzrNQ6fPDBXNmopaHC\n4cEHczaPssXuuy9fOKrMlCzQxhVWC69XDmU8UK9y9uw83KjjJi9UFZoSPYAkfKZPLyd7t2/P8wwM\nFeJAekY6Z8oE4vyFXacZFRZGsM/SsIIZKzwfainu2pULB1VQa9bkyVvNn7RBvyGvWbPyEFWv3Knu\nS1UwmjudNSsnAzGaFM2mTenedS88/ngy7pQZyQnz0dH43ebZPPZYbhRed13ucWoyfvFi4A/+II1/\n+lPg938/jTUcpsqh5G3wM9WweS9Sxni82MPmoZx8cp06/PjjaXzaaXlS9uijc3dZE/P699NPz2mm\n6tEsXFjOmyiFlSeFPZQmpogKbf4c3cxqnSrLRMNPOmmaP9FEmi4UTTqqW98G4/FyPPBGVwtUw1Y6\nZ5pHaQpHqofC4Zcmy5w3oW4wFkRK/2SFUlU5o0/zU3fckSsbDWdpMl+9nSeeyOfPy59wGA+I1isn\n/NUDsoRvE4aG8nUJNHsoEw176p7XkNeb35yHvRlNIS/2UNavz8kqQHofz0tV5ftS50TDXWrkPfNM\nyrdVVVz3HHbT8LkKf2V/6VjzufY9DCVlsEcFRC+JFVAbvGYeyowZ+WLQWpT583MLUum7/UIXT6/E\nPOdwNFzCYKG1e3dZoTDDS4WvWq+smFRoafJVY+Hqpo6M1O9dk60q/FThDA3VGTD6mb1ow22gVOJe\ntShKnlCml9YfMTSBqzUPJVjoAYjXxs8NSM9Fw5AcGnn44dw7uPJK4Ld/O401ubpyZbMHqgaIChON\nvXubXI2QjRvzNfD88/k6K0EZRQxVKIOGhryaygGaFArnJNnY4hA2kHtxOifqrbzySgo/bdkS5Y99\n/6pVeS5DIwiq7C1HynUyuv/UA9X8SRvloKxENWzbYGTkNVIoSr89+eTcYlRqsAfNg5x4Yp7r0MSs\nKjGGCqP589P1HHlkvhFKORRWKBqDZoGiybBt29JCVeGtYRRN3CmvXRNzGtcvsX34WnsVMwF1YTEI\nhdJvLYp6sTqHyvTSEAujrULhkBfnonTeOMTBngsAXH997nFs2ZI8lJGRPBZun2vzo8n8qopri8My\nmgPQuL8SQjw8/3ycD/7MfoWJCm3dR4PG7Nl5iKvJQwGaPST7G3t7/ByV5r9kSZrT/fvj+20OVZhf\nckkezrriitygUMLF5Zfn1GM1QHSPe3Rh3dOaQ/PCWSyzgLon3db7aOuJTkihaM7khBNyCqAHFQhv\nf3vcHAYNY2mBo1bYs0LSGLsyvRhMOVUPxUJL27blG5LBFogKZi1QVGtVLQ+tbVFrVBP4+n6Pu6/J\nvDbCRIWafXYTNOehXuOpp+YeirKCVECpkaIejFbXM9oqFGa3sULRsCWQNpJWL3M+Y/v2+N0GLUC9\n6aY8f6brYePGPNHqGQwvv5zH07U6upR8ZUGgwqQt+DOaqNpNGG9kopeH0va7TdDy3mEPtKryvJZX\nYGpGYFXlxbBm4JoxNjQUn5kJcvNuOPLx6KP5nGs+RmWA1kMBdaOWc0MlqGGp62qiGKhCUdonUBd2\n55yTC3lVEFqgqIl6VTAa5mI0KZRjjkkLtZRDaaKTcqzSY3BxMRt7PfY8bIF7C0XzKZs354JA3WmN\n1XrwFIpndTSxrDzMnZuzetRD6aeuwPteVTBnnJF7OIy2CoWT3awMHngghTyeeabe1sasXA1RXnVV\nzgTSglUV/uqBqrejr/esSBUMmkPzoO9p45Gq0B5P2FrXgNaBNaFtDqXtNfHc33NPMvzWrs2NCS4B\n2L07Khp739135+Ft9TZ0rN6lhryBeoHq0qV51ELJA0DdAPQMIoUaP1olP1FMWKHwZm/zGmV6aW5D\nk3QeK4S9IPVojjkmCRVVeEBaeCyoSjkUTqKpu8nQCndmSxk10KB0X+0bpbUhQF2wq3JQqqAHtU4H\nURAF1Kvl+6lmN/RTHe0VN9r6aKtQFPZ8WTGz9aoeg7JzuM/T6Gh8vVmn+/fH93Ll8+ho/nmaSL39\n9jyhv2pV7G3XhIceyoWJp4TUi/HWmod+E+8altY9q0SNI48sh7G1Fo3BCoULNpl8wP3vtLaLc07c\n4cCiFPbZl12WJ+PVe1m5Mp8fHWs+TGtXnnqqnsdTj0a9Wi+8pXkab341hKkU+oliwjkUFdiKM8/M\nLcrTT88FQluefEngeIl59nAYrNz42ks5FPZQuH+OLswnnkgCRRWIsn3U8lBr9vbbc+u1FKbg58ZV\nwIDP1upVIT1eaIPINh6J9m869th8HSmVmKEChtu7vOlNubfUL/iZsuLXEOYzzyQBoGtB49rabkfr\nibxQldK377wzXxP6nUBdKfE1lu6xjXU6nlCVepW92q+owmE0KRQg3Q+HqrmtChtb99yTIgdqIHJj\nyZ//HPj859PfuF3OK69EBch9vk47LV2HGXs2fvXV+FqObnDfMKDe323Hjvq8KGvTY64x69C7Rw/a\nTWOiGGjIy4MmYUtFhJy818S8hjmaBIeGyBicIG7joXDhHwsYZYMAea0Cu7fajoWTt0A9XKXWqLrX\nnrLQ+gQNwXkYlGWirK42UGqw9uzScVMBK1OH+7WkOQmvdSFsyXG9ySuv5O1QtPr9xhvrjf+UvcX5\nFE4EA/68qPLXkKeHQQmTJiZVCSoXeikUHbPHOm1anTZuf2NlxwqB68O40p3JMTfckGjb2uF5/foU\nNWBjEahTh73aEx5rct6KHRkqA266Kc/XeWFxpSx70LoZD4c15KWtVLQaGagn3bV2xINSQ9/97jyM\nddFFed5E/84CR0NsrHyYwsoKRXMoLDBsEpktwQpFaYjaaoF7SGnoROsIgHoxmzZ8bBP+ULqoh0Ep\nFA15edDq6V5dh9/2tpzppa9nKDW9H7DiXb06PVc9FmBoKClxTc7z/GgdgwlCm09tbQ/Ue0GpB+uF\nNpQ16L2mrUIZVMiLhbvXfqUp5KXjJi+XQ5ys7DhExAw5TjrzHmYhzkJfmXFajMjJ9NHRaExxeyUt\nNFTjU9leTGE3aDdhL9muHqi2bgLqeRfP21SjYiJlHUCfCqUp1mnQ4kbvPWpx9krMn3VWnlzXvIm+\nn3H22Uk4sYfCioc9FK3v8MCLV+tP2O3UfkJaCauCwUvW6fWodeoVq7Vhe7QVJr0WmIa8PPSiDrdp\nc89MLzZslKreCyx82evkTgWcP1GDgdk3e/fmrTO0O4KGqjRMYXkOFtjaO0qZPEA9fOkpjx078uSr\n14ZlUCEvrWZXo05rSzjPCfTXz4uNWFYoHPJSijeQU4TVG+WEOCuQ0dHc+9dO0dpr7d57czKOJcD5\nuSvBRtv1eHRhpSwbNOyteVQNgWkkw37Hikjf0y/6UiieR6KLVBWKh7POyi3Sc8+NWtmgSkjrWbR/\nkPb84hg9M73YKznqqPQZTb28gDrDi8ETqfFwtTY1n6KtOzQZq94SEIUAC4rHHssXOVAXON71t/FQ\n2limbRSKFjeqxarzO2NGnmNp6st29NHtGUNAnTJsgpifKyc3ubuwKgCtdtfxrbfmgkDDY0oH90Ib\nygjzBHyJ3dPUogNoHz/vtQ6UCq65Ve0yrUwxbbfi9fMylBqUcsjLDD4mJnDomRWBhoU4JH3HHXki\nXMNdmv+8/vq8c7B6N3quChBzI/z9XusfLb7VrtWAz/LrRebxXqfKtl/0pVD01EYPWi0P1G9MmV7K\n5PLQ1IJD619YiSxYkEJuJUqs99kcX2cLUE9N5CpSraxW1oVZtAZlcihTR5PBHtrQRTVvAwwu5OXN\n3YwZ+TpRD6TfmLyGQLR9Sz/gokbuQKuwa2Slo94ohzOs0tk2o9foTzn/alBoLB+oGxBeSx1dZx48\nYTIe2rAH70yUXrlVhvZdUwXDUA/F4J3WyK1nmAbMkQKmfGuomPNdBw/G77Xv2Lkz3rft/b1768l3\nVQTK7tLjCYB60bPHgPRyaOpJe+9T78iDGqP9om8PRRu16UV7IYiZM3OBoBanBz2cadGiXAlpUzkG\nM73aVOt7HgpvQI5p8mSqpcB0QbU2lZnjWQKaP9HW1h4jqE2+xBMm4+3b1eaEvRNPzBW8JmU9eKcx\nloRZk0JpOtUR8BtD8jwyq87CY5ycN4PBaNf2Nz1MTcMfXhhm//68UlsbiHq5EWWRAfWWK55HOt6i\nRqC3AaDHUij1v1+MJ+Tl9eBjY8t64VVVbgSywmaBv2dP/C57nVKFNfl+xRV1ZpiesaNeono0fJ38\nOVwACfheDNCbDtymuPmXzkPxFpMKAK8N+pw5OT3Qq5jnTsOaqGdl19QkkmET4HkoXIPCsXaOaasH\nwQtVG0lqHFTj615yTi1IbyGpMPHCZKWixvE09vPCngotbix9DysM9WzVoGA09WXjQjhPuXgKhT0N\nDjsyxdfao5jXqu1yrrsuFw6qHLT4Uc9OsWtjgeN5n228Ea8Ibjw9nID2HgqHvCaKppAXy6BSDsVg\nQp1DXxs2JC+Q+99VVZ5juOyy3JvQY3r1GWsOQ0959Aw/jUho0SHg509U6Xj5MS+vpjldD4fdQ2mq\nETFoLYoHrR9R+q9STbVCnsfqlbBVzYQApiaqh6JV8iZImMXDcVkNhWn+RFspKO1P8yde+40VK/LP\n8OAtJE+hjJfN0UahaG1Q6TWsdHoxvdgK1vmaPt0vbvQIIaZQ2GPk58qUUc6BqIHAVdZVlStyY3tx\nDyUNO153XX5UsPaVAur5E6B+sJrntWqXYcA3WBSeYGpDG9aQVxs0eZLaz4s7cKgM4sLGWbNSISmQ\njC02ApcsSSy9K68Efud34s9K72cChe1t+y4lSmhuRBtJAtGDse8C6v3dgDqDEKgbAt6+5S4OBu0y\nDDSfWmo47B6Ktpj34NE+WQgA9cO2tMeTd64C/109Gn7w3P/Jo5l6HkobIcytvpctyxkf2glUw1nK\nplBrVEMoQN0y8aiCnkLx4sttofc8Hg/FQy/qsBI3mo50Pumk9H2sULxrtVAAC93du3OhzGeWmJLg\n9igHD8bX2HyuW5czs5Tt9eKL9X5J3AsKqNerAHWPxesY7XkxnjXchr3I9GhGvyEvD7rH3vKWPM/S\nZJgy7byUQ7HP2Lq1Htrj3Bd7qNzQlcNPTz2VH/vsNYJk5aDhL60dA+r5M/VwgdyY8e4PqBdFAvU2\n+PY6voe2OOweSptqea1DAOoC4thj/bYZ/FC9Ajf7u0c9tDg+J+mVcQT4Hop5JRyPZu738HCu7bme\nQGPYukA911YtT+0HBtQ3tkcV1EKsEtp2KFC0USheCxQ1IFShaE6slwfD4OLGXgrFrFbvHHlrSQ7U\nKdw7dyZ2mIYrr7kmZ28p20tbtXiGgM635y1wrYzBY3h59NA28PJqgwp56f72ak9KbD1mEpZyKAbL\nVbEHal0E2LDinNLBg/GfGYc//SnwxS+m7+DK+6GhODf2nEZH62w59TS8ztB6XLjm64B6PQpQ7/MG\n+GSONt6I137psHsoXj+vJkaWQZleHjSOfsEFuRfT1HmWQ2DcEp0ZR1YkZdfrbR62ZDkhv3Jlns9g\nC1BzKxpPv/XWnMmhDA2gbp2ysDNozQJQP7irBL3XUkt8RRuF4hEh9DAs7ULc66Q/7ZbA7VnY6+yl\nUAx2TgSfiscsPe6/pNXNt92Ws/m4s7NXvKiJVS+0oR6rdpQG6nVMQPPpk/2ixPoahIfSb3Ejg9uv\n9FIoVtRoxhavc/YKONzFzCorVjTjYfny/JlffXW98SMbF0NDcf2zML/qqrzWxIxXjlhoSx6g3r8L\niF6MJum1RX1beLVoh91D8SiCbWLoGiMH6gwirYjXTsMXXZQn6rnAiosj2frhOhmrSzEPxXP5uRKa\nLQa2KPUkNe0uyie5AfXeW0o39hRbiSpYaq/fBO/z256F4hkVbZhfmgPrVfuj0FoGJnr046EYTIAz\n0aE011xJrXRg3eAaglLlAtSTsZ6xoDk1oN6vqy3azrf3uzYeSqngmd/rtV9hD7RJobCHouuv5KHY\nPPAZ7Fy8yhRaftYaelTv88EH8zymGgdaa1RV9VC11h8BUalxyBzwk+gayfBQmm81NL06pNfUQ/E2\nqKdQPOqwcs29+gU9XOm88/Iz55WffvbZcUMb2CvRtupeAt4Uinko1seLJ4RDI9xWnmOiWjHLQkSp\nn1VVXxTqEntJtlK9SS/r0auk9+ov2ioUz6gYT3GjB1VMTWdvcNdhLm7k1j9t6KvcysYsXq1W5ri1\negla6azJdu3t5YU2tBstUM+x2PUxPFqxJzi8My+8KvmSQum1xry/axjs+OPH389LPRRbf57gNOFt\n+RHzIHjfcURg374oi2yPsILYty/KLfMczXC0+92zJ28UCUTDjyMXelIn4CsP9Qw8xeFFHzxWp2d4\nlGrRdA1orVy/mHDISzuLAn61vIY4PChzq1erFw2taGLeY6x4125KxhQKt2Iwd5Jj2hrf5mpbPddg\nxYp8gWllrF2X5k80qaesIW8zeYwf79hXj6Pej0Lx1kAvhXLKKb17bmkIU09zBNJ9KxXdqzIvNZVk\nIWnKla1BVjLaPVqbO6qw1noIpXp7MXAlUnj5E++4V4+A4eXQPIKJZ52ON+TlQXOrnNME6sWqWpzM\nf+8n5KV49tkoezjfyOeVsEGwd2/cY6YgtBBx8WLgD/8wjbVy3vJufE1KJ1f6ORDvWxWAF/LUY4aB\n+lHUgL8uvKJGTw4c9pCXp1C8anmP6aW5Fq9+RCs+tZ01f4YuSrZwvWIpz0M5+ui8BsXAfXyYQqqC\nTHnjGm9XdseWLfWwlbJ0POWhYTSgPWW4ZJ22cXU9heI1iFRaqIasgHrXaI/pxR4rJ3Y9okevazVB\n7R1IxmEPpnyz92B0YBO62jNNw1/mHbKxoO3svd5Nehws4OdPvIR8aQ2oMPHi5/2EvHqFwjRy0eto\ngXnz8qQ9RyP6ScoDvkLmeeS9wxTgn/88Z2/xPHjJdz3LSFutjI5Gg5j3ldfcldv7GLz8SWkNaOSi\nrRzwjIqmNlNt0JdC0bAV4LeX9gTMmWfWrU2tTyklczkxr327zj03D3sxzjkn/Y2ZXibsbNGbh2LW\npbm2/HA5h8EHainfWw/X0aN4lbN+++29G7957ei9hPxEFlKpcr4NbdgLebU57EqNDD0/XpmATUyv\nNow08x4tjMm1Q8zS4znjokfNfWix4jXX5OEv7d0G1IX2qlX12iIvf+I1ifSKHL3GoF5RYz9GhT5b\nTw4oVKFxnrKuAAAaWUlEQVT0KqTV/l5sADIBQ2WE95kWMubcpiXamb1n547YZ7C1z+fQA1HA8z7V\nvwN1I0+ZXEC94h6oN50E/N5rXihLTw8F8oJNfl2bkBcwPo/U0Hf7+vG+zjtEx2N6qUeiYawLL4xC\n26AV9OxFcQiNFYp5MqUcinkozPDijcobVBsAcr3Jnj15fNOznLQ9PbOODF5C3qMUehXSXsuPtiGv\ntrRhL+TVhpihCkXzXlrgpgqHMXVqXcjptRojyhSKhqOmTMnpplZvYh6GtsrQPlsa/tLD0rz58eLp\nHtnCY/JolwSgngAG/DM42hoVnieiNWAePPZnP8W0XC3fq9OC/W8GoFntlj/hOWUjgKvZ1dBbvDhR\nh4HI1OTIwqWX5n/3an+83Njevfn8HDxYz6l5uTG+T4U+Hz1CHKjLIsAnjEwUEzpgqx94i0Lb0gN1\nQXTeeXkFvRZSLliQn6XCreyZaspEAaMOG93RKnFNoZjlqmdd2CTzvViMFqifUaCWpldHoJ/Hp8oZ\nmKli8KxJr/U0HxJmaJuQ9dA25FUqbtSC015NHvn1moebMiURLLi4sXStplDMqjOPg5O0S5emHMk9\n9+RsHK5jUCvwiSdyZTE6Whf4mrAH6jUjntAYGZmY1egpnn5CXvrdWnTsoU19moK/p6lBpF6b5RfN\nCn/kkRSdOPvsnLbPEQUmv7CCGBmJAtj2zc6dcd7NUBwdrRfDaquVkZH43NkI8A68Wr68nlPzvFEu\nwuT799DPWpnIuvIwEIXS9qL0AXhJfj3bpNfi1e/mkBi70fw6YxzZ7+x1lkMx5WGWLLdzYEtGPQ6t\nN9HDttSb8dxVbdFi36NMrX6gz6gkTMabQ/GOafU8FM1jqQfSCxru4Db28+cng8E8TrtWs/hModj8\nGk2cFf/q1Smsxew9DS9dcUVej3DVVflYq+WBvEAO8Bk62uLFvptj9YBPwOgHbWnDgN+lohd7zguD\n6uc0sfiaGkQaTKGYorcaFGZDhZAEND8zNv6qKu4J8zZvvjn3LC69NK+E1zyY9gED6qUAQCQAcIGr\nfZaGvL2Qp9cU1MujTvSArIliIAqlbeKuTd8v9Ui8z9IqamaDKT3Ruw5tp25Qho6FGTi2yvkTTZ5q\nSENdTy6YAupxWftO9jLMJWZ4LTj6WUieJ9O2+3Apj6Y1JZ6HUnruDD1PXg9jYmgtirHILARnCsWq\nhrkIkJ+X0bYt3GoGBj9nzY8orVcLDPUYV095eI3/tKU90D4Zq4y0ftE2Kd8m5OWtCQ1nK4kGSAbD\n7NnlUxsZplAsVG3hXcuVMOmB55AZWhoV0DnQIwO0xbxXV+IpCi8H4hlyXkiaDR2DtwY8NuBriUMW\n8vLYX6oIgDxsAfgHJrElCtQLHN/5zryCnhcth0KMRlo6lMnj8AN5DoOrnlkgKG9cmTXc98ugrCCl\nBgN+WwXPitUT+vpF25BXW2/UqwFpU4uycGFO3lCmF7PHehU3mkKxkM+OHdGKnDo1T1KaEmHDgbvS\nArkxoHkpLWwF6vFptXoBPzemXgyQqvoZHsOrzUmdTejHQ1GF0uvYAKBem6K1J23o5wzb5xZZeO65\nOL+zZ6cCRQ4zW68sO6rXhDv3aduxI3pXtpe1T9auXdEw4IiBtto5cCA+DzYmvWiE1yGDOyMz1NAE\n6sdHA35YjbugH2oMRKFwN1CDV4viUYf1PHkPyuzSTsSamOfiNk7MNzUYNJhiYIaXMqVsg7GVqnFP\njsUDPuODPwuonz8N+C0Z2jK8Ssk9DyVhMkgX2vNQNOyhiXdtEsmUdPZKWFiZQrEQmbWV2b8/buIF\nC1JCnrvwco0J5zs0YavNAbmuAfDZVytW1PNnGspsiomrUPdarnhdhvvxWjyjp62H0tSLy+AVN7at\nlvegIa+9e+NcLVqUajYsmc7t7VkZW/Gi7fXFi+uNILkWRanBQ0PxGljY694HYriLGYFAXREBfiFk\nSSF4nT28LsNe/uVQYSAKxYuXtq1F8ZheWpB2/vl5xbxW2Wu8lfMwXG3PTK8SrEZBmweGkCezVVhr\nXFUtBT1PwWNfcUt0g9c11GMLeQpl27b2rVk813vQCTtlcQHxPrSYUc+bZ6+GFY6GVUz4KWWZG+XZ\nvFrhGMe6mfnEeSZVGOpF6FibQVrYhdeL53V4uZJ+FIJ3xkqbU/oY3py38VDaNIjU4kbNs3khMIYq\nN1YoFqq2AkYLNxsphQ0EZnpxT6+qymvQhodTLs6gAvu66/JWK0AeFjd4tW3envWOLPAiEiUwScgw\nyOMremEgCsVjdHjtV7xqeW2fAsSN/vDDaey1aWnKxzDVmMNbbU6KtPCFTeLevWlBcaxU6wd27kyU\nUa9YTa0JL3/itciwe2V4rciVYQa0O6HNMN4THCcKVSAaKlPvt42XqQqFWU6mUCyPZN4e5zg0FMGK\nWfsreRRfJT3oWRuATyn18ifeoWolz9OzWMd7sJahrYfS5kwU7eelB2k1eSilNjoh5N7Hs8/GPTR7\ndp6LsHzD6Gg0nkxJc23RqlW5d6BexYYNdQ9Q53ZoKK6rXmfWWGhL97YnA/TsJKC8Xz3KsGdUHArK\nMDAgheL18/IWh5eo8/IZejYKUE/oafdZPvFRa17sfcwYM8qwtlU3S8J45VzZyj/z5tdwkca3vQWl\nR4R69Ewv+V6Ct5A8L6if978W8LzWJjSdtWLPTxP53ElZ63xMObOHyQ3+VLmot6JsLy9WvnRpnsQF\nfIvVE/5ensULa5QwUevUW5eeQmnbcZgViJ7c6h3rbdfqRUHsb3rez333xSS7ESOs4WEI0QOw5o+6\nP7RQVZW5tlqxiAU/nyVL6oWL2jDSrlGjEaWO3xpyBfyW+IA/t22LmweBgSkUTcD3202W4R2CpSE0\nPQJYOxMzvJ5O5i1pIpCr5GfMyDc0h07YkuADfIB6sZrnjQD5QvSqZbUyG6ifDNkEo1Ey2nQrHRQ4\nl1WCZ4wouL2OCrejjkqhFqYLszWrZ0NYGPGFFxKRgYXHk0+m0JEKGVUEzz6be4b6eqDeH8kTHCUh\n79UReQcqleC1XWlLETeMN+TlHarVRAXW9izM9PKo6RzyOuKINMdG11+7NkYZuDqdqf1MBzb5YAbc\npk3xudm9j4zEeWQWqOZT+LsZVhfD0JZMgG88AH4OzSNllNaQxxorVclPFIfMQ+kHyhDx4rhaMa9F\nkZq458W4aFHK09h3GeNI2WhKHbZqaHYRtd+N5j5U+2v/Hs/iUKUERHdaF5jXNK60kLyqe64aHiQ8\npe2x+lTJeHOt60mLVxncddgrbgSiUGMBbsrDWqNUVUp8avtu5vrrMboes06VuBfT9mjAGzfWDYrS\nvHqFrqV59TrUenVIJbQNeXkeihpr6pEodC1oP6+SQtmzJyqiWbOi4clhpxAS2WV4OJ1VokwvzqUA\n8ZRGbgSpNWZA/fRE69TL92GhKS/6oCErL/eya5d/2qrHLCzlTD0v01My/ZB4SjhkOZQS5s6tu66n\nnlr3SFTJaCt7LXLT0Bmzu/hnK4AzxpF37V6DNBYCeugOt8Xw2FLK0vA8Fq//jrdAPIXSdC6KV9TY\nj6vbNjzihSS84sY21OFexwFzm/sSdZhx4ECco5kzo3A1KrbF1R99NFmQLDhUIWu4S70RL7TpVcd7\nhWte/sTWqcKrfPcSvCX0a52O10OZaLU8k208uWGvHx2NhsTQUFS006enwlI+u509FQ05MUtzdLTe\nVkebNXr0XD07BfAbP3rn3wD1vn+Az/Lke2f0swZKrXcmUscEjFOhqJCZM6d3/NTQlumldGJvETed\nDMkeDVu4xvSyPIsXrrPELQsIbjvOloTSN9XK8PIY3mIEep9aCNSLI+3zNFlYgmeZlNAPy8tTKG2L\nG0OoJ961SSTXovB8nn66X9zIOHAgVcwvXJjyRabo+cwSrnDXYkZl12ndhxcr98KO3sb1cipeMraE\niQqTEjyDwute4SXl20YuSkYLJ+3VQ9Gc6qZN8b7274/zZ8WHXGDKe5jbyqsxp1XuO3dGGcfWu7Za\nAfwQtUeq8OjC3nG8pfeXntdEjYqJHq4FjEOheGeUlBhXXrt7T6Gcc05doXiJeSD/Hg5lAXmIhROA\nLBhNoVg1tIW8OLdgjRd5MlkIcKM1rZrVehNtV28bga/J44kzc6UXvMZ0JRyq2OlEPBTNj512Wm5M\naOiEmV5szHD7FcaBA3Gz7NuXFACHq6yeyHpm2TrgehJ9bnruDZAK5wzeBvWef0lAeIVrpZYrWs3d\nhH5DXm1ar3ghrzYKRY1RfhYc8mIPxYxLvrZNm6JnsmpVTMjbOTemlHftikovhGhcTJ2aBPjixXlr\nFWXg6d+rqs7G4s83eMcXAH73cI16GLzD8DwjBfDDYE3h8EEf/wuMQ6G0OVPc4BU3eg0BvfoEz2tR\nwaMFjVpBr6iqOqXUFIqFrY48MnkoFnvl6lWl6+lxvtoVVr0RjbsCfusGr8ixVEXrsTi8zsbAoWN3\neFXOXrcEz2tRI8M7xoBR6jpcoq4ePBg35q5d8ZksXJjyJzyf996bLFhtleK1qudaEy9WzueqGNTA\nAMqsLY95Vzq502unU7J6D0XIyzvvxvP+FboeOOqgIS89ZIsVyquvRoWya1e8tmOPTT34pkyJZ518\n/vPxtdx3bXQ0vseex8sv1yvh1TP1GrhqOBTwQ5slJpeXQ1UvzFBaA17PPy5nYHjdqw+Lh+I1dCzB\nUyhtjmUF/NoTTczrOeVaQT9zZkoAe8wxIDFLhofjZM+fnzantQLn+Cp3oNXGk8rAsgXPG9Lr3+Sx\ndvj8a0OptYZXl+KdpwD0F/LqB6VDtnRDeMqiX+pwmxYdTAc/cCB1lV63Lhoe1huJwxvsbd5wQ577\n0LCitsRp21rF63Lg5U+Umcbv17BKCVu21CnMwKFJynvwcih6zLN6sccdlww+JtZwQp+NWttbRryY\nMiXlLVh5M5GBTz7U47u1EaTnJXpnmnhFqXp8ARCNFlVGgO+JMNuQ4V0T4BuQ/dQheR5KyTAtYVwe\nStsOsV61fBPaNI7kxLxnObESuuCCFDbT44WB/HCw/fvjJvGSoFzlziEtFTJqkXoVrt7BPJ4lqt1L\ngfwY014o1aAcKpZXm3PlS/AMDyVlNNU5aD84IHqiJswOHIgex+zZqf7H8ifslTAjipvxaXjLa/So\njLzR0bqi5+aTDK9rbMkK9VqulPZNSZh4DKMmtPFQPOhRE0C9q0Wvs+U9sIdi2LEjzsl55yUKvu1b\nToJbnZDdk1L8VYhffnmeKxkerrda4UO7DKOjcY2pUeAZD6V6M6/DMOArn9Ia8GjjJXgeSr99wA65\nh+LFtEuv5QaQQJ3Z4XktSldlq5g9mkWLUvW9eSUsyIaH4wM94YS4CPhQIu4SykJHD91Rq1UXq9dK\nv5+KVa9HVMmCKMVZ+22F7y1UjwzR5kjeEjyvRddDU4W8Z7hwaPPgwaS0qyquKRMA9nt+tswMAurJ\n+euvz70X6wjNisLryaQ5FqAc1vCoxSWUvM5+OiWU4F2bnq5YgrcutVp+ogqFowMvvxyNPVufIcSf\nuSCR8yF6zomGrw8ciPudrXadeyAPpxlKORGPxs1GDcPb7yVoBwfDRD2UfgueD6mHUlI+Xs2ClzM5\n//x6K3vdgOyFAHn4RI8Rte88++xoxXCCeHg4WlN790bhZeEKXrSaKN+yJe9YCzQn2z0X1mvLUWpB\n7cXJSwtmEMKk5Op6ebS2QgZoZucZep0vz2uI83L22axQjDZsTB3rl8TeItN7NSmrikDZPF7/JW40\nafDCYk8+6YcxvXBVqU6gxO4ZlEIZZE837eelRqNXLa8wuWLXZqetjoxEb+9978uP1t66Ne7Tqsob\nump46/LL80aQXtW7p+j5oC6D11qn5C3oQW4Grwi55M2U1oB3WmcJnodSytuWMDAPpc2hO4Yzzsgb\nAgK+QvHCVF5inhPxOvaEnDG9mMJqRU8bN0ZvxoQGN4m0I0WBerGceide8s3Ln3jU0FLFrIcSw2v7\n9om1s29CP8QMrxO1V4CohoIqFG21c+aZacx/s5g8dxweGor5sJkz49xbDJ3pm1x3wEJDjYh9+6Iy\nYyGrNQqAfxSvZyh4IZASSsn7Em1cW5I0oelUyEEqFG2/op/dRqGwUcu9vEJIFej2XLnFkdZ0MEFm\n//64p1mgaqLc6pL4mj2BXVV+PZpHFwb8zgVe52eg3CiylF/1apZK+KXyUNowOgxvfWt+5oW9XxN4\nxx5bP1zp3e/OK+L1qFCtoOfEnlVpL1gQFRon6oeHge3bl/1CgZhraucqAHFBGo1Yhb7mTzyX96mn\n6rFyT/B4R4CWaMRN3PNegmDZsmXNLyigH4WiDQABvxZF14+yf5RFxAqHCxrtZz4TZd+++PxWr16G\n97wnWn7TpiUPkwW9Uj21m6wWqlVVfdN6zfhKlqXXOdrrcACUW654RbGGtnNcyqsMuiuthrwUs2bV\nDVO9Bg55GSU8hChQN2xYhhkz4t+mTcu9Di48VHLE1VfnYU3Li/Ae0l5eQN2rAcpC3/NkSk0etbuG\nQdeAze+mTX6I20NpTn+pPJR+qmL7ZfUw3v72vBOxQuPxfPiWvdeK2rhAa3gYeP75ZW7fJ9usLDjU\nstQT/DQM4oXESiEMz7opNYMrxU7bwBZjv0KjH4XSthZF10QvZcido3nOrbjRGkSaQnnlFWD9+mWY\nOzcKa2vPAcSCRAt3qbehcXU1JDx6t1cdz4aJwatJAsrxd6/lCtBsRbZVKE0HrA3SQ+Hea23gndzI\nCmXfvpi7Gx2NRsZLLy37hWdhx3rPmZPOLTGlrnkPJVVcckkeDgN8heCFnL3597pvAD4FGSgXtepx\nFryH27ZN8SrygV9CD8VTKF683OvLA/gJSqUZeol5bRzH7C1OzHshNHu9FTsBdatRz0IBcs/C4/sr\nE8MLV+gRwkBZuDcxvCa64b3T4JrQj0Lx6k7mz++tUErgI2K9dWjFjUYznTkzteQJIQr2D384/Q/k\nc8OGwJ49qUstkOZZu8tqGMNjbd19d11wPP20X4y4YkV+NIKhn/BFCaX19VocsAb0v1Y5AjFlSpx/\nVSjbtsW9tn59NDSM/s2hSyZWDA9HQWn3a7kmvjbtXOCxM72efIBf++HVmgHlg/e09YuhH8VRQqkO\nyTux9rB5KCWF0uvQHH2tusOWQGeo4tEw2KJFUVgDuYvNoRFWVPPmpdj49OnJCmW6IVuuXCkP1BeF\ndyiSF2fXQ7cA/1wTwG+LXkJJCJSKqpqsU++z+mH6lRpE6vv1oC2gbky0oaFrx9qZM+PYahRWroyh\nBNvgfO9G9bVntGRJas8BxBb32oZerbrSMdJeIVkpf9JPN+DSXDd5v95nD8JD8YzHNiSeJpIGE2ts\nbnn9WDjzpJMiweHNb46fdcQRsV7E5o89yWuvzcNbekrjQw/VQ86XXVYPd11xRZ3dVdq/nH9leKyv\nEkqMwCYjwfvsUnGzp6z6pQ2HqsEECSEM2D7p0KFDhw6TAVVV1UyNRoXSoUOHDh06tMVrdNRShw4d\nOnSY7OgUSocOHTp0GAgmnUIJIfwwhLAthLCWfjcvhHBzCGFjCOGmEMJc+tvFIYTHQggbQgifpN9f\nGEJYO/a379DvZ4QQLh37/X0hhEIFwGuDEML8EMJtIYR1IYSHQwj/dez3k/KeQwhHhhDuDyGsDiE8\nEkL4H2O/n5T3S9c0NYSwKoRwzdh4st/v0yGENWP3vHzsd5P2nkMIc0MIl4UQ1o+t6/e+Lu+3qqpJ\n9Q/AhwBcAGAt/e5/AvjK2M9fBfDtsZ/PBbAawHQACwA8jpRXWg7gorGfrwfw6bGf/zOA7439/EUA\niw/z/Z4A4J1jPx8F4FEA50zye5419v80APcB+OBkvt+x6/jvAP4fgKsn+5oeu46nAMyT303aewbw\nYwD/kdb1nNfj/R7WRXMIJ2cBcoWyAcDxYz+fAGDD2M8XA/gqvW4JgPcBOBHAevr9lwD8E73mvTTx\n2w/3/cq9Xwng42+EewYwC8ADABZN5vsFcAqApQA+CuCaN8KaRlQox8jvJuU9IyqPJ53fv+7ud9KF\nvAo4vqoqq4bZBsBq2k8CwP2QnwNwsvP7zWO/x9j/mwCgqqqDAHaFEA7BGYj9I4SwANE7ux+T+J5D\nCFNCCKsR7+u2qqrWYRLfL4D/BeDLALgn82S+XwCoACwNITwYQvhPY7+brPd8OoDtIYT/G0JYGUL4\nPyGE2Xgd3u8bRaH8AlVU0ZOOKx1COArA5QD+W1VVWXOLyXbPVVWNVlX1TkTL/cMhhI/K3yfN/YYQ\nPgvghaqqVgFwSwwn0/0SPlBV1QUAPgPgv4QQssY1k+yepwF4F2JI6l0A9gL4Gr/g9XK/bxSFsi2E\ncAIAhBBOBGD1+JsBcEu1UxA1/Oaxn/X39p5Txz5rGoA5VVW1bIt5aBBCmI6oTH5SVdWVY7+e1PcM\nAFVV7QJwHYALMXnv9/0APhdCeArAJQA+FkL4CSbv/QIAqqraOvb/dgBXALgIk/eenwPwXFVV1if9\nMkQF8/zr7X7fKArlagB/PPbzHyPmGez3XwohHBFCOB3AQgDLq6p6HsDuMaZFAPBHAK5yPusLAG55\nLW6ghLHr+wGAR6qq+t/0p0l5zyGEtxjbJYQwE8AnAKzCJL3fqqq+XlXV/KqqTkeMid9aVdUfYZLe\nLwCEEGaFEI4e+3k2gE8CWItJes9j17kphGCHEHwcwDoA1+D1dr+HKxF1qP4hWnFbAOxHjBn+BwDz\nEJOaGwHcBGAuvf7riCyJDQA+Rb+/EHERPw7gu/T7GQB+CuAxRIbRgsN8vx9EjK2vRhSsqwB8erLe\nM4DzAKwcu981AL489vtJeb9y7x9BYnlN2vtFzCmsHvv3MICL3wD3/A5EgslDAH6OmKh/3d1v13ql\nQ4cOHToMBG+UkFeHDh06dDjE6BRKhw4dOnQYCDqF0qFDhw4dBoJOoXTo0KFDh4GgUygdOnTo0GEg\n6BRKhw4dOnQYCDqF0qFDhw4dBoJOoXTo0KFDh4Hg/wN9srpwmG+W1wAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant locations with y axis inverted\n", "anhima.loc.plot_variant_locator(pos, step=10, flip=True);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAABUCAYAAAClWsGgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfWmwXtV15Tpv1AQa0QASSGJGYpCZR03MYDEb+JFKdVd1\nfriruqu6KnZjBygnTuxynOq0U5V04thdjqsLz8Z27NjY2IlJYoMAMQihASSBBpDEIIEkDBpu/zhv\n566z7r73u/e9Jwked1Wp9L77vvd9995z7tl7r732PiHLMrRo0aJFixZDRdeRPoEWLVq0aDEy0BqU\nFi1atGgxLGgNSosWLVq0GBa0BqVFixYtWgwLWoPSokWLFi2GBa1BadGiRYsWw4IRZ1BCCF8NIWwL\nITxLxyaFEH4eQlgbQngohDCBfndPCGFdCGF1COFqOn5uCOHZgd/9bzreH0L45sDx34YQTjh8V1dE\nCGFWCOFXIYTnQggrQwj/beD4iLzmEMKoEMKjIYSnQgirQgifGzg+Iq+Xzqk7hLAihPCjgdcj/Xo3\nhhCeGbjmxwaOjdhrDiFMCCF8J4Tw/MC8vvADeb1Zlo2ofwAuB7AAwLN07AsAPjHw8ycBfH7g5zMA\nPAWgF8BsAC8ACAO/ewzABQM//wTAtQM/fxzAXw/8fCeAbxzh650O4JyBn8cBWAPg9BF+zWMG/u8B\n8FsAl43k6x04j/8B4P8B+OFIn9MD57EBwCQ5NmKvGcDXAPxnmtfjP4jXe0QnzSEcnNlIDcpqANMG\nfp4OYPXAz/cA+CS976cALgIwA8DzdPwuAP+H3nMhDfyOI329cu0PArjyw3DNAMYAWA5g3ki+XgAz\nAfwCwGIAP/owzGlEgzJZjo3Ia0Y0Huud4x+46x1xlFcJpmVZtm3g520Apg38fCyAzfS+zQCOc45v\nGTiOgf83AUCWZfsB7AohTDpE590IIYTZiNHZoxjB1xxC6AohPIV4Xb/Ksuw5jODrBfC/APwhgIN0\nbCRfLwBkAH4RQng8hPBfBo6N1GueA2BHCOH/hhCeDCF8OYQwFh/A6/2wGJT/QBZN9IjrNxNCGAfg\nuwD+e5Zlb/PvRto1Z1l2MMuycxA99ytCCIvl9yPmekMINwLYnmXZCgDBe89Iul7CpVmWLQBwHYD/\nGkK4nH85wq65B8BHECmpjwDYA+B/8hs+KNf7YTEo20II0wEghDADwPaB41sAzKL3zUS08FsGftbj\n9jfHD3xWD4DxWZa9cehOvTNCCL2IxuTrWZY9OHB4RF8zAGRZtgvAjwGci5F7vZcAWBZC2ADgAQBL\nQghfx8i9XgBAlmWvDPy/A8D3AVyAkXvNmwFszrJs+cDr7yAamFc/aNf7YTEoPwTw+wM//z5insGO\n3xVC6AshzAFwMoDHsix7FcBbA0qLAOD3APzA+azbATx8OC6gDAPn9xUAq7Is+0v61Yi85hDCFFO7\nhBBGA7gKwAqM0OvNsuxTWZbNyrJsDiIn/sssy34PI/R6ASCEMCaEcNTAz2MBXA3gWYzQax44z00h\nhFMGDl0J4DkAP8IH7XqPVCLqUP1D9OK2AngPkTP8TwAmISY11wJ4CMAEev+nEFUSqwFcQ8fPRZzE\nLwD4Eh3vB/AtAOsQFUazj/D1XobIrT+FuLCuAHDtSL1mAGcCeHLgep8B8IcDx0fk9cq1L0Su8hqx\n14uYU3hq4N9KAPd8CK75bESBydMAvoeYqP/AXa9JzVq0aNGiRYsh4cNCebVo0aJFi0OM1qC0aNGi\nRYthQWtQWrRo0aLFsKA1KC1atGjRYljQGpQWLVq0aDEs6Kn6ZQihlYC1aNGiRYsCsiwrdG6ojFDu\nv7+oM167NsPXv148/tBDGR55pHj8/vszHDiQHjt4MMO99xbfe999/jH9+z/6o/T1Zz+bYe/e9PcH\nD8afN2zI8OUvx5/XrMnwD/8Qf7bvv+++DJ/5TIbFizNcdtn9uPfeDNu2ZfjSl9Lv2r07w5/+af4d\nn/pU/vOBAxk+/en0nPT11q0Z/uqv0mMvvpjhK18pXrN3b772tQzr1tV7b9m9tPHIf74fWZbhvfcy\n/MmfFN+7cWOGr361ePy7383w1FPVn23/Hnkkw89+lh7buTPDF7+YHluxIn5u1bXde28+rnqN9t7t\n2/P7fN99Gf7+7+N3nX9+hssvvx8PP5zh4YczPPpohgcfjO+75558jv3xH8extvvCc+3AgfhePqe/\n+Iv4nXzsnnvS88yyDJ/7XIZdu9Jj27Zl+Mu/LN6zL34xwxtv1B/rqjlgY9xpnDZvzvC3f1s8/pWv\nxHlQ5zPqHtO5qa/5b/7u7zJs2pQfP3gww6pVGR54IK5Dv/51hkWLMtxxR4ZvfSvDH/zB/fizP4v3\nes2aOP5ZFsf8pz/Nx/GTn0zn42c/m57DJz6Rvt6+PcOf/3l67PXXM3z+88Xr0zmSZRkefzzDt79d\nPP6FL2R4883i8U9/ujiHvLG28S173r3jL72U3xf+95nPZNi/Pz22dWuGv/kb97PdNkCNKa933gHG\njCkef/11YPLk4vEsA7rkW159FZg+vd736d+//TYwblz6nt/9Dhg9Oj0WBi73iSeAc8+NPz/9NHD2\n2cXvOHgQGDUKOHAgvp46FdixI/+cgweBsWOBvXvzv5k6Fdg20LatqyueJ2PKlPwzAGDGjHjdjDlz\ngPXr/etWnHkm8MwzxePd3cD+/fU+owy9vc0+4513ive7DK++Gq+d8fLLwAknpMdefBE48cTOn2fj\n+sYbwMSJxd9v2gTMnBnPcdSofL7aGJ13HvCb38Q58fjj8W/OOw948sn487XXAj/9afy5tzeOvc2L\nri6gpwd47738+5YtA374w/Qczj47zjXGkiXAr36VHuN5xrjwQuCxx4rHR49O5+BwYvToeK8UwV02\nBo+DB5t95o4dwDHHpOdj8y+EOLbTpsVjM2bEeXTDDcCPfwyccgqwbl38u0WLgH/+5/hzVxcwaRLw\n2mvx9fjxwO7d+TgDwFlnpWN4zDHxXPg5nzQJePPN4jkfdxyweXN6bMECYMWK4nuvuAL49a+Lx084\nIc5lRU8PsG9f8XgTTJ4cnx/FmDHFOdDTk96XTmhsUPbu9ReTN97wDYoHb/HQBbns2IoVwEc+Uv6e\nd98F+vry1889B8ybF39evRo47bRogPr7U+PU3x9v3NixcXIZ5s0DVq2KP48ZA+zZE3++8krgYWpe\ncOKJ8boMS5cCv/xlem5mnPi1h/HjgZ0702Onnw48/3zxvccfHxfow4kmBuWVV4rOw0svxfNmrF8P\nzJ2bvz5woOiIMNatiwsGEBf43t748+bNwKxZ0cGZMiWe69ixwNFHx3PeuDEe6+6OY5FlwNVXAz/7\nWfz7887LDQ0AXHYZ8G//lr9Ww3DSSem4A8A11+SfZ2ADxtA5AcT5/cQTxffOnw+sXFk83t8f5/RQ\nUGZQAP85HCx27vQdgbLveu+9eH0MNiijR8dnuKcnGoB3343O17MD2+sdeyywZUucS0cdlT9Xd98N\nPPBA/pnXXAM89FD++pZbgO9/P/3ec8/NHQ+DNyY33QT84AfpMXNo9PrOPdcf6wsvBH772+LxU08F\n1q4tHm9iaHgd63S8qcM6bBFKmcfo4cUX44PI2Lo1WnaGLQ6MJ55IDcpLL6Xe7jPPRO/CsG9fbmDs\n5xdeAE4+OS7EtrD19QFTpy76j4X7qKOAt96KC8ojj8T38OJyxhnRWBmuuiqdkPPn55PacOaZxcl3\nzDHA9u3pMW+SjRoVHxbFySfnXhijzkRYtGhR9RtK0MSgeJGrF6Hs3h3vedl79uxJ593atblB2bgR\nmD07/rx5c4xQXnstNyhjxsTPnjdvEZYvzyOX004D1qyJC5I9SCHEuWD3eulS4Be/yL/38svz+WAY\nNy46J4YJE4Bdu9L3mAFT8OJn8DxFIEY+XpTqGTXDwoWL/F8IRo3yjdJwRyjbt8fIrAxvv53OAw9s\nULq74/v374/R8Pz5i/DGG/ni+LGPAd/6Vvw7/nnWrBgB2AJ/xRXAv/xL/h1jxsT1gqPRj34U+NGP\n0nNZtqxoPGbOjEZMcc45xci1LAKYNy9dXwwaOdkzPGcOsGFD8f0eysZUWZiq8yvDsEUo+/fnXqKh\nzLPxPNRVq+IizXjqqSJFpZHQ8uXABRf4r8u+3yIVXrT6+4Fjj12EuXPjQNqiftxx0dgBqddgg2Lf\nMWtWGuaGUPRIFi7Mw27DpZcC//7v6bEyr8UDh/WM6dNzSq4Mh8OgeJTn9u0pjeGBDQaQRiRAfHjM\niKxfn0e89tlmUA4ejGN71FHAtGmLsHEjcNFFcRw5Mpk9O38gOfrs64tz24xBd3f8n40102QGpkQN\nc+cWF/6FC9OFjKHzd9YsPxotcyomTQLOPnuR/+GCKsMxnBHKtm2pQdFIdNu2znQ4z78si2O7Z0/8\nrI9/fBEeegi47jrgn/4pfQ7U8DLV2dUV1xWmIG+4AfjJT/LXo0dHI8Njf/TR0QjqPZo4MTpTDHU6\nDZMmFd9b5oDoWNszfPLJ0VFWmPNUB++rCMXDrl2RvlF4xmfVqkjrMFaujJ5+FZ57LjVEr7ySc/Yc\nvXAEtXZtHAA2bLbodHdHj9db1Pv7U4/FPFyD3vwTT0xzJNOmFaORs8+OhpNx1FGpx2vwBtyLcIDo\nJSmP2xRlC4mXs9q/P19s63xuFZ0FFA3KmjXp6wMHovcEpHSZfbYZFCDOte7u/N5dfHHMo/B43Hgj\n8I//GH++5JKU5rr44pR+uOKKNErx+HHj8RkeFebNCSAufrpAlC36J5/s0yDDMQfMMRoubN8er9lg\n1KShTn7V1iA7NzMGY8bEBXTlyshimLE4/fSctj7llPxeKTWlNJjNE8bVVxeNgkeFLVuWzyfDxIlF\nKhsodyp6e9P1BiiPGE46yXcq2CHuhPdVhOJhy5Y4qevAy8Go8XrrrWI4zAuL4oknohcCRFrBqDBb\nEHny9vfHG7plS/QMxo+P3wdEo/PSS/Fnzm9ceWVKh5x/foyQDJpnAYqceVlC3Fs85s8vhsFli0zZ\nYtLU4/A+/+DBokHYuTNSPYMB50AMO3akC80LLxRpUsOrr6aLFJAalJ6e3Bno7Y2v7cGxRPexx0Zn\nBIj3KIT8PqlnuXhxmkcxmowffs9jLKNCvDzKRRcBjz7qX68u8Ecdleb9+PvKDEpdI+GNf13nwUvA\nK+WlBqSOQbE1iI3d1Knx+f3Xf41z0yLT3/0uzYfcfjvw7W/Hn/v78yQ/kI+PfWYI0WHjSNOjPD3n\n4dRTU2fTMGFCMZF/zjlFpxIopzg9lCXaPYFAGd73EcrmzcW8yFCwYkX0Bg0aLqsCbOXKPCH/9NNp\nbgXIKamurji5+vp8ZQVPooULc1XGjBn5IgSkShLAT5h7k8QLS6dNK6rCypReHpSCM0ya5E+8oeLN\nN+vn0BSaBzPwYvTuu/E+Af5iqAsXR6S9vXFe7N0b58DKlflis3hxPmb8sHMUYvkr+15zAtgIqJEB\n4ndo3stLnnp8eVnUcfzx/hz1ULaYcM6oDjzhiybKPXhOxmuvpY5jlUEpM3qcQzFMmhS/b/36XFhh\niXampVTRdcstwIMP5p9z/vmpeOLuu4FvfCN/3dUV7585m0Ce81KnwFugr74a+PnP02Nl9FZZYt6j\nyMpQ5sR48CKU7u73WYSiBmXnziIN5qkfPHnwk0+mCfk1ayLtZGCJMJBSa6+9Fr0N/Z5du+Kk7++P\nD/vmzbnSy6S/7G2oAocXCM9TVO/Ty6NccEEa2QDFaAcoV3R5i9SECb7hKPNkhoo334wTneEptbxo\npK5k2OBJkYF0bDly7e2Ni0pXV4xyli/PIwB+aK+/PufMFy1KDYSOkdIhl14avWOGOhj2d7pIeJSH\nJ0UHiknZKkyfXnRKgGYLkheheAbFO1ePdTh4MI1u1KCwI/DWW/laYcpMoJhDsfO052zp0sgMXHBB\nfq95vFipp9HBRz+aysCPOy6NWgDgttuA7343vS5vXL3cmkdxA8UyA6A8D9okcqmKUHTMPAPYVJQx\nqAilrkHZujVSCQxv8di+vRjmPvts9MgZOkEfeywuvIbHH88prjLvRmm4N96ID1hvb7y2/ftzia4t\n6hxaK6eoNQOmDjMsWJBOIEsYMy66qMjVerUMZYPrKTzK3ttkMWkCT+W3Y0dR0ePRoDonVPqtY8n5\nFf6d8vGGnp7oJIwfHx+YLVvyPAmPJ0cKmoy/9to0/3HVVamn2dMT38vOg+ZigKKaCPBrlAA/ci2T\nDqvc3c7J8y6bOhV1IpTf/S6PIA2vv150MhRqUDi/xuINjnY4QrFzmzAhntf06dEZ3bMnfo7lIa67\nLl/clyxJqegTToh5UyBeV3d36qlrjkRzp0CMhtR4eHnYEHyK03M+ylDGVHg5l7Jo1FMSehFKUzQ2\nKO+9lz7sQHnNAEt2DZ5B8RReTz8dvYcqsLoHyKMNINIopgJiztcUXnbMDAovwLawsPcyeXK+ELM3\noQuETozFi4v1KF1d6YNuITnDo0sM+oCXKb088HUMJzzKiwUSBk/hp7y6yso1IjFRBZBHnkBe1Kjo\n7Y0Py6hR+X3ih+fEE2O+QyXD7CxYrsXu/ahR8VngsdC8h1Fj/J4ySTBQHFcv2V/2914SvwxNaE8v\nKe8ZlLffjvOYUVbszNizJ46FBy5q5PqVffviveVzmzUrPldnnhmpSst9WB6zry8u4vv2xeeeaau7\n7kpprdtuA773vfy1V1cyc2ZKPVodGz/XIfiLvDqZQIw8PSMxeXLRAZ0yxX+Gm0iHPceyrD6lCQbV\nHFK93yYJWS1gA3yDogVxXkI+y8o98ccfz+kv9mjNoBgd53lRc+ZEj2X06Fybf9llOaXBvLt6AFoY\nN3Gin4TT6MMSiQp9mDl5bCiTjXpo4p02Ufd4BsVLsHo1KEA6jp0UXkyl8nzy6paA+FBbVL1+fe75\nWx6F1VfsELCsGCguBFqseOWVRX7cE1J4i8QZZ+RKJMOFF5Yn5hVN5sBQKS+mnwxvvTU4g1IFjnDf\nfDNdY8zTN8yaFX+/b1+83xZRcpTICi1Ozlv+xZ4/rQ2yyIsdPK5pMXhV77xWGLx5UkZxVokzFENd\nB45IhOKhycTZu7fokZTVJfCE0YS80iL6GVwhzy1X7H22sHk8r4akWZYuJhpyHn107u14NINWMnu0\nh7egzJ1b9Dg8T6bM4/RC66aUV10O1aO8vAilLIpgcAQCRINy6qn5a1aZcZRqRY2MLMuVXe++GyMG\nizwtAmFpJTsO5lDYg3799amah2kUe/8776QLw7XXxnoIhqoDAT+PUmb8vRYsVcWNikNBeXkOX5Pu\nGR6U8vJEH1kWx3batLgerFwZ54flPU2xt29fOrZaca6dLebMSe+n5lamTo3nx/dGqTQgLYw2eIWv\nQHS+1Fk85xy/ZYunvhqqUzEcrX2GxaAMdeIA6cLlUWiacH/mmbTocfnyNJ/CCXn1eEPIqRejvAz6\nsFlFbU9PPoDqTSjNpQVol16a5ki8kJUnu8FLzHsV+GXwahzKalyGiv37i/SmJ+fVVhqetFSr5j1x\nh4ELHK1oThdAFS1YFMneqxl9cwjsMzgBOm5cPDf73ZgxKQ0GFCMSr9+TVx3fpF7Ay6OUUWFaxQ/E\nhbkJ5aVgxZ3Bo7xUgONFNgydC0x5aYRi52aFjTY2r78ev3PXrnzRXbQoGuuurnifLNdkVCdQNAZ3\n3QV885v5a2UegCjQ4cXe1huea15uDfDHwM6TUdUhQynOsjoXpdgBn/oui5KaYNgiFKWNOk2eKlhr\nFIaqiNSAPPVUnnPRm+LVqmzaFBd+naimpLIkJy/qzDGyekKpD03YXnppMTmrHoZXM+D17/LyLWXw\nPvNQVUR7f+stPAovivFg533gQKoS4mSwJXT37EkVgr296UNu3iXTlSwN54VCawxMdsyvOWL0ahKm\nTEkpLu2ywNBjnny8TOXjfZ4nG9UC3U4YbISiBaxa1KjQZ5ydDy9CYYOyd2++aFqLJKOLOPq7+eZc\nJsy0VXd3fLbM+E+YEJ8z+8wQ4jxlo+81Br3qqmL0aZ0ZGNo7DCi2c+LrVINUlnPx4KnFDpXa85BR\nXp7CyzMyqksH/JYrCl2IuHqbezsxWIps79cFyiKSM86Iizkrf5jTZumh0lxaIe/1SfISc4oyfboH\nr0FgWS2KB110DxeaSoa5/1oZuKgRSItHNQdlSXg2+qzY0YjuxhvTfk5qQDw6g+XIBm46ajjttJjj\nY3h5lLIWLB6GWi0/lKS8QsUXuh5UFTWWRShAOkZdXXHteOKJXPrd0xOv4cCBtJ/ehAnRENozpsl5\nrRnR33utWDyKyzMynlKrLELwVGWnnlqcK2XwpMOHSu15yAyKJw/dsKGYkH/++WJCXjlzz/th6CAw\nPcYGS+tWPFiEYt4CL0b8cJ9ySjrIqvpQr0IrZL08inVGZXjhqld3okbMzqlJlWzdnj8e6uRavIdF\nDcquXenCpHUrSl96UIPC98tab1gEesEFMY/S11deT8RjrRGiFzFyZwXAf/g1igX8PEoT+fjRRxeN\nWZNKaQ9161C8pLxi27Y0QtG+XVUGxVOMAnFO8RhY09d9+9I5zZQytyvipLmJcQxaM2aePs9jHcfu\n7jhf2bkzh5L/rkw+7K0BXoRT1mGDVYqG4YhS62JYDAoXIBm8KvkXXqgnGdZeX9qyXjuSskQYKE/I\nm8KLoQ+M0RM8ubw9UZS20ESrSj41z+Jx69oqHfBpL89j8RJynpqoDF5CrgkFVue9urcFkFOPhnXr\nUrpTDQ43ieSFjb/fi1DUoFivJ45M2BDwz0ptcG8oII4Hj5G+35OP2v4aDC+iZEOnqCMfP/bY+rkZ\nb4Hzvqcu5aWo03alihLzztcoLxPFzJ8fPX9T8hlNxTmS228HvvOd+LM+t1w4qvumAMW6M0/ZpWIN\nwO807tWpeB0XmuybNHdu8b1DdSqaYFgMiiff9RKpHr2xdWtnDl1b1utrLXBkg8TJ+/Xr4+B0WvyM\n+7WHixdrzn1wVKB1AEuWpKqRBQuKDeQ0jzJvXnHSeYl5L1z2WnU0aew33Bssed/rSYY1v6VdhTVa\n3b49NxZMbXKUbAZl//742ebNjRkTHZ8tW/IkK0cjLB/mVuXKN2sbczUgnrRbG0oC5Z1g9d55UarX\nNcFzKsqSuh68OVA3Qtm7t3NLJs+g8LNfp9OwnpsZlLffjjmW6dOj42j5E6tW50LVE07I7113d7xu\ni3BuvTXdB0Vbr3ATUSCOzdixaZTq0ZSqEAR8tZ9HezbZN8mbA5rDO5QYFoPiweNUVVFlUGWHVjtr\nso6bPgJplKMPI9dHmKHh8+D3l1WV8qLOxY5sNDSEVcVFd7e/uQ4bGS9n4hUreZOmSbLeQ9UGS4wy\nA6WT3ht/r6hR8eKLKS3qUVz2XWU1KLwXyujROeV1zDHR8GRZGiFa9DB7dh6VqOqKI5aJE1NqyVPX\naHR4+eVFT9ZrNMjdcA2epNxLzHveaROUzYG6vbw6OR9qdDRC4SinjiNk32dzf+bMaJT27ctpIpZT\nc3sU3kqAIxbLq5oRtnyVnU9fXzQifJ+0FYtt/MVrCTebNXjPbNVeJdoJ4ayzimrBJl2qDwUOmUEp\nQ6eLK9uml6GFlOzlliXkGWXFddo2xQaRF3DmYk84IeXKdYCVRtI8i7fIqMzTu19NW0orPM+4bg7F\nW0zefbfYn6tJUSNDZcUqIeb7oQbFcnYmVzWDYpSXRzMBaZ8upqbYA7z55rRaWtuF62KuUYx1I+aF\n8pJLinvheHkUr7jNi2aruit40EV7KBGKwuvZpvDqz+z7PBrdA1NeTBmOHZsv6OZQcNdoNiKnn55G\nBTfdlEac552XUlPaUFL/Hojjr+3rZ88uOoceLa05OKDYtBLw+7qV1bgcLhx2g8LwGkWqQenEzWrN\nCkcvPKm5/Yp5ylYgZ/+bx2BUlCm9WH3h1ZDY71TLrvyqtrP3PFtv8fASbU2gC4dXCFm3qMnr5Va3\n7YqOt6fcqYIaLqZQuGDSKEs7Vxs/Myg2vmYwuLEjL/K2SRNQ7Cyt7ThYjgr4+QxdeLz8CNMxBm+B\nGaqIwmtkOtQIhVGn2NkoSQ91NmFTysscNjMwNs62iR13jbZnwK6N8xRaBqBj7bWb1zyFV5Cohgrw\n5cNeXy+P+i7bS+dI4rAZFK9Y0VN4cT8uoNhhWBO7q1enG3Pxhlus6tq4MU4aIPeU7bvM2FjF+3HH\nxcWDJcOqrLKJyHSEGhtVZ3AhlUE/11RHDE9z7u2rYIV2DK8ZpVfUVJfyKjMoSmWWKXbY49WcWiea\nQ3u3ce5OFUTeuZpBMcGFJUWZemAnQAsQuSZE76EX/ShN4cmHPVWPXZthqJSFJykvcyq8CEXHpU6N\n2VCLnT0Bh8LOzXJkpsKyBLjVE7Gikj19ppzvvDOvSQkhrgHGJvT3F2kuZRv47+0zVHU3bVqxHsSL\nNMsobaXMmqJO3dNQ0cig1P1yr/pZ1TyAr/BSqEHRgsaqhDzvgcIKLyvEtP+NMrNBM9509uxc6cXG\nhXluVWXww6seqFfQpuG0ZxTOP983MnX4U0855BU1eYuJN95DiVAUalC08E3buWjCXs9VHZYyg2Lj\nxyobo7o4ia3qLPVUlcKwGiaDdij2oltvfwxvHL2Kdy+q9JwFrwq/rDngcFFedToNV4H7eHnriXdu\ntneHOWAWefKYckt5piWNLbAcpibjb701bRipv7ctn/mZUZqU38fXoKKLJoWvnjLPcyC8noLDYaQU\njQxKGX+u+nCvIlY7yAJFI+N9libyV6xIuxBzbYvebN7lb/XqVC0UQv7ZFqFYqGwKGqa6OOTkZKq2\nHldeXPvz6GZKXh5F1V8zZhQ9m7pKL68Wpe5iAhQfWs+geH286jQM1UahXlNIHjPu8VXHudFzNWrE\nDApHJtxagyuiOaehdQpKc918c6oQ0r1zgGJk4G3d6rXgKGvDU8e79WSjZU6FR3vqva6iqgxKeZUZ\nBf5MLjBmyosFHvv2pd/tzQMbP75Giyi4Xc7o0dFZsMWclXgW2dtirSIIcw74+7VY2RNTeJ2LvbYu\nnhjHo0M/vWzmAAAaTElEQVRPPLHYw62sJk0j4bJq+aFELY0Mirdboxfa1q1BUbnxqlV5/UgZdD9z\n/owNG3JaC0g38/GSfHbuHKG8/XbqaXJ7bDvG+ycAqYehPbm0SlbzLJ76o84GOh5/6nm2nkEpo7yG\nkkNRb1THds+e4tzRsexkULh4bseOzjULeq52Ph7dw+PGkYUqsbgHk46ByoW9+hPddMkrcOPI2OAl\nZb15UuZU6GLShPLy0ImGU4Oi64QuWrp3TlkfLx7TMlm8ntvBg2k+jO+bdpfmfMbChenY6yKvLZW8\nHIkm3dUpAfINwRhePYonR/bmgOek1K2WbyIz99DIoHi7NXqhrVeD4rViUWjLlV27UiPgJQc5otEG\nkp2gEYqFgKwQMR6cJ6lOWG5HrwOiVdLTphW5dl10vN5fHuWh8NpP657YwKHJoXTa/nfTps6S4Zdf\nTiPWKsUeRzc8L/jnOpvBTZ4cx4MNOxtm3dNk2bJUvaVtUHSXRI02PM7ciiwVfC1eU0+vBYtnjOou\nJnWT8nWgBkRrUDRXqkWNZX286hgUO5ZlMUpYvTqVhN9wQ67AYul/T0/8Tst76fYF2npFO0lbxMPs\ngta1AMXciieQ8CIUr9eXFtUCfpRaFqEM954ojSMUz6DUabuiHLe3xa8uINqyXvce1024OCejE9SM\nACcUOYfClBfDlF5AKrflCmRt0cAPq8eHaiJevU8v+vAWHY8/VWPn9QfyvJAmBkUjDY/e0vOoU4Oi\nc0S3i2VwvQrTohyJvvNOeXNKO7/zz89pKR4Xph1Z1aOFZx7txa81wgmhOP6eh+rRFkDnZH1vb1FS\n7i1YQ0nKV52TQZ3PTlXyVUWN3m6Ndm4Kqy8yA8L339SSfD96euI/yznw5lq206etCdb7y+5vb29c\nS3gB1mJFj5LyxBm6BbC3bnh1al4LFa++xWuNX7YnymE1KHUorzqSP1VnAcUFRSOO5cujCsrQKSFv\n0Q738+K8jYXSRnl5g8M5D15YeKLqwGiVvLZd0PYNXkO5LCsm7z2eVb3RwaKuQfGiVG2wCRQXGXUG\nmnaj3r073Udn48b887yiRqA6QjGKj+WdLJDgMVKaylqlA0WaS2uTrGaIFwKNQD3aU50U+2yvKeRg\nIghPslyX8tJjZbsu8vtUhadV8nUbQ+qY6rUbNW2JeaZ8OUfC+Q2uftccp6q3VPp/661pUaM3btOm\npXPE275ABRyAT2FbPVNTeL2/vCh1qJtsHRLKSw2DN+FV4eW9R1tWsyQYSL0a/ftnnskVXiwf5uI6\nMySaEGSwkocTo0pdsNepvObSpanXool4z6OYMyddmCZMKNaseIl53tq0Cep2NvYWaV1gbItWhqq+\nNN+lsnJ9rVsacGt8rkGpa1BsPwl+gC67LB8X61QLFMdHPUztl6W8uUrBvToDHTeP8rjwwmKTQK8F\nizojdTFYymswnYab9PFqQnmZAfHoIa43YoWXOmvsFOiirrkN3f/Ga8Vy++3FSnpVYml/OO+7AL9Z\nqMes1IGn/DrslJdGKHUKmHQyAamHCdTj2L19TQy6QPEkZ8kwUy+8x4Z+roV+TBlxTYdSSVxtrdv5\nai0Ie7gGLV7kyV8GT0XibbwDDG6BqSsb1vd5NSFKX6lkWMdfIxpvkzRDkwjFHmSvxQlXGaunxqIM\nLVpTmkuVPCoN9iLhJUtSI2OLJd9bpl8NXlLWq3eog7qyYUWdxpCq1lSDogINhkYotgZ5BsXGia/F\nDARHZeyghZDWdd15Z5or4cW+qyu+lz17lY+rxNiT7muXYk8+7EWknlPhNZ4s6xPH8DpuHPEIxTMy\nCq8ppC4wXkKevR71WHUSVyXk2YNlwYBNRi9C0ZoCD9yaQ6MO7hUEFD0STbKzQQLyLUwZyrPynvcG\nT+UzeXLRE6mLOrJhRRV9YdA50UnhpSpBXkhsp0ags0GxWhRWwjC/zJ2Jp0zJc1ns0YaQ7nGjC8bJ\nJ6dj720NrInT884ryoLnzk0pTauzYHhJfq9K39syVuHVLwDDE6Hwts2A3yWj7Pu4uzivQSqUOXiw\nWOcBpPQ0V8Tzs3LbbXkrFnte7Dx4/3kg1qA88ED++s47090dvb1uVEDhtaT3ipqB9H54+VWv8NnL\n3dTBEY9QFF64XWcTpZUrY/ho0Jb1updJVYNIVX8xBaet8QE/QmEqgSueWaHBNSeqwvF2beT6lIUL\nUwOkcmPP+/Kko4qh7ovSCXUorzpFjeqxqkFZuzY1KOzBav6FnRM2KLybo8EMChtjvq8sfuC2GHoP\ndXxVTXf00SlFqdsQKG3mUY4eH6+KQC+q8FQ+Wi/loSxfovNQX9fZC0XhdSgv+zx+fstyKBb1cxNW\ny5PxfjYsH2aZr0b2vJXE2LHxnpujYWuDfc+kSWlRJFBU1mkDya6u+I+NvOZeAb/7sEJ3krTr8TpP\nd4pajniEovDqPVS95X2OhrydWtQ//ngakbChqFOBz+gUoTDVwNXVyr3ygB1/fBrhKIWlhWoeDabt\nVbweQlo/4im4ygzKYGiwOpSXRihe2x0gXVBUyFFFpXIbHYVGrrpoeS1SuMkee7NqBNiZUA9TaS+V\nF2vbc4+W0gaTnpdZp0bJi64HuyfGYCivOrUhVdCaFEZZDoXzXHbMIgX+bu4aoLJZfk6uuSaVBN94\nY9p+/txzi/sdsTyctxcG/G141Yn0up17eZQyRRijrCbNa/PDeN9FKJ5kWHdbq7NzoibktY8TGy6d\nwKzw8qrvDVU5FE60svKDDUp/f6qUueQS4De/yV8zRaELvae60G6x+nleWxYvvFZ4BsWra6gDpSq9\nxUNzKHVqkID6i47OhSaf4xkUHhtu2KlFh6z2Mg/Txl8Vd5rv8BqB6n4WWvfgXYfHoatT4VFjZYtJ\nHaeiKeW1e3exJKCJ81KlEi0zKF5imqkgXoTZ8eP8xx135NRWb29cN+y+XnRR+iyqA6ERq9Uj8XV3\nSvADRcruuOP8HRy1wFFLCDypeB2n4n0XoXhV8gqNILyknhdmly0UmpDn16wO4vYP/LMXoXDCiovF\nVPHT15dTJ9p2RVtv6EY33GwQKFbCeju6Aem98ZReCv0ewC9qGgw8J0NVXtq2vixiKYMqCcuKGuuA\nDQrLJpmq5XwDS4k1OlCvVNVdqvbp1OvLS95qVwbP0/VasCjKnApddBR1IxQ2KHWEOgylMKsaQ5bV\noaiDlGXpvbv88jwa4HFjw6D39pZb8qLEENJcW39/PGaGqbs7Pgd8DlzwDBST9SYS4DmsEnX7bn6P\nx1Rw4WYZ3pcRSieD4lXJK3QTpWefjQujQSvkNcrgJCwQaSfOpzDnygov9pz5O8rUYzaIqtpiMHWl\nSXLerxooFrApZ6r1CZ4aSKk0rXsAIlXGdTGex1rWx4dRZ7H2quQ7FTXy3iVAcXy15kSbQvKujV7b\nlyqwuoeT13xfrUstkLbK0aJE3btE1V2aJ1m2LP29t5+FFj16+6Moyvad5znLLWMMXh2Ch07zQB1C\nrU3TjbW0F5hGtJ0oL/4sjlBsMbe5zfOQ79HFF+eOnya5WW6syW5Nxuu+KBzhAMUoxismVafRE+No\nbYy3742XmFd4EYoqyzwWpAmG1BzSm2hKb3jRhy7gWvGuHYZ5G1+gWOBYtWMjNxRkT5mTwhyhlKlh\neHIec0w+CZWWYvXW+PFpiwWNJlQr71WpahJWcy9eJbyXkFPUXUw6eaheY0jFyy+nBkVFGkpheQn6\nsl0b1Th1Onc+xtJhTsyzYk8jUv5dT09ctO2B1MSuepLa9TeEorFXOsMbSx27spyJfpeijlNRJ0LR\ndUEjSi0bUIOhVfJVlJeqvLwcivbiy7L0Pqtclhdajkr0s6ZPj+dm36njq8ZgzJh4b6q2+fZ6u+ka\n5NUtqZPr1d0oFerRWUrFelR8EzTeD4UnmMeV6o57dRRemsTt1LJeJcZsoLZuTSMkTvazp8xeFP89\nT0yesKz04kVdW0BzaA2kIaS2U/DaK7AcFShuE6zeiv09f4YnG1UMF+VVJ0JQr7KTZFgVXirq4DnI\nkvCmYEqEH3TeHhhIvXtV511ySZrTYCFFCMViQTUGF1+cOiS8d4d9hkI7LXjvqetUdDIoQL1Ilc9B\nKa86RY38WtcQBotvynIoTO8de2xO97Kzx/eHjYjRgLZY33VXKglWB5LpcKDYmFEr67WBpLfIK/Xt\nKfRU/eWVEHiJeUVdx7IuhrTBVh2uVBePqiS5QWkUlaFylbSXkLcKeQV7yuxFcYTCBoJzHqz40cI2\nDht17xLdtfHUU9M6kVNOSQdd1UMcngO+B6Ht7b0qayC9V3W80zrQsaqz+CgtylEkUKw5YYpEP5+L\nGpvmZthT1VYk7CVywtXqVOw8li5Nk7E33ZSqu665Js2T6NawixentKe3wGgNw4IFnZtJDleUWqeX\nl8IzKNp2pcqgMOqIaoCU8uKogqkgNgbc2FGjOZb4WzLf7sF116U0pjaM/NjHUtpLowtbX9QB5MjG\nS9YrhenVsSjKipwZw7UOGIZkUOrsyubtecHeZ9NFQCe3fj7nY7wCKUseMuWlEYoZFObVOTJQxRZX\nqurCpBTGlVemC5C2ZdHEvnrLQHFyKX/qGR1rbGfwOPXBQCmvOvugeI1COalb5aF6i5O9VmXgULoD\n8N4Wen/5NW8rCxQ9Ry1WU1rE8yx1odc8ivc32oLF2wlS581wUV4KjVo156kGpIri6jSfPMqLlVLz\n5+d9s9jZ0wQ8C1fUCWBDZA0hja2YODHS2uaYjB8fz8Ned3fH8WLxA+fogLRoFvCdCp2DWp0PFBPq\n3g6xivddhNKJ7ti9O82hqMJLE66akOcqWaCo6NIK+T17chpO6S9GWQ6FDQp7hhraAvlk1uaOvFhz\n5TVQbCaoyouypJh6nxzl8ENTBlX5eG0XBgNdPDSaHMpmPd7fq6CjrKixLniR5aJQzpUoNcnFcUCR\nlmT6JYRiWx198Ll6Gyh2rPUStWocVIHmOWnaAqiuU9F0DDXp7lFe7BR4xcYGdRIYHD1xlT/nFbnx\npjqCPA5Me/X1xb+35/b669Px1sr5pUvTKFNfq7qL2+cDPv2s3YG1yNUz9NpXzKtr0fzMcFHfhmGl\nvOrQWevWpTs3aj5E8yda4FiVkFewwkvBVA1HKLwYsNenSTBWVynHqXJSle2qkVGlhfLumhPRNi3c\ng8qgeZW61fJle4uUQUUX6n02jRr096+8koo8NCJl1DUoLHRguoMLVdV7Z5WQdgfWfTGU1lJ1H1fg\nA3GB4aI59Ua9xUPnhNeCRaEqH6/jsEK/2+ss3QlMUXuvq8B9vLxz03xkJ7BxZ5Wl9s3iKKW3N16z\nPZOnn55GoarW1Ne60I8dW2zFo8+nJuunTi3WT2n7Jk81qNDI7H1NedUpYPOKHJkC8wwKS4Kffz5t\ne8/GQPe/qDIo7EVphGKDpElLnrTelqys/uCQVhNzWpyme51oolar7Osk23QC1jEoKggAOj+oSl9p\nhKKS4ddfTxd9NUhqQDS/UrXpVl2DwrkxVnrpZkVs6LVwjRfncePSqHLevDRivOSSVA6uHRa8nR6B\ndMFRY6CJW68FizoqdeoQPPB51NlPvsnndfpdpwilDtiBYDWf0pFc/KgSfpUAM0vQ1RWdDFsrurri\nnGb1lHak1hyISs69HJgqArUVU53+gzqPdLOvoWLIBoUHW6vk60QsKiHWCaQcKntIOvm0H5juS62K\nNO8cmPLSCctKL1VbzZ2bc5o68NqaWhP13CodKBoUTbJ7UmENZXVCehvsKOruiVIFjVC0qFET7kp5\nrlmTvtbfs0Oi96CuQeHiRjYotpWBgWWhKq/UKITHOIR0rnR3p1X1Jg3l6EArvbWTtOZRvG1eFbrA\nDKanmyblm+5j0xRaJFkVoQDlxonP+7TTcsqQjbv1T7P38aKubei1SFFbq9x9d5qc131U9P0agaix\n9+TDKlNWp9YzsEqv6vd4a4nC28ivDEMyKMp9apW8epP79pW3n68DDbeV/uCWKwa7ybqwMcpyKApu\np6E5CN2Zj2ku9Tg16aaUlS7s3kTRAjhVj6lBUZrNw3AYFM15aYTiSYY5AtGIVecUT351PuoaFI7e\n1NAytamGn705TYLfeGOaWL3hhlQNpDSoGghtca8RkfYVq7N4aIt+rxOvIoR0Xuv3DEeEUgWtUamb\nQ1FMmZLnBphC1LoifqZPOintnaZRCQsdjj46PsO24GuBsfZhU1rc2jbxvZ44sVpSrnOmrH1SVbuX\nwfR08wqjyzAkg6JQOeiLL6b5Et0kSdUdO3em3om2YFC6SxPyWkDHqNqCliOUqr3btXiILbd6jEuW\npFGH8q66YZYu5pqc1x5U2m5Diya9jZc6waNNBgNehHSRVydDDYzSpkqpMbQGhSuyq7wuLkzVBZMN\ns7bK0WQrRxUawSjtqYuBGiuVAitl6RkQLVxTrl6dijreqKcgq6K8uIVRXVTVsumaUDeHopg5M4/O\n9NlkKvCaa9JIgeuI9H6qRPi661I12IIF6ZgrharPqNY0adSrkuOyfAffA41KdQ4MJgnfRMAzrAZF\nw1WlNzSBrhHFk08Wt/zlgsbHHuuckLfJqhRX3QhFHzqOBFTpxeobfaiUElPljtYfcLtsoCgfVk5X\n613UG/UmgZ6jJwaoarswWMUWf69Kgj2JcNkCpa1/uAZFUeVJew0iDZrf4Hukvdm0Fxertayo0e6n\nVtXb3u/mkHR1FaMD3Z9E25SrukyTslOnFvfO6IRO0bHeV1Vxdtp5c8+etMhVKVJ1Iqsk5FWGjBdW\n3RrgrLPyHBdHMkBsj8OFh7wtgRn5MnWnRjR33JHvseL9ntv6AEXJuefgKQ02e3bqOKraTyOUwcjA\n6+ylYxhWgwJUe6fKj2tLFU3Ir1iRtmTZsCH1bjmy8FquMHVStcVo1U6QrAJSY6McZpXna91HDerN\naAGjJu20X5Py4V7tSScDoF61Lib69++91zknNpzQxUjb1le1XVHjw1SWbn3KToNKss85J39AbfzL\nxBc335wmVtWDVadBjZeq95RG1ShHk8rqndZdPHicvUWsKkLxch5MUalwR9useAalrI9Xp3NnVOWL\n1BBwpKlMgvbn4nseQrGLMZA7AePGxXtpz+XYsfF37FQAKRWt7ZdmzEhzHtr7SwU+qvYbDtbhiEUo\nHnhSq1xQa0527UrDW29bUE8BAxTpLlV4MXWiYXrVnvJapcxQ66/7Gxx/fMqr8sKlXpM2fFPFlSaM\nrbV6FZQSUDpPi5q8HAp/R9nmRoOF/r2KOHQbYM2ZvfJKeYW1nivfe3UMWFCh8nBthcLJck2ua0Sg\nuzB68mE2OBrF6qZsmqhXw+jNB2+M+JjmFDRK1TnkGRSOULSUoFNRY5MiR0UV5aVRKHeC1oaQqsLU\n0gE2TLqvjTaM1ChEJeIqKdbv1mT9ddelr3VnWBX8DLX9PFDMzx6SCGWoi0fT79Dv0933OiXk16+v\n3oSJDVlVhKLdfVm9owuTRhGaR9GqeS1o0wWevwso7rimD79KAGfMqC6iVD61kzeji7Qm4NUg6Gul\nO9Qb1XyKKryqiho7nWuVKEHpQn6AdGG6+up0gdDebSwPDSGlrazuwwyW3m/Nm/X1pY1B6+RAgHLK\nFig2LFVpfFPKS7scaERSp+0K/76K4lJ4BoU7hPPvmObS9yqVqUaDK877+uJn25jafiVehwWgGGVq\nWxU1EBqBqFDAKvMN3k6fdcD3RhvQqlE6JBGK0h1KR6jXr7/XKEDlhyoPVu9UDYa2rFfJcFX1rXpR\nTSIU3TippycfDL3x2vpAE7Gq5FEDpB6qUiJKm+lDowk5pQE0ydcph6KLtFbJ62KxaVOa4/Da8FTt\nI6+/b7LPhmdQyrwspbnmzUvHmB9arTZX9ZZ6qKre0saOyptrVKt5Ey2S1d+rGEM3/lKnok6U2oTy\n6tRpWOdI071TGGrsqhwiTYhznyulMs88M32OdPtezbPoc8lj2NWVbifR3R0dQ6O1urvj3DR2IoTi\nAq9Ogb7WTgzqeHrMBD/3+lrVgockQtGOsWoA1NtUxdeGDelismpVfHANXodhTsDrFsCq+PLoMYPK\nlXU/86oIRVUvqvTiPl5AGlV4cmEeKK3O1QhHi9d01zilVLyHpsqg1FlMGLpIax3Sq6+mRY0qhKjT\nZZhf79hRlAGX0Xw6xlWUl0Jl25rb0Ps6Zkw+xhZF2BhrHkM5bo1w1BvWqnk1SJpH0ahXk7KDcSqa\nRChNKa/XXkt/X6XiU3hlB2zsvF0bDZp818iAnbEQ0jG2fnplEY324lI1mL7WFvlKc6m6Sx1Lzd1q\njYzu4qrUvMdU8Bw4LBGK7taoE0frBXTxUEVWp5Yrzz2Xvl9D5yojoNGR9vRSg1IVoShU6aWDq9px\nbZGhtBSriDRM17yK7rCnD4km/nXx0IaBSnk1NSheV2j2PjvVoLz8chrBcB82Q5kBUapDqZYmlJdC\nqUhdxJcsSR94lXCz3LirK41izUGx81fvcNasdIy8vm1VcnFNytZxKqoMClAdoSjlpeuCzpEqmrIT\ndEyV8vLqMsooQn2Odftl7d/FVFQIKZ3c1xevyRbhSZOiw21RrdaoqBOqAhxPjMHzTX+vTou27tE5\noLUo6lgetghFvVOeOF4NCi8emiTXegSvIr7MYOgk0R3+9LM1Ye9FKHUnuS76rHcHipSGyoN1F0dt\nn6BbvupDUlW16smC+Vy1d5Py6U1zKLpYaISiKizNXTVZXHTMVcChnq9GrGpQqh4SXai0JkTVVbr4\nLFuWeqyaeFWeXR9w3UOHx1zHVKk8HUP1PpUi8xYTjVCaJOWVArPPGAx0zDsZFBUYaP6LnTc9J1U8\nakfgG29MI0fdJ+W229KoQ5PvakS4oLZMYm5jrk0t9Tp1fuoaokWWnXKpyqQ0iVBCVpHlCyEchlR8\nixYtWrT4oCHLsoKrUGlQWrRo0aJFi7o45HUoLVq0aNHiw4HWoLRo0aJFi2FBa1BatGjRosWwoDUo\nLVq0aNFiWNAalBYtWrRoMSz4/ynXzfqjwCZCAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "loc = anhima.loc.locate_interval(pos, 11000, 20000)\n", "loc" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 5, "text": [ "slice(204, 315, None)" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "pos[loc.start-1], pos[loc.start]" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 6, "text": [ "(10999, 11044)" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "pos[loc.stop-1], pos[loc.stop]" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 7, "text": [ "(19918, 20017)" ] } ], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant counts\n", "anhima.loc.plot_windowed_variant_counts(pos, window_size=1000);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAbMAAACfCAYAAABz9Dd+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH6xJREFUeJzt3XmcFeWV8PHfYWkB2Wx2BQXBBZSoKKJGI7KoWU1eNTEz\n4XVMPia+JmbGTIySmUSJTnTIaNRXHT/jFqO+jkbjgoJhbY0OEXFjE3HDBeimoREakLXP+8epsi9N\nL7ebqlu36p7v59Mf7q27PQ/ddU89T516jqgqzjnnXJq1S7oBzjnn3L7yYOaccy71PJg555xLPQ9m\nzjnnUs+DmXPOudTzYOaccy71Yg1mInKEiLye87NRRH4qIuUiMktEVojITBHpGWc7nHPOZZsU6joz\nEWkHrAJOBC4D1qnqVBG5EjhAVa8qSEOcc85lTiGnGScA76rqx8A3gPuD7fcD3yxgO5xzzmVMIYPZ\nBcDDwe1+qloV3K4C+hWwHc455zKmINOMIlKGTTGOUNVqEdmgqgfkPF6jquUNXuPrbDnnnNuDqkpj\n2ws1Mvsy8KqqVgf3q0SkP4CIDADWNvYiVS2Jn6uvvjrxNnifvb/eX+9vsfe5OYUKZt+lfooR4Gng\nwuD2hcCTBWqHc865DIo9mInI/ljyx59zNt8ATBSRFcC44L5zzjnXJh3i/gBV3QL0brCtBgtwDhg7\ndmzSTSi4Uuuz9zfbSq2/UHx9Lth1Zq0lIlqsbXPOOVd4IoImnADinHPOxcaDmXPOudTzYOaccy71\nPJg551wbPPsszJiRdCtcyIOZczHZsiXpFrg4zZgBzz2XdCtcKPbUfOdK0Y4dcNBBsHYtlJUl3RoX\nh8pKaOfDgaJRiIume4rIYyLylogsE5ExXs/MZd3atbBxI6xZk3RLXFyqqvz3W0wKcVxxCzBdVYcD\nXwCWA1cBs1T1cGBOcN+5zKgKakKsWpVsO1x8Kis9mBWTuCtN9wBOU9V7AVR1l6puxOuZuYzzYJZ9\nVVUW0FxxiHtkNgSoFpH7ROQ1EbkrWKvR65m5TPNglm1btth5UVWorU26NQ7iTwDpAIwCfqKqr4jI\nzTSYUlRVbap22TXXXPP57bFjxxbdWmDONaWqCjp2hNWrk26Ji0NVFfTvDx062FRjt25JtyibKioq\nqKioyOu5sa7NGNQsm6+qQ4L7pwKTgUOBM1S1MqhnNk9Vj2zw2kyvzfjUU/bvOeck2w4Xj8svh4oK\nGDECHnoo6da4qM2fD//0T5apet11cPrpSbeoNCS2NqOqVgIfi8jhwaYJwFJgGiVez2zmTJg2LelW\nuLhUVsKoUT7NmFWVlTYyGzDAk0CKRSGuM7sMeEhEyoD3gIuA9sCjIvIDYCXw7QK0o6isXg0bNiTd\nCheXqio491x44YWkW+LiEE4zdurkwaxYFKKe2ZvA6EYeKul6ZmvWwCefJN0KF5eqqvqRmSpIoxMj\nLq0qK6FfP+jSxTMai4Vfv56Q1avti27r1qRb4uJQVQVDh1qCwKefJt0aF7VwmrF/fx+ZFQsPZgmo\nq7Od4ZBD4P33k26Ni9rOnbb6R69etqSVnzfLnqoqG5n5ObPi4cEsAevXWyrv0UfDu+8m3RoXtepq\nC2Tt21sw8/T87PEEkOLjCw0nYM0a2wmGDYP33ku6NS5q4VE7+Mgsq8IEkO7dPZgVCx+ZJWD1ajjw\nQAtmPjLLntxgduCBHsyyRrU+AaRXL1sBZPv2pFvlPJglIHdk5sEse3xklm21tVb6pWtX+7dv3/rl\ny1xyPJglIByZDR3qwSyLPJhlW+7vF/y8WbEoRD2zlSKySEReF5EFwbaSrmcWjswOOcQC244dSbfI\nRcmDWbaFyR8hD2bFoRAjMwXGqupxqnpisK2k65mFwaysDAYOhJUrk26Ri5IHs2wLkz9CHsyKQ6Gm\nGRuuf1DS9czCaUbw82ZZlBvM+vWzSzF27ky2TS46YfJHaMAAXwWkGBRqZDZbRBaKyMXBtpKuZxaO\nzMCDWRblBrMOHSxBwL/sssNHZsWpENeZfVFV14hIH2CWiCzPfbDU6pmpejDLuoYJAuFU46BBybXJ\nRaeyEk44of6+L2kVn9bUMyvEQsNrgn+rReQJ4ESgSkT659QzW9vYa3ODWVbU1NjipJ072/2hQ2HW\nrGTb5KKze7f9jvv0qd/m15pliyeAFE7DQcyUKVOafG6s04wi0kVEugW39wfOBBYDT1Oi9cxyz5eB\nj8yyZt066NnTphdDngSSLT7NWJziHpn1A54Qq3/RAXhIVWeKyEJKtJ5Z7hQjwKGHwkcfwa5de34B\nunRqOMUIHsyypmECSL9+th5nXZ1dRO2SEevXp6p+ABzbyPYaSrSeWcORWadOliDw8ccwZEhy7XLR\naCqYLVuWTHtctFT3/h2XlUGPHjYq79s3ubaVOj+OKLCGIzPwqcYsaTgFBT4yy5JPP7UD0PCcd8iT\nQJLnwazAGo7MwINZljQ1MvMyMNnQMPkj5OfNkufBrMAaG5n5Go3Z0dw5M230AhSXJo2NvMGDWTHw\nYFZgPs2YbY0Fs+7dQQQ2bUqmTS46DZM/Qr4KSPI8mBVYU9OMXqQzG5r6svNrzbLBR2bFy4NZATVc\n/SM0dCi8/76l9rp0a2xkBp4EkhVNHax4AkjyPJgV0IYNlgnVpcue27t2tdReTxJIPw9m2eYJIMWr\nEPXM2ge1zKYF90u2llljo7KQnzdLv7q6pq818mCWDT7NWLxaDGYiMiefbc34R2AZtno+lHAts8bO\nl4U8ozH91q+3ZI+OHfd+zINZNrSUAOIZq8lpMpiJSGcR6QX0CUZT4c9g4KB83lxEBgJfAe6mvqZZ\nydYya2lk5kkg6dbUFCP4tWZZ0dTIrGtXW8qqtrbwbXKmuZHZj4CFwBHAqzk/TwO35fn+vweuAHJT\nG0q2lllzIzOfZky/loKZj8zSra4O1q5teskqTwJJVpNrM6rqzcDNIvJTVb21tW8sIl8D1qrq6yIy\ntonPaLKWGWSvntmaNTB4cOOPeTBLv+aCmafmp19NjU0jl5U1/nh43uyIIwrbriyLtJ6Zqt4qIqcA\ng3Ofr6p/bOGlpwDfEJGvAJ2A7iLyAHnWMoPs1TNbswZOPrnxx8JzZqp2ga1Ln+aCWf/+trK6V0dI\nr6bOl4U8CSR6kdYzE5EHgf8ATgVG5/w0S1V/qaqDVHUIcAEwV1Un4bXMGnXAAbDffjaN4dKpuWDW\nsSP07u2rRKRZU2n5IV8FJFn5HCMeD4xQ3ec8nfD1N+C1zBoVjs6aO/pzxauqCg47rOnHw/NmAwcW\nrk0uOk0lf4R8ZJasfK4zWwI08xXcMlV9XlW/EdyuUdUJqnq4qp6pqp/uy3unhaqNzJoLZp7RmG7N\njczAk0DSzqcZi1s+I7M+wDIRWQBsD7ZpGJxcfjZutKmmrl2bfo4ngaRbPsHM0/PTq6WRmWczJiuf\nYHZN3I0oBc2dLwsNGwYzZhSmPS56PjLLtspKGDGi6cd9ZJasfLIZKwrQjsxr6XwZ+MgszVSbvwYJ\nLJjNnVu4NrloeQJIcWsxmInIZuqTN8qAjsBmVe0eZ8OyJp+RmS9plV4bNtgC0p06Nf0cv9Ys3Vqa\nZuzVCzZvhu3bLTPZFVaLCSCq2lVVu6lqN6Az8L+AO2JvWcbkMzLr08euQ6qpKUybXHRammIEn2ZM\nu5YSQETscR+dJaNVq+arap2qPgmcHVN7MiufkZmIZzSmlQezbNu92w4y+/Rp/nmeBJKcfKYZz825\n2w677uyz2FqUUWvWwJgxLT8vPG82usXL0l0xySeY9ehh6/tt2mTLIrn0qK62hQ1aWr3Fk0CSk082\n49epP2e2C7vQ+ZyWXiQinYDngf2wc21PqepkESkHHgEOCd7r26VwrVk+04xgwWzFivjb46KVTzAT\nqR+deTBLl5aSP0KeBJKcfLIZ/6Etb6yq20TkDFXdKiIdgBdF5FSsBMwsVZ0qIldi9cwyX9Msn2lG\nsLUbb7kl/va4aOUTzKD+WrPhw+Nvk4tOS8kfIR+ZJSeftRkHicgTIlId/Dwe1ClrkapuDW6WAe2B\nDZRgPTPV/EdmY8fCyy/D1q0tPtUVkXy/7Py8WTq1lPwR8mCWnHwSQO7DFgc+MPiZFmxrkYi0E5E3\nsLpl81R1KSVYz2zTJivc161by8/t1g2OPRb++tf42+Wi05qRmQez9Mn3YMUTQJKT13JWqpobvP4g\nIpfn8+aqWgccKyI9gL+IyBkNHi+Jemb5jspCEyfCrFlw1lnxtclFK98j9wMPhHfeib89LlqVlfmd\nJvCRWbQirWcGrBeRScD/AwQr57KuNQ1S1Y0i8iyWCVly9czyPV8WmjgRLrkkvva46LVmZJbnvumK\nSFUVjBrV8vM8ASRakdYzAy7CyrRUAmuA84NtzRKR3iLSM7jdGZgIvE4J1jNr7chs9Gj46CPbgVzx\nC5ey8mnG7Mo3m7FfP0vj3707/ja5PeUTzH4D/G9V7aOqfbBAdk0erxsAzA3Omb0MTFPVOVg9s4ki\nsgIYF9zPtNaOzDp0sESQ2bNja5KL0MaNUFYGnTu3/FwPZumU7zRyx47Qsyesa9XclYtCPtOMx6jq\nhvCOqtaISIsDblVdDOz1PFWtASa0qpUpt2ZN64IZ1J83+/u/j6dNLjr5TjGCjdCrq+Gzz/ILfq44\n5JsAAvVJIF5kt7DyGZlJcKFzeKccS7N3eWqpKGdjwmC2z/W9XexaE8w6doQJE+CRR+Jtk4vOzp02\n+u7VK7/nexJIMvIJZjcC80XkWhG5DpgP/C7eZmVLW0Zmw4bZF99bb8XTJhed1gQzgJ/8BG67zQ9U\n0mLtWluTsV2eK9l6Ekgy8lk1/4/YSvlrsSSQbwXbXJ5amwACtvRRODpzxa21wezss61kzIIF8bXJ\nRSff5I/QgAF+XjQJeR1rqOpSVf2/qnqbqi6Lu1FZ09oEkJAHs3RobTBr1w4uvdRGZ6745Zv8ERo3\nDh5+2BaVdoXTqhIwrvVqa206KZ/VPxoaPx5eeAF27Ii+XS46rQ1mABddBM88Y1NYrri1JvkD7Jxo\np072+3WF48EsZuGoTKT1r+3VC444AubPj75dLjptCWbl5XDuuXD33fG0yUWntSMzEZg8Ga6/3s+L\nFpIHs5i15XxZLp9qLH5tCWYAP/4x3HmnVRd3xWnHDnjoITjllNa97lvfsmKezz8fT7vc3mINZsGK\n+/NEZKmILBGRnwbby0VkloisEJGZ4UohWdTW82UhD2bFr63B7LjjYNAgmDYt+ja5aEydCoceCl//\neute1749XHmljc5cYcQ9MtsJXK6qRwEnAT8WkeFY/bJZqno4MIcM1zPb15HZKadYev6GDS0/1xWe\natuDGdjo7Pbbo22Ti8aKFXDzzfb7actpgu99D5Ytg1dfjb5tbm+xBjNVrVTVN4Lbm4G3gIMokZpm\nqjB9Ohx9dNvfY7/94ItfhLlzo2uXi87dd9vIu2vXtr3+vPNgyRK/nrDYqNpi3//yL3DwwW17j7Iy\n+NnP4IbML9hXHAp2zkxEBgPHYes0lkRNswcftHnzCy9s+bnN8anG4vTww3DNNfDcc21/j7IyuPhi\nuOOOyJrlInD//VaH8LLL9u19Lr7Yzpu9/XY07XJNEy1Auo2IdAWeB65V1SdFZIOqHpDzeI2qljd4\njV599dWf309bPbN162xE9swzcMIJ+/ZeixfDN78J770XTdvcvps2zb6oZs/et5E3wCefwBe+AB9+\n2LZLOFy0qqvtdzpjRn5lX1oyZYpVwbjnnn1/r1LTsJ7ZlClTUNVGJ31jD2Yi0hF4BpihqjcH25YD\nY3Nqms1T1SMbvE4LEWjjctFF0KOHzbnvK1WbynrpJTsZ7ZI1dy5ccAE8+6yV64nCeefZxbaXXhrN\n+7m2mzQJ+vaFG2+M5v1qamx5ukWLYODAaN6zVIlIk8Es7mxGAe4BloWBLJDpmmZz58KcOXDttdG8\nn4hdiDl9ejTv59rub3+D73wHHn00ukAG9es1eh2sZM2aBX/9q42molJebge3UQVH17hYR2Yicirw\nArAICD9oMrAAeBQ4GFgJfFtVP23w2lSOzLZtsymjG29sfTpvc156yaYab7gBfvCD6N7X5e/NN+HM\nM+G+++ArX4n2vVXhrLMs2eCuu9qWPef2zdattu/eemv0v99Vq2DkSMuQ7N072vcuJc2NzApyzqwt\n0hrMfvUry0x77LHo33v5cgtoY8fCLbdYpqMrjBUr6v/fzz8/ns/YvNmSfU46CW66yQNaIW3fbudA\nt2+PrzzPD39odezuuccSf1zrJTbNWGqWLbMVHW69NZ73P/JIW2m9qgrOOMMuyHbx+/BDCzLXXRdf\nIANL758+3aapo5zmcs1bvdoOVDZtslFxXH73O6uLNm6c1zuLgweziNTV2ZHXb36zbyt+tKR7d3j8\ncZsGGT0aXnwxvs9yti7fhAl2vdD3vx//5x1wAMycaWn/N90U/+eVuhdftP3oq1+FP//Z9q+49OgB\nTz5pU9WjR/uaq5FT1aL8saYVv23bVOfNU/3hD1VPPll19+7Cffb06ap9+qj+9reqW7cW7nNLxfr1\nqkcfrXrttYX/7I8+Uj3kENX/+q/Cf3YpqKtTvf121b59bT8qtGeesX33P//T2uLyE8SFRmOGnzNr\npbo6SwSYPdt+5s+H4cPt6P2SS2ytvUJ6/3244gp4+WX49a8ta6pjx8K2IYtqa60Ez9ix8O//nsz5\nq3fesc+/8Ua7FMBF47PP7BKIV1+FJ56AoUOTacc779iCxGPG2JJZnTol04408QSQCNTVWVLHr35l\n9ydOtAA2diz0LIJlkhcssLITH39slwScf37+Zd7dnj77DL78ZTtIueOOZBMxFi+2v7WBA+3vbcIE\nW96sc+fk2pQ2qnbQN2eOHYDOnWtTfXfdBfvvn2zbNm+26es337Rzsuee6/ttczyY7QNVO4cxebL9\nkV1/vX2hFGum2ezZ1tbdu+HnP7cvwj59km5V8Vm3zs49Nlb4dNo0+z974IHi+GLZscOubwtnAxYv\ntqP58ePtb3HUKFul3dVTtWvG/vQn+z/btq3+YGD8+OK6eDls6+TJ9r3y29/aflus3zFJ8mDWRvPn\n2x9YZWX9UVMa/sBU7WT2H/5glaoPPbT+i++005I/Gk1SbS38/veWcXrmmVYAtaG+feGqq4p3unbT\nJlvvb/ZsG22E2Xjhl/Vhh6Xj7zQuL71k++26dfCjH1lgGD68+P9P6ursAOtf/xUOOsgOnMeMSbpV\nxSWxYCYi9wJfBdaq6shgWznwCHAITVwwHTyv4MFs1y5YuNC+IP7yF1i5Eq6+2hYK7tChoE2JzM6d\n8Mor9Uf1r70Gxx9f/8U3enR6+9Ya27fbZRPhyHrKlOTOlURt9WqbOgt/x+3a1R+8jB8P/fsn3cLC\nWLTIVrlftMh+v5MmpXPEumuXHYhOmQJdujS+f/bta7/b8eNLZx+GZIPZacBm4I85wWwqsE5Vp4rI\nlcABqrpXPbMwmC1ebKnujR1B50vVVq2urt77sbo6m7aZPRsqKmwFhvCLfty47J2U3bzZlusJj+pX\nroQvfcn6e8wxjU+r9e5t17gV8sh240Y7/zd8eH5fSLW19iVWV7f3Y2+/bSPrkSPh3/7NVnnIqvBv\nPTw/VFFhU2rjx9u1ieXle7+mXTv7v4kzLT1fqrbowPr1+b9m50649976KfZLLsnGggLbtjW9uPiH\nH9rvONyHTz+9/nsr33113Tr7Tkxi316ypPF9tWtX2z+b2ucTnWYMSr9Mywlmy4HTVbVKRPoDFdpg\nkeHgeaqqXH65/aEOG7bnCfAuXZr/3E8+qd+h58yxK+6byjQ87DCbihg3ru1FFtNq7dr6tSSXL2/8\nOatWWVJE7tF+1Fmb27fD//xP/e9r6VL7XWzYYF/C4WcPG2Y73o4dlsEZPv/NN2HEiMZXVigvh1/8\nwv5uSs2uXTYanz3bppy3bNn7OTt22AX/I0fW/35POqlwAeGjj/bcV7t0sWm21pg4ES6/vDSrDoT7\ncDgy37lzz301/L/cssWuqwuf9/77tm+E+3b4mrbWb2vK9u17nvNdsgSOOqrxafx16+zn7LPtWtqz\nztrzAKzYgtnn5V+ChYhrNKccTM7rPp9m3LHDsvXC/4w33rChdWMXJ+/ebY+vW2fBKQyAvtr8vmmY\nDVZeDscdF800ztq19sc+YkT97+vkk21UvHr1nl907dtbQFu40A5CPMMvGp99Vn8wMXu2je5OPjne\ndQR377ZAu2HDnl++Q4bE95lZF2Zuhr/HefMsmalPH/u/HjVqz1MMHTvCBx/suW/36GGnIqLYt6ur\nbd8+8sj6zz3llOZnvD780MrvTJ9uMwsjR8LUqbaPF20wC+7vVcss2N5kPbPaWjvCqKlp/DNHjGh6\nysztu7o6m9JbujSa9+ve3RJTWrrEIZxCe/dd2yEamzJz0aipsX2stjbezzn6aPuy8n01HnV1dnBf\nXW3BoKWK6OFplyVLovn87t3h1FNtZZu2mDmzggceqKBXL/t+SLqe2WD2nmZstpZZ8LzEsxmdc84V\nj2JbaDjTtcycc84VXtzZjA8DpwO9gSrg18BTtFDLLHitj8ycc859zi+ads45l3rFNs3onHPORcqD\nmXPOudTzYOaccy71PJg555xLPQ9mzjnnUs+DmXPOudTzYOaccy71EgtmInK2iCwXkXeCUjAlq6Ki\nIukmFFyp9dn7m22l1l8ovj4nEsxEpD1wG3A2MAL4rogMT6ItxaDY/igKodT67P3NtlLrLxRfn5Ma\nmZ0IvKuqK1V1J/DfwDkJtcU551zKJRXMDgI+zrn/SbDNOeeca7VE1mYUkXOBs1X14uD+94AxqnpZ\nznN8YUbnnHN7aGptxg6FbkhgFTAo5/4gbHT2uaYa7JxzzjWU1DTjQuAwERksImXAd7A6Z84551yr\nJTIyU9VdIvIT4C9Ae+AeVX0ribY455xLv6KtZ+acc87ly1cAiYiI3CsiVSKyOGdbuYjMEpEVIjJT\nRHrmPDY5uGB8uYicmbP9eBFZHDx2S872/UTkkWD730TkkML1bm8iMkhE5onIUhFZIiI/DbZnuc+d\nRORlEXlDRJaJyPXB9sz2OWhTexF5XUSmBfcz218RWSkii4L+Lgi2Zbm/PUXkMRF5K/ibHpPa/qqq\n/0TwA5wGHAcsztk2FfhFcPtK4Ibg9gjgDaAjMBh4l/pR8gLgxOD2dCzrE+BS4I7g9neA/064v/2B\nY4PbXYG3geFZ7nPQji7Bvx2AvwGnlkCffwY8BDxdAn/XHwDlDbZlub/3A9/P+Zvukdb+JrqTZO0n\n+AXnBrPlQL/gdn9geXB7MnBlzvOeA04CBgBv5Wy/ALgz5zljcv7oqpPub4O+PwlMKJU+A12AV4Cj\nstxnYCAwGzgDmJb1v2ssmPVqsC2T/cUC1/uNbE9lf32aMV79VLUquF0F9AtuH8ielyKEF4033L6K\n+ovJP7/QXFV3ARtFpDymdreKiAzGRqUvk/E+i0g7EXkD69s8VV1Ktvv8e+AKoC5nW5b7q8BsEVko\nIhcH27La3yFAtYjcJyKvichdIrI/Ke2vB7MCUTs0yVy2jYh0BR4H/lFVa3Mfy2KfVbVOVY/FRixf\nEpEzGjyemT6LyNeAtar6OtDodZ9Z6m/gi6p6HPBl4MciclrugxnrbwdgFDYNOArYAlyV+4Q09deD\nWbyqRKQ/gIgMANYG2xteND4QO7JZFdxuuD18zcHBe3UAeqhqTXxNb5mIdMQC2QOq+mSwOdN9Dqnq\nRuBZ4Hiy2+dTgG+IyAfAw8A4EXmA7PYXVV0T/FsNPIGtI5vV/n4CfKKqrwT3H8OCW2Ua++vBLF5P\nAxcGty/EziuF2y8QkTIRGQIcBixQ1UpgU5BRJMAk4KlG3us8YE4hOtCUoH33AMtU9each7Lc595h\nZpeIdAYmAq+T0T6r6i9VdZCqDsHOg8xV1UlktL8i0kVEugW39wfOBBaT0f4G7fxYRA4PNk0AlgLT\nSGN/kzr5mLUf7Mh1NbADmyO+CCjHTp6vAGYCPXOe/0ssG2g5cFbO9uOxHehd4Nac7fsBjwLvYFl0\ngxPu76nYeZQ3sC/017GSPlnu80jgtaDPi4Argu2Z7XNOu06nPpsxk/3FziG9EfwsASZnub9Be47B\nEpneBP6MJYWksr9+0bRzzrnU82lG55xzqefBzDnnXOp5MHPOOZd6Hsycc86lngcz55xzqefBzDnn\nXOp5MHOuSInIj0RkUnD7H4LVGMLH7hKR4cm1zrni4teZOZcCIjIP+Lmqvpp0W5wrRj4ycy4GIjI4\nKGD4YFD08E8i0llExgcrlC8SkXtEpCx4/g1ihU7fFJGpwbZrROSfReRc4ATgoeC1nUSkQkSOD573\n3eD9FovIDTlt2Cwi14kVE50vIn2T+L9wrhA8mDkXn8OB21V1BLAJ+GfgPuDbqvoFbNXy/xOUxPim\nqh6lqscA1wWvV2zh8seBhcDfqeooVd0WPiYiBwI3YPXGjgVGi8g5weu7APPVVvl/AQhLmjiXOR7M\nnIvPx6o6P7j9IDAOK4b4brDtfuBLwEZgWzBS+xbwWRPv17AMiwCjgQpVXa+qu7GK0F8KHt+hqs8G\nt1/Fisc6l0kezJyLT+4JaQE+Zc+AJABBEDoRK8HxNaw6b0vv19Q2ydm2M2d7HTYSdC6TPJg5F5+D\nReSk4PbfYVOFg0VkaLBtElARlBvpqaozgJ9hK5mDBaYw+NUC3Ru8vwILgNNFpJeItMdKtTwfS2+c\nK2J+pOZcfN7GqhXfi9WJugkrg/GnoFDhAuBOoDfwpIh0woLX5cHrc6v8/gG4U0S2YkUz7QmqlSJy\nFTAveO0zqjot5/U08l7OZY6n5jsXAxEZDExT1ZEJN8W5kuDTjM7Fx48UnSsQH5k555xLPR+ZOeec\nSz0PZs4551LPg5lzzrnU82DmnHMu9TyYOeecS73/D5KyJSxwUuYqAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant counts\n", "anhima.loc.plot_windowed_variant_counts(pos, window_size=5000);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAbkAAACfCAYAAABk1qe3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xmc1vPawPHPVULFkYRCKT2yPfb1WE5ZonMecojCkexL\nyK5kSyeJk/TQSY8ocRI5SFmLmiNKpYVoopyTZKmolEqa5nr+uH6ju2lmmrmX3+93/+Z6v173yz2/\nuZfv18x03d/tukRVcc4555KoRtQNcM4553LFg5xzzrnE8iDnnHMusTzIOeecSywPcs455xLLg5xz\nzrnEylmQE5EhIrJYRGaX8b1bRKRYROqnXLtDROaJyFwROTVX7XLOOVd95HIkNxRoU/qiiDQGWgNf\npVzbH+gA7B88Z6CI+CjTOedcRnIWSFR1IrC8jG/1A24vde1MYISqrlfVBcB84Khctc0551z1EOpo\nSUTOBBap6ielvrUbsCjl60XA7qE1zDnnXCJtFdYbiUgdoDs2Vfnb5Qqeslm+MRHxHGTOOec2o6pl\nxpMwR3LNgabAxyLyH2APYLqI7Ap8AzROeewewbXNqGq1ud17772Rt8H76332/np/497fioQW5FR1\ntqruqqrNVLUZNiV5mKouBkYD54nI1iLSDNgbmBpW25xzziVTLo8QjAAmAS1E5GsRuaTUQ34Lv6o6\nBxgJzAHeBDrrlsKzc845twU5W5NT1fO38P29Sn3dG+idq/bko1atWkXdhFBVt/5C9euz9zfZ4thf\nyacBk4j4AM8559wmRASNwcYT55zLO7/8EnULXCY8yDnnXDlWroTGjWHEiKhb4tLlQc4558oxcCDs\nuy906QLz5kXdGpcOX5NzzrkyrF4Ne+0FEybY7amnYPJk2GabqFvmSvM1OedyZPhwOPtsWL8+6pa4\nbBs8GE44AfbfHzp3hmbN4Lbbom6VqyofyTmXpnXrYO+9oWFD+4dw6FCQihLVubzxyy/QvDm89hoc\neqhdW7HC7j/yCPz5z9G2z20qkpFcWfXkRORvIlIoIh+LyMsiskPK97yenMsrgwfDgQfaVNacOdCj\nR9Qtctny9NMW0EoCHEC9evD883DllfDVV+U+1cVMzkZyInIC8DPwjKoeGFxrDbyrqsUi0gdAVbsF\n9eSeA47Eqg+8A7RQ1eJSr+kjORcLa9fCf/0XjB4Nhx8OixfDscdC9+5w2WVRt85lYv16G6GPGAG/\n//3m3+/bF156Cd57D2rVCr99bnORjOS0jHpyqjouJXBNwRIxQ4Lrya1cGXULXC4MGgRHHmkBDmDX\nXeGNN+DOO+Gtt6Jtm8vM8OH2AaasAAdw881Qvz7cfXe47XLpiXLjyaXAG8H9RNaTW7oUmjSBa66B\nX3+NujUuW1avhgcfhJ49N72+zz72Cb9jR5g5M5q2ucxs2AC9e8Ndd5X/mBo1YNgwC4b+gSb+Qqsn\nl0pE7gR+VdXnKnhYmfOSPVIWPlq1ahXLXGklevSAc86BRYvg1FPhxRdh552jbpXL1IAB8Ic/wEEH\nbf69446zUd4ZZ8AHH8Cee4bfPpe+F1+EXXaBli0rflyDBvCPf8B558H06bDbbuG0z5mCggIKCgoq\n9dic7q4UkabAmJI1ueDaxcAVwMmq+ktwrRuAqvYJvn4LuFdVp5R6vbxZk5s717Yfz51rC9Z3321z\n/K++WvY/ji4/rFxpU1kFBbajsjz9+8MTT1ig23HH0JrnMlBcDAcfDH/7G7RpU7nn/PWvMH48vPMO\n1KyZ2/a58sXmnJyItAFuA84sCXCBxNWT69rVbjvtZL/8vXvb7eST4ZVXom6dS9ejj0Lr1hUHOIAb\nb4TTToOzzrKjBi7+Xn3VDnqfdlrln9O9u01f9uqVu3a5zORyd+UIoCXQAFgM3AvcAWwNLAseNllV\nOweP746t0xUBN6jq22W8Zl6M5AoK4JJLoLAQtt120+999JEdHr7iCpv393NV+WPFChvFTZoELVps\n+fHFxdC+ve3AGz7c/jF08aQKRxxhMy5VPQP33Xdw2GHw3HNw4om5aZ+rWEUjOT8MnmXFxbbr7rbb\nbL6+LN99Z4GucWM7QFy3brhtdOm55x74+mv7mVXW2rVwyilw/PG2WcXF05tvwu23w8cfp/dhZOxY\nuPRSmDHD1vRcuGIzXVkdPPecfXLv0KH8xzRqZAeI69Sxf/wWLgyvfS49P/4If/971beN165tZ+lG\njbLnu/hRtbW1O+9Mf7R96qm2q7ZTJ/ug6+LDg1wWrV1rc/R9+255GnLbbW1E0LEjHHMMvP9+OG10\n6enbF9q1s4S9VbXTTjZSuP9+C3guXgoK4Icf4NxzM3udnj1tY1LfvllplssSn67Moj59YNo0OytV\nFW+9BRddBA884Nky4mjJEthvPzv71qRJ+q8zbRr86U+WD/Hoo7PXPpeZk0+2D5sXX5z5ay1caMsV\no0aVf5jcZZ+vyYVgyRLbcffhh7Y5oao+/xzatrWdXf36wVaRnGB0ZbnlFtshOWBA5q/12mu26ej9\n9y0BsIvWpElwwQVWKy5bKbpGj7b6czNn+vGRsHiQC8G119ofSf/+6b/GihVw/vmWO2/kSEsd5KL1\n3XdwwAHw6afZO/A7aJB9kJk0yQ4Vu+j8z//Ywf2rr87u6954oyVxfvll30EdBg9yOZZ68HunnTJ7\nrQ0boFs3O0s3evSWz2O53OrSxc45PvJIdl/3jjvgX/+Cd9+1zSkufDNm2OzJ/PmbH/XJ1Lp1lv2m\nUye4/vrsvrbbnAe5HDvzTAtyt96avdd85hl7vSFD4PTTs/e6rvK+/toyYBQWWgLmbCoutnWgtWst\nlZRnywhfu3aWnu2GG3Lz+l9+aetyb71l5+hc7sSpnlx9ERknIl+IyFgRqZfyvbysJ1dQALNnZ//T\n2kUX2Uju6qvtfFUMY3vi3X+/rZ9lO8CBbVUfMgSWL7es9v7zDddnn1nKtSuuyN17NG8Ojz1mx4m8\nGkl0wq4n9xDwg6o+JCJdgR3zuZ5cycHv22+v+FxcJr75xjIwtGgBTz7pU1thWbDAyuh8/nlu181W\nrLBprcsus2DnwvGXv1gO2a5dc/9eV10Fq1ZZ1htfn8uN2NSTA9oCw4L7w4CSBDp5WU+u5OB3+/a5\ne4/dd7fijGBTK998k7v3chv99a9WIinXG0Pq1bMzdP362bSly7158yxDyTXXhPN+/fvbbM+QIeG8\nn9tU2IfBd1XVxcH9xUDJRFDe1ZMrOfj98MO5/3RWu7aV9TjnHDtfNWXKlp/j0jd/viXrveWWcN6v\nSRM7WtC5sycFCEOfPnDddfC734XzfrVr227pbt1sl64LV2SnsVRVRaSiucdY15Pr3x+OOsqmmsIg\nYlMrBxxgW5779rV1O5d9PXvarsowzzgdcoh9kGnXznZd7rtveO9dnXz1lR3Unjcv3Pfdbz8r4dOh\ngyUFqFMn3PdPmtjWkxORuUArVf1eRBoBE1R133yrJ5fpwe9MzZljW5/POss+lfrOvOwpOQ4yfz7s\nsEP47z90qAXZyZOhYcPw3z/prr3WRnAPPBD+e6vakYKtt7b1dZc9cUrQPBroFNzvBIxKuZ439eTu\nuw8uvDCaAAcWYKdOtYwKp59umxdcdvToYRtAoghwYCWaOnWyn+vq1dG0Iam+/dYKF990UzTvLwID\nB8LEibYJxYUjzHpy9wCvAiOBJsACoL2qrggenxf15AoLbQNINg5+Z6qoyNaN3n7bjhtUpsaZK9/s\n2VYQdf582G676NqharstlyyxqTVP8ZYdt9xiO6KzfbC/qmbNst+zDz7wv9ls8cPgWdS2LbRsGd6m\nhMp48knbBPPss1Wrauw21a4dHHtsPH6269fbaK5ZM3j8cd96nqmlS2GffeyDzO4x2NL2+OPwxBM2\nLZ3tbCvVkQe5LJkwwT5hFxbCNttE1owyvf++HWW49VabjvF/FKtm5kzLYzh/fnw2BaxcabMGHTpY\nGjCXvjvvhGXLLLjEgaqV9mnUyA6Mu8x4kMuCMA5+Z2rhQksxdvDBlgTYPyFW3hln2BRSly5Rt2RT\n335rqaF697YDzK7qli+39fPp06Fp06hbs9GKFZbuq29fOPvsqFuT3+K08SRvDR+e+4PfmWrSxEZ0\na9bAiSfamo7bsilTbJ3kyiujbsnmdtsN3njDRufjx0fdmvw0YIAtM8QpwIElAnj+eUvdt2BB1K1J\nLh/JVcLatTafP2JEeOfiMqEKd91lB5rHj4dddom6RfHWpo2lTst2uZVsmjDBZhDGj4f//u+oW5M/\nVq2yHJLvvx/fTR4PP2zZbiZOzF5Nu+rGR3IZCvvgd6ZEoFcvmwI5+WRbdHdl++AD2yl76aVRt6Ri\nJ55ov4d/+pOndquKQYPsbyCuAQ5slN6gga0buuzzkdwWRH3wOxOqcPfddrzg3Xdh552jblH8nHSS\nrXVddlnULamcPn1sRmHixPDSUuWrtWthr70sT+WBB0bdmor98IOtz/3f/8Ef/xh1a/KPbzzJQOfO\ntpMy6rM16SqZunztNQt0Xol6owkTrNRKYWH+TBOpWtaO+fPh9dfzp91ReOwxm9595ZWoW1I5Eyfa\nmv9HH8XjmEM+iV2QE5E7gAuBYmA2cAlQF3gB2JNSB8VTnhdqkIvTwe9MqNo5ujfftECXz33JFlX7\n2V55pRUvzSdFRZbSbfvt4emnLU2U29S6dTbz8sorcMQRUbem8nr1gnHj7O/UkwBUXkZrciLybmWu\nVaExTYErgMOCnJY1gfOAbsA4VW0BvBt8HamuXS1zeL4HBRHbgt6mja1P/Phj1C2K3rhxtlZ5wQVR\nt6TqttrKduWtXm2H/5cti7pF8fPMM7ZBJ58CHNh5yFq1rNSTy45yg5yI1BaRnYCdg4reJbemZFYG\nZyWwHqgjIlsBdYBvKb/WXCQmTLCyGNddF2UrskfEktKedhqcckr1/oexZK2yR4/8TW5dty68/LIV\ndj3mmPCz6sdZUZH9rt91V9QtqbqaNa0axeDBfmQkWyoayV0FfATsA0xPuY0GBqT7hqq6DHgYWIgF\ntxWqOo7ya82FrrjYUjv16RO/zCaZELE+nXJK9Q50b7xho6A4n3msjJo17SDxbbfB8cdDJSuPJN5z\nz8Gee+bPbujSGjaEYcNsGv2rr6JuTf4rd9ZXVfsD/UWki6o+mq03FJHmwI1AU+An4EURubDUe5db\nay6MenLDh9s6x7nnZv2lIycCDz1k/zC2bg3vvBNu3bSoqcI991gliRoJOUBzxRW2i7B9e/sQE/fj\nELm0YYNNzf/971G3JDOtW9vU5VFH2aiubduoWxQvWa8nJyLHYkHpt6Coqs+k0zgR6QC0VtXLg687\nAscAJwEnlq41V+q5Od94km8Hv9OlaqPV996z9anqEuheecXqtU2fnpwgV2LuXEvq3K6dTdclrX+V\nMXKk7YSeNCkZ+VsnT4bzz7dkBQ8+mKyZpWzKdOPJP4C+wPHAkSm3dM0FjgnW/AQ4BZgDjKHsWnOh\neuSR/Dr4nS4Ry7Rwwglw6qnVoyZdcbGN4nr2TGYA2HdfO885eTKcc071q0dXXGy7E++6KxkBDixv\n6cyZNm153HHw5ZdRtyj/bHEkJyKFwP7ZHEKJyO1YICsGZgCXA9tTTq25lOfldCS3eDEccEB+HvxO\nl6plXJg0yQ7N1qsXdYtyZ+RIW8OaMiU5/wiWZd06Oxrx6acwZozlv6wORo+Ge++FGTOS9/NVtRyc\nPXvaf+OaJD4qGZ2TE5EXsSKm3+aicVWR6yCX7we/06UKN95owX3s2OiqYufShg2W9aJfPztKkXSq\ntj73+OOWw/TQQ6NuUW6pwtFH27Gfdu2ibk3uzJhhAe6kkyzNW+3aUbcoHjLNXbkzMEdExorImOA2\nOrtNjF5hoSVJzcdtx5kS2Zif87TT4Kefom5R9o0YYeuO1aWorIhtXOjXz6ajX3016hbl1rhxNj17\n1llRtyS3DjvM1pNXrbK/18LCqFsUf5UZybUq67qqFuSgPRXK5UjujDOgVat4VIWOiipcf739Eb39\ndnJyIxYVwX77WV7Ak06KujXhmzbNNi7cdJP9fidtKg8se83VV+fn4f50qMKQIZas4m9/g06dkvlz\nrazYpfVKV66C3PjxcPnl8az4HbaS3IizZsFbbyUj0A0dahkwJkyIuiXRWbjQPsgddRQMHJisnJfv\nvWcJtgsLq18qrE8/tenLww+3n+t220XdomhkurvyZxFZFdzWiUixiKzMfjOjUVwMt96avIPf6RKx\nhe2DD7Zs6KtWRd2izPz6qy3W9+wZdUuiVVJQ9/vvbU1y+fKoW5Q9vXrZ1Gx1C3BgqcumTbNzvYcf\nbh9O3aa2GORUdTtV3V5VtwdqA2cDA3PespAMH27BLYkHv9NVo4Ydpj3wwPwPdEOHwt5721GJ6m77\n7WHUKPsAc8wxVskg302ZAp9/DhdeuOXHJlWdOvDkk7aztHVr22yURxN0OZfWdKWIzFLVQ3LQni29\nb1anK9essbNFzz8Pxx6btZdNjOJiuOYamDPHUmFtv33ULaqaX36xAPfPf9rOO7fRoEGWu3PkSFvP\nyldt29rItHPnqFsSD/PmWeab5s0t8CX5SFCqTKcr26XczhWRPsDarLcyAv372z9+HuDKVqOGfSrc\nd1+rSP3zz1G3qGqefNJGLR7gNnf11fDss3ZofNiwLT8+jmbNstpr1TmNWWl7723JABo1sp2YU6dG\n3aLoVWZ35dNAyYOKsIPag1V1SdpvKlIPeBI4IHjtS4B5hFhPrjoe/E5XcbEdLv7iCxvR5cPi9tq1\n9nMdPdrWKlzZCgstFViHDra2lU+ZYNq3t2nXm2+OuiXx9PLL9mGma1fbWZtPP9uqit3uShEZBvxL\nVYcE5XbqAncCP6jqQyLSFdhRVbuVel7Wgtw118C221a/g9/pKi62RMDz51ugq1s36hZVrF8/q7Sc\nL1Who7R0qZ0va9jQdqHWqRN1i7assNCO/Pz73/H/XYzSggVw3nnQoIEV2G3QIOoW5UamGU8aA49i\nuSsB3sMyoCxKszE7ADNVda9S1+cCLVV1sYg0BApylaB5zhxo2dIWrOvXz/jlqo3iYtuqvWABvPZa\nfP9xWb3a1iTGjoWDDoq6Nflh3To7RjN3ro1+GzWKukUVu+gim0bv3j3qlsTf+vVw552WEGH48Pxe\ngy1PphlPhmI15HYLbmOCa+lqBiwVkaEiMkNEBotIXUKsJ9e1q2059gBXNTVq2DrXnnvamas1a6Ju\nUdkGDLA/ZA9wlbfNNjaKO/NMW8OM81b0L7+02YRrr426JfmhVi0rr/XEEzbF26uXpbmrLiozkvtY\nVQ/e0rVKv6HIEcBk4FhVnSYi/YFVwHWqumPK45apav1Sz9V77733t6/TqSfnB78zt2GDLfYvWmQJ\ngOM0vbVypa3FFRTA/vtH3Zr8NHKkBZAhQ+zDTNxceaVNrVb3s4/p+OYbywpTq5ZVIG/YMOoWpad0\nPbn77ruv3JEcqlrhDRgPdARqYvXkLgTe3dLzKni9hsB/Ur4+HngdKAQaBtcaAXPLeK5mYsMG1UMP\nVX3hhYxexqlqUZFqx46qJ5+sunp11K3ZqGdP1b/8JepW5L8PP1Rt1Ej14YdVi4ujbs1GCxeq1q+v\n+sMPUbckf61fr3rPPfbzHTs26tZkRxAbyow5lRnJ7QkMwAqbAkwCrlfVhVUMvqmv+R5wuap+ISI9\ngJKxwI+q+qCIdAPqaZY3njzzjG2JT0pBxaht2AAXX2xZNEaPjj4j+ooVNoqbPNm2UrvMfPWVjeR+\n/3ubAo5DKrAuXWzD2EMPRd2S/Dd+PHTsaH/D992X3xljMt14Mgy4UVWXB1/XB/qqatqnU0TkYOwI\nwdbAl9gRgprksJ7cmjVW8fuFF/xcXDZt2GCbAJYutUz3UQa6e+6Br7+2LCcuO1autMrU69ZZlY4o\nK8h//71NQc+Zk7/TbHGzZIkFutWrbWNK48ZRtyg9mQa5zbKb5GPGk969rcLuiy9muVGOoiILdD/+\naGmjogh0P/4ILVrY4eBmzcJ//yQrKrLqBWPH2q7a5s2jacftt1sWm0cfjeb9k6q42EbGjzwCgwdb\nFpl8k2mQ+xg4UVWXBV/Xx864HZj1lm5BukFu7Vr7w5w4Mbo/0KQrKrJPhMuXW+mP9evt03/q7ddf\n07tWmcesWQOXXAKPPRb1/4nkGjjQNns89RTssUdufo4VPWbDBpg9O39HG3E3aZKN2s8+G666Kpol\nnX32Se95mQa5i7CD2iMBAc4F7lfVZ9JrTvoyGcktXx7tVEt1UFRkOQQ/+MB2rpbctt66al+n+5wd\ndvC11lx7+20bUUHufo7lPaZu3ejXfZNu2TLbWTtjRjTv//nn6T0v44wnInIAcBKWgmu8qs5JrymZ\nyWXRVOecc/kpdmm90uVBzjnnXGmZZjxxzjnn8pIHOeecc4nlQc4551xiRRbkRKSmiMwUkTHB1/VF\nZJyIfCEiY4Oac84551zaohzJ3QDMYWNB1m7AOFVtAbwbfO2cc86lLZIgJyJ7AH/CUnuV7IhpCwwL\n7g8D/hxB05xzziVIVCO5R4DbgOKUa6HVk3POOVc9hJ53WkROB5ao6kwRaVXWY1RVRaTMA3E9evT4\n7X469eScc87lt9L15CoS+mFwEemN1acrArYFfge8DBwJtFLV70WkETBBVfct9Vw/DO6cc24TsToM\nrqrdVbWxqjYDzsPShHUERgOdgod1AkaF3TbnnHPJEodzciVDsz5AaxH5AsuT2Se6JjnnnEsCz13p\nnHMur8VqutI555wLiwc555xzieVBzjnnXGJ5kHPOOZdYHuScc84llgc555xzieVBzjnnXGKFHuRE\npLGITBCRz0TkUxHpElz3enLOOeeyKorclQ2Bhqo6S0S2A6ZjZXUuAX5Q1YdEpCuwo6p2K/VcPwzu\nnHNuE7E6DK6q36vqrOD+z0AhsDteT84551yWRbomJyJNgUOBKXg9Oeecc1kWej25EsFU5UvADaq6\nSmTjSNPryTnnnCtPrOvJAYhILeA14E1V7R9cm4vXk3POOVdFsVqTExuyPQXMKQlwAa8n55xzLqui\n2F15PPAe8Akba8ndAUwFRgJNgAVAe1VdUeq5PpJzzjm3iYpGcl5PzjnnXF6L1XSlc845FxYPcs45\n5xLLg5xzzrnE8iDnnHMusTzIOeecSywPcs455xLLg5xzzrnEilWQE5E2IjJXROYF5XaqtcrmZkuK\n6tZfqH599v4mWxz7G5sgJyI1gQFAG2B/4HwR2S/aVkUrjr8wuVTd+gvVr8/e32SLY39jE+SAo4D5\nqrpAVdcDzwNnRtwm55xzeSxOQW534OuUrxcF15xzzrm0xCZ3pYi0A9qo6hXB1xcCR6vq9SmPiUdj\nnXPOxUp5uSsjK5pahm+AxilfN8ZGc78prxPOOedcWeI0XfkRsLeINBWRrYEOWI0555xzLi2xGcmp\napGIXAe8DdQEnlLVwoib5ZxzLo/FZk3OOeecy7Y4TVcmkogMEZHFIjI75Vp9ERknIl+IyFgRqZfy\nvTuCw/BzReTUlOuHi8js4Hv/m3J9GxF5Ibj+oYjsGV7vNicijUVkgoh8JiKfikiX4Hoi+ywi24rI\nFBGZJSJzROSB4Hoi+5vSppoiMlNExgRfJ72/C0Tkk6DPU4Nrie2ziNQTkX+KSGHwe3103vZXVf2W\nwxtwAnAoMDvl2kPA7cH9rkCf4P7+wCygFtAUmM/G0fZU4Kjg/hvYTlSAzsDA4H4H4PmI+9sQOCS4\nvx3wObBfwvtcJ/jvVsCHwPFJ7m/QjpuB4cDopP9OB+34D1C/1LXE9hkYBlya8nu9Q772N9JfnOpy\nC37wqUFuLrBrcL8hMDe4fwfQNeVxbwHHAI2AwpTr5wGDUh5zdMov49Ko+1uq76OAU6pDn4E6wDTg\ngCT3F9gDeAc4ERhTHX6nsSC3U6lriewzFtD+Xcb1vOyvT1dGY1dVXRzcXwzsGtzfjU2PTZQciC99\n/Rs2HpT/7RC9qhYBP4lI/Ry1u0pEpCk2ip1CgvssIjVEZBbWrwmq+hkJ7i/wCHAbUJxyLcn9BVDg\nHRH5SESuCK4ltc/NgKUiMlREZojIYBGpS57214NcxNQ+yiRu94+IbAe8BNygqqtSv5e0Pqtqsaoe\ngo1w/iAiJ5b6fmL6KyKnA0tUdSZQ5rnVJPU3xXGqeijwR+BaETkh9ZsJ6/NWwGHYdOJhwGqgW+oD\n8qm/HuSisVhEGgKISCNgSXC99IH4PbBPQt8E90tfL3lOk+C1tgJ2UNVluWv6lolILSzAPauqo4LL\nie4zgKr+BLwOHE5y+3ss0FZE/gOMAE4SkWdJbn8BUNXvgv8uBV7Bcu0mtc+LgEWqOi34+p9Y0Ps+\nH/vrQS4ao4FOwf1O2LpVyfXzRGRrEWkG7A1MVdXvgZXBDicBOgKvlvFa5wDvhtGB8gTtewqYo6r9\nU76VyD6LSIOSXWYiUhtoDcwkof1V1e6q2lhVm2FrLONVtSMJ7S+AiNQRke2D+3WBU4HZJLTPQTu/\nFpEWwaVTgM+AMeRjf6Na3KwuN+zT7rfAr9gc9CVAfWzh/gtgLFAv5fHdsd1Jc4HTUq4fjv1hzQce\nTbm+DTASmIft7GsacX+Px9ZqZmH/2M/Eyiclss/AgcCMoL+fALcF1xPZ31J9b8nG3ZWJ7S+2RjUr\nuH0K3FEN+nwwtonqY+BlbDNKXvbXD4M755xLLJ+udM45l1ge5JxzziWWBznnnHOJ5UHOOedcYnmQ\nc845l1ge5JxzziWWBznn8oyIXCUiHYP7FwfZJ0q+N1hE9ouudc7Fi5+Tcy6PicgE4FZVnR51W5yL\nIx/JORciEWkaFJb8R1CM8kURqS0iJwcZ3z8RkadEZOvg8X3ECtB+LCIPBdd6iMgtItIOOAIYHjx3\nWxEpEJHDg8edH7zebBHpk9KGn0Wkl1ih18kisksU/y+cC4MHOefC1wL4u6ruD6wEbgGGAu1V9SAs\nC/w1QemRP6vqAap6MNAreL5iieBfAj4CLlDVw1T1l5LvichuQB+s5tshwJEicmbw/DrAZLXKCe8B\nJaVjnEu5gjwBAAABQ0lEQVQcD3LOhe9rVZ0c3P8HcBJWpHJ+cG0Y8AfgJ+CXYGR3FrC2nNcrXfJG\ngCOBAlX9UVU3YFW8/xB8/1dVfT24Px0r6utcInmQcy58qQvhAqxg00AlAEFwOgordXI6Vk15S69X\n3jVJubY+5XoxNnJ0LpE8yDkXviYickxw/wJsyrGpiDQPrnUECoKyLvVU9U3gZiwzPFjAKgmKq4Df\nlXp9BaYCLUVkJxGpiZXF+VdOeuNcjPknOOfC9zlWXXoIVqerH1Zu5MWggORUYBDQABglIttiQe2m\n4PmpVZmfBgaJyBqsoKk9QPV7EekGTAie+5qqjkl5PmW8lnOJ40cInAuRiDQFxqjqgRE3xblqwacr\nnQuff7J0LiQ+knPOOZdYPpJzzjmXWB7knHPOJZYHOeecc4nlQc4551xieZBzzjmXWP8P/6nbnPLQ\nH80AAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant density\n", "anhima.loc.plot_windowed_variant_density(pos, window_size=1000);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAb0AAACgCAYAAACYEbK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuUFNW1+PHvBkQeIqiAiqAQowhGBYyAAoqCCmjUaG6U\nRBO9Xs1PRV1GjUpyr5BrguHG51Uj11c00Wg0akTxBYoaH7zkpYKKCCowM0iEEZHHwP79savT7dAz\n0z1T1dVdvT9rzZru6urqc6ane9c5dc7Zoqo455xz5aBZ3AVwzjnnCsWDnnPOubLhQc8551zZ8KDn\nnHOubHjQc845VzY86DnnnCsbkQY9ERkhIotF5EMRuaqOfW4NHp8vIn2DbT1FZG7GzzoRuSTKsjrn\nnEs+iWqenog0B94HhgMrgFnAaFVdlLHPKGCMqo4SkQHALao6sNZxmgXP76+qn0ZSWOecc2UhypZe\nf2CJqi5T1S3Aw8DJtfY5CbgfQFVnAB1EZPda+wwHPvKA55xzrqmiDHp7AZmB6rNgW0P7dK21zxnA\nQ6GXzjnnXNmJMujl2m8qdT1PRFoC3wMeDatQzjnnyleLCI+9AuiWcb8b1pKrb5+uwbaUkcAcVV2d\n7QVExBcOdc45tx1Vrd2gAqJt6c0G9hOR7kGL7XTgqVr7PAX8BEBEBgJrVbUy4/HRwF/qexFVLZuf\na6+9NvYyeH29vl5nr2+x17c+kbX0VLVGRMYAzwPNgXtUdZGI/Cx4fJKqThGRUSKyBPgKOCf1fBFp\niw1iOS+qMjrnnCsvUXZvoqrPAs/W2jap1v0xdTz3K6BjdKVzzjlXbnxFlhIydOjQuItQUF7f5Cu3\nOnt94xfZ5PRCEBEt5fI755wLn4igMQxkcc4554qKBz3nnHNlw4Oec865suFBzznnInTZZbB5c9yl\ncCk+kMW5GG3ZAtu2wY47xl0SF4XNm+29XboUevSIuzTlI7aBLI3Npxds7yAij4nIIhF5L1ixxblE\n+cMfYOzYuEvholJVZb9XrYq3HC4tsqAX5NO7DRgB9AZGi0ivWvuMAr6tqvsB5wN/yHj4FmCKqvYC\nDgYW4VzCfPKJtQJcMlUGiyp60CseRZlPT0TaA0NU9d7gsRpVXRdhWZ2LRWUlrFjR8H6uNFVU2G8P\nesWjWPPp9QBWi8h9IvK2iNwlIm0iLKtzsfCgl2ypll4q+Ln4FWs+vRZAP+AOVe2HLUZ9dYhlc64o\nVFbaF2JNTdwlcVGoqIDddvOWXjEp1nx6AnymqrOC7Y9RR9AbN27cv24PHTq0KNd6c64ulZXQvLkN\neOjSJe7SuLBVVECfPh70ojZ9+nSmT5+e076RTVkQkRbA+8AwYCUwExitqosy9hkFjFHVUcHozJtV\ndWDw2KvAf6jqByIyDmitqlfVeo3ET1m48EK48UZo1Srukriwbd1q72vv3nD33XDYYXGXyIXt9NNh\n771h6lSYOzfu0pSPWKYsqGoNkMqn9x7wSCqfXkZOvSnA0iCf3iTgwoxDXAw8KCLzsdGbv42qrMVq\n0yYb0r5kSdwlcVH4/HPo0MG+FP26XjJ5S6/4FHM+vflAWZ/7pi5+L1kC3/lOvGVx4aushN13h732\n8qCXVJWVcPDBsGaNXbdtEek3rsuFL0NWxFJnh97SSyYPeslXUWHv7267werVcZfGgQe9orZypf32\noJdMHvSSbeNG+Ppr2GUX2GMP7+IsFh70itiqVdC9uwe9pPKgl2yVldC5M4jAnnt60CsWHvSK2MqV\ncOSRHvSSKjPopVr1LjkqKqyFBx70iokHvSK2ahUcfrj93rQp7tK4sHlLL9kqKz3oFSMPekVs5Uob\nzt6tGyxbFndpXNhSQa9DB0tBs3593CVyYaqosPcXLOj5UmTFwYNeEVu1yj4s3/62d3EmUSroiXhr\nL4kyuzd9IEvxKOZ8estEZIGIzBWRmVGWs1itXGlLU+27rwe9JEoFPfCgl0SZ7693bxaPYs6np8BQ\nVe2rqv2jKmex2rwZ1q2DTp28pZdE27bZvK3One2+B73k8YEsxako8+llPJ517bRyUFlpAa9ZMwt6\nH30Ud4lcmP75T2jXDlq2tPse9JKn9kCWigpI+FLBJaEY8+ml9lFgqojMFpHzIitlkUp1bYK39JIo\ns+sLPOglUeZAltatbXHxtWvjLZMrznx6KYNVtS8wErhIRIaEU6zSkBrEAtCjByxf7jnXkiRb0PO5\nesmS2dID7+IsFsWaTw9VXRn8Xi0iT2Ddpa/VfpGk5tPLbOm1amUfnk8+gW99K95yuXB4Sy/Z1q+3\nk9R27dLbUiM4e/eOr1xJlU8+vSiD3mxgPxHpjuXTOx0YXWufp7D0Qw8H+fTWqmqliLQBmqvqlyLS\nFjgOGJ/tRTKDXpJktvQgPYLTg14y1A56Xbp40EuSVCtPMvqxvKUXndoNnvHjs4YLoHjz6e0BvCYi\n84AZwNOq+kJUZS1GmS098Ot6SZMt6FVWWmJZV/pqd22CB71iUZT59FR1KdAnyrIVu9otPR/BmSyV\nldZ6T2nZ0lZmqar65vvuSlPmIJYUD3rFwVdkKVLe0ku22i098Ot6SeItveLlQa9IZWvpedBLjmxf\nih70kiNbS8+XIisOHvSK0JYtNnk5tVoHWFfY0qW2kocrfXW19HzaQjJkrsaS4i294tBg0BORk0TE\ng2MBpVZjad48va1tW7vm4y2B0qdq1+4yT2rAW3pJ4t2bxSuXYHY6sEREJorIAVEXyG3ftZniXZzJ\nsHZteoWOTB70kiNb92YqhdSGDfGUyZkGg56q/hjoCywF/igib4rI+SLSroGnukaqPYglxUdwJkO2\nL0TwuXpJkq2lJ2LbPK9evHLqtlTVdcBjwCNAF+D7wFwRuSTCspUtb+klW7breeAtvaRQrfvExgez\nxC+Xa3onB8uATQd2AA5T1ZHAwcDPG3huo/PpBY81D/LpTc61QklQX0vPg17p86CXbNXV0KKFXYev\nza/rxS+Xlt6pwE2q+h1VnaiqVQCqugH4j7qeFEI+PYBLsdVcyiohh7f0kq2uoLfrrrBpE3z1VeHL\n5MKTrWszxYNe/HIJepWq+mrmBhH5HYCqTq3neU3KpyciXYFRwN2UWV69ulp6qfU3PSdXaasr6In4\ndb0kqKtrE9J59Vx8cgl6x2bZNiqH5zU1n95NwJVA2c1Mq6ul16GDjfirqip8mVx46gp64HP1ksBb\nesWtzqAnIheIyEKgp4gszPhZBizI4diNzacnInIiUKWqc7M8nnh1tfTAuziToKGg5y290lZfS88H\nssSvvgWnH8IWi74euIp08PlSVdfkcOym5NM7DTgpuObXCthZRB5Q1Z/UfpGk5dOrqYE1a7afuJyS\nCnqDBhW2XC489QU9794sfd7SK7yw8umpqi4TkYuo1WoTkV1V9Z8NHLux+fQqgLHBDyJyFHBFtoAH\nycunV1UFu+1mo7+y8ZZe6Wuopbd8eWHL48JVUQH9+2d/zINeNPLJp1df0PsLcAIwh+xdlT3qK4Sq\n1ohIKp9ec+CeVD694PFJqjpFREYF+fS+As6p63D1vVaS1Ne1CRb0nnmmcOVx4VJtOOi98UZhy+TC\nlW3dzZTOnW1d3Zqauk9sXbTq/LOr6gnB7+6NPXhj8+nVevwV4JXGlqHU1DWIJSU1gtOVpupq2GEH\naNMm++N+Ta/01de92bw5dOxoPTr1ndy66OQyOX2QiOwU3D5LRG4UkX2iL1p5yqWl50uRla76Wnng\nQS8J6hvIAt7FGbdcpizcCWwQkUOwFViWAg9EWqoy1lBLr2NH2LrVukhc6Wko6HXpYl+ankKqNKUy\naNT3HvsIznjlEvRqVHUbcApwu6reBvhi0xFpqKUn4oNZSllDQW/HHaF9e1i9unBlcuH54ovsGTQy\neUsvXrkEvS9FZCxwJvB0sLzYDtEWq3w11NIDD3qlrKGgB97FWcrqG8SS4quyxCvXfHqbgH8PphPs\nBfw+0lKVsYZaeuCDWUpZLkHP5+qVrvoGsaR4Sy9eueTTW6WqN6jqa8H9T1T1/uiLVp68pZds3tJL\ntoYGsYAHvbjlMnrztCD1T7WIfBn8VBeicOVm61a7ltPQh8ZHcJYuD3rJlktLzweyxCuX7s2JwEmq\nurOqtgt+ds7l4I3NpycirURkhojME5H3RGRC7lUqXVVVll5mhwaumHpLr3R50Es2b+kVv1yCXoWq\nLsr3wE3Jp6eqG4GjVbUPlqz2aBEZnG8ZSk0uXZtg+6xfD2vXRl8mFy4PesmW60CWykpPERaXXILe\nbBF5RERGB12dp4nIqTk8r0n59IIktQAtsWXMEj8zLZdBLGDTFgYOhH/8I/oyuXDlGvQ8vVBpyqV7\ns1Urm9bwxReFKZP7plyCXnvga+A44MTg53s5PK+x+fS6grUURWQeUAm8rKrv5fCaJS3Xlh7AscfC\niy9GWx4XrvXr7fdOO9W/n7f0Slcu3ZvgXZxxanDJU1U9u5HHbmw+PQ1edyvQR0TaA8+LyFBVnd7I\nspSEXFt6YEHvzDOjLY8LV6qVJw1kiNxtN9iwAb7+2loErnTk0tKD9GCWAw+MvkzumxoMeiLSE7gD\n2ENVDxSRg7GBLdc18NSm5NP7F1VdJyLPAN8Fptd+kSTl01u1Cg4+OLd9+/a1kZ6ffQZdu0ZbLheO\nXFsBItYSWLHCBi250rBtm30m68qFmclbeuEKK59eyl3AldganAALsbRDDQW9xubTqxSRjtjyZ2tF\npDVwLJA1QVKS8umtXAkjRuS2b7NmMGyYdXGeU1dCJldUcrmel5Lq4vSgVzrWrLEl5BoafQ2+KkvY\n8smnl8s1vTbBIBPAMssCWxp6kqrWYAHteeA94JFUPr2MnHpTgKVBPr1JwIXB0/cEXgqu6c0AJqvq\ntBzKWtLyuaYHfl2v1DQm6LnSkcvIzRRv6cUnl5beahH51/mmiPwAyOntamw+PVVdCPTL5TWSJJ9r\nemBBb+xY61Zplsvpi4uVB71ky7X7GizozZ4dbXlcdrl8VY7BujZ7ishK4DLggkhLVYa2bm04JUlt\n++wDHTrAggXRlcuFJ5+g17UrLFsWaXFcyHIdxAK+Kkuc6gx6InK5iFyOpRR6FvgNNnn8cSCXeXou\nD6tXWwBr2TK/53kXZ+nIJ+idcAI89hhs2hRtmVx48m3pedCLR30tvXbATsChwP8DdgU6AD+jDLse\no7ZqVX5dmyke9EpHPkGvZ0845BB49NFoy+TCk09LzweyxKfOoKeq41R1PDaloJ+qXq6ql2NBcJ9C\nFbBc5DuIJeXoo+HNN2HjxvDL5MKVT9ADuOgiuP326MrjwpXPQJb27WHz5vSCBa5wcrmm15lvjtbc\nEmxzIcp3EEtK+/Zw0EG+JFkpyDfonXii/V/MmRNdmVx48uneFIFjjoGHHoq2TG57uQS9B4CZIjJO\nRMZjUwg8n17IGtvSAzjuOHjhhXDL48K1YQNs2QI755SfxDRvDhdc4K29UpFP9ybA1VfDxIk2iM0V\nTi5JZH8DnAOsxRZ9PltVfxt1wcpNY1t64Nf1SkGuS5DVdu658MQTNvHZFbd8WnoAQ4ZYkHzssejK\n5LaX0+wuVZ2jqjer6i2qOjefF2hCTr1uIvKyiLwrIu+IyCX5vG6paUpLr39/+Phjm/LgilO+XZsp\nnTrBSSfBvfeGXyYXnmeesVZ8p075Pe+aa2DCBE8zVEiRTmluSk497NrhZap6IDAQuKj2c5OkKS29\nHXaAo46CaYlfs6Z0NTboAYwZA3fc4d1gxWr9erjwQpg0ybqk8zFqlAW8556Lpmxue1Gv49HonHqq\nWqGq84Lt64FFQCPDQvFrSksPvIuz2DUl6B12mLUgnn224X1d4f3Xf9lJ57Bh+T9XxK7tTZgQfrlc\ndlEHvSbl1EsJFq3uiw2iSZx582zKQRhBz7tJis+XX8I990CfPo0/xpgxcNtt4ZXJhWPOHHjwQbjh\nhsYf49/+zZace/318Mrl6hZ10GtSTj0AEdkJeAy4NGjxJcrWrXDeefC73+W/Gkum/fe3s8b33w+v\nbK7pNm6Ek0+2lFEXXdT44/zwh/D22/Dhh+GVzTVNTQ2cf76NwMz3Wl6mFi3gyiu9tVcouSw43RRN\nyqknIjsAfwP+rKpPZnuBUs+nd/vt0LZt09MDiaRbewccEE7ZXNNs2WJn8bvvDnfemf/IzUytWtlI\nzjvugJtuCq+MrvH+939t6cCf/KTpxzr7bPj1r20d3Vxzarq0fPLpiUbYHyYiLYD3gWFYTr2ZwGhV\nXZSxzyhgjKqOCnLq3ayqA0VEsGt9a1T1sjqOr1GWP2qffmrJYF9/3ZadaqqHH7bJrk891fRjuabZ\nutUy269fD48/nluOtYYsXw79+sEnn9iJkovP8uVw6KG2GtJ++4VzzIkT7VKHT1hvOhFBVbOeZkYa\n9IIXHwncDDQH7lHVCRn59CYF+6RGeH4FnKOqb4vIYOBVYAHp7s5rVPW5jGOXbNBThVNOsS+xa68N\n55irV1vS0VWroE2bcI7p8qdq3V4ffWRD2Vu3Du/Yp5wCxx9vk9ZdPFThe9+Dww+HX/4yvONWV8O3\nvgUzZsC++4Z33HIUa9CLUikHvccftw/MvHmw447hHffcc+2Yjz9uqYdcYanC5ZfDG29YV3O7duEe\nf/ZsG+b+97/bl64rvEcfhfHj7RprU67DZ/OrX8Hnn1t3uGu8+oKepx6Nwbp1cMkl8H//F27AA7j7\nbvjRj2DAAHjppXCP7Rr261/bfMkpU8IPeADf/S7cf7+1+ObNC//4rn7z59tnd9Kk8AMewKWXwl//\nCjlennKN4EEvBr/8JYwcacsQhU3EWhoPPWTB74YbfBpDodx0kw1ff+EF2HXX6F5n5EibvjBqlI/W\nLaQHH4Thw+19HjQomtfo1MmC3hln2Ov4Zzd83r1ZYG+9BaeeCu++C7vsEu1rLV8Op51mF9rvvtsH\nP0Tp7rvhuuvg1Vdh770L85r33WfXg197zbuyo7RlC/ziFzZA7IknCjO6ctky+57o1Qvuusuv0efL\nuzeLwIoV8MADNjT5xhujD3hgX4SvvWZdqIcfbteZXPgeecSCz4svFi7ggU1zueIKa314Fu5oVFXZ\nVKDFi2HWrMJNJ+je3UZ1N28ORxwBS5cW5nXLgQe9iKxbZ4MNLr4Yeve2D8vkybbk0OmnF64crVtb\ni+DKK2H0aFu8eOHCwr1+0j39tF2Hee658Iau5+OSS+CnP7X0Up6JIVyzZtkScIMH2/scZZd1Nq1b\n2/Xbc8+1k9bnny/s6yeVd2+GrKoKfvMb+OMfbTDJ8OH206cPNIv5FGPjRhsVNmGCfUmOH29DpF3j\nvPyyncA8/bRluoiLKowda+/t0Uen/+f2269pE+LLzaZNNu9u6lT7+fBD67b+/vfjLpl1m48ebe/r\n+PHWEnR18ykLBbBunQ0auf12m5Q8dmzjFxiOWnW1XSS/9Va7YH7mmXZG2yLq9XlK0GuvZR8luXEj\n/M//2KCDYlkEaNUqGzma+tIWSQfAYcPyS3BaLtautXVRX3jBuv97907/vY44wlbCKRbV1fYdc9tt\n8OMf24C4Yv2OiVvck9NHkJ6cfreq/i7LPrcCI4ENWJLaucH2e4ETgCpVPSjL82IPel9/bYFu4kQb\nTTduXOmcha1eDTffbBOoly2DI49Mf0n26lXerYQ5c+zEZckSGy2Z7W9x6qnWsipGqvDBB+kAOH06\n7LVX+v096qhoplSUig0bLHj8/vc22f+00+zkpUOHuEvWsKoq+O1v4U9/spRGV1wB7dvHXar4qaY/\np7EFvSCf3vvAcGw9zVnUvwzZAOAWVR0YPDYEWA88UExB7/PPrWtr6lS7TjdggI3cO/DAghclNFVV\nNq9v2jQbkLF5s53tps56u3Zt+BhJ8P778J//aYMIfvUru54SxXysQqupscnUU6faezxjBhxySDoI\nDhiQjHo2ZMsWS8j73/8NAwfa714lmqVz+XI7yX7qqeyt+GbN7LJKkj/D69bZtc7Jk61HZsECC3xx\nBr3DgWtVdURw/2oAVb0+Y587gZdV9ZHg/mJgqKpWBPe7A5PrC3qbN1uFDz00/ySOmb76yv5oNTXb\nP1ZdDa+8Yl8aS5akW0XHH1+6H5q6qNposVRX2UsvQceOVt9jjsm+oryIBf1CjErNLOcHH9jZeS7d\nPKr23lVUbP/Y1q02D+vJJ22e48UXJ3uKx4YN8I9/2Hs8bZr9HQcPtve4X7/sn6NddrHuv7ivTYOt\nabpgQX6JdT/+2ILcPvtYSynO67BhWrkSvvhi++1btsDMmenPcKdO6ZOcoUNzax1u2waLFlnas0IO\n5El9tquqspfp7bftWvqsWTbf+cQT7adbkLogzqD3A+B4VT0vuH8mMEBVL87YZzIwQVXfCO5PBa5S\n1TnB/e40EPSWLrUVKj791N7M1JlNz571d9HV1NgfLXX2O3u2BbBs/fitWqW/FPr3D2cB4VKxbZud\nVEybZt1k1dXb71NTY3MPDzgg/fcfNCj8ayIrVnzzulWLFpavrlu3dMv0yCPTXXcVFekv9qlTrS49\nemQ/9pAhNsq1kIG7WKxZY70X06bBO+9k36eiwq6BDRuW/lvX9bcMW+YX+LRp9qXXu3d+KxrtvDNc\ndpmVu9ykPsOpz82bb9pJaioIHn54+m/58cfp/V56yT5Ln39u36ep933QoHDXlAUL3qnP6bRpdnJV\n1/zTnj1t/dPhw7OfnMYZ9E4DRuQQ9K5X1deD+1OBX6jq28H97jQQ9FIqKr7ZRbdtm3VhZPtgrF1r\nXVjdu6ff+CFDkn12H7VNm2zyfeofd+FCO0EIYwBF6kNbVWWtzVRg3XdfO9ufMyf9YZk507ruqqvh\ns8/SIxqHDUvnHXSN88kn3/xiatvWlkaLchDUF1/YZ3XffdPv4+DB/lltio0b0yNVp02zE9bDDrNr\n+xs2fHMAVLdudrnjrbfS+8+fb5/tpiS+Ttm2zY5XUZH+bA8fbu93Yz+rcQa9gcC4jO7Na4BtmYNZ\ngu7N6ar6cHB/MXCUqlYG97tTT9C7NiNFQWY+vVRX1pw52btA2rSxs5XOncOpq9veunXWhbZ2bTjH\n69Urt6kfGzbYSLydd256l7erm6p9WS5YEO1yWW3bWpDr2DG61yh3a9faZ2affawF3VCwqa62z3a2\nbtXG6NnT0qw19rNaO5/e+PHjYwt6jc6nl/F4d3Js6TnnnHOxLUOmqjXAGOB54D3gEVVdJCI/y8ip\nNwVYKiJLgEnAhRkF/wvwBrC/iHwqIk3ML+6cc66c+eR055xzieILTjvnnHN40HPOOVdGPOg555wr\nGx70nHPOlQ0Pes4558qGBz3nnHNlw4Oec865shFp0BORESKyWEQ+FJGr6tjn1uDx+SLSN5/nOuec\nc/mILOgFufRuA0YAvYHRItKr1j6jgG+r6n7A+cAfcn1uOcpcW64ceH2Tr9zq7PWNX5Qtvf7AElVd\npqpbgIeBk2vtcxJwP4CqzgA6iMgeOT637BTjP1CUvL7JV2519vrGL8qgtxfwacb9z4JtuezTJYfn\nOuecc3mJMujluiimZzdzzjlXEJEtON2UXHpAj4aeG2z31aadc85tp64FpyPMd8xsYL8gH95K4HRg\ndK19nsJSDz0cBMm1qlopImtyeG6dlXLOOeeyiSzoqWqNiKRy6TUH7knl0gsen6SqU0RkVJBL7yvg\nnPqeG1VZnXPOlYeSzqfnnHPO5cNXZCkwEblXRCpFZGHGtl1F5EUR+UBEXhCRDhmPXRNM0F8sIsdl\nbD9URBYGj92SsX1HEXkk2P6WiOxTuNptT0S6icjLIvKuiLwjIpcE2xNZZxFpJSIzRGSeiLwnIhOC\n7Ymsb0aZmovIXBGZHNxPen2XiciCoM4zg22JrbOIdBCRx0RkUfB/PaBk66uq/lPAH2AI0BdYmLFt\nIvCL4PZVwPXB7d7APGAHoDuwhHTrfCbQP7g9BRgR3L4QuCO4fTrwcMz13QPoE9zeCXgf6JXwOrcJ\nfrcA3gIGJ7m+QTl+DjwIPJX0/+mgHB8Du9baltg6Y/Op/z3j/7p9qdY31n+ccv0J/hEyg95iYPfg\n9h7A4uD2NcBVGfs9BwwE9gQWZWw/A7gzY58Bwe0WwOq461ur7k8Cw8uhzkAbYBZwYJLrC3QFpgJH\nA5PL4X8aC3q71dqWyDpjAW5plu0lWV/v3iwOu6tqZXC7Etg9uN0Fm5ifkjl5P3P7CtKT9/814V9V\na4B1IrJrROXOi9ho3L7ADBJcZxFpJiLzsHq9rKrvkuD6AjcBVwLbMrYlub5g85CnishsETkv2JbU\nOvcAVovIfSLytojcJSJtKdH6etArMmqnOokbXSQiOwF/Ay5V1S8zH0tanVV1m6r2wVpAR4rI0bUe\nT0x9ReREoEpV51LHQhNJqm+GQaraFxgJXCQiQzIfTFidWwD9sO7HfthI+6szdyil+nrQKw6VYmuO\nIiJ7AlXB9hVAt4z9umJnSiuC27W3p56zd3CsFkB7Vf1ndEVvmIjsgAW8P6nqk8HmRNcZQFXXAc8A\nh5Lc+h4BnCQiHwN/AY4RkT+R3PoCoKqrgt+rgSew9YKTWufPgM9UdVZw/zEsCFaUYn096BWHp4Cf\nBrd/il33Sm0/Q0RaikgPYD9gpqpWANXBCCoBzgL+nuVYPwCmFaICdQnKdw/wnqrenPFQIussIh1T\no9hEpDVwLDCXhNZXVceqajdV7YFdo3lJVc8iofUFEJE2ItIuuN0WOA5YSELrHJTzUxHZP9g0HHgX\nmEwp1jeui6Pl+oOdDa8ENmN92OcAu2IDAT4AXgA6ZOw/Fhv9tBg4PmP7odgHbQlwa8b2HYG/Ah9i\nIwe7x1zfwdi1nnnYl/9cLGVUIusMHAS8HdR3AXBlsD2R9a1V96NIj95MbH2xa1zzgp93gGvKoM6H\nYIOy5gOPY4NbSrK+PjndOedc2fDuTeecc2XDg55zzrmy4UHPOedc2fCg55xzrmx40HPOOVc2POg5\n55wrGx70nCtxIvIzETkruH12sDpG6rG7RKRXfKVzrrj4PD3nEkREXgauUNU5cZfFuWLkLT3nYiQi\n3YNEm38UZEY/AAAB7klEQVQOknM+KiKtRWRYsKL9AhG5R0RaBvtfL5aQd76ITAy2jRORy0XkNOC7\nwIPBc1uJyHQROTTYb3RwvIUicn1GGdaLyHViiW/fFJHOcfwtnCsED3rOxW9/4HZV7Q1UA5cD9wE/\nVNWDsVXuLwhSrZyiqgeq6iHAdcHzFVvo/m/AbOBHqtpPVTemHhORLsD1WM67PsBhInJy8Pw2wJtq\nmSFeBVKpcpxLHA96zsXvU1V9M7j9Z+AYLGnnkmDb/cCRwDpgY9Dy+z7wdR3Hq53iR4DDgOmqukZV\nt2JZzo8MHt+sqs8Et+dgSY6dSyQPes7FL/PCugBr+WbgEoAgWPXHUruciGWbbuh4dW2TjG1bMrZv\nw1qWziWSBz3n4re3iAwMbv8I66LsLiL7BtvOAqYHaWw6qOqzwM+xle/BAlgqSH4J7Fzr+ArMBI4S\nkd1EpDmWBuiVSGrjXBHzMzrn4vc+ln37XixP2Y1YepVHg4SaM4E7gY7AkyLSCgtylwXPz8xa/Ufg\nThHZgCV4tR1UK0TkauDl4LlPq+rkjOeT5VjOJY5PWXAuRiLSHZisqgfFXBTnyoJ3bzoXPz/zdK5A\nvKXnnHOubHhLzznnXNnwoOecc65seNBzzjlXNjzoOeecKxse9JxzzpUND3rOOefKxv8HGB3vfS/R\nj4EAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot variant density\n", "anhima.loc.plot_windowed_variant_density(pos, window_size=5000);" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAcMAAACbCAYAAAAeLPDwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xm81XP+wPHXe0paRMqoKAohBhWSLbfCNI1JdkVjGWQp\nZKvwq2trs0yqQciuhVRTGLJ0yxhKWpS6LciUlMK00N7798f7e3W63eXce8/3fs/yfj4e38c953u+\n33Pen+7tvL+fz/eziKrinHPOZbLfRR2Ac845FzVPhs455zKeJ0PnnHMZz5Ohc865jOfJ0DnnXMbz\nZOiccy7jVYw6gLCIiI8Zcc45twtVlYL2p3XNUFUzYuvTp0/kMXh5vcxeXi9vspe3KGmdDJ1zzrl4\neDJ0zjmX8TwZpoGsrKyoQyhXmVZeyLwye3nTWzKWV4prR01VIqLpWjbnnHMlJyJoJnagcc658rJp\nU9QRuLIINRmKSFsRyRWRxSLSo5BjBgevzxGRpsG++iIyWUS+FJF5InJzzPHZIrJcRGYFW9swy+Cc\nc8XZvh2aNoX+/aOOxJVWaMlQRCoAQ4G2wFFARxFpnO+YdsBhqtoIuA54MnhpK9BdVY8GWgA3iciR\nwWsKPKaqTYPtnbDK4Jxz8Xj9dahcGQYNgo8/jjoaVxphDrpvDixR1aUAIjIKOBdYEHNMe+BFAFWd\nJiI1RKS2qq4EVgb7N4jIAuBAIDc4r8A2X+ecK287dsBDD8HDD8PWrdCpE8ycCbVqRR2ZK4kwm0kP\nBJbFPF8e7CvumHqxB4hIA6ApMC1md7egWXW4iNRIVMDOhSUnB7KyYO3aqCNxiTZhAuy5J/zxj/CX\nv8CFF8LVV4P330stYSbDeP8U8tfyfjtPRPYCxgC3qOqGYPeTQEOgCfA98GgZ43QuVKrQowf8+qt9\nUW7ZEnVELlFU4cEH4d57QYJvsn794PvvYfDgaGNzJRNmM+l3QP2Y5/Wxml9Rx9QL9iEiewBvAK+o\n6vi8A1T1h7zHIvIsMLGwALKzs397nJWVlZRjW1z6e+stS4QzZ8IFF8C118ILL+z88nSp6913YfNm\naN9+575KlWDUKGjRAk49FU44Ibr4Ml1OTg45OTlxHRvaOEMRqQgsBNoAK4DpQEdVXRBzTDugq6q2\nE5EWwCBVbSEigt1L/FFVu+d737qq+n3wuDtwoqp2KuDzfZyhi5wqHH+81RzOPx9++QVatYK2beH+\n+6OOzpWFKpx2GnTrBpdeuvvrr78OPXvaRdA++5R/fG53kYwzVNVtQFfgXWA+MFpVF4hIFxHpEhzz\nNvC1iCwBhgE3BqefClwOtCpgCMUAEflCROYAZwC7JMtUtX693Yh36WX8ePvS7NDBnlerBm++Ca++\nCs8+G21srmymTIHVq+Giiwp+/aKL7D7iddf5/cNU4DPQJIFNm+Coo+DYY+Hll6F69agjcomwYwcc\ndxz07WsdK2ItWgQtW8Lzz8Of/hRNfK5szjwTLr8crryy8GM2bYKTToKbbrKk6KLlM9AkuSFDLBnu\nvz+cfDJ8/XXUEblEeP11qFoVzjln99cOPxzGjoW//tWa0Vxq+eQTWLIELrus6OMqV4bRo+Gee+CL\nL8onNlc6XjOM2Jo10Lgx/Pvf9gX5j39Y77SRI+3ekktN27fDH/5gg7D/+MfCjxs71u45ffwxNGhQ\nbuG5MjrnHNuuvz6+419+2VoIZsywpnIXjaJqhp4MI3bLLfbFOXTozn0ffGBXnL17w403Fn6uS16v\nvAJPPmkXOcX1Gn38cRg2zBLivvuWT3yu9GbOtGbvr76yml+8rrzS/haefz600FwxPBkmqcWLrVl0\nwQL4/e93fe2rr6y7dsuW9mVZqVI0MbqS27bNavvDhkHr1vGdc9tt8PnnMGmSDeB2yeuCC+D00+HW\nW0t23oYNcOKJ0KuXNY+78ufJMEmdfz40b27drwuybp3doF+3zu4/5U+YLjk995w1i02eHP85O3bA\nJZdAhQowYgT8zu/mJ6Uvv7QLnK+/Ll1z5xdfQJs21mJwxBGJj88VzTvQJKGPPrKawC23FH7M3ntb\n1/xTT7Wk6Tfgk9+WLfDAAyUfQ/i731kCXb688IsjF72+faF799Lf9zv2WOsTcPHFsHFjYmNzZeM1\nwwjs2GGzU9xyS/G90fKMHGnHDxsG550Xbnyu9IYNgzfesObO0vjxR7v46dbNuuO75LF4MZxyit3C\n2Hvv0r+Pqg3Sr1ULnngicfG54hVVMwxzOjZXiNdes/8QHTvGf07HjtCokTWtzpu361yILjls2mRX\n/WPGlP49atWCf/3LEmK9enDuuYmLz5VN//52gVKWRAj2//bpp6FZM7v9UdigfVe+vGZYzjZtss4V\nL7wAZ5xR8vO//94SYr169h7eTTt5DBlic1W++WbZ32vGDBuM/+abNmjbRevbby15LV4MNWsm5j1n\nzIB27eDTT+GQQxLznq5ofs8wiQwZYrOSlCYRAtStax0z9trLag/ffpvY+FzpbNxoqxUkar7RE06w\nLvgdOtjgbhetgQNtgvVEJUKw3/Hdd1uTqa9kEj2vGZaj2AH2Ze1JpmoDugcOtKaW005LTIyudB57\nzDpFjRuX2PcdNgweeQT+8x/vTRyVFStsAoXcXJslKpHy5q097DB41BejC11kNUMRaSsiuSKyWER6\nFHLM4OD1OSLSNNhXX0Qmi8iXIjJPRG6OOb6miLwnIotEZFIqLe77wAPWfT4RXapFrFfbCy/YuCef\n9Dk6GzbYRcl99yX+vbt0sXtK7dvbMlCu/D36KFxxReITIewchD9mTGKa113phbmEUwVsCaczsTUK\nP6PoJZxOAh4PlnCqA9RR1dnBAr+fA+eqaq6IDATWqOrAIMHuq6q7dUZPtprhokXWE62gAfaJeO/2\n7eHss62GUtG7RZWr/v1h1iybgzIMqtC5sy3/NGaMjUV05WP1art4nTsXDjwwvM/5+GO7qJ0xw/oD\nuHBEVTNsDixR1aWquhUYBeTvG9ceW7cQVZ0G1BCR2qq6UlVnB/s3AAuAA/OfE/zsEGIZEqZnT7jz\nznCaug4/3G7CL15s6+T99FPiP8MVbN06uwCJWUc64URsIP/atdYakETXeGlv0CBrzQkzEYLd/7/l\nFus1vm1buJ/lChZmMjwQWBbzfDk7E1pRx+xyXSQiDYCmwLRgV21VXRU8XgXUTky44ckbYH/zzcUf\nW1o1algzS9OmNkB//vzwPsvt9PjjViNv3Djcz6lUySb1/vBDS74ufD//DE89BT0KvMGTeD16QJUq\n4TS3u+KF2aAW7/Vr/irrb+cFTaRjgFuCGuKuB6qqiCT1dfKOHXD77TZzRZUq4X5WhQrw8MNwzDGQ\nlWW1iYKWD3KJ8fPPlgw/+aR8Pq9GDXj7bWtur1/fZjFx4Rk61G4/lNdqInmzEDVrZr3NzzyzfD7X\nmTCT4XdA/Zjn9bGaX1HH1Av2ISJ7AG8Ar6jq+JhjVolIHVVdKSJ1gR8KCyA7pu0qKyuLrKyskpei\njEaPLvkA+7L661+t6fTCC20mk7vu8gH6Yfj73+3LslGj8vvMgw6yFoCzzrJhNqefXn6fnUnWr7dh\nUP/+d/l+bu3a8NJLdo941ix77kovJyeHnJycuI4NswNNRawDTRtgBTCdojvQtAAGBR1oBLsf+KOq\nds/3vgOD/QNEpCdQI1k70GzaBEceCS++WPpxhWXx3XfWbfvww623adg100zy44/27zpjBjRsWP6f\nP2mSfWHm5ITfRJuJHn7YlmoaOTKaz+/d21oc3n3XJ21PpEg60KjqNqAr8C4wHxitqgtEpIuIdAmO\neRv4WkSWAMOAvNX7TgUuB1qJyKxgaxu81h84S0QWAa2D50lpyBBo0iSaRAh203/qVHvcsqUlR5cY\nDz9sNe8oEiHYfcoBA2wGk5Uro4khXW3caPdl7747uhh694bNm62nsisfPug+JIkcYF9WqjYObvBg\nm0S6RYto40l1q1bZ73bOHLt3F6X77oOJE62GuNde0caSLoYMsY5KiZ5AoaSWL7dZal5/3ZvDE8XX\nM4zAzTdbEhoyJLIQdvPmm3D11VarueKKqKNJXbffbtNnJcPvVhWuucZqh//8p48xLavNm202mHHj\nLBFF7a234IYb7P5hrVpRR5P6PBmWszAH2JfV/PnW6aNDB2tm8wHcJZM3Nde8eXDAAVFHY7Zuhb/8\nBQ4+2IYCeGep0nvmGRvC8q9/RR3JTnfcAQsXwoQJ/rstK0+G5ez8822lgfIan1RSP/1k3fL32MM6\nCNRImQntotetm/27JdtYv/Xr7b7wRRdFe68rlW3bZp2iXn7ZBsEniy1brJn00ktt0gVXemXqQCMi\n7UXE+zPFaepU64VW1Ar2UatZE955x/7jn3SSr4oQr//+F159NTlXoq9e3ZrUnn4aXnkl6mhS08iR\nVrtOpkQINuHCqFG2Kspnn0UdTfqKJ8ldAiwRkYEicmTYAaWyHTusSaNvX6hcOepoilaxog0Y794d\nWrWyqdxc0fr2heuuC2fC5kQ44ABLiLffbh1AXPy2b7ff7733Rh1JwRo2hCeftKnh1q6NOpr0VGwy\nVNXLsOnQvgZeEJFPROQ6EakeenQpJm+A/aWXRh1J/K6/Hvr0gdatvYZYlG++sV59d9wRdSRFO/po\n+zu89FKbXNrFZ+xYu13QunXUkRTuggtswedrrvH5acMQV/Onqq7FpkUbDRwAnAfMil1aKdNt2gS9\netlyL6k2SPaaa2xcU+vW8NVXUUeTnB54AG68EfbbL+pIipeVZbX+P//Zuue7oqnCgw9arTDZO6g8\n+qi14gwbFnUk6afYjtgici5wJdAIeAk4UVV/EJGq2GD6waFGmCIGD7ZJslu2jDqS0rn2WvtSaNUK\nJk+GQw+NOqLksXix9eRLpabkjh3tHuef/2wTxe+9d9QRJa8337QL2Hbtoo6keJUrw2uv2X3NU06B\nY4+NOqL0UWxvUhF5ERiuqlMLeO1MVX0/rODKojx7k65ZY9Ou/ec/1ikllT31lN2onzwZDjkk6miS\nQ+fO9nv9v/+LOpKSUYWbbrIV2seP94RYEFWbhOKuu6wZMlW88orVZmfM8MkWSqKs07Gtyp8IRWQA\nQLImwvJ2//12JZ7qiRDsHmLPntZk+s03UUcTvQULbH7IZO4dXBgRmxigUSOrSXz7bdQRJZ/337dh\nKeedF3UkJXP55VYzvOmmqCNJH/Ekw7MK2JcCDQrlY9EiGDHC7rmlixtusCvlVq1g6dKoo4lWdjbc\ndlvq1qoqVLDa/tVXw8knw7RpxZ+TSR58EO65J/Xu84Nd6Hz2mS0E4Mqu0D8BEblBROYCR4jI3Jht\nKfBFPG8uIm1FJFdEFotIgUPQRWRw8PocEWkas/85EVkVxBB7fLaILC9gAu9I9OhhiSPZZpopqxtv\nhDvvzOyEOHcuTJkCXbtGHUnZiNgQmmHDbH3L116LOqLkMHWqTV5/ySVRR1I61arZ7/KOO2y6NldG\nqlrgBuwDNABGAQcHjxsAtQo7J9/5FYAlwTl7ALOBxvmOaQe8HTw+Cfg05rXTsSEdc/Od0we4LY7P\n17BNmaJ68MGqGzeG/lGRGTxYtWFD1aVLo46k/J13nuojj0QdRWLNmqVav77q/fer7tgRdTTROuss\n1WefjTqKshszRnW//VSfftp/p8UJ8kKBOaOoxgFV1aXATcB6YF2wqYjUjCPPNgeWqOpSVd0aJNVz\n8x3THlu3EFWdBtQQkTrB84+Anwt578g7QKfSAPuy6NYNbr3Vaoj//W/U0ZSfmTPh00+tyTidNGli\nTaUTJtgi0Js3Rx1RNKZNs45FnTtHHUnZXXCB9RgeOhQ6dYJ166KOKDUVlQzzlrX8vJCtOAcCy2Ke\nLw/2lfSYgnQLmlWHi0gkM2uOGpV6A+xL6+abbWvVCpYtK/74dNC7t40brVo16kgSr25da/7duBHa\ntIHVq6OOqPw99JDd4qhUKepIEuPII+3irUYNaNYMPo/nG9rtotBkqKp/Dn42UNWG+bc43jvecQ35\na3nFnfck0BBoAnwPPBrn5yTMpk02GXIqDrAvrVtvtXtnmZAQp02ztQqvvTbqSMJTtardb2rZ0uan\nnT8/6ojKz+zZNiTh6qujjiSxqlSxKdv69bOZagYP9plqSiKeQfenAnNUdYOIdMbu4z2uqsV11P4O\niF36tD5W8yvqmHrBvkKp6g8xsT0LTCzs2Ozs7N8eZ2VlkZWVVUzI8Un1Afal1b37zoH5OTlQr17U\nEYWjd2/rYZjOzd9gF3J9+9ri01lZNgn5WQX1HU8zffva/K1VqkQdSTguugiOP946Bn34ITz3nE3O\nn4lycnLIycmJ7+DCbibqzo4oc7Ea5HHALKArMCWO8yoCX2EdaCpRfAeaFsR0oAn2NWD3DjR1Yx53\nB0YU8vmJu+saY/Vq1Vq1VBcuDOXtU8Ijj6gedpjq8uVRR5J4H32k2qCB6ubNUUdSvqZMUa1dW/XJ\nJ6OOJFzz56vuv7/q+vVRRxK+zZtVu3dXPegg1Y8/jjqa5EARHWjiSYazdGcvzmuCxzOLOy847k/A\nQqxXaa9gXxegS8wxQ4PX5wDNYvaPBFYAm7H7ilcF+1/ChnbMAcYDtQv57FD+Mbt2tS3TDRyYngkx\nK0t1+PCoo4jG4sWqRxyheuutqtu2RR1NODp3Vn3ooaijKF8TJtgFQL9+qtu3Rx1NtIpKhvFMxzYV\neAe4ChvusBqYrarHFF3njFYY07EtXGgzeeTmpsaEzWEbOBCefdaaTJNl1feymDzZ7hPm5toSV5no\n55/hwgvtnuKIEbZOYrr46iu7P/rVV7DPPlFHU76WLbNZsqpVs8WLk3UZsrCVdTq2S7Da2dWquhLr\n7flIAuNLGT172gB7T4TmrrusE0KrVrBiRdTRlI2qzT3ap0/mJkKAffe1hZ/r1IHTTkuvzlIDBthk\nEpmWCAHq17eL1hNPtN6mvt7l7oqtGaaqRNcMp061cVm5uenfsaKk+vWzKaEmT7Zu+6no3Xetx+y8\neTaFWaZThcces23cOGjePOqIymbZMhtjuWgR1KoVdTTReu89uOIKawXp3Tuz/t7LVDMUkQuC6dLW\nicj6YMuoYZ07dljvs379PBEWpFcvu1Bo3RpWrow6mpLLqxVmZ2fWF0NRROxv/oknbBmoMWOijqhs\nHn4Y/vY3T4RgPYZnzoSPP7Zxpt8V2X8/c8TTTDoQaK+qe6tq9WBL0WmLS2fUKPtySNU5DMvD3XfD\nZZdZk2mqJcS33rIB6BddFHUkyefcc2HSJBtW07dvao5bW7nSljy67baoI0kedepYa8hZZ9kwjLff\njjqi6MXTgeZjVT21nOJJmEQ1k27aZLM7vPRS5o0rLI0HHoCRI63JtHbtqKMpnqp9Gdx7L5x/ftTR\nJK8VK6B9ezj6aHj6adhzz6gjit9dd9n/48G+DHmBPvrIpnHr2NFm5tljj6gjCk9RzaTxJMPHgTrY\nMIYtwW5V1bEJjTLBEpUMBwywaY7GjUtAUBnivvtsdpMPP0z+hDh2rC3j8/nnVvt3hfvlF5vLc80a\n+3dLhY5kP/5o64zOnm2dSFzB1qyBK6+0n6NGQYMGUUcUjrL2Jt0H2AicDZwTbH9JXHjJa/Vqu9cw\nYEDUkaSWPn2sybFNG/jhh+KPj8qOHRbr/fd7IoxHtWp27/CUU2x1+NzcqCMq3uOP20TWngiLtt9+\nMHEiXHyxdZYam9RVnXB4b9IidOtmX5LevFJyqpZoxo2zGmIyrvc4erT1lvz0U0+GJfX88zbUaMQI\nu+hJRmvXwqGHwvTpcMghUUeTOqZPtwUI2rWDRx5Jr06DZe1NeoSIfCAiXwbPjxWRexMdZLJZuNDu\nfaXTCvblScSaSzt0SM6VEbZvt96jXissnauusouJTp3sHmIy+sc/7AvdE2HJNG9uvU1XrYKTT7bh\nKJkgnmbSZ4C72Xm/cC7QMbSIkkTeCvapcF8kWYlYsmnf3hLimjVRR7TTiBHWzf7ss6OOJHVlZVnn\ni0cesWEY27dHHdFOv/xiTaS9ekUdSWqqUcPu+193nc269eqrUUcUvng60MxQ1RNEZJaqNg32zVbV\nJuUSYSmVpZl0yhQblOoD7BND1VaBeOst+OCD6C8wtm2zHsLPPGNDQVzZ/PST3Zfbe2/70txrr6gj\n2tn8/dprUUeS+ubMsXuJp51mt4yqVYs6otIrawea1SJyWMybXYitIxjPB7cVkdxg0H6PQo4ZHLw+\nR0Saxux/TkRWicjcfMfXFJH3RGSRiExK9OK+eSvY+wD7xBGxLtvt2sGZZ1oPvyi99JJ1qPBEmBg1\na9qYtf32S44p3DZtstrqPfdEG0e6OO446229ZYs1oc6bF3VE4YinZngoMAw4Bfgf8A1wmaouLea8\nCtiKFWdiaxR+BnRU1QUxx7QDuqpqOxE5CVsnsUXw2unABuCl2EnBRWQgsEZVBwYJdl9V7VnA55eq\nZjhiBAwaZFeVmbJwb3lRtU4XkybB++9HMxvIli22ft/LL9sXt0scVet9PXgwjB8PJ5wQTRxPPGHz\nq06YEM3npytVm3bxzjutsvC3vyXmfvu2bdasvX69bRs2FP+zX7/SrdFYqnGGInJ7vl2VsZrkr9g4\nw8eK+dCTgT6q2jZ43hM7sX/MMU8Bk1V1dPA8F8gKJgRHRBoAE/Mlw1zgDFVdJSJ1gBxVPbKAzy9x\nMlSFY46x1aJPP71Ep7o4qdr92Pffh+HD7fmWLbB5865bPPtKc97GjXDGGdaN3IVj3Di71/T3v8Ox\nx4bzeyzqmI0bbdKHVJ9PNVktWGDNpn/4g/Wr+PXX+BJYYT83b7am1+rVbdtrr+J/dupkzfIlVVQy\nLGp+/uqAAkcAJwJ511mdgelxfO6B2DqEeZYDJ8VxzIFAURN61VbVVcHjVUDChnWLWIeAffdN1Du6\n/ERs3GafPjZ925577twqVdr1eUH7KlWy/wwlPSf2eWn+E7n4nXceHHwwXH+9NVnG8zuJfV69esnP\nid1XtWpq39dKdo0b2/CLO++0OYmLSly1axef2KpWTY4e3YUmQ1XNBhCRj7BFd9cHz7OBeGayi7da\nlv+fIe7qnKqqiCR0oKQnwvDl9TK9//6oI3FhadbMvjBdeqpSBYYOjTqKxIpn5bb9ga0xz7cG+4rz\nHRA770N9rOZX1DH1gn1FWSUidVR1pYjUBQqd4yQ7O/u3x1lZWWRlZRUftXPOubSQk5NDTk5OXMfG\n04HmHmyB37FYLa4DMFpV+xZzXkWsA00bYAXWtFpUB5oWwKC8DjTB6w3Y/Z7hQOBHVR0Q3IeskcgO\nNM4559JTmSbqDt7geOB0rAlzqqrOivOD/wQMAioAw1W1n4h0AVDVYcExQ4G2wC/AVao6M9g/EjgD\nqIXV/nqr6vMiUhN4DTgIWApcrKr/K+CzPRk655z7TZmTYSryZOiccy5WWQfdO+ecc2nNk6FzzrmM\n58nQOedcxvNk6JxzLuN5MnTOOZfxPBk655zLeJ4MnXPOZTxPhs455zKeJ0PnnHMZz5Ohc865jOfJ\n0DnnXMYLNRmKSFsRyRWRxSLSo5BjBgevzxGRpsWdKyLZIrJcRGYFW9swy+Cccy79hZYMRaQCkLci\nxVFARxFpnO+YdsBhqtoIuA54Mo5zFXhMVZsG2zthlcE551xmCLNm2BxYoqpLVXUrMAo4N98x7YEX\nAVR1GlBDROrEcW6Bs44755xzpRFmMjwQWBbzfHmwL55jDijm3G5Bs+pwEamRuJCdc85lojCTYbyL\nCZa0lvck0BBoAnwPPFrC851zzrldVAzxvb8D6sc8r4/V8Io6pl5wzB6FnauqP+TtFJFngYmFBZCd\nnf3b46ysLLKyskoQvnPOuVSWk5NDTk5OXMeGttK9iFQEFgJtgBXAdKCjqi6IOaYd0FVV24lIC2CQ\nqrYo6lwRqauq3wfndwdOVNVOBXy+r3TvnHPuN0WtdB9azVBVt4lIV+BdoAIwPEhmXYLXh6nq2yLS\nTkSWAL8AVxV1bvDWA0SkCdYM+w3QJawyOOecywyh1Qyj5jVD55xzsYqqGfoMNM455zKeJ0PnnHMZ\nz5Ohc865jOfJ0DnnXMbzZOiccy7jeTJ0zjmX8TwZOuecy3ieDJ1zzmU8T4bOOecynidD55xzGc+T\noXPOuYwXajIUkbYikisii0WkRyHHDA5enyMiTYs7V0Rqish7IrJIRCb54r7OOefKKrRkKCIVgKFA\nW+AooKOINM53TDvgMFVtBFyHLdxb3Lk9gfdU9XDgg+B5Rot3va50kWnlhcwrs5c3vSVjecOsGTYH\nlqjqUlXdCowCzs13THvgRQBVnQbUEJE6xZz72znBzw4hliElJOMfVpgyrbyQeWX28qa3ZCxvmMnw\nQGBZzPPlwb54jjmgiHNrq+qq4PEqoHaiAnbOOZeZwkyG8S4mWODaUgUcs9v7BQsW+qKFzjnnykZV\nQ9mAFsA7Mc97AT3yHfMUcGnM81ysplfoucExdYLHdYHcQj5fffPNN9988y12KyxnVSQ8M4BGItIA\nWAFcAnTMd8wEoCswSkRaAP9T1VUi8mMR504ArgAGBD/HF/Thha1m7JxzzuUXWjJU1W0i0hV4F6gA\nDFfVBSLSJXh9mKq+LSLtRGQJ8AtwVVHnBm/dH3hNRP4GLAUuDqsMzjnnMoMETYrOOedcxvIZaJKE\niDwnIqtEZG7MvkInGBCRXsGEBLkicnbM/uNFZG7w2uMx+/cUkdHB/k9F5ODyK93uRKS+iEwWkS9F\nZJ6I3BzsT+cyVxaRaSIyW0Tmi0i/YH/aljmIqYKIzBKRicHztC2viCwVkS+C8k4P9qVzeWuIyBgR\nWRD8TZ+UsuUNqwONbyXucHQ60BSYG7NvIHBX8LgH0D94fBQwG9gDaAAsYWctfzrQPHj8NtA2eHwj\n8ETw+BJgVMTlrQM0CR7vBSwEGqdzmYM4qgY/KwKfAqdlQJlvA14FJmTA3/U3QM18+9K5vC8CV8f8\nTe+TquWN9D+Jb7v9YTVg12SYi42rBEseucHjXXrmAu9gPXDrAgti9l8KPBVzzEnB44rA6qjLm6/s\n44EzM6UCbLbdAAAEpUlEQVTMQFXgM+DodC4zUA94H2gFTEz3v2ssGdbKty8ty4slvq8L2J+S5fVm\n0uRW2AQDB2ATEeSJnawgdv937Jys4LcJDlR1G7BWRGqGFHeJiPUabgpMI83LLCK/E5HZWNkmq+qX\npHeZ/w7cCeyI2ZfO5VXgfRGZISLXBvvStbwNgdUi8ryIzBSRZ0SkGilaXk+GKULt0ijtejuJyF7A\nG8Atqro+9rV0LLOq7lDVJliNqaWItMr3etqUWUTOAX5Q1VkUMrlGOpU3cKqqNgX+BNwkIqfHvphm\n5a0INMOaMZthIwJ2mSs6lcrryTC5rRKbqxURqQv8EOz/Dqgfc1w97Mrqu+Bx/v155xwUvFdFYB9V\n/Sm80IsnIntgifBlVc0bL5rWZc6jqmuBt4DjSd8ynwK0F5FvgJFAaxF5mfQtL6r6ffBzNTAOm2c5\nXcu7HFiuqp8Fz8dgyXFlKpbXk2Fyy5tgAHadYGACcKmIVBKRhkAjYLqqrgTWBT26BOgM/LOA97oQ\nW/EjMkF8w4H5qjoo5qV0LvN+eT3rRKQKcBYwizQts6rerar1VbUhdh/oQ1XtTJqWV0Sqikj14HE1\n4GxgLmla3iDOZSJyeLDrTOBLYCKpWN6obr76tttN55HYbDtbsDbyq4CaWOeDRcAkoEbM8XdjvbFy\ngT/G7D8e+w+4BBgcs39P4DVgMdaLsUHE5T0Nu480G0sIs7Alu9K5zMcAM4MyfwHcGexP2zLHxHUG\nO3uTpmV5sXtos4NtHtArncsbxHMc1hFsDjAW61STkuX1QffOOecynjeTOuecy3ieDJ1zzmU8T4bO\nOecynidD55xzGc+ToXPOuYznydA551zG82ToXJoSkS4i0jl4fGUwG0jea8+ISOPoonMuufg4Q+cy\ngIhMBu5Q1c+jjsW5ZOQ1Q+eSkIg0CBZAfSVYNPV1EakiIm2CFQK+EJHhIlIpOL6/2ELJc0RkYLAv\nW0RuF5ELgBOAV4NzK4tIjogcHxzXMXi/uSLSPyaGDSLyoNhixJ+IyP5R/Fs4Vx48GTqXvA4H/qGq\nRwHrgNuB54GLVfVYbNWAG4IlbTqo6tGqehzwYHC+YgsHvAHMADqpajNV3ZT3mogcAPTH1htsApwo\nIucG51cFPlFbZWMqkLckkXNpx5Ohc8lrmap+Ejx+BWiNLaa6JNj3ItASWAtsCmqK5wEbC3m//Mso\nCXAikKOqP6rqdmxF+pbB61tU9a3g8efY4tPOpSVPhs4lr9gb+gL8j10TmgAESaw5toTOOdjq4MW9\nX2H7JGbf1pj9O7CaqHNpyZOhc8nrIBFpETzuhDV1NhCRQ4N9nYGcYLmgGqr6L+A2bCUBsMSWlzzX\nA3vne38FpgNniEgtEamALbU0JZTSOJfE/ErPueS1EFst/TlsnbjHsGVsXg8WOp0OPAXsB4wXkcpY\n8usenB+7yvgLwFMi8iu26K4doLpSRHoCk4Nz31TViTHnU8B7OZd2fGiFc0lIRBoAE1X1mIhDcS4j\neDOpc8nLr1SdKydeM3TOOZfxvGbonHMu43kydM45l/E8GTrnnMt4ngydc85lPE+GzjnnMp4nQ+ec\ncxnv/wEnXXVl15tkPgAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 11 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 12 } ], "metadata": {} } ] }