{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Dump and reload example\n",
    "This notebook illustrates how to dump a cancer simulation to disk, reload it, and continue with changed parameters. This feature allows to model different stages of tumour growth characterized by varying growth parameters (e.g. `death_probability` and `adv_mutant_division_probability`)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Import modules "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# The cancer simulation module.\n",
    "from casim import casim"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  3rd party modules.\n",
    "import os, shutil\n",
    "from wand.image import Image as WImage\n",
    "from glob import glob"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pickle\n",
    "from matplotlib import pyplot"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Define some utility functions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def unpickle(sim, pkl):\n",
    "\n",
    "    with open(os.path.join(sim._CancerSimulator__simdir, pkl),'rb') as fp:\n",
    "        return pickle.load(fp)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def clear_outdir(dir):\n",
    "    if os.path.isdir(dir):\n",
    "        shutil.rmtree(dir)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup parameters \n",
    "Initially, cells with advantageous mutations divide only seldomly, the tumour is dormant."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "parameters=casim.CancerSimulatorParameters(\n",
    "                                matrix_size=200,\n",
    "                                number_of_generations=50,\n",
    "                                division_probability=0.6,\n",
    "                                adv_mutant_division_probability=0.5,\n",
    "                                death_probability=0.3,\n",
    "                                adv_mutant_death_probability=0.05,\n",
    "                                mutation_probability=0.5,\n",
    "                                adv_mutant_mutation_probability=1.0,\n",
    "                                number_of_mutations_per_division=10,\n",
    "                                adv_mutation_wait_time=40,\n",
    "                                number_of_initial_mutations=2,\n",
    "                                sampling_fraction=0.9,\n",
    "                                plot_tumour_growth=True,\n",
    "                                export_tumour=True\n",
    "                               )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup the simulation engine. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "outdir = 'out'\n",
    "clear_outdir(outdir)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2020-08-19 12:29:15,983 INFO: Running in single tumour mode.\n",
      "2020-08-19 12:29:15,985 INFO: First cell at (100, 100).\n"
     ]
    }
   ],
   "source": [
    "cancer_sim = casim.CancerSimulator(parameters, seed=1, outdir='out/')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run the simulation "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2020-08-19 12:29:16,398 INFO: Ready to start CancerSim run with these parameters:\n",
      "2020-08-19 12:29:16,400 INFO: matrix_size = 200\n",
      "2020-08-19 12:29:16,401 INFO: number_of_generations = 50\n",
      "2020-08-19 12:29:16,403 INFO: division_probability = 0.6\n",
      "2020-08-19 12:29:16,404 INFO: adv_mutant_division_probability = 0.5\n",
      "2020-08-19 12:29:16,405 INFO: death_probability = 0.3\n",
      "2020-08-19 12:29:16,406 INFO: adv_mutant_death_probability = 0.05\n",
      "2020-08-19 12:29:16,407 INFO: mutation_probability = 0.5\n",
      "2020-08-19 12:29:16,409 INFO: adv_mutant_mutation_probability = 1.0\n",
      "2020-08-19 12:29:16,410 INFO: number_of_mutations_per_division = 10\n",
      "2020-08-19 12:29:16,411 INFO: adv_mutation_wait_time = 40\n",
      "2020-08-19 12:29:16,412 INFO: number_of_initial_mutations = 2\n",
      "2020-08-19 12:29:16,413 INFO: tumour_multiplicity = single\n",
      "2020-08-19 12:29:16,414 INFO: read_depth = 100\n",
      "2020-08-19 12:29:16,415 INFO: sampling_fraction = 0.9\n",
      "2020-08-19 12:29:16,416 INFO: plot_tumour_growth = True\n",
      "2020-08-19 12:29:16,418 INFO: export_tumour = True\n",
      "2020-08-19 12:29:16,419 INFO: Tumour growth in progress.\n",
      "2020-08-19 12:29:16,621 INFO: New beneficial mutation: 2076\n",
      "2020-08-19 12:29:16,653 INFO: Division of beneficial mutation carrier. Cell index = (108, 110), mutation index = 2076, place_to_divide=(109, 111)\n",
      "2020-08-19 12:29:16,675 INFO: Division of beneficial mutation carrier. Cell index = (108, 110), mutation index = 2076, place_to_divide=(108, 111)\n",
      "2020-08-19 12:29:16,688 INFO: Division of beneficial mutation carrier. Cell index = (108, 111), mutation index = 2076, place_to_divide=(109, 110)\n",
      "2020-08-19 12:29:16,738 INFO: Division of beneficial mutation carrier. Cell index = (109, 111), mutation index = 2076, place_to_divide=(110, 111)\n",
      "2020-08-19 12:29:16,741 INFO: Division of beneficial mutation carrier. Cell index = (108, 111), mutation index = 2076, place_to_divide=(108, 112)\n",
      "2020-08-19 12:29:16,757 INFO: Division of beneficial mutation carrier. Cell index = (110, 111), mutation index = 2076, place_to_divide=(110, 112)\n",
      "2020-08-19 12:29:16,771 INFO: Division of beneficial mutation carrier. Cell index = (108, 111), mutation index = 2076, place_to_divide=(108, 112)\n",
      "2020-08-19 12:29:16,772 INFO: Mutation  of beneficial mutation carrier. Cell index = (108, 111), mutation index = 2076, place to divide = (108, 112)\n",
      "2020-08-19 12:29:16,789 INFO: Division of beneficial mutation carrier. Cell index = (110, 111), mutation index = 2076, place_to_divide=(111, 112)\n",
      "2020-08-19 12:29:16,790 INFO: Mutation  of beneficial mutation carrier. Cell index = (110, 111), mutation index = 2076, place to divide = (111, 112)\n",
      "2020-08-19 12:29:16,792 INFO: Division of beneficial mutation carrier. Cell index = (110, 112), mutation index = 2076, place_to_divide=(111, 113)\n",
      "2020-08-19 12:29:16,793 INFO: Mutation  of beneficial mutation carrier. Cell index = (110, 112), mutation index = 2076, place to divide = (111, 113)\n",
      "2020-08-19 12:29:16,878 INFO: All generations finished. Starting tumour reconstruction.\n",
      "2020-08-19 12:29:16,884 INFO: Reconstruction done,  get statistics.\n",
      "2020-08-19 12:29:17,206 INFO: Exporting simulation data\n",
      "2020-08-19 12:29:17,443 INFO: Growth curve graph written to out/cancer_1/simOutput/growthCurve.pdf.\n",
      "2020-08-19 12:29:18,021 INFO: CancerSim run has finished.\n",
      "2020-08-19 12:29:18,023 INFO: Simulation output written to: out/cancer_1/simOutput.\n",
      "2020-08-19 12:29:18,023 INFO: Log files written to: out/cancer_1/log.\n",
      "2020-08-19 12:29:18,024 INFO: Consumed Wall time of this run: 1.601802 s.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cancer_sim.run()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualize the first stage of the run"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "pdfs = glob(os.path.join(cancer_sim.outdir,'cancer_1', 'simOutput',\"*.pdf\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "out/cancer_1/simOutput/sampleHistogram_97_107.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAJlQTFRF////ODg4AAAAiYmJ7u7uiIiId3d3mZmZ3d3dREREu7u7MzMzERERZmZmVVVVzMzMqqqqIiIi4e31w9vr0uTwiLbXH3e0W5vIICAgO4i9S5HD8Pb6l8DcS5LDT5TEpcnhXp3JQ43AUJTEV1dXP1VkDjhUHD9YKkhdW2ZuoqKiQEBAGBsdV2NsVmNsFxsdXWduYmpwR0dHwMDAuq7KgQAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHRdC9SyJAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAAC9tJREFUeNrtnQuDmzYWRhUtCCQERnUcN9ntJt02TfaR7uP//7iVBAOe+jFcfMEo8506lDEg2RwjJHERQgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA7eSMDJn3j1ZLkqMiFKpbL+HaUe/ZP5rtCGN72qUrYWpW207Y1BGCvcwqQRuRWqEsIW3TsQxgq3MK3zSonGS9K9KAhjhVuYstrmUVYUpvxZEsI4YRa2k6VopahrHGELwSys9OewQorCZZlsu7cgjBXuIrGR2peBWeV8XbEDwljhFiZ2pgz/M7unNyCMFXZhZ7wg7If9/u2jd0JKPFzY4d2746N3QkpAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJAWGKwh7kFKKHaEEaDP5BU24oSqg1hNBYoEm1BCdWGMBr8wgp/aBFCtSGMBr8w3Uf9RmGF1j9CGCfswkLsLyVUG8JosAtTzk8IodoQRoNdmMv9hBCqDWE02IWZrjY/OVQbwmigpyMxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwxICwx+KOmlDII1V4OBmFlLgrd9H5aWasGodrLwSAsr4VTVR+Y7aIehGovBoMwpUqZ9VpKqWW1G0O1EfnLDYOw1lZa1HmcN7IpazeGagscYcxwVDoKVQoVB44VmTReGkK1l4O7lujPYbRRtSGMBs8R5umOMF8+apkjVHs5GIQZX5MfhImMOKo2hNFgqSXesxjCaDAI27ny1mIIY4WlSAxcTQfCWGEQlplAdm0xhLHCUUsslVbXk4EwVjg6f61Wjbx6HoMwVjg6fxs/qa9qgTBWOIRpP2nya4shjBWOSod1urIoEteBp9IxdnScA2Gs3C1MSbTD1uRuYaVBO2xNOLqmQixAsbu2GMJY4eiaCkYaFInrcL8w2YFa4jrcL8wU6tYpDMJ44ajWZ6h0rAdLw9k6K/tKR9G1yRD5uxQsXVNGZVWfjq603iHydzlYQgSMHrR06SHydzE4AkldJvWTHm1tjUGaF4Sj0mHETvW+RGtaV2OQ5uVYoKcj14j8XQ7ung5fQ2w0In+Xg7unQ1bW7hD5uxzcPR2l6RJE5O9C4B7nxOCImlLq1iVnCGOFSRhiOtaCoy8xTOri2mIIY4XrHNYgLnEdOBrOWmuJInEluG6GwAXMleAQphHmth4c18MqXHFeDwZhRX1zMYSxwlEkuhoN59VguWXWVxM1AknXgSNEAEXiinAcYRhFYEUwikBiYBSBxMD1sMSAsMTgF2YQqr0kdwsrivJZEpnz9Q9aqPZ7OJvO/fc4N616/rcXRgvVhjAC99/j7ORptd7YMIosbZBmCCPAFUjak7lWeGG0QZohjADz4GC51GEcdFqoNoQRYB4czCilZF3SQrUhjMACg4P5IpEWqg1hBBYYHCz2UlFCtSGMwBYGB4MwAlsYHAzCCGyhLxHCCEBYYkBYYnC0w3Y3F0MYK1sIwoEwAlt4lAeEEdhCEA6EEdhCEA6EEWAZCadR4p5QbS/sw+FwePSuSAOOWqJU+q4HDbyPr3eP3hVpwPPAN31DC4SxwnG7kfbCNI6wdWDpra+ss/dVOiBsKixj/uZjFOI5EMYKh7C2ktXVEhHCeGFpOCuT2z6dVsWxLkmRvxBGgOuxwE9anFJW0SN/IWwqTDf0ldWYjm+VkSN/IWwqDI+jehb5Kwplc+IgzRBGgOFxVOZZX6LSsiAO0gxhBBa44pxb4iDNEEaAe+ii0htyxEGaIYwA99BF0jlf2yBH/kLYVLiHLsp6dcTIXwibyhaGLoIwAlsYugjCCGwkagrCprKFoYsgjMBGoqYgbCobiZqCsKls4ckQEEZgC0+GgDACW3gyBIQR2MKTISCMwBaeDAFhBLbwZAgII7CROzAhbCr3X3FWCtX6FWES5tDTsRIsRWLh7NVaPYTxwvJoe3sjUhvCeLm/SNQ3dUEYM3cLM8/jElsVh05EqPZS3C0sex6XaOtaFvNCteFsCtztMG+pnhmqDWFTWKDh7NQYqm2UegNhnPALq31ZSB6kGcKmwi6stjsxL1QbwqbALSz6mheqDWFTYBYWK/l6Xqg2hE2BWVis5IdjbEaoNoRNYUOXVyBsChCWGBCWGBCWGBCWGBCWGBCWGBCWGNsSBmcvAmGJAWGJAWGJAWGJAWGJAWGJAWGJAWGJwS2sbOONR3MifyFsCtzCdLzxaG7k7yBsv9//+dG7ZpvwF4lB2NzI30GYn9k/etdsk2WEzRqkGcKmsIywIfI3jLkNYZwsI2xu5C+EvcgywuZG/kLYi7DXEu+K/IWwF+EWtrsr8hfCXmRzPR1dC8zP/GW//+HRe2eDbE5Yd2x1M3jw7DkQlhgQlhgQlhgQlhhbFDbghR1va3t91X8ISwwISwwISwwIS4w0hB2Ox8va/Eo/HY9/PR4/Pmb3rU8awvz0eHHr01VfBxCWGCkJO8TpvlvWbT1f2PlvYJ/ECRHCBiCs4w5hnw6Hn56mHw6HTzzCPh4Oewi7zh3C3p+9OIRdOy2+amHzQ7Uh7DbLCLsjVPu2sJ+Pxw/dsmNsfg2r/i02xXx77VNssn281nATUZVf++cTYW99YoMwP//25ifu1n4Uywi7I1T7trDL63WrdofN+5P6yWWOJ4drxz4esb2wF4+z/btH3rKxjLDTUO1f3qhb/Pr582/xdYHfLr0+X1v1V6W+xJkvMdUv13L88pTSsMZXn0A3/apUN71Bt/aj+PHvi/wMTkbV/scvz3P859ev/7r0+vb0GufOV/r94nrfJib5+/Ssv93M+vKn/PftrL/dzPrqF/+Pl3Sy+/77vyWEnYZqn8m865dw3+/oVWY9hdNQ7Q199NeZ9RROQ7X/iCKlxLnxK816GmOoNgAAgAsUqu1mTjquJjNs0z+hbGbW/oRQUDceP245I+8ha9/EIZ8nTF78MZn1UNb/F7/2ScfVVMZtuieUzcxaiFxS61xj1q2sVTM361aqRpa0jY10+uwbrIav5+cuZn7ScTX5ew/bdE8om5l13PvEjces3YxdNn5rn6/NqZsr/cdvsBo7mflfTJg77biayrNtqDtuzNqnoajChqxLqWVFLNXGrI3MlSUeYYOwk2+wGiG/Lk89Q9jpNjW1PB2zzitBFjb2s8mmrIk/8zHrnauqaq6wMZn1KH0B3uV5q+PqGifbdE+8mpe1ddpRzwVD1pk05L02Zt34ZCpymdoLG5NZEZeH00FmbnZcXaPfxmQzfD1lXe6E8iWiI54+x6x9UVxQ99qQdTjx0k+CUVjIOiazLoWs/L72v5NbHVfX6LeRpn9C2Zysu18ruUgcshat1ZJabRiyNtY5R64lxgcu+1eXzMqU4dnBWWjIzOm46gYnyIYnlNGzLuNmJT3vPuvw4cknoTFr/8mp22bdA5djCobcdAUAAADAKhShx9YM9ehdeTp3oXqtNLlXGHCiQxdxNTSbxpZ6mDsXtpMtvQYIGNHOt1elDHKysiidVqVptMrinHdZ1tq3qVTe1HFU6dr5d/NaZbkObyitCuFXC3Xx0v+RCb+wCSMLNrowfvPd2n3n3z1a2VLVXljohtNZVZusKNqqjnNaZLZprRKyaqtwxSTLK/+uzdu6amsnclv4TY2OcbCuMVUttH/LisKqVhmZxQ4rwIlWtXLtk7BYEJa5aro5LcJOL1y4bNHGro8gJyyRuTHSxLlemP/T5DIszGRfslZ56FcErGi/p12w5c20naaqUfUoLPZ/h/7c58JqpVQZKiBOmKoTFgI4u4W9sMKtfXHqFeB3bZNHYa3oNXk7wxG2k7GH70xY6MsvRdP4FeJR6NQujsTfCwtduJkvT+mXIsELdMdCqHTYqtKh93SnbeX6uXCVprLtubDWry3Fzi+0/pB0rlJCSe1646L074SNJXr2uOkaXl1Xc9b33g5zYRq7U/2/+Eecxk0y03fz6jgty7Cm6RbG1MJiVDk2yI3LMBX9Wj9YnBuXYWZcZQEAAAAA+G75P+KdhneEegN5AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE5VDEwOjI5OjIzKzAwOjAwFVvQdwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xOVQxMDoyOToyMyswMDowMGQGaMsAAAAgdEVYdHBkZjpIaVJlc0JvdW5kaW5nQm94ADQzMngyODgrMCswENDwCAAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNCQxalcAAABKdEVYdHNpZ25hdHVyZQBhZTJkOGQ2NzgyOGQwYmE0YWU3ZTBlYjVhYmJiOGYzYTU5ZTU1YmE0YjgzNjI4NTQxMmMzMmU2Nzg4MmE5MDEzN+iLsQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<wand.image.Image: ae2d8d6 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "out/cancer_1/simOutput/growthCurve.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAASNQTFRF////ODg4AAAAiYmJd3d3REREu7u78Pb6pcnhH3e0lsDc7u7uiIiIzMzMmZmZqqqq3d3diLbXERERMzMzZmZmVVVVIiIiW5vI0uTwc6nQS5HDYZ/Kw9vrTJLDaqTNUJTETZPDx93sPYm+tNLmbKXO4e31osfggbPVTZPEl8Dcb6jPQozAea3SbKbOjLrZxt3sS5LDc6rQk77bLoC5yN7tV5nHVZjGdavRqszjoMbfR4/B1OXxZqPMSZDCXp3JQIu/Z6PMSJDCnMTeVpjGdKrQlL/brc7kO4i9aKLMZ6LMdKrRICAgutXoeq7SKn64nMPeXZzJZaHLmsLeU5fGhrbXfrHUdKvRYJ7JyN7seK7SpcjhgbLVW5zIV1dXoqKiR0dHwMDASvnegwAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHRdC9SyJAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAAFeNJREFUeNrtnYtj47hxxnEMH4KoFymJ1OooqUsdb2tXiWUzvsvlNtekj03a9JI2bfpO2///r+jMgNLJsmSJEGmT4vzW64ckmmt/C2AwGHwQgmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhmEYhimCzwymSH5UtmCm+db/Za4Kyy77DixYobBgNYMFqxksWM1gwWoGC1YzWLCawYJVGaf17CEWrMq0ZHv/oQIEc20Tv4nbMbvwoWf24f3ANN3saRZMH096+w8VIJjp+6jJcDgybPgC/oiBP7b8TDEWTJ9yBINvAprYhivGI2H0RScQ5lAIv6eeZMH0KVOwDmhkWl2QzTbE2FQPIiyYPp4M9x8qTDDTwr8gFgpmbQQzDcNgwbQJSxSsF6BgA2OAgo1G3MIuw5ng+zIFGxhdMTRF0MEBrBe4rtFXT7JgOngSFXs3/Xz/iSKiRFwHFWLkB0NX9Iyh3xXuMPBHm6dZMA08GcH7mZT7TxQg2MAG4GPXVl9ROG93N0+zYDrMJSY5ZCmCnYAF0yGkgJ4Fqw1bwaK9J1iwaqIEW/zZ4v3eEyxYNYkxoI/kgltYTZAbwZZ7T7Bg1SQT7ItnyUQWrJpIOWPB6gElpZLFYiHEl4sPLFjlmePUK6IZmCdbcr73NAtWNUIMDN8v/hxa2M3i/eJ272kWrGqQYJH8CznBhOKzdD0LVjVIsCXNwFZywoJVHomVUp78MQgWSjHbTyayYFVDYmC4FSxkwapOJthP4AMLVgekXKFgNyAY9IdzmTx9mgWrDhOcMt+t1/dC3K/v1g8CPvP2s78sWHUIU4EzMJx63cIM7IOAz27211dYsOqQZjmOGcb2AgJ6eOMWVmHWa0eInz5+9SjE14/i8Wt8+9n0m6cvYsGqA9UDePLn0NCglckwgrgj2s/+smCVwckEW6NgoYhTFqzaRJlgOGUGwULJglWbuymOVw/Tr6d3YvpOvJveTR+orT2BBasCITajiHIcoYTAMKI0h4cFHSxYFcFIXrQywVqyjYKpvpEFqyRU4OtRUoqGrh3B0vjpS1mwCvDtdCpw9IKRSw1dOHqpwexZ9pcFqwAwejnUwkJsURM5b8uWaMsUE78sWAVRAX1ISSnKSGFCStXhFCWY6XYtq/vkoYFJhg9s+5CfJQl2u/i4EJjwXdxiyvc9VbqJn+xlf/UE6wZiZI1Guw8N/LFp9Nj2QQfoDJe0P1blONJ4p4XtZ3/1BLOHIujb1pOHDCHGY7Z90MGjgD5GwSYQKYYSly0TKVNRlGCuAQ2oN37ykD/q+H22fdBBCSbDlZyoKTMNXBSDCG9vO4TmGIZuRb0nl7pWYPn2vu2DbZqfsWAvMMedsatMsG2OY/aDYPvJxMKixA7bPmigCkVD0mYrmKdGr41gT/3BdASzMwf1J0r0IMgYWWz7kI8dwWD0gukXzsA8NXodMOkQeoK5tmLw5MHAH/p9tn3IRyQx83S7WHyEMP4G3jCgv6GyDnj05sAVBU6cbTJ8YNuHPLSo80tjiDJUChGDwqXM2p134IrCusSjsGAv4KnRKgTBllvBIiXYvDDBDnaJR2HBXgAEmyjBHJDKkSHOwECwuTjovSe0u8QBxBe98y5lwV7Aw6SUQ/GhR5UBOANLKMp/NgNTaAoWjAPRGZ/1UhbsMJldG3WB3hPBVHm9aKXJgcu0U1OWeJqaOgoLdpAJza9C3FxEgi1X0DnG6Y5gh9ETbOC7ljC5hV2AymDc/gJidxXQ3y4wkv/uER59PBjQKzS7RAtmXQaPYReg3GHlX4JuqogD21a4Ceijo9fpzsN6ptk975Us2EEywR5V/Xwm2Lw8wc6HBTvIw/qX8H59v76HT+/wkx/jg7+6F2UI5hoDUGJ01mtZsIOEkgL4uZwJFR9i2/I2M7DJ0es0o0Rr8+40LNhBSDDVE8KsWchZJhiGIs5xvXSjRGphHCVeQIg5jo1g9OWciju8E9fpRom+xVGiJkvq8FJJaV4QLMkE8zbV2i9yQZQ4OO+VLNgeKkc4/cX0TjxM7+Dt8ynWj+LOh7vp9K9OXM1R4qujBJM/p7XlKLMnUqNXIl+IDxUs2KuzwhL6ifxcCdaWrUww5XtejmBujqtYsD1CZTb6jgraEog7MsFQqrIEs4PzX8uC7REqs1GPdllmymXV2qUJJvzzL2PB9ohxBqZ6QhAskX+dxYco1TOnomdoCsYlAvpQzfzN4v13jwJz84uvMDf/ngrpyxOMSwS0oBVJKZV3pTJPwS/VIRD43GR56lvo7l5BeOKckzZOmR06n2MuJ6F0ULAZ7YQQ8mTbUugLNt5Uip586Vv/nirDXLnDokLQvDy5RMHCLNxIz/semmE9vutxtj4nIZYF3E3/Zvog3k3Fw/RvMbvxboo7Y8XyZGeouGDi3OFsfU6oNHRJpwfEKXSQH5SX5QuLKc/RTf5als9dYl5SJdAnZdSWnfvgnTt6KS6IEgfnvZQF2yJVTiPaCPaI6ajXEcztWKPeeS9lwbZsBEvjRM5hyvzrLAjJ8z00BRsG5sg/TzEWDFhm3kQpDmSO8v3Kzt9TG1jORndTOsSJvb2gw6YlMnYROARWAWTbzGmJeSPYBPP2YZ7vpCkYWgTslWqbvjky2UXgMOFmBpaoGdhK2UpRZ/gagmGJgPU0ShyoigF2ETgIdX4P09+gs97fwSd/j3Ovd+RYJJZ5onr9oGM/M9UzLGPssovAYXaOe4AG1Yb4MKKy0fzfqbAVZ9MwB8PRvotAz7J+y4Jh/nB73AMF9AtlPpqrM1To5hKfLa/0DXYROEpEznq3i39Y3Hy7+LBxJbp5djjYGegGHf39mTOOYewicAS1CTaUE7naxIcCk1Sv1sIGB0oETH/o2+wicJBP0+kXQnwxFdN3ZOQ7pXAjwmaXF80ucWw9Xw8bsIvAMZRVCmakZm3M2NOx2iJJz8zQ76K5vOIPLcviTMcZzNHOdyvYTHrKGlGjM1Rw1VTJ0Gg1l9+TR4DyzH4DwVyumjoXqXxgMxuplfwlpqO8kyX0R+GqqXJJMnfYrbf5WrlWvrZgXDV1JjDl+iTE4uP7xc2nxZcw9foKZ2DLPYu2HHDVVLlEmRUizMB2fGGdlvY35Kqp0kgckZ0e4OxaZWuHGwqumiqNVG1XptyhBzMwdE5J3kgwgqumXsShleSV/F0snhn5XgBXTZVFlJ0ToJz1NoKl+vGhgqumyiLKZmBzmeDo5cnfKV/YtxEMRehyauolfr9e3wnxj//0h8efwpv4w+M/fy1wlfnhsm+r2yXC3w5PnF9CbYKVYVu2MNxoH7ahz42WYJaykOUW9hL/sl4/CGd9f7d+uF/TCeh/FC+b3JyFlmBde8Rj2ClCnIFBfAhT5lBuLTicC/ViF4HSCLczsJBO2TvtmXIWLFjxqIXJOIWJGJ4bkKrlSuXrcDEsWKG0MR2ljnsg/14M6DPXytlJg4Cz0BGsa4ozN0IgTRJM+cIu6dBYFIxmYHPVGYZvJ5gdDM7LShFNEszL/MuVsx5tPHc2nnqr1CniFpeE9TwP20eVhkJAn0AY//Cw/hkG9N9gfF8YF0ycz6RJgpHJjdpmrpaY43TrmVIQvKGvGBwPO7z79Vq9/z0a+d6tv1/f436iCgjGG/r28DIzNtwHRlYpPxj5iny7zk9Q4Ia+IzRDsFVmZ4ONScZyjkZtG5ObVGOTylGK29B3jGYIFqJlr/JM2Rq1Zb6wIjzTM+UsCtvQd5RmCEbZjPeLf1Unc3xciI90MsfiS3guSS7+9j9Q2Ia+ozRDMKk89T7Q6NXKbKS8QsMNBaemCiFZLBYJeuotbvG9+oQFqyYYAkaZRVsCEzGqCfAym9Eie0OiSM/fHm2SbZrtAy1Jfpr+2/QbdPxSb8oADF3NC79dgbtXuoYhGmj7EOKivydvYPRKYzTkSFP03thxrSyUAj1/h2OjibYPIWYKv1//+/peqJOKMLux/hVW4Kha0mIpbveKadmGaKDtA6123S7E4rtvF/8BEcd/wnvx3WLxXmTV2sVS2O6Vrj8ApZpn+5BIrPCFblD5EWWeemGh6ahdCosSx3h8TrNsHxKsBLhbr9cOdoZ/XP/X+hv48r+xi8zpqZeDws5x7pnmCDrJJtk+UO11W6YymkAQ70laU1bpxMoJdvAcZ+wPm2T7kCppVnQmPZ79Raeek2CqSqAMijzHmWZnzbF9cPAMZpgmL+UKi6OyvXrKoi26dJPKUfgcZ21UaSgeFzDzVGeodjtQZzgrpAjxAHyOszYtCughRpzhDpWNW4pXWmeo4HOctfGoSYFKK5mqzlBtuSyrM1RoCtYdG8E2W/gy1ygYrXDNQadEeTnQZktV4xYVU5F9FF1jlbHdGza3zI2K5cktZUmT5awzRKmS8ML9RCfQjhKF6De3poOWTWS8xGLRJYaLaMumCqfKRt+6yB01toVFkmTamG9cvnE5B1ql2kbDK3/pfGzl5RCqcKOQjSlnoSOYa9uNtS6Ks+L5EEt7IaCPVXqqpEzvAbhEIBdRZr6Rqp1Eq8x24/X0Ks6k+fhLr0Qwqs5oZemnlZzgZBmaWvFLlC9TmEnzUa5EsEmadYYYZMRt2d7ZuPyaFGfSfIwrEYzOM8dMYTuLDzEFlUh9Hz1NCjRpPsKVCOZRajfexocqc5gWX3h4AjZpPpNQTZZDig+9bdl88YWHJ2CT5jOJ1T6HVaziQ091htFrD2Fs0nwa5eJAW4eW0oPgcA7TrmUJVdhnwSbNJ8mqrlfQGULb8uQS48OkhCrss2CT5pPMMR5syXacYtk8NDKaPMdlVdmcgE2aTzKT5N0QQYAxk9C01D7YqNxVlKOwSfMLODTJyo6id7Kj6OXrT5Z3YZPmFyBzwwlMuTws7V3KOY5m4etPlndhk+YXoC0oMOuSM9xQmcjsnL03ig8VbNL8AhLDDRi9ZhIrAURKUrXfKj5UsEnzcSZku0GjFy1RqmRHzpPNi0ZLMNO89iixFTq0lgKDGI1eOwXzTsmFhyfQF8y65onzTG2rhD84ejlPzqxsv149wAF0gw4b9LriukSp4kGYdFGd4YzmykmqcWZl0WguYFrGaHDma2smWFtVg6JCVIWN+xzESrWtEjZU5kZLsNH5ctVOsBQDC0+mcjKRc1AqxtGr9aZz5SdoCXaozM3tmeZA1N72Qa0hY1FvG1VqqdFLeflWAi3BDpW52QGEIYPa2z6oLSlpnKg9RGoHc5UorMwN5QnM2ts+4AG+DnaGMznDthW/elnUCYqsSxwY/X3bB9s0P6uLYA52e3jQ8hJrRD1Kc2Q+iBWiQMHcofXM9kHUqIWp8xtmkVyt1OGHGMS33jRxeIDiBHOHQ1fU2PZhgnvKQS2RprE6mo2S8m+ZNzxEcWMY6iVqbPvQxlw8zrrmqv56VbW2pShMsKx4u762D3OZYg3ABANFOofoLVe9jlOYYIMs0K+h7YNDW/7TeCWX6DPvyEsPZSsT3r2CSd3UgSFstczKNWayMtPk57BgOFrByIWjV1zQoWxlwoJRUW8ost1DlQw0dmm2YEm8xPThfCYTTEFNZJU7Q0WzBaOa0DbmeWc0TQ7fdDH5LBorGI1TtMMBQnknLejAw/JpqmCYLqSSwxDXKTHuqFpK4whNFSzEUN6TLVxSXm02m9eBhgqmikLT1FGVUTCa1aRHbKpgbenJtIVtK65BZLhL8wRbxhNMZiRotrEs3y6vaJonWCxnAo9IgcgQt5k7qwonDg/QOMHQmtdrY4xIOytrR+MEW8l2CtF8gusntRq8MpohmErneqEjHIjf269vOFQcjRDMSdOEQvmZ01ZO2HUJ4p9z7YK1MKRYUUENhoUr1RlWPid/nCsXLJJ0BEAMOiUyTWZ17gwVVy4YKITpp2VLhlhhmKRVXv0/i+sWDHTC3FNICY0UgsKkvoNXxtUKFsnZBGTCIqg0eRNnw3K4QsEoRYgrXDKkqfGEYvpVWsdZ13OuSrBW2MJ9kqCUM5MeTJDrssiVgysRLGllm8hbuJ1hBgE8jFvJvO4RxgFqLBgdgOLM0gka8srUa8NfSaWF0L7i6+gAn1NfwZZp5p0hW6BXCP1fSmXWOFY50fX1hRl1FMxpQweI/Z/0VnLWlpjBEI4XYnDRCmucxTiHegnWSsOliGBOJeeqYcmZg/1hHddJNKm8YNByWrjMCFPgZI6tClRaebFM8dH5HIeq5Mob1RNKEqwQF4HIi6IWNqfYw1kVMJtEcxmjPtc7SJ2gHMH0XQQiUCJph6G3bM2USKtohaOVM/FKOwi0TpQj2GkXAQwQnAjxAPjQwo+h3CWEx+bYnBKvqe3pOeUINj6+Kd0LAXmMeOXNw7CdQHfoVXnTz9tRjmBbFwHzmVs6iRWHc9QNG5W3jKIlfGhja+OGdJJyBKuji0BNKEewOroI1IRyBKuji0BNKGviXCcXgVpR+UwH8xQWrGawYDWDBasZ5Qv2p89Mpjh++z9lC/ajP+3d8aJ/7mU/7BXc+n//r2zB9jnviJ0yLm7orS+lmb+1Ggt2UQxyWQDTzFszDMMwBeF2zK7elYM+HRWyU96T5+KO2cHL1EEjeembZk/71gBdrXfrnkknpGjf+mKGw5GheYpjgBfulvfkwByaQSCyg0ZyE5j755PkomNY2re2hpbV1b/1xdiGK8bnnU77HBRst7wnBy6eWyGyg0Z0MC3dW9NvW/vWKhOle+vL6QzpZ9cDBdst78lHzxDZQSMa15p+R/vWlgk/se6tLd8fuRf81JeCYl0kmKX7T+9Ch5QdNKLzzzZ6ureG/6OmpX3rvt0PRvo/9cX0ggsF2y3vyUMXSxayg0Z06Pi6t/YDK/DNS25taf/UlzMwumKoe2cUbLe8JwddVWKiDhrJ/6+G/2WB7q0hzLOC3gW3hhFQ89ZFMPKDoV64YxkGRFu75T05LzZEdtBIbowggCFf89ZC9Sm6tx76e0ejvDZd3QW4rm1TZY+tM4/Di/HG2UEjOXHt/fNJ8jHoat96YKvfl+6tGYZhGIZ5ZQYQXeOov83vmHsfX+BNckJNxzYsmL/2tQTjBd83gLIM40DY0NYgTjcHPUwhWD1UcDCyOiBLZzwa0Gt7Y6tnw+NdE15sjeiCjpU9ybwOJFjXAA1cbGe+sETX7/U78NH1x33fhClqfzjGl/Z8s2/issJoJAKrb8I1YjTsj/TS/IweKo9n2CalI8cmCmbgWqpFSz69ANc/+pTlVEm7Ycc1bHyKukSjA5+XXnPL/AAJ1jewe7MNaGlYOtazjKHYJtJBD3tHsF7QCTKZUbCRWvllXgv81Q+GFgUQwXCILct1oZFRS+titm8r2MjCBU/X9ztigK3KhWtwsXjw1j9Do7Ax+zt2hdop36Gu0LcCbGFi5A9hXNsKNgiCIT6K3aFpWD52o31/ONRbHGH0cOmQadVMXJyQdTGx2qWPKi8Lf90sw0q53hGlx/E16hqt3C3zagz9wVv/E5g8qAbJMAzDMAzDnOT/AfNLoKUWeOkBAAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE5VDEwOjI5OjIzKzAwOjAwFVvQdwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xOVQxMDoyOToyMyswMDowMGQGaMsAAAAgdEVYdHBkZjpIaVJlc0JvdW5kaW5nQm94ADQzMngyODgrMCswENDwCAAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNCQxalcAAABKdEVYdHNpZ25hdHVyZQA0ZmQ5ZTkxZjgxM2E4YTUwNDU0MGRiMWVmYWRhNmMyNzcxODAyNmRhMjRhY2E2YjRhNjE1ZWNkZGRlOTI0NzcxfRA+jwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<wand.image.Image: 4fd9e91 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "out/cancer_1/simOutput/wholeTumourVAFHistogram.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAKJQTFRF////ODg4AAAAiYmJd3d3REREu7u74e31w9vr0uTwiLbXH3e0W5vI3d3d7u7uqqqqZmZmVVVVERERIiIiMzMzmZmZzMzMICAgiIiIQ43AO4i9S5LDpcnh8Pb6UJTES5HDXp3JbKbOl8DcutXoc6rQuNTnV1dXP1VkDjhUHD9YKkhdW2ZuoqKiQEBAGBsdV2NsVmNsFxsdXWduYmpwR0dHwMDAw84SRwAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHRdC9SyJAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAAC29JREFUeNrtnQmDm7gZhllV0hDHXBI4NN0j3bab3aRHtsf//2uVBAZmbBPJfGCUeZ91vJhLYx4jhNAnJQkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAKDgOwYo+cPawjh/9E/mm0LItVOAMFIgLDIgLDIgLDIgLDIgLDIgLDIgLDIeLuwpTd88+iDExMOFHd6+PT76IMQEhEUGhEUGhEUGhEUGhEUGhEUGhEUGhEUGhEUGhEUGhEUGubCMc+ne834GhJFCLaxgJa+STFVC9cYgjBRqYdrp4XWSqKabA2GkEAvLmGD1KamMJGFFNUK8gzBKiIVJVmWldrJELwpnGCnEwnImjbSkLCFsJVa4hjUsaXSes6KbA2GkkJcSlWBtktdalf0cCCOF/D4sl5n9nzydZ0AYKajpiAwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwwIiwxqYQ3nPEN82HoQCMvapBFV70fUQpwQH7YeBMLaMtG8bqb7Q3zYahAIM3kgy89ahFJljviw9SAQVqhaJGXbf5CFLhEfth4UhY7GFDN4NnxsBeLD1oO6lGi8VQLxYetBc4Z1ZXkLq5U6IT5sPQiESVaOwjLZ7RDxYStBUkpcshjCwiAQdtLZ3GIII4UkS7Tc3A+EkUIgLJeW/NZiCCOFopSYccFv7wbCSKGo/FWCV+zmdQzCSKGo/K3MW3lTC4SRQiFMmLeqvbX468L+eDweHn0cooGi0KG0qNWCLPE9TjJ/aAodnN++FYMwUhYL42zxfRiEBbBYWCYX34dBWAAUVVO2LUBzurUYwkihqJqyRipkiduwXBjrQClxG5YLkw2fu4RBGC0UxfochY7tILlxVloxFDq2gaRqSvK8RqFjG0iaCEgxowXCSKFoSKpzJs4Nsy+BMFIoCh0yOfGbviCMFvqaDhkcbgRhAZDXdOSaycBwIwgLgLymg1dGWFi4EYQFQF3TIZUdL6cKil6BsACoRzfSRWKEDeFG3Jx7EEYJRaspPgZDtEzYEd/Cwo0gLAAiYec2HdJMszILCzeCsAAo6hLtWzneiJksMSzcCMICoLqGVaMWVwAJCTeCsAAobpyFuW6hmdtGUAVD4AHmRlAIE2jmth0Uz8NqPHHeDgJhTTm7GMJIocgSdTnXVhvCSCEJmTXFRIE2HdtA0USgnF8MYZRQnGFLexGAsAD20IsAhAWwh14EICyAPfT5C2EBQFhkQFhkLBbWNNn8LiCMlOUxzlWxtDc3CAtgeYyzZijWbwhVQ9LbQBgpe+gcDMIC2EPnYBAWwB46B4OwAPbQORiEBYDOwSKDunOwgrvYPsSHrQV11ZTmXPHw+DDzevvoQxEHK9QlchEeHwZhvpALa7hqw+PDIMwXivuwZ+1vuGBNeHwYhPmyQiOcVgUORwVhAVAP5WEmuQ4cjgrCAqBuhMO0NqWN4PgwCPOFuhFO3k8FxodBmC8kPeGYQuHCptoQ5gtFKZGZO6+FdYkQ5gvNgG+Le3ODMF8owo2EESZwhm0DSW19rbRa1vIXwnwh6fO3HevmL4EwUiiEFTWrb+aIEEYLyY0zl61a1swNwnyhGhYYpcSNIAroyxb2qg1hvhAMR0XR8hfCfCEYjkoSBPRBmC87CTeCMF920nURhPmyk66LIMyXnXRdBGG+7KTrIgjzZSddF0GYLzvpugjCfNlJ10UQ5stOui6CMF920nURhPlCPDLEFSCMFOKRIWx8mEzC48MgzBfqkSFUWbIG8WHrQT0yhLFUIj5sRVYYGUJzxIetB/3IEKXJC4f4sEaIdxBGCfnzsFKdEsSHrcfyJ86cT4v1zhfiw9aDSJjuazpcrYdAfNh6kGSJjT4XCd31zJ5jiA9bCZKh7dVMS20Io2V5lihmdUEYMYuFSbRL3JTFwnK0S9wUtEuMDAiLDAiLDAiLjL0JS81EmiTdO7gEwiIDwiJjR8IGIGwGCIsMCIuM/Qr70/H4/aOPzg7ZrzDD4dFHZ4dAWGRAWGRAWGRAWGRAWGRAWGRAWGRQC8sK1wb4nvgwCPOBWphwbYDvig+DMB/os0Qr7K74MAjzYR1hd8WHQZgP6wi7Kz7sEgi7ZB1hd8WHQZgP6wi7Kz4MwnwgLyXeHx8GYT5QCzvdHx8GYT7ssaYDwmaAsMiAsMiAsMiAsMiAsMiAsMiISdhTmr45T/+QpukPSdK9vyZiEmZWPZ6nX22UC4RFBoRFxq6FvX/u7HA18O+VOYOwyICwyIhG2PfH449nYW+Ox58gbC2IhE2D19PBE4TRA2Gk7F3Yj8fjBxc7u5owk9cel2y/MXsX9v58nq0m7PA2ql6TIAzCXkAh7MPh8OFS2J8Ph0th6eEQ1vnAIOzpcDg8rX00FhOHsOlrFDaQPt9f2CVpEJZGUXyBMAizLIoP+4qwNP35xXo/p6nJy96k6V/S9K9O2PTRWc/w6Kxbr8v8ntwGZpd/6/b6WoUtiw/7irCrmGLJ8bzG8ep5Nug4Tm7wDpeBMztnHWHL4sMgbIZ1hE3jw375js/x8ddff3OvK/x27XWdj5x/Oq/xye3104uEzEqf3US33sdJ6iOf+d559/dVfgZnYZLzf/zyPMV/fv78r2uvL+fXOHW50u9X1/viucvf/ZP+Mpv09b/y3/NJf5lN+uYX/4+RNDl8//3fGsKm8WEXMhf9Epb9jl5l0j5M48N29Ke/zqR9mMaHvYSH7Ih041eatB9jfBgAAIArNLwvh0wqrrwZtik4D793HJI2F4Tm7qTtVHDaQ9LmFif4OiHb5uVutoMr85/72pOKK1/GbVRZstBjPiSdJC0LLXONSRes5NW9SReMVywL21gyLS6+wWaYcn6rXeKTiivv7z1sYw5dGXrMh6Td0Q/ceExa33HIxm9t0lVt6OZcvPwGm3FiufnF2KlpxZUvz7YJPXBj0mYfPFTYkHTGBKsDc7UxaclargLPsEHY5Btshk2vS1PcIWy6TRman45Jt3USLGysZ2NVVgb+zMekT7qu63uFjbvZjsxk4F2acxVXt5hs0w3Gfl/SSgsdei0Yks6ZDD5qY9KV2U0dnKf2wsbdbIhu7eUgl7MVV7fot5H5Hb7OSWenhJscUQdePsekTVbchB61IWl74Q2/CDphNmm3m21pWG2OtfmdzFVc3aLfhkk3DHFortYl3f1ag7PEIemkUIKFFhuGpKXSWgeXEt2Ay+bV7WZjMjt2cG5vZO6puOp6RMrdMMTBm7ukM7dZFp52n7T944MvQmPS5i8P3TbvBlx2e5DBt64AAAAA2ITG1tjKoRx9yqZTV4rXXATXCgNKhK0irofbpvFO3U5dCjuxIrwECAgR2tyvMmbl5FmTacEzWQmeuynjMiuFuafibVW6oSxKbea2Jc9bYWdwwZvErGbL4pn5kCdmYWW7M65EI83mp63rzr95BFcZL40wWw0n8rqUedMUdemmRJKrqlA8YXVR2ycmeVubuaotyrooddKqxmwqhWsHqytZl4kws1TSKF5wyXJXYQUoEbzkujgLcxlh1vKqmxKJPeiNto8tClf1YeXYJayVkkk31QszH2XL7MKc9Tlr3dp6RUCKMEdaW1vGTNFpqitejsJc/betz30urOScZ7YAohNZd8JsA85uYS+s0Vs/nHoFmENbtU5YkfSajJ3hDDsxV8N3IczW5WdJVZkV3Fmo+ckN/9MLs1W4uclPwx9Fgq/QnQu20KHqWtja05NQte6n7FOaWhWXwgqzNktOZqEyp6TWNU84E7o3nmRmjt2YoWaPmu7Gq6tqzvva22HKvrvqVPPPfXDvbpNc9tW8wr1nmV1Tdgvd3uxiFDl2yMxjmDr8WT9YnZnHMHc8ZQEAAAAA+Gb5P1oVTpEBTrD3AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE5VDEwOjI5OjIzKzAwOjAwFVvQdwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xOVQxMDoyOToyMyswMDowMGQGaMsAAAAgdEVYdHBkZjpIaVJlc0JvdW5kaW5nQm94ADQzMngyODgrMCswENDwCAAAABN0RVh0cGRmOlZlcnNpb24AUERGLTEuNCQxalcAAABKdEVYdHNpZ25hdHVyZQA2NjliYTExODMyYjNiMTgwNzVmMTFlOGQ2MTQ5MzNlZTNmNWJjMDgwZmRhNjA3MzkzZjExYTkwY2YxOGZhYjQ0vswFQAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<wand.image.Image: 669ba11 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "for pdf in pdfs:\n",
    "    print(pdf)\n",
    "    display(WImage(filename=pdf))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Load the mutation matrix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "mtx=unpickle(cancer_sim,'mtx.p').toarray()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Plot the tumour as a 2D map color coding the mutation ID"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x7f3798e37390>"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAREAAADzCAYAAAChbyKzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3dd5xU5fX48c/ZpcMiEDoiIMVIXaUpoCBYQFBKBAVRVCLEICaWIJYoBlH8KRIbIChCYgHMV5AqVQQLCii9gyB1lyJdYJc9vz/uHRhhy7TdO7N73nnd18zcmbnzzJg9PP2IqmKMMaGK87oAxpjYZkHEGBMWCyLGmLBYEDHGhMWCiDEmLBZEjDFhyed1AYzJq6R0IeVMWmAvPpYyW1XbZm+JQmNBxBivpKQR16x8QC9Nm72zdDaXJmQWRIzxjCDx4nUhwmZBxBivCEicBRFjTIgECyLGmHAIxMXH/gCp599ARNqKyEYR2SIiA70uT0ZEZLuIrBaRFSKyzD1XSkTmishm97akx2UcKyLJIrLG71yGZRSRp9zffaOI3BIl5R0kIrvd33mFiNwaLeWNPEHiAjuimadBRETigXeAdkBtoLuI1PayTFm4QVUTVbWR+3ggMF9VawLz3cdeGgdcOAyYbhnd3/kuoI77nhHuf4+cNI6Lywsw3P2dE1V1JkRNeSPL7ROxIBKeJsAWVd2mqmeACUBHj8sUjI7AePf+eKCTh2VBVRcBhy44nVEZOwITVPW0qv4MbMH575FjMihvRjwvb6T5+kQsiISnErDT7/Eu91w0UmCOiCwXkT7uuXKquhfAvS3rWekyllEZo/m3f1hEVrnNHV/zK5rLGxoBiY8L6IhmXpcuvRAbrbskNVfVq3GaXv1E5HqvCxSmaP3tRwLVgURgLzDMPR+t5Q2DEBcX2BHNvA4iu4DKfo8vBfZ4VJZMqeoe9zYZmIxTlU4SkQoA7m2ydyXMUEZljMrfXlWTVPWsqqYBYzjfZInK8oZFQOIloCOaeR1ElgI1RaSaiBTA6Tib6nGZLiIiRUUkwXcfuBlYg1PWXu7LegGfe1PCTGVUxqnAXSJSUESqATWBHzwo3+/4Ap6rM87vDFFa3nA4fSJxAR0BX1MkXkR+EpHp7uOgR+dEpKE7ErlFRN4UkUyjmKfzRFQ1VUQeBmYD8cBYVV3rZZkyUA6Y7P6W+YCPVfULEVkKTBKR3sAvQFcPy4iIfAK0AkqLyC7geWAo6ZRRVdeKyCRgHZAK9FPVs1FQ3lYikojTVNkO9I2W8kZc9sxY/RuwHijuPvaNzg11p1AMBJ68YLSrIjBPRGq5v+lIoA+wBJiJMxo2K8OvYRs1G+ONfGWKaPEuNQJ67a+jVy/3m1qQLhG5FGcEbgjwmKp2EJGNQCtV3evW8haq6hUi8hSAqr7svnc2MAgncH+pqn90z3d33983w+8R0DcwxkSeEEynaWnfJEfXaFUdfcFr/g0MABL8zv1udE5E/Efnlvi9zjfaleLev/B8hiyIGOMRX59IgA5kVhMRkQ5AsqouF5FWAX78hTST8xmyIGKMV9zRmQhpDtzuLhMoBBQXkQ9xR+f8mjNZjc7tcu9feD5DXo/OGJOHRW7tjKo+paqXqmpVnA7TBarakyBH59ymzzERucYdlbmXLEYdsy2IBLuwzm8WaEyItfKClTka5cC096HATSKyGbjJfYw7Cuob7fqC3492PQS8h7O0YCuZjMxANo3OuAujNrmF3oUzH6S7qq7L5D3Lsup9jiaxVl6wMkebAuWLaumedQJ67d5hS7McnfFKdvWJnFtYByAivoV1GQYRY/KiaJ+NGojsCiLpLZZq6v8Ct5raByCO+IYFKUyClIiZSSuxVl6wMueE4xw5oKplAnqxRP8K3UBkVxDJcpjIHeMeDZAgJbQRLbKpKMbknIXM2BHoawWIj/IVuoHIriCS+xZLGRNpAvFBrIuJVtkVRM4trAN24ww59cimzzImJgliQSQjMbSwzhjv5JKNmrNtxqq7N+bM7Lq+MbHOaiLGmPAIxNsQrzEmVIJ1rBpjwiAiNsRrjAmP1USMMSETgTgLIsaY0FlzxhgTBrEZq8aY8FhNxBgTBqcmYvNEjDEhyi3zRGL/GxgTq0SIi48L6Mj6UlJIRH4QkZUislZEXnDPDxKR3SKywj1u9XtP7GfAMyYvi3BN5DTQWlWPi0h+4GsR8e2NOlxVX/vdZ0cwA57VRIzxijtjNZAjK+o47j7M7x6Z7QjXEZigqqdV9WecTZmbuGkliqvqd+pswPwfoFNmn21BxBiP+GoigRy4GfD8jot2wXeTea/AyS0zV1W/d596WERWichYv4Te6W1hWsk9LAOeMbFAJKjtETPNgAfgNkUSRaQETgL6ujhNk8E4tZLBwDDgASKYAc9qIsZ4RoiLiwvoCIaqHgYWAm1VNUlVz6pqGjAGJxMDWAY8Y2Kfb8ZqgM2ZLK4lZdwaCCJSGLgR2OD2cfh0Bta49yOWAc+aM8Z4RsgXuRmrFYDxbuK4OGCSqk4Xkf+KSCJOk2Q70BecDHgi4suAl8rFGfDGAYVxRmUyzYBnQcQYjzg1kfiIXEtVVwFXpXP+nkzeMwQYks75ZUDdQD/bgogxHnH2WI1MEPGSBRFjPJTPgogxJlQittu7MSYMgtVEjDHhEOsTMcaEwTpWjTFhseaMMSY81pwxxoTDqYnY6IwxJkRiNRFjTLgsiBhjQiaIdawaY0JnzRljTFhsiNcYEx5bO2OMCYezUbPVRIwxIcotHath1aVEZLubKWuFiCxzz5USkbkistm9LZnVdYzJk9ydzQI5srxUxhnwMvx7jFQGvEg0yG5Q1US/7ewHAvNVtSYw331sjLmAryYSyBEAXwa8BkAi0FZEriGDv8cLMuC1BUa4+7PC+Qx4Nd2jbWYfnB29Oh2B8e798WSRPcuYvCqSu71nkgEvo7/HqMmAp8AcEVnul5GrnLvtPO5t2fTeKCJ9fNm8UjgTZjGMiUVCvMQHdBB6BryM/h6jJgNec1XdIyJlgbkisiHQN6rqaGA0QIKUyDTDljG5kSDESbZnwMv449O5RCbnMxRWEFHVPe5tsohMxsmulSQiFVR1r1s1Sg7nM4zJzeIl8qMzqnpYRBbi9GVk9PfofQY8ESkqIgm++8DNONm1pgK93Jf1IovsWcbkVb6aSCBHltfKIAMeGf89RkUGvHI4VSbfdT5W1S9EZCkwSUR6A78AXcP4DGNyLyGY5kxWMsqA9x3p/D1GRQY8Vd0GNEjn/EGgTajXNSavECIXRDLJgJfh36NlwDMm5km29InkNAsixngkkjURL1kQMcYzQQ3xRi0LIsZ4RCR7hnhzmgURYzxjNRFjTJgsiBhjQhbktPeoZUHEGA9Zn4gxJmQiVhMxxoTJgogxJmSCEI81Z4wxYbCaiDEmZDbt3RgTJutYNcaEwaa9G2PCZDURY0yYckMQif1vYEyMEmfie0BHltcSqSwiX4rIejcD3t/c84NEZLebpXKFiNzq956IZMCzmogxHorL/O8zGKnA46r6o7uB+nIRmes+N1xVX/N/8QUZ8CoC80SklrvPqi8D3hJgJs6u8Rnus2o1EWM8FKmaiKruVdUf3fvHgPVknnQqajLgGWNC5DRmAvsfAWTAO3ddkao4mzZ/7556WERWichYv4TeEcuAZ0Ekl/u20zq+7bTO62KY9IizCC+QAzcDnt8xOt1LihQD/g/4u6oexWmaVMdJ8r0XGHb+0y+S8xnwTPRY1HYVcfnjaTGtzrlzi29dzXVT6nlYKpOVQJoqAV9LJD9OAPlIVT8DUNUkv+fHANPdh95nwDPR4/uuG5k66g1aTKvD17etAWDvv47z6uBH2PiPZDb+I5mfHtjOj722AbC0xxaW9tgCwNcd1nhWbkOAPSJZd766IyjvA+tV9XW/8xX8XtYZJ0slREkGPBMFOv/Yhs60oc/IoYzePpAVg58DoG/723jlv5/Q4porOXjXr3zXZT3XfnYlAAWKFaDB6Mv46paV51aRLmq3muKVipP4XhXPvkteIwgSuXkizYF7gNUissI99zTQXUQScZok24G+ENkMeFYTyQUefXo4DeteTs/B/yTx/SqkvCUkHT3Ik/d05/MPF1N7VnWKlUsAoMb0Kvzp9uuc+9dWI18h59+R62fVy6Lla7JDBEdnvlZVUdX6qproHjNV9R5Vreeev92tafjeM0RVq6vqFao6y+/8MlWt6z73sDtKkyELIjHu2eFO/9q7nV7lulb16PLTjQzvNhCApKMHqdG0Gu+Nmka9kZfy0Ja72NJhB6Dc+sN1VBpUHBSOv+78A5T4vtVCclqkmjNesiASY4qOLwbAsru3nDu3+P7prDm8kv6tenJz1TYMKPcyG/bsOPd8s89rA9Br4As8tOUuBrToR6daHfjx3m2sev9zTpw8xfHXz1Lnixp803Etl3x0Cbf+cF3OfrE8KXIzVr1kfSIxZE2/3aQuToU/Q6P3arBv8Am+qDWR1pO60O5sUw73OMKoX18G4Ms5K5h56FuaT6vDvsEnKPuH4tQ+cTm9Bg7iyqbV2bxiBxLv/J+z3LNFSXrxBBt+2Unzz+tw5PMjzGQxl3x0CUfuPuLlV87VBCLZJ+KZ2P8Gechvh06ScuIMhUsW5vtuG/lj9YrsbLiSa2+qx52Nb6Hv5m58tWsxp/6ttGhdlz3/XsrBl08RHydM7fIfypQqTt8+Hfnmr1NIHv0TKSfPwOqigBNIzvY5zb3rbueSjy7hrtXt+G7pBo+/ce6XG5ozVhOJITUSL2NBv0/49sAihlUdy/ST3/BSuwE826E0XZ59gu5dW/PxhHmoQrMbzs8XKfNMEa6P/xOlBhaCOTB//yyGVR3LrN0j6FPp4XOvG/Xry7z+7bu0q9eMCfVmUZsK6RXDRIxEfVMlEBZEYsA3t69178CfnryPx7c/wNcd1rDzle9Ze+JHALp3bc3C5avpcVcbft6XxBWXVuZfW4awoNsEXrhqCIoyptNUfluXxqiRU7hhbEMqbavFd53XA3Dm5BnqP9CJmW+9zR3P/4MiA4py9kwqV/7bAkl2ymKBbEywIBIDCiYU5PTR01S+qhI9XmpPi4JtiC8QT/c6vfmmozOlvdrVlan4fAIb2c4HR4ZT656baT61DlXubEaTibXY8HgSzabU5jtZT4t2DTj1wEkeGTuAayc7c0d6retI32eGcsfz/2D8s4N4aOSr1K9bjRROefnVczUB37qYmBb73yAPaN/pGv7vrdeoWK4UAOX7XU3Jqs795p/XZvaIkVR8PuHc630B5OYlzdkwZjY/3LmJBzvdyonXz3Jp/YpcVr4sC29awX3X3cbqh3bxl813AtBs8pV8/ep/GfvNZzSuX52UP1sAyV42OmNySFLXA3z06TTy9UvjY2ZwLU7tYU2/3RzaepB/VXoLgFZfN2buT8u44tjl/FBkExs+TeJNPqPJxFp80fQHij4Wz/qeB6jdoTr39+/AqJoTeWfz4/R/43XqvF2J25bdQPF8JTjS8yjrH91H8i1JtJzdwMuvnuvlhuZMliHOXT6cLCJr/M6VEpG5IrLZvS3p91y6uyWZ4DWcV5uG82qTPOQku5MO8e7hoUw/OfHcqtxDWw9y70PtAFj78G5urH49R46dBKBYuQSO7j7CH4eVA2D98p8p+UkJbm7fkMv/UJn/+2wxl0+9lFE1J3I86TgA0xp9CcCpN5TktecDSP05tXL0e+cVQW4FELUCKd04nJ2N/A0E5qtqTWC++/jC3ZLaAiNEcsF21h5ZfuM6lt+4jpSUVN7v+RJ9SwykVI96NJtSm4qflUXihO0ddwMw+ZnXebbCMEbd/yylny5M7Tcq0GRCLSZu+oAld2zgxP4THD32G4e6H6b3iy8CsO3289tGNFlQj4kbPuB46lF27NpP0/b1AbhlSXNWbXUW7q3ptzuHf4HcTwL8XzTLMoio6iLg0AWnOwLj3fvjOb/zUbq7JUWorHnOonarWNRuFZUGFadviYEsv2crzabU5tXkpzl6/CTFyhXjp/t/ZsCOB3m2wjCW37OVKeu+AGDAjgdJeVv4y7JHKVyyCDuGfs/ZPqcp/7/SlKxUgtRTqbRfej3JQ05yYOxKnm/8T/aX3UUaaVzxaln27T8MwNPPj+DQkROse2QPBzbu9/LnyH3EmWwWyBHNQi1dOd9CHve2rHs+o92STAiun1Wf2tfV5F97HqXh/DrE5XMqdQn5Ezh+73Fe6N+bB3t2oFmhG3hm11/Z+dr3fL9hAxseS+KBN17kp8Ub6d/yHp78cw+aTeoCwL47DjCzzX+YO/xdZjReRNlnitA+oRsbF/1Cm8EPAvDR0ZHs+mk3jRfUpeXsBpQaWJDylf7g2e+QW+Wl5kwwAt4VSUT6+LZ6S+FMhIuRe5R+ujDPVRzOex/O4KoPqvJjr200+fuddPqxNfdV+jPvjJ5MhUca0Th/CzqV68Eve/ZzPOkY+QrmY/XoaXR4rD8ARYsUYu+/jrOyzy881OoRqhS9HIAZxyYx49gkAI7sPEKjvt3In1SYFjPqsrT1Gpb22AxAqScL0mpeojc/Qi6lqgEf0SzUIJLk2+zEvU12z2e0W9JFVHW0b6u3/BQIsRh5w8nhaVz1QVVODk/j9PEzrH5zOlOuXkDDEe0Z9Fhvmk6sRYmeziK7v3e5A4ANT84DoF/fLkz94RuSdx0iedchTh87zQ/dN7P+yBo+PDqCSo82pdKjTXl0x0D+dPf13H3vjSyKm8vmgQcAKFC0ICsf/IVvOq715svndmkBHlEs1CAyFejl3u/F+Z2P0t0tKbwimiKPxrGo7SpWfr2JwiUK0+CxjqS+HUefWzvw3ZZVlPjoElrNTaTj8tYUzFeARh/VIFl20+bVB1jaeg0F+iupp1JJPZVK4+uvpMknNSlZ0BlQO330NL8dOsn23cns7pzED2s288Kk96k5tDQADcZcRpWa5Wn+eZ3MimhClCdqIiLyCfAdcIWI7BKR3sBQ4CYR2Qzc5D5GVdcCvt2SvuD3uyWZMKWcPEPi+1VIfL8KmzbvZlbTxRy++wgNqtTi9Jsw6ZsF/OOVt5l+eCLvLH2P/Vud2sTqv+wkLl8cCRUSkIdSeS35GWZv/ZJ9p3aTlnqWuHzxVBviTmQrU4Jd25JY/+g+rp5Xm6U9NrN5xY7MimVCpaBpGtARzbKcbKaq3TN4qk0Grx8CDAmnUCYd4nS0Aqz7214uu6wMC29cwZfjx/LEhNdJeDwfp9+Exh/XoFbl1tR6pQwjNgxg8BXjYR/MfnkkAANfeYUnyg5h61MHWVJ7NTd3bsLazTuZcWwSnSbeS8qfT3H/PGfuyY83rqPgXwpRf1TlDItlwhShpoqIVMbJEVPevepoVX1DREoBE4GqONsjdlPVX933PAX0Bs4Cj6jqbPd8Q85vjzgT+Ftmu5tFd7evOccXQABqv1GBZbPWUKBoAV6o9CYrF27k+LBUCj4Ci25ZRa1XyvDSpy/zwR8nc+LEKfp2v537xzxHq7/ez/6XfqPn2tsonlCYZ4o8A0CT+jXp9tkDfNr9Paac+JCJ8xYxdfESdj1/lILFCgKwuL1t6JwdItic8WXAuxK4BujnztsKZU6XLwNeTfe4cJ7Y71gQiTH5xhQk35iCFChagISKxblvfSdaTKvDjA7/BaD3I7fRb2sPnjs8CBmVn9RTqXwwfTabvnUmjBUrUogP60yj3B9KUP7GQpz47TQAm7bvpVPRnvx19mPsW7mXU2dSOLTP2ZDo1BvRXZ2OWRFszmSSAS+oOV2hZMCztTMxJvVB54++8/+uZ+KoeYy7cgoA137Uhcu4hM3td/DVC06t4WDSEdJSzlLnrYoU+mcBWjaqw5xvV8CTUONsee675O8AzO+/kuH9+7Ov/UEWrltGhcSK1Hj5DyztfojTR09R6LWyXEddb75wLqZAEH2mpUVkmd/j0ZkksKrK+Qx4v5vTJSL+c7qW+L3NN6crhSAz4FkQiVETR83jzr/cSNK8A7Rfej2F8hVkVZNNLF+/hWJSiDtatOTV9yfw26HfWPKnDfyn1ws8OX4kS5/9HysP/cj/qzIGgLUP76HOWxV5osCIc+tsmn1ai2QOknIqlWbuVgEmOygE3ml6QFUbZfWiCzPgZbLAL2IZ8Kw5E6OKVyrOrGnO6PmMxovoVrQHXRvcyqaVv1D8ifwAbB4+n2LlihGXP54RNT7h0O7DFI1LoEGpq/nm9rUU/7A4a4c6mQIW/3P8uWsndz0IYAEku2lkh3jTy4BH8HO6LANeXnH1+Mu5okEV9g0+zs/P/sqS/F/R5aknWDR8HACL16+gYs+m/Pvvf6NJq9pMPTSBt598nCr9W9B48B00n1rn3D85nxx7l+UHlnv2XfIyTQvsyEpGGfAIck5XKBnwLIjEsKSDR/jtdAqqyvBq41g0fBxd336Cwh84my83m1KbBwYMZsPmXXy+/1M6lunKVWOrAvD0zodoUr0O3Yr15qukBbQs19rDb5KHpWlgR9Z8GfBai8gK97iV0OZ0PQS8h9PZupUsMuBJNMyGS5AS2ogWXhcjJm175hADut7FP98ZQ4MxVUgbEU/cX8+y5akDpJ1N49SR09R/tzI/3fczV15VlTn3TOSakbezZ2syk154ha7PDeDSmuUBWPflJlrOtU2IwrGQGcsD6bsAuKr+Vbpg5oKArluqcqmAr5vTrGM1xl0+pBQfl5tDgzFV+LrDGjqXbun8+9LNeT7pxROA0zN29Phv3F7qLspShLu/607X5wZw1QfVANj/0m/pd6mZbKME1lSJdtacyQUKPAJr+++hxfS6TP7vV+zv5nSMfnXzSlo1rEfyiyc5kXycss8UYfm9znyRmcuWnAsgZ9+Jo8zThWk5x2ohOStvr+I1UabOWxUBaDHt/HyOlnMasL7dNlo2PH+u4X+cLQBWL9/G9902AbB25facK6g5T4lkn4hnrDmTB6xvt43rqPe7c1ePq3bufv3Rl+V0kYwryisZAbEgYoyHon2FbiAsiBjjFV9zJsZZEDHGQ9acMcaEzBnijf0oYkHEGM9E//BtICyIGOMVtZqIMSZMFkSMMSFThTQLIsaYcFgQMcaEIfrTQQTCgogxHlGFtLMWRIwxYcgNzRlbxWuMVyKYMkJExopIsois8Ts3SER2X7DTme+5p0Rki4hsFJFb/M43FJHV7nNvSiY7PftYEDHGI4qSlpYW0BGAcaSfZGq4qia6x0yIbOIqsCBijHfcId5AjiwvpboIOBTgJ0cscRVYEDHGUzmQ0PthEVnlNndKuucqATv9XuNLUFWJIBNXgQURYzyjBFUTKS0iy/yOPgF8xEigOpAI7AWGuecjlrgKbHTGGO8EN2M1oAx4v7u8apLvvoiMAaa7DyOWuAqsJmKMhwJryoTanPFlvnN1BnwjNxFLXAVWEzHGM5FcOyMinwCtcJo9u4DngVYikojTJNkO9HU+V9eKiC9xVSoXJ64aBxTGSVqVaeIqsCBijKciNWNVVbunc/r9TF4/BBiSzvllQN2L35ExCyLGeMX2EzHGhMM32SzWWRAxxkN5Yu1MpObkG2MuEMG1M14KpCYyDngbZwqsv+Gq+pr/iQvm5FcE5olILb+eX2OMyzfZLNZlWROJxJz8MMpnTO4VwbUzXgpnslkwc/IvIiJ9fFN4UzgTRjGMiVXZO9ksp4QaRIKdk3/xSdXRqtpIVRvlp0CIxTAmdjk7m6UFdESzkEZnQpiTb4xJR24Y4g2pJhLsnPzwimhMLqVK2tnAjmiWZU0kgnPyjTF+nNGZ2K+JZBlEIjUn3xhzAZv2bowJj0Z9p2kgLIgY4xFnKwALIsaYMET7RLJAWBAxxiOq1pwxxoQpN9REbI9VY7wSwRmrGay2LyUic0Vks3tb0u85y4BnTKxTFE1LC+gIwDguzlY3EJivqjWB+e5jy4BnTG6SzRnwOgLj3fvjOZ/NLqIZ8KxPxBivuM2ZbFTOTQOBqu4VkbLu+UrAEr/X+VbbpxBCBjwLIsZ4JMhNiUqLyDK/x6NVdXSIH20Z8IzJFTSojZqDzoAHJIlIBbcWUgFIds9bBjxjcgMl2/cTmQr0cu/34nw2O8uAZ0yuEMEFeBmsth8KTBKR3sAvQFewDHjG5CKRm7GawWp7gDYZvN4y4BkT62wBnjEmbLlh2rsFEWM8YgvwjDFhUwsixpiQKVG/CXMgLIgY4xFrzhhjwmbNGWNM6KwmYowJl1qfiDEmVKqQlmo1EWNMqFStT8QYEx7rEzHGhEytJmKMCYtNNjPGhMuaM8aYkFlzxhgTnuzf7T1H2B6rxnhIU9MCOgIhItvd7HUrfDvDh5IFL1gWRIzxiLMAL7AjCDeoaqLfzvChZMELigURY7yizgK8QI4wBJUFL5QPsCBijIcinDJCgTkislxE+rjnfpcFD/DPgrfT770BZbtLj3WsGuORIPcTCSQDXnNV3eOmy5wrIhsyuV5I2e7SY0HEGA8F0VTJMgOequ5xb5NFZDJO8yTYLHhBs+aMMV5xZ6xGomNVRIqKSILvPnAzsIYgs+CF8jWioiZynCMHFjLjBHDA67IEoTSxVV6wMueEKgG/MrKTzcoBk53sl+QDPlbVL0RkKcFnwQtKVAQRVS0jIstCSFjsmVgrL1iZo40vF29ErqW6DWiQzvmDBJkFL1hREUSMyZNUOZsW0j/+UcWCiDEeUSyNZqRdOFwV7WKtvGBljjJWE4modMa8o1qslReszNHGEnobY8KkpIU2IBJVLIgY4xEFzlpNxBgTMlXSrE/EGBMqq4kYY8JjNRFjTDhsnogxJkw2T8QYEwabJ2KMCZNy1uaJGGNCZX0ixpjw2CpeY0w4nJqIBRFjTKhUc0VzxvZYNcYjvhmrgRyBEJG2bja7LSIyMHtLf57VRIzxTORmrLrZ694BbsLZyX2piExV1XUR+YBMWBAxxiMa2Y7VJsAWd69VRGQCTpY7CyLG5GYR7BNJL6Nd00hdPDMWRIzxyDGOzF6g00oH+PJCWWTAi1hGu2BZEDHGI6raNoKXi1hGu2DZ6IwxucNSoKaIVBORAsBdOFnusp3VRIzJBVQ1VUQeBmYD8cBYVV2bE58tqjnSbDLG5FLWnDHGhOY9EswAAAAsSURBVMWCiDEmLBZEjDFhsSBijAmLBRFjTFgsiBhjwmJBxBgTFgsixpiw/H96fOetszyHSAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 288x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "cmap = pyplot.cm.get_cmap('PRGn', mtx.max()+1)\n",
    "\n",
    "pyplot.matshow(mtx, cmap=cmap)\n",
    "pyplot.colorbar()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Dump the simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "cancer_sim.dump()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Reload the simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "restart_with_mods = casim.load_cancer_simulation(cancer_sim.dumpfile)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Change adv. mutant properties and run for 20 more generations."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "restart_with_mods.parameters.adv_mutant_division_probability = 1.0\n",
    "restart_with_mods.parameters.death_probability = 0.5\n",
    "restart_with_mods.parameters.number_of_generations = 20"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "outdir = 'restart_with_mods'\n",
    "clear_outdir(outdir)\n",
    "\n",
    "restart_with_mods.outdir = outdir"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2020-08-19 12:29:32,477 INFO: Ready to start CancerSim run with these parameters:\n",
      "2020-08-19 12:29:32,479 INFO: matrix_size = 200\n",
      "2020-08-19 12:29:32,481 INFO: number_of_generations = 20\n",
      "2020-08-19 12:29:32,482 INFO: division_probability = 0.6\n",
      "2020-08-19 12:29:32,483 INFO: adv_mutant_division_probability = 1.0\n",
      "2020-08-19 12:29:32,485 INFO: death_probability = 0.5\n",
      "2020-08-19 12:29:32,486 INFO: adv_mutant_death_probability = 0.05\n",
      "2020-08-19 12:29:32,487 INFO: mutation_probability = 0.5\n",
      "2020-08-19 12:29:32,489 INFO: adv_mutant_mutation_probability = 1.0\n",
      "2020-08-19 12:29:32,490 INFO: number_of_mutations_per_division = 10\n",
      "2020-08-19 12:29:32,491 INFO: adv_mutation_wait_time = 40\n",
      "2020-08-19 12:29:32,492 INFO: number_of_initial_mutations = 2\n",
      "2020-08-19 12:29:32,494 INFO: tumour_multiplicity = single\n",
      "2020-08-19 12:29:32,495 INFO: read_depth = 100\n",
      "2020-08-19 12:29:32,496 INFO: sampling_fraction = 0.9\n",
      "2020-08-19 12:29:32,497 INFO: plot_tumour_growth = True\n",
      "2020-08-19 12:29:32,498 INFO: export_tumour = True\n",
      "2020-08-19 12:29:32,499 INFO: Tumour growth in progress.\n",
      "2020-08-19 12:29:33,198 INFO: All generations finished. Starting tumour reconstruction.\n",
      "2020-08-19 12:29:33,209 INFO: Reconstruction done,  get statistics.\n",
      "2020-08-19 12:29:34,007 INFO: Exporting simulation data\n",
      "2020-08-19 12:29:34,184 INFO: Growth curve graph written to restart_with_mods/cancer_1/simOutput/growthCurve.pdf.\n",
      "2020-08-19 12:29:36,200 INFO: CancerSim run has finished.\n",
      "2020-08-19 12:29:36,202 INFO: Simulation output written to: restart_with_mods/cancer_1/simOutput.\n",
      "2020-08-19 12:29:36,202 INFO: Log files written to: restart_with_mods/cancer_1/log.\n",
      "2020-08-19 12:29:36,203 INFO: Consumed Wall time of this run: 3.699930 s.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "restart_with_mods.run()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "restart_with_mods/cancer_1/simOutput/growthCurve.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAQhQTFRF////7u7uu7u7d3d3mZmZqqqq3d3dERERzMzMAAAAVVVVREREODg4iYmJMzMziIiIIiIiZmZm8Pb6lsDciLbXH3e0w9vrpcnhW5vItNLmTJLDb6jPUJTEl8Dc0uTw4e31PYm+LoC5xt3sS5HDk77bc6nQc6rQS5LDuNTnutXoSZDCO4i91OXxjLrZx93sea3SaqTN4u71rc7koMbfSJDCdKrQhrbXbKXOV5nHqszjgbPVNYS7ICAgZqPMVZjGmsLddKvRSZHCeK7SlL/beq7SdKrRirnYZ6LMbKbOXp3JT5TETZPDXZ3JVpjGirjYfrHUYZ/KZ6PMQo3AU5bFV1dXoqKiR0dHwMDA92jkDQAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHSdkLBwlAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAAGS9JREFUeNrtnQtjqzh2xwXYxphZbAPCwDg84uRmp5290Jlud7adbl/bdvqYPrft9/8m1ZEEBpvYMsHXPPS/14ntKLHhZ4mjo3OOEJKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpKSkpJ6lBRVU+mdmUa/zMnXhaYpj35fUu9I03UghZaGAQ/IP7TQV6bOiX1lSHWpX3SAzKTArBUBZszReoM0CyF9xnlqj/5EjUqm2sUfASaaqRqklymIfFtpBUUJrGN1BmypLwgp8h9uZgFM1bSvJLAu1RmwlW5ahrkwFgBsu5U97E7qDNhM07aGhjZruIDNNopizNkPJbBO1QEwDWwXuAPj4cyw9CVSrI2+LX4sgXWpDoAtVCK4o6j0ETXn1WXxYwmsU3UyJF6UBNapJLBey3ZOn5HAei0XeyfPSGC9lgQ2MElgA5OLdyfPSGC9lovdk2cksF7Ll8CGJQlsQAoCCWxQCjFCX++fTp6VwHorjGUPG4xg9rXfS2BDUYSdoof5Jz+SwPomOyYdDDoWAEsksN7LTWP0fHhB6HCA2+vJjyWwvsnHXtHDbHKTPazvAmA74ITJHQms33J3DJhLOEUSWP8FeD69PaOXt0/o+Y3ceSN36uoAGEuGUGaatkANyRASmLgAWNHDPNbDkpMmXYS50WQIdaOZxqIhGUICE1cDMHzSpLNAUsCz0WQyxEdkQ38K8a4EFt0RGNHCmMtkiLayib1BLluIXLZeimvYL727AlMs8ywZAqKCJTAROaRjeYAHkzmYS+6QB+5dgSmWpSCZDNFWEAzwzJwbL+iV3PmG3Hk+/NHhpF131zDghWQyRFsBsLKH+eSOS3vYH9+hh7FkCPpVk8kQbfUFgbFkCPp1IZMh2uodYPvTwETp6eiB4hihp/0T+ma//xbt979CX+/35In90+f9mwTWQ/kp62GkU3keuA95D8scHEpgPZRPRsKccKoC88gDn0ygJbAeyufBACWwVALro5w0ZndCbNeBkfmyl0tgfRNNTckjtlR5BsyXwPomlwYD5BeB5afpKxLY45RRYH4NGLgPwwqws/QVCexx8qm/twSWgHHoOuDvZcDSRALrjVj0hnsERgxDCszlwBzS5ci/QALrh2CpsgYMV4DFOMUuA+ZJYP1QkekAwGzSn2rAyJPYDVMJrEcqgDnExIClynNg2O8SmKYsTXMp1nQywOJIqFmEgyOwYm35zsCWG7Q1t1uhttMBlmNO431FQSXTATcA80+A2aeRpO2AqRbazFVTqO10gIFL8KxHNDexGbC8BszDIfYzBgtuGU4w6gaYYpi6MlsJtZ0SMI9ckwgw2648GwUNTQijmMZhF8Ai3rF8QOXiDOess3UEDEGw70zsV6cDLKQhoARYkhNohFOcBqxTUaVuESWaF1GiJbDMqwODPgbA/K6A3aDpAMMFMDjJLo7Y+FgCawjrZcACcqcGDObLXQJTeQV1MRJTA5Yzi8ItHIXpRWDEzojcc2AemTpnnQFTVKaFUOtpAEt2BTC/sNlLzy5rQR+EFWBxdW2ZA0uTAhiz8O8yJPKtPJS1BhOzaWav2IWHqQlYjGwPgCXHJt5pMAB0qR2z5OvA0qT+Sh0MiXwrD8vaGupUs1cAzZ+wmN1XFsD7evgOfQcZyofDM/qe3Cphva/kicN35IkDxPg+vxz+9PX58Gvy789eD785/PLww+H5EB5efgsZzq8nob+dDIkQ+asaClptp5q94hUeJp91n8KNUfY7pxp0mBYrX7joYaGD/xz/Bjy+0Nn22Ev/gjbwjzZLca7bDYmLrVkx6wHY2oLdPE6zV2am+eMUgH3e7+H2hL7ef00f0KJDT/u/hPoon2nQIX2WNSmCDskd8uDz0/53+6f9X+3/ev9EHn79ef/D/vPXf0Mb/G1XwDarDVqXE2dgo5l8+5Va9gqaQA+DQspe4bCA7uOxBUnmd4oKY78hrJfcSTEOMlhNYe4Nn9sfsc8auB0BUy1koqNriu4MsQFgU8xecan3ogSGC5u9dBQyYE4jMDAumKVRA4Y6BrbQFRNptR62MJbImuRWHi5NmcS2+x6wAOXsQexVgTkNwPIKMK9DYMjc6JbBf5Vv5bHVN5YyxeyVIsja5fFO9AH17GbHflfGYecoIc+Sn+QuGQnjI7CgsOTvAgz2ximWw4qtPJYqezS17JXsHWB+0d3csxi2qAw6PAKDuwH8jbgAFrNeWg9MlL7ED8u/BCxuBOadAIOFFODkucXSWH4cVrsAphDbAmlyAZPKJ8McAbDLOI2Yn/XSjXEdWBEqdQTmdguMGohyAZPJZ2iohVeede8MWNYEDGyNPEnfBbbrCNiC9jC5gEmXKsMjsNBuBsYfOG4TML9Yqqy25sC8rq5hJmygKBcwUR6igobr8yHO43OotMwaQhB+zSNsUps1SdL6UiWdQ3u7CrC4S2CIV5YS0KiBJXQyTE52irNTYGWf8U9CosIiwmaH8wqw0pIvFjY7BSauUQPDHBicdVjLOgO2awBWiWFjiym+wxswmzItgIXUh1JPX5HAPiQ2cDFgAIEua8EFKji6MQSAeZilQjC/iF8AK4uTVtTSrL/ht0YOjBXwOgLjDkG3CoxOhcMmYAl20oQC849ercoSTVfA1I142zEDs9llB2JmkhQg7IAG3K0BY+NjgvOcfKkBYwNofndgSBf/tdEC29l0BKNz36JbVQz1EDsEjwv2IDSJWBfMyLM8SjRNqle8GjC3ABafvWpLYDJqin70i4prDJhbA+YXCya8Jx1DonZFSFQdGFy27BNg5y/bNpBURk15OKPAdmz5sRFYhvPwHFgZYWOfAUNHYD55iYaw77bZK6BpT5zhlD4fDi8vh787/Prw8nt69+/Jg384vP4j+QJ3Xw6vh9dXGlLzzO5CExph8wIRNofD4fkVAnPgwcvBZyE7vydPPJ/tMMDVHtiqWKG82vTRp/Y+clmnyF38hH+g8y3Sa34qehL0pSfW2Xxy2fqnYw+DJoci6JD3MJ4XwYPi3KKHNamlWQ9fZtP21nNgPoDZU4+Gw8D8M3N6uPhXBTCg9ATWI7h5T4DZ3H3oHkMay3DvJn1g4ryerrceLi4u86nT0gzMBeXVojR4aCj2c3IpKwM24IpHr2FBZfXMuT8w0zT1CQ+JLnXy0ZMdQtIdgRWcA+OWRsVg5CFRZA7teGcrX84RmINQ+E6e2QesxAk7f31qy0FoBpkjMyev1wiMzaH985Ao7xyYx7IlmnbbrqglMGVtbmdiTccJDDx+2GalGeh5rwFjLihuUVwA5hfAkiqwHb6QLN0SmLXRtnqd2EKjaRCTSIbw2cn2OLCUAbO9ipPXrwCjvquLwNjUK2dby6L4wku3TUpXIAy7+tRCX2nGbCLJENAhEsqDDIopPe82i/KlD4IzYGVI1CmwFAJMAZhzLT2aqyUwoLGuhQhAkPZqNZFkiCKLHJYiuReQL40dE4iozV4tRkTX/8sYNruSOuHRi2J0R2AQImDWrURF3671qWzlwXk4PK3VZ6lf58DQ6eIzcE6PETYcWESBIc8WeOnWRsepZ0oxN6auniZDqJr21ciAJW4JzOXA8sIgPwPmNgCj820bnS2kiKmzFWcyEE4jGQKHKKaLkuDz9bnblgBzdkdg1VmxR+NxToExN29eAMvvDUw7W16Zkcva1pxAMgSzCIt1SK/0UbCflQGjbhWY/w4wtwDm3hnYUp+fzpyVjW7p8/EnQ3glMDAm3DqwtFrbJq0CIzZlHLwPLAhFrl+gloGkm4YnVZoGMfZkCJhvPb/9y9vz28vbv7692LDd13fkRvWJ7v31Bnt/vZAHP8GD52fYyPITe/Yb9DN5AE3QT3Svy39DNjwQV8shcWVOdT0MrHAH/w7/O96B0507lbhBXqvPVgTU2+8HbLhFnTBxtVxe0S3TNMV8UyMCBrvo/Zb0mZe3T+Tf888h6Udvb9+hb8mN6qSHfULfvL19iyo97Jk2eX77D2iLSFvSNYlJcoNk1NQNCrADfSZwWXVXiJYnVypiLhQXoKKH0RqjZV/CODuNsCndvBcWUprVtodNMmqK51q6UE4+ocC8+gwKaATQY2qDXy2VtgTmhvAJ+ELAJho1VQDLyfWL+Z1OgGXVwuZ5FZh7xFhZ//dg5SsQK2TKJaOmbpBb5ql4+Ohnj6oNKlsH1Go98OVkm+anHwM2RGdfR8moqRvEU/f8KjAnrzRoBpYKr/8LSEZNCQp20eO7N/h8zdI9M8hZkUS/WjoxrVUxgiYBGTpvGgVrklFTgsrpmac5DSHdHSVkpVJqsqOyzF507EuwPl0BdvswWJWMmhIUxLmnNE+FJxBxi+JMdlGGtALMeziw6UVNJWy7jZgDC98FRhf4WWiGeyz0W0697I/xklFTouL7o3gcGJtkhe+1hgHUdvl8uswaCgLRV3tfbY0OcltOyjXFl/cdmnfCihHZ7xt7djWOpinNq7XaDonktp7SxNmmSykJZJG7bF4lbp27aYwimCJ3oVbATFZCdjI9zIl4JhiNmT8CcwSt8xhdDl27Ra2ALdXtlK5hMTMwfAguDMs1y456zK2SVQSuy+PAWIk8GpSNvKSrLnOjJLDrggLm3x/Cn18OUAj7+RmS7h4mCey6dtQN9cZrGrISUA9TG2BLDQkmQoCGD8zlwH6gtZMjhHF7V+CH1Wqjgc2iySul0vJTI0uGCHIbgEWwPxSL5L01CKNjfcSsr5PQdG2rjS4ZAlx/Gc1qcEtg/uCAsYlzXQtWjW9syRAZW94PILOcRcWzxf2HqbOEvplhGitldMkQPBPMDTEPyn4kLFBnCX2aoS2s7WkyhCYc+tFT+ThD//lfby9vn2g86Pfk9lh1ltA3N0aZDJHA2td/w46UNlsCa7+43406S+iDa9gIkyFgZZnV1uMxvvmD31BnCX3EPrR0dVTJEAHdIBu2UoZqXjSWJk46WNL6kDpL6CvSIEaUDJGmKGZ5KmGxBPZ4SdfU++KJRQ4r1UZjaR4vCex9YRzxYog5xMpnj3QhlpI1f9+VR+vIpzikBRs85KQPWlGpSWavNMuj9WGdSuWvR78jLlnzt1EeTUOBcAC3qPz16LfEJbNXGhWwZOUso4vMUJNINAf53pLZK40qqlceK389+h0VklZio1ycwKpKwit/PXYJrKZO9nG+qEECy+jaF9ScJ99gB6Pk0e+oUCf7OF/UIIEVxRB5UDZy++DkoOpkH+eLGhwwyBmCbNg0ccu6RP1RJ/s4X9SggMGOexmO6AbmtPIXGIvZo99VVR3s43xFgwKW5ix3gZVqC1jt5N4Mh6AO9nG+oiEBs8HACGlZcoiip0HZnaUxdKOWwJYrY1MGsl3WkIAVSf8ezqCP0foo6IFBiA1qW1hlpc6scXk6IP8uYJUnc16q7bExvs1qbSUiNB+XleiwukQey1NxE1adfCTAoHSRsh1XD+NlbmgBUbesS9SHFZWaWoVqG8YInb85L3PjUAvRyyCC3nl0FOK52gBTVHWEzl+/3Bw7If8itz/uw5qk87dQAvtnJ7A5JU17cMYFTGsaEmcaGmj2ikt3qIciHLCiErAVleTR76pRnRVpJk8aBhpo9kqIYwho2+E8SWkUfULmy71yIZbqsEiztTKGmr1CjcMUajr4mAdlUy9VD9VdkWbNVA10mr0yM80f+w/MppvF+rDtclYUkkJxX4IC6uqsSPNSXxBSp9krqP89bEeLcGQurTLKNgMI++XvramzMLeVblrGALNXIL0hoDUdaMa505+g7GZ1VqR5pmlbYjcOLnsFKny51DkPpRAhPmrHytz0VF2GucF4OLjsFVYiICXXL5rAHNEY3/65EEt1GeZGA7iHlr0CJQJ87Jb5sLSmQ395ySLNkHEeYrp/VLkvb5816SLNuzQCj5SL/YilV2YE4E0ba3x5TbpIcwbBAGmZD2vnwUlpyh5q0kWafZzF1NULU68+ZOsJaLJFmiE0NMVhkQ/bj2w9AU22SDMLt8EOdnKw5ONehUZdUCtgmjZoK9HOIphvxR6dL/cpW09A7YGZgw0RCNimKB4UXU7p1pT99PM2qq3RoRJeQ41L5MlfPL0SOtuj39ENarmAaRpbwUtYD4HlOIV8ooz6e0O6M/Nw1ArYVhxX34B5Lqwv09KioY9ZLlHUr9jey2oFbMBhbjzTwUuTEPy9ffbLN6sVsAGHuaXYidim9Dm4D3f9ynQQ0JTC3GIbggFyD+eAjCZ/oYH4N46aEDCbAIKAKLZfgEOTv4anCQHzWHgomStHEDgf4cENh6AJASOwItgym+1eafcxl0hAEwEW+TFMvRwa1st3r+xhaoqAJgLMpSsoOEvTY+neXf7od9VGEwAWBSw1BcPMy2dBh+EQ7Q2qzoApM7qTRw+TIUJ6vSJjYeaznb92D35DH1JnwNSNZhqLHiZDkDEwIBYhmXw5LhiGUTYg3/y5uuth5LbRepgMAXHYDg0E6FPZw9bq8hq2MOZ928ojpoGHSUbjsG17UCtfzeoQmGKZZ8kQqqZ99UBgEc5QTlNgidmRIJT0PIZNQN0BUyxLQT1LhnBwihLYlwiiRPP+x7AJqLtrGPBCPUuGyCGrkkWJPn5Ph27UGTCe9dyfZIggjSGsNyGwkmFb8jV1BmzBV8h6kwzhw170PtTRI3bH8IfCQqP1dBCLMNnRFNi4b/W9PqRxAssyyKrEIUy9hm8Y1jRKYBAl6rL6KMgfrNewWaMDFucBDToMsZ32tNTGhzQ6YJBGmeA0Jd+yQcYAXNHogCWQ4JBntPhrPiJjo9CogHkOuXwlUIMoHoGbt1mjApbACsoONglA+ciMw1JjAga1RBMotZF/qVd8gEYDbOcRI8NnwaEjvHSVGgMwBwekd6VeSl3z9361B2vowBwXkflWAuU2YJ/lyB38CuUVDRyYR3fRS3GGQ1pVefwaODDYwxwn4Dfc2fkI/RrnGiiw2INU5XCHkyjFOwgDuPdh9EUDBLYLINYwi0nfglp5wC72x7GcLKBhAbM9Gzy7O0gmpylekxgFaxoGMC9GdphGZKaVgolBaMUhrcv7JU5RvzQIYATSjhiBKUSs4TTawbK/PXoDvll9Bma7WRyFOKSccOiQmx27xHYPxj47vqA7AftIMoTt+lB213dCeqVKgJOXpjEK3AkOgae6D7BWyRCxF6Eg9zN6jQILEOe7JHVQ7Exz8GvWfYCJJENEHpVNboEbpn6IC7mBjzM7cscTTNih7gPsQjKE6/s+PhfpUqGbp74HFqHU+7oPsDIZQjsrmAO0EgLNz1wqH74EctAT1X2A9SUZYoS6D7C+JEOMUPcB1odkiJHqXhPnxydDjFR99nRINUgCG5gksIFJAhuY7g/sD19pUt3px/+5N7Bf/OHkFW94dzewvuXPGndqe8N7uOVTXPuz//t/9wZ2KrEtdphuGE1v+bPGndre8B5uuU7ccmh3kAR245E9Gtgtb3Um3vReZ+BOn68bjuymQ5OSkpKSup+UtbYUaqiuZ6LNl5oGbSuBP+9rztqq2lrkTcw0wbbktenfZfVYr5+DmeChzcCSFzy0+8iytobIXF01NqZoc3Ol6dta4M/72mhQMXVuaBsBa2JpEBtRqK2qm7Bgy+qxXtFCN7WN4KFpprkxBA/tLlINBa0Ed6c1RZsrsEF7LfDn2h8mJ3dpXN86xloZSKytypCyeqxXtDVvOhOkkfChda+1xUgICJoJN99atcCfC5pp+hrpZOAS+XhDvU6htqpuWHPE67Fe0cbSNzPhQ1uQ1xY8tHsI3uINwESbz8hRmWJHpZnGjAK4CmGpL+DsC7VdzFXNKOqxXpGxWayNheihwegpeGj30GxzEzDB5jPCoBb4c1FrHZFP+HUIK920DFOsLUgv6rFekaXBHxQ9ExvthkPrXgtjSd+wgOBoxJpTXrXAnwtvgH5mtytiS1y7ihMLbWtoYm0XYKEU9ViviFzDyDgneCZUuHyKHdp9tNU3loi9o0KBU1WsOd0usBb4c6HtZkMu30vdMkQu4tBdhNpqpNEW8XqsV7QgF7GV6JkgHxckeGh30lJsAU6BAqeKWHO6XSCqBv5c/LsK/Sa0vaCiirbl+xXyeqzXGi+Fz8RyUfkNKSkpKSkpqf4LvLNgEJT2nnby/YIe4xSauFTDNDf6vBUwueL7AFEHxGqDVNLXiKWuLQi5tWnOgOBia64JlvVqy2zo2cqcqTNYtiGNzS39hbW5FdtnXKobUWBLgzBQoJ/pyCQz4tl8Tb4r+mqua8iw5hbMUskPtbkGzvTtFm3MubYE55A1324efQyTEnPxGSrzzK00AGbACqJJFzpmG1gamVPfHvPaWWvFUFXqkwLf35rcv3vMrdRRFNjcgOFNNUhPg2iomWlYqNzzjPBQK8Bmm/WGYwZgW41t8Sn1hQSnfmGZ1IDYWBb0LEUBhy30tCV4mktgsMCokIFSX9PFKHLBQwjWkSWvLynqTF4piGXKr+lQqJsb6GFoq1vkulYCW2w2FjwLw6FmmDoMo3Pdsm4JKZX6qJTCO0vvg1cZPLVL+p25bMlN4S5W6iDeUv84tGG/I+TWlXqYLF0OgYOS2PqLlJSUlJSUlJQU+n+hN9eBduOwwgAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyMC0wOC0xOVQxMDoyOTozOSswMDowMH2Bj6cAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjAtMDgtMTlUMTA6Mjk6MzkrMDA6MDAM3DcbAAAAIHRFWHRwZGY6SGlSZXNCb3VuZGluZ0JveAA0MzJ4Mjg4KzArMBDQ8AgAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjQkMWpXAAAASnRFWHRzaWduYXR1cmUANTEyOTEwNmMzYjM3YmJjOTE2NTNmOWFhMDFiMzRjNjNlZTgwOTdkZGE1NTY0MTJkYThkMzczMzA4ODVmZGRmOQsLXvsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<wand.image.Image: 5129106 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "restart_with_mods/cancer_1/simOutput/wholeTumourVAFHistogram.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAJxQTFRF////ODg4AAAAiYmJd3d3REREu7u74e31w9vr0uTwiLbXH3e0W5vI3d3d7u7uERERMzMzVVVVmZmZICAgZmZmIiIiqqqqiIiIzMzMQ43AO4i9U5bF8Pb6S5HDS5LDl8Dcpcnhc6nQXp3Jeq7SV1dXP1VkDjhUHD9YKkhdW2ZuoqKiQEBAGBsdV2NsVmNsFxsdXWduYmpwR0dHwMDAb26MpQAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHSdkLBwlAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAAC4NJREFUeNrtnQvDmygWhi0LNJMGUSua7d66s3Np99LZy///bwto1O9LYiUeDbTvM5nUL0ZIfKLi8QhZBgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAKHjDACW/21oY58/+yXxTCLl1DRBGCoQlBoQlBoQlBoQlBoQlBoQlBoQlxtOFvT0cfnj2SkiJpws7vnt3evZKSAkISwwISwwISwwISwwISwwISwwISwwISwwISwxqYdyRZYrzfHhldgEIC4NamBBCF5nSpdC9MQgjZYNdoq4yXvh/PBBGCr2wym5apZUknKhKiPcQRgm9MFF3skQvClsYKeTCFLMl1jWEbQS5MG7sU2XynDX9CxBGCbkw09qnvDC67l+AMFLIhcmuNS/PlxcgjBREOhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhIDwhKDPmuKc4lU7e2gFtawmpdI1d4OamHG60Gq9mYQC1NMsOI8pmpnEEYMsTDJSlWbMVUbufXUEAvLmbTSkKq9HRscwyqGVO3tIG8lasFapGpvB4Ew1dpDVXk57cqlcv8gVXsjCIS1td0RFtW92RBGCoEwzhXL72uBMFIIhDW6EFnd3psNYaRQNDoqrjKu7s2FMFIQrU8Mmi3Mgi1sHwiESVZD2G6QtBLXzIawMAiEnY2amw1hpJDsEh13y4EwUgiE5dKR35sNYaRQtBIVF/x+MRBGCkXwVwteMrQS94Ei+FvapxqxxH2gECbsU4lY4j5QNDq0EYXGLnEfaBodM4EOCKNltTDOcB62J6uFKYnzsD2hCE257ICqz+Gout0jcuu3giI05YyUfTmiEOKM3PrtWC+Mdahpecit34z1wmTFJ4cwoXU96QZdcv4GwiihaNbnk0ZHIxtToxv07SA5cdZGsyFx1EU+kFu/GSShKcnzoi/HHspKgdz67SBJEZBi0MIKrc/Ird8OikRSkzNxaRMq2RWI3PqNoGh0yOzM76bWQxgt1JGOayCMFOpIxzUQRgp1pOMaCCOFOtJxDYSRgpshEoMia4rjZoj9IBKGnI69oIgluqca9zjvA9UxrERe4j5QnDgLIRh2iTtBdTMEmvU7QSFMIM1tPyiuhxVIc9sPAmFVPTsbwkih2CUa3JS+HyS3zNpmosDllX2gSBHALnFHKLawl70ISKRqbwl5LwK5sZNI1d4M8l4EeGmFIVV7M6ivh0nt+mlGN+ibQd2rtmkyK2xI1eZ2ZwlhlBALa20TnxVnpGpvxmphVaUmRUh7Cs1qhVTtzVh/j3PZvDJid4lhqdq/P52Oz14PybD+HmfDXt2U7huMIanaH7CRLYcqkfQ+EEZKDJ2DQVgAMXQOBmEBxNA5GIQFEEPnYBAWQAydg0FYADF0DgZhAcRwMwSEBQBhiQFhiUFxHnaenQ1hpMSQhANhAcQwlAeEBRDDUB4QFkAMQ3lAWAAkPeGUfNWQihAWAEUrkXGBWOJe0Az4Jma0QBgpFLcbCStMYAvbB5JofaHNJTG74b5jt6DceggLgKTP33bUYzjXPDC3HsICoBDWFKyY7BFtEyQstx7CAiA5ceay1ZdyKq7bwNx6CAuAaljgQQsXrBpz6ysh3kMYJUQ39KliLKfVgd2gQ1gABMNRvcj8VfYVE9gNuhX2h8Ph8OxVkQYEw1HJaSyRGWNbG2G59R/8492zV0UaUF9xznt1gbn1ELaUGLougrAAYui6CMICiKHrIggLIIauiyAsgBi6LoKwACLJmoKwpdB3XfQaCCMlkqwpCFtKJFlTnbA/IkD1VWIYGWIQdsB29lViGBkCwgKIYWQICAsghpEhICyAGEaGgLAAYhgZ4oWw0+n0w7NXSsxEcgfmKMyCpv0M6684c07arIeweYiEGYJIx+n0Jwj7KiS7xMroNQNv98IGIGwGkqHt9ZCp7XLrXYmhufUQtpT1u0Qx0ZVluq5ZFZ5bD2FLWS1MvsxLtJbqB3LrIWwpq4XlL/MSM3cDy5hbvzRVG8KWQn8eVtt9oQi9GQLClkIurNbnLDy3HsKWQi3M+wrPrYewpRAL800Q8VhuPYQtgXrsFdf+cNvYA7n1ELaEiIK/ELYECEsMCEuMSIVB2z0gLDEgLDEgLDEgLDEgLDEgLDEgLDFiFPbn0wnC7hGjMMQ7ZoCwxICwxICwxIhX2F8Oh7fPXjsREq8wCwa8v4ZamGr8fSyrUrUhbAZqYcLfx7IuVfuesMPx+PHZ6+vp0O8SnbB1qdr3hB3RN/BWwoK7QYewpWwjbEjVdl04Ewg7nk5HCPNsI2xdqva1sJNXBWHZVsLWpWpD2AzkrUSCVG0Im4Fa2JkgVRvCZog60vHXS3QKwgaiFvbhspFB2ACEJUbswn48Hg8QNiF2Yd1GBmEDEJYYEJYYEJYYEJYYCQibPiAMwhIjMWEY3DQxYRijBcISI0Fhx+Pf/OP41i7tp2gSTt8eqUrakgSFDXT3kb2jyuo+kJW0JRA2AGEd2wnrbvzrJj761KoffXZVAB8nC0BYx3bCpnThkOAztWnw5LsWRpNbD2HXbCOMKLf+QWHHm+Ze6TjNHBZj1raNMKLcegi7Zhth09z6n97wOX7+5Zdf/eMGv9563OZnzj/52Z8uRX56VZF90+fJn59uFvn56n2x8f7vm/wMLsIk5//46WWN//z8+V+3Hl8uj3Hq+k2/3Xzfl4VF/ra86i+zVd/+lP+er/rLbNV3v/h/rKTJ6vvv/7YQNs2tv5K56pew7nf0XVa9hGlufUQf/fusegnT3PrX8KCSKBf+TqtexphbDwAA4AYV79shk8DVYoZl+jEAH6zaHhCq0IXHj6seqHuo2p7iBB8nZFu9LmY/uLb/+a89CVwtZVymGwPwwaqzrGWhba6x6obVvHy06obxkqmwhSUz4uob7IZt57fGV16MgavF33tYphsD8MGq/doPXHis2jywysZvbevVbejiXLz+BrtxZrn9xbipaeBqKS+WCV1xY9W2DB4qbKhaMcGKwL3aWLVkLdeBW9ggbPINdsPV19UpHhA2XaYO3Z+OVbdFFixsjLOxUtWBP/Ox6rMpiuJRYWMx+6HsDryrcy5wdY/JMt2Yco9VrY0woceCoeqcyeC1NlZd2mKK4H1qL2wsZkdM6w4HuZwNXN2jX0bmD/i6VK3OGbd7RBN4+ByrtrviKnStDVW7A2/4QdALc1X7YvalYoVd1/Z3Mhe4uke/DJP9GICPVN39WoN3iUPVWaMFC202DFVLbYwJbiX6AZftoytmZ5QbOzh3JzKPBK663iTyYQzA8KqVX0yF191X7T588EForNp+8tBl827AZV+CDD51BQAAAMAuVC5iK4d29FlNp240r7kIjgoDSoQLERfDadN4pu6mroWdWRPeAgSECGPPVxlzcnJVKSO4kqXguZ+yLlUt7DkVb8vadwNeG/tqW/O8Fe4FLniV2be5triyf+SZnVm6riBLUUm7+Hnv2Pk3j+Ba8doKc2E4kRe1zKuqKWo/JbJcl43mGSuawl0xydvCvqrbpi6a2mStruyiUvg8WFPKos6EfUlnleYNlyz3AStAieA1N81FmN8RqpaX3ZTI3EqvjLts0fjQh5Pj5rBWSib9VC/M/ilb5mbmrN+zFq2LKwJShF3TxtmyZppOU1HyehTm498unvtSWM05V64BYjJZdMJcAmc3sxdWmb0vTn0H2FVbtl5Yk/WarJ1hCzszH+G7EuZi+SorS/sGvxUafvZDJ/TCXAg3t/vT8EuR4Ct024JrdOiiEC56eha6MP2Uu0pT6OZaWGPfzbKznantJmlMwTPOhOmNZ8q+4hZmiOxR0514daHmvI/eDlPu2YdT7f/+D//sF8llH+YV/lkp907ZzfSludlockTIzGWYIvxaP9icmcswD1xlAQAAAAD4Zvk/wYZM6yrgIgMAAAAldEVYdGRhdGU6Y3JlYXRlADIwMjAtMDgtMTlUMTA6Mjk6MzkrMDA6MDB9gY+nAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDIwLTA4LTE5VDEwOjI5OjM5KzAwOjAwDNw3GwAAACB0RVh0cGRmOkhpUmVzQm91bmRpbmdCb3gANDMyeDI4OCswKzAQ0PAIAAAAE3RFWHRwZGY6VmVyc2lvbgBQREYtMS40JDFqVwAAAEp0RVh0c2lnbmF0dXJlADQ5ZDliZjc5NjI4YTc1NzE3M2I2MTEwOGUxNTNhNDgwMTNlYzNkYWQ1MTNmMjU3NzRlYTNmMjIyNGM4ZmRkMWQXxtbCAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<wand.image.Image: 49d9bf7 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "restart_with_mods/cancer_1/simOutput/sampleHistogram_86_78.pdf\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEgCAMAAADi9tZbAAAJJmlDQ1BpY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpNzTVQAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAJ9QTFRF////ODg4AAAAiYmJd3d3REREu7u74e31w9vr0uTwiLbXH3e0W5vIQ43AO4i9mZmZ3d3d7u7uMzMzzMzMERERVVVVZmZmqqqqIiIiICAgiIiIS5HDuNTnl8Dc8Pb6pcnhc6nQc6rQS5LDXp3JT5TEV1dXP1VkDjhUHD9YW2ZuKkhdoqKiQEBAGBsdV2NsVmNsFxsdXWduYmpwR0dHwMDA3TaNUAAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAAHdElNRQfkCBMKHSdkLBwlAAAKeHpUWHRSYXcgcHJvZmlsZSB0eXBlIGljYwAAWIWdl22SZKkNRf+zCi+BTwktBwSK8P434ENWd0+33TNjOyuIzHoPhJCu7hXpn+7pH3zqUEn5fVbJPz7167ccudq1jtq115rHHDZWzX/2SVKkadPcy8gjd//TmX/xCXZ9Hv1w57R6/h9DH4/+x/lLugxt0r758u0E6omDZa3aP/8XnV8v6lQlQvn78/XNTulSZf/xfPfvzxPh/ITx63+fPxboz8+P/2Ho5+eRfzbUycyXqzV/7TCzY+j3z/9kfvr1zN8/tfbyDiwmwvtJ+puECMdV4Y2MmrV9h0a33lJvTCJKbxrLXMbvo/x3ptN/2v5vTf/6+dv06zv6/JYlPh0/yJqxkYkKb9j+efTXcWi15VYakP1diUQbD8zlu0eliPvf1dL3z+/mSaz6OVqb8RWHZr+fWM3e99b5mVfmWf8+72Oo9m/IjfmJxRYPED/Ikvxi8Uek8jP4FsUDI8MwVC6m2isLBkVL0jJ1k9v+WtlZ9HbqLBo8GHg3WPOwJ/MRDil5R1N9RQc8CdrEg4mBdxLDgGHAMLAwsHi4MLDrOySDNc4aZ41vDD3mOCw6GGBevvy+++M1TMPY5OX9KeOQmsYwRuRSB4P3DY9Km4zLUXkIsRWyXnC/YKMIi4V3yju8LhMjeFyMOXhboNaCp2UXDG1+4GJxvg/fh+/L9+U7WBCL4mwMh4Y741AvwghCO8lUYXA0qpnBS3avykNlIdmr8+ZqTCTHdWFks5gNq29yMnJ9OSIEFei0l/6WN+AVklXyo9rGLtQbI3KDd5rwTvFJL4Djf+N/jDcC3zb/u+Z2Goaw3K7nFka2hcJpmfphHApr594nCEAXSHfH447BPp36XqCCd3javafcDxOIyYNJjwvUTh7F8yAboy2gA9zHzIOjD6AygMjAq7EYG+lxxhkJbPGDNH/+OKJUzY/IBU+E7ImsLLrBnmexk2VFFn84LFluo9DgnKwpK5hQdtd24IzIVD4Y7VnZWakxJdC6eX4gLjbVmFDrBr+RJ1Uwu+Q5VgLMN084ZOLuXAtg8z+L5tU8AaMBXgN4xjGNjUx6NrVsk98g3gi4eaRs7GIsWKXkxbEWni0gsTjSomwWEFhkaBGLhZqseHnmD0Ld0MWGk7ZQtJu620ze+5UP3wR+k0EvQLCu7EDBh2cH3Q62fGn2V2YA1zF63l9Fsk9/pbbyIS6HiQfIH2fC4TfxuMDhgr5L9i7Huhr52qYcJV9CcO+lLPEoOH8A84AaAlQHsYrdUOPIcV95E6VKBjqMK5xfcdk2bvP86FtYKOTE4LsHfHtKmV7KIlpupdzJ4bRQV6X2Uar0QumUulqpzriQ+SP0ykDXCuIIATAWmPYBEQxKU0qn8Ho3RHqVPnfp60AOlz0hh1LLaHRCQwqyAVnsVMY+hVO9ait0CEVYLOJFZhTZFUd5Fqso1KC9FJVBr2FF1y1gq2homQVDFHqZvJxzlbkCYuc3Cz+Uw5FMdjFOahvonkNj0suqqyxCs1Sho1uARiqLgOJ42W2XzTE3Bjee7LPKYyAgUHzwrbs48XH34gT4QFqHKj76KMwSHUsrB2O3SLl4d4nJtV4ugLrXSpCNaLeE8JvnsaPEXfVDpcSewqvAPIE6SAOyI1UQ4OTQbL+Ipt/Kqlqr1jpGrZOfK2o9B81ZFd6qcFVt1mvzmmqLx5ZRez90Eo7G7drPetVVB5OHMJD64YxAyetTc8bU17xVuZP84pF2q6pUGQb0OOp26mxB8wdsFo6cXu2JLUYJPKJ7KmxC8eAgbcxio0X6oeOARGrdTaBlq5uJIKI+avNm1eVWx6AfhTO9HuJyVOph43PBJaC53VPFMzhcKzVTOSBcvmpYqcFRImCuNmAvim9RvWdTB0C5kz5CVDbfURu+pValtWob3u+Nma1Bzk2jtT1bI2UdX+mRWrfb+pl0Mq0N+HlM+jOvbcShODQ1UYK/bpNriEVv+kTDvOnRNktvNCBtTm/T52tWPkkyNrLNwQO6w8zSnhpHRVmiceK2BViu1fadZFQbbV9zjuS3tVNro1oaOG0wTLso0mXTiyLBJIn8lBZMoFlqcSvK2KjZ/ijykQ+hBYVCRS8HpRd/UCpcr3sQUCUe7KSHrhaJ6shhpx3tc3Uq/JEGUkZDDSmPc+nSa389oazdJZA2oqS6gR0Sh2BNJLtTyH1Cj0blmBDTZZ1OhrxoX3o6jvQN/Dfx3hjeeE39dZLafa8OpDqzUj9GMo73SxNw5Xag8KWVtMrEssd5Qg9hKxex/ageqkAKoYNBYQ5AMCqXGlCnA1ob5BFhXYOAjd6xSmPZz6bK5hjKQZ1qgVcFaZVlgy55EIyhVBIqnsYEglPPmL6HwTImBuEheVnHYtlajBhjE7VtjIvNxoDE/Mg4eHt0pnHcBtQ0rvi4+wwoHwUvAwGg1cIJLqwIG844/MubBY3iWCWi1bjkoOCPswV0SUNb+ku6denXQA9bGUV+VYTflKBQ5YKsixoYZg6FLaizzOvyLjVitsTiIWVy9KBHUNnsvBffEfip4otrK+J+6DHONqFW5cqW66CBiAdHk4DTaccQevqWS24AfLGh9AgkmGpeOEIH2YgE9QdC+9fd0skSZEPnrsQmvXOpwOwSXD9pgnQ3BAah4Lo+mWx1qU3ahgtrcbEksTQ5XeF33dQRvKo+MeRPVbjfUEP6+tcLBV4mwA50MF3j0mV1LrtrvpZiolGz+IFEMkwHAUeHEjRNqhT9PBOsz34pdhaNtemOXnQrgeGW9c5kMbE4pxhkcKdB2mb4GndSlmkuXxOpn8Rw7vDpAmPw7EBdhzUnYt5Pcu6MhmwafTO9G+0a3QbSQvNZ1kyGfEDay9DyVywGl0A59FSToqNOxggbbp8yJL1GB2UE04iDze42N47VnvAum4UDgmnrAGq4fq8wZNCcOR5qB4ShQobu2V0XtBwOui2CFk9ob89MdAiKtAr0zjBZEDSFz0ApO1VFmVOAc43FXrQqBGCBGVB2F16tiZBM2uMFwTLFaGZ8LUQfRVmbMtvXkHRfTid4Or0IWn7RjovsP/zi0X53O0qSrmulTRuyy0GwOorvMH0j9utyQurUqOTS9piL/gy/1TbEBujmxhtKm/I+3Gbgo20shqX32gNLlx8PZ2W77dfw7ENrywmgcTgtUH6UNIKmklYyXzoKURqHlmCZQPWQBIikHS4DtP3QrY++ORlo6Fz9nRtHfw0J+GjH53ZHP9jLaFCmE4vksIVvbrFYcg7iKJbDZwiH+H2326YeHIDbzMmbtq05h6ENbXG4LR3Y/iA3iTgafkBE/Z5xiNYYRw4sjj3icKYgixdsCg0xeSddZ8Um9jS/3EJ8LtqvnA4zkHA/tDwnaA9icbNBLvPmcee64/Q3Axk7GyfbhbsuMnJ7OFUIzedzxSRd+OICACSRNmA7PRbYPyQUUl0X0oRcNvGGWi997z3mdAnzktcbKF84ffSYie57RKFfKBH0MoSkWEBJ0REQdAe2hnvPDZET8pJGozmZMwEdrQ4loAGzpFi08ls1yCeFMomgxaFGbt9xj8ORlG1E+hftkQTIS62KtQAACyVJREFUeNrtnQ2/m7YVh6kmKa7jFYEFqHRb99amyV7SvXz/zzYJuOBc20wyB8zJ/T9x/OOCQbYeg8VBB2UZAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACg4BsBKPnN2sKkfPZX5qtC6bVLgDBSIIwZEMYMCGMGhDEDwpgBYcyAMGY8Xdi7w+HbZ1cCJ54u7Pj+/enZlcAJCGMGhDEDwpgBYcyAMGZAGDMgjBkQxgwIYwaEMQPCmAFhzIAwZuxB2G9Pp+Oz64ENexD2HXayeCCMGRDGDAhjBoQxA8KYQS0sl1KaLDNSFsMcCCOFWpgqlTpnxlbKDsYgjBRyYd32ZJllNu/nQBgp5MKsrYus8pLUIArCSKEW1ujG1Z2sTliu1PcQRskKrcRWZXWNPWwlqIWZLKtUlruiEE0/B8JIoRYmSmvPWVE6Ww9zIIwU8j1M9xvU55c5EEYKIh3MgDBmQBgzIIwZEMYMCGMGhDEDwpgBYcyAMGZAGDMgjBkQxgwIYwaEMYNAmGmzXFXFvcUQRgqBsLbOnCzze4shjBQCYVIaUdzXAmGkEAhrbKmyur23GMJIoWh05NJkoUP9bSCMFLQSmUGzh/UpK7eBMFIIhGlRQ9hmkLQSlyyGsDQIhJ2dmVsMYaSQHBIDd7cDYaQQCCt0AKGpbaBoJRqp5LQZjRznNaEI/lolK2GGvwrnj47IcV4NiuBv5Z/qFy3enUaO83pQCFP+qRpiidoWXhhynFeDotFhnSqt6addk3lhyHFeDZpGxxjoaIVSojwjx3k1FguT4vI8THt3ojbIcV6NxcKMfn0e5t0hx3k1KEJToTmYjznNWecOOc4rQRGaCkYqhKa2Ybkw0WPuLYcwUpYL07mcCyVCGC0UzfoCwd/tIDlxts6K873FEEYKSWhKy6JEo2MbSLoIaJWhI+lGUHQkdYVQFl21t4Gi0aGzs7zrC8JoWSHS8QoIIwWRDmYg0sEMRDqYgWQIZlD0mpJIhtgOImEvfTpuAGGkUMQSw1ONE+dtoPoNqxCa2gaKE2ellMAhcSOokiHQrN8ICmEK6UbbQXE9rMQV5+0gEJbXs4shjBSKQ6JDUvp2kKTM+maiWnZ55Yfj8fjsquABRRcBgkOif7x/dlXwgGIPI7iLAITFspO7CEBYLNR3EWhk170jNSkdwmKhvh7mpLQyPSkdwmJZ4QKmVOlJ6RAWC7mwXNo2PSkdwmJZLCzPzRebkErkU1K6lvIbCKNkeY5z1bwy0tr0gbchLJblOc5OXDbrjTfk0gfehrBYqDqSviCc862N5KR0CIuF+uZgxXBGlpiUDmGxkN8c7AoII4X85mBXQBgp1DcHuwbCSCG+OdgNIIwU4puD3QDCSNlJMgSExQJhzIAwZlCch51nF0MYKeiEw4ydDOUBYbGgEw4zdjKUB4TFQnInnEouHVIRwmKhaCUKqRBL3AqaAd8W380NwmKhSDdSXphavocdDr87HA7PrpC9QxKtL62zyxsdPc+ukL1Dcs/fduqYfQ2EkUIhrClFefeICGG0kJw4S93a5SfOEBYD1bDABK1ECIuBKKHPLLyrNoTFQjAclSCKJUJYDATDUWmiWCKExbCjK84QFsNObl0EYbFQ37qokV03+0dynCEsBupbF9m6FvljOc4QFgP1rYu8pfrBHGcIi2GFWxc5+ViOM4TFQN9rqvbHwuSBtyEsFvJbF9XW72sP5ThDWAzUvaY6X4/lOENYDMS9pjp56rEcZwiLgXhkiE5e2MceyHGGsBh2MjIEhMWyk5EhICyWnYwMAWGx7GRkCAiLZScjQ0BYLLgexozlV5ylJBjwbRIWugA/u1L2DJEwR3YBE8LmITkk5s4uHHgbwmIhGdrezvTUfkDY7w+HPzy7XnbL8kOimtX1iDDsZDMsFqaJ+yVC2DyLhRXE/RIhbJ7dnYeNwo7HI7Rds19h/hkjVF0DYcyAMGbsWli/ANougTBmQBgzIIwZEMYMCGMGhDEDwpgBYcyAMGZQCzNN1x3n8RxnCJuHWpjquuMsyHG+5Ds4u4L+kBiELchxhrB51hG2IMcZwuZZR9iCHOd7wg7H44/Prqw9sI6wBTnO94T5152eXVl7YB1hC3KcIWwe8lbi0hxnCJuHWth5aY4zhM2z00gHhN0DwpgBYcyAMGZAGDMgjBn8hL07Ht9yaJGfsMP7N32rAQhjBoQxg5Mw3IAlgzB2QBgzeAj74+HwJwjr4CHs8gFhKwNhpEAYMyCMGRDGDAhjBoQxg62w9yR3EON3lQ3CIOwVEEYKc2HH4/HdkjfnS//h8gr20u2tD3NhS/ezvvSpBUOz364JhEHYKyCMlJWELU1KjxL259Opr+DT6fSX7pGQ8nc8ncarbO/9Bvq5/Vb/ejp9++Wrf+xevQfWEUadlH5b2OV+lpxae+rah8dxk/3cey3Q/bQm1xFGnZQOYSPrCLtMSv/pGznHzx8+/NI9bvDLrcdrPko5vvRnGc2nDx8+vZQe6OdebvXVu/wUv+01+f5vq3wNXoRpKf/+05cl/uPjx3/eenx+eUxT1y/69ebrPkdu8tf4oj/PFn37Xf5rvujPs0Xf/eD/9pIuqu8//11D2GVS+pXMRd+EZd+jN1l0DJdJ6Tt662+z6Bguk9JfI5O2RLnyGy06jikpHQAAwA1yObRDLgJX0YzrNFKmnzuORfsfhDx15entmgfKHov2pzjJvxO6zV9vZjuk9f+6j30RuIplWsfWtUit87HoLGtFaptrKroRtaweLboRshImbWUtnLr6BJvh2/mt6wq/CFxFf+5xHV91dWqdj0V3tZ+48lS0e6DKpk/ty7Vt6upSvf4Em3EWhf/GhKnLwFUsX6yTWnFT0X4bMlXYWLQRSpSJR7WpaC1aaU1qrQ3v9uITbEYory9TPSDscp069Xg6Fd2WWbKwKc4mKlMnfs2nos+uLEuTWmvDu502sx3GH8D7MucCV/e4WKe2qb/dU9HWKZf6WzAWXQidXGtT0ZXfTJl8TB2ETZvZENeGn4NCzwau7jGso4sHfL0Ubc6Z9EdEl/jzORXtD8V5aq2NRYcf3vQfwU5YKLrbzLbkovR17b8nc4GrewzrCN0NQ5x6VOuL7r+tyYfEseissUqkNhvGorV1zpm0lXU/4LJ/9JvZGBPGDi7Cicwjgav+NoxFNwxx8upd0aZbzaSXPRQd3rx57FOHQv07T1236Adc7ragk09dAQAAALAJeYjY6rEdfTaXUzea11IlR4UBJSqEiMvxtGk6Uw9T18LOojHPfstvG+X8+aoQQU5hcuOUNLpSsuimvEtTK39OJduqNln4y/m5bS2LVoUZUsk88y8LbXHj/ygyv7AKYyhUKtd+9fPWsfOvHiWtkbUXFsJwqihrXeR5U9bdlMoKWzVWZqJsynDFpGhLP9e2TV02tctam/tVter6wbpKl3Wm/Cyb5VY2UouiC1gBSpSspWtehHUHQtPKqp9SWaj03IXLFk0X+ghywhLRai10NzUI83/qVoSFhRiOrGUb4oqAFOVr2gVb3kzTayorWU/Cuvh3iOd+KayWUprQAHGZLnthoQNnv3AQlrutL069AXzVVm0nrMkGTd7OuIedRRfhuxIWYvkmqyr/gm4vdPLcjTk4CAsh3MIfT9MvRYL/Q78vhEaHLUsVoqdnZUs3TIWrNKVtroU1/tUiO/uF1u+SzpUyk0K5wXhm/JywskBkj5r+xKsPNRdD9HacCs9dONX/7/7onrtVCj2EeVX3bEx4pe4XdlsLi9Hk2CEzl2HK9Gv9YHVmLsM8cJUFAAAAAOCr5X9A5GoDI5EdUAAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAyMC0wOC0xOVQxMDoyOTozOSswMDowMH2Bj6cAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMjAtMDgtMTlUMTA6Mjk6MzkrMDA6MDAM3DcbAAAAIHRFWHRwZGY6SGlSZXNCb3VuZGluZ0JveAA0MzJ4Mjg4KzArMBDQ8AgAAAATdEVYdHBkZjpWZXJzaW9uAFBERi0xLjQkMWpXAAAASnRFWHRzaWduYXR1cmUAYTEyNjQ0MGIyNTdiZjZiMzFmZDRmNDExNTc4YTk4NzNhOTNhNjYxNzUzYTUwZDRjYWE2N2U5M2VjZmRiMmFiNQwTnBoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<wand.image.Image: a126440 'PDF' (432x288)>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pdfs = glob(os.path.join(restart_with_mods.outdir,'cancer_1', 'simOutput',\"*.pdf\"))\n",
    "for pdf in pdfs:\n",
    "    print(pdf)\n",
    "    display(WImage(filename=pdf))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the mutation matrix\n",
    "mtx=unpickle(restart_with_mods, 'mtx.p').toarray()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x7f37961dc668>"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARgAAADzCAYAAABdegl5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deXxU9d3o8c93JisYBAQRcFdA1qQV2Qn05QIiSsUNu0irrdZbb9tre2+1z/NcvbePt8tT6+tpfbTF5aVt3XBBUUGgthDCKmDIwr6Dsq9hScIk3/vHOTNMICGTyRzOzOT79nVemTlzzskvAb7+9q+oKsYY44WA3wUwxqQvCzDGGM9YgDHGeMYCjDHGMxZgjDGesQBjjPFMht8FMKa1kk45Sk1dbBdXnpylqmO9LVHiWYAxxi8n6wgMuyimS+tmbe/kcWk8YQHGGN8IEhS/C+EpCzDG+EVAAhZgjDEeEEAC6T3OYgHGGL8IBNK8ieR7+BSRsSKyVkQ2iMhjfpenMSKyRUTKRKRERJa55zqKyBwRWe9+7eBzGV8WkT0iUh51rtEyisjj7u99rYiMSZLyPikiX7i/5xIRGZcs5U08QQKxHanK1wAjIkHgv4CbgT7AvSLSx88yNeFrqlqgqgPd948Bn6pqD+BT972fXgFOH8pssIzu73kS0Ne95zn3z+NceoUzywvwjPt7LlDVGZA05U0stw/GAox3BgEbVHWTqtYAbwITfC5Tc0wAXnVfvwp83ceyoKpFwIHTTjdWxgnAm6paraqbgQ04fx7nTCPlbYzv5U00ASQYiOlIVX6XvDuwPer9DvdcMlJgtogsF5EH3XNdVHUngPv1Qt9K17jGypjMv/tHRKTUbUKFm3TJXN74WA3Gcw395pJ1B6zhqvpVnObcD0Wk0O8CtVCy/u6fB64CCoCdwNPu+WQtbwsIgUBsR6ryO8DsAC6Jen8x8KVPZTkrVf3S/boHmIZTPd8tIl0B3K97/CthoxorY1L+7lV1t6rWqmod8AKnmkFJWd4WEZCgxHSkKr8DzGdADxG5QkSycDrxpvtcpjOISFsRyQu/Bm4CynHKOtm9bDLwgT8lPKvGyjgdmCQi2SJyBdADWOpD+eoJB0PX7Ti/Z0jS8rZEeB5MLEeq8nUejKqGROQRYBYQBF5W1Qo/y9SILsA0EQHnd/a6qn4iIp8BU0XkAWAbcJePZURE3gBGA51EZAfwBPBrGiijqlaIyFRgFRACfqiqtUlQ3tEiUoDT/NkCPJQs5U24VjCTV2zTb2P8kdG5jbabeHVM1x6cUrY8anpEyrCZvMb4RUjpDtxYWIAxxie2FskY4x13FCmdWYAxxjepPYkuFp7Vz5q7iDFqdmxKSLXygpU5GdlM3jjEuYgx1f4ipVp5wcqcVEQgEAzEdDT9LLlERP4pIqtFpEJEfuyeb/ZqehG51t05YIOI/EHc+RnuHKS33PNLROTypsrlVQ0m1RcxGnNOJLAGEwJ+qqq9gSE4y1n6EN9q+udxAnsP9wiveH8AOKiqVwPPAL9pqlBe9cE0tDBtcPQFbtX3QYAAwWuzySVP2qfMpJxUKy9Ymc+Foxzep6qdY7pYEtf8cReyhhe1VorIapx/hxNwJjOCs5p+LvBzolanA5tFZAMwSES2AO1UdZFTRPkLzgr8me49T7rPegd4VkREzzKZzqsA0+TCNFWdAkwByJP2OpARHhXFmHNnLh9vjfVaAYIebMXgNl2+AizhtNX0IhK9mn5x1G3h1ekn3dennw/fs919VkhEDgMXAPsaK4tXASb9FqYZk2gCwdjnwXQK76TomuL+T7r+I0XOA94FfqKqR9zuk0a++xn0LOfPdk+jvAowkUWMwBc4bb1vePS9jElJgjQnwOxraqmAiGTiBJfXVPU99/RuEenq1l5iWU2/w319+vnoe3aISAZwPk1sGOZJJ6+qhoDwIsbVwNQkXcRojH/EaSLFcjT5KKeq8hKwWlV/H/VRs1bTu82pShEZ4j7zvtPuCT/rTuAfZ+t/AQ8n2rl7qc7w6vnGpDpBCCRuqcBw4NtAmYiUuOd+QXyr6R/G2S85F6dzd6Z7/iXgr26H8AGclslZ2UxeY/wiEEzQUgFVLabhPhKA6xu55yngqQbOLwP6NXC+imZuSWIBxhifCM3q5E1JFmCM8YmIeDJMnUwswBjjI6vBGGM8IUIiO3mTkgUYY3xjTSRjjEekeTN5U5IFGGN8YzUYY4xHnBpM6m4mFQsLMMb4xObBGGO8IxLTbnWpzAKMMT6xGowxxjs2k9cY4xWrwRhjPCPizZaZycQCjDG+adaOdinJAowxPmkNa5HS+6czJqkJGcFATEeTTxJ5WUT2iEh51Lm3RKTEPbaEd7oTkctF5ETUZ3+KuidhSdfAajDG+MaZyRts+sLYvAI8C/wlfEJV7zn1veRp4HDU9RtVtaCB54STri3G2fJ2LM6WmZGkayIyCSfp2j0N3F+P1WCM8YmTVSAY09EUVS2ikR3+3VrI3cAbZy2Pk3WgnaoucjfzDiddAyfp2qvu63eA68O1m7OxAGOMjzICwZgO3LxIUUdzcnaPBHar6vqoc1eIyOciMk9ERrrnuhNj0jWc2tAFTf58zSikMSaBRBKbF+ks7qV+7WUncKmq7heRa4H3RaQvCU66BhZgjPGNQLh24t33cBKkTQSuDZ9z81FXu6+Xi8hGoCcJTroG1kQyxj+SuD6Ys7gBWKOqkaaPiHQWkaD7+kqcpGubEp10DawGY4xvwp28CXmWyBvAaJy+mh3AE6r6Ek5ytNM7dwuB/ysiIaAW+IGqhmsjCUu6BhZgjPFNIptIqnpvI+e/08C5d3FyWDd0fcKSroEFGGP8I4mrwSQrCzDG+MSpwaR3N6gFGGN8IlaDMcZ4yQKMMcYTgng+D8ZvFmCM8Yk1kYwxnrFOXmOMd6wGY4zxirPptwUYY4wHrJO3CSKyBajEWc8QUtWBItIReAu4HNgC3K2qB1tWTGPSUGJ3tEtKiehh+pqqFkTtVfEY8Kmq9gA+dd8bY04TrsHEuOFUSvKiCzt6a71XObXlnjEmirMnbyCmI1W1tA9GgdkiosCfVXUK0MXdVwJV3SkiFzZ0o7vl34MA2eS2sBjGpCIhKKlbO4lFSwPMcFX90g0ic0RkTaw3usFoCkCetG9y4xpj0o0gBCR1ayexaFGAUdUv3a97RGQaMAjYLSJd3dpLV2BPAsppTFpK9xpM3OFTRNqKSF74NXATUE79rfUmc2rLPWNMlHANJpajyWc1nHjtSRH5IirB2riozx53k6itFZExUecTmnitJfWzLkCxiKwElgIfq+onwK+BG0VkPXCj+94YczohYQEGZ5vLsQ2cf8Yd5S1Q1RkAItIHZ8vLvu49z4X36OVU4rUe7hF+ZiTxGvAMTuK1JsXdRFLVTUB+A+f3A9fH+1xjWguBhPXBqGpRrLUKnJHeN93sApvdfXYHufPa2qnqIgARCSdem+ne86R7/zvAsyIiTW38nd49TMYkNWcUKZaD+BOvPSIipW4TqoN7LpJEzRVOsGaJ14xJF82swcSTeO154Jc400l+CTwN3E/jSdQs8Zox6cPbYWpV3R35TiIvAB+5b8NJ1MLCCdYs8Zox6UKE5jSR4ni+dI16ezvOKC84I72T3JGhK3A6c5da4jVj0kriajANJV4DRotIAU5TZgvwEICqVojIVGAVEAJ+qKq17qMs8Zox6SKBo0gNJV576SzXPwU81cB5S7xmTDqwpQLGGE+l+1IBCzDG+ETEajDGGA9ZgDHGeEIQglgTyRjjEavBGGM8kcjFjsnKAowxvkn/Tt70/ulMk5Z/e2O99wtuq/CpJK2P10sFkoHVYFq5a/96Vb33w6f39akkrZHVYEwKK7q5jPm3lDF/XFm985W/C1E0tjTyfuk96wAoHl9O0ZhSzLmTwB3tklLqltw0qXBmf7ROGTmjP5v/5dTK+ryfZVD4yQAAFt6+muP7j7P0nnXU1SqFswb4VdxWRxACMf6XqlK35CYmHa/oSNHNpVzxVMczPjvw6yqGTetNzvk5DHqrJ4Uz+/tQwtZNRGI6UpUFmDSx7ud7zzg376aV9HvuYgpn1q+VLL7DSV/V8bEcAFSVZd/c4H0hzRmEQExHqkrdkpt6dpfuOuPcqNn5LL5jDcXjnX2Glty9FoAh715T77qh7/Vm4GtX1xtBCvfbLLl7nVdFbvWsiWRSwuK71pKRW39AcOX3tjJvzEr6D+/BiI+c7T0GT+3FgltPBZEld62NvF4wYRXtup/PZ99wajIjZ/Rn/rgyBk/tCUDZf4veC9okhFgTyaSAIW/3irze/sRhAI7srOT8S9rT9lFnDsXanzkJNod/2JeDv6lmyd1rGRx13/AP+tD/+Yu57vWrmT+ujIUTVpGZm8nCr68C4OjOynP147QqiWoiNZJ47T9EZI2bVWCaiLR3z18uIieiErL9KeqepEm8ZpLA59/dzJK711J7so6SB7Zwyf85H4C6k7UUvHgZK+7bxLwxK9m/fh8LJzjBosPPs7mkX7dIk2nBbRUsuXsdRTeXMn9cGcHMIIGMACdPnCS7XQ4Lv76Kuro6Ft0Rc+pxE6PYwktMNZhXODPx2hygn6oOANYBj0d9tjEqIdsPos4nNPGaBZgUV3OshsFTe5F3UR6Vu45SfGsFc28siXz+1b9cyahZ+Zw8cZLak7WR893+dx41R2sAyGmfC0DhzAHUheoY9kEfhrx7DSNn9Of4vmOoQrtu7chtn3Nuf7g0JwgigZiOpqhqEaft8q+qs90cRgCLqZ8x4MzyOJuEt1PVRe6G3uHEa+AkXnvVff0OcL3E0HazmbwpbvBUp5mT/8Kl9c6vmLyJ0oe2c3zfMYa8ew29Rl7NRf/WFnD6Xga/3YtgltN8qjtZy+CpPVnxnc3knsiNXAOQRRY1x2o4sOkAGdn21yXRmjFC1ElElkW9n6KqU5rxre4H3op6f4WIfA4cAf5VVefTjMRrIhJOvLbvbN/U/sakiOX3beLav1wJwIrJmzlx8Hi9af3zb3GaNm06taXgxcvomX8Z5QvWE84ssWPtLrJ+3Z2Oj+XU63vZ+Iv9XPf/egBQdfAEwaxgJABFmzdmpS0j8ECMzR+IL/EaACLyLzjZA15zT+0ELlXV/SJyLfC+iPTFg8Rr1kRKYovvXMPiO91+D1WW37eJz797KrgsnLAqMuW/e0F3hr3fh4IXL2P+uDLOezSIqjL0vd7MG7OSgX+7ilULnIWN4eHoYe/3YXdFJDcXwz7ow+C3exGqdmrVS+5ay/YnDrP9icOMmpXPottXn8OfvjUQz+fBiMhkYDzwzXAeI1WtdnPIo6rLgY1AT2JLvIYlXksTQ965JvKP/dq/XkXVoRO07XxevZpE4ScDKL61giuf6si6/7WX/b+qIiM7g83/coAOVzizdwfd3J/iWysY8WFfFtxWQSDD+WNffNdaAm4zaf64Mkq+t5U1j+6OPD9UHeKLip1sX/lFpF8nPIxtWk4gYX0wDT5fZCzwc+A2VT0edb6ziLNEW0SuxOnM3eRF4jULMEluxIenUtQMn96XE4dOsODWCorGlhLICrLm0d0Egk7ttedvO7Nr5wGGTutNzckQh7YeBGDpjDIysjMoGlvK8Ol9CQSdP/Yhb/fisv7dWTChgpEz+lPw4mUc2ubcs/jOtQyf3pe6WmXY+33IzM1k6LTeHNt7tMFZwyY+iRpFchOvLQJ6icgOEXkAeBbIA+acNhxdCJSKyEqcDtsfqGq4NvIw8CKwAadmE5147QI38dqjwGOx/HwSQxDyXJ6014GM8LsYSa3k/i1U7qpk5AxnvdDK723l8I7D1IXqyGybyfAPovtjyhn5sROYFt+xhpz2OdQcraHqcFVkkePiO9cy5J1erPz+No58eZiRH59ah7T4zrWEqk/S/tIO1ByrobqymtqaUL1gN39cWaQs5pS5fLw81r6Sgq8W6Kfz58T03E7nXRjzc5OJdfKmiMpdleRdlMe6n+/lxP7jHNp2iFGz8wGYe0NJvWtrq0OsmLyJr756JaGaEIFggEFvOTNy5964ktFz8hnyTi9W/Wgn+X+4lPm3lDl9OQKZuZkIAUZ82I/iW8vrBZXqP0D2j5zXGbmZLLpjDZm5GQz829Xn5peQhlJ5lm4srAaTQpa4+7YMfqtnpAYx98YSRs8piFxT9vAOamtqObT1IBm5GYz4sB9rHt3FoW2HCFWFIssGGjN/XBkZORlIIEDN0Wouyu9Kz990rnfNggmrGP5Bn8j708vQmjWnBvOVrxboP4v/EdNzO7S9ICVrMNYHkwLC21oGAkJWmyzAWSs076aVjJ5TwLybVrJo4hqWfWsj1ZXVFLx0GR0u74AEnD/ea35/EaGqEBk5GZHazqKJZ44IzR9XRr9RPcntkMuQd3pR+MkAev6mc2TUadHE1ZQ9vKNecAFiGKw0DfN+FMlvqVvyVuTo7qOn3ggU3ewMTYdrn6Nm59OuWx4nj9cw8G9X8fn9WwhVh+jUsxML3OUBIz7qRzDLaREf+d1JaqtrKb6lPDLM/fl3tzByRn/Wr9hKMDsjsgL70G9r0Drn++R2bBPpOJ43ZiVb/+0QSyfZauuWaPWLHRtZRNVRROaIyHr3a4eozx53F0StFZExXhW8NRk1O58V39lMXfgfuju1f/ScgshixENbDzH0vd4AnDxWQ+XOSiq/PMLwD/pQdHMZn31jPdWV1eR1zWP9ym1kts0ks20mGTlO0Dm29yhzbyzh+IHjHNi4P7KsYPXiTZFmVfXhqkjHbna7HLYs2cqJg1XkdWtHua22bjbbrsHxCmcuonoM+FRVewCfuu8RkT7AJKCve89z4fF20zI1ldUMfqsn2XnZXPdGj0izJVQVoujmMtp1bwfA0nvXR+7J6ZBL+SNfUDizP6HqWgJB4dq/XoXWKdWV1VRXVjPs/T6RIDV6TgGj5xSQ2z6X3Pa5dHvvQoa+dw0LJ6yK1GhKH9zGoomrGfJ2L6eTWZWju4/S77mzLnMxjZAY/0tVTY4iqWpRA0uzJwCj3devAnNxJvRMAN5U1WpgsztmPghnfN60QHiTqD7/2ZUFt1UwfHpfyh7eQeHz9XerG/RGDz77xgYy22bR71lnGcln39zA0NdObTI18DVn1Kd4fDlzbyghWB1kZNR2mVWHq8hsk8mXE/cw94YSRn9QQNHY0jOWDyyauIZR7+XbKut4CXFPoksV8f50XdxZf7hfL3TPRxZEuaIXS5kWKr6lnIr//mVkpm3/50/VGsoe3sHSSU7tpeZYDcGMACu/vw2A6147cxh52bc2EqoKMfrvBfVWWc+/pYzCTwYw9L3eLJq4mtF/d0aHCj8ZwPL7NtV7Rs2xauaNWcnQ03bIM7FpDU2kRM+DiXlBlIg8iLPvBNnkJrgY6SmQFeTwjsOUPFBDqCrEwNeuZvEdaxjy7jX1gk0wI0AgM0jlrkqKb60gEBQycjMZ9EaPyDVHd1VGgkcgM8DC21eTmZtJ3b46isaW0v6y9rSpa1Pv+1dXVrNo4upIX8+oWfnn4KdOX6pKMkwT8VK8AWa3iHRV1Z3uHhJ73PORBVGu6MVS9bhLzaeAMw8mznK0KuE1RMf2HoucG/LuNSydtJ5Bb/aIrBeSo8KoWfnsfeo4FXPXNThHJRxciseX06X3RfT6D6cSuvDrqxj2fh9KvreVghcvY+Htq6l110O1yWsbWdH9+f1boE6pPlpD1aETkeeZZqrzuwDeijfAhBc+/dr9Gr0g6nUR+T3QDWcR1dKWFtI4hrzdK9L/AlD6g+0c3VUJVc7nGTkZ1NbUEswMUnxLOefvaM/oOQWRiXALJ6xi2Ad9IoFo9JwC2nU/n31r9rLv1r10z+/GsH935riEg8uwab3rlWHNT3dzeMdhqFQGv92LuTeURDasMs3X6msw7iKq0Tgb3uwAnsAJLFPdBVXbgLsAVLVCRKYCq3D2n/ihqtY2+GDTbM6cFqF4fDkiQh7tyMjJZNCbPSh5YCsZhzOora5lxEf9KB5fzuFth4CLye2Qy/L7NpF5MpP5t5QhJ4VRs/MJ/DmLE/84wUfP/ZGvP/o/uPzfO0SyCAye2pNgpjMAWDSmlC79u9Drd1245ukuQBfA2c5h9NtWc4mbEpljlK5iGUW6t5GPrm/k+qeAp1pSKNOw7POyIiNAAEsnrSczN4NdvzzGic+PM/S93iy+cw0bHt/PiF+dWhIweGovim8td7I8fuzMAC4eX067FecjwK8u+RN6V/hqZfDUXiy5e10ko0DhrAEw68zynD6qZOKQ5k2k1O2eboXCwWXh7atZNHE1g97sQdsLz+PY8SrOv6Q9ADdOuI6rf3UBgT87SwrCWyuM+LAfIz/u73TgXu7Mi6w6dCLyurYqxOI71zB4ai+WfWsjg6f2ZN5NKyPfe8Xk+iNIJjHCHb1NHanKAkwKGjatN3W1zl+6fWv3kpOd6cyPubWCj19fAEDR9BWUPLCFnr/pzMLbV0eGmNtf3oGaozW0634+Oe1z0do6Sh7YyvAP+zLknWtYfOcaBv7tKoDIau3l9zkrs+FUVkiTAG4TKZYjVdl2DSkoPPkNnG0veR+KxpZS+KEz6a54fDl5XdtFRp1CVSHCy1kG/MkZ5Atvz9CuWzsKXroccBYzDn2nfqfuoomr6dC5I/NuWsmo2flnZIU08VMghSsnMbEaTIqZe0PJGUPCSyetI7NNZr0V0vkvXEpdyOlfz8zJQN0aT9GYUj67dz2FnwygcOYAjnxxJLLCOjy/BaD0oe0UjSklVF3LNU93idRmTCIp1MV4NCFRawYt8VorFw4u0cFk0Js9Gfpe70iACC9OzH/hMubeWEJ1ZTUSEEp/sJ3CWQO47o0ekaZO4ScDaHNBG5beU39V9KEtB8k+P4c2F9SfbGcSSBPaB/MKiVkzaInXDNQcOxnJOLD2Z3to++p5kfktS+9Zx7wbV9bbCGrotN6R5hFQr6lTV1sX2fEOnCZW4awBDJ7ak1BVCOMdrYvtaPI5DSReo36ytFepn0TtTTe7wGac/XcHeZF4zQJMisrOy2LIO06QOLTtIMcmH41MqsvKy2bUnHxyO7Zh6b3rCWYHyX7ZTah296mayvJvb2T+LWVktc2i7GFnu4XiW8ojNaDi8eUMntqTVT/eyZK71kaSsZkEir2J1ElElkUdD8bw9OauGYw58RoQTrx2VtbJm6LCwQVOZXfc/K8HkaXCkR2HKRpbSocLOtL/uYtZcFsF1fefcDpqp+az+ic72V2+m+CBDDpeeQF9/9iNoptLKR5/iA6XR5rpkUDT5z+7Al3P6c/XKmizZvLGnXitAY2tGbTEa6ZhZQ/voHLvUUbNziczN5PsdtmRz4ZP7xsZBQLYu2YvgYwAtdUhDmzcz/ybyyicOYARH/VzlgGYc0JJXBOpEbvdZk8473RTawYt8Zpp2IGN+zmwaT/g9LecOHiC/s9dzEJ3y8xRs/OZf0sZcKpjV9x8SiNn9mf+uLLIQse1/3NPw9/EJFhsHbwtmGgXnSzt9DWDk9yRoStw1wx6kXjNmkhpInoYefEda8ghx+nk/cBZ5BjMDpJFVuSaUHWI7LwcqiurzshxFF5ZbTymxDQEHYsErhl8GGdEKhcn6Vp04rW/upvIHcAZhWq6XMkwDdnSlvhj7g0l5HbMjfThmJZrVuK1/gX69w8+jem5na/qlJJpS6wG04rZHi7+S+VlALGwAGOMn5KgBeElCzDG+EVbNEKUEizAGOMTZ5jaajDGGE+k9l4vsbAAY4xfbMtMY4yXLMAYYzyhSiTfeLqyAGOMjyzAGGM8ktr77cbCAowxPlElsnl7urIAY4yPrIlkjPGGDVMbY7yiKHV16b1WwAKMMX6xYWpjjJfSvYlkW2Ya4xPFqcHEcjRFRHqJSEnUcUREfiIiT4rIF1Hnx0Xd06zka/GwGowxfklgE0lV1wIFAG4StS+AacB3gWdU9XfR15+WfK0b8HcR6elunRlOvrYYmIGTfG0mcbAajDG+iS3xfRzNqOuBjaq69SzXxJN8rdkswBjjk/BapBibSM1JvDYJeCPq/SMiUurmrw4nvoon+VqzWYAxxkd1tRrTgZt4LeqY0tDzRCQLuA142z31PHAVTvNpJ/B0+NIGbm8q+VqzWR+MMX7xZqLdzcAKVd0NEP4KICIvAB+5b+NJvtZsVoMxxifhiXaxHM1wL1HNo3BmR9ftQLn7Op7ka81mNRhjfJTIiXYi0ga4EXgo6vRvRaQAp5mzJfxZnMnXmq3JACMiLwPjgT2q2s899yTwfWCve9kvVHWG+9njwANALfAjVZ0Vb+GMSWsJbiKp6nHggtPOffss1z8FPNXA+WVAv0SUKZYazCvAszjDVdGaO7ZujIkSnmiXzprsg1HVIpxctLFocGy9BeUzJn01b5g6JbWkk7c5Y+tnEJEHw2P6J6lpQTGMSVWeTbRLGvEGmOaOrZ95UnVKeEw/k6w4i2FM6nJ2tKuL6UhVcY0ixTG2boxpQCo3f2IRVw2muWPrLSuiMWlK1WowIvIGMBpnLcQO4AlgdBxj68aYKM4oUuoGj1g0GWBU9d4GTr90lusbHFs3xpzG9uQ1xnhHU7r5EwsLMMb4xNmuwQKMMcYj6T6KZAHGGJ+oWhPJGOMhq8EYY7zhzuRNZxZgjPGJoqh18hpjvJLuTSTbMtMYvyR4saOIbHETppWIyDL3XEcRmSMi692vHaKu9zzxmgUYY3ySyMyOUb6mqgWqOtB9/xjwqar2AD5135++OdxY4Dk3YRucSrzWwz3GxvszWoAxxi/qyabfp5sAvOq+fpVTSdTOSeI164MxxidKs0aROoWbPa4pDeRGUmC2iCjwZ/fzLm6mAFR1p4hc6F7bHSc1bFh4c7iTJDDxmgUYY/zSvMWO+6KaPY0ZrqpfukFkjoisOcu1lnjNmPSW2Jm8qvql+3WPiEzD2Q97t4h0dWsvXYE97uWWeM2YdBZe7JiIPhgRaSsieeHXwE04G8FNBya7l03mVBI1S7xmTLpL4DyYLsA0d0Q5A3hdVT8Rkc+AqSLyALANuAuSKPGaMcYbiVzsqKqbgPwGzu8Hrm/knqRIvGaM8YjaWiRjjCcU6mrTe9Y3lCcAAANwSURBVKmABRhj/KJqNRhjjDeaOdEuJVmAMcYvtqOdMcZLan0wxhgvqO1oZ4zxjCoasgBjjPGI1WCMMZ5QG6Y2xnjGJtoZY7xkTSRjjCesiWSM8Y4NUxtjvJTuw9S2o50xPnH2g4ntaIqIXCIi/xSR1SJSISI/ds8/KSJfuLmSSkRkXNQ9nudFshqMMX7RhO4HEwJ+qqor3K0zl4vIHPezZ1T1d9EXn5YXqRvwdxHp6e5qF86LtBiYgZMXKa5d7SzAGOOjBO5otxMIpyepFJHVnD3dSCQvErBZRMJ5kbbg5kUCEJFwXqS4Aow1kYzxSXjLzESljg0TkcuBrwBL3FOPiEipiLwclTq2O7A96rZw/qPuJDAvkgUYY3yktXUxHbiJ16KOBxt6noicB7wL/ERVj+A0d64CCnBqOE+HL22oOGc5HxdrIhnjl+bN5G0y8ZqIZOIEl9dU9T0AVd0d9fkLwEfu23OSFykpAsxRDu+by8fHgH1+l6UZOpFa5QUr87lwWcxXJnA1tTvS8xKwWlV/H3W+azh1LHA7Tq4kcPIivS4iv8fp5A3nRaoVkUoRGYLTxLoP+GO85UqKAKOqnUVkWQypMZNGqpUXrMzJRqGlie2jDQe+DZSJSIl77hfAvSJS4H67LcBDYHmRjEl/qtTW1TZ9XUyP0mIa7j+ZcZZ7LC+SMekqwTWYpJRMAWaK3wVoplQrL1iZk0ziajDJKmkCjKqm1F+kVCsvWJmTjarVYIwxnlHq1GowxhgPKFBrNRhjjCdUqbM+GGOMF6wGY4zxjtVgjDFesXkwxhgP2TwYY4xHbB6MMcZDSq3NgzHGeMH6YIwx3kngaupkZQHGGJ9YDcYY451WMA/GNv02xifhmbyxHLEQkbFuErUNIvKYt6WPjdVgjPFN4mowIhIE/gu4EWfj7s9EZLqqrkrIN4iTBRhjfKKJ7eQdBGxQ1U0AIvImTnI1CzDGtFYJ7ORtKJHa4EQ9PF4WYIzxSSWHZ/1DP+wU4+U5IrIs6v2U03b7S2jCtESxAGOMT1R1bAIf11giNV/ZKJIx6eEzoIeIXCEiWcAknORqvrIajDFpQFVDIvIIMAsIAi+raoXPxUJUfW+mGWPSlDWRjDGesQBjjPGMBRhjjGcswBhjPGMBxhjjGQswxhjPWIAxxnjGAowxxjP/H8qBNii2QTqUAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 288x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot the tumour as a 2D map color coding the mutation ID\n",
    "cmap = pyplot.cm.get_cmap('PRGn', mtx.max()+1)\n",
    "\n",
    "pyplot.matshow(mtx, cmap=cmap)\n",
    "pyplot.colorbar()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}