covid-venn.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "%precision 3\n",
    "\n",
    "from pandas.plotting import register_matplotlib_converters\n",
    "register_matplotlib_converters()\n",
    "\n",
    "import matplotlib\n",
    "import math\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "plt.rcParams['figure.figsize'] = [9.5, 7]\n",
    "import datetime\n",
    "import pandas\n",
    "from matplotlib_venn import venn3\n",
    "from matplotlib_venn import venn2\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "prefix='/home/ahu/git/antonie'\n",
    "data=pandas.read_csv(prefix+\"/combined-codons-s-protein.csv\")\n",
    "data.head()\n",
    "data=data.reset_index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "def amalg(row, nam):\n",
    "    return str(row[\"index\"])+\" \"+row[nam]\n",
    "\n",
    "data[\"indorig\"] = data.apply(lambda row: amalg(row, \"orig\"), axis=1)\n",
    "data[\"indcvac\"] = data.apply(lambda row: amalg(row, \"cvac\"), axis=1)        \n",
    "data[\"indbntc\"] = data.apply(lambda row: amalg(row, \"bntc\"), axis=1)        \n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>index</th>\n",
       "      <th>orig</th>\n",
       "      <th>cvac</th>\n",
       "      <th>bntc</th>\n",
       "      <th>indorig</th>\n",
       "      <th>indcvac</th>\n",
       "      <th>indbntc</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>AUG</td>\n",
       "      <td>AUG</td>\n",
       "      <td>AUG</td>\n",
       "      <td>0 AUG</td>\n",
       "      <td>0 AUG</td>\n",
       "      <td>0 AUG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>UUU</td>\n",
       "      <td>UUC</td>\n",
       "      <td>UUC</td>\n",
       "      <td>1 UUU</td>\n",
       "      <td>1 UUC</td>\n",
       "      <td>1 UUC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>GUU</td>\n",
       "      <td>GUG</td>\n",
       "      <td>GUG</td>\n",
       "      <td>2 GUU</td>\n",
       "      <td>2 GUG</td>\n",
       "      <td>2 GUG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>UUU</td>\n",
       "      <td>UUC</td>\n",
       "      <td>UUC</td>\n",
       "      <td>3 UUU</td>\n",
       "      <td>3 UUC</td>\n",
       "      <td>3 UUC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>CUU</td>\n",
       "      <td>CUG</td>\n",
       "      <td>CUG</td>\n",
       "      <td>4 CUU</td>\n",
       "      <td>4 CUG</td>\n",
       "      <td>4 CUG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1269</th>\n",
       "      <td>1269</td>\n",
       "      <td>UUA</td>\n",
       "      <td>CUC</td>\n",
       "      <td>CUG</td>\n",
       "      <td>1269 UUA</td>\n",
       "      <td>1269 CUC</td>\n",
       "      <td>1269 CUG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1270</th>\n",
       "      <td>1270</td>\n",
       "      <td>CAU</td>\n",
       "      <td>CAC</td>\n",
       "      <td>CAC</td>\n",
       "      <td>1270 CAU</td>\n",
       "      <td>1270 CAC</td>\n",
       "      <td>1270 CAC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1271</th>\n",
       "      <td>1271</td>\n",
       "      <td>UAC</td>\n",
       "      <td>UAC</td>\n",
       "      <td>UAC</td>\n",
       "      <td>1271 UAC</td>\n",
       "      <td>1271 UAC</td>\n",
       "      <td>1271 UAC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1272</th>\n",
       "      <td>1272</td>\n",
       "      <td>ACA</td>\n",
       "      <td>ACG</td>\n",
       "      <td>ACA</td>\n",
       "      <td>1272 ACA</td>\n",
       "      <td>1272 ACG</td>\n",
       "      <td>1272 ACA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1273</th>\n",
       "      <td>1273</td>\n",
       "      <td>UAA</td>\n",
       "      <td>UGA</td>\n",
       "      <td>UGA</td>\n",
       "      <td>1273 UAA</td>\n",
       "      <td>1273 UGA</td>\n",
       "      <td>1273 UGA</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1274 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      index orig cvac bntc   indorig   indcvac   indbntc\n",
       "0         0  AUG  AUG  AUG     0 AUG     0 AUG     0 AUG\n",
       "1         1  UUU  UUC  UUC     1 UUU     1 UUC     1 UUC\n",
       "2         2  GUU  GUG  GUG     2 GUU     2 GUG     2 GUG\n",
       "3         3  UUU  UUC  UUC     3 UUU     3 UUC     3 UUC\n",
       "4         4  CUU  CUG  CUG     4 CUU     4 CUG     4 CUG\n",
       "...     ...  ...  ...  ...       ...       ...       ...\n",
       "1269   1269  UUA  CUC  CUG  1269 UUA  1269 CUC  1269 CUG\n",
       "1270   1270  CAU  CAC  CAC  1270 CAU  1270 CAC  1270 CAC\n",
       "1271   1271  UAC  UAC  UAC  1271 UAC  1271 UAC  1271 UAC\n",
       "1272   1272  ACA  ACG  ACA  1272 ACA  1272 ACG  1272 ACA\n",
       "1273   1273  UAA  UGA  UGA  1273 UAA  1273 UGA  1273 UGA\n",
       "\n",
       "[1274 rows x 7 columns]"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'1134 AAC',\n",
       " '431 UGC',\n",
       " '107 ACU',\n",
       " '325 AUU',\n",
       " '1231 AUU',\n",
       " '383 CCU',\n",
       " '660 GAG',\n",
       " '807 GAU',\n",
       " '188 CUU',\n",
       " '625 GCA',\n",
       " '652 GCU',\n",
       " '1250 GGA',\n",
       " '3 UUU',\n",
       " '115 UCC',\n",
       " '1062 UUG',\n",
       " '103 UGG',\n",
       " '1223 UUG',\n",
       " '211 UUA',\n",
       " '342 AAC',\n",
       " '1220 AUA',\n",
       " '363 GAU',\n",
       " '232 AUU',\n",
       " '2 GUU',\n",
       " '533 GUU',\n",
       " '313 CAA',\n",
       " '872 UAC',\n",
       " '509 GUA',\n",
       " '494 UAU',\n",
       " '121 AAC',\n",
       " '1156 AAG',\n",
       " '724 GAA',\n",
       " '93 UCC',\n",
       " '388 GAU',\n",
       " '689 CAA',\n",
       " '1005 ACA',\n",
       " '1018 AGA',\n",
       " '826 ACA',\n",
       " '1025 GCU',\n",
       " '977 AAU',\n",
       " '840 CUU',\n",
       " '1006 UAU',\n",
       " '497 CAA',\n",
       " '346 UUU',\n",
       " '6 UUA',\n",
       " '843 AUU',\n",
       " '935 GAC',\n",
       " '371 GCA',\n",
       " '766 UUA',\n",
       " '98 AAC',\n",
       " '965 CUU',\n",
       " '133 CAA',\n",
       " '246 AGU',\n",
       " '476 AGC',\n",
       " '505 CAA',\n",
       " '14 UGU',\n",
       " '490 CCU',\n",
       " '1153 AAA',\n",
       " '570 GAC',\n",
       " '171 UCU',\n",
       " '348 UCU',\n",
       " '13 CAG',\n",
       " '992 AUU',\n",
       " '328 UUU',\n",
       " '140 UUG',\n",
       " '1126 GAU',\n",
       " '925 CAA',\n",
       " '536 AAA',\n",
       " '1164 GAU',\n",
       " '308 GAA',\n",
       " '864 CUC',\n",
       " '945 GGA',\n",
       " '17 CUU',\n",
       " '937 CUU',\n",
       " '1016 GAA',\n",
       " '538 GUC',\n",
       " '113 ACC',\n",
       " '915 CUC',\n",
       " '1255 UUU',\n",
       " '77 AGG',\n",
       " '511 GUA',\n",
       " '95 GAG',\n",
       " '445 GGU',\n",
       " '1248 UCU',\n",
       " '355 AAG',\n",
       " '474 GCC',\n",
       " '1252 UGC',\n",
       " '748 UGC',\n",
       " '713 AUA',\n",
       " '831 GGC',\n",
       " '1004 CAG',\n",
       " '531 AAU',\n",
       " '731 ACC',\n",
       " '296 UCA',\n",
       " '921 UUG',\n",
       " '330 AAU',\n",
       " '716 AAU',\n",
       " '927 AAU',\n",
       " '290 UGU',\n",
       " '722 ACC',\n",
       " '183 GGU',\n",
       " '508 AGA',\n",
       " '1246 UGU',\n",
       " '670 UGC',\n",
       " '407 AGA',\n",
       " '958 UUA',\n",
       " '866 GAU',\n",
       " '34 GGU',\n",
       " '100 AUA',\n",
       " '550 GUU',\n",
       " '762 UUA',\n",
       " '830 GCU',\n",
       " '847 GAC',\n",
       " '467 AUU',\n",
       " '791 CCA',\n",
       " '219 UUU',\n",
       " '1030 GAG',\n",
       " '828 GCA',\n",
       " '909 GGA',\n",
       " '70 UCU',\n",
       " '1225 GCC',\n",
       " '257 UGG',\n",
       " '778 CAA',\n",
       " '215 CUC',\n",
       " '626 GAU',\n",
       " '229 CCA',\n",
       " '263 GCU',\n",
       " '853 AAG',\n",
       " '390 UGC',\n",
       " '478 CCU',\n",
       " '801 UUU',\n",
       " '1138 GAU',\n",
       " '1261 GAG',\n",
       " '936 UCA',\n",
       " '1084 GGA',\n",
       " '545 UUA',\n",
       " '910 GUU',\n",
       " '405 GAA',\n",
       " '594 GUC',\n",
       " '964 CAA',\n",
       " '67 AUA',\n",
       " '1201 GAA',\n",
       " '1055 GCA',\n",
       " '74 GGU',\n",
       " '1159 ACA',\n",
       " '1124 AAC',\n",
       " '860 UUG',\n",
       " '345 AGA',\n",
       " '1263 GUG',\n",
       " '880 ACA',\n",
       " '777 ACC',\n",
       " '911 ACA',\n",
       " '486 AAU',\n",
       " '1061 UUC',\n",
       " '316 AAC',\n",
       " '581 CUU',\n",
       " '253 UCU',\n",
       " '1234 UGC',\n",
       " '1208 UAU',\n",
       " '408 CAA',\n",
       " '666 GGU',\n",
       " '137 GAU',\n",
       " '822 UUC',\n",
       " '193 UUU',\n",
       " '1051 UUC',\n",
       " '601 ACA',\n",
       " '286 GAU',\n",
       " '1149 GAG',\n",
       " '1010 CAA',\n",
       " '920 AAA',\n",
       " '380 GGA',\n",
       " '695 ACU',\n",
       " '1191 AAU',\n",
       " '827 CUU',\n",
       " '1039 GUU',\n",
       " '282 GGA',\n",
       " '426 GAU',\n",
       " '871 CAA',\n",
       " '1116 ACA',\n",
       " '158 GUU',\n",
       " '635 UAU',\n",
       " '503 GGU',\n",
       " '341 UUU',\n",
       " '589 UGU',\n",
       " '1152 GAU',\n",
       " '1260 UCU',\n",
       " '1117 GAC',\n",
       " '1031 UGU',\n",
       " '146 AAA',\n",
       " '1221 GCU',\n",
       " '294 CCU',\n",
       " '952 AAC',\n",
       " '723 ACA',\n",
       " '354 AGG',\n",
       " '795 GAU',\n",
       " '142 GUU',\n",
       " '185 UUC',\n",
       " '404 GAU',\n",
       " '998 GGC',\n",
       " '1119 ACA',\n",
       " '1163 GUU',\n",
       " '440 CUU',\n",
       " '1141 CAA',\n",
       " '903 UAU',\n",
       " '815 UCA',\n",
       " '197 GAU',\n",
       " '329 CCU',\n",
       " '747 GAA',\n",
       " '848 CUC',\n",
       " '1046 UAU',\n",
       " '502 GUU',\n",
       " '718 ACU',\n",
       " '218 GGU',\n",
       " '559 CUG',\n",
       " '347 GCA',\n",
       " '699 GGU',\n",
       " '1072 AAG',\n",
       " '37 UAC',\n",
       " '288 GUA',\n",
       " '606 CAG',\n",
       " '214 GAU',\n",
       " '391 UUU',\n",
       " '844 GCU',\n",
       " '134 UUU',\n",
       " '1093 GUC',\n",
       " '305 UUC',\n",
       " '62 ACU',\n",
       " '72 ACC',\n",
       " '969 UUU',\n",
       " '892 GCA',\n",
       " '654 CAU',\n",
       " '534 AAA',\n",
       " '779 GAA',\n",
       " '261 GCU',\n",
       " '427 GAU',\n",
       " '315 UCU',\n",
       " '616 UGC',\n",
       " '744 GAU',\n",
       " '571 ACU',\n",
       " '1076 ACU',\n",
       " '504 UAC',\n",
       " '887 UUU',\n",
       " '1054 UCA',\n",
       " '39 GAC',\n",
       " '543 AAU',\n",
       " '259 GCU',\n",
       " '499 ACU',\n",
       " '1214 UAC',\n",
       " '628 CUU',\n",
       " '157 AGA',\n",
       " '399 UUU',\n",
       " '1135 ACA',\n",
       " '974 AGU',\n",
       " '378 UGU',\n",
       " '1019 GCU',\n",
       " '849 AUU',\n",
       " '818 GAA',\n",
       " '857 CUU',\n",
       " '1176 GUA',\n",
       " '258 ACA',\n",
       " '738 ACA',\n",
       " '236 AGG',\n",
       " '468 UCA',\n",
       " '275 CUA',\n",
       " '401 AUU',\n",
       " '636 UCU',\n",
       " '1129 AUA',\n",
       " '639 UCU',\n",
       " '1179 CAA',\n",
       " '841 GGU',\n",
       " '1036 UCA',\n",
       " '1015 GCA',\n",
       " '33 CGU',\n",
       " '271 CCU',\n",
       " '1177 AAC',\n",
       " '755 UAU',\n",
       " '1057 CAU',\n",
       " '36 UAU',\n",
       " '58 UUU',\n",
       " '882 ACU',\n",
       " '640 AAU',\n",
       " '335 UGC',\n",
       " '552 ACU',\n",
       " '979 AUC',\n",
       " '16 AAU',\n",
       " '4 CUU',\n",
       " '618 GAA',\n",
       " '117 CUU',\n",
       " '240 UUA',\n",
       " '765 GCU',\n",
       " '242 GCU',\n",
       " '333 AAC',\n",
       " '451 CUG',\n",
       " '1079 GCC',\n",
       " '206 CAC',\n",
       " '1168 AUC',\n",
       " '302 UUG',\n",
       " '1125 UGU',\n",
       " '1206 GAG',\n",
       " '1151 UUA',\n",
       " '1273 UAA',\n",
       " '1012 AUU',\n",
       " '1236 AUG',\n",
       " '251 GGU',\n",
       " '76 AAG',\n",
       " '752 UUG',\n",
       " '753 UUG',\n",
       " '957 GCU',\n",
       " '637 ACA',\n",
       " '528 AAG',\n",
       " '829 GAU',\n",
       " '577 GAU',\n",
       " '297 GAA',\n",
       " '500 AAU',\n",
       " '548 ACA',\n",
       " '498 CCC',\n",
       " '162 GCG',\n",
       " '168 GAA',\n",
       " '367 CUA',\n",
       " '1091 GAA',\n",
       " '466 GAU',\n",
       " '1103 GUA',\n",
       " '379 UAU',\n",
       " '521 GCA',\n",
       " '816 UUU',\n",
       " '1139 CCU',\n",
       " '262 GCA',\n",
       " '303 AAA',\n",
       " '513 UCU',\n",
       " '825 GUG',\n",
       " '274 UUU',\n",
       " '1228 AUG',\n",
       " '1245 GGC',\n",
       " '441 GAU',\n",
       " '272 AGG',\n",
       " '374 UCC',\n",
       " '737 UGU',\n",
       " '1184 CGC',\n",
       " '47 UUA',\n",
       " '558 UUU',\n",
       " '685 AGU',\n",
       " '526 CCU',\n",
       " '873 ACU',\n",
       " '483 GAA',\n",
       " '141 GGU',\n",
       " '84 CCA',\n",
       " '858 ACU',\n",
       " '1262 CCA',\n",
       " '529 UCU',\n",
       " '671 GCU',\n",
       " '403 GGU',\n",
       " '812 AGC',\n",
       " '852 CAA',\n",
       " '248 UUG',\n",
       " '484 GGU',\n",
       " '805 UUA',\n",
       " '1217 CUA',\n",
       " '1209 AUA',\n",
       " '94 ACU',\n",
       " '97 UCU',\n",
       " '309 AAA',\n",
       " '1241 AGU',\n",
       " '926 UUU',\n",
       " '464 GAG',\n",
       " '999 AGA',\n",
       " '754 CAA',\n",
       " '1142 CCU',\n",
       " '567 GAC',\n",
       " '537 UGU',\n",
       " '1048 CUU',\n",
       " '951 GUC',\n",
       " '54 UUC',\n",
       " '1131 AUU',\n",
       " '166 ACU',\n",
       " '471 AUC',\n",
       " '470 GAA',\n",
       " '900 CAA',\n",
       " '340 GUU',\n",
       " '1017 AUC',\n",
       " '733 ACA',\n",
       " '694 UAC',\n",
       " '353 AAC',\n",
       " '838 GAU',\n",
       " '938 UCU',\n",
       " '145 CAC',\n",
       " '879 GGU',\n",
       " '78 UUU',\n",
       " '1268 AAA',\n",
       " '57 UUC',\n",
       " '312 UAU',\n",
       " '576 CGU',\n",
       " '1066 UAU',\n",
       " '147 AAC',\n",
       " '894 CAA',\n",
       " '270 CAA',\n",
       " '1097 AAU',\n",
       " '99 AUA',\n",
       " '1056 CCU',\n",
       " '712 GCC',\n",
       " '686 GUA',\n",
       " '135 UGU',\n",
       " '1013 AGA',\n",
       " '821 CUU',\n",
       " '15 GUU',\n",
       " '1205 UAU',\n",
       " '811 CCA',\n",
       " '1022 AAU',\n",
       " '55 UUA',\n",
       " '845 GCU',\n",
       " '327 AGA',\n",
       " '641 GUU',\n",
       " '306 ACU',\n",
       " '1146 UCA',\n",
       " '781 UUU',\n",
       " '1162 GAU',\n",
       " '155 GAG',\n",
       " '60 AAU',\n",
       " '417 AUU',\n",
       " '1096 UCA',\n",
       " '1267 GUC',\n",
       " '1107 AAU',\n",
       " '1002 AGU',\n",
       " '622 GCU',\n",
       " '373 UUU',\n",
       " '91 UUU',\n",
       " '48 CAU',\n",
       " '291 GCA',\n",
       " '1239 UGC',\n",
       " '177 GAC',\n",
       " '1264 CUC',\n",
       " '1130 GGA',\n",
       " '1186 AAU',\n",
       " '1037 AAA',\n",
       " '322 ACA',\n",
       " '163 AAU',\n",
       " '350 UAU',\n",
       " '908 AUU',\n",
       " '1258 GAC',\n",
       " '984 GAC',\n",
       " '456 AGG',\n",
       " '750 AAU',\n",
       " '987 GAG',\n",
       " '784 GUC',\n",
       " '1021 GCU',\n",
       " '1069 GCA',\n",
       " '120 AAU',\n",
       " '745 UCA',\n",
       " '439 AAU',\n",
       " '607 GUU',\n",
       " '25 CCU',\n",
       " '799 UUU',\n",
       " '800 AAU',\n",
       " '359 AAC',\n",
       " '289 GAC',\n",
       " '482 GUU',\n",
       " '610 CUU',\n",
       " '243 UUA',\n",
       " '23 UUA',\n",
       " '1227 GUA',\n",
       " '1088 UUU',\n",
       " '409 AUC',\n",
       " '891 GCU',\n",
       " '863 UUG',\n",
       " '855 AAC',\n",
       " '681 CGG',\n",
       " '1009 CAA',\n",
       " '150 AGU',\n",
       " '1192 UUA',\n",
       " '835 CAA',\n",
       " '813 AAG',\n",
       " '532 UUG',\n",
       " '170 GUC',\n",
       " '975 GUU',\n",
       " '1183 GAC',\n",
       " '592 GGU',\n",
       " '375 ACU',\n",
       " '976 UUA',\n",
       " '143 UAU',\n",
       " '457 AAG',\n",
       " '488 UAC',\n",
       " '293 GAC',\n",
       " '479 UGU',\n",
       " '1251 UCC',\n",
       " '181 AAA',\n",
       " '808 CCA',\n",
       " '1032 GUA',\n",
       " '890 GGU',\n",
       " '874 UCU',\n",
       " '1155 UUU',\n",
       " '611 UAU',\n",
       " '991 CAA',\n",
       " '541 AAC',\n",
       " '972 AUU',\n",
       " '475 GGU',\n",
       " '264 UAU',\n",
       " '110 GAU',\n",
       " '127 AUU',\n",
       " '352 UGG',\n",
       " '1029 UCA',\n",
       " '889 GCA',\n",
       " '959 AAC',\n",
       " '988 GCU',\n",
       " '1104 ACA',\n",
       " '361 GUU',\n",
       " '1207 CAG',\n",
       " '593 GGU',\n",
       " '71 GGG',\n",
       " '237 UUU',\n",
       " '292 CUU',\n",
       " '644 ACA',\n",
       " '454 UUG',\n",
       " '525 GGA',\n",
       " '881 AUC',\n",
       " '1080 AUU',\n",
       " '824 AAA',\n",
       " '919 CAA',\n",
       " '1218 GGU',\n",
       " '64 UUC',\n",
       " '204 UCU',\n",
       " '763 AAC',\n",
       " '929 GCU',\n",
       " '9 CUA',\n",
       " '823 AAC',\n",
       " '485 UUU',\n",
       " '46 GUU',\n",
       " '30 UCU',\n",
       " '397 GAU',\n",
       " '599 CCA',\n",
       " '487 UGU',\n",
       " '933 AUU',\n",
       " '512 CUU',\n",
       " '83 CUA',\n",
       " '648 UGU',\n",
       " '1224 AUU',\n",
       " '372 UCA',\n",
       " '751 CUU',\n",
       " '591 UUU',\n",
       " '986 GUU',\n",
       " '961 CUU',\n",
       " '967 UCC',\n",
       " '65 CAU',\n",
       " '650 AUA',\n",
       " '793 AUU',\n",
       " '1020 UCU',\n",
       " '674 CAG',\n",
       " '851 GCA',\n",
       " '190 GAA',\n",
       " '798 GGU',\n",
       " '651 GGG',\n",
       " '73 AAU',\n",
       " '714 CCC',\n",
       " '480 AAU',\n",
       " '81 CCU',\n",
       " '75 ACU',\n",
       " '1058 GGU',\n",
       " '35 GUU',\n",
       " '250 CCU',\n",
       " '868 AUG',\n",
       " '1095 GUU',\n",
       " '861 CCA',\n",
       " '705 GCU',\n",
       " '398 UCA',\n",
       " '31 UUC',\n",
       " '715 ACA',\n",
       " '394 GUC',\n",
       " '656 AAC',\n",
       " '1167 GAC',\n",
       " '563 CAA',\n",
       " '772 GAA',\n",
       " '870 GCU',\n",
       " '85 UUU',\n",
       " '632 UGG',\n",
       " '249 ACU',\n",
       " '627 CAA',\n",
       " '102 GGC',\n",
       " '944 CUU',\n",
       " '943 GCA',\n",
       " '814 AGG',\n",
       " '898 GCU',\n",
       " '960 ACG',\n",
       " '970 GGU',\n",
       " '680 CCU',\n",
       " '1180 AAA',\n",
       " '1068 CCU',\n",
       " '948 CAA',\n",
       " '331 AUU',\n",
       " '463 UUU',\n",
       " '859 GUU',\n",
       " '940 ACA',\n",
       " '49 UCA',\n",
       " '389 CUC',\n",
       " '116 CUA',\n",
       " '743 GGU',\n",
       " '438 AAC',\n",
       " '179 GAA',\n",
       " '1173 GCU',\n",
       " '207 ACG',\n",
       " '239 ACU',\n",
       " '92 GCU',\n",
       " '1008 ACU',\n",
       " '69 GUC',\n",
       " '1077 GCU',\n",
       " '732 AAG',\n",
       " '410 GCU',\n",
       " '280 GAA',\n",
       " '1200 CAA',\n",
       " '668 GGU',\n",
       " '1110 GAA',\n",
       " '1064 GUG',\n",
       " '702 AAU',\n",
       " '1098 GGC',\n",
       " '797 GGU',\n",
       " '834 AAA',\n",
       " '510 GUA',\n",
       " '787 AUU',\n",
       " '918 AAC',\n",
       " '1266 GGA',\n",
       " '1003 UUG',\n",
       " '1100 CAC',\n",
       " '619 GUC',\n",
       " '697 UCA',\n",
       " '63 UGG',\n",
       " '1082 CAU',\n",
       " '803 CAA',\n",
       " '323 GAA',\n",
       " '105 UUU',\n",
       " '336 CCU',\n",
       " '783 CAA',\n",
       " '338 GGU',\n",
       " '138 CCA',\n",
       " '1011 UUA',\n",
       " '287 GCU',\n",
       " '719 AUU',\n",
       " '1007 GUG',\n",
       " '1023 CUU',\n",
       " '130 UGU',\n",
       " '245 AGA',\n",
       " '230 AUA',\n",
       " '124 AAU',\n",
       " '520 CCA',\n",
       " '1169 UCU',\n",
       " '191 UUU',\n",
       " '216 CCU',\n",
       " '573 GAU',\n",
       " '196 AUU',\n",
       " '175 CUU',\n",
       " '42 UUC',\n",
       " '703 UCA',\n",
       " '1147 UUC',\n",
       " '1171 AUU',\n",
       " '449 AAU',\n",
       " '523 GUU',\n",
       " '1235 UGU',\n",
       " '56 CCU',\n",
       " '429 ACA',\n",
       " '428 UUU',\n",
       " '128 AAA',\n",
       " '459 AAU',\n",
       " '1188 GUU',\n",
       " '518 CAU',\n",
       " '642 UUU',\n",
       " '1144 UUA',\n",
       " '953 CAA',\n",
       " '922 AUU',\n",
       " '337 UUU',\n",
       " '435 UGG',\n",
       " '598 ACA',\n",
       " '195 AAU',\n",
       " '32 ACA',\n",
       " '26 GCA',\n",
       " '804 AUA',\n",
       " '728 GUG',\n",
       " '912 CAG',\n",
       " '1196 CUC',\n",
       " '1133 AAC',\n",
       " '701 GAA',\n",
       " '414 ACU',\n",
       " '51 CAG',\n",
       " '764 CGU',\n",
       " '700 GCA',\n",
       " '709 AAC',\n",
       " '604 UCU',\n",
       " '1172 AAU',\n",
       " '165 UGC',\n",
       " '44 UCC',\n",
       " '87 GAU',\n",
       " '711 AUU',\n",
       " '29 AAU',\n",
       " '757 AGU',\n",
       " '634 GUU',\n",
       " '1140 UUG',\n",
       " '1060 GUC',\n",
       " '66 GCU',\n",
       " '119 GUU',\n",
       " '673 UAU',\n",
       " '12 AGU',\n",
       " '544 GGU',\n",
       " '672 AGU',\n",
       " '213 CGU',\n",
       " '1238 AGU',\n",
       " '1187 GAG',\n",
       " '299 AAG',\n",
       " '580 ACA',\n",
       " '368 UAU',\n",
       " '1086 GCA',\n",
       " '339 GAA',\n",
       " '365 UCU',\n",
       " '400 GUA',\n",
       " '366 GUC',\n",
       " '247 UAU',\n",
       " '1052 CCU',\n",
       " '108 ACU',\n",
       " '433 AUA',\n",
       " '717 UUU',\n",
       " '790 ACA',\n",
       " '421 AAU',\n",
       " '924 AAC',\n",
       " '643 CAA',\n",
       " '973 UCA',\n",
       " '314 ACU',\n",
       " '620 CCU',\n",
       " '721 GUU',\n",
       " '995 UUG',\n",
       " '549 GGU',\n",
       " '669 AUA',\n",
       " '678 AAU',\n",
       " '402 AGA',\n",
       " '568 AUU',\n",
       " '833 AUC',\n",
       " '729 UCU',\n",
       " '1148 AAG',\n",
       " '1226 AUA',\n",
       " '539 AAU',\n",
       " '1083 GAU',\n",
       " '768 GGA',\n",
       " '1085 AAA',\n",
       " '420 UAU',\n",
       " '61 GUU',\n",
       " '1035 CAA',\n",
       " '1190 AAG',\n",
       " '126 GUU',\n",
       " '554 UCU',\n",
       " '447 AAU',\n",
       " '182 CAG',\n",
       " '19 ACC',\n",
       " '43 AGA',\n",
       " '1045 GGC',\n",
       " '50 ACU',\n",
       " '932 AAA',\n",
       " '994 AGG',\n",
       " '1089 CCU',\n",
       " '645 CGU',\n",
       " '52 GAC',\n",
       " '901 AUG',\n",
       " '836 UAU',\n",
       " '1143 GAA',\n",
       " '450 UAC',\n",
       " '1027 AAA',\n",
       " '612 CAG',\n",
       " '437 UCU',\n",
       " '583 AUU',\n",
       " '1189 GCC',\n",
       " '1059 GUA',\n",
       " '1154 UAU',\n",
       " '865 ACA',\n",
       " '1033 CUU',\n",
       " '82 GUC',\n",
       " '956 CAA',\n",
       " '524 UGU',\n",
       " '981 UCA',\n",
       " '125 GUU',\n",
       " '597 AUA',\n",
       " '7 UUG',\n",
       " '996 AUC',\n",
       " '904 AGG',\n",
       " '377 AAG',\n",
       " '897 UUU',\n",
       " '129 GUC',\n",
       " '432 GUU',\n",
       " '867 GAA',\n",
       " '806 CCA',\n",
       " '1265 AAA',\n",
       " '905 UUU',\n",
       " '1063 CAU',\n",
       " '238 CAA',\n",
       " '20 AGA',\n",
       " '200 UUU',\n",
       " '655 GUC',\n",
       " '1212 CCA',\n",
       " '273 ACU',\n",
       " '22 CAA',\n",
       " '1211 UGG',\n",
       " '792 CCA',\n",
       " '159 UAU',\n",
       " '862 CCU',\n",
       " '519 GCA',\n",
       " '966 AGC',\n",
       " '914 GUU',\n",
       " '572 ACU',\n",
       " '1094 UUU',\n",
       " '208 CCU',\n",
       " '472 UAU',\n",
       " '1120 UUU',\n",
       " '386 UUA',\n",
       " '684 CGU',\n",
       " '687 GCU',\n",
       " '575 GUC',\n",
       " '358 AGC',\n",
       " '5 GUU',\n",
       " '630 CCU',\n",
       " '587 ACA',\n",
       " '122 GCU',\n",
       " '982 CGU',\n",
       " '356 AGA',\n",
       " '1127 GUU',\n",
       " '1050 UCC',\n",
       " '324 UCU',\n",
       " '351 GCU',\n",
       " '776 AAC',\n",
       " '739 AUG',\n",
       " '584 CUU',\n",
       " '771 GUU',\n",
       " '461 AAA',\n",
       " '540 UUC',\n",
       " '1044 AAG',\n",
       " '1123 GGU',\n",
       " '564 UUU',\n",
       " '662 GAC',\n",
       " '1136 GUU',\n",
       " '906 AAU',\n",
       " '895 AUA',\n",
       " '220 UCG',\n",
       " '761 CAA',\n",
       " '186 AAA',\n",
       " '1203 GGA',\n",
       " '203 UAU',\n",
       " '192 GUG',\n",
       " '614 GUU',\n",
       " '983 CUU',\n",
       " '281 AAU',\n",
       " '886 ACC',\n",
       " '978 GAU',\n",
       " '88 GGU',\n",
       " '492 CAA',\n",
       " '949 GAU',\n",
       " '1170 GGC',\n",
       " '144 UAC',\n",
       " '527 AAA',\n",
       " '542 UUC',\n",
       " '893 UUA',\n",
       " '1090 CGU',\n",
       " '679 UCU',\n",
       " '819 GAU',\n",
       " '1216 UGG',\n",
       " '741 AUU',\n",
       " '114 CAG',\n",
       " '663 AUA',\n",
       " '785 AAA',\n",
       " '1121 GUG',\n",
       " '566 AGA',\n",
       " '735 GUA',\n",
       " '332 ACA',\n",
       " '677 ACU',\n",
       " '164 AAU',\n",
       " '1256 GAU',\n",
       " '773 CAA',\n",
       " '41 GUU',\n",
       " '423 AAA',\n",
       " '0 AUG',\n",
       " '968 AAU',\n",
       " '1219 UUU',\n",
       " '1157 AAU',\n",
       " '304 UCC',\n",
       " '786 CAA',\n",
       " '139 UUU',\n",
       " '794 AAA',\n",
       " '698 CUU',\n",
       " '1115 ACU',\n",
       " '907 GGU',\n",
       " '1081 UGU',\n",
       " '649 UUA',\n",
       " '1113 AUC',\n",
       " '664 CCC',\n",
       " '820 CUA',\n",
       " '465 AGA',\n",
       " '1215 AUU',\n",
       " '278 UAU',\n",
       " '675 ACU',\n",
       " '1269 UUA',\n",
       " '370 UCC',\n",
       " '756 GGC',\n",
       " '688 AGU',\n",
       " '319 GUC',\n",
       " '406 GUC',\n",
       " '516 CUU',\n",
       " '980 CUU',\n",
       " '413 CAA',\n",
       " '954 AAU',\n",
       " '1202 CUU',\n",
       " '942 AGU',\n",
       " '1272 ACA',\n",
       " '710 UCU',\n",
       " '493 UCA',\n",
       " '1041 UUU',\n",
       " '223 GAA',\n",
       " '154 AGU',\n",
       " '310 GGA',\n",
       " '1065 ACU',\n",
       " '726 CUA',\n",
       " '1229 GUG',\n",
       " '382 UCU',\n",
       " '856 GGC',\n",
       " '950 GUG',\n",
       " '96 AAG',\n",
       " '491 UUA',\n",
       " '28 ACU',\n",
       " '376 UUU',\n",
       " '495 GGU',\n",
       " '646 GCA',\n",
       " '1165 UUA',\n",
       " '507 UAC',\n",
       " '45 UCA',\n",
       " '86 AAU',\n",
       " '746 ACU',\n",
       " '706 UAC',\n",
       " '1233 CUU',\n",
       " '276 UUA',\n",
       " '360 UGU',\n",
       " '1213 UGG',\n",
       " '877 UUA',\n",
       " '106 GGU',\n",
       " '963 AAA',\n",
       " '136 AAU',\n",
       " '411 CCA',\n",
       " '101 AGA',\n",
       " '1198 GAU',\n",
       " '228 UUG',\n",
       " '522 ACU',\n",
       " '556 AAA',\n",
       " '1043 GGA',\n",
       " '317 UUU',\n",
       " '18 ACA',\n",
       " '1222 GGC',\n",
       " '209 AUU',\n",
       " '1028 AUG',\n",
       " '21 ACU',\n",
       " '172 CAG',\n",
       " '586 AUU',\n",
       " '1 UUU',\n",
       " '869 AUU',\n",
       " '1128 GUA',\n",
       " '514 UUU',\n",
       " '565 GGC',\n",
       " '854 UUU',\n",
       " '1111 CCA',\n",
       " '917 GAG',\n",
       " '231 GGU',\n",
       " '506 CCA',\n",
       " '596 GUU',\n",
       " '1102 UUU',\n",
       " '1174 UCA',\n",
       " '1271 UAC',\n",
       " '1122 UCU',\n",
       " '260 GGU',\n",
       " '759 UGU',\n",
       " '658 UCA',\n",
       " '1166 GGU',\n",
       " '300 UGU',\n",
       " '477 ACA',\n",
       " '1075 ACA',\n",
       " '453 AGA',\n",
       " '1114 AUU',\n",
       " '255 UCA',\n",
       " '609 GUU',\n",
       " '725 AUU',\n",
       " '387 AAU',\n",
       " '631 ACU',\n",
       " '962 GUU',\n",
       " '446 GGU',\n",
       " '267 GGU',\n",
       " '481 GGU',\n",
       " '151 UGG',\n",
       " '579 CAG',\n",
       " '913 AAU',\n",
       " '452 UAU',\n",
       " '318 AGA',\n",
       " '1259 GAC',\n",
       " '603 ACU',\n",
       " '1034 GGA',\n",
       " '344 ACC',\n",
       " '1247 UGU',\n",
       " '415 GGA',\n",
       " '602 AAU',\n",
       " '434 GCU',\n",
       " '782 GCA',\n",
       " '224 CCA',\n",
       " '1049 AUG',\n",
       " '720 AGU',\n",
       " '1092 GGU',\n",
       " '295 CUC',\n",
       " '574 GCU',\n",
       " ...}"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "set(data.indorig)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.ratio });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "    if (this.ratio !== 1) {\n",
       "        fig.send_message('set_dpi_ratio', { dpi_ratio: this.ratio });\n",
       "    }\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    var resizeObserver = new ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    resizeObserver.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = 'image/png';\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.which === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.which;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which !== 17) {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    if (event.altKey && event.which !== 18) {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    if (event.shiftKey && event.which !== 16) {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data']);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.one(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA7YAAAK8CAYAAADbFiNKAAAgAElEQVR4nOzdd3icd522/XEc99hxekIIOCGkLiEFshuSAA/FG/JCyhMCKRAIIQXSWLKbBQKpmCykk+LMjEZlxrbc5V5kFcuSm0ausizLauNuuXdbVrmeP/RaG1mSLckz+s59z/k5jvPYA2kspBCtfpdndN8eAQAAAADgYB7rTwAAAAAAgFPBsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWwBAAAAAI7GsAUAAAAAOBrDFgAAAADgaAxbAAAAAICjMWzhWCtXrtR9992nL3zhC+rTp48+97nP6Xvf+57++c9/tvn4+++/Xx6PRy+++GKb78/NzZXH42nutNNO03nnnaf77rtPq1evbvPPTJkyRd/85jd13nnnqV+/frr00kt1//33a+bMmR3+OrZu3aoXXnhBV155pfr166f+/fvrxhtv1BtvvKHdu3d3+OMcPXpU55xzjm699dZ2H9PY2KjPf/7zuuGGG9p9TGFhoZ5++mldc8016t+/vy655BLdf//9Kisr6/DnAgCAE1RUVOiJJ57QpZdeqj59+mjgwIH6xje+offff1+HDh3qts8jWj/DgUTGsIUjzZ8/X71799bll1+uN954Q36/Xy+//LKGDh2qL33pS60ev3fvXvXt21dDhgzRJZdcosbGxlaPOTZsn3vuOYVCISUnJ+t3v/ud+vbtq3POOUdbtmxp8fi33npLHo9H3/rWt/Tuu+/q008/1X/+53/q+uuv1y9+8YsOfR2FhYU699xz1bdvX/3617/W8OHDNXz4cD322GMaMGCAvv/973fqn8tTTz2lHj16KBKJtPn+uXPnyuPx6J133mn3Y9x333268MIL9eyzz8rv9+uNN97QBRdcoAEDBqi4uLhTnw8AAPFq2rRp6tevnwYPHqznnntOPp9PH330kR544AH16tVLjz/+eLd+PtH4GQ4kMoYtHOnOO+/Ueeed1+YzmjU1Na3elpycrF69eiknJ0cej0dz585t9Zhjw3bcuHEt3j58+HB5PB79/e9/b35bXV2dBg0a1O7wbOtzON7u3bt18cUX64ILLlBpaWmr92/dulVvvPHGST/OZ+Xn58vj8ejNN99s8/1PPPGETjvtNG3atKndjzF//nzV1ta2eNvatWvVp08fPfzww536fAAAiEdVVVU644wzdNVVV2nz5s2t3l9eXq7333//lP97Dhw40OHHRuNnOJDIGLZwpCuvvFLf/va3O/z47373u7rzzjslSVdffXWbfwvb3rBdtWqVPB6Pnnjiiea3bdmyRR6PR6+++moXvwLpf/7nf+TxeDRy5MgO/5mPP/5Y11xzjXr37q2LLrpIv/3tb1uM+8bGRg0ZMkRf+cpXWv3Zo0eP6uyzz9Z3v/vdLn2+N954o2688cYu/VkAAOLJU089JY/Ho/nz55/wcdXV1fJ4PEpJSWn1Po/Ho1deeaX5P7/yyivyeDwqKSnRgw8+qMGDB+v6669vfn8oFNKNN96ovn376qyzztJPf/pTrV+/vvn9nfkZvmLFCv3iF79ofgn1BRdcoEcffVQ7duxo9Wc3btyoX/3qV7rooovUu3dvDRkyRE899VSrv8QGnI5hC0caOnSoBg4c2KGXxm7atEmnnXaaQqGQJOn111/XWWed1er/obc3bKdNmyaPx6P//u//bn5bQ0OD+vXrp5tuukk7d+7s0tfwjW98Q/369evwD5ZjPzC/973v6cMPP9Qzzzyjnj176utf/7qOHj3a/Lg//elP8ng8WrVqVYs/P2XKFHk8HiUnJ3f6c21sbNTFF1+soUOHdvrPAgAQby6++GJddtllJ31cV4btNddco7vvvluffPKJPv74Y0nSX//6V/Xo0UM//elP9cknn+i1117TueeeqyFDhrT4C+qO/gx/++23dfvtt+v111+Xz+fT888/r379+unmm29u8etWmzZt0uc+9zn1799fv/vd7/Tpp5/qL3/5i66++upOXccDcAKGLRwpMzNTPXv2VM+ePXXLLbfoxRdf1OzZs1sMvGPefvtt9evXT/v27ZPU9LJaj8ejjIyMFo87NmyTk5O1fft2bd68WbNmzdLll1+uHj16qLCwsMXjX375ZXk8Hg0YMEA/+MEPNGzYMC1ZsqTDX8NZZ52lr371qx167LZt29S7d28NHTpUDQ0NzW//6KOPWo3VkpISeTwe/fGPf2zxMR544AH17dtXe/fu7fDneEwoFJLH41EgEOj0nwUAIJ7s3btXHo9Hd99990kf25Vh++CDD7Z4XCQSUc+ePTVs2LAWby8uLtbpp5/e4u0d/Rne1oWt0tPT5fF4NG/evOa3PfLIIzrttNMUDodbPb6t640ATsawhWMVFhbq3nvvVf/+/ZuvZHzeeedp8uTJLR53ww036Mc//nGLt910002t3nb8VZE/+zFHjBjR5ucwatQo3XbbbTrttNOaH3/DDTe0exXlz+rZs6duu+22Dn2to0aNksfj0YwZM1q8vba2VoMGDdJ9993X4u033HCDLr300ub/fODAAfXv37/V19wRpaWlGjRokG655RbV19d3+s8DABBPNmzYII/Ho5/97GcnfWxXhm1eXl6Lx7377rvq0aOHysvLtX379hZdffXV+t73vtfi8Z39GX748GFt3769+XM99rvBDQ0NGjRoUIcGPOAGDFs4Xm1trQoLC/XHP/5Rffv2Va9evVRSUiJJWr16tTwej9577z2Vl5c398ILL7R69vLYsH355Zc1Z84cZWRk6JFHHlHv3r01atSoE34Oe/fuVWZmph566CF5PB596Utf0uHDhyU1/T7uZzv2t6ydecb2zTfflMfjUWVlZav3XX/99fra177W4m1vv/12i98dGjFiRItnqevr61t9Xm29JHrLli267LLLdMkll3CxCgCAK8T6GdvP/t6sJP3mN79p8y/Oj3Xddde1ePzJfoZL0s6dO/Xcc8/p/PPPb/XxXnvtNUlNF6H0eDx66aWXOvqPBnA0hi1cJSUlpcVFnY79rkp7ffYlvO39ju2Pf/xjDRgwoNUPqvb84he/aHHl5eP/O4/9cLzllls6/Du2nR22x36v+Omnn5bUdBXpz/5e8bEf1J8tNze3xcfYs2ePrr/+ep199tnNf1EAAIAbfO5zn2vz9oDHi0QibQ7b+vr6doft9u3bWzz2ySefVI8ePTRr1izNmTOnVQsXLmzx+JP9DJeaLorZr18/vfzyy5o4caIyMzM1a9asFp8TwxaJhmELVykuLpbH49GTTz7ZfHXB//N//o/GjRvXquuuu67FFYLbG7YVFRXq2bOnnnzyyQ59Dh9++KE8Ho/S09MlqdUPsGO3Ffjb3/4mj8dz0meDpRO/FPnMM89s9VJkqemH3vnnn6+tW7e2uh/f4cOHW31eu3btavH+22+/Xf3799eCBQs69HUDAOAUTzzxhDwez0l/xh17dve9995r8fbKysoOD9t//OMf8ng8Kisr6/Dnd6Kf4bt27WrxzOwxx64hcuxz4qXISDQMWzhSTk5Omxc9+Pvf/y6Px6N33323+X5wwWCwzY8xbNiwFveDa2/YStJPfvIT9enTR1u2bJEkHTx4sN0fhnfeeac8Ho+WLl16wq9h165duuiii3TRRRe1+cOupqam+T62xy4edccdd7T4uj/55JN2r3ScnJzc/FKr9u7d25b6+nrdddddOv300zV9+vQO/RkAAJykoqJCAwYM0DXXXKOtW7e2+f5jv6t67rnn6t57723x/hdeeKHDw/bYX5A/9NBDrc4ujY2Nbd6i50Q/w4+N7eNvOfjb3/621efExaOQSBi2cKRrr71Wl156qX7/+9/L5/Ppo48+0kMPPaSePXs2Xzr/qaeeUs+ePdu9Hc+xZ3ffeecdSScetuFwuMUtf7Zv3y6Px6N/+7d/06uvvqpAIKB33nlHt99+uzwej+65554OfR2LFi3S2WefrX79+unxxx/Xp59+qk8//VRPPPGEBg4c2OL2Osd+YA4dOlQfffSRnn322TZv93PM3r171bdvX3k8Hl1yySUd/gH2/PPPy+Px6Ec/+pFCoVCrAABwg8mTJzffU/b555+X3+/Xxx9/rIcffli9e/duvn/9H/7wB3k8Hj322GMaPny4HnzwQd10000dHrbS//5K0Te+8Q394x//0PDhw/Xiiy/qy1/+st56661Wjz/Zz/BvfvOb6t+/v1566SV98sknuueee/TVr3611ee0ceNGXXjhhc23+/F6vXr11Vd17bXXcrsfuA7DFo40c+ZM/epXv9JVV12lM844Q71799bll1+uZ599VjU1NTp69KjOOecc3X777Sf8OJdeeqluuOEGSScetpL07W9/W4MGDdKePXtUV1cnv9+ve+65R1/84hfVp08f9e/fXzfccIPeeuutTt30fPPmzfqP//gPXXHFFerbt6/69++vm266ScOGDWt1a56PPvpIV111lXr16qULLrhAv/nNb074g+n++++Xx+PRiy++2OHP51vf+tYJfy8ZUVJfLx06JB04IO3dK+3eLe3YIdXUSFu2SJs2SevXS9XVUmWlVF7e9H/XrWt6X02NtHNn0589eFCqrZX423cgbjQ0NuhI3REdqD2gvUf2auehndp2YJu27N+iDXs3qHp3tSp3VapiZ4Uqd1WqclelqnZVKbI7ovV71mvj3o3asn+Lth3Ypj2H9+jQ0UNqaGw4+X8xOmXt2rV6/PHHNWTIEPXu3VsDBw7Urbfeqg8//FBHjhyR1HRrnccee0xnnnmmBg4cqJ/85Cfatm1bp4atJE2YMEG33XabBgwYoAEDBuiqq67S008/3e5LlE/0M3zjxo269957NXjwYJ155pm6//77tXnz5lafkyStW7dOjzzyiM477zz16dNHl112mZ5++ulOnVUAJ+CUCgDRUlcn7dvXNDojEam0VFqyRCookObMkSZPlkaPllJSJK83Nvn9TR9/xAhp4kRp9mwpP19aulQqK5M2bmwa0RxogE5rbGzU/tr92rRvk9buWKsVW1do8cbFyovkaVb5LE0qnaT04nSlLEuRt8gbkwJLAwqtCGnMqjGaVDpJM8tnKrc6V4UbC1WyrUTr9qzTzkM7daTuiPU/LgDoVgxbAOis/fulDRuk4uKm0Th1qhQKxW6sxqrkZGnMmKbPf+5cacWKpq/rwAHrf8KAmbqGOm0/uF2Vuyq1bMsyzYvM0/S105VenC7/En/MBmusRvCYVWM0rWya8tflq7imWBv2btCBWr7HAbgPwxYA2tLYKO3aJVVUSEVFUlaWNH68FAjYD9LuKDW16Rnm/Hxp1Spp82bpCM8AwV32HdmnyO6IlmxeojmVczS6eLR8RT7zQdodpSxLUUZphnKrc7VsyzKt27NOB48etP6fBAC6jGELAFLTS3PXr5fCYWnatNi+XNjJhULS9OnS4sVNL7c+fNj6fzmgQw7UHlDlrkrNXz9fk9dMjunLhZ3cyJUjlVmRqeVblmvTvk06Wt/64oQAEI8YtgAS065dTb8DO3du08txrQejkxszpumfY2lp08WsAGONjY3acXCHSraVKLsqWyNXjjQfjE7NV+TT2FVjNbd6rtZsX6P9tfut/+cFgDYxbAEkhh07pGXLmp5t5NnY2DZihJSdzdBFt9p1aJdWbF2h6Wun82xsjEsvTldeJE8VOyt0uI5XbQCIDwxbAO505EjT78fm5jrzwk5uauTIpitDb9woNXC7EkRHbX2tqnZVKS+SpxErRpiPvURuXMk4LVi/QOv3rFd9Q731vxoAEhTDFoB77NkjLV/edNEjn89+0FHrUlKans2trJSO8rt76Jydh3Zq6ealmrxmsuOuUJwoJS9N1pzKOSrfWa7a+lrrf2UAJBCGLQBn27pVWriw6f6w1qONOpffL82YIZWUSAe5GivatvPQToU3hTVm1Rjz0Uady7/Er+lrp6tkWwlXXAYQcwxbAM6zd2/TLXjS0+3HGUWviROlpUub7hOMhMaYdWcZpRlauXWlDh09ZP2vGAAXYtgCcIbaWmn16qaXGVsPMIp906ZJ5eVSPb+vlygO1B7Qks1LNHbVWPMBRrHNv8SvmeUzVbmrkt/JBRA1DFsA8auhoeleqZmZUlKS/dii7i8lRcrPl7Zts/63ETHQ0Nigyl2VmrF2hnxFPvPBRd1fyrIU5UXytGX/Fut/HQE4HMMWQPzZvr3pKrppafbDiuKnceOklSulw9xexOl2HdqlhRsWKm15mvmwovgpvThdy7Ys4xZCALqEYQsgPjQ2SlVVvNSYTp7f3/Qs/vr11v/WohPqG+pVtqNMGaUZ5gOK4rukJUnKrc7VtgO8UgNAxzFsAdiqrW16Fm7UKPvBRM5rzBiptJTfxY1jh+sOq2hTkYLLg+aDiZxXRmmG1u5Yy+/iAjgphi0AG3v3SvPnS8nJ9uOInF8oJC1b1vQXJYgLuw7tUl4kT0lLkszHETm/4PKgCjcWctsgAO1i2ALoXlu2SLNnSz6f/Rgi95WcLC1YwC2DDG3Yu0HT1043H0LkzpKWJCl/Xb72Hdln/a86gDjDsAXQPSorpQkT7IcPJUZ+v5SdLe3YYf1vfsKo3l2t8SXjzYcPJUb+JX7lVudq9+Hd1v/qA4gTDFsAsbVunTR+vP3QocRt+nRp0ybr7wTXYtCSZb4in+ZUztGOg/wlFpDoGLYAYmPTJmnSJPtRQ3Ss6dN5BjeKqndXa8LqCebDhuhYM8tnavvB7dbfGgCMMGwBRFdNjTRtmv2IIWqv7GxpH7+f11WR3REGLcV12VXZ/A4ukIAYtgCiY+dOadYs+9FC1JH8/qarch8+bP2d4xg1B2o0qXSS+Wgh6kj+JX7NXz9fh+v4HgcSBcMWwKnZs0fKyrIfKkRdKSVFWrJEqquz/k6KW/tr9yu7Ktt8qBB1peSlyVqyeYnqGvgeB9yOYQuga+rqpEWLmp75sh4nRKdaKCStXi01NFh/Z8WNo/VHVbixkPvQkisKrQhp9bbVamxstP7WAhAjDFsAnVdZKY0YYT9GiKLdmDHSxo3W32GmGhsbVbq9VKEVIfMxQhTtJqyeoJoDNdbfZgBigGELoOP27Gm6sqz1+CCKddnZCfn7tzUHarh1D7k+X5FPeZE8Hak7Yv0tByCKGLYATq6uTgqHpaQk+8FB1F2lpkqlpdbffd2itr5WBesK5CvymY8Oou4qbXmaSrcnxvc4kAgYtgBObN06adQo+5FBZNWUKdLu3dbfiTFTtauKlx1TQpdRmsH9bwEXYNgCaNu+fdy+h+hYfn/Tqxbq662/M6Nmf+1+zSqfZT4qiOIhX5FPCzcsVH2De77HgUTDsAXQWkmJFAjYjwmieGv0aGnTJuvv0FPS2Nio4ppiJS9NNh8TRPHWmFVjuLgU4FAMWwD/6+BBacYM+/FAFO/l5kq1tdbfsZ2278g+TV4z2Xw8EMVzviKfFm9czLO3gMMwbAE0qa6W0tLsBwORUxo5Utq61fo7t8PWbF/Ds7REnWjsqrHadmCb9bcugA5i2AKJ7ujRpmefrEcCkRPz+aQlS6TGRuvv5HYdrjus2RWzzUcCkRPzL/GrcGOhGhobrL+VAZwEwxZIZFu2cMVjomg0ebK0f7/1d3Qr6/es54rHRFFowuoJ2ntkr/W3NIATYNgCiaihQVq8uOnZJutBQOSWUlKkykrr725JUl1DnfLX5ZuPASI3lbIsRRU7K6y/vQG0g2ELJJrdu6UJE+xHAJFbmztXqqsz+xbfc3iPxpWMMx8BRG5tXmQeF5YC4hDDFkgkVVVScrL9wZ/I7Y0ZI23f3v3f4ruquEAUUTc0rmScdh/e3e3f4wDax7AFEkFjY9NLj60P+0SJlN8vFRd3y7d4Q2ODFm5YaH7YJ0qkAksDKttR1i3f4wBOjmELuF1tLfemJbIsL6/p99pj5ODRg9yblsiwvEgeV00G4gDDFnCzXbuk9HT7gz1Rojd5snT4cNS/xTfv26zg8qD5wZ4o0Zu8ZrIO10X/exxAxzFsAbfi92mJ4quRI6WdO6P2LV5cUyxfkc/8QE9ETY1cOVI7Du6I2vc4gM5h2AJuw+/TEsVvgYBUXX2K3+KN3MqHKE4LLA2oaldVdH6eA+gUhi3gJvw+LZEzWrKkS9/iR+uPasbaGeaHdyI6cUWbiqL8Ax7AyTBsAbfYt08aPdr+wE5EHSsrq1P3u91fu5/70xI5qMyKTNU12N3TGkg0DFvADXbulEIh+4M6EXWuCROkAwdO+i2+7cA2hVaEzA/qRNS5JpVO0pG6I91wEADAsAWcbssWKSXF/oBORF1rxIimK5i3o2pXlQJLA+YHdCLqWmNWjdH+2v3ddy4AEhTDFnCySERKSrI/mBPRqZWaKtXUtPoWX1WzyvxQTkSn3ogVI7TrUPt/gQXg1DFsAadas0by+ewP5EQUnQIBacOG5m/xok1F5odxIopeqctStXX/VsODA+BuDFvAiZYtsz+EE1H08/uligrNXz/f/BBORNEvsDSgyO6I9SkCcCWGLeA0CxbYH76JKCY1er3aPiuonMKg+QGciGKTr8insh1l1qcJwHUYtoBTNDRIOTnmB28iik2NPq+2zw4qkuVVJMurPMYtkatbvW219ckCcBWGLeAEDQ3SrFnmB28iik2NPq+2ZaY1j9pjFSzmFj9Ebm5VzSrrEwbgGgxbIN41NEiZmeYHbyKKTY1+n2raGLXHWsC4JXJ1K7eutD5pAK7AsAXiWWOjNGeO+cGbiGJTg9+nrXPaH7XHWrSIcUvk5opriq1PHIDjMWyBeNXYKGVnmx+8iSg2Nfq8J3ymlpclEyVWvCwZODUMWyBe5eaaH7yJKHZ99kJRHY2rJRO5u5JtJdanD8CxGLZAPCooMD90E1Hs2jUj1OlRe6xZ4TTzwzcRxa6KnRXWpxDAkRi2QLwpLDQ/dBNR7NozteujNpLlVXW2T1PCqeaHbyKKTf4lfm3Yu8H6NAI4DsMWiCfLl5sfuokodu2f1PmXH7dVVY5fE4pSzA/gRBSbkpcma9uBbdanEsBRGLZAvFi92vzQTUSx6+D4jl8oqiNV5iZpdDjZ/ABORLEpbXma9hzeY306ARyDYQvEg/XrJZ/P/OBNRLHp8JhUrcvyRXXYRrK8qpgbUKgoyfwATkSxaeTKkTp49KD1KQVwBIYtYG3XLiklxfzgTUSxqXZkitZn+aM+ao9VNi9FqWG/+QGciGLTuJJxqq2vtT6tAHGPYQtYOnxYGjXK/OBNRLHpaDCgDVlJMRu1xyrNT1WgyGd+ACei2DRlzRQ1NDZYn1qAuMawBazU10uTJpkfvIkoNtWnJGljViDmo/ZYq+anyRe2P4ATUWzKi+RZn1yAuMawBaxkZ5sfvIkoNjX6vNo6J7XbRu2xFi0Kmh++iSh2FdcUW59egLjFsAUsLFlifvAmoti1c0Z0buvTlTILGbdEbs1X5OMet0A7GLZAd6usND90E1HsOjDRbtRGsryqzvFxj1siF5eyLIXbAAFtYNgC3ammRgoEzA/eRBSbakcmx+S2Pp2tYm5AaUVcKZnIrY1ZNYYrJQPHYdgC3eXAASkUMj94E1Fsagj4u/ViUSdrdX6q/HFwACei2DR97XQ1NjZan26AuMGwBbpDQ4M0ebL5wZuIYldNZpr5mD2+wkUh88M3EcWuxRsXW59wgLjBsAW6Q2Gh+aGbiGLX7ukh8xHbXrPDaeaHbyKKXev3rLc+5QBxgWELxNrGjZLPZ37wJqLYdGhc/D1T+9mqc/waHU42P3wTUWxKWwyEgQUAACAASURBVJ6mA7UHrE87gDmGLRBLhw5JwaD5wZuIYtPRYEDrs/zm4/Vkrc1LVqDIZ34AJ6LYNHnNZH7fFgmPYQvE0vTp5gdvIopNDX6fNmWlmI/WjrZ8Afe3JXJzhRsLrU89gCmGLRAry5aZH7yJKHbtmGV7v9quNDvMuCVya74inzbs3WB9+gHMMGyBWNi6VfL7zQ/eRBSbDo6P79+rba/K3CSFipLMD+BEFJuCy4M6ePSg9SkIMMGwBaLtyBFp5EjzgzcRxab6ZL82ZCWZj9Sutmo+V0kmcnPT1063PgkBJhi2QLTNnm1+8Cai2LV9tvNegnx8OYW8JJnIzZVsK7E+DQHdjmELRFNpqfmhm4hi18EJznwJ8vFV5fg1sihgfvgmotiUvDRZ+47ssz4VAd2KYQtEy8GDUkqK+cGbiGKT01+CfHyrC1LND99EseyhPzykiy+/WH0H9FXfAX116Vcu1bMfPNviMS8mv6grv3alevftrb4D+uryGy7XhwUftnjMM+8/oyHXDlGvPr3Uf2B/ffVbXzX/2jrS1LKp1icjoFsxbIFo4SXIRK5umwtegnx8eYt5STK5t9+++1s98/4zen3i63p9wuv6waM/UM/Te+qVMa/IW9Q0avsO6Kt7n7lXr4x5Ra+Nf02Pv/m4PlrwUfPHePLvT6r/oP566A8P6fUJr+uVsa/o8TcfN//aOlpxTbH16QjoNgxbIBoqK80P3UQUuw5MdN+ojWR5VZ3j0+hwsvnhm6i76j+ov37+55/LW+TVpf9yqe587M52H/vJok80+PzBzY93YoGlAe05vMf6lAR0C4YtcKpqa6VQyPzgTUSxqT4lyVUvQT6+0vxU+cL2B3CiWDZ88XD9etivdXqv0/XK2Ff0VuZb8ng8+ul//lSXXXeZBp49UF++8cv6r6T/av4zf0j9gzwejx55+RFdcsUlGnTOIF37jWv18uiXzb+ezjSpdJIaGxutT0tAzDFsgVOVl2d+8Cai2LUt0x0XjDpRBYtD5odvolj08uiX1adfH53W8zT1O6Ofnnn/GXmLvPrvlP+Wx+PRgDMH6JGXH9FLI17Sdx/8rk7vdbpen/i6vEVe/XrYr+XxeHT2hWfryb8/qT+F/qSv//vXNeDMAXo3+13zr60zrdy60vq0BMQcwxY4FZs2mR+6iSh27c9w50uQW5Xt1VjTlyT7FSgaqeTweKWGpyoYnq1gYa5CiwsUWrRYwYVLlTZ/pVLyS5Qyr0yB3Aol5VbKn1Mpf07VZ6qWL7tKvpxqebOr5cuuljc7Im92REm5lQrklStlXplSC1Yrbf5KBRcuU2hRWKHFCxUqzFewMEfB8GylhqcqEB4lb5HffJDQqfXxwo/1RsYb+lPoT7rjl3fojMFn6JWxr+jFwIvyeDy649E7Wjz+4ssv1h2/bHrbY399TB6PRw//6eHm93+04COdMfgMPfzHh7vta4hGyUuTdfDoQetTExBTDFugq+rqpPR084M3EcWm+pQkrc/y24/ObirWV0n2F6UopWhC02hdXKDgwiVKKyhRYG6FvFmRuC0wt0Kp+aVKm79MoUWLFCzMUVp4mpLDo+ULJ5kPFupcV918lW6/93YNmzxMHo9Hj77+aIv33/T9m3TzHTfLW+TV7z/9vTweT4uXJ3uLvBpy7RD94NEfmH8tnS2rMsv65ATEFMMW6KpFi8wP3kQUu7bPcv9LkI8vs/DUr5KcFB6h1MKZCi1eqOCCFUqZt0b+nErzgRqrknIrlVZQotCixUorzFRy0Rh5wz7zEUNtd+XXrtQtP7xFn4Y/1eDzBre6eNQlV1zS/Czu+3Pf1+m9T29x8ahPFn2igWcPbPEsrpPauHej9ekJiBmGLdAV27dLPp/5wZuIYtOR9BTzkWlRRW5AyZ0YZUlFIaWFm0Zs2vxVSsp174DtTL6caqXkr1FwYZGChblKKcqQryhgPmoSrTt+eYde8L2gYVOG6eXRL+uOX96hHj166PmPnpe3yKv7f3+/+g7oqyf+5wm9kfGG7nzsTvXq00t/nfTX5o/xnQe/o8HnD9bzHz2v18a/plvvvlUDzx6od3Oc9Tu2xxpdPFoNjQ3WpyggJhi2QFdkZJgfvIkodm2ec/Jh+/zPfyiPx9Oiyy65oPn9/3rdFa3e/9APv2k+Xk/W4kVtX0jKH05VWngGI/YUSs5bq+DCZQoW5ihQlG4+ctzerXfdqnMuOken9zpdA88aqKtuvqp51B7r3mfu1VkXnKXefXvrsusua/Wy408WfaLv/+z7Gnj2QPUd0FdX33x1831wndrSzUutT1FATDBsgc4qLzc/dBNR7No/qWMXjHr+5z/UFUM+p8Kx/2hu6YR3WgzbB++8rcX7iye/bz5cT1q2V2OKkuULJystPF3BhYuVMq/MfBS6sUBuhYILlitYmMvQpW4rsDSg/bX7rU9TQNQxbIHOqK+XRo0yP3gTUWxqCPg7fM/a53/+Q139pc+3+/5/ve4KPfp/v2M/VDtRefZorZozU/lZC5quJhwH4y+RCsytUHDhUgXDs+UvSjMfQOTeZlfMtj5RAVHHsAU6Y8UK84M3EcWuXTNCHR6Bz//8h+rXt7fOP+dMXXLhubr7Ozdr/sg3Wwzbs888Q2cNGqArhnxOv3ngDpVO+9B8vH626qwUlWZN1rJZWSqYlKe8Cf/b+Kwy86GX6KXMW6Pg4gIlh0ebDyFyX5v2bbI+VQFRxbAFOurIESklxfzgTUSxqS4Y0LpOjMKUvz2rj//yhGb6/qK0N5/TjddcpovPP1urpnygSJZXf/vdw0p78znN8r2s9//4K1147mD9+23Xm4/Zquxkrc6aovD0nBZD9vjmTl0kb1a1+bijYyO3TMHFC5QcHms+iMgdTVw90fpkBUQVwxboqPnzzQ/eRBS7ts0+tdv7rJj0ngb276u/v/DzNt8/6q3/kMfjUV7wrwbPzCapNGuSlszM0bwTjNnjm5q12nzQUeuS89YqtGiRksPjzccRObuKnRXWpysgahi2QEfs3Sv5/eYHbyKKTdG6vc91V3xRv3ngjjbft3rqP+XxeJT25nPdNGb9WpOdoaUzs5U/seNjtkWT5iuJZ23jukBeuUKLFislPNF8JJHzSi9O5/Y/cA2GLdARmZnmB28iil1b5qSe8pAsmfqBzhzYX688/dM23z/+/f+Sx+PRTN9fYjhofSrLGq9ls7KUn9HFMXtcszJXmo836lgp88oULMyTvyjFfDCRcyquKbY+ZQFRwbAFTmbrVvNDNxHFroMTuvYS5Mfv/75Gv/OC8kcM04QPXtRtN16ts888Q0vGv6284F/1+1/epamf/En5I4bJ//pv9YWLztXN1305JoO2MitNxZkzNX9SdMZsiybmK5hdZT7aqOP5cqoVXLhMqeFJ5qOJ4r+05Wk6Wn/U+rQFnDKGLXAyGRnmB28iik2NXq82ZSV3aUz+8Ntf0/nnnKnevU7XhecO1g+//bXm359dMOpN3XzdlzV44AD17nW6hlx8vp78ydCo38d2bfYYLZ2Z3anfm+1Kc2YvNx9r1LVS5pUpxLO4dJLCm8LWpy3glDFsgROpqjI/eBNR7DowIRjDlwXHpuqsJK3OmqLFU+fGdMwe/6xtKs/aOrpjz+Km8CwutVFgaUCHjh6yPnUBp4RhC5zIhAnmB28iil2bs6Jz0ajuqDI7qJWZs1UQpd+d7WyzM4vNxxlFp5R5ZQqG58gb9pkPKoqfFm1YZH3qAk4JwxZoz4YN5oduIopdh8ae+gWjuqOyrPFaMvPE95ztljLyFcjiWVs3FcgrV7AwV75wkvmoIvuSlybrSN0R69MX0GUMW6A9U6eaH7yJKHZtyYzvZ2vLsiZ078uNO9DMOavMxxhFv6TcSoUK8+UrCpiPK7KN37WFkzFsgbbU1Jgfuokodh0eHb/P1pZljVfhtFzzEdtmkwq4r62L8+dUKbR4ofzhVPOBRTalLkvlCslwLIYt0JZZs8wP3kQUu2oyu3aLn1i2NnucCqfH6aD9TNOySswHGMV64FYrtKhQSUUh86FF3d+yLcusT2FAlzBsgePt2mV+6Cai2FU7Mr5eglyePVZhBwza5iYvkC8OxhfFPl92tYKLCnkGN8EKLg+qrqHO+jQGdBrDFjheTo75wZuIYte22fHxbG159hiFp8fBRaG60OSsUvPRRd1X00uU58lb5DcfXdQ9FdcUW5/GgE5j2AKftW+f5PebH7yJKDYdDSXHwaBNV9EMZw7aY82dssh8bFH3F8grV1p4tvnootg3YsUINTQ2WJ/KgE5h2AKfVVBgfvAmoti1Y1bQbNBWZaVo+ew55qM0Wk3MLjMfWmRTan6pUooyzMcXxbbKXZXWpzKgUxi2wDGHDkmBgPnBm4hiU11qktaZjFqfVs+ZroIM+zEa1Wdtpy02H1hkW3DBCgXCo8wHGMWmKWumWJ/MgE5h2ALHLFlifvAmoti1a0ao20dtWdZ4LZoSX/eijWbjsteajyuyzZddreDCxfIXpZgPMYp+Ow/ttD6dAR3GsAUkqbFRGjnS/OBNRLGp0efVhqyk7nvZcXaqls7MNh+esS5nRth8WFF8FMitUFp4pvkQo+iWF8mzPqEBHcawBSRp/XrzgzcRxa4DE7rnd2urs7wqyZqmfJe97PhEhbIqzUcVxU/BBSuUFA6aDzKKToGlAR2pO2J9SgM6hGELSNLs2eYHbyKKXTWZsb/FT3nWaC2e6t6XHbfXzMxi8zFF8ZU/p0pphdnmo4yi04qtK6xPaUCHMGyBgwe5xQ+Ri6tLC8T4WVq/VmbONh+YZk1ZaD6kKD5LKyjh4lIuaNTKUWpsbLQ+rQEnxbAFli0zP3gTUezaMy12F42qyB6lRQn4LO3xjc3iIlLUdr6caoUK8+UN+8wHGnW9dXvWWZ/WgJNi2ALp6eYHbyKKXRuzov+MbXWWV6uypmteHIzKeChr9jLzAUXxXcq8NUouGmc+0KhrzamcY31aA06KYYvEtnGj+aGbiGLXobGpUR+1VdlBhafnmo/JuCojX0lZ1ebjieK87IhCiwvMRxp1vqQlSaqtr7U+tQEnxLBFYsvKMj94E1Hs2j47uldDLs2arIIEuuJxZ5qcVWo/nMgRpRWskr8ozXysUedavW219akNOCGGLRLX4cNcNIrIxdUn+7UuyxelZ2mTE+K+tKdSLve0pU6UlFup1PBU87FGHW9S6STrkxtwQgxbJK4VK8wP3kQUu/ZOic5Fo8qyxmv+ZPvh6ISC2VXmg4mcVWjxQi4s5aD2HN5jfXoD2sWwReIaP9784E1EsWtTVsopjlpfYt/GpwvNnLPKfCiR80rNX62k8Ajz0UYnL7wpbH16A9rFsEVi2rvX/NBNRLHrSPqpjdrqrBSFZ3CBqE43ZYH5SCJnlpRbqbTwDPPhRidu5MqR1ic4oF0MWyQmXoZM5Op2zej6y5Arskdp4ZQ4GIkObUw297SlrhdaVChvkd98wFH7bd632foUB7SJYYvElJFhfvAmotjV1XvXrsnKUP5E+3Ho5ObMXm4+jsjZpeaXctXkOG7++vnWpzigTQxbJJ6DB80P3UQUu46Gkrvw0mOvVmbOMh+FriijgHva0ikXyCtXcnis+Yij1vFyZMQrhi0Sz6pV5gdvIopde6Z27t611VkBLZmZYz8IXdTErDLzYUTOz59TpbTwNPMhR63bfnC79WkOaIVhi8Qzdar5wZuIYteWOR2/cFRl9ggtnDLXfAi6rczZK8xHEbmnYGGe+ZCjlnF1ZMQjhi0Sy+HDks9nfvAmothUn+zvxP1pJ6ogw34EurIpC83HELmr4MIi7ncbR41dNdb6RAe0wrBFYlmzxvzgTUSxa39GWodGbcmc6fbjz+WNyK40H0PkrtIKSuQvSjEfddTUnsN7rE91QAsMWySWmTPND95EFLu2ZZ582K6czUWiuqOpc0rMhxC5r+S8tQqER5mPOvJq+Zbl1qc6oAWGLRJHba2UlGR+8Cai2NTo92ldlu+EVz5ePnuO+eBLlHJmLjEfQeTOknIrlVKUYT7sEr2M0gzrkx3QAsMWiaOy0vzgTUSx69C4Ez1b69PSmdnmYy+hmpgvP7f9oRjly6lWaniq+bhL5HxFPtXW11qf7oBmDFskjrw884M3EcWuHbPavs1PdVaSimZwOx+LxmWvNR9A5N582dVKC88yH3iJXNWuKuvTHdCMYYvEkZ5ufvAmoti1ISup1aityk5W4fRc84GXqM3OXGk+fsj9BcPZ5gMvUctfl299ugOaMWyRGPbtMz90E1HsOpLe+t61VdmpWjSVe9RaNnfqIvPRQ4kR97q1aXTxaOsTHtCMYYvEUFpqfvAmoti1Z2qoxaitzA5qwRT7YUd5Ssvitj/UPYUWzzMfeonYgdoD1qc8QBLDFokiK8v84E1EsavmM7f5qcgeqfmT7QcdNTUla7X54KHEKbS4wHzoJVplO8qsT3mAJIYtEkUoZH7wJqLY1Oj1an2W//9/pnaECibZjzn637JnLjUfO5RYhRbPNx97iVR2Vbb1KQ+QxLBFItizx/zgTUSxq3Zk0+/XVmaFeKY2HssokJfb/lA3F1q8wHzwJUqhFSHrkx4giWGLRLBmjfnBm4hi177JQVVlB7WAURu3jckuNx86lHiFuKBUt7Xn8B7r0x7AsEUCyM01P3gTUezaMjtdC6dw9eN4blpWifnIocQsGJ5jPvoSobU71lqf9gCGLRIA968lcm31/lQVcvXjuC9r9jLzgUMJWnZEaeHp5sPP7RWsK7A+7QEMW7jcwYPmB28iik0NvoA2jlhlPtro5OVyP1syzJdTrdTwJPPx5+Ymrp5ofeIDGLZwuUjE/PBNRNGv0etXTdoqrRm9xny0UUeaJz8XkCLD/DlVSg6PNR+Abs2/xK/6hnrrUx8SHMMW7rZkifkBnIii3/bgckW8ERVPLI6D0UYdKZ0LSJFxgdwKBYpGmo9At1ZzoMb61IcEx7CFu2Vmmh/AiSi67QqFFfFGFPFGFJ4UNh9s1LGmZK02HzZEyXnl8helmY9AN1ZcU2x96kOCY9jC3bhwFJGr2hMqaB611d5q5U/MNx9s1LEyZ68wHzVE3qyIUuaVyRdONh+Cbiu7Ktv61IcEx7CFe9XVmR/CiSh6HUyb1TxqI96IKtMqzccadbzcaYXmg4boWGkFJfKGfeZj0E2lF6dbn/yQ4Bi2cK+aGvODOBFFp6PJY7TeX9Vi2HLhKIc1MV9eLiBFcVRoUaH5GHRTviKf6hrqrE9/SGAMW7jX6tXmh/HurH74cP35zjs15Jxz1LdXL1127rl6/a671Pjppy0et/rVV/Wj667ToL591b93b33ti1/UujfflLxeVQ8bJo/H02Zjn3jC/GukxKzBn6JNyWtbjNqIN6JVE7jVj9MakV1hPmaIPlswPMd8ELqp7Qe3W5/+kMAYtnCv/HzzA3l3Nuzuu3XOgAGa9swzqh42TOOeeEJn9OmjD3760+bHVPz1rzp7wAD919ChWvrSS6r46181+be/Vc1bb0nepnG85R//aNFrP/qRzujTR/s/+MD8a6TEbFvaylajNuKNaFnGMvOhRp1rUvYa8yFD9Nl8OdVKKZpgPgjd0toda61Pf0hgDFu416RJ5gfy7uz/+8pX9Ktbb23xtv97ww16+Oabm//zT7/2Nf3sX/+1Ux/3+ksuafVxibqrPaEFbY7aiDeixRmLzYcada7ZmSvNhwzR8QXmVnCl5Ci1eONi69MfEhjDFu6VnGx+KO/Oht19t754zjkqe/11yevV8j//WecPHKgRv/qV5PWqYfhwndGnj16/6y4NveYanTdwoG4eMkQZv/lNux+z6E9/ksfj0fwXXzT/+ijxOpQ2o91RG/FGuCKyA8udETYfMURtlVqwWt4iv/kwdHqzK2Zbn/6QwBi2cKe9e80P5d1dw/Dh+u9//3f16NFDp592mnr06KG/3XNP8/u3/OMf8ng86t+7t969/34t+/Of9ea996pHjx6a+8ILbX7M33zrW7r6oovMvzZKvI4mj251sagWV0RO4YrIjiyjwHzAELVXaFHYfBg6vdHFo61PgEhgDFu4U3W1+cG8u0v/9a/1+bPOUvqvf62VL7+s4KOP6uwBA5T6y19KXq82/f3v8ng8evDrX2/x53503XV64Li3yevVoQ8/1Jn9+untH//Y/GujxKrBl9zmxaI+W9nIMvuRRl0qmF1lPmCI2isYzjYfh07OV+RTfUO99SkQCYphC3cqLjY/nHd3nz/rLH30wAMt3vbGXXfpygsukLxe1X78sU4/7TS9cdddLR7z4tCh+saXvtTq4wUffVS9evbUtrffNv/aKLHaFlxxwlEb8UZUMrbEfKBR1xqTvdZ8vBC1ly+7WilFGeYD0cntPLTT+hSIBMWwhTstXGh+OO/uzh4wQJ889FCLt/3tnnv05fPPb/7Pt1x2WauLR91z/fWtnsWV16tvXXGF7rvxRvOvixKrPaGCk47aiDeiFRNXmA806loZWaXm44XoRAVyK+QPp5oPRKdWuavS+hSIBMWwhTvNmWN+QO/ufnHLLbp48ODm2/1MfOopnXvGGXpx6NDmx0x86in16tlTvp/9TOVvvKEPH3hAPU87Tfn/9V8tPlb5G2+oR48emvnss+ZfFyVOh9KmdWjURrwRhSeFzQcada3pc0rMhwvRyUqbv9J8IDq1FVtXWJ8CkaAYtnCnjAzzQ3p3t++DD/T8d76jL5x9tvr26qXLzj1XL/3gB6r9+OMWjws88oguP/989e3VS1/9/Oc1qY2rIv/xjjt0yVlnqWH4cPOvixKj+qTQCS8WdXzzM+abDzTqWtzyh5wSv2/bteavn299CkSCYtjCnUIh84M6EXW8mrRVHR61EV/EfJxR18uaudR8sBB1JH9OtQJF6eZD0WllVmRanwKRoBi2cJ+GBvNDOhF1vP1p2R0ftdzqx/HlTi80HyxEHS01v1TesM98LDqpiasnWp8EkaAYtnCfBLyHLZFTq0sa0amXIEe8EZWHys3HGZ1CUxaajxWizhRavMB8LDqp4PKg9UkQCYphC/fZtMn8sE5EHWtrWkmnRm3EG1FZOvewdXbzzIcKUafKjiglPNF8MDqphsYG69MgEhDDFu5TVmZ+WCeik7cvLafTozbijah0bGkcjDM6lVKyquzHClEnSs4rl68oYD4YndLeI3utT4NIQAxbuM+SJeYHdiI6cXWBkVrvr+7SsC0eX2w+zOjUSs8uNx8qRJ0tuHCJ+WB0Spv2bbI+DSIBMWzhPnl55od2IjpxW1NXd2nURrwRrZi4wnyY0ak1MavMfKQQdaW08Ezz0eiEKnZWWJ8GkYAYtnCfWbPMD+1E1H77gnO7PGoj3oiWZiw1H2Z0ak3NWm0+UIi6UlJupXzhZPPhGO+VbCuxPg0iATFs4T5Tp5of3Imo7Y4G0rXO17WXIB8rnBE2H2Z0as3MLDYfKERdLbQobD4c471lW5ZZnwaRgBi2cJ8JE8wP70TUdqfyEuRjLcpYZD7M6NTKmr3MfJwQnUop4fHm4zGeW7RhkfVpEAmIYQv3GT3a/PBORK3bG8o75VEb8UZUMLHAfJjRqZUzo8h8mBCdSqn5pebjMZ7Li+RZnwaRgBi2cJ9QyPwAT0QtqwuMPOWXIEe8EUV8EfNRRlEYttMKzYcJ0akWLMwxH5Dx2pzKOdanQSQghi3cJznZ/BBPRC3bHlwelWdrq5OqzUcZnXq5DFtyQf6cKvmL0sxHZDw2rWya9WkQCYhhC/eJg0M8Ef1vR1IyojJqI96IqpKrzEcZnXpzpy42HyVE0ShtwVLzERmPTVg9wfo0iATEsIW7HD1qfognopZtSS2N2rCtTK40H2UUjWG7yHyQEEWr1KLJ5kMy3hq1cpT1iRAJiGELdzl40PwQT0T/24G02VEbtRFvRJWpDFs3NHfKQvMxQhStUuaVyVvkNx+T8VRwedD6RIgExLCFu+zebX6QJ6KmGn1J2hioiOqwrQhWmI8yisKwnbzAfIwQRbPQ4gLzMRlPpS5LtT4RIgExbOEu27aZH+aJqKk9wYKojtqIN6LyULn5KKMoNHm++RAhimb+nCr5w6nmgzJeCiwNWJ8IkYAYtnCXTZvMD/NE5FW9P03r/VXRH7YjGLauaFKB+RAhinahRYvNB2W85F/itz4RIgExbOEuGzaYH+iJyKudwSVRH7URb0TlIxm2riiDYUvuy5dTraSikPmojJeA7sawhbvwjC2RebXJY2MyaiPeiMpGldmPMjr1Js4zHyFEsSi4sMh8UMZLjY2N1qdCJBiGLdxl82bzQz1RoleTtip2wzadYeuWrAcIUUzKjihQlG4+KuOhhsYG61MhEgzDFu6yZYv5oZ4okTuUNi1mozbijWjN6DXmg4yiNWyr7UcIUQwKLlhuPirjofqGeutTIRIMwxbuUlNjfrAnStQavT5tSimL7bAdw7B1S/451eYDhChWpYTHmw9L6+oa6qxPhUgwDFu4C7f7ITLrQDAzpqM24o1oTTrD1i0FeMaWXFza/FXmw9I6XoqM7sawhbvs2GF+uCdK1DbH+NnaiDeispH8jq1bSmHYkstLLZpsPi6t4nY/sMCwhbvs3Gl+uCdKxA6lTo35qI14I1o7Yq35IKPolMSwJZeXml9qPjCtSl6abH0iRAJi2MJddu82P+ATJWJbU2N3JeTPVh7kPrbuiNv9UGKUGp5qPjItSlueZn0iRAJi2MJd9uwxP+ATJVq1KRO6ZdRGvBFVplXGwSijU25SgfngIOqOEvV3bUesGGF9IkQCYtjCXfbuNT/kEyVa24Mrum3YVqVU2Y8yOvWmLDAfHETdVXLRGPOh2d2lF6dbnwiRgBi2cJf9+80P+USJVF1gVLeN2og3oupAtf0oo1Nu7tTF5mODqLsKLlxiPjS7u3El46xPhEhADFu4S22t+UGfKJHaGSzq3mHr5ZpEPwAAIABJREFUZ9i6oZzpYfOxQdRd+bKr5S9KMx+b3dnE1ROtT4RIQAxbuI/fb37YJ0qE6v1pWuer7tZhG/FGzEcZnXrZM5eajw2i7iy0eIH52OzOpqyZYn0aRAJi2MJ9QiHzAz9RIrQ7uKDbR23EG1H+xHzzYUan1pzZy82HBlF3lpRbKV84yXxwdleZFZnWp0EkIIYt3GfcOPMDP5Hba/AFtCGp0mTYFkwsMB9m0ejxhx/XlV+6Uv369tPgQYN129dvU+ifoRaPyUzP1D133KNBZwxSv7799M1//aYyAhnN7w+8E9B3bv2OzjvnPPXu3VtfuPgLevbRZ82/tpM1O3OF+dAg6u6ChTnmg7O7mheZZ30aRAJi2MJ9pk41P/QTub29wTyTURvxRrQgY4H5MItGN19/s/7w9B+U+l6qAu8E9G83/psuOPcCzRo5q/kxdw29S+efe77ee/U9+f7h0zVXXKN/ufJfmt//4m9f1H133qcPXv9A6R+n66XnXlKf3n30/GPPm399J2rGnGLzkUHU3SXnrTUfnN1V4cZC69MgEhDDFu4zZ475oZ/IzTV6/doYqDAbtoUZhebDLBZNTp4sj8ejf77+T+VNyNP04HSdfvrpeu0/X2t+TPCDoDwejz5585N2P849d9yjG/7lBvOv50RNyyoxHxlEFqWFp5uPzu5o5daV1qdBJCCGLdwnP9/84E/k5g6mzTQbtRFvREszlpoPs1g08qOR8ng8SnkvRXkT8vTeq+/J4/FoWnBai8ddcO4FevqXT7f7cb7/ze/rW//2LfOv50RNzio1HxhEFqUVlJiPzu6ofGe59WkQCYhhC/cJh80P/kRublvaStNhu3LiSvNhFu1yx+Xqlptu0Veu+krz2/78uz+r1+m9Wj32qsuv0oP3PNjmx/l42Mfq2bOn3v7L2+Zf04mamF1mPjCIrAqER5kPz1i3Ye8G69MgEhDDFu5TXGx+8Cdya/X+VK0zHLURb0Ql40vMh1m0u2voXbrwvAs1zjeuy8M25b0UnTnoTD324GPmX8/JGpu11nxcEFkVWrzQfHjGuu0Ht1ufBpGAGLZwn/Jy88M/kVvbF5xrOmoj3ojWjF5jPsyi2b133KvzzjlPoz8Z3eLtnXkpctr7aTrrzLP0s//7M/OvpyONyq4wHxdEVgXyys2HZ6zbX7vf+jSIBMSwhfts2GB++Cdya1tS1pgP2/KR5ebDLBrNHT9X995xr849+1yN/HBkq/cfu3jU6//5evPbQv8Mtbp4VOp7qTrrzLP0wN0PmH9NHS01u8p8XBBZlhqeYj4+Y1l9Q731aRAJiGEL99m50/zwT+TGjiaPMR+1EW9ElamV5sMsGt3973frjP5n6IPXP9DEpInNZY7KbH7MXUPv0gXnXtB8u59rr7xW1155bfP7U95L0eBBg/X9b36/xceYnDzZ/Otrt4n55qOCyLrggqXm4zNWhVaErE+CSFAMW7jP0aPmA4DIje0OLTIftRFvRNX+avtxFoU8Hk+b/eHpPzQ/JjM9U/fccY8GnjFQffv01e3/ersmJk1sfv8vf/LLNj/GheddaP71tdfcqYvMRwWRdf6cavmKAuYjNBZNXjPZ+iSIBMWwhTuFQuYjgMhtWd679vgKJhaYDzTqWtkzisxHBVE8FCzMMh+hsSinKsf6FIgExbCFO02ZYj4CiNzU4dSp5mP2sy3KWGQ+0KhrZWauMB8URPGQW+9pG94Utj4FIkExbOFOeXnmQ4DITW0PLjMfs5+tKKPIfKBR15o2p8R8UBDFS4GikeZDNNqt2b7G+hSIBMWwhTstW2Y+BIjcUoMvoPX+avMx+9lWTFxhPtCoa03MKjMfE0TxUmjxfPMhGu027dtkfQpEgmLYwp2qqszHAJFb2p+WZT5kj2/1uNXmA4261kjuYUvUXLIL72nLPWxhhWELd+KWP0RRa2taifmQPT633Ms2EfNnVZuPCaJ4KrlonPkYjVb+JX41NjZanwKRoBi2cKe6OvMxQOSG6pJGmI/YtqoOuOOWP4nW3CkLzUcEUbwVWrzAfJBGq7GrxlqfAJHAGLZwrxEjzEcBkdPbH8wxH7HtNT9jvvlQo86VM3OJ+YggirdS5pWZD9JolVWZZX36QwJj2MK9pk41HwVETm9b2krzAdte4Ulh86FGnWtW5krzEUEUj7nl6shLNy+1Pv0hgTFs4V7z5pmPAiIn1+j1a72/ynzAthdXRnZek7NKzQcEUTwWKswzH6XRaN2eddanPyQwhi3cq7jYfBgQObnDqZPNx+uJWj2WKyM7rdFZ5eYDgige+3/t3Xl83Hd95/EJJDY5TBJCAykJISnQNqENpbCPZVlykAWyBRK6AZaj2wDtgwJZFnZZulDSBpJmExKgDZBNNKO5LV+yJMuybMsa+ZR86PAtW7Iky7Isy7Z8y5Ksa977x9SKFR3WMTOf32/m9Xw8Xo9HH7ZQfWn0fWdmfr9wZb35KE1GXBEZlhi2yFzHj5sPAyI3dya61Xy8TlbTfK6M7LZyuSIy0fhVHJK3JmA+TGdTcEfQ+uSHLMewReYaGpJ8PvNxQOTWjgYazMfrZLXmcmVkV8UVkYkmLVyz2nyczqbihmLrkx+yHMMWma2gwHwcELmxIV/YfLhOJa6M7J5iq3eYDwciJxfZst18nM6mTW2brE99yHIMW2S2TZvMBwKRG+sOl5uP1qlUV1RnPthoapXE9pkPByIn51/XbD5OZ9O+E/usT33IcgxbZLbGRvOBQOTGusI7zEfrVNpbsNd8sNHUyqtoNh8ORE4vWLPUfKDOtFO9p6xPfchyDFtktjNnzAcCkRtrz20xH61TqSmPC0i5oqJK88FA5Iai2zaaD9SZFNwRVDwetz71IcsxbJH5gkHzkUDkpi4GCswH61Rr9bZqY8FG++FGk1axstZ8MBC5oXDVLvOROpNKD5Ran/YAhi2yQEmJ+VAgclNnI1Xmg3U6VRdVmw83mryV5XvNB4Mbe2HRVv27hz+v6+fdpGvmzNXvv+cP9ZNXlo/8/Ec/9bg8Hs+o7vnw/aM+x2Pf+N+6+54P6Zq5b9G1188z/z3R5OWuazEfqTOp7mid9WkPYNgiC2zbZj4UiNxUZ2i/+VidTrsLdpsPN5q8xRUHzAeD2/p10S7d8o536aOf/oJ+/Ltlei66Sd9/IaJ/jmwY+ZiPfupx3fuRB/TikuqRfl20a9Tn+cx/+x/64ref0n/6wt8ybF2Sv3ah+VCdbkfOHbE+7QEMW2SB1lbzoUDkloa9AbU5YKxOp8aFjebDjSapcJN8sVbzseC2Pv1fv633fuAjk37MRz/1uO77D5+c0ud74kcvMWxdUrg6Zj5Up5O31quBoQHr0x7AsEUW6OkxHwtEbqkvtNx8qE631txW+/FGE7autNp8KLix2979Xj38X76pD93/F5p30y264w/u0V/9z+dHfcxHP/W4rr1+nubddIvecfvduv+zX9OvCsa/XzDD1j1FttSZj9XptLR+qfVJD5DEsEW2WLzYfDAQuaFz0Y3mQ3UmbS3aaj7gaPxWr9ltPhTc2NXXzNHV18zRI1/5rn766gp97Qf/V9fMmauv//0vRz7mb3/6G333GZ/+ybta3/l5jm5793v1nj+8T6+WtYz5fAxb9xTYcMB8rE6nyrZK61MeIIlhi2xRWWk+GIjcUFdkp/lInUm7CneZDzgav4KKRvOh4MbefPU1uvueD436sYc+/3Xd9cd/NuH/5p8jG+XxePSDF/PG/BzD1l35akLmg3WqHTx90PqUB0hi2CJb8D5boinVEThgPlJnUsPiBvMBR+MXqDhoPhLc2NtufZc+9p//66gf++r/eFY33fKOSf93N9z4Nn3tB8+N+XGGrbsK16w0H6xTyVfnU/9Qv/UpD5DEsEW26O+XfD7z0UDk5Ia9fvOBOtNagi3mA47Gtr5kq/lAcGsfeejRMRePevi/fHPMs7iX98LCLbrqqqv03Wd8Y36OYeuuolu3mY/WqVTcUGx9wgNGMGyRPZYtMx8ORE7uYmiZ+UCdTbzP1nmVl+00Hwhu7Se/K9ab3ny1Hvvmj/RseL3+5icva85brtU3f/Kvyokd0ssl9frkF7+l//ObQj03f5N+8GKe3v2+D+jWd92l3618/eXfz+dV6anXSvW5J/6n5l57vZ56rVRPvVaql0vqzX+PNHGhTfvNR+tU4v61cBKGLbJHba35cCBycucjG8zH6WzaXcj9bJ1WUUWD+UBwc0/+s1+//54/1NXXzNE73/0Ho66K/NvSBt3z5x/XvJtu0Zuvvka3vONd+o9/8RW9lF8z6nN89FOPy+PxjOl//XKh+e+PJs67tlU5NV7z4XqlTlw4YX26A0YwbJE9TpwwHw5ETu5kZLv5OJ1NB/IOmA85uqzCTcrl/rVEM85fu9B8uE5WeGdY8Xjc+nQHjGDYIrtEIubjgcipdQQazcfprPIeUmVhpf2gI20o2KCKVdvNhwGRmwtVrzIfr5NVcbDC+lQHjMKwRXZZt858PBA5sbg3V23WwzQJbS/abj7oKNEyXoZMNKui2yrNx+tkNZ5stD7VAaMwbJFdWlrMBwSRE7sYLDQfpclo/+L95oOONmhDwUb5eRky0ayKbN5pPl4nylvrVe9Ar/WpDhiFYYvswm1/iMbtfGSd+ShNRgf9B7WxYKMDhl12t3ZVnfkoIHJ7wY2N5gN2opY3LLc+0QFjMGyRfUpLzUcEkdM6GakzH6XJqmZZjfmwy/aWV+w3HwVEbs/JV0bee3yv9WkOGINhi+zT2Gg+Ioic1tFgg/kgTVb78veZD7vsbqMCFQfNRwFRJuTEKyN7a73qGeixPs0BYzBskX0GBqTcXPMhQeSkDvtazQdpsjoYOOiAcZe9rS2tMR8DRJlSqGal+ZB9Y8UNxdYnOWBcDFtkp/Jy8yFB5JSGfGHzMZrsqouqzQdetlYS22c+BogyJSdeGXnP8T3WpzhgXAxbZKfWVvMxQeSULgYKzIdosqvPrzcfeNlaMMbLkImSVXjzDvMh+8Z4GTKcimGL7DQ8LAWD5oOCyAldiJSZD9Fk1+pv5erIBq3jZchESc1pV0bmZchwMoYtsteGDeaDgsgJnYtWmg/RVFRXVGc+9LKtFbF68yFAlEnlrmsxH7O8DBluwbBF9uroMB8URE7oVAbd6ufyGhc2mg+9bCscazEfAkQZV63PfNDm1ObIV+dT32Cf9ekNmBDDFtkrHpeiUfNRQWTd8fBe8xGakryHVFVYZT72sqV1pdX2A4AoA/PX5pmP2pzaHJW3lFuf3IBJMWyR3bZsMR8VRNZ1BA7Yj9AUtbtwt/ngy5aKY/vNBwBRJhasLTIftTm1OTp89rD1qQ2YFMMW2a2ry3xUEFl32HfQfICmquZIs/ngy4qWVcoXazUfAESZWLhmtfmozdudp3g8bn1qAybFsAXy882HBZFVcW+u+fhMddzTNvWtXrPb/PBPlKlFqtebD9uajhrr0xpwRQxbYN8+83FBZNWgf4H58Ex1+/L3mQ+/TC9Swb1riVJVZOsW82F7/uJ569MacEUMW2BwkHvaUtbWF1xmPjxTXWtuqzYVbjIff5na2lV15gd/okwusnm76agtaSyxPqkBU8KwBSQuIkVZ24VImfnwTEc7inaYD8BMbWms0fzgT5TJhSv3mg7bplNN1qc0YEoYtoAknT8veb3mI4Mo3Z2LbjAfnemoKa/JfABmYuuXbzE/9BNlesGNjWajNrorquH4sPUpDZgShi1wyerV5iODKN2diW5N6oBc/MPFevhPHtatN94qj8ejnO/kjPr5xz/6uDwez6juv+f+UR+z89c79di/e0w3vOUGzbt2nr70sS+p/uX6Wf/athVtMx+CmdaKWL35oZ8o0/OtbTEbttuPbrc+nQFTxrAFLjlyxHxkEKW709GapA7b4PeC+u//+b/rtW+/NuGwfeDeB1T9YvVIu369a9THPHDvA/rj2/9YRT8uUv6P8vWeW9+jRz/y6Kx/bQ2LG8yHYEZVuEl+bvFDlJ5qvGkftf7tfvUN9lmfzoApY9gCl1uyxHxoEKWzU5G6lL38d6Jh+8n7Pjnh/yb2s5g8Ho+W/2T5yI+FvhfSVVddpW2/2DarX09rTqu2FG2xH4QZUnnZTvvDPlGW5K31p33Ybjy00fpUBkwLwxa4XH29+dAgSmddkR1pH7bzrp2nW+bdorvfcbe+dv/XtONXr/8aXvzrF/XW69466n/T/P+a9eY3vVmvffu1Wf+a6vPrzQdhprSgotn8sE+ULflqw2kftmf6zlifyoBpYdgCl+PWP5RldUV2zXosTmfY/uZvfyPfd31a/U+rlfOdHL33tvfqvvfcp5ZXW3Qo55B+9NiPdPc77h7zuW6Zd4ue/eqzs/41tfpaVVlYaT4K3d7aFdXmB32ibMpfm5fWUVt6oNT6RAZMG8MWeKPNm83HBlG6OhHendZh+8Y2/vNGeTwe5f0gLy3D9lDOIe0p2GM+DN1eUazB/KBPlE0Fahalddi2n2u3Po0B08awBd6ou1vy+cwHB1E6Ohae/dWGZzNsD+Uc0ttueJue+9pzOpST+pciH8o5pIOBg9pUuMl8HLq24ip5HXDQJ8qmgjVL0zZql+xdYn0SA2aEYQuMZ+NG88FBlI6OhfaZDtstL2zRVVddJd93fTqU8/rFo0r+oWTkYyLfjyTl4lGXt7Nop/1AdGklsX3mh3yibCtYW5S2Ydt0qsn6FAbMCMMWGE93t5Sbaz46iFLd0eD+pI7Z+pfrVfpUqUqfKpXH49FTX3xKpU+Vqur5KtW/XK9vffJbKvw/hdr03Cbl/SBPH3j3B3TXrXep8XeNI5/jgXsf0L133KtlP16mpT9aqrtuvSspt/u5vJZwi/lAdGXLN/NsLZFBoZqStIzaxXsXKx6PW5/CgBlh2AITqaw0Hx1Eqa4j2JjUwbjwfy2Ux+MZ0+MffVwNv23Qx+/5uG6Zd4uuefM1etct79JX/uNXVPPS6Hvp7vz1Tj36kUd1/dzrNe8t8/TF//BF1b+c/JdM1xXV2Q9Fl1Uc229+wCfKxsI1K3m2FrgChi0wkQsXeNaWMr4j/gNJH4xuqSnaZD4U3dT6kq3KibWaH/CJsrFITRnP1gJXwLAFJlNVZT48iFJZu7/ZfGBatr1ou/lgdEuFFY3mh3uibC1SHePZWuAKGLbAZHp6JL/ffHwQpar23BbzcWlZS6hFGws2mo9Gp7d+xTbzgz1RNhepXseztcAVMGyBK+G+tpTBtXlbzceldbsKd5kPR6e3NMaztUSWRbat49la4AoYtsCV9PbyrC1lbNaj0glxX9vJW1daY36oJ8r2IjUVKRu1S/Yu4dlaZASGLTAVW7aYDxCiVNSWwzO2h3IOaU/BHvMB6dQWVxwwP9QTZXuRmvKUDdvDZw9bn7KApGDYAlPR1ycFAuYjhCjZHfZl93tsL9Wa26rKwkrzEem0KlbWmh/oieiQwjWrUzJqSw+UWp+wgKRh2AJTtXOn+QghSnZHsvyqyJdXn19vPiSd1sKKJvMDPRGl5j623lqvTvWesj5dAUnDsAWmanhYWrzYfIgQJbOOQPbex/aNtXpbtblos/mYdEoVq7abH+aJKFGopiTpw3bjoY3WJysgqRi2wHS0t5sPEaJk1hlsMB+UTmr/kv3mg9Ipza9oMT/ME1GiUG1xUkdtYHtAvQO91qcqIKkYtsB0lZWZjxGiZHUstM98TDqp1pxWbSvaZj4qrVuzZpf5QZ6IXi9YU5jUYbv96Hbr0xSQdAxbYLrOn5dyc80HCVEyOh7eYz4mnVZTtMl8WJpWvFm5sVbzgzwRvV6gNj9pozZvd56GhoesT1NA0jFsgZmoqTEfJETJqCuy03xIOrFdhbvsB6ZRhbFG80M8EY0uULMoacP24OmD1qcoICUYtsBMDA5KeXnmo4Rotp2M1JmPSCfWmtuqqqIq85GZ7mKr6swP8EQ0Nn9tXlJG7eqm1dYnKCBlGLbATB08aD5KiGbb6Ui1+Yh0ag2LGsyHZlorqlSo4qD5AZ6IxpZbE0nKBaO6+7utT09AyjBsgdlYscJ8mBDNpjPRLeYD0snVFdXZD840tTxWb354J6IJqvXNetjuOrbL+tQEpBTDFpiNM2ckn898nBDNtHPRSvPx6ORaQi3aVLjJfHSmunWl1fYHdyIaN9/allmP2oJ9BYrH49anJiClGLbAbNXVmY8Topl2PrLefDw6vb35e82HZ2rbqLyKZvPDOxGNX3Bj46xGrbfWq66eLuvTEpByDFtgtuJxqajIfKAQzaTucMx8ODo+7yFtLdrqgAGamlaX7zE/uBPRxIUq985q2FYdrrI+KQFpwbAFkuHMGcnvNx8pRNOtJ7zKfji6oKb5mXlv2/XLt8jHPWuJHF1k885Z3bN2YGjA+pQEpAXDFkiW3bvNRwrRdOsP5JuPRreUife2XVLRZH5oJ6LJi26tnvFLkI+eP2p9OgLShmELJFNJiflQIZpOw76g+WB0S63e1ox6SXL56p3mB3YiunLR6k0zGrZb2rdYn4qAtGLYAsl0/rwUCJiPFaLpdNh30Hw0uqXmSHNmXCV5WZUC3LOWyBWFq8unPWrz6/M1NDxkfSoC0ophCyTb/v3mQ4VoOh0NNpoPRjdVn19vP0xn2dJYo/lhnYimVqimZFqj1lfn06neU9anISDtGLZAKqxaZT5WiKbaifBu87HotuqK6szH6UzjKshE7ipQmz+tYbuzc6f1KQgwwbAFUqGnRwqFzAcL0VQ6Ha02H4pu62DgoKoKq8xH6nRbt7JGOVwFmchV+WrDUx61xQ3Fisfj1qcgwATDFkiVlhbzwUI0lc5HNpgPRTd2YOEB86E6rYqrFIzxvloiN+WtaJ3yqA3uCOr8xfPWpx/ADMMWSKWqKvPRQnSlesMrzUeiW3PTLYDyKw6YH9KJaHoFNzZOedi2nmm1PvUAphi2QCoND0vFxebDhWiy+gNLzAeiW3PLLYBWlu81P6AT0fSLbN7FrX2AKWLYAqnW0yNFo+bjhWiihr1+84Ho5px+C6C1q2rND+dENLOi27ZM6X21w/Fh69MOYI5hC6RDZ6fk85kPGKKJas9tMR+Ibq5hcYP5gB234ioFuVgUkWsLV5dNOmqju6LqGeixPuUAjsCwBdJl927z8UI0UZ3BBvNx6PZ2F+y2H7Kj2qjFFU3mB3MimnmT3erHW+vV0fNHrU83gGMwbIF0isXMBwzReHVFdpkPQ7fXmtPqqPvblpbXmx/KiWh2eWtyJxy2Ozp3WJ9qAEdh2ALpNDgoLVliPmKI3tjp6DbzYZgJteY642JSa1fVmR/IiWh2+dc3TThqy5rLrE80gOMwbIF0O3NGCgbNhwzR5V2IlJmPwkypJdSiysJKu2G7fIsCFdyvlsjthav2jjtqC/YVaHB40Po0AzgOwxaw0NpqPmSILm/Qv8B8EGZSjQsabUZtUaUWVDSbH8iJaPZFt1aPe7GoC/0XrE8xgCMxbAErO3eajxmiy+PKyMmtPr8+vaO2cJOWcLEooowpUr121Kj1b/erq6fL+vQCOBbDFrBUVWU+ZogudTy813wMZlo7i3ambdgWVTSYH8SJKHmFaotHXQG59Uyr9akFcDSGLWApHpfKy80HDZFycnQmstl8CGZc3kOqLqpO+agtie0zP4QTUXLz1gRGhu2uY7usTyyA4zFsAWtDQ9Ly5eajhqg3vNJ+CGZgrf7UXil51Zo95gdwIkpuwY2NI6N2w6EN1icVwBUYtoAT9PdL+fnmw4ayu6HciPkIzNRagi3aXLQ56aO2vGyn+QGciJJfZMt25dTmaOWBlRqOD1ufUgBXYNgCTnHhgjR/vvm4oezuiL/ZfARmai2hFlUVViVt1K5dVSevAw7gRJT8ItVrtWz/Mm7rA0wDwxZwktOnucctmXYisst8AGZyzZHmpNzjdt2KauXGWs0P30SUmor2xnRx8KL1qQRwFYYt4DRHj0q5ueYDh7Kzc9FK8/GX6TXNb9Kmwk0zHrXrl29VkFFLlLGF1x9RT3+P9WkEcB2GLeBEBw9KXq/5yKHsqy+43Hz4ZUONCxq1sWDj9IdtcZUiFS3mB28iSk3+WLtWx3j5MTATDFvAqZqaGLeU9oa9AfPRly01LGqY3qgtqtSCimbzgzcRpabcWLsC0QHV1VkfQAB3YtgCTsa4JYM6AgfMR1+2tG/JvqmN2sJNWlLRZH7wJqLU5IsdVjCvXzk5UkeH9eEDcCeGLeB0Bw4wbimtdUV2mA++bGrf0iuM26JKLWbUEmVsubH2kVHr80mDvBIZmBGGLeAGjY2MW0pb3ZGY+djLtvblTzBuiyq1iFFLlLHlxtoVmN8/8hBcVGR94ADci2ELuEVTU+I/5Tpg+FBmN+QLmw+9bGzMy5KX8Z5aokzO/2/vqb38IXjrVuvDBuBeDFvATQ4eZNxSWuoM7TcfetlYw+KGxNWSi6uUx6glytj8sSPyRwbHPPy2t1sfNAD3YtgCbtPWxn1uKeWdjVaZj7xsbe/iJuWt5ZY+RJmaP3ZE/vDYURsISMPD1ocMwL0YtoAbdXRIfr/5+KHMrT+41HzgZWMHAh3K8w8qvKRX3lib+QGciJJbINah3NDQuA+95eXWhwvA3Ri2gFudOCFFIuYDiDK3dn+z+dDLphqCRxX2vX7gDS68KF/ssPlBnIiSU6i8U77A+KM2JydxEwQAM8ewBdysu1tassR8AFFmdjJSZz72sqV9oWMKeIfH/DUE5vfLH2s3P5AT0ewKrzkhry8+4UOu1ytdvGh9qADcjWELuF1/v7RihfkIosyrJ7zKfPBlQ7siJ5TrnfjAmxsaUjB21PxgTkQzK7rytHJyJv4az8mRli+3PkwA7sewBTLB8LC0fr35EKLMatjrV5u31Xz4ZXJbomem9Nfh9Q0rvOaE+QGdiKZXpPjclL7Gd+2yPkgA7sewBTKqoC8eAAAfvElEQVTJ9u3mY4gyq2PhevPxl4m1ettUFrkwzb+OeOKZHwcc1onoSrUpXNAz5a/vc+esDxCA+zFsgUzT1MTtgChpnYtsMB+BmVazv10Fwf4Z/7VEis8rhysmEzm23Fi7QosuTvlrevFi64MDkBkYtkAm6uyUwmHzUUTub9C/wHwIZlINwU5Fcye+KupUCy/p5YrJRA4sGDsqf2TsPWona9s260MDkBkYtkCmOntWWrTIfBiR++sIHDAfhJnQzkjXpBeJmm6JKyYfMT/IE1GiSFmXvL6xVze/UseOWR8YgMzAsAUy2cCAFIuZDyNyd6ej28xHodvbHD2bkr8en5+LShHZ16bo8qldJOqNRSJSPG59WAAyA8MWyAb19bzvlmZcfyDffBi6tYO+wyoLT/0CMjMtWnJOvO+WKP0l3k/bN+Ov3fXrrQ8IQOZg2ALZoqtLWrjQfCSROzsabDQfiW6rIXhUC/0DaftrCi3qU26s3fygT5QtBcs7lRua3vtp39jhw9aHAyBzMGyBbNLfL61ZYz6SyH2dj6w3H4puqjp6Wr6c5L2fdqrlBocUKj9mfuAnyvQiq04pZ5bvmc/L42XIQDIxbIFstGeP5POZjyVyT0O+kNq8reaD0em15LarNNxr/NfF/W6JUlVurF3hpcl5e0FNjfVhAMgsDFsgWx0/nvjPxQ4YTeSOuiK7zIejk6sPH0vKrXySVXhJLy9NJkpi4TXHlBtKzte41yt1d1sfBIDMwrAFstnFi9Lq1fYncHJFveEV5uPRmbWpMkVXPZ5tvsCQImVcNZlods38qscTVVpqfQAAMg/DFoDU2CiFQvancHJ08RyvjvibHTAknVOz/4iWhS5a/9VcsUjhBflihx0wEIjcVSDWoWBef9K/JltarL/xA5mHYQsgoadHKiuzP4GTozsT3WI+Jp3S9shJBX3D1n8lUy43NKgwF5YimnKRVafkTcHXeDgsDQ9bf9MHMg/DFsBozc2J77rWp3ByZIP+BeaD0rom/xGVhGZ+30rrIsXn5OWet0QTFoh1KLQ4dV/jmzdbf6MHMhPDFsBYfX1Sebn9CZwc2bHQPvNxaVObtkVPK3eWt/hwQoHoALcFIhpTm6IrT8vrS+3X+OnT1t/kgczEsAUwsYMHpUjE/hROjqo7EnPAyExvjcGjKggm/3121kWKurlyMlHskIKxoyl5L+0bKyqy/sYOZC6GLYDJXbworV1rfwInxzTs9euw76D52ExHrd42bYqckzfH/c/STpQvMKTw6i7zYUFkkTfWpmjJOeWk6Wt8/37rb+pA5mLYApiaw4elxYvtT+HkiE5G6sxHZ6qrDx/TQv+A9R912got7lMg1mE+NIjSVXjNcfkj6fsaDwSkgQHrb+ZA5mLYApi64WFpzx5uDUTqDxaYD89U1exvVyzSbf1HbJM3ruiKs1xcijK6QKxD4fzetH99rV9v/U0cyGwMWwDTd/GiVFUl+XwOOImTVcdCe81HaDJr9bZpS/SM/F733MInVfnDg4qU8fJkyqx8scOJlx0bXQDu1Cnrb95AZmPYApi5M2ekVavsT+FkUl+o2HyMJqudkS7l+Qet/0gdV3DhRa6eTBlRZNVJ+QJDZl9LZWXW37CBzMewBTB77e3SkiX2p3BKe52h/eajdDbtCx1TUfCi9R+j4wsv7eH9t+TKQuXH0nK14yt1/Lj1N2og8zFsASRHPC7V10vhsP0pnNJWb7jUfJzOpAOBDq0M9Vj/8bkrb1zR5ee4PRC5omD5UZP30Y7XihXW36CB7MCwBZBc/f1SXR0DN4s6Gmw0H6pTrdl/RGsj3Rl9+55U5/MPK1p6Rr7YYfPxQvTGgrEOhZc66z9adXRYf2MGsgPDFkBqDAxIO3dKkYj9qYJS2oVImflgvVJNgQ6tjXTLx6BNWj7/sKIrzjJwyREFYh2KFFww/7p4Y8uWWX8zBrIHwxZAag0OJm4RNH++/QmDUlI8x6uOwAHz8TpejcEOrYlc4BnaFObNTQxcXqJMFvljRxQp6laOQ7/G29qsvwkD2YNhCyA9hoelffukBQvsTxqU9LrDFeYj9vIagp1aFXbWyxEzPW/usKIlDFxKT8HYUUUKL8ipgzYnR1q61PobL5BdGLYA0mt4WGpslBYtsj91UNKK5/h0xN9sPmj3hY6pJNRn/ceR1Xl9w4oUn+cqypSSQuXHFF7ijItCXanmZutvuEB2YdgCsBGPSwcPSiUl9qcPSkrnIxuMBm2bdkVOcNseBxbO71V4DffBpdnWpsjqLkfctmeqLVqU+DYHIH0YtgDsnT0rbd4shUL2pxGacXFvrtpzW9I2aJv9R1QVPato7pD1b52uUCA6oMiqk/LG2hwwksgt5cbaFS09o9yQ+77G9++3/sYKZB+GLQDnGBpKvEy5qMj+VEIz6my0MuWDdk/ouFaGe7kglAvzBYYULTknf+yI+Wgi5xZeczxxyx6vO7/G8/IS77oBkF4MWwDOdPKktHGjFAjYn1Joyg37gil51rYlt11bomeU5x+0/i1Skgot7lOkrItncUk5scTVjaOlZ+QPu/9rnGdrARsMWwDONjAg1ddLhYX2pxWaUsl7r22b6sPHVBbu4f6zGZwvMKRI8XkFyzvNxxWluzaF15xQOL9XTr668XRaupT31gJWGLYA3OPcOWnHDqmgwP70QhMWz/HN6r629eFjWhs+r7DPfe+ro9kViA4knrXjpcoZXbj8mCLLuuULZN7X+NGj1t8ogezFsAXgToxcR9cTWjmtMbsvdEzrIoxZer3ggn5FV57mtkEZUqj8mCLF55UbzNyv8TVrrL8xAtmNYQvA/Ri5juxYuH7yZ2aDx7Q+cl4RrmpMVygwv1/RFWcVZOS6qlB5p6LLz7nyqsbTLTdXOn/e+pshkN0YtgAyy6WRW1Ii+Xz2p50srj+4dNSQPeg7rN3hE1oXOc8temjG+SMDipacVXDNMeVw4SlHlRtrV2R1lyKFFzLyZcaTtW2b9Tc/AAxbAJlrYEBqa5MqK6VFi+xPPtmW16vOpQ3aEjmj4lAfF4CipOf1DSuc36voytM8m2tUqLwz8Wz6gn7zfw9WRaOJbzcAbDFsAWSP7u7EfXLXr5cWLLA/DWVi+flSVZXU2ir196u3VwoG7X9ZlB3lhoYUKepWpKxLubF289GXebUpWN6p6MrTChf0ZN2zshPV0GD9zQ2AxLAFkM0uDd1Nm6SiIsnvtz8hualwWCotlaqrpYMHpd7ecf+Y9+yx/6VSduYPDypScCHxjG55J/fMnWa+2GGF1xxPPCO7qE9e37D536nTKihI8/ctABNi2ALAJfG4dOpUYuxWVUnLl/N046XeOGK7u6f1x5qfb/9bIMrxxhVYcFGR4vOKlHX92xWXGbs5scT7Y0PlxxRZeUqRZd0K5PUrU+4tm8o6O1P3LQnA9DBsAeBKzp6VmpsTVwcpLZUWLszcC1MFg4mnINasmdGInUhnp/1vjWjcvHEF5vcnntktPaNI2Qn5yzN38HpjbQrGOhQp61K05KzC+b3yBQft/x5cWCw2+28vAJKHYQsAMxGPJwZfZ2fiGd7aWmndusSzvHl5ktdrf+oaL6838f7ikhJpw4bEFaSbm6UTJ6SLF1P6R7Z2rf1vn2jK/dvgDS/tUaT4vKIrTytS1qVQeaf8Dn//7sizr6u7FF1xVpGibgUX9mXFbXfSld+flP/mByCJGLYAkArDw4lbDx09mriQUmNj4s2mdXXS5s2JC1iVlSUGZkFBYmyGQonTkt8vBQKvFwyOLRRKlJcnLV2a+Dxr1iTG6taticG6b5/U0iK1tyeG6/nziV+XkZ4eXtlNGZQ3Ln94UKFFfYnxu6xb0ZJziWd9V51SeHWXwmuOKVjeqUCsQ/7YEeXG2uWLtcsXOyxf7LC8sTa9NmaYtskba5M3dli5sXb5Y0cUjHUoWN6ZGKtlJxRZdUqRFWcSL6kuvKDQ4j4F8voTw9XLy4fT0a5dZg+lACbAsAUApE1jo/2BlMiRMUhdU3Fx4kU7AJyFYQsASKuyMvuDKRHRTPL7Ey/GAeA8DFsAQFr19iYusmx9QCUimm719daPoAAmwrAFAKRda6v9AZWIaDqtWGH9yAlgMgxbAIAJrpJMRG4pGOQqyIDTMWwBACb6+xMXdbY+sBIRXamGButHTABXwrAFAJg5csT+wEpENFmrVlk/UgKYCoYtAMBUZaX9wZWIaLxCocQ9uAE4H8MWAGBqcFBatMj+AEtE9MaamqwfIQFMFcMWAGDu+HHJ67U/xBIRXWrNGutHRgDTwbAFADhCdbX9QZaIKCdHWrBA6uuzflQEMB0MWwCAIwwPSyUl9gdaIsru/H6pq8v6ERHAdDFsAQCO0dfHLYCIyLYDB6wfCQHMBMMWAOAoXV1Sbq794ZaIsq+qKutHQAAzxbAFADjOgQP2B1wiyq6KixNviQDgTgxbAIAjVVXZH3SJKDuKRqXeXutHPQCzwbAFADgSF5MionTk8yVuOQbA3Ri2AADH4mJSRJTq9u2zfqQDkAwMWwCAo3V1JW6/YX34JaLMa8MG60c4AMnCsAUAOF5Tk/0BmIgyq8JCaWjI+tENQLIwbAEArrB5s/1BmIgyo0hE6u62flQDkEwMWwCAK8TjUmmp/YGYiNxdMJh4iwOAzMKwBQC4xuCgVFRkfzAmIneWmysdPWr9SAYgFRi2AABXuXhRys+3PyATkbvyeqVDh6wfwQCkCsMWAOA6PT3SwoX2B2Uick+NjdaPXABSiWELAHCl8+elaNT+sExEzm/XLutHLACpxrAFALjW6dNSKGR/aCYi57Ztm/UjFYB0YNgCAFzt+HEpELA/PBOR89qwwfoRCkC6MGwBAK7X0ZG42qn1IZqInNOaNYnbhAHIDgxbAEBGaG1NXPXU+jBNRPatWCEND1s/KgFIJ4YtACBjNDbaH6iJyLaiImlgwPrRCEC6MWwBABmlvt7+YE1ENhUXS/391o9CACwwbAEAGaepSfL57A/ZRJS+VqyQBgetH30AWGHYAgAyUlub5PfbH7aJKPWtWiUNDVk/6gCwxLAFAGSszk4pGLQ/dBNR6iov50JRABi2AIAMd/KkFInYH76JKPmtW8ctfQAkMGwBABnv3DlpwQL7QzjZ9/nP/1/deeeHNXfuDZo37/d0332P6ec/bxj5+eeea5XH4xm3b31rycjHPfTQ9/Tud39IV189R7fffp/57ysb27LF+pEFgJMwbAEAWaGnR1qyxP4wTrbdc8+n9cQTQT399F499dROfeADf6G3ve3d+s1vLignR3r11SG9+GLnqD73uZ9r7twb9PLL3SOf56GHvqcvf/l3+vf//r8xbA3atcv6EQWA0zBsAQBZ4+LFxD0urQ/l5Jx++csT8ng8+uEPN0z4MXfc8UF97GPfHPfnPvvZpxm2acznkw4csH4kAeBEDFsAQFYZHJRKS+0P6OSMnn22SR6PR//0T3vG/fl/+IdaeTwe/f3fV4378wzb9OX3J652DgDjYdgCALLO8LBUUWF/UCfbXn11WH/yJ5/RH/zBxyb8mAce+I5uu+2PJ/x5hm16mj9fOnHC+pEDgJMxbAEAWWvnTsnrtT+0k0333/9t3XLLnXrhhfZxf/63v+3VtdfeqC984ZcTfg6Gbepbvlzq7bV+tADgdAxbAEBWO3JECoftD++U3h588EndfPPteu65gxN+zDe+EdGb33yNfvnLExN+DMM2tVVWco9aAFPDsAUAZL3z56WCAvtDPKW+116L68EHn9RNN/2+nnnmwKQf+/73P6APfejxST+GYZuacnOlxkbrRwYAbsKwBQBAiYtKrV1rf6Cn1PbAA9/RtdfeqB/+cP2oW/r89re9oz7u2WebdNVVV+l731s17ud59tkmPfXUDt1//9/pHe94v556aoeeemqHXnml3/z36PZ4Py2AmWDYAgBwmT17ErcUsT7cU2ryeDzj9sQTwVEf98gjP9HNN9+hV18dHvfzvP/9D4z7eZ57rtX89+jmiot5Py2AmWHYAgDwBkePSpGI/SGfKJvatIn30wKYOYYtAADjuHBBKiqyP+wTZXq5uVJDg/VXPAC3Y9gCADCBoSFpwwb7gz9RppaXJx0/bv2VDiATMGwBALiCpiYpFLIfAUSZ1Pr1Un+/9Vc3gEzBsAUAYAp6eqRVq+zHAJHbi0SkQ4esv6IBZBqGLQAA07B/vxQI2I8DIje2Zo3U12f9VQwgEzFsAQCYpu5uqaTEfiQQuaVgUDpwwPorF0AmY9gCADBDe/ZIfr/9aCBycqWliauMA0AqMWwBAJiFs2elZcvsxwOR0/L7pfp6669QANmCYQsAwCzF49KOHYn7cVqPCSIntGyZdO6c9VcmgGzCsAUAIElOnZIKC+1HBZFVwaC0e3fiP/YAQDoxbAEASLLGRikatR8ZROnK65U2bOCKxwDsMGwBAEiBgQFp2zZenkyZX3Gx1NVl/RUHINsxbAEASKFz5xL37rQeH0TJLi9Pam62/goDgASGLQAAadDRIeXn248Rotnm90s1NdLgoPVXFQC8jmELAECaxOPSvn1SOGw/TohmUnm51N1t/ZUEAGMxbAEASLP+fmnzZsnnsx8qRFMpPz/xqgMAcCqGLQAARs6c4f235OwWL5YOHOD2PQCcj2ELAICxU6cSL/H0eu2HDFFOjrR0qdTSYv2VAQBTx7AFAMAhTp+WKioYuGRXUZHU1mb9lQAA08ewBQDAYc6ckdat4z24lL6WL5eOHLH+lw8AM8ewBQDAobq7ExeZ8vvthw9lZqWlUmen9b90AJg9hi0AAA538aJUV8dtgih5rV4tnThh/S8bAJKHYQsAgEsMDibug1tQYD+MyH2Fw9LWrdLZs9b/kgEg+Ri2AAC40MmTUmWlFAzaDyZydqWliSscDw9b/6sFgNRh2AIA4GJDQ1JTk1RSYj+gyDlFo1J1tXTunPW/UABID4YtAAAZ4tw5adu2xKixHlaU/rxeadUqqbWVZ2cBZB+GLQAAGSYelw4dSlwgiFsGZX55eVJtbeIq2gCQrRi2AABksJ4eaffuxH1KvV77EUbJaeFCacuWxK164nHrf2UAYI9hCwBAlrh4UTpwQCovlwIB+3FG06uwUNq+XTp92vpfEgA4D8MWAIAsNDwstbdLmzYlXspqPdpobD5f4orG9fXShQvW/2IAwNkYtgAAQF1difdpco9c2wKBxDPqTU1Sf7/1vwoAcA+GLQAAGOXCBamhQVq/Xlq82H7sZXLhsFRWJu3aJR0/ztWMAWCmGLYAAGBSFy9KbW2J+6KWlPD+3Nm0eLG0YUPiPxycPWv9NwsAmYNhCwAApiUeT7x0ee9eqaJCWrDAfjA6sdxcqbhY2ro1cfulvj7rvzkAyFwMWwAAMGs9PYnxtnNn4hnJ5culaNR+XKYjrzdx+53Vq6Vt2xJXnu7qkoaGrP9WACB7MGwBAEDKDAwkRl5TU+LiVBUVidvWBIP2g3S6RaOJZ2DXrUvcdqe5OfF7Gxy0/lMGADBsAQCAid5eqbNTammR9u1LjMUtWxLDcfXqxIhcvDhxgSWfL/lD1e9P3OqooCBxW52KCqmqKjHA9+5NDNcjR6RTpxIDHQDgXAxbAADgCv390rlz0okTUkdHYnR2dIzf0aOv19n5eidPJq76zLOsAJBZGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAA0sLj8ejpp5+2/mUkVWtrqzwej1566SXrXwoAZDWGLQAAmLFgMCiPxzOq3/u939ODDz6olStXjvrYmQ7bdevWjXzu2traMT//xBNP6Prrr5/w1zNed95550x+u2MwbAHAGRi2AABgxi4NyWeeeUbRaFSRSEQvvfSS7r33Xnk8HpWUlIx8bF9fnwYHB6f9/+PyYfvZz352zM9fPmxbWloUjUZHNXfuXH384x8f9WNFRUUz/01fhmELAM7AsAUAADN2adjW1NSM+vHTp0/rmmuu0Ve/+tVZ//+4NGw/+MEPyuPxqK6ubtTPXz5sx3P99dfriSeemPWvYzwMWwBwBoYtAACYsYmGbTwe11vf+lb99V//9ciPjfdS5O3bt+uRRx7RvHnzdP311+sTn/iEtmzZMupjLg1bv9+vm2++WZ/73OdG/fxMhu2RI0f0jW98Q7feeqvmzJmje+65R36/f8z/tq+vT08//bTe9773ae7cuXrnO9+pv/zLv1Rzc7Ok0cM2JydHd999t+bMmaMPf/jDqq6unvDXBABILoYtAACYsUvDNhaLqaurSydOnNDevXv1d3/3d3rTm96kNWvWjHzsG4ft3r17df311+u2227Ts88+qxdeeEF33XWX5s6dq61bt4583KVhm5+fr2eeeWbMs7bTHbbHjh3T7bffrjvuuEPPPPOMXn31VT366KPyeDz6l3/5l5GPGxoa0sMPPyyPx6Mvf/nL+t3vfqfnn39en/jEJ7Rs2TJJrw/bP/uzP9N73/te/eIXv9CLL76ot7/97br99ts1MDAwmz9eAMAUMWwBAMCMTXSxprlz5yoUCo362DcO289//vOaM2eOWlpaRn7s6NGjmjdvnu6///6RH7t82J49e1Y333yzHn300ZGfn+6w/Zu/+RvddtttOnny5KiP+/KXv6wbb7xRvb29kqRAICCPx6Nf//rXYz5nPB6X9PqwveWWW3T69OmRny8uLh7zHmMAQOowbAEAwIxdGravvPKKysvLVV5ervnz5+uRRx7R1VdfrYKCgpGPvXzYDg0N6brrrtOXvvSlMZ/z0rO9586dkzR62ErSz3/+c3k8Hm3fvl3S9IZtPB7XTTfdpG9961vq6uoa1aXfS2VlpSTpM5/5jN7+9rdPesGrS8P2u9/97qgfP336tDwej15++eUr/AkCAJKBYQsAAGZsovfYDg8P60//9E912223qb+/X9LoYdvZ2SmPx6N//Md/HPM5//Vf/1Uej0d79+6VNHbYnj17VjfddNPIs7bTGbbHjx+/4q2ACgsLJUl/9Ed/pI997GOT/v4vDdsXXnhhzM95PB797Gc/m/R/DwBIDoYtAACYsYmGrSR9//vfHzVQkzVsJelnP/vZyLO20xm2l/7//tVf/dXIM8xv7Pjx45KmN2zHuyryTO/bCwCYPoYtAACYscmG7ZNPPimPxzNyIaipvhT529/+9qQvRZZef9b2sccem9awHRoa0rx58/SVr3zlir+3Sy9FnuwCUAxbAHAGhi0AAJixiYbtwMCA3ve+92nOnDkjA3W8i0fNnTtXra2tIz927NgxvfWtb53w4lGXu/Ss7X333Teti0d9/etf15w5c7Rnz54xH3vixImR/3s6F49i2AKALYYtAACYsUvD9plnnlE0GlU0GtWvfvUr/fmf/7k8Ho9+/OMfj3zsRLf7ede73qXnnntOv/jFL3T33XdPerufy505c0Y33nijPB7PtG/3c+edd+q6667T97//feXk5Oj555/XF7/4Rd18880jHzc0NKQHH3xw5HY/r7zyil588UV96lOfGnO7H4YtANhi2AIAgBkb73Y/b3nLW/TBD35Qr7766sgzm9L4Q2/79u369Kc/rRtuuEHXXXedHnroIW3evHnUx0w0bCXp6aefnvawlRIXkXryySd1xx136JprrtE73/lOPfzww/J6vaM+rre3Vz/96U911113jXzcF77whZFbFDFsAcAZGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXI1hCwAAAABwNYYtAAAAAMDVGLYAAAAAAFdj2AIAAAAAXO3/A/i5yWGrtaMhAAAAAElFTkSuQmCC\" width=\"950\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib_venn._common.VennDiagram at 0x7fa13ab9feb0>"
      ]
     },
     "execution_count": 50,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plt.figure()\n",
    "venn3([set(data.indorig), set(data.indcvac), set(data.indbntc)], set_labels = ('SARS-CoV-2', 'CureVac', 'BioNTech'))\n"
   ]
  }
 ],
 "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}