{ "metadata": { "name": "", "signature": "sha256:33c2ef1b18f3bedffe864251f1930f63055b0c6238517f2ac805ee9b33a159bc" }, "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(1)\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, '../src')\n", "# %load_ext autoreload\n", "# %autoreload 1\n", "# %aimport anhima.ld\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+/AAAIABJREFUeJztfWmsXdWZ5XrveQBspmDANhgMhClMZh4MZjLBQAhQlaj7\nT6vULfWfbqlbaqm6QpKqpGvqVJdKPZQqFSkiqVTUqiRABjKRhICZScA2Y4IHZtuYwdjgebz9Y7+v\n99rrrHPvuW+wyctZkuW373v33nPO3vsb1jdsoEWLFi1atGjRokWLFi1atGjRokWLFi1atGjRokWL\nFi1atGjRokWLFi1atGjRokWLFr/PuOWWL3QAtP/af+2/9l/7r/3H/74Ig0H3YmDevC+i0+mg0+ng\nzTc7+PKXO/9//NJLHXzjG3n83HMd3HlnHj/+eAc/+lEe//KXHTzwQB7fc08HTzyRx9/6VgfPPpvH\nd9zRwcqVefz3f9/BG2/k8d/+bQfvvJPHn/98Bzt35vFnPpN/7nQ6uP32/PPu3R187nN5vHFjB1/6\nUh6/+GIH//RPeXzffR3cf38ef/3r5bX99V93sGlTHn/ucx3s3ZvHf/qn5bXweOfODv78z8vnyt/9\n85938NBDefyVr3SwenUe/9mflZ/9hS+MfLxrV3ktGzd28Hd/l8dr16bvj/GKFR388z/n8dKlHdx9\ndx4/8kgH996bx/feW97LnXd2sGRJ+VyXLy/neP368l55jm+/vXzOn/1seW883rWrfO4bNnTwN3+T\nxytXpjUX44ce6uAnPynX57Jl9XPebY737i3nafXq8jnef3+5N/7xH9N+q5vjfsc8x3v3luNNm9Jz\njvGbb6bv5+fyzW/m8bJlHXzve3n88MNpjcb4pz/t4NFH8/jb3y739T/8Qwdr1uTxX/1VB1u2lM9t\nz548/vzny3vhfavjvXvLOd++vYMvfjGP168v73X58rTmYnz33WkNx/gv/7KDbduaPdedO9Pfx/ip\npzr4wQ/y+Gtf6+Dll8vP4rWr98Hj3bvL7962rYO/+Is8fu+9cp++9loHX/1qHj//fFq/MX7ssTRP\nMf7Zz8o5+853kjyP8R13pM/kex+RQmEceCCwbVsef+QjwPr1eXzMMcDq1Xl88snAypV5fNZZwLPP\n5vHFFwO/+lUeX3klsHhxHi9aBNx7bx7feivw/e/n8R/+IXD33Xm8cCHwy1/m8WmnAb/9bR4ffDDw\nwQfp56EhoNMB9u5N40MPzb8DgFNOAVasyOP584FHH83jq64CHnggjy+9FHj88Tw+6STg5ZfzeGgI\n2L0bFpMnl7/T5zhzJvDWW3l8+OHAxo15PDgI7NnjP7tfTJpUftbBBwObNuXxUUeV13L88cCrr+bx\naacBL76Yx/PmAcuW5bE+p2uuAe6/P49vvBH4yU/yeNEi4Gc/y+PLLy/n4fTTgeXL81jX5OTJwM6d\n+d52707zDgCHHQa8/37+249+FFi1qrzWxx7L4+uuA37xizy+7LLyXvS7GQMD5XjGDODdd/N42jRg\n61b/3l7Yu7f8fB3v2ZPWSGD37vRcAro2p08HNm+uv9Zjjy3Xpz63M88Enn8+j88/H3jqqTy+6CLg\nySfL3y9Zksenn17u25kzgTffzONp08rr4+uP+445njo1zz+Q5mjDhjw++eTqPn/kkTz+2MfKa+HP\nVuhz3Lo1XWtg3Tpg1qzyPXG9Omcvvpj2UuCFF9JzDTz7LHD22Xn8/PPl71esSDIs8Npraa/WXct7\n76VnE9iypbz23bvT/mmCESsUFcKHHVYKuiOOSBcaOPLIcmEedRTwzjt5rIJz9uxyIR17LLBmTR6f\ncEIpzC6/HHj44Tz+5CeBe+7JYxVW558PLF1aXm9c38BA+hcK54ADgO3b89/OnZsmKaDK8bLLSmF0\n+umloO2mYHQTHH10WgCBww8vN8VBB5Xz0g11m6EOg4Ple4aG8jMBgClTgF278vjAA8vnpILykEPK\nNaMbXNfEGWeUwmnBAuDBB/N40SLgpz8tx2yEXH55KSB0HlQJTJqU7yeUa9y/XutFF5VzfsEFpWBU\nwcfQOT7ooLSJGU3n6p130t4KbNyY1khAhcOuXaVwUCNCr+XQQ6v7mvexrk81iE48sTSuzj4beOaZ\nPNa9o8bbggXlvtbf63NXo2bWLGDt2jw+5JBsSAwMlHN+9NHA22+X18pGsMowhu6VLVvSswxs357k\niMPKlUm5BR54ALj66jx+9NEkUwJLlqT7Djz3XDLYA06hzJ2bx2++meRtoJdC0TXTDY0VyuTJpfAY\nGCgfoFphIwELcSBNAAvLGTNKgXPccVmwDw2lm96xI43VktfF8fGPAz//eR6rNXzeeaXC0Y0F5PtX\nQaoeznnnlYv+5JNLq44/S3HEEaXQO+ywqkJpat3u2JGE2WjQr1JS9Jpjfs4DA6XyVUGsG/xjHwN+\n85s8VmGkCsd5QKyATj21nEcWRjrn6o2p8uqGadNKId7PXnrjDWDOnDx+5520TwKbNiVPM9DLQ1HB\nqNcy2rEaIeopHndcuqeAGhUXXQT8+td5rEbGwoXAffflse5z/b16mkC+/xNOKJWheuTdoEKZsWNH\neg6BJ58ELrwwj9euTYq523j27Dxet65UEOvWJeUYeOutZKwF3n23XCPvvdfdCNE10w2NFcpYKAyl\nZ9TaOeuscvFcfXVJgynt9elPA3fdlccqID760Uy7DQykSQyFo6796aeXwujaa0sKTa/llFNKSo/p\nFVW2J50EvPRSHistcOSRpaJkqFegirIfhbJtWxKE4wm2+IB0b2z16b1fcUUpxFXoz59fens8p0B3\nBTRlSm8F9MILeXzlld09IKVVp07NSmX69FIp6Hrq5rH0Q3lt3VrOoVMo7LFs3lwqFLU2VYGMB3RN\nqFHR6dRfg17fAQfkPQxUmZJZs0qPSan3c88tFb9S7cx86HcrM9ENSnkxXnopyYTA8uXJeAGq9FdT\ndFPknU5Je+7dm/ZJQBWGyolxobxGAhayQBIGLFgvvri0NnRyL7mke2xC+d1LLgGeeCKPb7kF+MEP\n8liVxJw5wOuvp58HBkqFpwv13HNLj0WvVV3v6dNz/EEn+9RTS+tVFz1Qv8Ec5bUvFYrey4EHlt+v\nnLpy5OoFaKzirLOSCx/Q56zU5fXXl0aExmlOPLFcc0p58Jw7upGVoXqaun4ZJ5wAvPJKHncTRv1Q\nXu++WyqM1au7K5RNm9JaDPRjbTaFKoijjiqf2xlnlMpVjbHjj8/7MN7P9Ld65aqc2VAESgUW9HXM\ncSiJOnZB1ydDPZRu8UulvBisQIBS4L/wQnpeAfU2tm0rqbPxMAZUAe1XhcI36PhTpp1USOsG1kAq\nUKWeONg/MJAedgi4GTNKukhjG5/4BPCjH+Wx8rkayOdNo5tAuV0VNt0ElyqUblyt3v++VihKkWgs\n64wzSqt/3jzg6afz+JhjSk5bn4V6GWqR9uK51ZNUBaTBdZ1zFRosvFQYqTJkQaixiW50ic5hNwtV\n6Yp3302ePo/78VBGCt6TGt9UL/SCC6qBed4bl15aGoK6l3Ss+1iNCP18jZfqmpkzJ9NsGqg/+ug8\np4ccUiapHHpoSdcx2ENRoc8xE52PxYtT0k/gscfK+MmzzwLnnJPHq1en5x/YuXPsDYZxiaE0gQpC\njRVo5pUKC6Bq7ag1c/PNwA9/mMfqhdx0UylA5s7NlmK4ebHR58wpA4iadaRjDTgCZdYQC/oLLywX\n9amnlhlJjOOOK5WTjhkqqDRZohvGQqGoQFc6T++zSdKA0mIqME46qVxHrOjVy9AECqVAzj23VHA6\nxzfc0J3mYoV58MGlpaxeKuP446seSqydfmgnjZEApQJyMZSx9lDUUFOvVOMeqkyVtlR6UAW+xjN7\nBfKVutQ51rjKDTdkmaHKXLNTGR/5SJl4xGAPReNYW7bkOVGP5K23yniIZnQ99VRSkAENyL/8ckmn\nbd8++rjpfvNQVODqwtGAHFBVILqYlFPXwOu0aUlgxWcoLaFxFxVWTE2p8FO6RikzFZ58L8qpq5XE\n/LsqCadQ6gTOvvBQ+Ls1a0a9UGckKPT+9LledVVvmosFggoYTtYASi+DUzWB6pxrUPiCC8o0V/1u\nplfOOadUVjzHSgNpxltTqAei2LGjpERUmI3UQ+E14PY5KxQ1LFVI98oe1P2gsTBdYzqHuveUvtY4\npNKRnP2nFCxDM/8Ye/bk5/zmm9WU4QAH5F38hD8HqKb8qkLSDC+NsY1VjKYOo1IoGmxT17fbAw8o\n3aPWhaYTDwxUOVTmPQcGSiXBLitQFV5qrbDAHBoqU0fVHVbBp64+KxhVrqed1t1j4QWumTCM8VYo\n+vmzZ5cKpVv9RUA5ceWp9bmpwNA1oBlVaoHedBPw4x/nsc65WryzZpXrljPPomYp1oDG7dhQ6Jbt\np0K0WyIG0D2Goh5KNyjlNRIPRakerUWZPbt8fg4aT9W4BytfoJr9p7E6jTupTJg+vVQiHDsDqtQx\nXw8ngqjyAfLcHH54vYfC6KZQuEbkmWdKOquJd7FpU7q3wIoVZQqypgyvX9/f+ukXo1IoqjAcJaRg\nugKoWpdaiwBUF4/SEsqLaw0KC+jg7GNRqHDRz2baLgJ88V6lzK64okxTZaWhAkWFKFAuVH6umuvP\nGG+Fotei9UHOclFDw9UKMPWpBWlAfzSXWqhq2Fx4Yell6Hq58cZyzq+6qjRqNC7EBovWPzFcYVyg\nm0LpZg1u2JAMDKBatOiglJfzUHrRbb2KG+sou27r/ZxzqrQWGwkad7nkkmqtF8dNNFtQ51CNjuuu\nK9OHr7wSeOih9LPGdhlcX9fEYAZKr2L9+rLmA8jz/eCDZfxE6000w8/hgw+Ssgy8+mpZ1Kg1KGMd\n1B+1QmHrVN16B00NdhlOmj1x9dVlZboqgcmT06SEBaSL97bbStqLabXBwTILRDPJ1IPpVj2rwlaD\nheyxqMCcNSu/VwWKWuiMfa1Q1Ntw0DlVKsgJIK35UC9DqSatuOaMPSDdZxghGkzX2MeJJ5YZWRps\nv/76kna95pq8HtWCZQ9OPVpGLw8loM+JM3CUb3fYvHn0MRRVKE32uSatqFGhMcaLLy4D85qx2Wus\nhaZqRKjMUC+Xq+S1ewUniXA8qB8PJeaJM7w0Y0u9z8cfT/cZUA9m164y/dfB1aywt7RtW3022kgw\naoXCD1S5UKDq2mqMxFljar1oxtTAQNVF5ljLwEDS0iEINTCrnLwKfnbPNSvpmmtK+sQVsNUVRnGC\ngfLETrEGPkweiptjjYOpRe8SBzQV9NprS4tRg/8qAJTG0jm99trSItVsM1WMkyfndRoCJeZRuz6o\nMGQDiC1rjRnymm1KeW3fXj9nyo/XfQ57Mc5D6cWPq0Jxa4CvGagadarwdU/qM9a4h9bqaPKFPmvd\nX0pfq1fs5irWLBs7HG9pqlA2bsxeJSuUp59O6xJI602VA78PqAbkV60q4yUOGoNxVfJc1DhajEqh\n6CJw0KCtBnWBKl+q1a9hIfCC/fjHS9pClYJ6JczzKies7rD2jGIFo96C9vVia1c3qiYMMLoplA+T\nhwJUrWYVDi6jTZ+FFiyq4ga614zofWswXVNHXWYPex0LFjSvkg+BGs+BjaRuWUGcjNCU8tq6td6C\nfOONMmW0CZp4KGogaAGyg64BLRp1Sktf03Wl16FUqsZV1EjR2I5mjKqRwazBRRdlw4Dnlz2UKVPq\nWygp4l5ffz2tA6AMyC9blqi2gKOi3n67VAbPPVf28FJ600HjORs2VCm40WBMKS8HFZRuYak76uow\ntBGgFj1qMF4Fwm23Ad/7Xh6z96DB1CuuyHwqUE1LZBdYBaHGUXiRa0sWvl7lpVnJqiXEG2vq1GpW\nVSxGzegYK4WiUKNBa0uAapM/p1xVYGjRonqqxx5bKhEWMCr01eJVOkpjIVolrwYMx8cWLMjrRedD\n6ZKwbpv2YOulUNhD6dbuI9Aky0uVtaYJO2iml9KKQJWt0EC5xmDVMNEqd42jaJKMxknUqNBYGhsZ\nLF9OOSVfhybI9BuD4Cr1t9/ORYsPPpjiOIGXXiqD6w5aILlyZW+PRT1e7eM1Woy7h+L6VgHlRKg1\nCVTjKMpjDw6miWEBpEqDXXUN1GrgnrOwVCCqklB+lu9Hs7k4mKhUAVcQa1yBA/76O6550TYv4doD\nVXd3rBSKGgUad3JQDt0pHV0H55xTNhJUulHjLOppKr3KveC0MaBSHr0oN1Z2LpGE10PU6bhalF7o\nplCUEnE1KgrnoahQ1LRbFfzxHn6fKhSHs88uU3B1HzUpcGQvUo1KVThqNGiMU5UD/55ZDKXPxgPa\n1FMbQrrK+127yp5gmjLsaDTFh0qhaMaWe00L/oBq3rdTTPPmlcLEbVqtitYaldtuA7773TzmWIQK\nDBVOXBCpaaya2szfq8JW23ezG+8yvQKaOsxwqYwBVihakNREoagid/USWi3vrlXTRLXLLFCNg3Gm\nDdA7m0tpLs0YVA/H9QljgaVtWrRKHshzp+noQBawfF1ctKcp190Qn8UKxdUQaFFjtxoVoJmHon3u\nHKtw2GHlutAYCVBdJ1rXpWPdD1oQqXJCmQU12DQrE6h6SUqTcTKHUm6jBWflvf9+Tvd1c7JqVVLS\nAY631GHVqrKocc2a3pSoKpTRdlLoS6G4Jm9q3RxzTO+c9G48c0DjKEA119z9DW90zd5RD4YDs6pg\nNH2Y/1bdeXW1WXHpoufKchXEvIA1S4bRjYZShcIWinosDqo8XUaWPivtRA1U+WpX8HjhhaWyVUEG\nVGNj+syU5uLAq1qg2tpcDYNeVfJaZc/fzed4cGCeqbGmlm5dDEWtWEWvokegWQylWyPLgNJTLlCv\nXou+R+dbn08oBMbgYPk9GkfRBBat2tfYndJenMF35pk5DqR1K02pLu4s/OqruSZkyZIcYH/qqTI9\nOMD3zn8PeI9l586ybkVThh16dRruF30pFJetoxOu1JKDZgEBSWDzQnC8rWbuTJqUJpYX4ac/Ddx5\nZx4zd69CWmkvtlb0PpR/5doYdaW1cI+tbC6E02fHyrhbsZhaVYxulFcTaIKEg0usUGhQ1kHTc4Gq\nh8sbHPBFi7wmNENQz7nRdvis5HQee9WvsIfMcRQWYt2SJnoFpYFSoagHou8fSw9Fm1UqNObnoCUC\nztNR6kw9Wz2yQg/gUppM14tmA2qsTGlV9nKZYmO6jhmCXhlynNrNMQ9WEA8/nGj1wMaNZS0JkJ4B\nt6DXCnkHPVjLQRXIPlUoLpCom8AtNKU23Odo4RJQpV+0DQZQpS30+//gD0raixeoejBa4MbCVa1b\njauwd6FcMVs62mmYwQkMmgvPnz9SyqsJ1KNycNSNbixNQHDv00A5UPX2ND1UFb3WHyjNpdlc2npn\n7tySVmFDIazhusA+10JwzMu1GOoFNdYc5cV1Co4fH2kMRefOeYqKOoXC8sCtdZ3POXPK+VSqW+dL\nE2Y0jqKGjFLl6inrYXrs5XLCTV2mV6+CzroalPffz/EvLUZ84omy/sRBe3i563BJGy4ex/O/zz2U\nXimqTtg0CdqqpQBUg2yauQNULRCgdLW1mv3WW0vaiy1idY/ZAtWFp3w9e13q/vOmULeeM+W6pQ4z\nzdUP5TUSD6WXQuEizABzz4APums8Cajy1FohDVTreQ46KAs8XRMas9OqfE3z1voVF9hn44DnK7ph\nM7jOgX+Oe2TamPcT11nUUV6sUNaurbbz+OCDsg2HQxMPpQnl5bxUzc50a0ATJTQRQxWEpvYee2z5\nvWqcOo9BvR6tQdLvcMXLTKVzDNixNkw9cZquoyx37KgqeK032by5KuTV+1APRq8D8McQK/Y75aXQ\nRm5ANfMJqFIrzqrTIC1QDbxrKxWget48p5eqcmOvZGCgXHzdFrder1rWrKicNRfXy95Lt6At16Io\n5cWbqFsMpQmaKBRXLa/enoMKZ6BKf7o4g9Jc111XGhFq1XbL5lIL1Xk8Wr/CHo9Sn2y8cPo0t/Xn\n7C6mNbkWxZ2JAtRTXnVFjXqwkqJJDKWO8lJqSpVpk0wvNSrOPLNcE+pRuO9RmkxliQbzNQ7XK32Y\nM0pdsgF7tc64Y0/A9fF6++08j7/+dTJOGTt2lMkzTz9d1qgEeK41w8tBixodPnQKxcFZ3nroDlCl\nuFx7aKU0gCrXHZMci04VDB+rq9lHTGXpeSzqDbHg0uJDddUZLHw5s8Vl8HBWUV2gn7EvYiiTJlUD\nsFqLAvgGgBqYdwcaaX2JFs0p9alCX6nLXjRXt/oV7dekApCFDxtA7GmxxcvKhRWKHgMcUA8lzj5p\nUiXv2mo06eXlKC939oeuV6dQlO7WOJVTGAr9DM0SdXEU9jJdjz2m4pTKZE+Li5xjP7JB52QUn4Xi\nPAcuaHzkkfQdgd27q/3Z1GNxaKJQfic8FOeRKBxXrq4vUJ1ooJp66IL1mg4KlMFvFfZ8hooeDayZ\nY1xQqbELDQCyAlAlx8KDi/q6tVXhnkndquXr6lDGK4ai3wl4WtPFUfS9LolDC9KA0uOLzLN41rq+\nlLvXbC6lubRXHHuNQGmFaioz8+ycdszpwvwzC0NVKEwnuxgKKwg9qdHBVU6PNMurSXGjO0TMtSZS\nqCGpe0LTizkBAqgmd6hx42gxpVq1eDjS5ZkyjVM4WZE6hcIeyp49aU/yXHAKMJ+NAiRjhc8/Aaod\ngtnDCbz+eu829R86D8UF07sdMtMNbvFpXAKodiEFqplAkUnDgkpPY+TNrimt7HkElRWfpQWVXN2u\n1ipntag1zm53t7oV3lycldZN8TD6USiOEnGUVxOvxfU1ctl8TdLKtXEmUKUpLrig7K7ANFy/NJcm\nhKgXrN+t6c7xzHgumbpjhccUSh3lpd0NuIFgQNvSO+hZKMDIKuWBaj8vB+dtuForNUz1KAfXCJIp\nK2U86uIo3U561biJdkIIo4bT5h1lW0d5qVBmDyIMA9eeXk9o1PsAqgF5oHpO/Lp1VeXRpIhxvwfl\nnUJRqsN+salv0OArUM3qAHwgXif7wAPTNcSm1hoUdluVEmELVTldtmY1J17PP2dXnRewixNwUVwI\nHP7ZFYg69BtDUYXmlEeT9isOjup0QXflqZ0FqQkaCxd2b0Gv7cy5UBXoXr+ia1ppUW1kyRYz11q4\nwjiOvXXzUALa3LEfOKXTxENxe9MplLq4D8MZB5qA06tDgjIDrv5NDSFVUuqFqgzR2rhzz81eUSgj\nlwrfi/IKRIYXX/fjjydDk6ExF2cUOIWicCnDrlGoyof9TnlxPCLgGv05QaXBNZeh4c4cUBcYqHok\nQJmlpS77Jz+ZjxLWo4G1wI29ErVO2b3Wgj9use6oQdf2noUwF/F1q5BmsLIaSQzFUV5OoSgV6eAS\nLZzFqs04Ad+Vmo8Z0OJSVV5Ki/TK5lJvirtVB81VdxonF8txHIXbrgBVIXjwweVpofyZznvsFnDf\ntauqKBzlNZJuw4BXKE1qzlzBowbmlaJyno4qEU0MUa/m6qvL+dV+gU5JsQHARx5EzzD2dEPJOPnk\n0nNffjndJzeHdN6IQptGAkmBRSwN8OfiNKlBcVln+12hOA3tDtpyHPvs2VUaxwUEVXjp+RZAtYMw\nUBUazK1rgRQ3mtNFwpw+ewBAtZqaF6Vrix8CloOznK7MgtHFngJ1cRN+VmNFeTmF4mI66rE5uFRS\nPVMc8F6oUlNa/MkZeko5qpHTq03L9deXHpD2ImOBxt/L64oFZ13tTswBB+XrBHy36mwO2AfqKK9e\nvbwcXMdht8+btHR3p7AqdM+rwaAeqMZRdL6dYmNDEagqneiM4NLdY/+69lMhlDmOEX23OCCvSRNr\n1pRnlwDVCnkHVlKBJgrFdRr+UCoUZ7k4+sMF5p33oeexA/6sZ20GNzSUNk8ISD1fnheDNh5kgaGd\ncZma01gQ35MqOXbl+Z7YOm5yeJGCOeR9RXm5WhTN4AHq6U8WYo5icf3dtCpevQ5NM9d0dfYCNJtL\nv68XzXXDDaUCCkqNlQQbMJzp5aiwOsrLCVtu5xFwbVfqKK+R9Gpibyrg0uHVK4v31hlGAc3kUtmg\nBcz6e/f89Dlr2x6lwRYuLOf4uutS7IyNoFhDkVzh6LegvNxxzc8/n7x01/Tz8cfTfTLUG3HB9hUr\nyq7DfA3d4GIqGovpF6MOytdRXm6hqULRoDbgO/m6OIoGSoEyeyvANQzR2ykWwM03Z5os3Nn4Hb9P\nPSK1ZnnhakYRB9k1MB9KhD2iJvRDfG58JwdR+0kbdpZpU8rLWdzaXBHwZ6PMmVOmBQO+nYx6S2rQ\nRNZNQNtqaFxFvUml35Tm4nlVgcotzYEy0/DII5NRwFYxe2EzZ1Y9817xiB07cgB3zZq0xxiu7Yqj\nvEYak3Hrsm6fa8xEs+aA6v1qnzSlsJQWU6YAqColjcVo9qDKGjUAuRdcGEZxL5p6zgjKK+IhrNQi\nhvXoo1XloeebOLzySlr3jCYpw5FtxhjrTsPAGATl1W0E/FkIjpJylI429AO8QHPcpbNqtWaBg/dn\nnllaOWz1aNt9XpxafT9vXn6fBvH5M7nadiTtOeoEXJ1CaWKNOg+liUJxHoqrRXGZXkofAdWzLYCq\nVQpUPUL2qIJSDEWpgltPWtQMPj1mWBMIjjyyzCRjWoY7DC9YkBVbGBRc26B1FEA1s0qVPVu7LmXY\nKZQmmWCjgWv46RJrXNxMFYjOjRosTqFpJ2ylMd3hd6qU1FPUM55iHIxC7OduZ0GFdxB1H2vXJm+O\nvQvXXkWF/gcfNAvIN2m541KN97tCcdlJriPoeMB1tHXnNLgOvpwhdfPNORg/MJAmPgSzUihs8ahH\nxGmpGgBkr0StZgZ7PdpPKAQVW3LsDfI9uYwloFQoTavmnUJxRW2uWp6z0gJ1Ra0qYBxPfdllVYXi\nmkOy0NCD2Fz9Sl3RotIoSoEozaUFdRG3YbpMD3oDfMo8Cze3n1ihuKJGJ1SanOA31nAxMi1UBarU\ndq8jkYGqENfsTB3PmFE1MJUGU69a9/miRWkcMZrISOwm89hDmTkzrceoyQpPQtuiuKMlli5NdDvj\nN79J7Ia2oW8VAAAX5ElEQVSiV4cEV9Q41sf/AqM8DyXQtJVzkz47QLJO1JV2sRVtmQBU4yRAahAZ\nlfIR9I9r/sQnshJRL4o5e/WcmL5S2o+9Gc38YsuLPR22WrltB1u+HMxkr6GbhxJKxC1YN28usO4s\nOWdcuNdcxo7zzpwgchawBlu1RuGGG0qjQOsLtP6AhZTSXJrNpfSdWsBR9c/Py6VJs+JtkmIPlB7I\n6tXVcy6UaweapQjvC7h930SBaKdhDbwrZeVaxuiRG9qNWulr9abC+ImssCbGczz3oLwiZTgC8k7R\nc7A+4ALyeuKiw/r11bXgihr1gLaxwJgolKY44IAqZab54oDvPKztE4CqcACqbTqAaiEjUy76GVyv\noAuYM1g0LZUFg8ZRgCy8WRGxkuSFzMkF7HFxcWMThcIxlLrDtcbau+znSFT92ylTqt6Ryy5imlWD\n65oefu65ZUKHs0BZoKiQ0oJJ9h41psNtZCImxLHDEJAcQK4rGOxGebmA7miDqSPFaNaPvlczxFQO\naBaWa+WjcRRN89fkG00gUqOCOyPw603WeQj/t95Kezc8FG23AiRPV3t69Tr/Jr5DiyNdhpfzUMZj\nzYybQtF2CkB5uFTAZXrpOQqAj6249GGgen45UKaYMrU1MFDSSrfcks9IUWHF1owmKFx2WaY/VKix\nQuMCKlYiHKRnmoir5bVBZFBeTWIoIzn+d6ygmxzwB4hdeGHVC9XMLqDqmarVyQF+PWlR43Ya49F0\nZa1J0lRwXlcs4KK2hjPYNFUd8JZ6HeXV7ayTfs82HwnqCjVVqLusrkhUYGiLF67bAqrFrM6zVYpW\niyQ1ocfFWZUlUTYkjIwwfsNLjvtski4f2WCDg94b0Rb2Di6776WXypMdgaRQ5s4tX1u3Lp9fP54Y\nN4WiVA/gz5d32V+uLXi8rgtKO80CPtuLD96KTKZYBHxmhh7tylaRusPsks+fX3pQrFC14WS8zjQP\nW6286LlNeDcPpS5tuJtCGQ8h5OqIXI82FzNxB2659HBN49bmkM7r4O/SbsQ8V0qbqHfr2t+zwol4\nAV9jeDWuLXoT6gcYH767G5xB6Iw6V8TsUoddppcqDPUmHQ2qBcLqUWoih0sGUgNHY2FqVIRXEx5o\nGMFRi+JijAouKNWz4DudqhHBRwQHfvvbZJQwXIaXO6lRYzbjhXFTKE1TCh0/H6+rUnF8dOSKM1wW\njVJhvIhUuDHNwV6JCh+2lNVj4SwWTTFkS9Z5WAx+PpzVMhYxlLinkcLNkWtj7zpLu+yvOi9Uz8AZ\nHCytZW1Lo+no2qZFvQ5NUVUPV88wYWGkmUNhnLA1HYpECzE7neoZInVwaZ/jCdfPq2lxY13qsBoG\n6gm4WhJVbEp7qRHi4iiqQFTpKH2tSig+M2JfcS8hZ7r1M9y2LX1/1J84ReGUwtKl1fjJ889XM7xW\nrEiGOoMP79rXGBOF4rhMV9zoiqPq4ASOi6M4IQT4tD7uHMpW7cBAmX7KbVn0mtnade064lq4a7EK\nbd4UbOH1amfCn8NxhNHEUEYD9eYAX4vi+jm5rD3A0ypOQCm10S3W4YLr7DFpgZsqHD2HXAv62Grm\nhIqgdMIQ4vkLWrTOQ+l0SoMg3ruvYmFOMLuOw6640a0BpzhdDRtQGg+6xzSTy3kgSsPpe1Sh1LWU\n4utldiDKIsJQdGn1YSiuWpUEftBcDz2UZAPDtWBZsqSa4bVqVdpLDO1WPNbol8UYE4XStP1Kt/fr\nwnKdh+tcS1cUp+fFAynbKxpETp2aNmwIX65BUKuT0x45qKddTTkYr1wtU1ecmcQV80wJuuw3Bhen\njTSGMlrKy1XGu1qUbmde6DXoAWqAz+bTyvWbbqrGOtjrYO8gDlILhabepVI2eqSBZpJpL7LwaKK2\nxtVmROqwE0asPCL4Hs/JZXj14u9HAkddugQCt89dAXRTaGxN59HVn6khpsaGKhSXPagZZVqTFJ4n\n7+N4D8u/mKfNm5Oyiyr2aI/iMreizxfDZWDtj8SLfjMFx0ShOEujSbvrgGsL7U5KA7yg1dRQwOf/\nB2UUk84FaHo2OQcWufGkWm6cVKDU2eBgXujsXbGyqQvSs2B2fDaQ76NJHcpoLVsXfHXV8nUeowNb\n8wHt0QR4b3fatFJoab2Leh1aY6RKQGsleJ2pANLr1sBvCLBYU/yMo3Yq+HfXdibgsrlcDcpYpH/q\nNdSdiaJ72hUx9/M9WpyoysDVsChNrBXxzgPRpBA1YtVo0c7CEUuLdcMGSWSn8d6LNN2VK3PQPJqb\nugPOxrOWr9e59wHXZLLfVj3j5qHUPSDHzWrr6G7vv+iiahxF+y4FXLYJpyKytzE4mBZJCA6mPdRq\nYveXi9/0JEM+y8VReJ1OKayYimEB6Vp1MJgj5qB/kxhKUzStlq+D6xHmYmKOQgE8HThnThkrY6Wn\n2T8a21GvQxWO0qtaMMvfpcHj+GwWZCGA4sApbg7q0OmUCiX2g1MovTLAesHVD9VRXk2NRAeuqQpo\nYF4FuZMDmi2qRohTctrRo1e3aY2XAuk5nXJKMvr44LDwmvg5Rg1KrPkDD/SJFe61DRuqBsLGjdUs\nMPeaq653TULda64x5IdGodTB8eF1gUm1XgBfee7OzwCSZ8FV1UCpKMKVC0HFVopaPUwBcFZYN2qL\ng31qibKSitdZMNWlDvMij82mZ8rH/fSKofRDeTmF4qrl6+AyvRy9VQfN1AKqVfPaHJJrioBS8Wrc\nTwW8KhTXiJJ/z9XvcRYPf2cIrLjnbp5czGcoFJ6n1aub9fHqB64vlqO8nNdSB9ebzLVg0cC8O176\nsMPKPaYeiItDKZPBbARQTTV3h3KpwXPZZcl72ru3ZCbYQ4m9x+eaxAmNDz6Y1g3jiSeqZ6IsXZq+\nmxFBfcbKldVgvksZdufaN+003ORANsa4UV51cK046uAaRTolA1QtHcB7PlOnps0T7ipbKsy16tHA\nXBDnToyLzcMph8r31h3q5epmWOBwZTUrX3eWzHhRXk6hOMsW8BSd89Dq3u9ojquuKg8vA6o91ebP\nLy1VVQL6GWrpshel6am9vksD+TNnJsPi4ouTEArh5NquMHhOQ6FwjYJL/3QKxdUs1MGl6Tvl0Q8t\n446rcAdU6bkkQEkVA/kZBupkCD87rWnRXmGudkRr4tSLiQP0pk9PxtSqVek57dqV5oXXc5yY2Onk\ngDwf/RtwMZWlS6tnoLgeXk1ThuvarjRRKPslhtLNQ3Fnj6uHAngh5CgRwB9WpX2XgLQBXENL7jrL\n1MfgYBLI8dlcVOcsqVBKLKg0g4itFras9NyMCMy7IkCOJ/U6Cni0rVfq0M+pja42wRW11l2HO3Cr\nLkOQBZ82h9Tr0OCsO1t+NAqHlWDUwkRaa1Bm3WImStOFQuE+Xb36fAUcpVGHOg+l12mMAaeQHBPh\nMrKA6vM444yyXkf3nnsG2hZHFb5LDNEEEpUhSp3F/Fx6aVISe/eWnu3QUJYJGzakn6NZ7FFH+dRv\n7iIdcEWOzst4+eVq52FXJe/arvSjUMbdQ2kSwAO89TF9uj/3oe4cBfe5an0A9YdRcVvxAFNRQ0Pp\ngYUnsnBh6bGEYIvah9h4nErKGWlBv8Uz4s3A6cJsITG9xtysew5cLV93MNF4tF7pR6G4TK86b8T9\nrbNkAd/ZWluSa/GgKpzdu/PcuNokroJv4uGw0GRjIKzomLtuzzmMKa7D4BiKUxiM0XYadgrBUV51\ncO2O6gxHB83e1MC8S+vV9aiUlmMyegXiXRW9Zn/NmZP2YFwf17xNmpTmPzxDPqfEdfvtl07SNaQF\nkoBfC6PxUMad8nLeQd1m6Sd12LVgAXzhX10HX5d+rO3rgbSBDzggW2BMZ/Fn69HAHNBfsCALFxWW\nbC1xBpC2jgiqjIO33InXBbLZQ3G/r6O8nCXUD0arUOrgWtnXWfGuR5p6sUo9OYXDcRtWODGPrHCY\nilSFw7E0ICkcvr7wTsObDQo1vObI+ArKONaDxlB6Bd137EhrgdFPp2GX8NBPvMTt8zpj0H2XHqCn\nCRBAlaLSuXBJDhpHUYPAdcfW2I82j7zhhpQ9uH17UjbTp2cKc2goKbFIpFi+PCmfQw9N33vVVeV3\nPfNMlQJzQfp+0/tVHrvgvfue/UJ5OWuhDnVnTrsHxEWHDLdQ6harO3RrcNCnvHLxIi/OoaE0ISH4\n+Xdc+6CppMyNcxzl8MNL64uVgAs0M83iaKK6fl6Bbueh6ELrZ6G674rvUwHhYiCApx/r5j3O7GY4\nzzSUTwgbjedpKveiRdVDt1jhaL8tVTjs4WgdjlJqYYDwWRrPPZd7QsXajqAyHwMMVCmvzZv9KXxu\nHvv1UNR77HefNzUcHROhrXScgapJHa77uEIzQs85p9pMVql2VTq6HsIjmTQpGSfr1uX9OGlSOUer\nV6ffn3++j4E89lg1IO9iKu5AtdEqmQ0bmiuUcfVQxmKh6cYB6g+ccp2H4+9VSbgsEsCfT86LKzyW\nuK+rr86WJr9X752VHVsyKmg43sL3w1RZNM/jze2Cj0xFsNcQC2y0RwDHZ6iSqKOsXLV83d+6TK+6\nv3UHbtX9rSpm9lS1OaSekaEejlKk2pxSz1thJak8fQin8Hrju+L0xhBOrFC6UV4uZbgO/XgojvKq\n8xJdNmU/CsXtUVewrMaHeiSuO4fGsZTJGBqq3hOn9gPVgmqXPjx1arqPTqesnxsaSs+di1GXLUsK\nwp2U+fbb1YaNy5Y1C8g7GqsfuDb4H3qF4rKXAB+wq0NdaqnyrEB9+rA7+W9goGyvwNXPkc0BVI8G\nZm64W60KZ6owrcNV8szzqtsPlM/J1dQ4GkpjKCOprHWHbNXBVcvXwWV6BXSj1yVkzJpV9Xq1V9eN\nN5ZKQZ8tt0nXjtXaXkQzDdXD0UwgbnYaCiJiiVGtH55neLWhUOoor9jorkoe8EbBaIPydXBxUdep\nAvAJJpq0EujWngiopoAD1WaxmurtkgA4rR+oFsA6o0XXwJVXpvXywgulcRIeypQpmSb64IM0fxoU\nr8OmTdVeXy5l2GV4jbbWbL/EUPppq1BnUdal/TmhWdc8ks+eYDhOPmgstcI+9Sngrrvy+0KATZqU\nFkxcO1fyc22LcrC8wPk6Lr88x1FcA7xOp/SYIvjM1J4rCmMayqUN8+sOdW6zqyyug6uWr0Ndppej\nQZwFCvhsPk304PPb4z2scNQLUSuV6TY9wkBrptSi1c+OY4M5BhCeTBQKOg8FqBY29uOhjDYoXwdX\n3Fi3xlyjUNf2BKgqJU3zdd+hpQLueABVOtdcUyoQR+VqarcaDfPnJ3nAhkeUI2zZkuI2M2ZkQbx4\ncVJcjH69zSYNJV97LcmJkWK/xVBcllY/qMsAceegAL5tfV1ao/Z4CmhHUaBsoRHHAYcA56wR9l60\nap2teRYu7OVoERtTStG/i4WJ687rAoi9PJQmcBu1Hw+lrlreZWTVGQeulT3gac04SlVxzDFZ4UZz\nyFgzGnNT2kU7VmssThNG2PCJ512XORY9ocIoYWo1nr0G5QO7dpWbee3apMAVzjAYbVC+Dq7jcB36\nYSLUi3Q9u/TIb41ZcXZlQOlQF7NTmaaBeJU1TG0ODCRlMn16WqubN6fPGhxM93DUUf7YXtcQ8t13\nm5/x7gpcXcpwPx7Gh4rychlHdXBWOlCf6VX3et3hPi5gr65tIDJKFi9eXJw3z21ZNCWZrSnO7OEs\nJH0PZ5ywZcXccGwGxzO7zJeRKpTFWiEo6Gcu66rlXWFbHVxsBfDepqtoBqpV81HLEvfKVEfw4iEU\nNMVUuySogtHMMa2bYMUQ9xbZfnF+D6+H8OQ5tuiUhLMW2Yvhee2H8qpjEhxcx+E61BmOjgrTs1Ac\nuMBx8eLFtR2M+TWNozgDSqkyd/ibHhB20klp3RxySG40u21bkj2bN6d53LEjH6Sl3xtdExhLllQr\n5HfvBl5/fXHlml1MxlXJx2mRTeDW17hTXnUKpZ/2K4CfWOU3A67zMFAfsD/iCN+KYc+eqoUcB28t\nXry4UFzx97HRmMPldFE+j0EVJVs+bCkxXcf0X2wGVh4hnNgzCo/NZUL1o1DGgvLqp74E8C056q7Z\nHbgFeE9WhUAo97hXrSnRoL+2v+eaI40bqPC7/vqS5mJaJZTXjBlJIIU35mjfWD/9VKRz3QErlH7S\nxPvxUOr6ebnYpYu3AP7Ybxfv0BoQpqX5XnkdKxXuAv7a/VkD/kCVFdFD2xYtyt7I9u3p39at6d+m\nTcnQ27Ah3YOLe7nOwS4gv3IlsGnTYjTB2rXVQP1og/f7LW3YFQXF3zelyOo2Uh1Pz4WHjEWLqunD\nQPUgJaDkxAcG0sKOTcDFi3w0MLd00epsTlvlpnV13Ya1CHLp0vI5BG3AWTfBuYfSA7JQ2NeUVx1c\ny3bA89xxHU0FksvaA0plpTEY5fNdN2IOtmsgVjPHNJDPQlY7V7MgjKy4KMDk/eGUbS906+PVVDH1\n66E4heJOaK1DXTamQtPEHbOhnZ9dCyY9e0djIi7DVNu9aDLB7NnJw9y2LSmsTZvSz1u25I7mnU5S\nVBo/2bKl2kUaSApKacrnnkvGDqMue9N5La5Kvh+MO+VVF5R3RYVAfS2K49gDTXOs6+IoGpQNuDb3\nQFIQYZ3eems+PpgDwJweq+dp8DGx/B5WYM4Ci/sMTt5tBj2CFihrUQIxL/0IhzrUUV69zmhh1HUu\nqMv0cvcJeEVTdygVt9QBkheggdNQlOqF6fer0lJeXeeK+X0N5IcAC2Xa6aTvi1qU8Kwc5dWLchht\np2Ggv6B8ndKrSx12e7mOClMZMm9edT/oZ2oGp1v/Sp2ecEI1WUApLS2idenD06al8dSp2TvZujUn\nAQE+eP7kk0lZNsHy5dUuCY7aqoPzUPrp89avQullw3wRwBeaf1yLFi1atPg9wH9D0g8tWrRo0aJF\nixYtWrRo0aJFixYtWrT43cHXALwFgE+o/wiAXwBYAeDnAPhwzdsBrATwIoCP0+vnD3/GSgD/m16f\nCuDbw68/AUDKiPYp5gB4AMALAJ4H8J+GX5+I93sAgF8BeBrAbwD89+HXJ+K9BoYALAMwXBk1Ye/1\nVQDPIt1r5FVN1Hs9DMBdAH6LtI4vxsS91wmBKwCci1Kh/A8A/3X45z8B8KXhnz+GJKAmA5gLYBVy\nksKvAUQexk8ALBr++T8A+PLwz/8KwLfG9Or7w0wA0fx6OoDlAE7HxL3fSLSchLRZLsfEvVcA+C8A\n/i+A4WT1CXuvryAJVcZEvddvAPh3wz9PAnAoJu69ThjMRalQXgQQtaIzh8dA0v5/Qn93L4BLAMxC\nsiAC/xrAV+hvLh7+eRIAk7y63/B9AAsx8e/3IABPAjgDE/dejwVwH4CrkT2UiXqvrwA4Ql6biPd6\nKADXO2JC3euYHAH8IcfRSDQYhv+PyZsNgDPnVwM4xry+Zvh1DP8fLRp3A3gfVetqf2Aukmf2K0zc\n+x1EstjeQqb6Juq9/k8AfwyAK3Am6r12kJTnUwD+/fBrE/FeT0AS8F8HsBTAVwFMwwS7198HhcLo\nDP+bSJgO4G4A/xmAlopOpPvdi0TxHQtgAZL1zpgo9/oJAG8jxRTq6sQmyr0CwHwkY+gGAP8RibZm\nTJR7nQTgPCRK6jwAWwB8Rv7md/5efx8UyltIriSQ3MWoh12DFNgOHIuk+dcM/6yvx3uiQXRwoH10\nMBtzTEZSJt9EoryAiX2/QLK6fowUmJyI93oZgE8iUUH/AuAapPmdiPcKANGr+h0A30OKDUzEe109\n/C+aydyFpFjWYQLd6++DQrkHwB8N//xHyIL3HiT+cQqSO3oyUrBrHYAPkLjIAQD/BsAPzGd9CoDp\nKLXPMADgDqRskf9Fr0/E+52BnP1yIIDrkCz4iXivn0USJCcg3cP9SNc5Ee/1IADRE3kaUibTc5iY\n97oOiY6KRiwLkWjbH2Li3euEwb8AWAtgJ9Lk/VskDvE++LS8zyJlT7wI4Hp6PdLyVgH4P/T6VADf\nQU7LmzsO99AUlyPRQE8jCddlSNkeE/F+z0LinZ9GSjH94+HXJ+K9Mq5EzvKaiPd6AtKcPo2U+n77\n8OsT8V4B4BwkD+UZAN9F8iAm6r22aNGiRYsWLVq0aNGiRYsWLVq0aNGiRYsWLVq0aNGiRYsWLVq0\naNGiRYsWLVq0aNGiRYsWLVq0aNGiRYsWLVqML/4fl6yKp5t1qcAAAAAASUVORK5CYII=\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+/AAAIABJREFUeJztfWmwXdWZ3XoghCaQQExGAgmBwMyYwRaDbdntsdttenBX\n8ifVlVTnT1KVVKXKM22c7q6K2z86Dl2VSnXKSbnTbtsxYJAASUxiFGhCEohJDGKQ0MCgWWh4eic/\n9vty1l5nnXvPfYOEn8+qevXuOe++e885e+9vWN+wgRYtWrRo0aJFixYtWrRo0aJFixYtWrRo0aJF\nixYtWrRo0aJFixYtWrRo0aJFixYtWrT4rcL/ArANwHN07lQADwDYAOB+ANPob98F8AqAlwB8ic5f\nM/gZrwD4b3T+RAC/Gjz/NIBZI3v5PeEcAEsBPA9gPYD/MHh+LN7vBADLAawF8AKA/zJ4fizea+B4\nAGsALBw8Hqv3+gaAZ5HudcXgubF6r9MA3A7gRaR5/CmM3XsdE/g0gE8gVyg/BvCtwdffBvCjwdeX\nIAmoEwDMBvAqgL7Bv60A8MnB1/cB+Mrg638H4L8Pvv4XAH45olffG84CcNXg6ykAXgZwMcbu/U4a\n/D0OabHchLF7rwDwnwD8HMCCweOxeq8bkYQqY6ze688A/JvB1+MATMXYvdcxg9nIFcpLAM4cfH3W\n4DGQtP+36X2LAcwD8DEkCyLwLwH8D3rPpwZfjwPw7khd9AjgLgBfwNi/30kAVgK4FGP3XmcCeBDA\n51B6KGP1XjcCmC7nxuK9TgXwujk/pu71uKP9hccAZyLRYBj8HYN3NoBN9L5NAGaY85sHz2Pw99uD\nr/sB7ELVujoWmI3kmS3H2L3f45Astm0oqb6xeq//FcA3AQzQubF6rwWS8lwF4N8OnhuL93oekoD/\n3wCeAfA/AUzGGLvX3wWFwigGf8YSpgC4A8B/BLBH/jaW7ncAieKbCeAzSNY7Y6zc69cAbEeKKfTV\nvGes3CsA3IhkDH0VwL9Hoq0ZY+VexwG4GomSuhrAPgDfkff81t/r74JC2YbkSgLJXdw++HozUmA7\nMBNJ828efK3n43/OHXwdHOgHI3/JjXECkjL5P0iUFzC27xdIVte9SIHJsXivNwD4OhIV9AsAn0ca\n37F4rwCwZfD3uwB+gxQbGIv3umnwZ+Xg8e1IimUrxtC9/i4olAUA/nzw9Z+jFLwLkPjH8Uju6Fyk\nYNdWALuRuMg+AP8KwN3ms74B4KFRvvZO6APwU6RskZ/Q+bF4v6ehzH6ZCOCLSBb8WLzX7yEJkvOQ\n7uFhpOsci/c6CcBJg68nI2UyPYexea9bkeioCwePv4BE2y7E2LvXMYNfAHgHwCGkwfvXSBzig/Bp\ned9Dyp54CcCX6Xyk5b0K4DY6fyKA/4syLW/2KNxDU9yERAOtRRKua5CyPcbi/V6OxDuvRUox/ebg\n+bF4r4zPoszyGov3eh7SmK5FSn3/7uD5sXivAHAlkoeyDsCdSB7EWL3XFi1atGjRokWLFi1atGjR\nokWLFi1atGjRokWLFi1atGjRokWLFi1atGjx24Wbb741Kjfbn/an/Wl/2p/2J35+CIOOhY1XXfVD\nFEWR/axaVeDuu4vK+V/+ssALL1TP//VfFzh4sHr+Bz+oniuKAn/5l9VzAwMFbrmlev7ddwv83d9V\nzz/1VIEFC6rnf/Yzf43f/naBW2+99f9f14ED6fw//VOB9evT6w0bCvz0p+n14cP59Xzve+Xrv/mb\nAnv3ptcLFqRrKYoCW7YUuO228n3f/3719aJFBR5/PH8OP/95gRdfzJ/ZbbcV2L49vb711vT7hz8s\ncORIgX37CvzoR/nfqq/Tve7bV+Bv/7b6PBYsSOPcdMz4s+Nn+/YCf//31fO/+EX5TPnnllvSOOt5\nfk7x8/rrBf7hH/xnHD6cj0vc68BAge9+t/zb++8X+PGPy+MVKwrceWc+V+K5F0WBv/qr9LzimD9L\n58D+/en9cfyP/5jm3fe/X87l5csLLFxYro9bby2wY0eaz//8zwVeeql8rvHc3fNfs6a87rjXujGp\nO3/XXQWeeabZe3fu9GvuJz9Jz7TJnFmzpsAdd/jx4+MHHijwyCP+PXGvt91WYOvW8u/79qU1qHOI\n59a6dQV+9av68Yt5um5dedzfn1/fwECSG9/5ToFvfSt9xw9+kH5/85tpftxySxqblSvzzz50yD9b\nJ+MWLCjwF39xa3ZuYMA/V3du0aICy5Y1e6+7pgceKGUS/wxJoTh8+CEwcWL1/PvvA6eaNmSHDwPj\nxzf77IEBoM90L3r1VeCCC6rnH30U+Oxnq+cfegj4/Oer5196Cbj44vzca68Bc+ak10eOAP39wIkn\npuPnnwcuuSS9XrgQ+MM/TK+fegq4/vr0et8+YNKk8vP27wcmT06vV64Err02vX7kEWD+/PT6zTeB\ncwcbJPT3A8cfn16vXQtcdVV+7vXX0/UdPgyMG5fOvf8+MH06cPBg+WyLAjjuuHQP8b9NkOZGjgMH\ngAkTmn+GQ1y34pVXgLlz/f/o2G/eDJx9dvV9DzwAfPGL+bmiSPcez2j3buCkk8q/r1uXnm1gyRLg\nK18pjx98EPjCF8rjl18GLrqoPP7ww3KcDx3K5/TGjcDs2eXxo4+WYw2keXfRReka33gDOO884IUX\n0lzk9cFzpylOPx14t6anrBtbB55v3fDee8Bpp1XPf/ABcMopzb5//Xrgssu6f9fq1cA115THTsZs\n3w6ceWZ5vHw58KlPoQKeW0uX5uOzezcwZUr+/ueeAy6/vDxetQq47rryeMOG9D8nnwxMm5ZeT5yY\nxm/q1LR+xo1Lc+7RR/PPPuGE9MwV06en58u4/HJg27b8nJORQPq+w4fzcx/7GLB1q39/E4wb56+1\nDiOmUD74wCuUXrBxoxdCdZNk3Trgyiur593C3LvXL9bbbwe+8Y30+vHHgU8Ptqbbsye9PwbvvffS\n4gWAhx8uFdZDDwG/93vle6ZTI24WcM8/D1x6afV7nnsOuOKK9HrfvjQxX321FLohcDZvBmYM9hQN\n5bFzZ3UR96JQ6ibmhx82Vyis1Bh1CsUZGAcOlEqc8eSTwE03Vc+/+WYuwIGqAFiyBPjyl/PjL30p\nfz8Ltb17SwUUyjuez9tvAzOpexKPHwAsWgR89avl8bJlpcFx5EgaqzfeSM9j1apkZLzxRvUewjhp\nqgiAJNydQpk4MY1jE/A87YY6hQJU55MK+8CrrwLnn5+f27YNOOOM/NzevbmgX74cmDevPC6K6rN6\n+uny2QPAli3AWWfl73n33fy71ACIz+T7USN10aL0zE44Ic3dCRNKhXLSSUnJnHJKWgc7dlSfwfjx\nae0wrr02KVHGrFnArl3V/3c45xxg06b83FlnpWfQFPo8nZLqhBFTKP396eE2wY4d6YErnn22FK6M\nDRuACy+snh8YqApPNzEB4L77gN///fxcUZSKcP78+ZlyuO8+4GtfS6/feSdp+sDBg6XAXbmytFzY\nytWJHEoASBbrxz+eXrPwCTgL7q230uRiuOfYRKHM59VjcOCAH2OHrVurCxaoWu2dUDfu69eXSjgQ\nAlqxeHHucYS3F/e6a1f5rMIiD4GxZUsu+JYtA264oTxetCifO6pQNm0qFU5RpM8PAR0K5LHHkjcd\nXi/P3VDK+/aVRk9cGy/wCROqSuLEE5PHBOTjetJJyShqgl49lOm6e0kN6pgFt25Xry69ecAri2ee\nAT7xifR6/vz5dk3wMwSqysLdK69hwK+/MPYC27enc7t3p999fckYmDw5/UycmObp6tVeKF9zTbof\nxlVXAWvW5Of6+oBZs+ZDcfzxVc9h1qxkqDBOO63q9dRh/PhyLgXqvKk6DEmhMMUzFNRNtKABFCyM\nA+++6y0lpTICIWAYa9aULvVNN81HUZRKcd26UsgtWAB8/evp9YYNpfdQFPniYCv5oYdK5aIKCciF\n2dln5+58UCSMt94qabLAjh1D81BY8AyX8orrV7DS5c91nogKk4ATPGvWAFdfXX0ve2uhdPr60r2q\nkaFW7KJFwB/8QXmsVKp6E0x5qYH1yiv52IU3G1Y5ewNB8YbVr541KwsgKS21QIFyLqlC2bu3+l6H\nXjyU99+v91AUdevcYe3anG3YvDlZ3AyeU/Pnz8eTTwI33lj+nWnhANPWQBLiTKMB1ftXSnX79pKd\nAKpsR39/Ut7jx5fU16FDaVz37k3fp57HvHmJOmdMmZKUk2L27PmVc+edV1Ues2Yl751x/PFpnjXB\n5MnV7x91ymv//ubWa50AqePRexFCankEXn656s18+GH6DHXJFy4svRCmsYJ6iPezJXTvvaXwef75\n0pJRN/nFF0svhK+1zsLjiX7oULpeVjJKuwC5EI3vZ4XS19d5Mo0E5fXOO16hOGhcotNn7NqV+GnF\n0qXA52QnFF3wSpGqh/Hww/lnvPZaTsHwPNT5xzE3IMXG+LMWLcoNGqbSiiI983i+W7YkQ+Pdd9P1\nK+U1aVJab/G/55yT5kETTJnSm4cyHMqrLvbJscKAxrYCaqiuWAF88pPlsfNM2dMHvLGhxuijjwKf\n+Ux57OJ0qjSVPn3oobTuWeH196dnMGVKMkrfe68cx898JnmojOnTEzvSBKefnuY448ILk3HLqDM4\nmqJOoRwTystZuu+8U/L+jDp+3WHt2tLNZTiaRKmMwAMP5Pw5kCbowYPlJGardPHikhPfuTMF2QLM\nvy5eXH4uu8mqXNhKeuKJMi6wdWtJtaxZU/WgWPEGDcXPuRvldfzx6bhX9OKhOO+rDhxH6gb1IgI7\nd1bv+b77cg9D+W6lTJnS0xjQG2/kNIp6K6qcnnoq5/U5bhAGQQit115LFnso1ogFsUJhyisUyrRp\nScE6jhzwa2+0KC+nEDQmEXCezwsvNJsDGuN68cXc0wCqyuKJJ3KPhddXQOePJmM4g1kV18qV6dnO\nnZvLt/7+pFAOHCgNuTPOSMe7d/v71LFzyuPyy9PzYMydW1Uo48f3Jvz1u51COSqUV1MPhbllhgvM\n1gXNly/PLZXAwEDVYqkL3rvz7DUcOpQWVCyq1atLS4c9Es1k+eCD0ttgN3nDhqolHsplzZpSQT71\nVMnV79mTLHK2iJ0nx14JU17x+axQxo3rrlCcMHIUY13wPTLOmrzXUSB13pDy2kC9d+wUBmdkcWxP\nFcZjj+UWq8balFZhb0ozy3QO339/MjgifhLxlBCOcS2hUPbvz630UCgReJ8xo7kFOlzKqy5WBVQN\ntrp17uDiEy5wrjFZVdxOsfF6BKosxsGDVbZD544aEEyfAiXNvWFDmruXXlo+pyNH0jWFVzluXFrr\nK1eme9H4xJw5KdbIcPSYUyinnuqD/U3h6LVjQnnVKRTn8m7a5D0UB6aPGJxdFVCKI+DokEOH0kPR\nxcFWbCx8oBRw7F2ERcXWqVo67CY/+GAZ3FeXmhcvx2kCbA25TCC2upy1rgolJoOjv+ooL6dkXLwm\n3qvPNlJjFS5x49lnfaae85JU+ANV5aU0SwjzgCoM9hiBNF4aYI9r1hRxpfCYNgVK4RmxwYiNvfJK\nUqwRF1MPRSmv009P62DChGpmUEDHrFfKSz0UVW6d4Na5C6oD6Z41NqIxNKfMNIaoGV+Odnvuudwb\nUq/XZXM9/XT+uUq9RVZmUaTxnT69vPfwUEIon3deMhLXrUufsWJFfn033JASQBhXXZVYGUYv9JiD\nW/unnFL9zGOiUJxQqEtz5FTXbqjL9HGoi5+4hAHluIEkrMaNKwU7Z/UwX+qUS7jdTKOp5cxUH3+/\nC+adcEJOXbFijb/z8922rbTmnKdXp1AmTEj3PVTUKRSHXijNVauqQVL1KgK62IGqRckeJZAyslhh\nbNpUCrSiSM+QFQaPo3qaS5fmCkMzy9irCgHX11cq3UgyiHGN+RqGQR3l1S1Tx1mbvVBezkPRrKZO\nYCUc+OCD5unFali9/HJOMbn/4xovICntXikxZ8SqfFMjYcmSpBymTUvr+b33SqNPPZSLLkrG1cGD\naQ4+8UT+XRdemO6VMXly+f9DQdSiMU47LRm8DOfhHJMYisPu3XmcIcAByYArgAJ80Lkuk+vZZ/Oa\nA8C7zUBVoAA55x4B+5h4bJE8+GBJY4WQi0nNWSks1NS64uA8F0RyptAzz5QUm8vR37GjDM47Xpih\nMZRQKHV1Cc4YcJ7LcBVKXXLF1q3VGIzLxAG8Au3mYfBzVg9DFcbDD5eeJVANsCt1yrEDzfgLajO8\nj04UUgi+bpRXHVxxYy+UlwvKaw1IJzjGoJcML01TV69A45hA1XjUVO5t26rrRP+HmQTAe0+aebdj\nR/I0wwBlCrW/PxluEaCPeBmQnqWOx3HHNa85cl7GxIlV5fOxj1XrTlwtyqmnNvNQRj2G4lBXJe/A\nRXsKFWR1MRFHs6i1CKRFXBTVxbJiRWlJclFa0CdxHU8/XX4/0y1sfcb7wnJ2KbDxPr4fjqVwdkpY\nizzR2N1nD8WhLobiFEod5eXQi0JxAkYDmwG3oJ58Mq8DAbyAUEpq794qJcVp6JzKDVQLElWQcYA9\n4iXxbPfsyQWuBpsj2BuCLpQXe0QKrcBXyqsOzoMZLuVVF9N0cOvRKRQ3L9z4ayadygF3veolKYuh\nxgRQTSbQbC5NNw8mYfXqNNbjx+dKiD29009PY71nT5nOe+KJVZbgpJOqAfszzqhWt593XjXeMndu\nes4MlzrsquWbKpRRp7wceqmSD/6YUaelXe1IXeGiU1RqvQKltRPClGkKprH6+9MiiYXy+OOlQlEL\nmoXaI4+U3op6XWyJcZBx166qBbZ9e3mfrFDqPDyXNsyToZfKaTceriq/LvgOVJWVy/CqC7KzRxZQ\na9J9JnuUQPJEVWFwoL+Th6FBX6VHVDktXpwLo/Daw5ONgHwI2ro5HxQZUCqUSZPKhT55ctXSdR5K\nL5ZlU8qrlwp+V9jqKKa6tHOeP9qCRbsiuHiNZok98UTuwThaVWN/qmBCPhw6lO7v4ovT64GBNI9Z\n0bFXcN11Sc7Mm5eMT8a8eUlhMlzFvAvMu9Th2bOrCmU4Hsoxobw0u6ITNm6sBmyZ12a4uhTuiRWo\nS3tU/hMA7rwT+JM/Sa/DCovJy0KelZFy7Szc1E3mimy+Vp3ALtmArV7O8NJAswprDkiOBOXl4DyU\nuip5B2exrltXTQmvux7HkauHEUI7wM9YFYbWuaiQUuWkwkVjP7wGOCYW3xmGVNyHG38Gx1B4vF3q\ncKd+Xk1Q56GoQuG53Q2OnnQZXjpmLhlDqUDNvHNJIOo1MdsQxxygd89AC4zXrk1e9vjxecw3DCBW\nzEw9XXJJGvcbb6wG4a+7rhqsv/LKtDYY8RmMCy5I84px7rneQ1GF0jQo/5GhvOqEvJtoLiBfVyzl\nJiWn3/L/u+9i3vyee8rCRs0G44wxtax27y6FEQua/fvzxcBZQGplBViwMW3CCqVT/yQg93xGgvJy\n511GmauSr1MIzgp2AXkXkHXZOEDuyQS/HOOnNTpaz6S1SaowOHUcyNOjwyKO71LlFMqIuf94fxTe\navoy32eACxsDM2dWixu7xVi6waXgO8rLzcNevJaIJzG4lx1QHScnB5T21fiJo9aUTtT/UbpTa9ri\nOp55JilApuWCNWCZF0J86tT0LI8c8eM5aVJ1XbpzEydW6bLhnKu7Fj33kaG8HOddB9ezqa6q2vG1\nTDMFVq6s1q9oCi8H1jkgGx5JTEC2hHWysnWkzQyjMjr+FlbV22+XXg3HT1hZappkp3gHBxtHi/Jy\n2X2OrqgrcnNwzQN1oQMpsUFb8qhwUxpSW/BwajhQLSRlSk+VkypTTXVWbybmVcwHFjSh/N98MykU\nZ3i5SvmAq5bvJaOrDjq/HOXlFIqjJ3uB1oYsX57TknUZX4xutSTOu1avR5mPFSvyjMKYX9F7j+Oo\n8QxYMZ91VvLgL7ooXV8E1V0g3WVmjTacPHFJAh8ZheJSCevgJq8LyNdlcrlmhmpxAsCvfw382Z+l\n1xG3iAfLKYUcMAfyCblkSalcVBhwFpfSVCyw+H3c8I7/x/WxYnDWEC+OXhVKLxamwlXJa0A1rq9p\nMaxr16HxCqBaT6L9lzQTkONUmlShHraOvyqj++/Pv4uVE3tT4YG5Ku9QKOz58CLv6/PCZ8aMtLYY\nfC8jBUd5OYXiihpdVlZdAa1CqUBNFXfZoEB+/+r1aIDeXZ/KIG3P88ADaQ4y7RceOsu/uI7wCiI1\n+KKL0m9Xe3LZZcmYZLi4h+tQPNpwTSg7YUQUiks7dMVOvWSDaX8lwNefuF48ReGVDKeo3n03cPPN\n6fXhw2kihCBmj0Q9LeZV1U3mmM+jj5YTUt129sjUwuskGPhzOFVbFUqMRbcYSi+Ul1v8rkr+9der\n4+asTJd1U/f9Llaj/DYLQKUrNJFDq/Bd8JUVhm6ToD3U2Gtm+ja83FAsbATE+DENpM/YWYx1nPZw\nDAOHOspLx9ut82gxw2CvnP/XxU4ZmoK+bFkeP3FJOkrhsdEGVAugXaKLphhHfLOvryxniMLjukQZ\noMzOisD89ddXm0I6JeMC81Ecyzj11GqNCTBy86FXQ2VEFIr7UuehuFYivWYKqaXHfbcCLiajAoit\nRhb+QG6tcKFcZH7FdbGbrJOKuxJr3Keb98GT4eDBkmpi5cMxAlYo7DV1i6E0+f5u73NV8hobcOPm\nMvicx+M6yGpNkBoVrpswzxHepiCujylXVtbqzXzwQU5/6cZdoYxYIIfX5oLHrtbKPX/XumM04Tpi\nuE7Dbp27BIz166u1Yy6xwQlmlgM6lzR7a/t237yS1xvXfQHVfl5a6xU0d2zQFx5QzPVOtF+swQsu\nSM/FdTuYObPqdV5xRTUwf9llPtNLA/PDTdIYDkZEodS5rmrNOIXi6IA6isQJMFdEp9XSQJ7d9f77\n+cTlDqSa1syLQwUVKx6togbKhcBxAe5dxVld7L1FF1ogb+PNljrTcPy8jibl5RCdkhnO23RJCm5D\nLVcdry1Y7r23czdh5dh53NRI0MCxXuf99+fejNJfIVx0rvT1+Uw1/j6lvBiu2nk04daa1t4AnoZ2\na9IZFWpoadzTfZ9el1bZa6udOg+GFYwaoFrTFvG4yPwMOjXinO46FTyeLvVb4QLpoZQYTVOHOfV8\nNDEiCqUOuihcbyrnTbhK6abWK5AGmGsIiiK3Olm59PfnxY+81a82LmRXWa0p5tHVSudsKC58ZK6e\nM8nqUoaZfmPaZagxlJHm3evgMrxc8N5l8Ln+XRq4V8uf6U71ZnQeqcLS6ni1XjW9mDl1rl2JGCCn\nm7vizk6UF6Nb+5WjNZb6PW5sXepvND9lqPGh1JRSyo6G66Yc1NDTGJ1SlkDV8I1uF0HbRYo4z7M6\nmjgScwYG0nM6dMi3YYkgfie4vU3cplru3HC3Am6KUVUoCpcp5IK4OpEAHz/RfHTAZ4etWJEH+F9/\nvfQ6uGARqG71y9QITyAVNDwpuX+X26I0BBgrThamr7xSWtS8F4t6KI7y6rUOxeFoCSeF2xxMW19o\nBp5SM1pQp+OrwXwVOCr02Thx9BcbFZw8ENfFhkJ8FgeFXRzKgVODTz65ui1sCKvfFjjlqWOtsS5N\n1HEdh1XBaAaprlv1GjX+VrclRjcURRqnPXuSAnvrrWSIrF/fvFHk2WcnA6gTXEdxVy3vgvyjQaP2\npFBGmh6Jz1TX2qX5OZfZWa+8u2KAhYi6wOx1qPvOhY6bNuU8PVtCavVwPy4t5OQUVBYmzNtyoJNf\ns4eilFedh8IxFNccUsfUpbE6T9DFvly/LkdfOkvTCQfXQSB47IAWr7rmkOzNbNyYGzBsJOhe8jpX\nVq3KBZzSX1G7ws9BC/eAMsML8Ba+W2fcfsWlDjsP5lhkBTnUdQ9mL97ds2ZeqaGpAXpHs3PqPlBN\nGtEEDPVYo6txxOk69WQLxFa6UYsSGV4RmNedOAEfM7nmmjR/GFOnJhnSCaHIGM5DOeWU4bXAd+hJ\nobg+RK7wqGmKYK/Q73GV9Bog6+9PP7HA77ijpLuOHEnXH/fEymhgIF0zZ36FUlI3mS1TLcJjC7nb\nM+A9F+KzmRZgy5YpryZ1KNHt1n0fw9EVTavkXdDZ1ZC4TdPchlpuOwJtp6ICQDd144I2VYLa5UAL\nZDXhQ9u/cHsPnhMcCwhjg8c+NtZiqNJ2MZROCsX1+xqJ+pRe4YwFl/WliraJglEPRlvpaPxEqWcn\nq5R61XqUmIMR34sWT3WF10DZMTiEeNSicMq3epnO07jiimS4MlyKcbedWYH69ivHVKG49vCuHYOr\nS3EZIi5g5oJ8LjXY5aO7DCOtYeCJqxYOL3TNBGFhqW4y03YqVDnIz9YxtyrvRdnGJOa02yYxlKZw\n6dZOobgqeUdfuh5erkJehQNQrVYH8nmk7VRUmLmCNxYWSn+pMNHtpDk1Wekv7nLA2yEA6T2cIuvm\nqQpLBaeH1imU4XQcHgqcV+mUh4uNKRWl3odmijnPQGMwui/90qW5AaAFqW5zN5VxcRzXG8xEpyat\nEQAPIT59elL2rIA+/enqtsDaPNJtOeF6es2Y0Z0ac7273LnhoieF4uiLpkWNrimkFiABvqDR9e/S\nVFAgry0JcAaX5syzxanFTlzMppYtu8nhUsdkWbq0eq3xNxY0nDbLwrmT5aOI9zWJodRBlZlbZE6h\nuCQJV4PiMn5cgdzBg9W5pXy41jF1q1hXb0a7GL/9dk5Vsnfm6C9OJVX6a8mScr5ErITnDXtqLm7S\nzcDg+hPXm8kpFNdx2HmqTeDmZdOUYd7CIaDxVKW01WPV9keOXtX5ogaObg2tlLkas5s3l3M8vM8o\nP3AJRgH2ULZsyZ9bFKtee22VznK1J/HdATf2LtNL4WqajrlCcbnprjGkK3ZyKcNuXxOX9eVoE0cb\nKA2zf3+adDHJ7rwT+NM/Ta8HBtKkjgWv+5LzXi7Kw7ObrNYXC1DNZ2eLl61vtta48lgnQJ2i4R5i\nvWwB3JTycp2GnYfiOOxuNTeAj9G4ink1IhYvzivYddMlbT/PlK0qTv0+9V4d/cXCKTwnvhc2mGLX\nxoA+e227Jge2AAAYDUlEQVQFwvSnwmX7uH5ejvLqtXtswPVyc+vc7eej1LS7J40lcR0XUE3AUY9G\njQ3X/01lFXesAKrpwkFzsycWCSHhZTrPKTyUadOq8Y5PfCJld7oC1XnzqkWP55yT16i4Nes26nLf\nrXAKZajzIzBsheKq353l4jS6E1SuoSOQP0gXEObajcDChXmAnntP6RaiLNSVTuPOwyqI2OrhNEH9\nv/h7TD5WSkwLdWoKWUeN8eeOBOXVxENxygPo7l1pSjdQpReBarAUqFKNTLdq/y1t2KjBfPV6lf7S\nflAq4NjgYKqNe7NxnGDHDq9s43rZQxkKnHBwlNcJJwxNYDQtauxG3QHVRJW6fmY8l7RvnCoDjZ/o\n5mnOo1FZo0ZvGIcuvheJFY7yj669fP2hZNgz0TE75ZSqEnCBeSCXBS7l2KUOK5p2HO4FI+KhqEJx\ncRCXzaJwaaPOCnL7nNx1F/BHf5Sf49oQTSflgrTYtTEmwD33lLUoQD7x1E1mxaCUD3tbnXbA40QC\nntQc73EdYR1GgvJqEkNx2XkK116FBW7Abail2VjhXdRRUFq7pN2E1fpUgbRpU2f6i1NHea8aIM80\nY2+WY391e8yEMBpqTC3gvBZHefXajjzgaEoX21Q4Qa6b0L34Yh5PaaJcta5Fm31qcat6nCpXtH0U\nbzcR6ctspERsxVX3u218585Na5sLVD/72WS4KHj862Im3bwWlzqsmDq1mn5+1BWKCgjHBzdRHs4q\ncdSW2xPeKRTlw8NziofNdFdR5G64cu/cd0h5WBZEStXotTLdwy56p2aJ7Glxx+Gm9QqjQXk12a3R\nCUE3ni7QrlSEU5465upRuIaNnEmme6PwHO1Gf2n2l3Yy5noJfVa8YRaQK1kuaqyjvDp1ou2meJyH\n0oTScDSOUyiuK4LC9XHTeInSV5oBpvE6JzvUGNWmlUqZacNPHWNu6RJrwrWPcW1XXFV6pA4zgv5i\nXHBBkjkBl2LslIzCKZTjj8/nkjNCjqpCcYJQ951wcBNfXVKgugkOUO3BBVQnsqNfbr8d+MY3ymNe\nvFowxamn2kZBeXv2VpQ2Yy9IM2B40XAw2QU745gFW9PNrEaD8qqjIRku68VleKnyUGoD8H2+lIJS\neoIpKM3A0oQLVXSaSNEt+4vnJG8NzIKYlRR7sex18pyss8p1b4xYS00CqkP1UNy1dNuTB/AJHS7D\nS70BVeDLl+cKRpWBNut0BonOKZ3Xmo2mW2A89lhSKGw0uphvnYeiQpnbpkQtiRPorsARyGWoSx3W\nIsXp06utek4/PXnXnXDMKa8m3ojbWMcF37ngK6ACx6VcumJGpkw0fZTpD93qV6kS5m7ffjt3k2PS\nxXUCeUYXLwKmvLg/lEsuCLC1yxRPp2D9cBWKG+MmcNlcbndOhTMslK4Ij5LrSdigUAtWDQa1RrW9\nitJf6u2yN6uV02yBr19fzmluGsmxH56/qlDCQ+Hx5T1R2ONwqcOKoQbl6/ZC6eYhuyw/DdI7A0LP\ncZskoKqUtCEkr0Ogmv2ncV7ngWm6fBitTNEGe8Dj4zwUpryiuJQ9DY6jaDJFUGMMHWtXuDhnTr7n\nvGMfmrRfOeYKpQlchpfm+Af4QegmU4BPF9b95NXi+c1vgD/+4/S6KPL70IaETMkoDaGV0bzwNF2S\n6yrUYmJvgxeLm+jxf/w/nTjmXmMoCuehKFxShFMomuHl6hYef7xKX+oi3bAhp000AKsehRYgqkWr\ntSyd6C9N0FCvlFvfs2DjWhtWKGw0KeXFYxrjzgqFhU8ThXLiidVK+SYeiov3NYlruJRhNTi1UajG\n2RybofNIDdRu7em1vY5SapweDOSGkMb3+vry7EEXQ2bK68wzkzHIuPrqdM1A8nwfeST/fH0GdenE\njCapw664UfGRUyguEKsCH6hONFdh6+pPdA+FXbuqQkrpLrb016zJg7c8+UKg8zbAPBHZAlNqTK+V\naTlnhQeYT3YeWqCuMSSQT8BeYigOTRRKXZV83bUHXDaX/p+jTVRh8GZoQNUyZkUfYxoCWila7Xbt\nih/ZgFF6jIULjxFTgJw+zs+OaaROlFcolG7tV7Ri2lmpTTwUtxb189w4uZRhhQbkNUlDC4ObKhyW\nJZqRpgaFKhxNCWejNeg0/l7OWO0WlNe6kaLIPU1XDa/tVVzbejUW587trlA+kh5K3aZIAZdK6GoK\nFFooBlS9mAMHqpbxvfeW+8IHOtWJsIAoirw1h7bx4JiOeg9qJTN1pZktTKmowGbL1HlyAbbuO8Wt\nRiOGonA1KBpncULJtawHckGldAZQ7cjL9ST6rDdvzmlJbYmuXqZmf2mVu1KSvAY4puioHHePdW11\n6jYcUw8lFMoZZ1T58Ca9mZrGULq1Y3frXKkq19Zdixy1Q4Km6Or4afNXFeiRkcfQJA/NEFODJAL6\nde1y2AhyDW85NZsVikvvdQWH11+fnkNg/PiqEaDU2NSp6fkzNDngzDM/YgqlyVaubktQ9UacQFTe\nG6hubfr449W941VhcH8lIFXPR3ylKHJLUDOBtKaEhbNeHy8EnRD6Xs4Uc7GjACsUZwGylV2XdXU0\nYiiuSl7h9v5Qz8rVDmk8Q+9VU0w1s06bQ2p1tD5/9hKKIg8Ya6xGqRH2SjmWosq0SecD/q66GApT\nXk4QuX5eiqYeSjeF4ta5gtvRBJwsYMpIFY52ztDu0Vp0vHp1vvaUitUEDc3U5HXHXb+1GWw3bzzA\nXgFnerFwVyqqSe3JUDK9XMaYerUfOcrL1Y0odE8JIGlXthqcV6PZHmHV8YK9++6yHqUocmvkuedy\n11fjMRyYU8GlbjJbJjzxgKqVC5TXyLSPpj+yJd5pe1RXEBqo2wI40C3V1Hko+j/OQ1E4YaJwG2qp\nEHfVy3ysHLcW0LIBofSXChdNcVUvVKkR9mi5xQ9nqe3YUe2YHOiUXOFiKN0URpOd+prGULrFS5qu\nc92ptJty1fomVW4qzJVG1Tqx2As+oLt16lplA4VT1ZmabuLBBTizihVKdB4GqnEU1xl81qwUUw64\nNGRFk1qUk0/OPZujqlA0OOaCyEo5OAGmPbyiwpzh4ifa7ViFzcBAPtjKz95zT269alCd3d7Fi3Nq\nhNNStahOeXVWYiq0OGvMxZZcyrCiE+VVtwUw0KxKWikkl/OvGT/Om9JC0p07qxssqcCJrVY7vYfr\ncboF1FXoqefYrfixW6yGv5trlzggzzGaTttd10GzvJTWYLj2K4qhZnkpdJ27e9OkGp3vuv2xjl8T\nBaQej3ZI0C7GykpoOvqyZaWHzHO4Cd3vwGuQ9ze58sqyXf2llyaFxVBjUL0WR4PpTpCzZ3dXKNpx\n+KgqFIVrO6Dpd1pVDFQFsrPGXRsETUFVTl5TdVWBsFLQjBSNxTAVogJc9+DgNFNVspoVxFy7y9EP\n6EJgDJXyarLJllqIroeTvker2t173C6cWoym1qNWqGtxpzYMVPpL+7Pp52u9i85VNmCUGmEDRI0h\nLqLkmhdWOp1QR3m5eiWmK0bSQ2GF4gxHzUpzWxfoNWuXaY2P6BxR2lS3G9B1qcaQXrfbnZFTwtVA\n0esfLvizmH5y33HFFTmldeml1doThe4v3yRmovVMrjamFwxLobg+XopOgeaA6zCsQU6lp8Ii4vdw\nCmekB8ek18m5YEHeXoWVnFo5Koi475d6V5oDzymKasVxgZwKLG5roYKXg58acORjtXJ62bUx0KRK\n3qUMK5SacD29VMEuW5ZTYi4bRwPqHH/RgDoLQS1+1JRm7Y6t9Bx7M2+8UX//LASZrmm6Nw4rFIUq\nkJGMobCy6ESxBlzKsEID6touRQubleLWeIkea7W7MhT6/erBcjYmb63BntNIbzIYnzdjRrqewPXX\n5wWOziPRoLumDrs4G38nMPIdh4elUFynYYUqFOfGqiBxjSR14yTdyjO0fViUqqRUgbAVqZbOgw/m\nVAhbslrAqLy71lVw3ElbN3CWj94zx2g0756Vhot5xLXpcx6KQmkiTFShuOwe9XQ0m8d1h1Xa8/nn\n8wp3Duiqdemyvzigrtaytt9R5cUxEiCPI/HfNBYIlGPFCsUpVIdOCkVTh13diaKJh6Kx0iZFjapQ\nlM4CqoaRshk615S50AafaqSoAtK4pxa4uvhcjDkbEBwPZLkxFHBF++zZZQPHz30uXW/g7LN9zQgr\nA6XKzj+/rMavgxa7fuQUSjcPRXlUR4+om6nUhePo1RrVdFC1Xnk7WW1eybsxArnlpG6ya9lQV5eg\nLndd2izQ2ZNTipChxZedMFoeinqqLsNLlZumi7paHQ7Q6/4kKuQc/dWpnYoaJCqctm7NFRArf1VW\nTH+pVcwCgFNtXfcI9z/dFApbtU3QtD05j5frNKxQWaDjoZ60Ws4ujqpQ1kJjqk5BsRGj7YHUEOJy\nA26zwhlebPg5ueTA1OT555e9ujgw3yTIPmdO3kE49qgPuCwuhdbFuI7Dw8GwKS9e1NyhM6CBfE3b\n1L4+QBKuHPzUwFm0S+HP5WyNI0fSZIvB1n5g2qqFKazY+jcm/0sv5R6IVmGzFaULQntZcdW1CqVu\nCqWuj9fRVihujIF8oet9uyQCXfwq4PV5PPlkTm80ob+Y/lBDhuNpajRoTYgmSGj2F1Nn7PkodcRz\nvZNCic8EOu8LP3Nm92p5pT2H0m1Y+3jVzQGGsg5qMGzcWK39YG9E9xJSA0KTXdQr1hIHXW96DzpW\nbORy3Qm/dm1XHLh+iBXHxRenGBvg4yhnn513FdbA/Lnndg+66/zR4saR3ld+RD0UV1eg0I2GdOI5\n6kNrE1TB7NmTBEAIBO1MqwqEuXVti677naubzAF4VSAqBJXT5ftS953pQ9eq5aPioajl7uBSO9kz\ncxtq6X4X995bbdDIlJPOo070l6Yia6W/Buc16ULjd6ystIU70zYvv1ztthzoplACKmhYQbCgqoN2\nHB7KBkqqUJqkjWuCg1KM2mFY17gm2OiafuKJfD5oQaw75v/X+BwbNK4639GWTRgaIBfiHOfQLEzd\nw+SGG/INty65JM37QJMYidKiWvMyYUJ3mrQXDEuhqOWp1bPuZnX3No11uACnVqNqZ1CuPQGq1u7O\nneXAa8BdA3u6sVMnCku71upncVYPV+QDnTO8dDF+lCgvFfx1E1ozvFhJa+qmUzD83IGcclL6Sz1l\n106F54sqCDUaVNhps0imWrTugaFKj9FUoSimTy8piiYZSNpxeCQ8FC1qdOOn0AatyhqoV6u1ajpn\nnIJiClUVlMZXdNw4eYYzM5WK45IA13bFgYX4lCl5IJ23J9A4irZcaTJ2mpihqcNN2q8MB8NSKEpn\nadZEk5brOjm1pkOzbQYG0o/uRx3Cub8/TYAQ3lpfopXUGmzlWIu6yeq9KE/fqbuuLphOfY80ZbhT\nXcBoKxQ1ALRKvonHoi05um3nqty3CqNOGVeAL37k71NPkse8KPJ5re1+tPcYKy81sHg+dCpiVXSK\nJXTzSlQAahB2KB6KzjFd5zpfmxgZrqaNj5UK1zirGplqsLljXsvamYBjPNw+f+vWem+sKeXl9oEP\nMAXG8RWgvhcfP19tSa+ZXur1aJ+wkcawFIpCPRQtZGpSgasbWnXbaU+L4VQh3X03cPPN5TF7BtEX\nKiae8rTqcWhbdc08A0qrUbNcNOjLC0rTiV2X5TqMtkJRKN2h4+U6Cis0O1BjEpz+DVQbNKqHqhs2\nafEjexTKneuYa4W/tvtR74YDxUrtsVBUz67bVrl13ke34kW1gJXyGuqOjQy3ztno03iIevbq0ej8\n12C3yg1lGTQpRGWCGkHaQVoLrXls3B4ogU4eCmdzdVIoHJjnfm+BCRPyNTtnTkomCOjeKKpQmuzu\n2A1a79QJI6pQNFCm3oVmfujEiMWp8QOeDNpK4de/zjsLq9fAGSq61a8KD92JTy1b5sfVClMuXdNi\nO/X/0WyTTgpFBWIvWwqMhEJRAazV45rhpYvOWa/qgWniRif6K4R2XZ8zDQard6Njrh0StFqehaUa\nDW43ykCn8XfdCOqglIbu6Kh/V8prKB6KQtetKhSlczWxReWAZsZpRwM1IjU+oh2iH344T5xx++Gw\nUbBkSfl3HQum3nTudvJQuOOwtlI59dRy8yv1SvT4uuvyQLy2steeXueck7doGUrHcUUvPQFHVKEA\nuTLQymAVFFopq4JVMzV0614g3wv80KG0wFx2BpAKFDk+ohNLhSF7L0p/OUuWJ7nGDdiS1YyUThSM\nQnP1j7aHotSFZkAptaep0tr3bPfu/FlEBl9dvy01UjSbS70bHWPlzjt5N3E9ocDVktYxZ35dEwE6\nbaSmret7obzOOivfb0M9GKW8RsJDAfJ1rh4BZ00C1Zb1uiOjFjRq12+Ns65cmX+eJlWoQtP1pWPB\nMVYtTGa6Vj3rTt642wY4wF6EegwaR5k3L+88fPHFeWCe96gHhl/p7tDLnBlxhcJQnlRdTRU26l1o\nPYpSCq+9lscstDfTXXflwXqdqJyloY0D1ePQrJOlS6u0Cy8iVkbq4qsyatJNIKA8/dFWKArlstW7\n0mCqFn7qmGtm3P335x5Et+15NZjLtQVALuiDt48xV0WucS71bnSfDYbGfbRtCMPRX3XUxPTpuULR\n1GGtnu+W5TVSld98vdoeXp+jxmC0WFD/X1PMtf6E5UzEkOJ6jhzJWQ/18HU9qUfKcAXXdePEHopC\na07qih2BNB937SqPR8Lj6Ab1ekfFQxmNiadUhwa11XJRa/OOO3K6SykqpkrigYR1rZlUagWpm6yf\nrbEEttxVgWjWkVpQGrRm6GcppXOsFYpCPRitHNdFqbvnaXXzunW5V6s1OeyxqlGgFqT2g1MDo1uz\nSFWGLMiUCtQxV8+aUbcXioMWr2laqCqUblleWnA4GtA5AdQLYk2ZP3gwVwA6n7RIWT1gNUI1fVjH\nnA0QXXtNNpELdPJQtAuwa0XvCj8ZupFap2OtRVKjQuM0Wkw7KgpFXXgNtOlE6PZANBtFLQutJdC9\nTIoi9zAOHEgTL65JvSG1LjXbSwPuygtzixO9V7Vy1G1WT6zJhmMBzQYaDuXVpNvw0QA/O83e4efc\nREGwclKjQFPA1btR5aXUiRoNHLtRKk5pnQ0bmnudvVBeCqdQ2IPpluWlz1/X9XCPNSaxZ09+r5qs\noPEW3d9EFYLuqqoxVq0pUpqSDRZ97tpWqVMcTP+3U9depaU4MA+k73zppfJYmz5ecEEeZ1Gv5txz\n8zkxY0bZ5RiobkvcreNwL3G3xgpFg79qKavg0yJHzfxQS1VbmihloDTR2rV5IE/bp/zmNzndpRNP\naSZuJNlt21gNKqr3om1hXDpxnYWm1qpa5HotTds/dPrOgC4KFQ5KF2iMS4WsWu1aoKlKX/cb18QG\nba+jCkKFh8axtmzJqRZWxrHA4xnpvXUzULTRoe5gydDP7oXyUkyblqeBaot7tTbVQ1GBr6yBBp51\nPuq61iQNrdV65plc8aqXqNSy7nejAXut9dHMS/X+2SjQBqHaTUFjvq67ekA9PaW81FNkqMDXOMqN\nN6bnELjmms6B+W6pw1qLov28VKGMSgxFFYoGqDSVUAW2Do5akxovUbrhnnvy9vILF+bHWsvA1qVu\n9asBdk0rVDdZr015e/VAXKuUOgGhCkGL6FQIqyJ3dMJQofEQjSloxp0uQE1qUKtdhbA2ZNTqePUa\nNR6jCoL3rAgFUddORffz0M9WA0Tno2Z0daoV0iCpGl+9UF4K19K+07HSH+qhaOv6bkWN3TK8dA5o\nQF3lgsbgVIiromZWQ5Wj0mMax9J1q96PJoB0gj5Hpbw0dVjjFIyZM/NU3zlzco/k4x8vW7YA1Qp6\np1CYYtNq+W4KZVQoL6dQ2HLRQJvWoGjmh1ou6rFoqihPrIGBRMHp3t4xsdQ70oJEFRZaCKeWrgob\nVRidrFG18jVnXq1yXVBaJe/2JxkpaOdipddU8WpWnnqR+tw0p1/vfdu2/PN5zqk1qTSRehBqFLgx\n53icdldwHRzYy3SCLaCCTeeyBqG71aTo5w8HqmD0WrXmQ71Mt85Z6OoWvlqnpApCCxhZQbi94DvF\nXFetyr0XXcfa/0378qkCUSXRCfocVSirQlGaShVQX59vQwVUA/OTJuUxEC12PPfcPJVYr0UVitYy\n9UJ5dXOsfwjg1mYf1aJFixYtfkfwn5H0Q4sWLVq0aNGiRYsWLVq0aNGiRYsWLVq0aNGiRYsWLVq0\naNGiRYsWLVq0+K3B/wMOaIK3ghEEcgAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "loc = anhima.loc.locate_region(pos, 11000, 20000)\n", "loc" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 5, "text": [ "slice(118, 305, 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": [ "(10972, 11026)" ] } ], "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": [ "(19990, 20088)" ] } ], "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": "iVBORw0KGgoAAAANSUhEUgAAAbMAAACgCAYAAACG2IIGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGpFJREFUeJzt3XeYVPW9x/H3wIqAqIhiRVmwIRYUI7tSBFFijIgiXGOM\nNdGUaxIevUbFmwgY84hejEqMMRgL9q4RjcSCqxRZUClLU0CWJmVxpUiHnfvH90xmmJ2dPVPOnJkz\nn9fz7LMzZ86c/f3YZb7n174/EBERERERERERERERERERERERERERERGRPDYUmAtUAc8BewNtgPeA\nL4F3gda+lU5ERKQRpcBXWAADeBG4GrgXuMU5diswMuclExERcakN8AVwAFACjAP6AQuAQ5xzDnWe\ni4iI5K2fA5uAtcDTzrFvY14PxT0XERHJK0cD84ADsZbZ68AV1A9etTkul4iIBEyJh9f+HjAF+MZ5\n/hpwJrAa615cDRyGtdrq6dKlS3jWrFkeFk9ERArMLODURC808fCHLgDKgRZYd+K5WEttHDYRBOf7\nG4nePGvWLMLhcOC/hg0b5nsZVFfVVfVVXQuhrkCXhgKOly2zWcBTwKdAHfA5MAbYF3gJ+BlQDVzq\nYRlERKQIeBnMwKbh3xt3rBZrpYmIiGSFl92M4kKfPn38LkLOqK7BVUz1VV3zU8jvAiQRdvpIRURE\nCIVC0EDcUstMREQKnoKZiIgUPAUzEREpeF4Hs+OBGTFfG4Dfosz5IiKSRbmcANIEWAl0A34DrMOm\n7d+KJSO+Le58TQAREZH/yJcJIOcCi4DlwABgrHN8LHBxDsshIiIBk8tgdhnwvPP4EGCN83gN0S1h\nREQkDbfcAhs3+l0K/+Sqm7EZ1sXYGajBMucfEPN6LTaOFkvdjCIiLoTD0LIlfPIJnJowDW8wJOtm\n9DqdVcT5wGdYIANrjTWaOX/48OH/edynT5+CWo0uIpIra9fCtm1QG7ANtSoqKqioqHB1bq5aZi8A\n7xAdJ7sX2xrmHmziR2s0AUREJC2VlVBeDi+/DIMH+10a7/g9AWQfbPLHazHHRgL9sKn5fZ3nIiKS\nhupq+x60llkqctHNuBk4KO6YMueLiGSJgpkygIiIFLzqamjXTsFMREQKWHU1dO2qYCYiIgVMwUzB\nTESkoIXDsHQpnHaagpmIiBSomhpbMN2+vYKZ11oDrwDzgXlAGcqaLyKSFdXVUFoKBx6oYOa1B4F/\nAScApwALsAXS7wHHAR9Qf8G0iIi4EAlmbdoomHlpf6AX8LjzfBe2p5my5ouIZEEkmLVoAXV1sHWr\n3yXyh9fBrAOWj/EJ4HPgUSwjiLLmi4hkQSSYhULF3TrzOpiVAF2Bh53vm0mQg9H5EhGRFFVX2+QP\nKO5g5nU6qxXO13Tn+SvAUCxbvrLmi4hkKNIyg+AFs3zLmv8xcB02c3E40NI5rqz5IiIZCIehVStY\ntQr22w8uvhiuvhoGDvS7ZN7wez+z3wDPYht0LgauBZoCLwE/A6qBS3NQDhGRQFm3Dpo3t0AGwWuZ\npSIXwWwWcEaC48qaLyKSgdguRijuYKYMICIiBUrBLErBTESkQCmYRSmYiYgUKAWzKAUzEZECpWAW\npWAmIlKgFMyicjGbsRrYCOwGdgLdsKz5LwLtiU7NX5+DsoiIBEI4vGf2DyjuzPm5aJmFgT7AaVgg\nA2XNFxHJyLp1sPfe0TVmUNwts1x1M8av2FbWfBGRDCxdumcXI1g2kG3bYMcOX4rkq1y1zN4HPgWu\nd44pa76ISAbix8uguDPn52LMrAewCmiLdS0uiHtdWfNFRFKUKJhBNJgdemiuS+SvXASzVc73GuB1\nbNxsDcqaLyKStupqOO64+seD1DLLp6z5LbGkwpuwTTnfBUZgeRmVNV9EJE39+8PPfw4DBux5/MIL\n4frr6x8PgmRZ892MmX3g8lgihwATgZlAJfAWFtBGAv2wbWH6Os9FRMSlxroZi02ybsYWWMuqLbYu\nLGI/4AiX118CnJrgeC3Kmi8ikpZEa8wiFMzq+wUwBDgc+Czm+CbgIS8LJSIiDfvmG2jWDPbfv/5r\nCmb1PeB8/RYYnZviiIhIY5YuTdwqAwtmc+fmtjz5wM1sxtFAd6A07vynvCiQiIgk19B4Gahllswz\nQEdsEsfumOMKZiIiPlAwq89NMDsd6IwWNouI5IXqajjmmMSvFWswczM1fw62sDldTYEZwDjneRss\nE8iX2DT91hlcW0Sk6CRrmRVr5nw3wawtMA8LPOOcrzdT+BlDnPdHWnbKmC8ikgF1M9bnJgNInwaO\nV7h4bzvgSeBPwE3AhVhuxt5EU1pVAJ0SvFcZQERE4oTDtu3L8uXQOkG/Vl2dTdvftg1KcpGwMIeS\nZQBxU9WKDH72/cDvsIXWEcqYLyKSptpaC1KJAhlAkyb22vr1cNBBuS2bn9x0M36HLZTeBGwH6rCd\noxvTH0sgPIOGW4DKmC8ikoJkXYwRbdrYwupi4qZl1irmcRNsY81yF+/r7pz7Q6A51jp7GpcZ80FZ\n80VE4iXalDNeUMbNcpE1fyaJcy42pDdwMzZmdi+NZ8wHjZmJiNTz5z/DsmXwwAMNn/PDH8INN8AF\nF+SuXLmQ6ZjZoJjHTbB1Z1vTKEckMo0EXgJ+BlQDl6ZxLRGRolRdDR07Jj8nKC2zVLgJZhcSDUS7\nsAB0UYo/5yPnC5QxX0QkbdXV0Ldv8nMUzBK7xutCiIiIO19+CUcfnfycYgxmbmYzHgm8DtQ4X69i\n68dEpMCsX2/Z1g86qP7XOefYGibJX1u32gSQ449Pfp6CWWJPYBk/Dne+xjnHRKTAvPcedOoECxbs\n+TV/PixZAlVVfpdQkpk/33IyNmuW/DwFs8TaYsFrp/P1JHCwh2USEY+MHw/9+9dvlbVtC4MGwSuv\n+F1CSaaqCk45pfHzFMwS+wa4EksYXAJcAazzslAikn3hsAWzH/wg8euDBsGrr+a2TJKaqio4+eTG\nz1MwS+xabPr8amAV8F/OscY0ByqxNWnzgLud48qaL+KDqipo0QKOPTbx6926wcaN1pUl+Wn2bAWz\nhrgJZncCV2HdjW2xQDbcxfu2AWdji6tPcR73RFnzRXyRrFUGltPvkkvUOstnbltmxbgNjJtg1gX4\nNuZ5LdDV5fW3ON+bYd2U32IprsY6x8cCF7u8lohkoLFgBho3y2fr1tlsxiOPbPzc1q1hwwbLoF8s\n3ASzENY1GNEGC0xurz8Ty8f4ITAXZc0XyblNm2D6dDj77OTn9egBq1fD4sW5KZe4V1UFJ50EIRdJ\nCJs2hX33tYBWLNwEs/uAT4A/Anc5j//P5fXrsG7GdsBZWFdjLGXNF8mBCROgvBz22Sf5eU2bwsCB\n6mrMR25nMkYUW+Z8NxlAngI+A/pigWcgNqEjFRuAt7G8jsqaL5JjbroYIwYPhqFD4ZZbvC2TpGb2\nbDj9dPfnB2ESSC6y5rtxEJbLcT3QAvg3MAI4D2XNF8mZcNgS0779NnTu3Pj5u3bBYYfBp59athDJ\nD2VlcN990LOnu/PPOw9uvNH9TUwhSJY13003Y7oOAyZgY2aVWOaQD7Cs+f2wqfl9neci4pEvvoDd\nu+GEE9ydX1ICAwbAa695Wy5xr64O5s61MTO3gtAyS4WbbsZ0VZF41qOy5ovkUKSL0c3EgYjBg+FP\nf7I7e/HfkiUWnFqnsCq32IKZly0zEckD48fD+een9p5zzoF58+Drr70pk6TG7fqyWApmIhIYW7bA\n5MmN738Vr1kzy+H4+uvelEtSM3t2ajMZQcFMJCU7dsDUqTBqlLUAJL989BF07Qr775/6e7WAOn+o\nZdY4L8fMpIBt2WILZxNNKF29GiZNgokTbcbbMcfAqafC6NGwaFHj21NI7qTTxRjx/e/D1VdDTY1l\n1Rf/VFXBHXek9h4FMylK69ZZd9TEiRaoqqpsWnZJgr+QNm0sU8Stt8KZZ0bv+vv1g+eeg2uuyWnR\nJYnx4+HFF9N7b4sWNnHkjTfg+uuzWy5xz+2GnPEUzLLrSGzR9cHYgusxwGgsJdaLQHugGsvKv97j\nskgCU6fCddfBsmUWmHr1gpEjLYN6y5apXWvoULjhBrjqKktaK/766ivLgt+lS/rXGDQI/vEPBTM/\nud2QM16xBTOvP3J2AjcCJwLlwA3ACShzfl6YMAEuvBBGjLA/+n//G37/e+jTJ/VABpb3b9994Z//\nzHpRJQ3pTMmPd/75dsNTTB+K+SadyR9QfJnzvQ5mq7FF0wDfAfOBI1DmfN+NGweXXWYD/IMGJe5O\nTFUoZK2zu+9OPNYmuZVKCquGtGpl0/TffDM7ZZLUpTP5A+CAA+Dbb4vn/2IuO4NKgdOwbCDKnO+j\n55+3rsW33oLevbN77YsusgztH36Y3etKaubMsbHPc7OQnmDwYCUe9lO6wWyvvWzcc9Om7JcpH+Uq\nmLUCXgWGAPH/tMqcn0NjxsDNN8P779u4WLY1aWITQ+6+u/FzxRtr1lj38V/+Yl1Nmerf36b4b9yY\n+bUkdekGMyiuzPm5mM24FxbIngbecI65ypxfbFnz1661u6kDDmj83N27YeFCOOood+Nb27fDgw/C\n3/5mH0zHHJN5eRty+eXwhz/YtP3vfc+7nyP1bd0KF19sk3B+8pPsXHO//eCss6wlf/nl2bmmuFNT\n435DzkQik0A6dMhuuXIlX7LmR64/FsuSH5vl7V4az5xfVFnzd+6EE0+09EEdO1pm7J49bXbhkUfa\nH/S0adGp8598Yn+oa9fa4HDk/J497W58/XqYMiV6/owZtnj22WfT/4+RigcesJ+rRbe5Ew5Hg81z\nz2U28SPek0/auJmSD+fWhAm2vmzSpPTef+651lPSr192y+WXZFnzvQ5mPYGPgdlEuxKHAtOAl4Cj\naHhqflEFszFj4KWX4J13YObM6KLkSZNscsaGDdbV0KuXBazu3W0h65YtFuQi50+danfS69dbN2Ik\nIJaV2UzDXNm82e4GJ05MfX2MpGfYMHj3XfsAbNEiu9eO3N1//XXjG3xK9jz4oO168PDD6b3/0ktt\ngtePfpTdcvnFz2CWiaIJZlu3wrHH2l1v/DhWOAzV1XDIIe66E3ftsvPbt7cuSz+NGGHr1x57zN9y\nFINnn7VlFVOn2t+KF847z9abDR7szfWlvuuusw05f/Wr9N7/y1/aOsN0359v/NrPTFx66CFrOSWa\nkBEK2R2x23VfJSU2HuZ3IAP49a8tUe2KFX6XJD+Fw9mZNj15sm3VMm6cd4EMlKvRD5lM/oDiWjit\nlpnP1q+H446zSRluN08sJDffbKmynnzS75L4b/t2mD492iU8ZYrdpMSOj558MjRtmvw6tbV7jodW\nVVnKqnRzMLq1dq39ra5alf1uTKmvrs6GDFasSG0fs1ijRtnv6777sls2vyRrmSk3o89GjbKpz0EM\nZADDh8Npp9k6pUGD/C6Ntz7+GB5/vP7xcNhSS82YAZ06WdD66U/t3M2bo0Hpr3+1JM7l5XDoofWv\ns2sXzJplefrKyiwA/vGP9jgX41gHH2y/y3fftfWE4q10NuSM16aN7VBdDBTMfLRmjU2VnzHD75J4\np1UreOYZGDDAPqSPOMLvEnlj7lwbS7r99sQfPldcYfVPNAmnY0fLTg82FXvKFMvcEC8UgiFDbIcC\nv7qRIwuoiymY1dZa66ZTp8ZbzQ0Jh+0mBGw8281M09mzM+tihOLqZlQw89Fdd9l6oKOO8rsk3ior\ns/Gza66x/I9BS0JcU2OLlEeNst9nJtq2ze9AMXCgrSHcsSO4W/0sXWot5Uh38LJl1iqtqbFk3JEu\n4W7dGu5u3b3bsrBEWt2TJtmxcNjGtSPX6NkTTjopcZCsqkovJ2OsYgpmXo+ZPQ5cgC2KjtxjuM2Y\nH+gxsyVLbEHxggXFsVfUrl2WOmvwYJusEBTbtlnuwrPPtpuTYtCjh82c9HqMLtdeeAFuu81mF0cC\nTa9eNhuwpMSC2eTJ0eBUVWVrNhPdnH39tXUVxwatjh3ttcWL9wyWa9bA4YfXv8aqVfDII5lNq58z\nx94flK5GP6fm98ISDD9FNJjdC6xzvt8KHEDirPmBDmZXXWV/3DFJTgLvq6+slfbBB5nfceaDcNi6\nD3futA/CoLU4G3L//fYhGZQlFxs3Ws9BZaVNVCovd9cNuHlztOswXtu27m9S1661SVLxQiGbcJNu\n1yZYQOza1b4Hgd/rzEqBcUSD2QKgN9GUVhVApwTvC2wwmzPH7uYXLrTZSsVk7Fjrjps+HZo39+7n\n7N6d2YeAG3feCW+/DRUVxTW7b+lSW/u0alX9sbtw2GbhZfpvv3u33Ry4zWKS7u978mS48krbVfu+\n+4K3IHzbNvuM+eab7GaEyZVQaM/fSb7NZiy6jPmbN9tdX6T/vLIS7rmn+AIZWIv0rbdsq5j778/e\ndVes2DNryty5lnkktpvH7cC7Gy+8YC2TysriCmRg/44dO9pykt694fPPo//2kyfbEoTy8ui/e1lZ\n4+skN260FG2R39+0aTaRJvb3d8op0YC1Zk20q27SJMuaU1oaPbdnT0tE0NDve+dOmwn66KPw97/b\nBKUgat7c/l0OO8zvkqRn+3ZbdnLJJY2f60fL7FusazGiFhtHixcOh8PccQcsWuRtAb0SmZI9d671\nu0f+Y3bvbgOzxaq21v49ysszn5W3Ywd89hl891399Vrz50c/ZCdOtJ91xhmZtwjDYesqff/9YHSX\npuOee2D0aAtCRx8d/Xfv0cOCe+w6uNmzbZJDhw71g0s4bOmavvzSxpAjwai83P5OYidQrFxp5yxf\nbuNX3btHf27XrvZ/LXL+xIn2QVhWZjNq482bZx/wTzyReBmE5IfRo60na8wYe56P3Yx9iGbM/5AG\nuhmHDRvGokU2IHvCCX3o3LmP96XNsnbt7AO02O7eG7NwoXU1ZqppUwuMxx+fvNUVubH4/HO7K8/U\nySdnPm26kG3caKmzunVrfB1UJEn2ypWJX+/QwbotG5sdWVNjfzPt2llwbGyMctkyO3/79vqvtW5t\nE1gKseutmDzxRAU33VTBkCH2fMSIEZBHwcxNxnwI8JiZiIg0Lhy2lnNlpXWX+pmb8XlgCnA8sBy4\nFhgJ9AO+BPo6z0VERPYQCtmyFzc71+dzI1stMxGRIvfoo5Yq7umn/R8zS5eCmYhIkVu82Cb5rFwJ\nTZpoCxgRESlAHTvaTOQvvkh+noKZiIjkrVAI+va1HdSTUTATEZG85iaYacxMRETy2sqVlqCgtjY/\nx8x+gC2gXoglHBYREanniCMaT9zsVzBrCjyEBbTOwI+BgO61nFxFRYXfRcgZ1TW4iqm+qqs/+vZN\n/rpfwawbsAjbz2wn8AKQx1sSeief/li8proGVzHVV3X1x5VXJn/dr2B2BJYRJGKFc0xERKSeM89M\n/rpfwUwzO0REJGv8ms1YDgzHxswAhgJ1WPLhiJlAl9wWS0RE8tgs4FS/CxGrBFiMZdRvhgWuopwA\nIiIihe184AtsIshQn8siIiIiIiIihe5xYA1QFXOsDfAetm/bu9gmpBFDscXiC4Dvxxw/3bnGQuDB\nmON7Ay86x6cC7bNb/JQcie0OPheYA/zWOR7E+jYHKrFu8HnA3c7xINY1oikwA9tQF4Jb12pgNlbX\nac6xoNYVrC6vAPOxv+Uygl1fSVMv4DT2DGb3Arc4j28luglpZ+zDcS9szHAR0Yk407A1eAD/IjpB\n5r+Bh53HP8LW5fnlUKIDsK2wruITCG59WzrfS7D/pD0Jbl0BbgKeBd50nge1rkuwD/NYQa0rwFjg\np87jEmB/gl1fyUApewazBcAhzuNDnedgdzyx6bvGY7M7D8PumiIuAx6JOafMeVwC1GSr0FnwBnAu\nwa9vS2A6cCLBrWs74H3gbKIts6DWdQlwYNyxoNZ1f+CrBMcDVV9lzffOIVjXI873yB/N4dgi8YjI\ngvH44yuJLiSPXWS+C9hA/btKP5RiLdJKglvfJthd6hqi3atBrev9wO+wZTIRQa1rGAvcnwLXO8eC\nWtcOWHB5AvgceBTYh4DVV8EsN8IEb6F4K+BVYAiwKe61INW3DutWbQechbVaYgWlrv2BtdgYUkPr\nT4NSV4Ae2I3Y+cAN2FBBrCDVtQToinUDdgU2A7fFnVPw9VUw884arOkO1jxf6zxeiU2iiGiH3e2s\ndB7HH4+85yjncaS/uzb7RXZtLyyQPY11M0Kw6wt2p/k2NgAexLp2BwZg3W/PA32x328Q6wqwyvle\nA7yOjQMFta4rnK/pzvNXsKC2mgDVV8HMO28CVzuPryb6of8m1tfcDGv+H4sNqq4GNmL9ziHgSuCf\nCa41GPjA47InEwIew2ZEPRBzPIj1PYjoDK8WQD+s5RLEut6OfYB1wOowAStnEOvaEtjXebwPNluv\nimDWFaycy4HjnOfnYt3l4whmfSUDzwNfAzuwP5prsf7i90k87fV2bIbQAuC8mOORaa+LgNExx/cG\nXiI67bXUgzq41RPrepuJfbDPwGY0BbG+J2NjDDOxady/c44Hsa6xehOdzRjEunbAfqczseUlkaQN\nQaxrRBesZTYLeA1rOQW5viIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiJB8QtsYSrANViWhohH0e7s\nIiJSYD7EFqyKiIjkTCmWPeEZLPXXy1hKrHOwrCKzsbRgzZzzR2IphmZh+0wBDAf+BxiEJXNe4Ly3\nOVBBNLj92LleFdE9qQC+A+7CMl18AhycxfqJiEgRKMXSfp3pPH8M+D2wDDjGOTYW23WgDdG9pAD2\nc74PwzbLBGuZdY05J/L8cGAptjdXUywn3kXOOXXABc7je4D/zaxKIvlLiYZFvLMcaxGBtdD6Ypsk\nLnKOjcW2ldkAbMMC3kBgawPXi9+aJQScgbXSvgF2Y7tEn+W8vgPL9A/wGcqXJwGmYCbindj9oULA\nevYMSJHHu7EtSF7B9hUb7+J6DR0LxRzbGXO8DtuaQySQFMxEvHMUtt08wOXYrsalwNHOsSuxVtU+\nWMbyd7BuxS7O6yGiAW8T0e7HiDC2NUdvot2MlwEfZbUWIgVAd2oi3vkC28X4cWxyx5+x7TFexv7v\nTQMewfZNewOb2BECbnTeH7v775POuVuwjTQjVmO7Bn/ovPctbJ8q2LPVVvA7CYuISO6VYrMLRSQH\n1M0o4h21hERERERERERERERERERERERERERERERERKQQ/D8XOo/7UFLFfAAAAABJRU5ErkJggg==\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": "iVBORw0KGgoAAAANSUhEUgAAAbkAAACgCAYAAACR+hLPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH+dJREFUeJzt3Xl4VOUVx/EvshZBBRfAokJBFKkFXIpYFcRS3OqCtuKG\nQqvUBRSRzQ0EF4oSQeuCIApFUFREEVEECVhkx7CGJbIIWnYVEVlz+8e50xnCJCSz3Tt3fp/nyZOZ\nm8nMe0nImfu+5z0HRERERERERERERERERERERERERERERERE4jIM2AQsLnC8I5ALLAH+GXG8J7AK\nWA78KRUDFBERidWFQGMODnIXA58BZd37x7ufzwBy3OO1gDzgiJSMUkREJEa1ODjIjQFaRHlcT6B7\nxP1PgPOSNywREckEqb5aOhW4CJgFZAPnuMdPBDZEPG4D8OuUjkxERAKnjAevVwW7SjsXu7L7TSGP\ndVI1KBERCaZUB7kNwFj39lwgHzgO+BY4KeJxNd1jB2nYsKGzcOHCZI9RRETSy0KgUbQvpHq6chzh\nNbl6QDlgK/Ah0Ma9Xxub1pxT8JsXLlyI4zgZ8dGrVy/Px6Bz1bnqXHWu6XCuQMPCgk4yr+RGA82A\nY4H1wGPYtoJhWDLKXqCt+9hl2NTlMmA/cDearhQRkTglM8jdWMjxWws5/pT7ISIikhDai+ZTzZs3\n93oIKaNzDSadazCl27mW8noAJeS4868iIiIAlCpVCgqJZ7qSExGRwFKQExGRwFKQExGRwEpmkCus\nCwFAF2wjeNWIY+pCICIiCZXMIPc6cGmU4ycBLYF1EcfOAG5wP18KvJTksYmISAZIZiD5Avg+yvEs\noFuBY1djm8f3AWuxVju/T+LYREQkA6T6aulqrH7logLH1YVARCSFMmU3ViqDXEXgIaBXxLGi9ull\nyI9ARCS1Ro2CP/whMwJdKrsQ1MGaqIbaCNQE5gNNKGYXAoDevXv//3bz5s3Tbve9iIiXfv4ZunWD\nUqVg/Hi46iqvR1Ry2dnZZGdnF+uxya54UgsYD5wZ5WtrgLOB7VjCyShsHe7XwGSgLodezaniiYhI\nHHr3hpUr4frr4amnYO5cC3jpzKuKJ6OBL7GWOuuBdgW+HhmtIrsQTERdCEREEm7DBvjXv6BfP7jm\nGtizByZO9HpUyZVu8VtXciIiMWrbFk4+GZ54wu6PGQNZWTBzZnpfzal2pYhIhps3DyZPhu7dw8eu\nuw527IDPPvNuXMmmICciEnCOA507Q9++ULly+Hjp0vDII/D448HNtFSQExEJuLFj4aef4PbbD/3a\nDTfA1q0wdWrKh5US6TYLqzU5EZES2LMHzjgDXn0VLrkk+mNGjIBhw6CYWfm+ozU5EZEM9fzz0KBB\n4QEO4KabLPNy2rTUjStVUt2F4BkgF9sQPhY4OuJr6kIg4jP/+Q9MmeL1KCRWW7ZA//7w7LNFP65M\nGXjoIVuzC5pkTldeCOwERhDeDN4SmIK12ennHutBeDP4uYQ3g9dzHxdJ05UiKbJokb37L1sWli+H\no47yekRSUnffDeXKwcCBh3/svn1Qrx6MHGklv9KJV9OV0boQfEY4cM3GyneBuhCI+Mp338GVV9rG\n4ZYt4emnvR6RlNTSpfDuu/DYY8V7fNmydjXXp09yx5VqXq7JtQc+dm+rC4GIT+zcaQHurrss8+7p\npy1pYfVqr0cmJfHgg/Dww1C16uEfG3LbbXbVPmtW8saVal4FuYeBvdgUZWE0LymSYvv3Q5s2cPbZ\n0KOHHTvxRNtj1a1gF0jxrU8+ga+/tjcqJVGunP3cg7Q2l8ouBCG3A5cDkbk+6kIg4jHHgfvus7WZ\nl146uMxTly5Qv75l3zVr5t0Y5fD277ef17PPWtAqqfbt4cknrULKOeckfnyJ4OcuBJcCA4BmwNaI\nx6kLgYjHsrLg9dcto/Loow/9+ttvW2HfefOsUob408svwzvvWFZsrPUoX3jBSoB98EFix5YsRSWe\nJDPIjcaC2XHYVoJe2DaBclh7HYCZWMcBsIaq7YH9wH3Ap1GeU0FOJAnGjoVOneDLL62AbzSOAxdd\nZOs2f/97ascnxfPjj3DaaTZd2ahR7M/zyy9Qpw5MmACNGydufMniVZBLBgU5kQSbPdsSTT79FM46\nq+jHzp9vj12xQlsK/KhbN9i+HYYOjf+5nnvOrurfey/+50o2BTkRiWrNGtsT9eqrFryKo107OOEE\n+Oc/kzs2KZnVq+H3v4fFi6FGjfifb9cu+M1vrEPBmdHaXvuIynpJ0jiOzf1fc41NeUn6+P57uPxy\nSzMvboAD6yY9dKhl74l/dO9uWbCJCHAAFStaAkuo91y60pWcxGTvXnjrLUtW2LfPat8NGgQLFkDN\nmof/fvHW3r3QqpWtt2Rllfz7n3rKElD0xsYfvvgCbrnF9rj96leJe96dO21tbupUK/LsV5qulIT5\n/nsYPNiyr+rXt3d6rVrBEUfYO77sbJg0ye6LPzmOJY/s3GlZeLFkSu7ebT//YcPg4osTP0Ypvvx8\naNLEruJuuinxz9+vn5V4G1XUrmaPabpS4vb119Cxo72ry821rKvJk+Gyy8IBrUcP+PlnKwUl/tWn\nj73jHzky9q0AFSpY4d/774cDBxI7PimZN9+0n+ONNybn+e+5x9blli9PzvMnW6q7EFTF6leuBCYB\nx0R8TV0IfOjLL+G66+ydYqVKtqg9fHj09OQyZeDf/7ZqCbm5qR+rHN6IEfDGGzB+vK25xOP6620/\n3WuvJWRoEoOff7Z6k1lZse+JO5zKla1IwFNPJef5ky3VXQj6Y5vA+wPdgSqoC4Hv7N8P48bBgAGw\nebO9W2/XzoJccQwebNl6M2fGVnFBkiM722pRJnJ9ZcECS15ZsSL6BnJJrj59YNkyWx9Pph9/tFmc\nWbOgbt3kvlYsvFyTq8XBFU+WYxvENwHVgWzgdOwqLh8IJSV/AvQGCpYJVZBLop9+sjWWQYMsQ6tL\nF7j66pJPaTmOZeuddVawauCls9xcaN4cRo+GFi0S+9x/+5sVAX7mmcQ+rxTt22/hd7+zvYu1aiX/\n9Xr3hm++sb8RfuOnNblqWIDD/VzNva0uBB7asMHSj2vXts2fo0bBjBnQunVsazalStkU1pAhdjUn\n3tq0Ca64wtbQEh3gwOocvv46rFqV+OeWwj3yCNx5Z2oCHNiU5Qcf2N7KdOJl4olD0Z0GdMmWZDk5\ncOut9m5w926YO9ey7c47L/7nrl7divy2bWtZfOKNXbvgqqvs53zbbcl5jerVra1L167JeX451IIF\nVrqrZ8/UvWaVKtbVIN16C6a6C0FomnIjUAPY7B5XF4IUyc+HiRNtvW3lSqtX+MILcMwxh//ekmrd\nGj780P4AvvJK4p9finbggO2dqlfPppqS6f77bZ1vyhTrJi7J4zi2XeDxx1NfWq1zZzj1VCsgcMop\nqX3tSH7uQtAf2IatvfXAsisjE0/UhSBJdu+2zMfnnoPy5W297a9/TX5iyI8/QsOG8OKLNmUmqdOl\ni63XfPqp/cyT7b33LJh+9ZVl2kpyvP8+9Opl/85edIPo0QN27LCZGr8oak2uOKYU81hBo4HvsOao\n64F22BaCyUTfQvAQkIclp7Qq5DkdKZlNmxynd2/HqVbNcS6/3HGmTHGc/PzUjiE723Fq1HCczZtT\n+7qZ7MUXHee00xxn+/bUvWZ+vuM0a+Y4L7+cutfMNLt3O06dOo4zaZJ3Y9i82XGqVHGc9eu9G0NB\nFLG8VVTk+xVQEZgKNI84fhSW/Xh6Ed+bLO75yOEsX257Z955B/7yF5tmqF/fu/F06wZ5efZuP1n7\necRMmGCtcGbMsAK7qZSTA5dear9/yZgCz3RZWfD55/DRR96Oo2tX2LMHnn/e23GExLqF4H6sr9uJ\n2BVZyE/Aq4AXdS0U5A7j22/hH/+AOXNskfjuu61ivNf27IFzz7UptGQlQIglJLRqZZu9E5FAFIs7\n7rC1ogEDvHn9oNq61d6oTp/u7RtWsIzd+vVhyRI48URvxwLx75PrBPgkXivIFWX9eksRv+UWu3JK\nZKHWRFi0yJIS5s5NXdpzJlm/Hpo2tX2O113n3Tg2bYIGDaxaTr163o0jaO6912ZBXnjB65GYzp3t\n83PPeTsOSMxm8POxJJLI5eQRcY0qNgpyhfjmGyuUe8898MADXo+mcM88Y1cZU6d6s2geVDt2wAUX\n2JaNBx/0ejS2J+8//7HsWolfbq51Zc/NheOO83o05rvv4Le/tYor1at7O5Z4g9xI4DdADhBZirVj\n3CMrOQW5KNatswDXqZOlcvvZgQN2tXnlldpXlSj79tm/Z506lsXqhzXPPXtsS8Err0DLll6PJv1d\ncYXNgvjtDWzHjlas2+tqN/EGuVwsxd8P0UVBroC1ay3Ade5sQS4drF1r63NTpthGdImd40CHDrYW\n+8EH/krdf/99ePRRS0bx07jSzaRJNkOzdKn/asFu2GD/h1esgOOP924c8Zb1WoJt3E6knsBSrEPB\nKKA8RXcokCjWrLEA16VL+gQ4sPW4Z56xtcM9e7weTXrr39/WON96y3+B5JprLOnp1Ve9Hkn62r/f\n/n/37++/AAfWIPmGG2JrvJsqxbmSywYaAXOA0J8kB7gqxtesBXwO1Hef723gY6AB0TsURNKVnGv1\nagtw3btbBmW6cRxLjqhb1/4DS8m9/batv82aBb/2aaXXhQvhT3+yLQVVqng9mvQzeLAV1Z461R/T\n0NGsW2fF2FeuhGOP9WYM8U5XNi/keHZsw6EqMBM4D9uO8D6WvfkC0TsURFKQwxqYtmhhdev+8Q+v\nRxO7LVusGsro0dCsmdejSS8zZsC111ozy4YNvR5N0Tp0sN51fsjCSyc7dlh26scfWxDxszvusOQT\nr7qOeNlqpzB3AgOAX4BPgVuB77Grt9C4tkfcD8n4IJeXZwEuVIE83U2YYKnRCxemvg5fusrLs0zK\nN96wjdd+t3mzJaHMmAGnneb1aNJHjx62HeP1170eyeGtXm3r7Hl53lyxxxvkdhJOOikHlHWPxfon\nqQ5Wz/JC4EfgHeA97Eou8p9nO3bVFymjg9zKlZZh1auXVbQIig4dYO/e9PjP7LVt22wvXJcu9u+W\nLp591pq2el2pI12sWQPnnAOLF/tjs3VxtGtn6+29eqX+tYsKcsVZqo7sB30EthYXTy2Fc4AvsULN\nAGOBplhngmgdCg6SqV0IVqywANenD7Rv7/VoEmvAAGjUCMaOtc4FEt3u3ZbMce216RXgwBKjBg+2\nYtGtCqtMK//Xo4f1b0uXAAfWmaBpU9vGlOwu8anoQpCDJaPEoiHwJnAusBt4A0tqOYXoHQoiZeSV\n3PLl8Mc/whNPwO23ez2a5Jg50/545+R4v7HUj3butA7cjmOZlEd42QkyRh98AA89ZFPTfssE9ZMZ\nM6BNG3tjW7Gi16MpmbZtbR3xkUdS+7rxTldGFgg6AjgbSxBpGseYugG3AfnAAuDvQGVgDHAysBb4\nK/BDge/LuCCXm2sB7umn7RcoyB591NqHjB/v30yyVNu92zZU9+tnvwdDhvivXFtxOY5tDL/2Wtv3\nJYfKz7eroXvvtUa36WbFClsvXr0aKldO3evGG+TeILwmtx8LQEMoZDoxyTIqyC1bZn/Y+ve3PWVB\nt2+fFRXu0CEYSTXx2LfP1ij79rXMur59g7FxfvFi+53OzYWqBVfchTffhIEDYfbs9LxaB7jpJvtd\n7VFwHi6J/JhdGauMCXJLltj+omeegZtv9no0qbNsmdXomzXL9tBlmgMHbEtF795Qu7ZNUTdp4vWo\nEuuuu2xj86BBXo/EX3btgtNPt0B34YVejyZ2S5daBvjXX0OlSod/fCLEW/HkJGwv2xb34z2gZqIG\nJ4davNimdQYMyKwAB5Zq/thjNlWzf7/Xo0kdx7HEm9/9Dl5+GYYOtT1wQQtwYMlTo0bZ1ZyEZWXZ\nzzudAxxYB4pmzWya3Q+KcyU3GUsUGenev9n98KLsauCv5BYutL1PAwdauZxMlJ9vGXjNmqV+ATvV\nHAc++cTO03Hsyu2yy4K/JpmVBZMn20ZnsYr+Z55pJdpS3eg2GRYtspmo1atTkzwT73TlQiwj8nDH\nUiHQQS7UVfmFF6ybdybbsMHWoiZOhLPP9no0yTFtmgW3bdvs6qZ16/RdhympvXutTcugQRbUM93f\n/mYlsYJU4q51a1t6SEVnlKKCXHF8jlUkKY3tq7sFmJKQkZWcE1Tz5ztOtWqO8+67Xo/EP0aNcpzT\nT3ecXbu8HklizZnjOC1bOk7t2o4zfLjj7N/v9Yi88eGH9vPdu9frkXhrwQL7v//DD16PJLEWLHCc\nGjVS8/+XIrrkFOd9YzssnX8j8F/gL+6xeBwDvIu18VkGNCGDuxDMn2/vZl96yduOzn5z4422STyV\nWVrJtHhxeDN369a2/7Ft28xtHnvllXDSSf5Zu0mV/HxYtQrGjLF9g23aWJWQZG+gTrXGja1qy2uv\neTuO4lzeDQfux2pLggWjZ4F46m4MB6YBw7CrwyOBh8nALgTz5llDxMGD7Q+gHGz7ditAPGxY+jbf\nXLXK/ohNmWJdI+66K333uiXakiWWiZeb610F+2TavdvOMSfHPr76ytarjj3W3sA1bmzT8ZdfHsyp\n6nnz7E1dXh6UL5+814l3TS5adZN4Kp4cDXyFdRuPtJwM60Iwd669mx0yBK6KtXFRBvjsMytltmhR\nerVr+eYbW2sbN87WJe67L7UbZNPFPffY1ezzz3s9kvhs3x4OZqGAlpcHp55qwaxRo/BHOv0ex+uK\nK+DPf05ux5REJJ5cjBVMBruSmwacGeN4GgGDsWnKhsB87EpxAxnUhWD2bAtsr71mgU6K1qkTbN1q\nqed+t3EjPPWU7Xfq0MF6vmnjc+G2boX69S0R54wzvB7N4TmO9VCLDGY5OfD99zbrEApkjRvb+VSo\n4PWIvTVrlmWKr1qVvMav8Qa5tthU4hj38X8BngRGxDiec7B+cucDc4GBWF+5e8mQLgQzZ8LVV1tF\niyuu8Ho06WHXLpvWeewxW6vzo+3bLTvu1Vdtra1nT6hWzetRpYeBA61488SJXo/kYPv22VRqZDDL\nybHp5sirs8aNLfU/iFOOidCqFVx/vfWdS4ZEVDxpALTAMlg+x67CYlUdC3K13fsXAD2x6cuLCXch\nmEqU6cpeEX0c0rELwZdf2trb8OFKnS6pefNs7WLBAqjpo3IEO3bYH+nnn7eEkkcftYQKKb59+2yf\nWFaW/Yy9sH27VdyJDGi5uXDKKQcHs4YN9ealpGbMsNKEK1dC2bLxP1/BLgSPP/44+Kys13SsKPNK\noDcQ2i4Y6C4EoW7OI0akR7NLP+rbF774wjZQe/2u+Zdf4MUXrfRay5ZWiisTS5ElyoQJ1idv8eLE\n/CEsyHFsajQvL/yxalX49oEDVlYrFMwaNbLAe+SRiR9LJrrkEgt07eLNzY/Cj7UrGwJDsSasX2Nb\nEkoT4C4EX3xh2wNGjrRKABKb/futyvnNN0PHjt6M4ZdfbKr5ySetDFOfPraxWeLjODa7cdlllqQT\n63Ns3hw9iOXlWSWZU0+1NyN16x58+7jjgl9pxkvTptmm9+XLE99qyY9BLlZpGeSmT7f56FGjrAK7\nxGfVKjj/fPt3rV8/Mc+Zn2/TVf/9ryWOFPyIPL5zp70r7dvX9gFJ4ixbBs2b2+fjjov+GMexn0Nh\ngaxcuUMDWOi+EoC8dfHFloSS6ExLBTkPZWdbia633rI/jJIYgwdbgsfMmUVnbO3aVbzAtXmzpffX\nqGFNWyM/Ch6rUsX7qdIg69jRpg4ffrjwQFaxYvQgVqdOZqXnp5tFi+yN/pIlcMIJiXteBbkEOnDA\n6u4V52PdOltjGDPG3sFI4jhOuGJGkybRA9fGjZbQEC1QFTx2wgnJ3awqxbdtmyV3HDhwaBCrW9cC\nWdCqg2SSrl1h0ybLTUiUjA9yy5ZZgkAo+OzZU/xAVfAjP9+uHKJ9lC9/8P0KFaxsT7NmSfiXEDZu\ntP1zFSsWHryOOkrrLOnIcfRzC6qdO23/4PDhiXvzn/FBbt06+Oij4genoj5Kl9Z/PhGReIwbZ/tI\nFy5MzAbxjA9yIiLiH45jFZ+aNrXZrnj5NciVBuZh5bz+jFU3eRs4hQBuIRARkbC1ay07ec6c+BvF\nFhXkvMwRuw+rnBKKWj2wVjv1sH51AWmwIiIiBdWqZXVd773XruySxasgVxO4HNsQHoq+V2EteHA/\nq/GMiEiAPfCA5UyMHZu81/AqyD0HdAXyI45Vw9rs4H5WdTgRkQArVw5eftlaUf30U3Jew4sgdyWw\nGespV9iaYJHtzEVEJBguusg2iEfU3k+oBFcQK5bzsanJy4EKwFHAvwk3Sw11Idgc7Zt79+79/9vp\n2IVAREQO1r+/1X9t29YKYx9OwS4ERfF6C0Ez4EEsu7I/Ae9CICIi0Q0dak2kZ8woedk8v2ZXhoSi\nVj+gJdZ+p4V7X0REMkD79lZoY+jQxD6v11dyJaUrORGRgIq1gLNfN4PHQkFORCTAHnwQtmyx2pbF\npSAnIiJpIVTAecQI6y1YHH5fkxMREQGgUiUYNAjuuss6v8RLQU5ERHzlmmusd+Czz8b/XJquFBER\n3ylJAWe/TVeeBEwFlgJLgE7u8apYgeaVwCRsn5yIiGSgUAHnjh3jK+DsRZDbB3QGGgDnAfcA9VEX\nAhERifDAA3ZF9/77sT+HH6YrxwH/cj+aES7vlQ2cXuCxmq4UEckg06fDzTfDsmVQuXL0x/h5C0Et\nYBrwW+AboIp7vBSwPeJ+iIKciEiGuf12qFoVsrKif92vQa4SFuD6Yldz33NwUNuOrdNFUpATEckw\nW7ZAgwYwaVL0As5FBTkvuhAAlAXew7oPjHOPqQuBiIgc4vjj4cknbe/cjBkwfbq/uxCUwjp/b8MS\nUELUhUBERKLKz4cLLrCpyzvvPPhrfpuuvACYDiwi3IGgJzAHGAOcDKwF/gr8UOB7FeRERDJUYQWc\n/Rbk4qEgJyKSwaIVcFaQExGRQAgVcH7zTbjwQjumICciIoGxdCnUqQMVKth9BTkREQksv9WuFBER\nSQkFORERCSy/BblLgeXAKqC7x2MREZE056cgVxor0nwpcAZwI9adICMVdzd/EOhcg0nnGkzpdq5+\nCnK/B/KwjeD7gLeAq70ckJfS7RcpHjrXYNK5BlO6naufgtyvgfUR9ze4x0RERGLipyCnvQEiIpJQ\nftondx7QG1uTA6tnmY8VbA7JARqmdlgiIuJzC4EoTXj8pQzwNdZItRwW0DI28URERILnMmAFloDS\n0+OxiIiIiIiISCYahnU8XxxxrCrwGbASmIQ1hw3piW2EXw78KeL42e5zrAIGRRwvD7ztHp8FnJLY\n4ZfIScBUYCmwBOjkHg/i+VYAZmNT6suAp93jQTzXkNLAV8B4935Qz3Ut1uvyK6zHJQT3XI8B3gVy\nsd/jJgT3XCVJLgQac3CQ6w90c293B/q5t8/A/miWxdYl8wgnBs3B9hECfEw4Oedu4CX39g3Y3kKv\nVCe88FsJm3auT3DPt6L7uQz2H/gCgnuuAA8AbwIfuveDeq5rsD/0kYJ6rsOB9u7tMsDRBPdcJYlq\ncXCQWw5Uc29Xd++DvUuKLGX2CZZxWgN7pxXSBngl4jFN3NtlgC2JGnQCjAP+SPDPtyIwF2hAcM+1\nJjAZuJjwlVxQz3UNcGyBY0E816OB1VGOB+5c/bRPLlNUw6YwcT+HfqFOxDbAh4Q2wxc8/i3hTfKR\nG+j3Az9y6LtQL9TCrmBnE9zzPQJ7Z7uJ8DRtUM/1OaArtqUnJKjn6mABfR5wh3ssiOdaGws6rwML\ngCHAkQTwXBXkvOUQvE3wlYD3gPuAnwp8LUjnm49Nz9YELsKuciIF5VyvBDZja1SF7asNyrkC/AF7\ng3YZcA+25BApKOdaBjgLm048C/gZ6FHgMYE4VwW51NuETQOAXepvdm9/iyVvhNTE3iF9694ueDz0\nPSe7t0Nz6tsTP+RiK4sFuH9j05UQ7PMFe3c6AVt8D+K5ng9chU3jjQZaYD/fIJ4rwH/dz1uA97G1\npiCe6wb3Y657/10s2G0kYOeqIJd6HwK3ubdvIxwMPsTms8thUwmnYgu6G4Ed2Nx2KeBW4IMoz3U9\nMCXJYy9KKeA1LEtrYMTxIJ7vcYSzzn4FtMSudIJ4rg9hf9xqY+fwOTbOIJ5rRaCye/tILINwMcE8\n143YVGI99/4fsSn38QTvXCWJRgPfAXuxX6h22Jz0ZKKn6D6EZS0tB1pFHA+l6OYBz0ccLw+MIZyi\nWysJ51BcF2BTeDnYH/yvsCyrIJ7vmdg6Rg6Wbt7VPR7Ec43UjHB2ZRDPtTb2M83BtsGEClIE8VzB\nSiTOxUpijcWutIJ6riIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiLiLx2wTbcAt2OVKUKGYN0fRERE\n0t5UbDOuiIiI52phFSNGYiXQ3sFKg12CVVFZhJVHK+c+vh9Wbmkh1usLoDfQBbgOK4K93P3eCkA2\n4aB3o/t8iwn3BQPYCTyBVfaYCZyQwPMTEZEMVgsrf9bUvf8a8AjwDVDXPTYc6+JQlXA/L4Cj3M+9\nsCamYFdyZ0U8JnT/RGAd1hutNFY38Gr3MfnAFe7tfwIPx3dKIv6lAs0iqbceu4ICu6JrgTWwzHOP\nDcfa9/wI7MYC4bXAL4U8X8EWOKWAc7Grum3AAayr90Xu1/dinRMA5qOaghJgCnIiqRfZo6sU8AMH\nB6rQ7QNYq5d3sb5unxTj+Qo7Viri2L6I4/lYGxSRQFKQE0m9k4Hz3Ns3YV2oawF13GO3YldhR2JV\n4Cdi05MN3a+XIhwIfyI8jRniYG1QmhGermwDTEvoWYikAb2DE0m9FVjX6WFYUkkW1orkHez/5Bzg\nFaxv3TgsoaQU0Nn9/siOzW+4j92FNTgN2Yh1ep7qfu9HWK8wOPgqLxDdn0VExB9qYdmOIpICmq4U\nST1dOYmIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIgkzv8AyFJ2NH9haz0AAAAASUVORK5CYII=\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+/AAAHO5JREFUeJzt3XucVHX9x/HXwoJyEZTECwIuihdIxTSFn4KgogJqJVCm\nppl5SbO0fmmmFV5+5uXx+2mapaCWWoakkHiDigyDQgXlsorcBRTiEgquXHbZ3fP743NOc3Z2Lmd2\n5uzMnPN+Ph7z2Jkz58x8v3v7zPf6ARERERERERERERERERERERERERERERERkbIzAlgCLAd+mOac\nB93nFwKf8x3/EfAuUA38HtgjvGKKiIjkpy2wAqgC2gELgH5J54wCXnHvDwRed+9XAatIBLpJwNfD\nK6qIiMRBmxBf+0Qs6K0GdgPPAF9MOucLwJPu/TeAvYH9gU/cazoCle7XdSGWVUREYiDMoHcQ8IHv\n8YfusSDnfAT8H7AWWA9sBWaEVlIREYmFMIOeE/C8ihTHDgWux7o5ewCdgYsKUywREYmryhBfex3Q\ny/e4F9aSy3ROT/fYMOCfwBb3+BTgJOBp/8UDBgxwFi5cWLgSi4hIFCwEjk31RJgtvXnAYVhrrT1w\nPvBC0jkvAJe49wdh3ZgbgaXu4w5YS3A4sDj5DRYuXIjjOLG4jRs3ruhlUF1VV9VVdS2HugID0gWm\nMFt69cC1wJ+wmZyPA+8BV7nPj8dmbo7CJrxsB77hPrcAeAoLnI3A28CEEMsqIiIxEGbQA5jm3vzG\nJz2+Ns2197o3ERGRggize1MKaNiwYcUuQqtRXaNJdY2mcqtrqpmT5cRx+29FREQAqKiogDTxTS09\nERGJDQU9ERGJDQU9ERGJjbCDXkuzLBwBzPfdtgHfDbWkIiISeWFOZGmLLTIfju2yMhe4AFur5xmF\nLVkYhWVZeABblO7Xxr3+RJru0wmayCIiIkmKNZElnywLfsOBlTQPeCIiIjkpxSwLPZPO+SqWRFZE\nRPKweDE8+GCxS1FcpZhlwX9de+Bc4NmClEhEJMZmzYIpU4pdiuIq1SwLnpHAW8DmdG9y6623/uf+\nsGHDym53ABGR1rJ6NXz0UbFLUXgzZ85k5syZgc4NcyJLJTaR5XQsEeybZJ7IMgj4OU0nsjyD7d35\nJKlpIouISEAXXGCtvQ+Tmx8Rk2kiS6lmWQDohE1iuSLEMoqIxEZUW3q50N6bIiIxceCBsGED7NgB\nHToUuzTh0d6bIiIxt3OntfK6d4ePPy52aYpHQU9EJAbWroVevSzoxbmLM+wksiIiUgJWr4aqKqit\njXfQU0tPRCQG1qyxoNetW7yDnlp6IiIx4LX0GhriHfRKNcsC2D6cz2HLHBbTfCNqEREJyAt6n/mM\ngl5Y2gIPYYGvP7YwvV/SOaOAvsBhwJXAw77nHsDW8fUDjqHponYREcmBF/Ti3r1ZqlkWugJDgF+7\nz9VjOfVERKQFFPRMqWZZ6IPtt/kb4G3gUaBjaCUVEYmwXbtgyxZbnN6tm92Pq1LNslAJHAf8yv26\nHbipcEUTEYmPtWuhZ09o21YtvVLNslDhnjvXPf4caYKesiyIiGTmdW1CNINeVLIs/B24HFgG3Ap0\noPkMUO29KSKSxYQJ8MYb8PjjFgCHDrV1e1FVrlkWvgM8jSWSXZn0nIiIBOQtTIdotvRyoSwLIiIR\nd9FFMGIEXHwxOA60bw/bt9vXKFKWBRGRGPOP6VVUwD77xDfTgoKeiEjE+YMexLuLU0FPRCTCamvh\n3/+GHj0SxxT0REQkkvxr9DwKeiIiEknJXZugoBemfLIsrAYWAfOxNX4iIpKj1avh4IObHotzpoUw\n1+l5WRaGY7uszAVeoPnidC/LwkAsy4K3ON0BhgEx/dGIiORPLb2mgrT0vhDwvGT5ZFnwlPs6QhGR\novIvTPco6GV2Pha87gWOzOG1W5plwTvHAWYA84ArcnhfERFxpWvpxTXTQpDuzYuw/HYXAE9gweg3\nwESgJsN1Lc2y4BmM7dnZHfgLNjY4K+BriogI6t5MFnRMbxuW6aADcD1wHnAjNgnlwTTX5JNlASzg\ngeXV+yPWXdos6CnLgohIarW1sHlz0zV6EL2gV+gsC18ELsUmmzyFtfY2YUldFwNVaa7LJ8tCR2wi\nTA3QCfgzcJv71U97b4qIpLFiBZx5Jqxa1fT4ypVwxhnNj0dFvlkWRgP3Y6l+/HZgqX/SySfLwgHA\nFF8Zn6Z5wBMRkQxSdW1C9Fp6uQgS9DbSPODdg627m5Hl2mnuzW980uNrU1y3Cjg2QNlERCSNdEGv\na1f49FOor4fKMBeulaAgszfPSHFsVKELIiIihZUu6LVpY4Fv69bWLlHxZQp6VwPVwBHuV++2Gtsp\nRUREStiaNc13Y/HEtYszU8P291jX5N1YV6Y3KFgDxHSFh4hI+UjX0gMFvVQcrFX3bZqvueuGtgcT\nESlpCnrNZQp6E4GzgbdIvdC8TyglEhGRvNXVwaZNcFDyPliuuAa9TGN6Z7tfq7AAl3wLIp8sC2BL\nHeYDLwZ8PxERAT74wBalp5udGddMC0Fmb54MdHbvXwzcB6QZGm3Cy7IwAuiPLUzvl3SOP8vClViW\nBb/rsAXwWoEuIpKDTF2boJZeJo9gC9EHAN/H1tA9FeC6fLMs9MSC4mMo24KISE4U9FILEvTqgUbg\nS8AvsdbbXgGuyzfLwv3ADe57i4hIDhT0UgsS9GqAm4GvAS9h3ZbtAlzX0iwLFcA52P6e81M8LyIi\nWQQJenFMLxRkA5rzgQuBy4ANQG/gfwNcl0+WhTFY1+coYE+gC9aleknymyjLgohIc6mSx/pFqaVX\n6CwLLZVPlgW/ocAPgHNTvIeyLIiIpNC7N/z97+kD37JlcPbZsHx5qxarVeSbZWEMtivL/r4XcbDW\nVyb5ZFlIpsgmIhJQXR1s3Ag9e6Y/J0otvVwEaemtxMbY3st2YhGopScikmTVKjjtNBvXS6e+Hvbc\n0wJkmyCzO8pIppZekKpuoDQDnoiIpLBsGRx6aOZzKiuhc2fYtq11ylQqggS9ecAkbDxujHsbHWah\nRKS4HnoI9t039W3SpGKXTrKproajj85+Xhy7OIOM6XUFdgJnJh2fkuJcEYmAiRNhwgQ45ZSmx//8\nZxg/Hs4/vzjlkmCqq2Ho0OzneUEvW6swSoIEvUvDLoSIlI6PP7Z/mqNG2ZiP33nnwTXXwObN0L17\nccon2VVXw3e+k/28OLb0gnRvHgH8FXjXfXwM8OPQSiQiRTVjBgwZ0jzgAXToACNGwPPPt365JJj6\neli6FD772eznKuil9ii2I0ud+7gaG98LoqVZFvbE9uJcgG04fVfA9xORPE2bBiNHpn9+7Fh47rnW\nK4/kZtkySyfUsWP2cxX0UuuIBSCPg20gnU0+WRZ2AacCx2Ity1OBwQHeU0Ty4Dgwfbq15tIZORJe\nfz1+/yzLRdBJLBDP9EJBgt5mLDB5xgL/CnBdvlkWdrhf22MBNGY/GpHWV11tLYS+fdOf06kTnH46\nvPBC65VLgssl6Kmll9q1WHqhI7DtxL4HXB3gupZmWfD2EGiLdW9uBP6GdXOKSIiydW16xo6FyZPD\nL4/kTkEvs0yzN//bd38aFnjaYC2w0Vgy2UxammXBu64B697sim1lNgyYGfA1RaQFpk+HH/wg+3nn\nnANXXw2ffAJdsm1IKK2quhqOOSbYuXHMtJAp6O2FBaAjgBMArzPjYmzz6GzyybLgtw14Gfg8KYKe\nsiyIFEZNDcybB0H+hLp0sTV8L70EF14YetEkoJoa23Mz6Lq7qLT0Cp1lYRY24aTGfbwXtlH0kCzX\n5ZNlYV9sw+qtQAespXcbtnTCT3tvihTI1Knwy1/aAvQgnngCXnxR3ZylZM4cW583b16w8997z9Ze\nLlkSbrlaW75ZFvaj6WzN3e6xbPLJsnAgNsGljXv7Lc0DnogU0LRpmWdtJvvCF+C662D7dpvcIsWX\ny3geRKell4sgQe8prJU2BYucXyIx4zKbae7Nb3zS42tTXFcNHBfwPUQkT95She9+N/g13brBoEEW\nLMeODa9sElyuQW+ffWwHHseBijCzq5aQILM378RaYFuxZQOXAj8LsUwi0sqWLrV/fP2SV9JmMWaM\nFqqXklyDXvv2tvNOTU32c6MiSEsP4C33JiIR5HVt5vpp/0tfghtvhJ07bYsyKR7HgUWLgs/c9Hhd\nnHGZhRux1IFSqjZuhClT4Cc/sWnuUlqy7cKSzn77wec+F3zyi4Rn/XrLkbf//tnP9YvbuF7Qlp5I\nM44DH35oYwLJ6uvtU+esWXbbvBlOOsn+uDp1gptuav3ySmo7dsA//wnPPtuy68eMsRmcX0zeb0la\nVa5dmx4FPZE0Ghvh3XctiM2ebV9ra+GAA5qfW1EB/fvbbv3XXQdHHQVt2sA778AZZ9gxdYeVhtde\ng+OPb3n31ujR8NOfQl2djRFJcSjoBdMaQW8Etv6uLfAYcE+Kcx4ERmK7vVwKzMcWrT+FLY9wgAnu\nedLKGhutZfboo5Y5e8gQGD4cbr0VDjsst3Ggo46CE06wNV5XB9nMTkKX61KFZD162ASYv/412BZm\nEo6giWOTKegVlpdpYTi208pcbGeX5AXqXqaFgVimhUHYesDvYftvdsYm0vwl6VoJWX09XHYZvP++\ntdIOSt49tQVuugm+9jW44gobg5Dimj695V2bHm8vTgW94lm0KFji2GRxy7QQ9kSWfDItbMACHsCn\nWLDrEW5xxa+2Fr7yFdi0Cf70p8IEPLCxvV694A9/KMzrScutXGnT1XOd8Zds9Gjb0aW+vjDlktzs\n3h08cWyyuLX0wg56+WZa8FRhCWbfQFrF9u1w7rk2Djd1arCElLm46Sa4+26bDCPF88AD9nPOd2Hy\nwQdDnz42Piitb/ly6NmzZX+nCnqFlW+mBbCuzeeA67AWn4Rs61Y46ywbq3nmGdhjj8K/x4gRFlBf\neaXwry3BTJhgLfh7Uo2yt4AWqhdPSyexQPQyLdxyS+bnwx5RyTfTQjtgMvA74PlUbxC3LAuOY11S\nffpA27bZz//4Y/sUd8ghwT7Nr1ljG9CefLK1AtqE9LGooiLR2jv77HDeQ9KbMcNmXM6ebVtRFcKY\nMTB4MDz0ULDfTSmcfINeubf0vCwLu3bBfdmS3oWsEliJdU+2x8bokjc6GoVtPA02geV1934FNnvz\n/gyv78TN1KmOs+eejtOli+OcdZbj3HGH48yc6Tg7dtjza9c6ztNPO87VVzvOUUc5TufOjnPQQY6z\n336OM3q049x3n+PMnes4u3c7TmOj4yxe7DgTJjjOJZc4Tp8+jtOtm+Pcdps9F7bdux3nkEMcZ9as\n8N9LEhYvdpzu3e33ptAGDHCc114r/OtKZuee6zjPPtuya6urHad//8KWp1imTnWc4cMdhwy9jK2x\nxehIEksWHgfuommmBbAZniNIZFp4GxgM/B1YRKICPwKm+17bcWI0KNTQAAMGWOto4EBbUOytmauu\nhq5dbSLBkCH2iXvwYDj2WGjXDtauTaytmz3bWnTt28Nee9l53jVHHhle6y6V8eMtPc1LL7Xee8bZ\nv/9tvzs/+QlcemnhX/+OO+w9Hnig8K8t6fXpY13Vhx+e+7Xr19s6zX/9q/Dlam3XX2870tx8c/rU\nQuW+r3asgt5vfwuPPGJBK7mrcscO2+qrqipYN+ZHH9l+iYWakdlSu3ZZ1+v06fnPIJTMamttfeXg\nwXDXXeG8x3vvwZln2oeq1vzwFGc1NbZBxCeftKxbeedO2Htv+1ss90wLRx8Njz8OAwemD3r6tSwT\ndXU2BnP33al/MTt2tE97QX9pu3UrfsAD2+H9+usLN5kiigrxuc5x4PLL7VPwnXfm/3rp9OtnvQdv\nvhnee0hT77xj3/eWjqN26GDX7thR2HK1to0b4YMP4LgsSem0NLhMTJhgv9hDsuWrL0Pf+hb07WvZ\nnj//+WKXpvjWr2/aFb14sbWC/d3Q+2VJ49zYaNmw/VvGHXig7ZoSdgvMW6g+aFC47yMmn0ksHm8y\nSzknA545E045JfuGFwp6ZeDTT+3T+bTkdLwR0aUL/OIXcNFF8Pbb5f2Hl01DA9x+u3X/Jdu1C+bO\ntSUj3pjsww/bHqYLF1rweuwx2yFn//1tHCbVXpdbtsCcOTbGO3iwbU11yy1wxBGt0301ZoylHLr3\n3vLvLisHhQx6vXplP7dUvfoqnHZa9vPK/VcyFmN6d95pXRgTJxa7JOG65BLrpn3kkWKXJDw33ACv\nvw7f/Gbz5yorLU1Pv36ZW2MNDfb7sGiR3U+2117WyipW97Xj2ISKSZOydzVFRWMjLFtmyz9yTe3j\nV1trLfS+fYN/+Bs6FH78Y9vIvaWGDYNx4+DUU1v+GsXWty/88Y/2AaCiIv2Ynlp6JW7LFrj/fvtH\nGXUPPWSzU1980XYJiZrHHrPdbebMsf0OW6ptW/s+DRhQuLIVUkVFYqF6VINeXR289Vai6/gf/7Ae\ni61boXv3pjOo+/ZN3+LdutVmYc+ebbe337adVT74wLYU879O9+7Nr3cca+nlOwms3NfqrVljE3mC\nbMPWGi29lmZZAPg1cDawCUjVgI98S+/GG+2HGeXWj9/s2fDlL8OCBfl9Yi41r74KF1xg/yBbMq28\n3Myda93VS5dGq4tz61ZLizV5sv0cvXHWk0+2HYwaG60V7gWxWbMsQKb6kFNXZ5MvTjghEdwGDbKW\n+s6d9j30XmPOHGtFJqfjamyEbdvyX25w+eW2lOWKK/J7nWJ54gkb/pk0yR4Xs6WXT5YFgN8Av8AW\nqcfOunU2/ba6utglaT2DB1vX32WX2dq9KPzDXLrUAt4zz8Qj4IFNSKqttQCQ73hTqXjtNeuCP/dc\n+9vs2rX5OW3aWKvrmGPgmmsSiZZralKfe+ihto42WYcONinjlFPscUMDrFplG0snK8SOOuWeaSHo\neB6EH/T8WRYgkWXBH/RSZVk4AMuyMAvbzSWWbr/dPoH1iFluiXHjLBPDww/bP46wOI59Ug5zy6wt\nW+Ccc+BnPyvv8ZJceV2ckyenDnoNDYX5vufyOi39edfV2e/kk09aF/WoUcGvragozOSQtm0td2VY\nunWDDRts0ly5cRwLej/9abDzww56qTIoDAxwzkFY0IsNx7F+af8U808+sckKcdOuHTz9tHUZnXqq\nTewohN27Yf78xPd49mybMTloUKJ7aeDAws0erauzlDvnnZd64krUjR0LV15pyYY3bbJxL+97v2AB\n9O7ddMzq8MMzt+wdB1asaPo3snq1jW36X8frFm9osL8ff1fj1q1w4omJc086KXPG+KVLrZv2wAOt\nzNmWipSrww+33XQefbTYJWmZww6zVnMQYXcejcHG9Lye4q9hQc+f6vBF4G7gH+7jGcCN2FZkYC29\nF8kwprd8uX0SK1e7dtli3oaGpmuxjjkm3klWJ0ywae8nnpj/a23caN/jQw5pOg7TqVPT7dwWLrTB\n8KAbdGeyZo39k5w8OZ67kzQ2WiunY0fYvNkCjPe7ffzxlpjYH8B27rQPIJ07N3+t7dvhjTdsiUby\nJBFvQsns2RZYu3e39503z3pJ/Ofvu6+Nj3nvOW+e/cM88sjmP+/GRlvXeMcdcNVV0ehqj4tMY3ph\n/xgHAbdigQ9s78xGmk5meQSYiXV9AiwBhgIb3cdVZAh648aNY8cOyzzQr98w+vcfVsDit4527WyW\nWyH+0UaJ48DLL1uLN1/77GP/ULONf+zcacFx3brM5wVRWWldm4XORVhO5s+3gH/UUdm7Fdeutckb\ntbXNn9tjD5vw0bt35tdobIR337XZjyeckHrGo583C/P991M/f+KJFliltHlZFjy33XYbFCnoVQJL\ngdOB9cCbwAU0n8hyrft1EDbT07+XQxVZWnoiIiKeTC29sDtd6rGA9idgMTAJC3hXkci08AqwCpvw\nMh7wT12YCPwTOBwb9/tGyOUVEZEIK/fONLX0RESkiWK29EREREqGgp6IiMSGgp6IiMSGgp6IiMSG\ngp6IiMRG2EFvBLbYfDnwwzTnPOg+vxD4XI7XioiIBBZm0PMyLIwA+mOL0pN3UfRnWLgSy7AQ9NpY\n8e82EHWqazSprtFUbnUNM+j5MyzsJpFhwS9dhoUg18ZKuf1i5UN1jSbVNZrKra5hBr102ROCnNMj\nwLUiIiI5CTPoBd0qpdx3hREREWEQMN33+Ec0n5DyCPBV3+MlwP4BrwVYgAVX3XTTTTfddPNuCyiC\nSmAlliWhvVuIVBNZXnHvDwJez+FaERGRkjISSy20AmutQdMMC2CzNFdgSxaOy3KtiIiIiIiISGn5\nNZYRvtp3rBvwF2AZ8Gds2YbnR9ji/CXAmb7jx7uvsRx4wHd8Dyxn4XKsq/jgwhY/J72AvwHvAu8A\n33WPR7G+e2JLbhZgeSPvco9Hsa6etsB8LMEzRLeuq4FFWF3fdI9Fta57A89hOU8XAwOJbl2llQzB\ndp3xB717gRvd+z8E7nbv98f+ibbDxjZXkJjp+ia2lhFsTHSEe/8a4Ffu/fOx9Y3FcgBwrHu/M9ZV\n3Y/o1rej+7US+4MeTHTrCvB94GngBfdxVOv6PvaP3y+qdX0SuMy9Xwl0Jbp1lVZURdOg581YBQsU\nS9z7yTNWp2OTfQ7EPol5vorNgvXOGejerwQ2F6rQBfA8MJzo17cjMBf4LNGta09gBnAqiZZeVOv6\nPvCZpGNRrGtXYFWK45GrqzacLr79sS5P3K/eL1gPbFG+x79w3398HYmF+/7F/vXANpp/Si2GKqyF\n+wbRrW8b7JPvRhLdulGt6/3ADUCj71hU6+pgAX4ecIV7LIp17YMFod8AbwOPAp2IYF0V9EqLt8Yk\nSjoDk4HrgJqk56JU30asO7cncArWCvKLSl3PATZhY1zpNpaISl0BTsY+sI0Evo0NUfhFpa6V2Oz5\nX7lftwM3JZ0Tiboq6BXfRqzbAKxrYJN7fx02GcTTE/sEtc69n3zcu6a3e9/rk/+o8EUOrB0W8H6L\ndW9CtOsL9un1ZWwwP4p1PQnbM/d9YCJwGvbzjWJdAf7lft0M/BEbq4piXT90b3Pdx89hwW8DEaur\ngl7xvQB83b3/dRLB4QWsP7w91vVwGDZAvAH4BOsbrwAuBqameK2xwF9DLnsmFcDj2Cywn/uOR7G+\n+5KY1dYBOANrCUWxrjdj/+z6YHV4FStnFOvaEdjLvd8Jm6FYTTTrugHrejzcfTwc66J/kejVVVrR\nRGA9UIf9gn0D69OeQeopwTdjs6KWAGf5jntTgldg+Qg9ewB/IDEluCqEOgQ1GOvyW4AFgPnYLK4o\n1vdobBxkATa9/Qb3eBTr6jeUxOzNKNa1D/YzXYAtu/E2yYhiXQEGYC29hcAUrCUW1bqKiIiIiIiI\niIiIiIiIiIiIiIiIiIiIiIiIiPhdhS0CBrgU2znD8yiW3UJERCRy/oYtDhYRESk5VdiOFr/Dtmx7\nFtvK7HRsl5dF2HZu7d3z78a2h1qI5ToDuBX4b2AMtqn3EvfaPYGZJILgBe7rVZPIiwbwKfA/2M4j\nc4D9Clg/ERGR/6jCtmv7L/fx48CPgbVAX/fYk1iWim4k8pkBdHG/jsOSuoK19I7zneM97gGswXLD\ntcX2Pfyie04jcLZ7/x7glvyqJFK6tOG0SPF9gLWwwFp8p2EJPVe4x57E0hVtA3ZhgfE8YGea10tO\n+VMBnIC1+rYADVjW81Pc5+uwzBAAb6E9ESXCFPREis+fo6wC2ErTwOXdb8BS2zyH5bWbHuD10h2r\n8B3b7TveiKV9EYkkBT2R4usNDHLvX4hl6a4CDnWPXYy10jphu9xPw7ozB7jPV5AIjDUkuj09Dpb2\nZSiJ7s2vAq8VtBYiZUCf6ESKbymWlfvX2CSV+7DUK89if6NvAo9gefuexyaoVADfc6/3Z7R+wj13\nB5bw1bMBy4T9N/fal7BcadC0FRiJ7NgiIlKaqrDZlCLSCtS9KVJ8almJiIiIiIiIiIiIiIiIiIiI\niIiIiIiIiIiIlK//Bw+JtujWJRx6AAAAAElFTkSuQmCC\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": "iVBORw0KGgoAAAANSUhEUgAAAcMAAACgCAYAAABwkQeeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYE1XWx/FvCyLjrjiyK6CCKIqICoiOLS5s44aOguMy\n4oL7Liio4Iyyqiyj4i64MoIOoDCgKC3OKCgCDQg0+6uyC7KvTdf7x6mYdJPuTndSqaTy+zxPniSV\nquQWTffJvXXvOSAiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIeaQMsABYB3YrZZ4j7ei7QxN1WGZgG\nzALmAX0i9j8S+BxYCHwGHJ7wVouIiCRIBWAxUAfYHwtsDYvs0w4Y7z5uBkyNeO1A976iu72l+7w/\n0NV93A3om8hGi4iIJFILYELE80fdW6SXgWsini8AqhbZ50Dge+CkKPtUc5+LiIiU234evndN4OeI\n57+420rbp5b7uALWm1wDTMaGS8EC4Rr38Rr2DZ4iIiJl4mUwdGLcL6uY4/YCp2HB8U9AdjGfEevn\niIiIRFXRw/deAdSOeF4b6/mVtE8td1ukTcA4oCmQg/UGqwGrgerA2mgf3rhxYyc3N7ecTRcRkQDK\nxTpZ+/CyZzgdOAGbQFMJuzY4tsg+Y4Eb3MfNgY1YsDuK8CzRPwAXYUOmoWNudB/fCIyO9uG5ubk4\njhP4W8+ePX1vg85V56pz1bmmw/kCjYsLWF72DPOBu4GJ2PW/N4D5QBf39VewmaTtsFmn24Cb3Neq\nA8OxYL0f8A7whftaX+BD4GZgOXC1h+cgIiIZwMtgCPAf9xbplSLP745y3Bzg9GLecwNwYZztEhER\n+Z2Xw6SSBNnZ2X43IWl0rsGkcw2udDrfojM5g8Rxx4hFRETIysqCYuKeeoYiIhLVb7/B0KF+tyI5\nFAxFRCSqvn3hzjth1qzS9013GiYVEZF9rFgBp54K110Hv/4K773nd4viV9IwqYKhiIjs4/bb4ZBD\n4PHHoV49+OEHqFPH71bFR8FQRERitmgRtGgBeXlQpQp06wY7d8LgwX63LD4KhiIiErNOnaBRI+jR\nw56vXGnPFy2y4JiuFAxFRCQms2ZB27YW+A4+OLz9llvgmGPgySf9a1u8FAxFRCQm7drZ7e4iucEW\nLIDzzoNly+DAA6Mfm+r8XGfYBiu+uwirSh/NEPf1XKCJu602VsPwR2AucG/E/r2w6hcz3VubRDda\nRCQTTZkC8+fDbbft+9qJJ9p1xGHDkt6spPCyZ1gByMPyiK7AqtV3wpJ1h7TDcpO2A5oBg7HqFdXc\n2yzgYOAH4DIssPYEtgDPl/L56hmKiMTIceCcc2wW6fXXR9/nm2/stbw8qOh1ZmsP+NUzPAurRrEc\n2AOMwAJapEux6hQA07CyTVWxWoWhZZ5bsQBaM+K4IA/viogk3aefwqZNcO21xe9z9tlQowZ89FHy\n2pUsXgbDmsDPEc9/oXBAK26fWkX2qYMNn06L2HYPNqz6BuG6hyLik1WrYOpUv1sh5bV3L3TvDr17\nQ4UKJe/btSv07289ySDxMhjG+k9VtJcXedzBwCjgPqyHCDAUqItVK14FPBdHG0UkTtu3w5//bLdN\nm/xujZTHBx/YAvtLLil93/btYccO+PJL79uVTF6O+q7AJsKE1MZ6fiXtU8vdBrA/8BHwLoWr2a+N\nePw68ElxDejVq9fvj7Ozs9OqnIhIOnAcuPVWaNgQTjkFBgyAp5/2u1VSFrt323KJt96CrBguQO23\nHzzyiPUOL7jA+/bFIycnh5ycnJj29fLaW0VsAs0FwErgO0qeQNMcGOTeZ2HXEtcDDxR53+pYjxD3\ntTOBaKPcmkAj4rEBA+Bf/4Kvv4Z166BJE5g7F6pX97tlEqsXX7Trhf8pWoa9BLt2wXHH2XGnneZd\n2xLNz3WGbbEAVwG7vtcH6OK+Fqp4/wK2PGIbcBMwAzgHmALMJjxs+hgwAXgbGyJ1gGXu+62J8tkK\nhiIemjABOneGadOgtju+89BDNoT20kv+tk1is20bnHACjBtnX2TK4tlnYcYMeP99b9rmBS26F5GE\nWrQIWraEjz+26fgh69dDgwbw7bf2R1ZSW+/eMHs2jBhR9mM3b4a6dWH6dLtPBwqGIpIwmzdD8+Zw\n//3RF2f37g25uTZ8Kqlrwwb74vLNN+X/4vLoozaBasiQxLbNKwqGIpIQBQVw+eVQs2bxFdC3bYP6\n9WHMGDjjjOS2T2LXrRts3AivvFL6vsVZtQpOPhkWLoSjjkpc27yiYCgiCfHkk5CTA5MmQaVKxe/3\nyiswcqTtJ6lnxQqb/Ttnjn2xicett0KtWtCzZ2La5iUFQxGJ20cfwYMPwvffw9FHl7zvnj3WY3jx\nRbjoouS0T2J3++1w6KG2PCJeeXlw7rmwfHnqJ/BWMJSU4jixrWeS1DF7tq0pmzABmjaN7ZiRI6FP\nH5tgsZ/XJQEkZkUL9ybCFVfAhRfCXXcl5v284mfVChHAegrDhlnm+zvv9Ls1Uhbr19t1wsGDYw+E\nAFddZcmcP/zQu7ZJ2T35JDzwQGKL9HbtCs89B/n5iXvPZFMwFE/t2GFDZccfD++8Y2uTxo0LXiqn\noMrPh6uvtsBWUgLnaLKyoG9fePxxy3Ii/ps506753ndfYt+3RQu79pjOCbwVDMUTmzfb9Yh69WDi\nRJtm/8UXlr9y6FC76L5tm9+tlNI8/LBNlOnTp3zHt2plX4Reey2x7ZLy6dHDbpEV7BOlWzfo1y99\nE3grGEpCrV9vs8rq1bNvoRMnwtixti4tpH17+yb5xBP+tVNKN2wYjB9vGUZKq2RQkj59LF/p1q2l\n7yvemTLFqtVHWxuaCO3aWZq2L77w5v29lo6V7o8EPgcWAp+hEk4pYeVK60WccIJN2/72W8uEf+qp\n0fcfNMheV9mf1DRtml0HGj0ajjgivvdq0gTOPx+eL60ct3jGceCxx+Cpp0peEhOPyATeUlgFrLhv\nHawCxSygYZF92gHj3cfNgNCfxmpY/lGwMk55wInu8/5AV/dxN6BvMZ/viPeWLnWc2293nCOOcJz7\n7nOcn36K/dgRIxznpJMcZ+dO79onZbdihePUrOk4Y8Yk7j2XLHGcKlUcZ+3axL2nxG7sWMdp1Mhx\n8vO9/Zxdu+z/zowZ3n5OeVFCacF0rHQfecxw4PLEN11KM28e3HCDZRg58kgbfhk0KJywORZXX209\nyWee8a6dUja7dsGVV9o6tEsvTdz71qsHnTrpZ+2HUOHeZ56Jb7g7FpUq2UzVdOwdpmOl+6qEq1Ss\ncZ9LkkyfDh062LDXiSfCkiX2S1baIuxosrKsusHLL1suS/GX48Add9iswB49Ev/+jz9uM4qXLUv8\ne0vxylK4NxFuvRU+/zz9fs7pWOm+6L5pOncpfTiOXXxv3drWm513Hixdat82D4/zim2NGjbB4uab\n03uNUhC88IJ92Rk2zJukCFWrwj332Do3SY5Q4d6+fZOX6OLQQy0gpts14nSsdL8Gu6a4Giv0u7a4\nBqjSfXwcxzKOPPMMrF5tGerHjoUDDkjs53TubCVknn/eJm1I8k2ebD/nb7/1Ztp9yEMP2dB4bi40\nbuzd54h57TWrTPGnPyX3c++919Lx9ezpbwLvslS691JFYAk2zFmJ0ifQNCc8gSYLK+I7MMr79ic8\nM/VRNIEm4fLzHWfkSMdp0sQuur//vuPs2ePtZy5dahMs8vK8/RzZ19KljlO1quNMmpSczxsyxHHa\ntk3OZ2WyrVsdp1o1/yaz3HKL4/Tq5c9nF4c4RxIvpfzDqW2xmaCLsUr1YJXpu0Ts84L7ei5wurvt\nHKAAC6Az3Vsb97UjgUmUvrTC73/3tLN7t+MMG+Y4DRo4zlln2WzCvXuT9/mDBjnOuecm9zMz3dat\njnPqqfZvnyy7djlO3bqOM3ly8j4zEz3zjONcc41/n79ggeP88Y/2fyxVUEIwjGUU+T2gBXbt7k1s\n3WA6cM9dSrNjB7z1ls0AO+44uxbYqlXyk2nv3WvZ76+7TvlLk8Fx4Jpr4KCD4M03k/vzfu89Kwg7\ndaqStnthwwarKfntt+Uv3JsIHTrY35K77/avDZESUbXiMKAT8Dcssr4FfABsib95nlEwLMWWLZYa\nbeBAOPNMC4KRmWL8MH++BcQZM+CYY/xtS9D17m0FeL/6CipXTu5nFxTA6afb5I4OHZL72ZkgEYV7\nE2HqVFtSs2iRJW33WyKqVmzCeob/AmoAV2BDl/eWdJCkrunTLWdkcSnT/NKwoa1T6tIlfXMcpoNP\nP7VlLf/+d/IDIVi2kr597QuYZhEn1ooV8PrrqTFrt3lzW3s8apTfLSldLMHwMuDfQA42w/NM7Frg\nqcCDnrVMPDN3riXMfu21klOm+aVrV1i1ytakSeLNn28zeEeNsqUtfmnd2j7/rbf8a0MQ/f3vtlQp\n3gr2iZIuCbxjGSYdDrwBTIny2oXYZJZUpGHSKBYvhuxsGDDAhi9S1YwZ0LatFZWtqrQKCbNxIzRr\nZn+gOnf2uzXw3Xc2TLpwYepXSU8HocK9CxdaZqhUUFBgX7gHDoSLLvK3LfEOk65h30DYz71P1UAo\nUfz8s/1n7NkztQMh2PWkzp1T58J7EOzdazUJW7dOjUAIcNZZNpT2z3/63ZJgeOIJePDB1AmEkD4J\nvGPpGc4kXE0iZA5wSuKbk1DqGUZYu9YW3t52m/2ypIOdO21hdp8+mmSRCI8+aj2xiRNh//39bk1Y\nXh6cc47dp9If8XQzc6aVUVq82GYIp5Lduy0/7dix9kXXL+XtGd6BBb0G7n3othyYndAWiqd++w0u\nvhg6dkyfQAg2seONN6x3+NtvfrcmvY0YYQWWP/wwtQIhWIaUDh1sQo2UX/fulv811QIhhBN4Dxjg\nd0uKV1LP8DDgCCzDS7eIfbcA6z1uVyKoZ4gVVL3oIruO8Nxz6bmm65577Dw00aJ8Zs60L0OTJqVu\nCrSVK+GUU2DWrLJVPhHz1Vdw001WPcareoXx2rIF6ta10Yl69fxpQ3nXGR4KbAaqEH3V/oa4W+at\njA+GO3daVfl69eDVV9MzEIIFwkaNbM1U69Z+tya9rF1r1+UGDIC//MXv1pSse3dYs8ZGAyR2jgMt\nW1qiiuuu87s1JeveHTZvtqTwfihvMBwHtMeGRaNFlbrxNsxjGR0M9+yxunQHHmjZPryuY+a1zz6z\n651z5lg5Gindnj1w4YV2PS4d6ghu3GhZU3Jy4KST/G5N+vjkEwsys2al/u/56tX2s83Lgz/+Mfmf\nX95rhu3d+zpY4Ct6i0UbLH3bIsLJtYsa4r6eS+GJOm9iM1nnFNm/F1b9omjOUnHt3Qs33mjfGN95\nJ/V/QWJx8cVWQ7F7d79bkj7uv9/K6fzjH363JDaHH25LPvQzjl2ocG/v3unxe16tGlx1Fbz4ot8t\nKZ+WWF1BgOuB54FjYziuApaAuw62WL+0qhXNCFetADgXC45Fg2FPYlvs70ceWN8VFDjOrbc6zvnn\nO8727X63JrE2bHCcGjUc5+uv/W5J6nv1VUu4vnGj3y0pmx07HKd2bcf53//8bkl6ePttx2nRwn7v\n00Venn8JvCkhUXcs6wxfBrYDjd0gtBQrr1Sas7BguBzYA4zAstlEuhRb1A9Wyf5wrFYhwNdAcXMI\n0/Tql7ccBx5+2BaqjxkDf/iD3y1KrCOOsPVoN99s10Mluv/9zyrVjxkDhx3md2vKpnJleOopWwaS\nwVc5YuJH4d5EqF/f8g+/+abfLSkslmCYj5VTuhx4ESu5FMtVm5rAzxHPf3G3lXWfaO7BhlXfoPgS\nThnnH/+wGYPjxwf3ulqHDpbN4qmn/G5J6tm1y7J8XH65Vatv0MDvFpXPDTdY1YVx4/xuSWp79VXL\n5Zvswr2J0LWrzW5Ppby0sQTDLUB34DrgU2z4M5aVSrF+ryv6naa044Zi1yxPA1YBz8X4OYE2cCC8\n+65NNAn6wuUXXrBvlTNm+N2S1OA4MHKkTUz44gubgNKund+tKr8KFewa2GOP2TUx2dfWrTYpKh0m\nRkXTrBkce6z9v00VsRTVuAa4FugMrAaOAZ6N4bgVQOSKodpYz6+kfWq520qyNuLx68Anxe3Yq1ev\n3x9nZ2eTnZ1dylunp9dfh8GDYcqUzMjjWbWqLRXo3Bm+/z71FpEn0zffwEMP2bDxq6/CBRf43aLE\nuOQSS9/17rs2GUwKGzzYcgw3KZobLI1062bD+R07ejfMm5OTQ05OjjdvXgYVgSXYBJpKlD6BpjmF\nJ9DgHlt0Ak31iMcPAO8X8/nJvzrrgxEjbFLJwoV+tyS5Cgocp00bx3n6ab9b4o9Fixznyisdp1Yt\nxxk+3HH27vW7RYn33/86zjHH2KQaCfv1V8epUiX9f+cLChynUSPHmTgxeZ9JnBNorsSWPmzGhky3\nuI9Lkw/cDUwE5mG1EOcDXdwbWCBcik20eQWIrG/+AfANUB+7rniTu70flg4uFzgPC4gZ6dNP4b77\nLNekn9Ws/ZCVZYvwBw2CefP8bk3yrF9vSyaaN4emTa06wQ03WDLkoGnZ0jLmDB3qd0tSS79+tjwh\n3X/ns7JSK4F3LJ3TJcCfsUCWTtwvAsH05Zc2vDBunFWpz1QvvWRDaV9/nR7rrMpr5067VtqvH1x9\ntVUeOfpov1vlvblzbeh34cL0mxnrhRUrbALZ7NmpU68wHrt3w3HHwejR9uXOa/GWcFpN+gXCQJs6\n1QLhyJGZHQgBbr/drhn6ld7JawUFVoC5YUP4738t6L/4YmYEQrA0fO3apXaC52T6+9/hlluCEQgh\ntRJ4x9IzHIyt/RsN7Ha3OcDHXjUqQQLZM8zNtWwsw4ZZ8VuxXsPZZ/ubANgLU6bYulHHgWefhfPO\n87tF/vjpJ5soMncuVK9e+v5BFfp/nkqFexMhmQm8y5ubNGSYe180stxEagtcMMzLs5RkgwenftLl\nZOvf35aVfP55ei1AjmbhQptpN3OmLTHo2DGY1wTL4uGHYdu2zLp+WFBgtQmnT4cffrD1w9dfH8x0\ndT16WG5ar9O0xRsM01WgguH//Z9lbXjqKSvVIoXl59ukkjvusAw16WjdOhsGGzHCJhbce69lZBGb\nONSggS0lqV/f79YkXkEBLFliQS8U/GbMsKxLZ5xht6ZNoVWrYF4bT1YC73iDYQPgJWyo9GTgVCyN\n2tMJap9XAhMMV62yLBP33GN/ICW63Fyr0pCbCzVq+N2a2O3YYb39Z5+Fv/4VnngCjjrK71alnj59\nrLf84Yd+tyQ+jgNLlxYOfD/8YBOEmjYNB76mTTPr/0GXLpbI28vsUvEGwynAI1iO0ibuMXOxwJjK\nAhEM16+3xbUdO9pQgpTsySctGI4enfrDpQUF8P779nM94wzLMZnu0+W9tH27/fuMGWP/XunAcWD5\n8n0D30EHhQNfKPj5UdIolSxaZMtpli2zfx8vxBsMpwNnYOWSQvkOZmHp0FJZ2gfDzZutp3P++emX\njNcvu3bB6adbULzmGr9bU7zJk+062P77W47Gli39blF6eOUV6xlOmpR6vw+OY5N9QoEvFPwqVw4H\nvNB9JmSKKo+rrrJRMK9GwOINhv/BEmOPxILhVcDNQKrPZUzrYLh9u80WPekkW0uXar/4qWzqVEtW\nPXdu6g0zzZ9vk2PmzrUvOH/5i362ZbFnjy23+Oc/bVa1XxwHfvklHPBC9xUrFu7tNW2a2TNgy2ra\nNPsSu2hR+dIsOg6sWQM//mi/Yz/+CFdcEZ55H28wPA7LDnM2sBFYBvwVK82UytI2GO7ebX/Mq1SB\n4cM1k7A8HnwQ1q61BfmpYM0a6NULRo2yBNR33QUHHOB3q9LTqFE2y3b6dO9/N/LzbZgzL89m+ebl\n2S2U9Sgy8J1xRnpdq05V2dlw221w7bUl77dunQW70C0U/LKy7AvTySfb7eKL4fjj7ZjyBsOHijyv\njC3S344ts3g+hvPyU1oGw/x86NTJ7keOtG+aUnbbtlmmjiFDoH17/9qxfbtVFBk40BJO9+gRrDVi\nfnAcmzl8//32u5KI9/v113Cgiwx6y5ZZz65+fZvN2qCBPW7Y0Ba+q1efeOPH2/KRmTPt33fDhsJB\nLxT49uyxYBcZ+E4+2RJSFPdzKW8w7IUFvQbAmcBYd/slwHdYSafStAEGYWWfXsfyihY1BBty3Q78\nDbs2CfAm0B6rUnFKxP5HYnlOj8V6p1djPdai0i4YFhTYsoAVK+CTT9RziNeXX1oAmjvX21RejgOb\nNsHKlTbzd+XK8O3jj22hdO/elnZKEmPyZMvEMn++ZTGJxY4dtm6vaMDLy7M/nqFAFxn0jj8+eEWy\nU53j2BfZQw+1LyNbt4YDXWTgq1697F9G4h0m/RqrLrHFfX4IlmD73FKOqwDkARdiZZm+BzpROLVb\nOyyZdzugGZbtprn72rnAVuBtCgfD/sCv7n034Ajg0Sifn1bB0HHsm+706bZ43KvZVJnmtttsKO3l\nl8t+rONYdoxQYCsa6CK3VaxoQ2ShW/Xqdn/OOUqZ55U2bazU0113hbcVFNi1vGi9vDVrLNNJtKB3\n1FHq5aWSefNsbXWjRlCrVuJ+NvEGwzygMbDTfV4ZqxhRWh3tFkBPrHcI4YDVN2Kfl4HJWE8PYAGQ\njeVDBSvh9AmFg+ECrFrFGmztYw5wYpTPT6tg+PjjNjzw5Zdw+OF+tyY4Nm2yX6h33rFrESFbt0YP\nakW3gQ2HhYJb0WAXenzwwb6cXkabNcsmRtx8czjoLV5svz+RgS70+Nhjddkh05UUDGP5r/E2Niz6\nsfsmlwPDYziuJlZ6KeQXrPdX2j41CQfDaKpigRD33tNJygUFNqFl1y67j3xc3H1ZX1u3zn6Zp0xR\nIEy0ww6z2bjXXmt/EEOBrqAgenBr2rTw80MO8fsMpDinnWaTklavthmDDRrYOkT9zKQ8YgmGzwAT\nsGFLh8LX9UoSa7esaJQuS3euxGKNZbV4sQ1rRQaq/Hy7dlep0r730bbFsu9BB+27T+vWWnTrlUsu\nCQ87h4LfIYdoWCwIunQpfR+RWMQ6aPCDeyuLFUDtiOe1sZ5fSfvUcreVJDQ8uhqrer+2uB179er1\n++Ps7GyyI8fJojj2WJvBFBmkKlbUH80gaNXK7xaISLLl5OSQk5MT075e/pmviF1vvABYiQ21ljSB\npjk287R5xOt12PeaYX9gPTYz9VHgcAIwgUZERLwVb3Hf8srHAt1EYB42SWY+0MW9gc1KXQosxhb2\n3xlx/AfAN0B97LpiqFZDX+AiYCHQisITckRERMosyAOA6hmKiMjv/OoZioiIpAUFQxERyXgKhiIi\nkvEUDEVEJOMpGIqISMZTMBQRkYynYCgiIhlPwVBERDKegqGIiGQ8r4NhG6z+4CKsEG80Q9zXc4Em\nMRzbC0v4PdO9tUFERCQOXpa6rAC8QOFK92PZN1H38cAJWK3DoVii7pKOdYDn3ZuIiEjcvOwZnoUl\n4F4O7AFGAJcV2edSwoWCp2EVKKrFcGyQc6qKiEiSeRkMi6tiH8s+NUo59h5sWPUNLICKiIiUm5fB\nsLyV7kszFKgLnAasAp4r4/EiIiKFeHnNsLyV7n8B9i/h2MjK9q9jxX+jKmulexERCY4gVLov6djq\nWI8Q4AHgTODaKJ+veoYiIvK7kuoZetkzjKx0XwG7vheqdA9W2X48FggXA9sIV7Mv7liAftgQqQMs\ni3g/ERGRcgnyrEz1DEVE5HeqdC8iIlICBUMREcl4CoYiIpLxFAxFRCTjKRiKiEjGUzAUEZGMp2Ao\nIiIZT8FQREQynoKhiIhkPAVDERHJeF4HwzbAAmAR0K2YfYa4r+cCTWI49kjgc2Ah8BmqZygiInHy\nMhhWAF7AgtpJWNWJhkX2aQccD5wA3IbVKizt2EexYFgf+MJ9nrFiLU8SBDrXYNK5Blc6na+XwfAs\nrBrFcmAPMAK4rMg+lwLD3cfTsF5etVKOjTxmOHC5F41PF+n0ny1eOtdg0rkGVzqdr5fBsCbwc8Tz\nX9xtsexTo4RjqwJr3Mdr3OciIiLl5mUwjLV+UixlpLKKeT+nDJ8jIiKSdM2BCRHPH2PfSTQvAx0j\nni/AenolHbsAG0oFq3q/oJjPn0U4WOqmm2666abbLHxQEVgC1AEquY2INoFmvPu4OTA1hmP7Ew6M\njwJ9E95yERGRBGoL5GGTYR5zt3VxbyEvuK/nAqeXcizY0opJaGmFiIiIiIiIBMmb2MzYORHbSkou\n8BiWjGABcHHE9qbueywCBkdsPwD4l7t9KnBsYptfJrWBycCPwFzgXnd7EM+3MrZkaBYwD+jjbg/i\nuYZUAGYCn7jPg3quy4HZ2Ll+524L6rkeDowC5mP/j5sR3HMVn52LZd+JDIb9ga7u426Er42ehP1x\n3R+7prqY8Izc77A1mmDXYtu4j+8EXnIfX4Ot2/RLNeA09/HB2FB4Q4J7vge69xWxX/RzCO65AjwI\nvAeMdZ8H9VyXYQEhUlDPdTjQ2X1cETiM4J6rpIA6FA6GoZm1YAEkNGu26KzcCdjko+rYN7eQjths\n3dA+zdzHFYF1iWp0AowGLiT453sg8D1wMsE911rY9fzzCfcMg3quy4AqRbYF8VwPA5ZG2R64c1Wi\n7tRVXHKBGlgSgpDIRAWR21cQTlQQmdwgH9jEvt9q/VAH6xFPI7jnux/2TXkN4eHhoJ7rQOARoCBi\nW1DP1cEC/3TgVndbEM+1Lhac3gJmAK8BBxHAc1UwTA+hNTJBcjDwEXAfsKXIa0E63wJsWLgW8Ces\n1xQpKOf6Z2Atdg2tuEQaQTlXgJbYF7m2wF3YpY5IQTnXitgs/5fc+23smw86EOeqYJi61lA4ucBa\n9/EKbBJKSC3sG9cK93HR7aFjjnEfh8b8NyS+yTHbHwuE72DDpBDs8wX7tjsOm0QQxHM9G8sbvAz4\nAGiF/XyDeK4Aq9z7dcC/sWthQTzXX9zb9+7zUVhQXE3AzlXBMHWNBW50H99IOGiMxcbbK2FDGCdg\nF6ZXA5uxsfcs4HpgTJT3ugqr9uGXLOANbFbaoIjtQTzfowjPsvsDcBHWcwriuXbH/gjWxc7hS6yd\nQTzXA4F/AllmAAACZklEQVRD3McHYTMm5xDMc12NDWHWd59fiA31f0LwzlVSwAfASmA39h/vJkpO\nLtAdm6W1AGgdsT00dXkxVicy5ADgQ8JTl+t4cA6xOgcbOpyFBYaZ2KyyIJ7vKdh1llnYNPxH3O1B\nPNdI5xGeTRrEc62L/UxnYcuDQklBgniuAI2xnmEu8DHWcwvquYqIiIiIiIiIiIiIiIiIiIiIiIiI\niIiIiIiIBEkXbHEzwN+wTCEhr2HVQkRERDLGZGzRs4iISNqog2XweBdLXTcSS+l2AZbVZjaW1q6S\nu39fLE1WLlZrDqAX8BBwJZYMfYF7bGUgh3Bw7OS+3xzCdekAtgJPY5lWvgWOTuD5iYiIlKoOlrau\nhfv8DeBx4CfgeHfbcKzqx5GE68kBHOre98SK7YL1DE+P2Cf0vAbwf1htvgpYXsjL3H0KgPbu435A\nj/hOSSR1KVG3SOr6GeuRgfUQW2GFVhe724ZjZaE2ATuxgHkFsKOY9ytaWikLOBPrJa4H9mJV6v/k\nvr4bq7QB8APKGSkBpmAokroia8RlARspHNBCj/diJYRGYXUFJ8TwfsVty4rYtidiewFWXkckkBQM\nRVLXMUBz9/G1WFX1OsBx7rbrsV7dQVjVgP9gw6KN3dezCAfMLYSHT0McrLzOeYSHSTsCXyX0LETS\ngL7piaSuPKyK+pvY5JjnsRI3I7Hf3e+Al7G6iaOxiTFZwAPu8ZEVyIe5+27HCvGGrMYql092j/0U\nq1UHhXuNgahmLiIi6aUONrtTRJJAw6QiqUs9MRERERERERERERERERERERERERERERERERHx2v8D\nCdQE/Bo1/cIAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 11 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 12 } ], "metadata": {} } ] }