tb_profiler_data_visualization.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 151,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 152,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create URL to JSON file \n",
    "url = 'https://fe4qpw.by.files.1drv.com/y4mxQxTWOjdBgxzDzcQS5G50gBvRmeocJMV1rKibOTYGFmHKcEJ_UtY59CoqV2nHifh5RaNqhZcET6Jan1RTmOi6jNbCJ7ghFQrFJgpKbsKu-pE95l8YDU5DKdQvdKQk1UVTlyEtQPSkG7bXQXsJxJw25zVqUmscf0-7e8y0cCmN1ktbRJ_RBCOIz4CtuAidGui6mJEXfDgrK0wEqoxzXnLcQ'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 153,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the JSON file into a data frame\n",
    "df = pd.read_json(url, orient='index')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 154,
   "metadata": {
    "scrolled": true
   },
   "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>rifampicin</th>\n",
       "      <th>isoniazid</th>\n",
       "      <th>pyrazinamide</th>\n",
       "      <th>ethambutol</th>\n",
       "      <th>streptomycin</th>\n",
       "      <th>fluoroquinolones</th>\n",
       "      <th>moxifloxacin</th>\n",
       "      <th>ofloxacin</th>\n",
       "      <th>levofloxacin</th>\n",
       "      <th>ciprofloxacin</th>\n",
       "      <th>...</th>\n",
       "      <th>cycloserine</th>\n",
       "      <th>linezolid</th>\n",
       "      <th>bedaquiline</th>\n",
       "      <th>clofazimine</th>\n",
       "      <th>delamanid</th>\n",
       "      <th>main_lin</th>\n",
       "      <th>sublin</th>\n",
       "      <th>drtype</th>\n",
       "      <th>MDR</th>\n",
       "      <th>XDR</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>ERR760889</th>\n",
       "      <td>rpoB_p.Ser450Leu</td>\n",
       "      <td>katG_p.Ser315Thr</td>\n",
       "      <td>pncA_p.Gln10Pro</td>\n",
       "      <td>embB_p.Gly406Ala</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>...</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>lineage4</td>\n",
       "      <td>lineage4.1.2.1</td>\n",
       "      <td>MDR</td>\n",
       "      <td>R</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR10851621</th>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>...</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>lineage4</td>\n",
       "      <td>lineage4.1.2</td>\n",
       "      <td>Sensitive</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR751514</th>\n",
       "      <td>-</td>\n",
       "      <td>katG_p.Ser315Thr</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>rrs_r.517c&gt;t</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>...</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>lineage4</td>\n",
       "      <td>lineage4.1</td>\n",
       "      <td>Drug-resistant</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR760825</th>\n",
       "      <td>rpoB_p.Ser450Leu</td>\n",
       "      <td>katG_p.Ser315Thr</td>\n",
       "      <td>pncA_p.Gln10Pro</td>\n",
       "      <td>embB_p.Gly406Ala</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>...</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>lineage4</td>\n",
       "      <td>lineage4.1.2.1</td>\n",
       "      <td>MDR</td>\n",
       "      <td>R</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR1765879</th>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>gyrA_p.Asp94Asn</td>\n",
       "      <td>gyrA_p.Asp94Asn</td>\n",
       "      <td>gyrA_p.Asp94Asn</td>\n",
       "      <td>gyrA_p.Asp94Asn</td>\n",
       "      <td>gyrA_p.Asp94Asn</td>\n",
       "      <td>...</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>-</td>\n",
       "      <td>lineage2</td>\n",
       "      <td>lineage2.2.1</td>\n",
       "      <td>Drug-resistant</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 26 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                   rifampicin         isoniazid     pyrazinamide  \\\n",
       "ERR760889    rpoB_p.Ser450Leu  katG_p.Ser315Thr  pncA_p.Gln10Pro   \n",
       "SRR10851621                 -                 -                -   \n",
       "ERR751514                   -  katG_p.Ser315Thr                -   \n",
       "ERR760825    rpoB_p.Ser450Leu  katG_p.Ser315Thr  pncA_p.Gln10Pro   \n",
       "SRR1765879                  -                 -                -   \n",
       "\n",
       "                   ethambutol  streptomycin fluoroquinolones     moxifloxacin  \\\n",
       "ERR760889    embB_p.Gly406Ala             -                -                -   \n",
       "SRR10851621                 -             -                -                -   \n",
       "ERR751514                   -  rrs_r.517c>t                -                -   \n",
       "ERR760825    embB_p.Gly406Ala             -                -                -   \n",
       "SRR1765879                  -             -  gyrA_p.Asp94Asn  gyrA_p.Asp94Asn   \n",
       "\n",
       "                   ofloxacin     levofloxacin    ciprofloxacin  ...  \\\n",
       "ERR760889                  -                -                -  ...   \n",
       "SRR10851621                -                -                -  ...   \n",
       "ERR751514                  -                -                -  ...   \n",
       "ERR760825                  -                -                -  ...   \n",
       "SRR1765879   gyrA_p.Asp94Asn  gyrA_p.Asp94Asn  gyrA_p.Asp94Asn  ...   \n",
       "\n",
       "            cycloserine linezolid bedaquiline clofazimine delamanid  main_lin  \\\n",
       "ERR760889             -         -           -           -         -  lineage4   \n",
       "SRR10851621           -         -           -           -         -  lineage4   \n",
       "ERR751514             -         -           -           -         -  lineage4   \n",
       "ERR760825             -         -           -           -         -  lineage4   \n",
       "SRR1765879            -         -           -           -         -  lineage2   \n",
       "\n",
       "                     sublin          drtype MDR XDR  \n",
       "ERR760889    lineage4.1.2.1             MDR   R      \n",
       "SRR10851621    lineage4.1.2       Sensitive          \n",
       "ERR751514        lineage4.1  Drug-resistant          \n",
       "ERR760825    lineage4.1.2.1             MDR   R      \n",
       "SRR1765879     lineage2.2.1  Drug-resistant          \n",
       "\n",
       "[5 rows x 26 columns]"
      ]
     },
     "execution_count": 154,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# A Glimpse of the dataset.\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 155,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "object    26\n",
       "dtype: int64"
      ]
     },
     "execution_count": 155,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# What's the count of different types of objects.\n",
    "\n",
    "df.dtypes.value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 156,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "Index: 735 entries, ERR760889 to ERR775312\n",
      "Data columns (total 26 columns):\n",
      " #   Column                    Non-Null Count  Dtype \n",
      "---  ------                    --------------  ----- \n",
      " 0   rifampicin                735 non-null    object\n",
      " 1   isoniazid                 735 non-null    object\n",
      " 2   pyrazinamide              735 non-null    object\n",
      " 3   ethambutol                735 non-null    object\n",
      " 4   streptomycin              735 non-null    object\n",
      " 5   fluoroquinolones          735 non-null    object\n",
      " 6   moxifloxacin              735 non-null    object\n",
      " 7   ofloxacin                 735 non-null    object\n",
      " 8   levofloxacin              735 non-null    object\n",
      " 9   ciprofloxacin             735 non-null    object\n",
      " 10  aminoglycosides           735 non-null    object\n",
      " 11  amikacin                  735 non-null    object\n",
      " 12  kanamycin                 735 non-null    object\n",
      " 13  capreomycin               735 non-null    object\n",
      " 14  ethionamide               735 non-null    object\n",
      " 15  para-aminosalicylic_acid  735 non-null    object\n",
      " 16  cycloserine               735 non-null    object\n",
      " 17  linezolid                 735 non-null    object\n",
      " 18  bedaquiline               735 non-null    object\n",
      " 19  clofazimine               735 non-null    object\n",
      " 20  delamanid                 735 non-null    object\n",
      " 21  main_lin                  735 non-null    object\n",
      " 22  sublin                    735 non-null    object\n",
      " 23  drtype                    735 non-null    object\n",
      " 24  MDR                       735 non-null    object\n",
      " 25  XDR                       735 non-null    object\n",
      "dtypes: object(26)\n",
      "memory usage: 155.0+ KB\n"
     ]
    }
   ],
   "source": [
    "# Let's see information about the variables. \n",
    "df.info()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 157,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(735, 26)"
      ]
     },
     "execution_count": 157,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Taking Description of the Datasets\n",
    "df.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 158,
   "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>count</th>\n",
       "      <th>unique</th>\n",
       "      <th>top</th>\n",
       "      <th>freq</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>rifampicin</th>\n",
       "      <td>735</td>\n",
       "      <td>34</td>\n",
       "      <td>rpoB_p.Ser450Leu</td>\n",
       "      <td>333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>isoniazid</th>\n",
       "      <td>735</td>\n",
       "      <td>32</td>\n",
       "      <td>katG_p.Ser315Thr</td>\n",
       "      <td>414</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>pyrazinamide</th>\n",
       "      <td>735</td>\n",
       "      <td>67</td>\n",
       "      <td>-</td>\n",
       "      <td>400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ethambutol</th>\n",
       "      <td>735</td>\n",
       "      <td>29</td>\n",
       "      <td>-</td>\n",
       "      <td>409</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>streptomycin</th>\n",
       "      <td>735</td>\n",
       "      <td>18</td>\n",
       "      <td>-</td>\n",
       "      <td>650</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>fluoroquinolones</th>\n",
       "      <td>735</td>\n",
       "      <td>21</td>\n",
       "      <td>-</td>\n",
       "      <td>693</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>moxifloxacin</th>\n",
       "      <td>735</td>\n",
       "      <td>21</td>\n",
       "      <td>-</td>\n",
       "      <td>693</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ofloxacin</th>\n",
       "      <td>735</td>\n",
       "      <td>21</td>\n",
       "      <td>-</td>\n",
       "      <td>693</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>levofloxacin</th>\n",
       "      <td>735</td>\n",
       "      <td>21</td>\n",
       "      <td>-</td>\n",
       "      <td>693</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ciprofloxacin</th>\n",
       "      <td>735</td>\n",
       "      <td>21</td>\n",
       "      <td>-</td>\n",
       "      <td>693</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>aminoglycosides</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td>-</td>\n",
       "      <td>572</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>amikacin</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td>-</td>\n",
       "      <td>572</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>kanamycin</th>\n",
       "      <td>735</td>\n",
       "      <td>6</td>\n",
       "      <td>-</td>\n",
       "      <td>564</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>capreomycin</th>\n",
       "      <td>735</td>\n",
       "      <td>5</td>\n",
       "      <td>-</td>\n",
       "      <td>568</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ethionamide</th>\n",
       "      <td>735</td>\n",
       "      <td>33</td>\n",
       "      <td>-</td>\n",
       "      <td>627</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>para-aminosalicylic_acid</th>\n",
       "      <td>735</td>\n",
       "      <td>10</td>\n",
       "      <td>-</td>\n",
       "      <td>716</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>cycloserine</th>\n",
       "      <td>735</td>\n",
       "      <td>1</td>\n",
       "      <td>-</td>\n",
       "      <td>735</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>linezolid</th>\n",
       "      <td>735</td>\n",
       "      <td>1</td>\n",
       "      <td>-</td>\n",
       "      <td>735</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>bedaquiline</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td>-</td>\n",
       "      <td>725</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>clofazimine</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td>-</td>\n",
       "      <td>725</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>delamanid</th>\n",
       "      <td>735</td>\n",
       "      <td>1</td>\n",
       "      <td>-</td>\n",
       "      <td>735</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>main_lin</th>\n",
       "      <td>735</td>\n",
       "      <td>7</td>\n",
       "      <td>lineage4</td>\n",
       "      <td>698</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>sublin</th>\n",
       "      <td>735</td>\n",
       "      <td>32</td>\n",
       "      <td>lineage4.1.2.1</td>\n",
       "      <td>248</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>drtype</th>\n",
       "      <td>735</td>\n",
       "      <td>4</td>\n",
       "      <td>MDR</td>\n",
       "      <td>429</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MDR</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td>R</td>\n",
       "      <td>449</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>XDR</th>\n",
       "      <td>735</td>\n",
       "      <td>2</td>\n",
       "      <td></td>\n",
       "      <td>715</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                         count unique               top freq\n",
       "rifampicin                 735     34  rpoB_p.Ser450Leu  333\n",
       "isoniazid                  735     32  katG_p.Ser315Thr  414\n",
       "pyrazinamide               735     67                 -  400\n",
       "ethambutol                 735     29                 -  409\n",
       "streptomycin               735     18                 -  650\n",
       "fluoroquinolones           735     21                 -  693\n",
       "moxifloxacin               735     21                 -  693\n",
       "ofloxacin                  735     21                 -  693\n",
       "levofloxacin               735     21                 -  693\n",
       "ciprofloxacin              735     21                 -  693\n",
       "aminoglycosides            735      2                 -  572\n",
       "amikacin                   735      2                 -  572\n",
       "kanamycin                  735      6                 -  564\n",
       "capreomycin                735      5                 -  568\n",
       "ethionamide                735     33                 -  627\n",
       "para-aminosalicylic_acid   735     10                 -  716\n",
       "cycloserine                735      1                 -  735\n",
       "linezolid                  735      1                 -  735\n",
       "bedaquiline                735      2                 -  725\n",
       "clofazimine                735      2                 -  725\n",
       "delamanid                  735      1                 -  735\n",
       "main_lin                   735      7          lineage4  698\n",
       "sublin                     735     32    lineage4.1.2.1  248\n",
       "drtype                     735      4               MDR  429\n",
       "MDR                        735      2                 R  449\n",
       "XDR                        735      2                    715"
      ]
     },
     "execution_count": 158,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's see statistical information about the variables. \n",
    "df.describe().T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 159,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x265247f4708>"
      ]
     },
     "execution_count": 159,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD4CAYAAADLhBA1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAYEElEQVR4nO3df5BdZX3H8fdHwg9hIRt+uEQSXSxRsUQiWTGIdnaJVX5YEisZYVACk87SihRHUULrWKydGlSKAsp0BUpwIisiNGlEKBPZItUA2YAJASEBA8TEREyILiAM+O0f59ly2SzZm3uf3T17+bxm7txznvPcc57vvbn7uee5P6KIwMzMrF6vG+0BmJlZY3CgmJlZFg4UMzPLwoFiZmZZOFDMzCyLcaM9AIADDzwwWltba7rtM888wz777JN3QKPEtZRPo9QBrqWs6qmlt7f3qYg4KPOQalaKQGltbWXFihU13banp4f29va8AxolrqV8GqUOcC1lVU8tkh7PO5r6eMrLzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLIshA0XS2yTdX3H5vaRPS9pf0u2S1qbrCam/JF0maZ2kVZKOGv4yzMxstA0ZKBHxcERMi4hpwHTgWeBmYD6wLCKmAMvSOsAJwJR06QSuHI6Bm5lZuezqlNdM4NGIeByYBSxM7QuB2Wl5FnBdFJYDzZImZhmtmZmVlnbl/0ORdA2wMiKukPR0RDRXbNsWERMkLQUWRMRdqX0ZcEFErBiwr06KMxhaWlqmd3d311RAX18fTU1NNd22bFxL+TRKHeBayqqeWjo6Onojoi3zkGpW9TflJe0BnAxcOFTXQdp2SK2I6AK6ANra2qLWb4pevmgxl9z1TE23rdf6BSdl3Z+//Vs+jVIHuJayaqRadmXK6wSKs5PNaX1z/1RWut6S2jcAkytuNwnYWO9Azcys3HYlUE4Drq9YXwLMTctzgcUV7WekT3vNALZHxKa6R2pmZqVW1ZSXpL2BvwTOrmheANwgaR7wBDAntd8CnAiso/hE2FnZRmtmZqVVVaBExLPAAQPafkfxqa+BfQM4J8vozMxszPA35c3MLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyyqCpQJDVLulHSLyU9JOkYSftLul3S2nQ9IfWVpMskrZO0StJRw1uCmZmVQbVnKN8Ebo2ItwNHAg8B84FlETEFWJbWAU4ApqRLJ3Bl1hGbmVkpDRkokvYD/gK4GiAiXoiIp4FZwMLUbSEwOy3PAq6LwnKgWdLE7CM3M7NSUUTsvIM0DegCHqQ4O+kFzgN+HRHNFf22RcQESUuBBRFxV2pfBlwQESsG7LeT4gyGlpaW6d3d3TUVsGXrdjY/V9NN6zb1kPFZ99fX10dTU1PWfY6WRqmlUeoA11JW9dTS0dHRGxFtmYdUs3FV9jkKODci7pb0TV6e3hqMBmnbIbUioosiqGhra4v29vYqhrKjyxct5pLV1ZSR3/rT27Pur6enh1rvh7JplFoapQ5wLWXVSLVU8x7KBmBDRNyd1m+kCJjN/VNZ6XpLRf/JFbefBGzMM1wzMyurIQMlIn4DPCnpbalpJsX01xJgbmqbCyxOy0uAM9KnvWYA2yNiU95hm5lZ2VQ7V3QusEjSHsBjwFkUYXSDpHnAE8Cc1PcW4ERgHfBs6mtmZg2uqkCJiPuBwd74mTlI3wDOqXNcZmY2xvib8mZmloUDxczMsnCgmJlZFg4UMzPLwoFiZmZZOFDMzCwLB4qZmWXhQDEzsywcKGZmloUDxczMsnCgmJlZFg4UMzPLwoFiZmZZOFDMzCwLB4qZmWXhQDEzsywcKGZmloUDxczMsnCgmJlZFg4UMzPLoqpAkbRe0mpJ90takdr2l3S7pLXpekJql6TLJK2TtErSUcNZgJmZlcOunKF0RMS0iGhL6/OBZRExBViW1gFOAKakSydwZa7BmplZedUz5TULWJiWFwKzK9qvi8JyoFnSxDqOY2ZmY4AiYuhO0q+AbUAA/x4RXZKejojmij7bImKCpKXAgoi4K7UvAy6IiBUD9tlJcQZDS0vL9O7u7poK2LJ1O5ufq+mmdZt6yPis++vr66OpqSnrPkdLo9TSKHWAaymremrp6OjorZg1GnXjqux3bERslPQG4HZJv9xJXw3StkNqRUQX0AXQ1tYW7e3tVQ7llS5ftJhLVldbRl7rT2/Pur+enh5qvR/KplFqaZQ6wLWUVSPVUtWUV0RsTNdbgJuBo4HN/VNZ6XpL6r4BmFxx80nAxlwDNjOzchoyUCTtI2nf/mXgg8ADwBJgbuo2F1iclpcAZ6RPe80AtkfEpuwjNzOzUqlmrqgFuFlSf//vRcStku4FbpA0D3gCmJP63wKcCKwDngXOyj5qMzMrnSEDJSIeA44cpP13wMxB2gM4J8vozMxszPA35c3MLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyyqDpQJO0m6T5JS9P6oZLulrRW0vcl7ZHa90zr69L21uEZupmZlcmunKGcBzxUsX4xcGlETAG2AfNS+zxgW0QcBlya+pmZWYOrKlAkTQJOAq5K6wKOA25MXRYCs9PyrLRO2j4z9TczswamiBi6k3Qj8BVgX+B84ExgeToLQdJk4McRcYSkB4DjI2JD2vYo8J6IeGrAPjuBToCWlpbp3d3dNRWwZet2Nj9X003rNvWQ8Vn319fXR1NTU9Z9jpZGqaVR6gDXUlb11NLR0dEbEW2Zh1SzcUN1kPRhYEtE9Epq728epGtUse3lhoguoAugra0t2tvbB3apyuWLFnPJ6iHLGBbrT2/Pur+enh5qvR/KplFqaZQ6wLWUVSPVUs1f4mOBkyWdCOwF7Ad8A2iWNC4iXgQmARtT/w3AZGCDpHHAeGBr9pGbmVmpDPkeSkRcGBGTIqIVOBX4SUScDtwBnJK6zQUWp+UlaZ20/SdRzbyamZmNafV8D+UC4DOS1gEHAFen9quBA1L7Z4D59Q3RzMzGgl168yEieoCetPwYcPQgff4IzMkwNjMzG0P8TXkzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLIYMFEl7SbpH0i8krZH0pdR+qKS7Ja2V9H1Je6T2PdP6urS9dXhLMDOzMqjmDOV54LiIOBKYBhwvaQZwMXBpREwBtgHzUv95wLaIOAy4NPUzM7MGN2SgRKEvre6eLgEcB9yY2hcCs9PyrLRO2j5TkrKN2MzMSkkRMXQnaTegFzgM+BbwNWB5OgtB0mTgxxFxhKQHgOMjYkPa9ijwnoh4asA+O4FOgJaWlund3d01FbBl63Y2P1fTTes29ZDxWffX19dHU1NT1n2OlkappVHqANdSVvXU0tHR0RsRbZmHVLNx1XSKiJeAaZKagZuBwwfrlq4HOxvZIbUiogvoAmhra4v29vZqhrKDyxct5pLVVZWR3frT27Pur6enh1rvh7JplFoapQ5wLWXVSLXs0qe8IuJpoAeYATRL6v9LPgnYmJY3AJMB0vbxwNYcgzUzs/Kq5lNeB6UzEyS9HvgA8BBwB3BK6jYXWJyWl6R10vafRDXzamZmNqZVM1c0EViY3kd5HXBDRCyV9CDQLelfgPuAq1P/q4HvSlpHcWZy6jCM28zMSmbIQImIVcC7Bml/DDh6kPY/AnOyjM7MzMYMf1PezMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZeFAMTOzLBwoZmaWhQPFzMyycKCYmVkWDhQzM8vCgWJmZlkMGSiSJku6Q9JDktZIOi+17y/pdklr0/WE1C5Jl0laJ2mVpKOGuwgzMxt91ZyhvAh8NiIOB2YA50h6BzAfWBYRU4BlaR3gBGBKunQCV2YftZmZlc6QgRIRmyJiZVr+A/AQcAgwC1iYui0EZqflWcB1UVgONEuamH3kZmZWKoqI6jtLrcCdwBHAExHRXLFtW0RMkLQUWBARd6X2ZcAFEbFiwL46Kc5gaGlpmd7d3V1TAVu2bmfzczXdtG5TDxmfdX99fX00NTVl3edoaZRaGqUOcC1lVU8tHR0dvRHRlnlINRtXbUdJTcAPgU9HxO8lvWrXQdp2SK2I6AK6ANra2qK9vb3aobzC5YsWc8nqqsvIav3p7Vn319PTQ633Q9k0Si2NUge4lrJqpFqq+pSXpN0pwmRRRNyUmjf3T2Wl6y2pfQMwueLmk4CNeYZrZmZlVc2nvARcDTwUEf9WsWkJMDctzwUWV7SfkT7tNQPYHhGbMo7ZzMxKqJq5omOBTwCrJd2f2v4BWADcIGke8AQwJ227BTgRWAc8C5yVdcRmZlZKQwZKenP91d4wmTlI/wDOqXNcZmY2xvib8mZmloUDxczMsnCgmJlZFg4UMzPLwoFiZmZZOFDMzCwLB4qZmWXhQDEzsywcKGZmloUDxczMsnCgmJlZFg4UMzPLwoFiZmZZOFDMzCwLB4qZmWXhQDEzsywcKGZmloUDxczMsnCgmJlZFg4UMzPLwoFiZmZZDBkokq6RtEXSAxVt+0u6XdLadD0htUvSZZLWSVol6ajhHLyZmZVHNWco1wLHD2ibDyyLiCnAsrQOcAIwJV06gSvzDNPMzMpuyECJiDuBrQOaZwEL0/JCYHZF+3VRWA40S5qYa7BmZlZeioihO0mtwNKIOCKtPx0RzRXbt0XEBElLgQURcVdqXwZcEBErBtlnJ8VZDC0tLdO7u7trKmDL1u1sfq6mm9Zt6iHjs+6vr6+PpqamrPscLY1SS6PUAa6lrOqppaOjozci2jIPqWbjMu9Pg7QNmlgR0QV0AbS1tUV7e3tNB7x80WIuWZ27jOqsP7096/56enqo9X4om0appVHqANdSVo1US62f8trcP5WVrrek9g3A5Ip+k4CNtQ/PzMzGiloDZQkwNy3PBRZXtJ+RPu01A9geEZvqHKOZmY0BQ84VSboeaAcOlLQB+CdgAXCDpHnAE8Cc1P0W4ERgHfAscNYwjNnMzEpoyECJiNNeZdPMQfoGcE69gzIzs7HH35Q3M7MsHChmZpaFA8XMzLJwoJiZWRYOFDMzy8KBYmZmWThQzMwsCweKmZll4UAxM7MsHChmZpaFA8XMzLJwoJiZWRYOFDMzy8KBYmZmWYzO/53bIFrn/yjr/j479UXOrGKf6xeclPW4ZmY5+AzFzMyycKCYmVkWDhQzM8vCgWJmZlk4UMzMLAsHipmZZTEsgSLpeEkPS1onaf5wHMPMzMole6BI2g34FnAC8A7gNEnvyH0cMzMrl+E4QzkaWBcRj0XEC0A3MGsYjmNmZiUyHN+UPwR4smJ9A/CegZ0kdQKdabVP0sM1Hu9A4Kkab1sqf19lLbp4BAZTv0Z5XBqlDnAtZVVPLW/OOZB6DUegaJC22KEhogvoqvtg0oqIaKt3P2XgWsqnUeoA11JWjVTLcEx5bQAmV6xPAjYOw3HMzKxEhiNQ7gWmSDpU0h7AqcCSYTiOmZmVSPYpr4h4UdKngNuA3YBrImJN7uNUqHvarERcS/k0Sh3gWsqqYWpRxA5vb5iZme0yf1PezMyycKCYmVkWIx4okvrS9Rsl3TjSx381kvaT9GtJV1TZv1R1SJom6eeS1khaJeljA7aXbbytkp6TdL+kX0j6maS3VWx/n6R7JP0yXTpTe7uknw/Y1zhJmyVNLGGdb5bUm+pcI+lvR3tMI22kH5NqjifpVklPS1o6AuPZ6XOzot9Fks5Py/8s6QPDPbZqSTpfUkg6cKf9Rvo9FEl9EdE0ogetgqRvAgcBWyPiU1X0L1Udkt4KRESslfRGoBc4PCKeTtvLNt5WYGlEHJHWzwbeGxFzJR0M3APMjoiV6R/xbcAXgR8DjwPvj4j16bbHA5+LiJklrHMPiufZ85KagAco6nzNfJR+pB+Tao4naSawN3B2RHx4mMez0+dmRb+LgL6I+PpwjmdXSZoMXAW8HZgeEa/6JcxRm/JKr1AfSMtnSropvWpYK+mrFf0+mNJ9paQfpCclkr4o6V5JD0jqkqTU/u70KuDnkr5WcYzd0vq9afvZFceYDrQA/z1W64iIRyJibVreCGyhCMhSjncQ+wHb0vI5wLURsTLV8xTweWB+RPwJ+AFQ+SrvVOD6MtYZES9ExPPpcHvyGp5mHqnHBHh9an+LpN+m2zyZ9nerpLXAh4A/5DhePc9NFWciJw9yX10r6ZS0vF7Sl9L4Vkt6e2rfR9I16Xj3SZpVcT//NPVfKem9qf11kr6t4kxpqaRbKo4xXdL/qDibvk3SxIrhXErx/Bv67CMiRvRCkcAArcADaflM4DFgPLAXxSvQyRQ/SXAnsE/qdwHwxbS8f8U+vwv8VVrufwUIsKDiGJ3AF9LynsAK4FCKJ3hPOt6ZwBVjsY4BYzsaeAh4XVnHm8bxHHA/8CiwCXhT6ncTMGtATeMpzh4B3g3cV7HPLcCEMtaZ1icDq4BngXNG+jk32pdReExeSssXApvT8t8AzwPvrDjeHIoz31F9blZsuwg4Py1fC5ySltcD56blTwJXpeV/BT6elpuBR4B9KM689krtU4AVafkU4BaKv3kHU7yAOwXYHfgZcFDq9zGKr3sAnAx8s2IcB+7ssR6On16p1bKI2A4g6UGK36hppvjF4v9NLw72APrnzzskfZ7iztsfWCPpp8C+EfGz1Od7QP/p7AeBd/YnMsU/5CnAScAtEfFkOsZYreNX6ZgTKf7xz43i1XxZx/sI8GhETEvH/hjF5/GPp/j5nsFeDQVARNwrqUnFey6HA8sjYtsg/ctQ568i4sm07Y3Af0q6MSI2DzHe14Lhekw+l5bfDzRLuj/1fx44OCJWpeMdDEwA3lLn8YbruVnppnTdC/x1xfFOVnrfhSIo30TxyyRXSJoGvAS8NW1/H/CDdOzfSLojtb8NOAK4Pd0HuwGbJO0N/GM6TlXKFCjPVyy/RDE2AbdHxGmVHSXtBXwbaEtBcBHFnbmzRBBFyt82YF9zgfdL+iTQBOyhYg621v/HZVTqSPvbD/gRxauk5WUer4r3UCotAf4jLa8B2njlLyxMBx6sWO+mmOo6nAHTXa9i1B6XfhGxUdIaij90o/6BgRIYicdkU0RMk3Rmum3/tPZLFH84yXy8nM/NSv33Vf/91H+8j0bEK35YN411M3AkxdnIHyv6v9q410TEMQP2M5ViNuEXKWgmASslHR0RvxlsR2Wfz10OHCvpMABJe6t4g2uvtP2pNN95CkB6lfoHSTPS9lMr9nUb8HeSdk/7equkfSLi9Ih4U0S0AucD19URJqNWh4o3f29O4/9B2cc7yDHfRzH1BcX/p3NmeoWFpAOAi4GvVvS/Hvg4cBy1/7TPSDwukyT1z+lPAI4Fav1l7deCnI/JncD+/Y8JxdnKwH97T2c83i4/NyV9RdJHqr97XuE24Fzp/9/beVdqH08RpH8CPsHLwXkX8FEV76W0AO2p/WHgIEnHpP3sLunPI2J1RLwhIlrT38cNwFGvFiZQrjOUHUTEb9Mri+sl7ZmavxARj0j6DrCaYl7v3oqbzQO+I+kZivdGtqf2qyjmcFemB+C3wOzhrgFGrI6PAH8BHJCOBXBmRNxf0vEC/FmajhDwAsU8NxGxSdLH0/72Tdu/ERH/VTHGByU9C/RGxDO7WuMI1nk4cImkSHV8PSJW1zLe14JMj0m/buAzwEqKqS14+W/eMRQvYvYE+oAfSeo/CxjJ5+ZUan9B9GXgG8CqdLz1FNNv3wZ+KGkOcAfQ//z4ITCT4j2gR4C7ge0R8YKKabrLJI1P99E3KGYKdknD/fSKpKaI6P8c+nxgYkScN8rD2mVjrY6xNt5avVbqHEtG+jHJeTxJt0XEh7IOcOfHa4qIvnTWfw9w7M7OOHZVqc9QanSSpAspanuc4tMkY9FYq2OsjbdWr5U6x5KRfkyyHW8kwyRZKqmZ4sMHX84ZJtCAZyhmZjY6yv6mvJmZjREOFDMzy8KBYmZmWThQzMwsCweKmZll8X904HyYxalxzAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# What's the dristribution of the main lineages\n",
    "df[\"main_lin\"].hist()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 160,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x26523e23ac8>"
      ]
     },
     "execution_count": 160,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAASz0lEQVR4nO3df5BdZX3H8fcXgoBAQX64g0l0+ZGqKIo1oyA4rKIjP5zCTKFoqRIHzTiFQYSpjbVjR4fpQC2iBKWmRYkMQhB1giAiArcWRZAfgYhRk0KElChFfsgi0ga//eM8C5fNLnt3997szZP3a+bOnvOc557z3PPcfO5znvsjkZlIkuqy1Uw3QJLUfYa7JFXIcJekChnuklQhw12SKjRrphsAsPvuu+fg4OCU7vvkk0+yww47dLdBmjb7pf/YJ/1pOv1y++23P5yZe4y1rS/CfXBwkNtuu21K9221WgwNDXW3QZo2+6X/2Cf9aTr9EhG/Gm+b0zKSVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklShvviG6nSs/O/HWbDo6hk59tqzjpqR40rSRBy5S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SapQx+EeEVtHxJ0RcVVZ3ysibomI1RGxLCJeVMq3LetryvbB3jRdkjSeyYzcPwKsals/Gzg3M+cBjwInlfKTgEczc1/g3FJPkrQJdRTuETEHOAr497IewNuBK0qVpcAxZfnosk7ZflipL0naRDr9VcjPAR8DdirruwGPZeaGsr4OmF2WZwMPAGTmhoh4vNR/uH2HEbEQWAgwMDBAq9Wa0gMY2B7O2H/DxBV7YKpt3hIMDw97fvqMfdKfetUvE4Z7RLwbeCgzb4+IoZHiMapmB9ueK8hcAiwBmD9/fg4NDY2u0pHFlyznnJUz88vFa08YmpHjbg5arRZT7VP1hn3Sn3rVL52k4sHAn0fEkcB2wJ/QjOR3iYhZZfQ+B3iw1F8HzAXWRcQsYGfgka63XJI0rgnn3DPz45k5JzMHgfcAN2TmCcCNwLGl2onA8rJ8ZVmnbL8hMzcauUuSemc6n3P/O+D0iFhDM6d+YSm/ENitlJ8OLJpeEyVJkzWpyerMbAGtsnwv8KYx6vwBOK4LbZMkTZHfUJWkChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalCE4Z7RGwXEbdGxF0RcU9EfKqU7xURt0TE6ohYFhEvKuXblvU1Zftgbx+CJGm0TkbuTwNvz8zXAwcAh0fEgcDZwLmZOQ94FDip1D8JeDQz9wXOLfUkSZvQhOGejeGyuk25JfB24IpSvhQ4piwfXdYp2w+LiOhaiyVJE5rVSaWI2Bq4HdgX+ALwX8BjmbmhVFkHzC7Ls4EHADJzQ0Q8DuwGPDxqnwuBhQADAwO0Wq0pPYCB7eGM/TdMXLEHptrmLcHw8LDnp8/YJ/2pV/3SUbhn5jPAARGxC/At4NVjVSt/xxql50YFmUuAJQDz58/PoaGhTpqykcWXLOeclR09jK5be8LQjBx3c9BqtZhqn6o37JP+1Kt+mdSnZTLzMaAFHAjsEhEjqToHeLAsrwPmApTtOwOPdKOxkqTOdPJpmT3KiJ2I2B54B7AKuBE4tlQ7EVhelq8s65TtN2TmRiN3SVLvdDKfsSewtMy7bwVcnplXRcTPgMsi4kzgTuDCUv9C4OKIWEMzYn9PD9otSXoBE4Z7Zt4NvGGM8nuBN41R/gfguK60TpI0JX5DVZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpApNGO4RMTciboyIVRFxT0R8pJTvGhHXRcTq8vclpTwi4ryIWBMRd0fEn/X6QUiSnq+TkfsG4IzMfDVwIHByROwHLAKuz8x5wPVlHeAIYF65LQQu6HqrJUkvaMJwz8z1mXlHWX4CWAXMBo4GlpZqS4FjyvLRwFez8WNgl4jYs+stlySNa9ZkKkfEIPAG4BZgIDPXQ/MCEBEvLdVmAw+03W1dKVs/al8LaUb2DAwM0Gq1Jt96YGB7OGP/DVO673RNtc1bguHhYc9Pn7FP+lOv+qXjcI+IHYFvAKdl5u8iYtyqY5TlRgWZS4AlAPPnz8+hoaFOm/I8iy9ZzjkrJ/Ua1TVrTxiakeNuDlqtFlPtU/WGfdKfetUvHX1aJiK2oQn2SzLzm6X4NyPTLeXvQ6V8HTC37e5zgAe701xJUic6+bRMABcCqzLzs22brgROLMsnAsvbyt9fPjVzIPD4yPSNJGnT6GQ+42DgfcDKiFhRyv4eOAu4PCJOAu4HjivbvgMcCawBfg98oKstliRNaMJwz8ybGHseHeCwMeoncPI02yVJmga/oSpJFTLcJalChrskVchwl6QKzcy3f7TZGlx0dUf1zth/Aws6rNuJtWcd1bV9SVsCR+6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKjRrogoR8WXg3cBDmfnaUrYrsAwYBNYCf5mZj0ZEAJ8HjgR+DyzIzDt603SpboOLru7q/s7YfwMLOtjn2rOO6upxNTM6GblfBBw+qmwRcH1mzgOuL+sARwDzym0hcEF3milJmowJwz0zfwA8Mqr4aGBpWV4KHNNW/tVs/BjYJSL27FZjJUmdmeqc+0Bmrgcof19aymcDD7TVW1fKJEmb0IRz7pMUY5TlmBUjFtJM3TAwMECr1ZrSAQe2b+YSZ8JU27w56/Rcd7tfPNfT12mfbInneiYNDw/35JxPNdx/ExF7Zub6Mu3yUClfB8xtqzcHeHCsHWTmEmAJwPz583NoaGhKDVl8yXLOWdnt16jOrD1haEaOO5M6eUMOmhDpZr94rqev0z7ZEs/1TGq1Wkw1/17IVKdlrgROLMsnAsvbyt8fjQOBx0embyRJm04nH4W8FBgCdo+IdcA/AmcBl0fEScD9wHGl+ndoPga5huajkB/oQZslSROYMNwz873jbDpsjLoJnDzdRkmSpsdvqEpShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoVmzXQDJGmmDS66esaOfdHhO/Rkv47cJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoV6Eu4RcXhE/CIi1kTEol4cQ5I0vq6He0RsDXwBOALYD3hvROzX7eNIksbXi5H7m4A1mXlvZv4vcBlwdA+OI0kaR2Rmd3cYcSxweGZ+sKy/D3hzZp4yqt5CYGFZfSXwiykecnfg4SneV71jv/Qf+6Q/TadfXpGZe4y1oRc/PxBjlG30CpKZS4Al0z5YxG2ZOX+6+1F32S/9xz7pT73ql15My6wD5ratzwEe7MFxJEnj6EW4/wSYFxF7RcSLgPcAV/bgOJKkcXR9WiYzN0TEKcC1wNbAlzPznm4fp820p3bUE/ZL/7FP+lNP+qXrb6hKkmae31CVpAoZ7pJUob4N94jIiLi4bX1WRPxPRFxV1heU9TsjYnVEXBsRb2mrf1FE3BcRKyLirog4bCYex+YgIj4REfdExN3lfL25i/v+Ufk7GBF/1VY+PyLO69Zx+lVEPFPO6T3leXh6RPTlv7uJ+mR0H07xGKdFxIuns4/NQUTMLfmza1l/SVk/NCKeKrm1KiJujYgT2+43kmsrIuLnEfHRKTciM/vyBgwDdwLbl/UjgBXAVWV9AXB+W/23Ab8GXl3WLwKObdu2eqYfUz/egIOAm4Fty/ruwMt6cJyhkb7bkm7AcNvyS4HvA58ao96sHhx7637rQ2AtsPtM98sm6vuPAUvK8peAjwODwE/b6uxdcu0DZf3ZXAN2o/ly09ypHL8vRxBtrgGOKsvvBS4dr2Jm3kjzrvPCMTbfDMzueuvqsCfwcGY+DZCZD2fmgxHxxoj4j4i4vVwV7QkQEa2IOLuMOH4ZEW8t5a8pZSvKFcC8Uj5cjnMW8Nay/aMRMRQRV0XEVhGxNiJ2GWlQ+cG5gYjYIyK+ERE/KbeDN+mZ6bLMfIjm+XlKNBZExNcj4tvA90bOyUj9iDg/IhaU5SPLSO6miDivvV67ci4/GRE3AcdFxD4R8d3Sj/8ZEa8q9Y6LiJ+Wq4kflLJnj19GmCvK7c6I2ImN+3Cw7POOcntL235aEXFFafMl5fGeCrwMuDEibuzRae4n5wIHRsRpwCHAOaMrZOa9wOnAqWNs+y2whubf6OTN9KvbC7zqDQOvA64AtqN5dRtinJF7KTsGuKYsX8RzI/djgK/N9GPqxxuwYzm3vwS+CBwKbAP8CNij1Dme5iOtAC3gnLJ8JPD9srwYOKEsv4jnrriGy99n+270OvB5nhu5vLltn18DDinLLwdWzfT5msL5HR6j7FFgoDyH1wG7jnOOzi91tgMeAPYq5ZcyzgiaZmT8sbb164F5bef2hrK8EphdlncZo0++DRzc9hyZNUb7XgxsV5bnAbe17edxmi8wbkUzuDqkrX1bxMi9PN530XxD/51lfZC2kfvI+QeeKssLeG7k/nKaf5vbTeXYvfj5ga7JzLsjYpBm1P6dDu4y+qcPPhMR/0xzOXxgd1tXh8wcjog3Am+lmb5aBpwJvBa4LiKg+b7C+ra7fbP8vZ3myQrNP+BPRMQc4JuZuXoSzVgGfBL4Cs2X3paV8ncA+5U2APxJROyUmU9MYt/9qP15el1mPjJB/VcB92bmfWX9Usa+Qh2xDCAidgTeAny97RxuW/7+ELgoIi7nuf5s90PgsxFxCU1/rmvbx4htgPMj4gDgGeBP27bdmpnrSjtW0DxPbprgcdboCJp/O68FrhunzugTe3xEvI3mN7c+lJl/mMqB+zrciyuBf6EZDew2Qd03AKva1v+W5ol7KrAUeGMP2rfZy8xnaEbkrYhYCZwM3JOZB41zl6fL32coz6HM/FpE3EIzjXZtRHwwM2/osAk3A/tGxB40V1lnlvKtgIMy86nJPqZ+FRF705y3h0rRk22bN/D8DzlsN3K3F9jftTRXAbdl+bG+tn1uBTyWmQeMvl9mfjiaN86PAlaUgG7fflZEXE1zdfbjiHjHGIf/KPAb4PXlWO0h9HTb8rPPky1JOafvpBlY3hQRl41TdXRuLcvMUyLiIODqiLgmM3892eP3+5w7wJeBT2fmyheqFBGH0oxm/q29PDP/SHPZv1VEvKtnrdxMRcQrR+bHiwNonmh7lCcXEbFNRLxmgv3sTTO6PI/mBfl1o6o8Aew01n2zuQb9FvBZmqmX35ZN3wOe/TXR0QG0uSkvXv9Kc9k91rcHf0VzpbJtROwMjHzC6+fA3uUqFpppMgAy812ZeUBbsNO27XfAfRFxXDl+RMTry/I+mXlLZn6S8qbdqLbuk5krM/Ns4Daaq4fRfbgzsL78G3sfzRXeRMZ9HtQkmsucC4DTMvN+4DM0g9TR9QZL+eLR2zLzZuBi4CNTaUPfv5qWS7vPj7P5+Ig4hGbu7z7gLzJz1ehKmZkRcSbNu9fX9qyxm6cdgcXlDc0NNG/gLKR5c/q8EjKzgM8BL/QzEscDfx0R/0fzqaVPj9p+N7AhIu6ieT/kzlHbl9H8LtGCtrJTgS9ExN2lDT8APjzJxzfTti/TEtvQnN+LaV7ENpKZD5RpkruB1ZRzlJlPRcTfAN+NiIeBWydx/BOACyLiH0obLgPuopmynEdzVXB9KTu07X6nlamBZ4Cf0Xy44Y88vw+/CHyjvHjcyPOvQsazBLgmItZn5tsm8Tg2Nx8C7s/MkamYL9I8t18B7BMRd9JcmT0BLM7Mr4yzn7OBOyLinyY7HenPD0ibgYjYsbw/EjT/09nqzDx3ptul/rU5TMtIgg+VK4B7aKZDvjTD7VGfc+QuSRVy5C5JFTLcJalChrskVchwl6QKGe6SVKH/B71jx1xlY1i5AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot frequencies of drug resistance profiles\n",
    "df[\"drtype\"].hist()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 161,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "drtype          main_lin         \n",
       "Drug-resistant  lineage2               4\n",
       "                lineage4             103\n",
       "MDR             lineage2              17\n",
       "                lineage3               1\n",
       "                lineage4             411\n",
       "Sensitive                              1\n",
       "                lineage1               1\n",
       "                lineage2               9\n",
       "                lineage2;lineage4      1\n",
       "                lineage4             166\n",
       "                lineageBOV             1\n",
       "XDR             lineage2               2\n",
       "                lineage4              18\n",
       "dtype: int64"
      ]
     },
     "execution_count": 161,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Find the relationship between drtype and main_lin variables\n",
    "grouped = df.groupby(['drtype','main_lin'])\n",
    "grouped.size()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 162,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "drtype          sublin        \n",
       "Drug-resistant  lineage2.2.1       4\n",
       "                lineage4           2\n",
       "                lineage4.1         2\n",
       "                lineage4.1.1       4\n",
       "                lineage4.1.1.1     1\n",
       "                                  ..\n",
       "XDR             lineage4.1.1       1\n",
       "                lineage4.1.1.3     2\n",
       "                lineage4.1.2.1    11\n",
       "                lineage4.3.3       3\n",
       "                lineage4.3.4.2     1\n",
       "Length: 61, dtype: int64"
      ]
     },
     "execution_count": 162,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Find the relationship between drtype and sublin variables\n",
    "\n",
    "grouped = df.groupby(['drtype','sublin'])\n",
    "grouped.size()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 163,
   "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>drtype</th>\n",
       "      <th>Drug-resistant</th>\n",
       "      <th>MDR</th>\n",
       "      <th>Sensitive</th>\n",
       "      <th>XDR</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>row_0</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>ERR3087504</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR3129939</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR3148147</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR3148148</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ERR3148149</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR9224998</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR9225000</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR9225001</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR9225002</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SRR9225003</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>735 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "drtype      Drug-resistant  MDR  Sensitive  XDR\n",
       "row_0                                          \n",
       "ERR3087504               0    0          1    0\n",
       "ERR3129939               0    0          1    0\n",
       "ERR3148147               1    0          0    0\n",
       "ERR3148148               0    0          1    0\n",
       "ERR3148149               1    0          0    0\n",
       "...                    ...  ...        ...  ...\n",
       "SRR9224998               0    0          1    0\n",
       "SRR9225000               0    1          0    0\n",
       "SRR9225001               0    1          0    0\n",
       "SRR9225002               0    1          0    0\n",
       "SRR9225003               0    1          0    0\n",
       "\n",
       "[735 rows x 4 columns]"
      ]
     },
     "execution_count": 163,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## Find the relationship between samples and drtype variables\n",
    "\n",
    "drtype_sample_table = pd.crosstab(index=df.index, columns=df[\"drtype\"])\n",
    "drtype_sample_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 144,
   "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>drtype</th>\n",
       "      <th>Drug-resistant</th>\n",
       "      <th>MDR</th>\n",
       "      <th>Sensitive</th>\n",
       "      <th>XDR</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>main_lin</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage2</th>\n",
       "      <td>4</td>\n",
       "      <td>17</td>\n",
       "      <td>9</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage2;lineage4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage3</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4</th>\n",
       "      <td>103</td>\n",
       "      <td>411</td>\n",
       "      <td>166</td>\n",
       "      <td>18</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineageBOV</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "drtype             Drug-resistant  MDR  Sensitive  XDR\n",
       "main_lin                                              \n",
       "                                0    0          1    0\n",
       "lineage1                        0    0          1    0\n",
       "lineage2                        4   17          9    2\n",
       "lineage2;lineage4               0    0          1    0\n",
       "lineage3                        0    1          0    0\n",
       "lineage4                      103  411        166   18\n",
       "lineageBOV                      0    0          1    0"
      ]
     },
     "execution_count": 144,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## Find the relationship between main lineages and drtype variables\n",
    "\n",
    "drtype_mainlin_table = pd.crosstab(index=df['main_lin'], columns=df[\"drtype\"])\n",
    "drtype_mainlin_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 164,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x26523e539c8>"
      ]
     },
     "execution_count": 164,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAIwCAYAAABTHMhgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3de7hVZbn///cNmJjHPEagAUpqJgIiaB7ykMcsy0PYEU3D1DxgmWx330qr77b0UuJrOzdaou1Ms7JM3bk959lAyLM/wE1KmpIl4kZM4v79MQe0gAUscK01n7XG+3Vd65pzPGPMOe/HgeuzxphjPE9kJpIkqQw9ml2AJEn6J4NZkqSCGMySJBXEYJYkqSAGsyRJBenV7AIANt100+zfv3+zy5AkqVNMmTLlL5m5WWvrigjm/v37M3ny5GaXIUlSp4iIP65onaeyJUkqiMEsSVJBDGZJkgpSxHfMrXnzzTeZPXs2CxYsaHYpWoXevXvTr18/1lprrWaXIkldXrHBPHv2bNZff3369+9PRDS7HK1AZvLyyy8ze/ZsBgwY0OxyJKnLK/ZU9oIFC9hkk00M5cJFBJtssolnNiSpnRQbzICh3EW4nySp/RQdzJIk1Y3B/BZ94xvf4IILLliu/Ve/+hVPPPFEEyqSJHVlBnMHWLhwocEsSVojBvMa+Pa3v822227LBz/4QZ5++mkA9t57b84++2w+8IEP8J3vfIfrr7+eM888kyFDhjBz5kyGDRu25PXTp09n5513BhrDkZ511lmMGDGCESNGMGPGDADmzJnDEUccwS677MIuu+zCvffe2/kdlSR1umJvlyrVlClTuPrqq5k6dSoLFy5k2LBhS0L2lVde4a677gIa4XvooYdy5JFHArDhhhsybdo0hgwZwuWXX84xxxyz5D032GADHnroIa688kpOP/10brjhBk477TTGjh3LHnvswbPPPsuBBx7Ik08+2en9lSR1LoN5Nd1999187GMf4+1vfzsAH/nIR5asGzVq1Apfd/zxx3P55Zdz4YUXcs011/DQQw8tWfeJT3xiyePYsWMBuPXWW5c6Ff7qq68yb9481l9//XbtjySpLAbzGljR7UHrrrvuCl9zxBFHcM4557Dvvvuy8847s8kmm7T6foufL1q0iPvvv5911lmnnaqWJHUFfse8mvbaay+uu+46Xn/9debNm8dvfvObVrdbf/31mTdv3pLl3r17c+CBB3LiiSdy7LHHLrXtNddcs+Rxt912A+CAAw7g4osvXrLNtGnT2rsrkqQCrTKYI2LbiJjW4ufViDg9IjaOiFsiYnr1+I5q+4iICRExIyIeiYhhq/qMrmTYsGGMGjWKIUOGcMQRR7Dnnnu2ut3RRx/N+eefz9ChQ5k5cyYAn/rUp4gIDjjggKW2feONNxg5ciTf+973uOiiiwCYMGECkydPZvDgwbz3ve/lkksu6diOSZKKEJnZ9o0jegJ/AkYCJwN/zczzImIc8I7MPCsiDgFOAQ6ptvteZo5c2fsOHz48J0+evFTbk08+yfbbb79anSndBRdcwNy5c/nmN7+5pK1///5MnjyZTTfdtImVvXXdcX9JUkeJiCmZOby1dav7HfN+wMzM/GNEHAbsXbVfAdwJnAUcBlyZjcR/ICI2iog+mfnCGlXfTXzsYx9j5syZ3H777c0uRZJUsNUN5qOBn1bPt1gctpn5QkRsXrX3BZ5r8ZrZVdtSwRwRY4AxAFtttdVqltH1XHfdda22z5o1q3MLkSQVrc3BHBFvAz4C/MuqNm2lbbnz5Zk5EZgIjVPZba1DktTw5HbN/fpo+6ccW6EjrM5V2QcDD2fmi9XyixHRB6B6fKlqnw1s2eJ1/YDn32qhkiTVweoE8yf452lsgOuB0dXz0cCvW7R/tro6e1dgbt2/X5Ykqa3adCo7It4O7A+c0KL5POBnEXEc8CxwVNV+E40rsmcA84Glb9qVJEkr1KZgzsz5wCbLtL1M4yrtZbdNGrdStav+425s1/ebdd6HVrlNz5492XHHHXnzzTfp1asXo0eP5vTTT6dHj/LGZZk8eTJXXnklEyZMaHX9rFmzuO+++/jkJz+5xp8xfvx4xowZs2Q4UklS+ysvYQqyzjrrMG3aNB5//HFuueUWbrrpJs4555zltlu4cGG7f/Y//vGP1dp++PDhKwxlaATzVVdd9ZZqGj9+PPPnz39L7yFJWjmDuY0233xzJk6cyMUXX0xmMmnSJI466ig+/OEPc8ABB3DnnXdy6KGHLtn+i1/8IpMmTQLgpptuYrvttmOPPfbg1FNPXWq7lvr378+5557LHnvswbXXXsvMmTM56KCD2Hnnndlzzz156qmnALj22mt53/vex0477cRee+0FsNTn33XXXQwZMoQhQ4YwdOhQ5s2bx7hx47j77rsZMmQIF110EbNmzWLPPfdk2LBhDBs2jPvuu2/J++y9994ceeSRbLfddnzqU58iM5kwYQLPP/88++yzD/vss09H/WeWpNpzEovVMHDgQBYtWsRLLzUuQL///vt55JFH2Hjjjbnzzjtbfc2CBQs44YQT+N3vfseAAQOWzCS1Ir179+aee+4BYL/99uOSSy5h0KBBPPjgg5x00kncfvvtnHvuudx888307duXV155Zbn3uOCCC/j+97/P7rvvzmuvvUbv3r0577zzuOCCC7jhhhsAmD9/Prfccgu9e/dm+vTpfOITn2Dx6GtTp07l8ccf513vehe777479957L6eeeioXXnghd9xxR5cfpUySSuYR82pqOYTp/vvvz8Ybb7zS7Z966ikGDhzIgAEDAFYZzIunjnzttde47777OOqooxgyZAgnnHACL7zQuLh9991355hjjuHSSy9t9ZT37rvvzhlnnMGECRN45ZVX6NVr+b+/3nzzTT7/+c+z4447ctRRRy01xeSIESPo168fPXr0YMiQIQ6CIkmdyCPm1fDMM8/Qs2dPNt+8MchZy2kee/XqxaJFi5YsL1iwAFg6yJd14IEH8uKLLzJ8+HAuu+yypd5z0aJFbLTRRq3OKnXJJZfw4IMPcuONNzJkyJDlthk3bhwf+tCHuOmmm9h111259dZbl3uPiy66iC222II//OEPLFq0iN69ey9Zt/baay953rNnzw75Dl2S1DqPmNtozpw5fOELX+CLX/xiq/Mxv/vd7+aJJ57gjTfeYO7cudx2220AbLfddjzzzDNLjjoXT/EIcPPNNzNt2rQlodzSBhtswIABA7j22muBRsD/4Q9/AGDmzJmMHDmSc889l0033ZTnnntuqdfOnDmTHXfckbPOOovhw4fz1FNPLTcN5dy5c+nTpw89evTgxz/+cZsuNlv2PSRJ7a/LHDG35fam9vb6668zZMiQJbdLfeYzn+GMM85oddstt9ySj3/84wwePJhBgwYxdOhQoHFl97//+79z0EEHsemmmzJixIg2f/5PfvITTjzxRL71rW/x5ptvcvTRR7PTTjtx5plnMn36dDKT/fbbj5122om77rpryevGjx/PHXfcQc+ePXnve9/LwQcfTI8ePejVqxc77bQTxxxzDCeddBJHHHEE1157Lfvss89SR/8rMmbMGA4++GD69OnDHXfc0eZ+SOoYH/+X5v4Kf7Spn959rda0jx2lu0/7+Nprr7HeeuuRmZx88skMGjSIsWPHNrusdtWd9pfUVex4xY5N/fxHRxvNa2pl0z56KrsTXHrppQwZMoQddtiBuXPncsIJJ6z6RZKkWuoyp7K7srFjx3a7I2RJUsfwiFmSpIIYzJIkFcRgliSpIAazJEkF6ToXf31jw3Z+v7mr3CQi+PSnP82Pf/xjoDGLVJ8+fRg5ciQ33HADkyZN4swzz6Rfv3689tprDBw4kK9//eu8//3vB+CYY47hrrvuYsMNNyQzufDCC9lvv+VmypQkaQmPmFdi3XXX5bHHHuP1118H4JZbbqFv375LbTNq1CimTp3K9OnTGTduHIcffjhPPvnkkvXnn38+06ZNY/z48XzhC1/o1PolSV2PwbwKBx98MDfeeCMAP/3pT1c6CcU+++zDmDFjmDhx4nLrdtttN/70pz91WJ2SpO7BYF6Fo48+mquvvpoFCxbwyCOPMHLkyJVuP2zYsCXzJrf029/+lo9+9KMdVaYkqZvoOt8xN8ngwYOZNWsWP/3pTznkkENWuf2yQ5yeeeaZfOUrX+Gll17igQce6KgyJUndhEfMbfCRj3yEL3/5y6ucSxlg6tSpS40Zff755zNjxgy+9a1vMXr06I4sU5LUDRjMbfC5z32Or33ta+y448oHjL/rrruYOHEin//855dq79GjB6eddhqLFi3i5ptv7shSJUldXNc5ld2G25s6Sr9+/TjttNNaXXfNNddwzz33MH/+fAYMGMAvfvGLVmdZigi++tWv8t3vfpcDDzywo0uWJHVRTvuoduH+kjqf0z52XU77KElSF2EwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBeky9zG3920Bbb3M/9vf/jZXXXUVPXv2pEePHvzHf/zHKsfLbqv3v//93HfffcyaNYv77ruPT37ykwBMnjyZK6+8kgkTJrTL50iSuo4uE8zNcP/993PDDTfw8MMPs/baa/OXv/yFv//97+32/vfddx8As2bN4qqrrloSzMOHD2f48FZvb5MkdXOeyl6JF154gU033ZS1114bgE033ZR3vetdTJkyhQ984APsvPPOHHjggbzwwgsA7L333px11lmMGDGC97znPdx9990APP7444wYMYIhQ4YwePBgpk+fDsB6660HwLhx47j77rsZMmQIF110EXfeeSeHHnooixYton///rzyyitLatpmm2148cUXmTNnDkcccQS77LILu+yyC/fee29n/qeRJHUQg3klDjjgAJ577jne8573cNJJJ3HXXXfx5ptvcsopp/Dzn/+cKVOm8LnPfY5//dd/XfKahQsX8tBDDzF+/HjOOeccAC655BJOO+00pk2bxuTJk+nXr99Sn3Peeeex5557Mm3aNMaOHbukvUePHhx22GFcd911ADz44IP079+fLbbYgtNOO42xY8fy+9//nl/84hccf/zxnfBfRJLU0TyVvRLrrbceU6ZM4e677+aOO+5g1KhRfPWrX+Wxxx5j//33B+Af//gHffr0WfKaww8/HICdd96ZWbNmAbDbbrvx7W9/m9mzZ3P44YczaNCgNtcwatQozj33XI499liuvvpqRo0aBcCtt97KE088sWS7V199lXnz5rH++uu/1W5LkprIYF6Fnj17svfee7P33nuz44478v3vf58ddtiB+++/v9XtF5/27tmzJwsXLgTgk5/8JCNHjuTGG2/kwAMP5LLLLmPfffdt0+fvtttuzJgxgzlz5vCrX/2Kr371qwAsWrSI+++/n3XWWacdeilJKoWnslfi6aefXvJ9MMC0adPYfvvtmTNnzpJgfvPNN3n88cdX+j7PPPMMAwcO5NRTT+UjH/kIjzzyyFLr119/febNm9fqayOCj33sY5xxxhlsv/32bLLJJkDjNPvFF1+8VG2SpK6vyxwxN2MWk9dee41TTjmFV155hV69erHNNtswceJExowZw6mnnsrcuXNZuHAhp59+OjvssMMK3+eaa67hP//zP1lrrbV45zvfyde+9rWl1g8ePJhevXqx0047ccwxxzB06NCl1o8aNYpddtmFSZMmLWmbMGECJ598MoMHD2bhwoXstddeXHLJJe3af0lS53PaR7UL95fU+Zz2sety2kdJkroIg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSpIl7mP+cnt2vdWnO2fenKl65977jn22msvpkyZwsYbb8zf/vY3hg0bxqRJkzjooIPYbrvtWLBgAeuvvz4nn3wyo0ePBmDSpEmceeaZ9O3blwULFnDCCScsNf61JEkr02WCubNtueWWnHjiiYwbN46JEycybtw4xowZw7vf/W623nprpk6dCjRG9Tr88MNZtGgRxx57LNAYEOTiiy/m5ZdfZtttt+XII49kyy23bGZ3JEldhKeyV2Ls2LE88MADjB8/nnvuuYcvfelLy20zcOBALrzwQiZMmLDcuk022YRtttlmybSQkiStikfMK7HWWmtx/vnnc9BBB/Hf//3fvO1tb2t1u2HDhvHUU08t1/7ss8+yYMECBg8e3NGlSpK6CY+YV+G//uu/6NOnD4899tgKt1l2WNNrrrmGHXbYgYEDB3LaaafRu3fvji5TktRNGMwrMW3aNG655RYeeOABLrroohWekp46depS40SPGjWKxx9/nLvvvpsvfelL/PnPf+6skiVJXVybgjkiNoqIn0fEUxHxZETsFhEbR8QtETG9enxHtW1ExISImBERj0TEsI7tQsfITE488UTGjx/PVlttxZlnnsmXv/zl5babNWsWX/7ylznllFOWW7fbbrvxmc98hu9973udUbIkqRto63fM3wN+m5lHRsTbgLcDZwO3ZeZ5ETEOGAecBRwMDKp+RgI/qB7fklXd3tTeLr30Urbaaiv2339/AE466SQmTZrEH//4R2bOnMnQoUOX3C51yimnLLkie1lnnXUWw4YN4+yzz2b99dfvzC5IkrqgVU77GBEbAH8ABmaLjSPiaWDvzHwhIvoAd2bmthHxH9Xzny673Yo+w2kfuz73l9T5nPax63qr0z4OBOYAl0fE1Ii4LCLWBbZYHLbV4+bV9n2B51q8fnbVtmxRYyJickRMnjNnzmp0R5Kk7qstwdwLGAb8IDOHAv9L47T1ikQrbcsdlmfmxMwcnpnDN9tsszYVK0lSd9eWYJ4NzM7MB6vln9MI6herU9hUjy+12L7lMFf9gOfXpLhVnWZXGdxPktR+VhnMmfln4LmI2LZq2g94ArgeGF21jQZ+XT2/HvhsdXX2rsDclX2/vCK9e/fm5Zdf9pd+4TKTl19+2Xu1JamdtPWq7FOAn1RXZD8DHEsj1H8WEccBzwJHVdveBBwCzADmV9uutn79+jF79mz8/rl8vXv3pl+/fs0uQ5K6hTYFc2ZOA1q7emy/VrZN4OS3WBdrrbUWAwYMeKtvI0lSl+LIX5IkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQXo1uwBJ0pp59H+ebXYJ6gAeMUuSVBCDWZKkghjMkiQVxGCWJKkgbQrmiJgVEY9GxLSImFy1bRwRt0TE9OrxHVV7RMSEiJgREY9ExLCO7IAkSd3J6hwx75OZQzJzeLU8DrgtMwcBt1XLAAcDg6qfMcAP2qtYSZK6u7dyKvsw4Irq+RXAR1u0X5kNDwAbRUSft/A5kiTVRluDOYH/jogpETGmatsiM18AqB43r9r7As+1eO3sqk2SJK1CWwcY2T0zn4+IzYFbIuKplWwbrbTlchs1An4MwFZbbdXGMiRJ6t7adMScmc9Xjy8B1wEjgBcXn6KuHl+qNp8NbNni5f2A51t5z4mZOTwzh2+22WZr3gNJkrqRVQZzRKwbEesvfg4cADwGXA+MrjYbDfy6en498Nnq6uxdgbmLT3lLkqSVa8up7C2A6yJi8fZXZeZvI+L3wM8i4jjgWeCoavubgEOAGcB84Nh2r1qSpG5qlcGcmc8AO7XS/jKwXyvtCZzcLtVJklQzjvwlSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVJA2B3NE9IyIqRFxQ7U8ICIejIjpEXFNRLytal+7Wp5Rre/fMaVLktT9rM4R82nAky2WvwNclJmDgL8Bx1XtxwF/y8xtgIuq7SRJUhu0KZgjoh/wIeCyajmAfYGfV5tcAXy0en5YtUy1fr9qe0mStAptPWIeD3wFWFQtbwK8kpkLq+XZQN/qeV/gOYBq/dxq+6VExJiImBwRk+fMmbOG5UuS1L2sMpgj4lDgpcyc0rK5lU2zDev+2ZA5MTOHZ+bwzTbbrE3FSpLU3fVqwza7Ax+JiEOA3sAGNI6gN4qIXtVRcT/g+Wr72cCWwOyI6AVsCPy13SuXJKkbWuURc2b+S2b2y8z+wNHA7Zn5KeAO4Mhqs9HAr6vn11fLVOtvz8zljpglSdLy3sp9zGcBZ0TEDBrfIf+wav8hsEnVfgYw7q2VKElSfbTlVPYSmXkncGf1/BlgRCvbLACOaofaJEmqHUf+kiSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFWWUwR0TviHgoIv4QEY9HxDlV+4CIeDAipkfENRHxtqp97Wp5RrW+f8d2QZKk7qMtR8xvAPtm5k7AEOCgiNgV+A5wUWYOAv4GHFdtfxzwt8zcBrio2k6SJLXBKoM5G16rFteqfhLYF/h51X4F8NHq+WHVMtX6/SIi2q1iSZK6sTZ9xxwRPSNiGvAScAswE3glMxdWm8wG+lbP+wLPAVTr5wKbtPKeYyJickRMnjNnzlvrhSRJ3USbgjkz/5GZQ4B+wAhg+9Y2qx5bOzrO5RoyJ2bm8Mwcvtlmm7W1XkmSurXVuio7M18B7gR2BTaKiF7Vqn7A89Xz2cCWANX6DYG/tkexkiR1d225KnuziNioer4O8EHgSeAO4Mhqs9HAr6vn11fLVOtvz8zljpglSdLyeq16E/oAV0RETxpB/rPMvCEingCujohvAVOBH1bb/xD4cUTMoHGkfHQH1C1JUre0ymDOzEeAoa20P0Pj++Zl2xcAR7VLdZIk1Ywjf0mSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkgqwymCNiy4i4IyKejIjHI+K0qn3jiLglIqZXj++o2iMiJkTEjIh4JCKGdXQnJEnqLtpyxLwQ+FJmbg/sCpwcEe8FxgG3ZeYg4LZqGeBgYFD1Mwb4QbtXLUlSN7XKYM7MFzLz4er5POBJoC9wGHBFtdkVwEer54cBV2bDA8BGEdGn3SuXJKkbWq3vmCOiPzAUeBDYIjNfgEZ4A5tXm/UFnmvxstlVmyRJWoU2B3NErAf8Ajg9M19d2aattGUr7zcmIiZHxOQ5c+a0tQxJkrq1NgVzRKxFI5R/kpm/rJpfXHyKunp8qWqfDWzZ4uX9gOeXfc/MnJiZwzNz+Gabbbam9UuS1K205arsAH4IPJmZF7ZYdT0wuno+Gvh1i/bPVldn7wrMXXzKW5IkrVyvNmyzO/AZ4NGImFa1nQ2cB/wsIo4DngWOqtbdBBwCzADmA8e2a8WSJHVjqwzmzLyH1r83Btivle0TOPkt1iVJUi058pckSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQVYZzBHxo4h4KSIea9G2cUTcEhHTq8d3VO0RERMiYkZEPBIRwzqyeEmSuptebdhmEnAxcGWLtnHAbZl5XkSMq5bPAg4GBlU/I4EfVI+SpHbWf8FVTf38WU399O5rlUfMmfk74K/LNB8GXFE9vwL4aIv2K7PhAWCjiOjTXsVKktTdrel3zFtk5gsA1ePmVXtf4LkW282u2pYTEWMiYnJETJ4zZ84aliFJUvfS3hd/RStt2dqGmTkxM4dn5vDNNtusncuQJKlrWtNgfnHxKerq8SmtEHAAABJwSURBVKWqfTawZYvt+gHPr3l5kiTVy5oG8/XA6Or5aODXLdo/W12dvSswd/Epb0mStGqrvCo7In4K7A1sGhGzga8D5wE/i4jjgGeBo6rNbwIOAWYA84FjO6BmSZK6rVUGc2Z+YgWr9mtl2wROfqtFSZJUV478JUlSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEYJYkqSAGsyRJBTGYJUkqiMEsSVJBDGZJkgpiMEuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIKYjBLklQQg1mSpIIYzJIkFcRgliSpIAazJEkFMZglSSqIwSxJUkEMZkmSCtKr2QVIb9WT223ftM/e/qknm/bZkrong1ld3sf/pXn/jB9t2idL6q48lS1JUkEMZkmSCmIwS5JUEINZkqSCGMySJBXEq7LV5T36P882uwRJajceMUuSVBCDWZKkghjMkiQVxGCWJKkgBrMkSQUxmCVJKojBLElSQQxmSZIK0iHBHBEHRcTTETEjIsZ1xGdIktQdtfvIXxHRE/g+sD8wG/h9RFyfmU+092dJAP0XXNW0z57VtE+W1F11xBHzCGBGZj6TmX8HrgYO64DPkSSp2+mIsbL7As+1WJ4NjFx2o4gYA4ypFl+LiKc7oJa22hT4SxM/v9ns/xr2P77TzpV0Pve9/V/j/vvv/y1594pWdEQwRyttuVxD5kRgYgd8/mqLiMmZObzZdTSL/a9v/+vcd7D/9r/M/nfEqezZwJYtlvsBz3fA50iS1O10RDD/HhgUEQMi4m3A0cD1HfA5kiR1O+1+KjszF0bEF4GbgZ7AjzLz8fb+nHZWxCn1JrL/9VXnvoP9t/8Fiszlvv6VJElN4shfkiQVxGCWJKkgBrMkSQUxmCVJ3V5EbNHsGtrKYNYSEbF/s2voDBGxQURs3Ur74GbU00wRsXmza+hMEfHOiHhn9XyziDg8InZodl3NEhH/t9k1dKI/RMQtEfG5iNiw2cWsjFdla4mIeDYzt2p2HR0pIj4OjAdeAtYCjsnM31frHs7MYc2sryNFxMbLNgFTgKE0fhf8tfOr6jwRcQIwjka/vwMcAzwO7A58NzN/2LzqOl5ETFi2CfgMcCVAZp7a6UV1omqCpQ/SGFvjEOB+4KfA9Zn5ejNrW5bBXDMRsaLBXgLYNzPX7cx6OltETAMOzswXImIEjV9KZ2fmLyNiamYObXKJHSYiFgF/XKa5H43R+jIzB3Z+VZ0nIh6lMW7/OjT+O2yTmX+OiHcAd2TmkKYW2MEiYjZwJ/Df/HPo5AuALwNk5hXNqazzVYNfHUwjpPcBbsvMTzW3qn/qiLGyVbY9gU8Dry3THjRmBuvuembmCwCZ+VBE7APcEBH9aGVM927mKzSOGM7MzEcBIuJ/MnNAc8vqNG9m5nxgfkTMzMw/A2Tm3yKiu+97gO2BbwIH0fg38KeI+HqdAnmxzPx7RDwBPAnsDLy3ySUtxWCunweA+Zl517IrmjzDV2eZFxFbZ+ZMgOrIeW/gV0C3/q4xMy+IiKuBiyLiOeDrdP8/RlpaFBFrZeabwIcWN0ZEb2pwvU1mzgNOj4idgf+MiBupQb9bioitgFHAJ4B1qaYlzswnm1rYMjyVrVqJiJ2A/83MGcu0rwV8PDN/0pzKOldEfBj4V6B/Zr6z2fV0huqX8vOZuXCZ9r7A9pl5a3Mq63wREcBJwG6Z+elm19MZIuI+GtMS/xz4aWZObnJJK2Qw11hEvBsYlJm3RsQ6QK/qr+pasP+xDrB1Zj7W7Fo62zL7/u00vuKo076v3b/9iPgA8LvsAqFXq9MY+qeI+DyNvxz/o2rqR+N0bi3Uuf8RsUVE/BD4ZWY+FhHvjYjjml1XZ2ll3/elJvse6vtvv/r67qCI+F1E/CUi5kTEXRFxSLNrW5bBXF8n07hN5FWAzJwO1Ome1jr3fxKN2d/eVS3/f8DpTaum89V530NN+1/9QfJN4BvAQGBr4BzgGxExpomlLcdgrq83MvPvixciohf1uhCozv3fNDN/BiyCxlStwD+aW1KnqvO+h/r2fyxwQGbenpmvVj+307htamyTa1uKwVxfd0XE2cA61Yhf1wK/aXJNnanO/f/fiNiE6pdxROwKzG1uSZ2qzvse6tv/VgfRycyXm1HMynjxV01FRA/gOOAAGvcw3wxc1hUujGgPde5/RAwD/h/wPuAxYDPgyMx8pKmFdZI673uob/8j4kFgTGb+YZn2nYBLM7OYcRwMZqmGqtOX29L4xfx0dW+v1G1FxB7AT4DLaQxFm8AuwGjg05l5TxPLW4rBXFPV8ITL7vy5wGTgWyWe3mlPde5/RBzeSvNc4NHMfKmz6+lsdd73UO/+VxOYnERjMKGgMVb69xePAlcKg7mmIuK7NC74uapqOrp6fBXYIzM/3JTCOkmd+1+N+LQbcEfVtDeNEeHeA5ybmT9uUmmdos77Hux/V2Aw11RE3JuZu7fWFhGPZuaOzaqtM9S5/xHxG+D4zHyxWt4C+AFwPI0BGN7XzPo6Wp33PdS3/xExiMZod38FLgQupTF3wEwa/z/8vonlLcWrsutrvYgYuXihmmlpvWpxYesv6Vbq3P/+i0O58hLwnuqK1Tp811znfQ/17f/lwH3A88CDwI+ATWnMrnVxE+tajpNY1NfxwI8iYj0a37W8ChwfEesC/9bUyjpHnft/d0TcQOM2GYAjgN9VfX+leWV1mjrve6hv/9fLzIkAEfGFzFz87/+WiDi/iXUtx1PZNRcRG9L4d1CHX8jLqWP/qwkMjqAx+lMA9wC/6O63yyyrjvu+pbr1PyIezsxhyz5vbbnZDOYai4gP0bg6sffitsw8t3kVda6697/O6r7v69j/iJgPzKDxx+jW1XOq5YGZuW6zaluWp7JrKiIuAd4O7ANcBhwJPNTUojpRnftfjfT1/4DtgbcBPWlMhblBUwvrJHXe91Dr/m/f7ALayiPmmoqIRzJzcIvH9WjMNnRAs2vrDHXuf0RMpnGLzLXAcOCzwDaZ+a9NLayT1Hnfg/2H8qe99Krs+nq9epwfEe+icTXugCbW09lq3f/MnEFjDuJ/ZOblNI6e6qLW+56a978rTHvpqez6uiEiNgLOBx6mMRLQZc0tqVPVuf/zI+JtwLRqsIkXgGK+X+sEdd73YP9PBkbQuGWKzJweEUVNe+mpbBERawO9M7NOMwwtUbf+V6fxXqTx/fJYYEPg36uj6Fqp275fVh37HxEPZubIiJiamUOrceMfzszBza5tMU9l11REvD0i/k9EXJqZbwCbR8Shza6rs9S5/5n5RxpXovbJzHMy84w6hXKd9z3Yf7rAtJcGc31dDrxBY8xkgNnAt5pXTqerbf8j4sPANOC31fKQiLi+uVV1qtru+0rd+z8OmAM8CpwA3AR8takVLcNgrq+tM/O7VEMwZubrNI6i6qLO/f8Gje/YXgHIzGlA/ybW09nqvO+h5v3PzEWZeWlmHpWZR1bPi/pO14u/6uvv1W0CCRARW9P4K7ou6tz/hZk5tzEAWC3Ved9DzfvfFaa9NJjr6+s0TmVuGRE/oTE84zFNrahz1bn/j0XEJ4Ge1Yw7p9IY3L8u6rzvwf7/Fyue9nIS0PRpL70qu8YiYhNgVxqnsR7IzL80uaROVdf+R8TbaUx/dwCNvt8MfDMzFzS1sE5U132/WJ373xWmvfSIud56A3+j8e/gvRFBZv6uyTV1plr2PzPn0wjmWoz0tQK13Pct1Ln/60XEyMx8EMqc9tJgrqmI+A4wCngcWFQ1J1CL/znr3P+IeA+NOWj70+J3QGbu26yaOlOd9z3Yf7rAtJeeyq6piHgaGFzdx1g7de5/RPwBuASYQuO7NgAyc0rTiupEdd73YP8XK3naS4+Y6+sZYC1qdDXmMurc/4WZ+YNmF9FEdd73YP+XmvZy8d0JJU17aTDX13waYyXfRov/QTPz1OaV1Knq3P/fRMRJwHUs3fe/Nq+kTlXnfQ81739XmPbSU9k1FRGjW2vPzCs6u5ZmqHP/I+J/WmnOzBzY6cU0QZ33Pdj/rjDtpcEsSaqNFpNYPAAcDrwMPJaZg5pc2hKeyq6ZiPhZZn58BaPfUNIMKx2hzv2PiH0z8/aIOLy19Zn5y86uqTPVed+D/W+h+GkvPWKumYjok5kvVFP/LaeaeajbqnP/I+KczPx6RFzeyurMzM91elGdqM77Hux/a0qd9tJgliTVRjXy3ZeArTLz89WwtNtm5g1NLm0JT2XXTETMo5XTWDRutM/M3KCTS+pUde5/RJyxsvWZeWFn1dIMdd73YP9buJzGPfwtp728FjCY1RyZuX6za2immve/zn2v+76vff9b2DozR0XEJ6Ax7WUUNtWawSzVRGae0+wapAIUP+1lj2YXIKlzRcShyywfFhEjm1WP1MmWnfbyNuArzS1paR4xS/WzC0t/nzYS2DEiemXmwU2qSeoUmXlLRDzMP6e9PK20aS+9KluSVCsR0Rd4N0vPrlbM7FoeMUs1EhEbAJtl5sxl2gdn5iNNKkvqNF1h2kuPmKWaiIiPA+OBl2jMLnRMZv6+WvdwZg5rZn1SZ+gK01568ZdUH2cDO2fmEOBY4Mcthucs6nYRqQMtnvayWJ7KluqjZ2a+AJCZD0XEPjTGDe5H6wNPSN1R8dNeGsxSfcyLiK0Xf79cjZu8N/ArGpPGS3VwffVTLL9jlmoiInYC/jczZyzTvhbw8cz8SXMqk9SSwSzVUDXD0KDMvLUaBalXZs5rdl1SR+lK014azFLNRMTngTHAxpm5dTW7ziWZuV+TS5M6TFea9tJglmomIqYBI4AHM3No1fZoZu7Y3MokgRd/SXX0Rmb+ffGEOhHRC6/KVjfXlaa9NJil+rkrIs4G1omI/YGTgN80uSapQ3WlaS89lS3VTET0AI4DDqBxtHAzcFn6y0AqgsEsSVJBPJUt1cwKbheZC0wGvpWZL3d+VZIWM5il+vkv4B/AVdXy0dXjq8Ak4MNNqElSxVPZUs1ExL2ZuXtrbd42JTWfs0tJ9bNeRIxcvBARI4D1qsWFzSlJ0mKeypbq53jgRxGxHo2rsl8Fjo+IdYF/a2plkjyVLdVVRGxI43fAK82uRdI/GcxSDUXEh2hM9dh7cVtmntu8iiQt5nfMUs1ExCXAKOAUGqeyjwJaHdhfUufziFmqmYh4JDMHt3hcD/hlZh7Q7NokecQs1dHr1eP8iHgX8CYwoIn1SGrBq7Kl+rkhIjYCzgcepjEK2GXNLUnSYp7KlmosItYGemfm3GbXIqnBU9lSzUTE2yPi/0TEpZn5BrB5RBza7LokNRjMUv1cDrwB7FYtzwa+1bxyJLVkMEv1s3VmfpfGRV9k5us0bpuSVACDWaqfv0fEOlRTP0bE1jSOoCUVwKuypfr5OvBbYMuI+AmwO3BMUyuStIRXZUs1FBGbALvSOIX9QGb+pcklSaoYzFINRURfGsNwLjlrlpm/a15FkhbzVLZUMxHxHRpjZT8OLKqaEzCYpQJ4xCzVTEQ8DQyu7mGWVBivypbq5xlgrWYXIal1nsqW6mc+MC0ibqPFbVKZeWrzSpK0mMEs1c/11Y+kAvkdsyRJBfGIWaqJiPhZZn48Ih6lGvWrpcwc3ISyJC3DI2apJiKiT2a+EBHvbm19Zv6xs2uStDyDWZKkgngqW6qJiJhHK6ewaQzLmZm5QSeXJKkVHjFLklQQBxiRJKkgBrMkSQUxmCVJKojBLNVERAyPiAlr+NrXqsd3RcTP27cySS158ZekVYqI1zJzvWbXIdWBR8xSFxIR/SPiqYi4LCIei4ifRMQHI+LeiJgeESOqn/siYmr1uG312r0j4obq+Tci4kcRcWdEPBMRbZrAovr8x6rnx0TELyPit9Vnf7fjei7Vh8EsdT3bAN8DBgPbAZ8E9gC+DJwNPAXslZlDga8B/3cF77MdcCAwAvh6RKzJVJBDgFHAjsCoiNhyDd5DUgsOMCJ1Pf+TmY8CRMTjwG2ZmdUY2P2BDYErImIQjQFFVhS4N2bmG8AbEfESsAUwezVruS0z51a1PAG8G3hudTsk6Z88Ypa6njdaPF/UYnkRjT+2vwnckZnvAz4M9G7D+/yDNftDvT3eQ1ILBrPU/WwI/Kl6fkwT65C0Bgxmqfv5LvBvEXEv0LPZxUhaPd4uJUlSQTxiliSpIF6oIQmAiNgEuK2VVftl5sudXY9UV57KliSpIJ7KliSpIAazJEkFMZglSSqIwSxJUkH+fw6hTZpuUqVzAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Stacked Column Chart: visual form of the two-way table\n",
    "\n",
    "drtype_mainlin_table.plot(kind=\"bar\", \n",
    "                 figsize=(8,8),\n",
    "                 stacked=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 165,
   "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>drtype</th>\n",
       "      <th>Drug-resistant</th>\n",
       "      <th>MDR</th>\n",
       "      <th>Sensitive</th>\n",
       "      <th>XDR</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>sublin</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage1.2.1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage2.2.1</th>\n",
       "      <td>4</td>\n",
       "      <td>16</td>\n",
       "      <td>9</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage2.2.2</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage3</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4</th>\n",
       "      <td>2</td>\n",
       "      <td>18</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.1</th>\n",
       "      <td>4</td>\n",
       "      <td>5</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.1.1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.1.3</th>\n",
       "      <td>10</td>\n",
       "      <td>16</td>\n",
       "      <td>4</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.2</th>\n",
       "      <td>4</td>\n",
       "      <td>5</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.2.1</th>\n",
       "      <td>28</td>\n",
       "      <td>171</td>\n",
       "      <td>38</td>\n",
       "      <td>11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.1.2.1;lineage4.1.1</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.2</th>\n",
       "      <td>5</td>\n",
       "      <td>34</td>\n",
       "      <td>15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.3</th>\n",
       "      <td>32</td>\n",
       "      <td>110</td>\n",
       "      <td>44</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.3;lineage4.1.1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.3;lineage4.1.1.3</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.3;lineage4.1.2.1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.4.1</th>\n",
       "      <td>5</td>\n",
       "      <td>9</td>\n",
       "      <td>12</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.4.2</th>\n",
       "      <td>6</td>\n",
       "      <td>37</td>\n",
       "      <td>20</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.4.2.1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.3.4.2;lineage4.3.3</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.4.1.1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.4.2</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.7</th>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.8</th>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.8;lineage4.1.2.1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4.9</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineage4;lineage2.2.1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lineageBOV</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "drtype                       Drug-resistant  MDR  Sensitive  XDR\n",
       "sublin                                                          \n",
       "                                          0    0          1    0\n",
       "lineage1.2.1                              0    0          1    0\n",
       "lineage2.2.1                              4   16          9    2\n",
       "lineage2.2.2                              0    1          0    0\n",
       "lineage3                                  0    1          0    0\n",
       "lineage4                                  2   18          2    0\n",
       "lineage4.1                                2    0          1    0\n",
       "lineage4.1.1                              4    5          3    1\n",
       "lineage4.1.1.1                            1    0          0    0\n",
       "lineage4.1.1.3                           10   16          4    2\n",
       "lineage4.1.2                              4    5          6    0\n",
       "lineage4.1.2.1                           28  171         38   11\n",
       "lineage4.1.2.1;lineage4.1.1               0    1          0    0\n",
       "lineage4.3                                1    0          0    0\n",
       "lineage4.3.2                              5   34         15    0\n",
       "lineage4.3.3                             32  110         44    3\n",
       "lineage4.3.3;lineage4.1.1                 1    0          0    0\n",
       "lineage4.3.3;lineage4.1.1.3               0    0          1    0\n",
       "lineage4.3.3;lineage4.1.2.1               1    0          0    0\n",
       "lineage4.3.4                              0    0          3    0\n",
       "lineage4.3.4.1                            5    9         12    0\n",
       "lineage4.3.4.2                            6   37         20    1\n",
       "lineage4.3.4.2.1                          0    0          1    0\n",
       "lineage4.3.4.2;lineage4.3.3               0    0          1    0\n",
       "lineage4.4.1.1                            1    0          5    0\n",
       "lineage4.4.2                              0    0          1    0\n",
       "lineage4.7                                0    2          0    0\n",
       "lineage4.8                                0    3          6    0\n",
       "lineage4.8;lineage4.1.2.1                 0    0          1    0\n",
       "lineage4.9                                0    0          2    0\n",
       "lineage4;lineage2.2.1                     0    0          1    0\n",
       "lineageBOV                                0    0          1    0"
      ]
     },
     "execution_count": 165,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## Find the relationship between sub lineages and drtype variables\n",
    "\n",
    "drtype_sublin_table = pd.crosstab(index=df['sublin'], columns=df[\"drtype\"])\n",
    "drtype_sublin_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 147,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x2651df78088>"
      ]
     },
     "execution_count": 147,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAJfCAYAAACqvnXvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzde7xc873/8dcniQoSVIKmCSKkEppkSyKh7lV3pa5B67iUqGtEqVQdt3JOWw6Rqmq0Fdq6VJVqUEKJS1yaSCQinARBSEkdIn5EE/n8/lhrb5OduX5nZu3v7Lyfj8c89p416zPr+5211nxm3T7L3B0RERGJQ4e2boCIiIh8TolZREQkIkrMIiIiEVFiFhERiYgSs4iISEQ6tXUDALp37+69e/du62aIiIhkZtq0af9y9w1bD48iMffu3ZupU6e2dTNEREQyY2av5xuuXdkiIiIRUWIWERGJiBKziIhIRKI4xiwiIuVbtmwZCxYsYOnSpW3dFClD586d6dWrF2ussUZZ4ysxi4g0mAULFtC1a1d69+6NmbV1c6QId+e9995jwYIFbL755mXFaFe2iEiDWbp0Kd26dVNSbgBmRrdu3Srau6HELCLSgJSUG0el80qJWUREJCJKzCIiUjcXX3wxV1555SrD7777bl588cU2aFH8lJhFRCRTy5cvV2IuQolZRERq6vLLL2errbbiG9/4Bi+//DIAu+22G+effz677rorP/3pT7nnnns499xzaWpq4pVXXmHw4MEt8XPnzmXIkCFAUrL5vPPOY9iwYQwbNox58+YBsGjRIg499FC22247tttuO5588snsO1onulxKRERqZtq0adx2221Mnz6d5cuXM3jw4JYk+8EHHzB58mQgSb4HHHAAhx12GADrrbceM2bMoKmpiRtvvJHjjjuu5T3XXXddnn32WW6++WbOOussJk6cyKhRoxg9ejQ77bQTb7zxBnvvvTdz5szJvL/1oMQsIiI18/jjj3PwwQez9tprA3DggQe2vDZixIiCcSeeeCI33ngjV111FbfffjvPPvtsy2tHHXVUy9/Ro0cD8NBDD620K/zDDz9kyZIldO3atab9aQsld2Wb2SZm9oiZzTGz2WY2Kh1+sZm9ZWYz0sd+OTE/NLN5Zvayme1dzw6IiEhcCl0etM466xSMOfTQQ7n//vuZOHEiQ4YMoVu3bnnfr/n/FStW8NRTTzFjxgxmzJjBW2+91S6SMpR3jHk58H137w9sD5xmZlunr13t7k3p4z6A9LUjgW2AfYDrzKxjHdouIiKR2WWXXbjrrrv45JNPWLJkCX/961/zjte1a1eWLFnS8rxz587svffenHLKKRx//PErjXv77be3/N1hhx0A2Guvvbj22mtbxpkxY0atu9JmSiZmd1/o7s+l/y8B5gA9i4QcBNzm7p+6+2vAPGBYLRorIiJxGzx4MCNGjKCpqYlDDz2UnXfeOe94Rx55JFdccQXbbrstr7zyCgDf/va3MTP22muvlcb99NNPGT58ONdccw1XX301AOPGjWPq1KkMHDiQrbfemuuvv76+HcuQuXv5I5v1Bh4DvgqcDRwHfAhMJdmqft/MrgWedvffpzG/Ae539z+1eq+RwEiATTfddMjrr+e9X7SIiLQyZ84c+vfv39bNqLkrr7ySxYsX8+Mf/7hlWO/evZk6dSrdu3dvw5ZVL988M7Np7j609bhlXy5lZl2AO4Gz3P1D4JfAFkATsBD4n+ZR84Svkv3dfby7D3X3oRtuuGG5zRARkXbo4IMP5uabb2bUqFFt3ZQ2V9ZZ2Wa2BklS/oO7/xnA3d/Jef0GYGL6dAGwSU54L+DtmrRWRETapbvuuivv8Pnz52fbkAiUTMyWnAL3G2COu1+VM7yHuy9Mnx4MvJD+fw9wi5ldBXwZ6As8i0hE5vQrvBuw/0vt41pIEWlM5Wwx7wgcA8wys+bT3s4HjjKzJpLd1POBkwHcfbaZ/RF4keSM7tPc/bNaN1xERKQ9KpmY3f0J8h83vq9IzOXA5VW0S0REZLWkWtkiIiIRUUlOEZEG13vMvTV9v/k/2b/kOB07dmTAgAEsW7aMTp06ceyxx3LWWWfRoUN823tTp07l5ptvZty4cXlfnz9/PlOmTOHoo48OnsbYsWMZOXJkSynSasT3CYqISPTWWmstZsyYwezZs5k0aRL33Xcfl1xyySrjLV++vObT/uyzyk5bGjp0aMGkDElivuWWW6pq09ixY/n444+reo9mSswiIlKVjTbaiPHjx3Pttdfi7kyYMIHDDz+cb37zm+y11148+uijHHDAAS3jn3766UyYMAGA++67j379+rHTTjtx5plnrjRert69e3PppZey0047cccdd/DKK6+wzz77MGTIEHbeeWdeeuklAO644w6++tWvMmjQIHbZZReAlaY/efJkmpqaaGpqYtttt2XJkiWMGTOGxx9/nKamJq6++mrmz5/PzjvvzODBgxk8eDBTpkxpeZ/ddtuNww47jH79+vHtb38bd2fcuHG8/fbb7L777uy+++5Vf57alS0iIlXr06cPK1as4N133wXgqaeeYubMmWywwQY8+uijeWOWLl3KySefzGOPPcbmm2/echepQjp37swTTzwBwB577MH1119P3759eeaZZzj11FP5+9//zqWXXsoDDzxAz549+eCDD1Z5jyuvvJJf/OIX7Ljjjnz00Ud07tyZn/zkJ1x55ZVMnJiU4/j444+ZNGkSnTt3Zu7cuRx11FFMnToVgOnTpzN79my+/OUvs+OOO/Lkk09y5plnctVVV/HII4/UpEKZtphFRKQmcks877nnnmywwQZFx3/ppZfo06cPm2++OUDJxNx828iPPvqIKVOmcPjhh9PU1MTJJ5/MwoVJWY0dd9yR4447jhtuuCHvLu8dd9yRs88+m3HjxvHBBx/QqdOq26fLli3jpJNOYsCAARx++OEr3V5y2LBh9OrViw4dOtDU1FSXAijaYhYRkaq9+uqrdOzYkY022ghY+RaPnTp1YsWKFS3Ply5dCqycyFvbe++9eeeddxg6dCi//vWvV3rPFStWsP766+e9o9T111/PM888w7333ktTU9Mq44wZM4b999+f++67j+23356HHnpolfe4+uqr2XjjjXn++edZsWIFnTt3bnltzTXXbPm/Y8eOdTmGri1mERGpyqJFi/je977H6aefnvdezJttthkvvvgin376KYsXL+bhhx8GoF+/frz66qstW53Nt3cEeOCBB5gxY0ZLUs617rrrsvnmm3PHHXcASYJ//vnnAXjllVcYPnw4l156Kd27d+fNN99cKfaVV15hwIABnHfeeQwdOpSXXnpplVtQLl68mB49etChQwd+97vflXWyWev3qIa2mEVEGlw5lzfV2ieffEJTU1PL5VLHHHMMZ599dt5xN9lkE4444ggGDhxI37592XbbbYHkzO7rrruOffbZh+7duzNsWPl3CP7DH/7AKaecwmWXXcayZcs48sgjGTRoEOeeey5z587F3dljjz0YNGgQkydPbokbO3YsjzzyCB07dmTrrbdm3333pUOHDnTq1IlBgwZx3HHHceqpp3LooYdyxx13sPvuu6+09V/IyJEj2XfffenRowePPPJI2f3Ip6LbPtbL0KFDvfnAukgWVCtbGll7uu3jRx99RJcuXXB3TjvtNPr27cvo0aPbulk1V5fbPoqIiNTaDTfcQFNTE9tssw2LFy/m5JNPbusmtTntyhYRkTYzevTodrmFXA1tMYuIiEREiVlERCQiSswiIiIR0TFmWS0d8cPCi/6sDNshItKaErOISKO7eL0av9/ikqOYGd/5znf43e9+ByR3kerRowfDhw9n4sSJTJgwgXPPPZdevXrx0Ucf0adPHy666CK+9rWvAXDccccxefJk1ltvPdydq666ij322KO2/WhQ2pUtIiIVW2eddXjhhRf45JNPAJg0aRI9e/ZcaZwRI0Ywffp05s6dy5gxYzjkkEOYM+fzOgFXXHEFM2bMYOzYsXzve9/LtP0xU2IWEZEg++67L/feey8At956a9GbUOy+++6MHDmS8ePHr/LaDjvswFtvvVW3djYaJWYREQly5JFHctttt7F06VJmzpzJ8OHDi44/ePDglvsm5/rb3/7Gt771rXo1s+HoGLOIiAQZOHAg8+fP59Zbb2W//fYrOX7rEtDnnnsuP/jBD3j33Xd5+umn69XMhqMtZhERCXbggQdyzjnnlLyXMsD06dNXqhd9xRVXMG/ePC677DKOPfbYejazoSgxi4hIsBNOOIELL7yQAQMGFB1v8uTJjB8/npNOOmml4R06dGDUqFGsWLGCBx54oJ5NbRjalS0i0ujKuLypXnr16sWoUaPyvnb77bfzxBNP8PHHH7P55ptz55135r0rlplxwQUX8LOf/Yy999673k2Onm77KKulATcV/nU/61iVGJG4tafbPq4udNtHERGRBqXELCIiEhElZhERkYgoMYuIiEREiVlERCQiSswiIiIR0XXMIiINrtjlfyHKvWTw8ssv55ZbbqFjx4506NCBX/3qVyXrZZfra1/7GlOmTGH+/PlMmTKFo48+GoCpU6dy8803M27cuJpMJ0ZKzCIiUrGnnnqKiRMn8txzz7Hmmmvyr3/9i3//+981e/8pU6YAMH/+fG655ZaWxDx06FCGDl3l0t92RbuyRUSkYgsXLqR79+6sueaaAHTv3p0vf/nLTJs2jV133ZUhQ4aw9957s3DhQgB22203zjvvPIYNG8ZXvvIVHn/8cQBmz57NsGHDaGpqYuDAgcydOxeALl26ADBmzBgef/xxmpqauPrqq3n00Uc54IADWLFiBb179+aDDz5oadOWW27JO++8w6JFizj00EPZbrvt2G677XjyySez/GiqpsQsIiIV22uvvXjzzTf5yle+wqmnnsrkyZNZtmwZZ5xxBn/605+YNm0aJ5xwAj/60Y9aYpYvX86zzz7L2LFjueSSSwC4/vrrGTVqFDNmzGDq1Kn06tVrpen85Cc/Yeedd2bGjBmMHj26ZXiHDh046KCDuOuuuwB45pln6N27NxtvvDGjRo1i9OjR/OMf/+DOO+/kxBNPzOATqR3tyhYRkYp16dKFadOm8fjjj/PII48wYsQILrjgAl544QX23HNPAD777DN69OjREnPIIYcAMGTIEObPnw/ADjvswOWXX86CBQs45JBD6Nu3b9ltGDFiBJdeeinHH388t912GyNGjADgoYce4sUXX2wZ78MPP2TJkiV07dq12m5nQolZRESCdOzYkd12243ddtuNAQMG8Itf/IJtttmGp556Ku/4zbu9O3bsyPLlywE4+uijGT58OPfeey977703v/71r/n6179e1vR32GEH5s2bx6JFi7j77ru54IILAFixYgVPPfUUa621Vg16mT3tyhYRkYq9/PLLLceDAWbMmEH//v1ZtGhRS2JetmwZs2fPLvo+r776Kn369OHMM8/kwAMPZObMmSu93rVrV5YsWZI31sw4+OCDOfvss+nfvz/dunUDkt3s11577UptayTaYhYRaXBtcUe0jz76iDPOOIMPPviATp06seWWWzJ+/HhGjhzJmWeeyeLFi1m+fDlnnXUW22yzTcH3uf322/n973/PGmuswZe+9CUuvPDClV4fOHAgnTp1YtCgQRx33HFsu+22K70+YsQItttuOyZMmNAybNy4cZx22mkMHDiQ5cuXs8suu3D99dfXtP/1pNs+ympJt32URqbbPjYe3fZRRESkQSkxi4iIRESJWUREJCJKzCIiIhFRYhYREYmIErOIiEhEdB2ziEiDm9OvtpdO9X9pTslx3nzzTXbZZRemTZvGBhtswPvvv8/gwYOZMGEC++yzD/369WPp0qV07dqV0047jWOPPRaACRMmcO6559KzZ0+WLl3KySefvFINbFFiFhGRAJtssgmnnHIKY8aMYfz48YwZM4aRI0ey2WabscUWWzB9+nQgqex1yCGHsGLFCo4//nggKQpy7bXX8t5777HVVltx2GGHsckmm7Rld6KiXdkiIhJk9OjRPP3004wdO5YnnniC73//+6uM06dPH6666irGjRu3ymvdunVjyy23bLk1pCS0xSwiIkHWWGMNrrjiCvbZZx8efPBBvvCFL+Qdb/Dgwbz00kurDH/jjTdYunQpAwcOrHdTG4q2mEVEJNj9999Pjx49eOGFFwqO07r08+23384222xDnz59GDVqFJ07d653MxuKErOIiASZMWMGkyZN4umnn+bqq68uuEt6+vTpK9WJHjFiBLNnz+bxxx/n+9//Pv/85z+zanJDUGIWEZGKuTunnHIKY8eOZdNNN+Xcc8/lnHPOWWW8+fPnc84553DGGWes8toOO+zAMcccwzXXXJNFkxuGjjGLiDS4ci5vqrUbbriBTTfdlD333BOAU089lQkTJvD666/zyiuvsO2227ZcLnXGGWe0nJHd2nnnncfgwYM5//zz6dq1a5ZdiJZu+yirJd32URqZbvvYeHTbRxERkQalxCwiIhIRJWYRkQYUw2FIKU+l80qJWUSkwXTu3Jn33ntPybkBuDvvvfdeRddq66xsEZEG06tXLxYsWMCiRYvauilShs6dO9OrV6+yx1diFhFpMGussQabb755WzdD6kS7skVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIiUTs5ltYmaPmNkcM5ttZqPS4RuY2SQzm5v+/WI63MxsnJnNM7OZZja43p0QERFpL8rZYl4OfN/d+wPbA6eZ2dbAGOBhd+8LPJw+B9gX6Js+RgK/rHmrRURE2qmSidndF7r7c+n/S4A5QE/gIOCmdLSbgG+l/x8E3OyJp4H1zaxHzVsuIiLSDlV0jNnMegPbAs8AG7v7QkiSN7BROlpP4M2csAXpsNbvNdLMpprZ1EWLFlXechERkXao7MRsZl2AO4Gz3P3DYqPmGearDHAf7+5D3X3ohhtuWG4zRERE2rVO5YxkZmuQJOU/uPuf08HvmFkPd1+Y7qp+Nx2+ANgkJ7wX8HatGixSC7Nee6OtmyAiklc5Z2Ub8BtgjrtflfPSPcCx6f/HAn/JGf4f6dnZ2wOLm3d5i4iISHHlbDHvCBwDzDKzGemw84GfAH80s+8CbwCHp6/dB+wHzAM+Bo6vaYtFRETasZKJ2d2fIP9xY4A98ozvwGlVtktERGS1pMpfIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhKRcm77KCINZk6//gVf6//SnAxbIiKV0haziIhIRJSYRUREIqLELCIiEhElZhERkYjo5C+RduiIHxZetWdl2A4RqZy2mEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISERKJmYz+62ZvWtmL+QMu9jM3jKzGeljv5zXfmhm88zsZTPbu14NFxERaY/K2WKeAOyTZ/jV7t6UPu4DMLOtgSOBbdKY68ysY60aKyIi0t6VTMzu/hjwf2W+30HAbe7+qbu/BswDhlXRPhERkdVKNceYTzezmemu7i+mw3oCb+aMsyAdtgozG2lmU81s6qJFi6pohoiISPsRmph/CWwBNAELgf9Jh1uecT3fG7j7eHcf6u5DN9xww8BmiIiItC9Bidnd33H3z9x9BXADn++uXgBskjNqL+Dt6pooIiKy+ghKzGbWI+fpwUDzGdv3AEea2ZpmtjnQF3i2uiaKiIisPjqVGsHMbgV2A7qb2QLgImA3M2si2U09HzgZwN1nm9kfgReB5cBp7v5ZfZouIiLS/pRMzO5+VJ7Bvyky/uXA5dU0SkREZHWlyl8iIiIRUWIWERGJiBKziIhIRJSYRUREIqLELCIiEhElZhERkYgoMYuIiESk5HXMItJ4Zr32Rls3QUQCaYtZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiESmZmM3st2b2rpm9kDNsAzObZGZz079fTIebmY0zs3lmNtPMBtez8SIiIu1NOVvME4B9Wg0bAzzs7n2Bh9PnAPsCfdPHSOCXtWmmiIjI6qFkYnb3x4D/azX4IOCm9P+bgG/lDL/ZE08D65tZj1o1VkREpL0LPca8sbsvBEj/bpQO7wm8mTPegnTYKsxspJlNNbOpixYtCmyGiIhI+1Lrk78szzDPN6K7j3f3oe4+dMMNN6xxM0RERBpTaGJ+p3kXdfr33XT4AmCTnPF6AW+HN09ERGT1EpqY7wGOTf8/FvhLzvD/SM/O3h5Y3LzLW0RERErrVGoEM7sV2A3obmYLgIuAnwB/NLPvAm8Ah6ej3wfsB8wDPgaOr0ObRURE2q2Sidndjyrw0h55xnXgtGobJSIisrpS5S8REZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESUmEVERCKixCwiIhIRJWYREZGIKDGLiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISEQ6tXUDRERiN6df/4Kv9X9pToYtkdWBtphFREQioi1mEZESjvhh4a/KWRm2Q1YP2mIWERGJiBKziIhIRJSYRUREIqLELCIiEhElZhERkYgoMYuIiEREiVlERCQiSswiIiIRUWIWERGJiBKziIhIRJSYRUREIqLELCIiEhElZhERkYgoMYuIiEREt30UidiAmwYUfG3WsbrhoEh7pC1mERGRiCgxi4iIRESJWUREJCJVHWM2s/nAEuAzYLm7DzWzDYDbgd7AfOAId3+/umaKiIisHmqxxby7uze5+9D0+RjgYXfvCzycPhcREZEy1GNX9kHATen/NwHfqsM0RERE2qVqE7MDD5rZNDMbmQ7b2N0XAqR/N8oXaGYjzWyqmU1dtGhRlc0QERFpH6q9jnlHd3/bzDYCJpnZS+UGuvt4YDzA0KFDvcp2iIiItAtVbTG7+9vp33eBu4BhwDtm1gMg/ftutY0UERFZXQQnZjNbx8y6Nv8P7AW8ANwDHJuOdizwl2obKSIisrqoZlf2xsBdZtb8Pre4+9/M7B/AH83su8AbwOHVN1NERGT1EJyY3f1VYFCe4e8Be1TTKBERkdWVKn+JiIhERIlZREQkIkrMIiIiEVFiFhERiYgSs4iISESqrfwlItLuzXrtjbZugqxGtMUsIiISESVmERGRiCgxi4iIRESJWUREJCJKzCIiIhFRYhYREYmILpcSiZgu0xFZ/WiLWUREJCJKzCIiIhFRYhYREYmIErOIiEhElJhFREQiosQsIiISESVmERGRiCgxi4iIRESJWUREJCJKzCIiIhFRYhYREYmIamVLFOb061/wtf4vzcmwJe1D76W3FHxtfnbNEJEASswoKYiISDy0K1tERCQi2mKWKBzxw8KL4qwM2yEi0ta0xSwiIhIRbTGjrTUREYmHErOslnTWsojESruyRUREIqLELCIiEhElZhERkYjoGLNEYdZrb7R1E2Q1MOCmAQVfm3WsTvWUOGiLWUREJCJKzCIiIhFRYhYREYmIErOIiEhElJhFREQiorOyG4TOJhURWT1oi1lERCQiSswiIiIRUWIWERGJiI4x0xhVpxqhjSIiUj1tMYuIiEREiVlERCQiSswiIiIRUWIWERGJiBKziIhIRJSYRUREIqLELCIiEhFdxywNTTXERaS90RaziIhIRLTFLCKrDVXQk0agxCx5zenXv6OanO0AACAASURBVOBr/V+ak2FLRERWL9qVLSIiEhFtMUteR/yw8KKhU6pEROpHW8wiIiIRUWIWERGJiHZlSxR6L72l4Gvzs2uGiEib0xaziIhIRLTFLHk1yvWejdJOEUmoWl9p2mIWERGJiBKziIhIRJSYRUREIqLELCIiEhGd/CUiUkIjXM6nk6raD20xi4iIRERbzCIRa4QttUaiz1MagRIzWlnzaZTPpFHaKSJSLu3KFhERiYi2mNs5nRAi7ZGW68alan2laYtZREQkInXbYjazfYBrgI7Ar939J/Wa1uog9Fjqkjn62EWkDi5er8hri7NrRwmNuHelLlvMZtYR+AWwL7A1cJSZbV2PaYmIiLQn9dpiHgbMc/dXAczsNuAg4MVygkN/4WT+yyjkF2OD/MoM1XvMvQVfm/+T/TNsSZ0Ezr9G+VxC2xm67jXK5xIqy++y0L1jjbJFmXU72zIPmbuXNWIlzOwwYB93PzF9fgww3N1PzxlnJDAyfboV8HKBt+sO/CugGYpTXHuLa4Q2Kk5xiis/bjN333CVoe5e8wdwOMlx5ebnxwA/D3yvqYpTnOIao42KU5ziqo+r11nZC4BNcp73At6u07RERETajXol5n8Afc1sczP7AnAkcE+dpiUiItJu1OXkL3dfbmanAw+QXC71W3efHfh24xWnOMVlPi3FKU5xbRRXl5O/REREJIwqf4mIiEREiVlERCQiSswiIiIRUWIWERGpATPbuBbvo8ScMTP7kpl9Kf1/QzM7xMy2CXif/6p968qa7vGBcV0qHH+DMsfrZ2Z7tH7/9CYqFau0nfVkZuua2RZ5hg8sEVeTZayCdg4zs+3S/7c2s7PNbL8q3q+iZayt1oVQpfrXBvMvaDmTvJ43s0lmdoKZFanfW0JIJZNGegAX1iMO2Bv4LtC71fATisScDLxGckOoU4BngN+SlCP9bpG4ca0ePwc+aH4e2L/7A+PeqHUcsCMwB5gNDAcmAa8CbwI7FIk7M/3s7k4/04NyXnsu4/6NL/Lal4BfktzYpRtwMTAL+CPQo0DMESRFeWakn8t25fQtdBkLXVaAi4CnganAfwN/By4EHgN+VIdlJdN1AVg37dfvgKNbvXZdHfpXj/k3q8hrocvZc8AFwBYVtmUo8Ajwe5IiVJOAxSS1L7atQ/+6AJemfVsMLEqX1+NCplXG9DqS5IYbgXfS76YRwFqVTKNut32MyIkkM6Zmcekv9J1IFs7zzWysu/88ffl0khUpn9OBbYC1gNeBLd39n2b2RZKF9TcF4g4BHgUeBCwddiQwrVgHzGxwoZeApiJxM4vEFdxVY2ZnF4krtiV6NckXRBfgXuBb7v5E2v6fkyTufE4Chrj7R2bWG/iTmfV292v4/HOqWTuLbMUbUGwLcQJJv9Yhmc9/APYnubHL9enf1s4n6dtCMxsG/M7Mznf3P1OkbwQuY6HLCnBY+vqawD+BXu7+oZldQZJULi8wvaBljIzXBZIv2LnAncAJZnYoSYL+FNi+yPRC+xc6/w4pMr0vFZle6HL2RWB94BEz+ydwK3C7u5eq8HgdyY+59YEpwGh339PM9khf2yFvJ8L79wfgLpJkeQTJOngbcIGZfcXdz6/l9Nz9M5L6HQ+kxbX2JVk+rzGzh93920Xa2qJdJGYz+7DQSyQLeE3jgG+S/LpbbmYXA7eYWR93H03xhXmZu38MfGxmr7j7PwHc/X0zK3ZBeX/gx8A+wLnu/paZXeTuNxWJgeRX6OQCbVq/SNzGJAvy+62GG8nKVMh/AVcAy/O8VuywyRruPgvAzBa5+xMA7v6cmRWbDx3d/aN03PlmthtJct6M4vMhtJ2LSL4sc9/b0+cbFYnbuPmHm5md6u4/TYf/3My+WyCmo7svBHD3Z81sd2CimfVKp1lI6DIWuqwsT7+Mmqf3YTq9T8xsRZG40GUs63VhC3c/NP3/bjP7EfB3MzuwxPRC+xc6/24nSUL5xulcJC50OXvf3c8BzjGznYGjgOfMbA5wq7sXKqqxhrvfD2BmP3X3P6XTftjMrqxD/3q7+4T0/6vM7B/u/uP0cMKLJD9Majm9Fu7+bzN7kWRv4BCSWyCXJ3RzPqYH8AbJl1++196sQ9ycVs87kvySvQOYXSRuKsmCCcmWRfPwzsDzZfRzCMmv5nOA+WWM/wLQN6B/vwF2KvDaLUXippD8+q50es/n/P+t1n0oEvd3oKnVsE7AzcBndWjnXGDTKvt3WavXZhZp4xathnUFHgY+rfUyVsWy8gywdvp/h5zh61F8V2jQMpYzTlbrwpzcfqXDjiXZNfp6rftXxfybBnw1oH+hy9kq85bke3Af4MYicU8Be5Hc6Oj15vUd2JUiN3uosn87pf9/E3gg57WXaz299PVNgXNJ9qi+DFwC9C+1jK70HpWMHOsDuAwYVuC1n9YhbiKwa4H3W1FihnXKM7wn8I0y+2rAacDvyxj3MGCrAq99q5zpVTgftgK6F3gt7w+g9LUDSb/cWw3fAvhBkbhewJcKvLZjHdp5GjCowGtnFIm7FOiSZ/iWwJ8KxAwi2Y3ZevgawLdrvYyFLivAmgWGdwcG1HoZazWNuq8LwM/yfW4kCWhuHfoUOv92pvCPxqFF4kKXs9sC+zeIZFfv/UA/4BqScwRml1hnQ/s3EHg2ncYTwFfS4RsCZ9ZhelNIfnD8T7HxSj1UkjNA8+5Vd/8kz2s93f2tMt5jM5Jf8A+Z2doku5SWVBi3FslKXDJORETqy8x2BR7zKhOrLpcK4O6f5EvK6WvlJOWTgD8Bv0oH9SQ5e6/SuF7lxBV4r0Inw5SKmxgYNzLjuKzbWewElGJxFc8HMwsqpl9FXNbLSmhc1v3L+nMJ7d8BGU8v688ztH81n567Twb2MbPHzOxfZrbIzCZXevlgu0/MZvZchHGnkZxp3HyizFyKnzxUbVw+pwTGnRQYV+xkrHrEZd3OQmfUlxIyH35VepSaxmW9rITGZd2/rD+X0P5tl/H0sv48Q/tX8+mlG08/Jrkksg/J4bhLgIsr+dGvXdltwMyecffhZjbd3bc1s04kJ1OUKhwRFCciUi4z28jd323rdjSi9Czsndz9/1oN7wY84e79y3mfdr/FHKnJZnY+sJaZ7UlyNvdf6xjXwqqskmRm9wfGhVYMC40r2k5r44phpeaDJdWffmlmvzCzbmZ2sZnNMrM/mlmPWseVaEu/wLiC88CSalP/bWa/M7OjW712XZG4mvUvdF0oJy6G/uW8555FXtug1aMb8KyZfdHKrL6X8151+zzT8WpSoazO07PWSRnA3d+roIntY4vZzDYhuS61J8nZfle4+7L0tbvd/VsxxOXEdyCpGrYXya7TB4BflzphoNI4MxvXehBwDMnlRLj7mQXiihVjmOjuFX9BmNkb7r5pLeNC22lmZ5IcFphDUlxilLv/JX3tOXcPOe5brJ0Vzwcz+xufFyU5muSayltJipF8w93zFSUJjquib6Hz4E6Sy8+eBk4AlpEW7ig2D6r4XELXhdC4TPtXTIn5t4LkLOJcvYAFSfe8T4G4rD/PI4CxwLskZ4wf5+7/SF8r9nlmPb1ngJHu/nyr4YOAG9x9WL641tpFgRGSSlt3kqwE3yXZsvxm+itls4jiAHD3FcAN6aNsAXFBVZIILMZg4RXDQqskhRaNyLRiGGHzIaQoSXBcni+wlpco/llmXrgj8HMJXRdC4zLtn5ndU+glkjKwhfwA+AZJsZbmIj+vufvmJdqZ9ecZWqEs6+l9H7jHzG5Mp+Ekx6SPBb5TYpqfC73OKqYHMKPV8++QXBe3BcWLHGQalzP+LGBmq8fjJKUpu9UqjqRQwFjgFqBnOuzVMtoXWozhHZIt0M1aPXoDb9chLrSdL7Z63gX4G3BV63nbarylJCd2XJTn8UGRuIrnAwFFSaqMWwKMJPkCaf34Vx3mQWjhjtD+ha4LoXFZ9+99kjKvu7Z67Aa8U6KtvUgOi12V9reen0to3KxWz3uQJL4zKf5dnen00nG/RFK74E7gzyTfGXnrLRR8j0pGjvWRLuydWw37BjAPWBhLXM64PyMpjD8gfVyePs4D/lqHuEqrJIUWYwitdhQaF9rOTCuGhcwHAoqSVBn3d+BrBV57rQ7zIKhwR2j/QuZBNXFZ94/kkNruBV57rMw+fpNkL+A/6/W5VPF5BlUoa6vpVfuo65tn9QBGk78S17bApFjicsZ7stAwit+5JCgufb3sKknt/UHGFcMaYT4AG5Cn+lp7fITOg1jnXY0/m7UoUIqyLT9PAiuUZT09oC/JTWuuSr9n7gc+Ap4n585dJaff1gvC6vhIZ9LwnOfDSHdhAdNrHZe+vhnpL/h05esa2PbBbf35NXI7azEfQvuW9WeSdTtjneeN2r9ypxe6TNfqO6mC/tR9eiRlP0eSbJm/RVITvDOwJ/BM2e+T5YxuiwdwQGxxJCcDzOLz+67OTJPsOsARdYg7ieQknVfS532BhwP7d0Ng3MSM47Ju58gyxqnJfKiib6FxF2c8vazjCt5Lu05xjdK/ktMLXaZr/J1Usn9ZTY+cc1SAeYVeKzmNkIY10gO4JNY4kjvwrB8wjYriSG6C/gVytqopseu7DvOhR5ZxbdDOkxthPgT27Ztt3YY69y/veQP1imuU/pX53kHLdC3XhXL6l9X0yDkxjFYnibV+XuzR7guMuPtFMcaZ2f7AycCZZnahmV1Yx7hP3f3fOe/RieL3Ws033aoKk3h6z9d6x7VhO8spX1jVfKh3AYdC3L2iIjY50w0tTFJRXA36V+rSmZrEZT3/zGwjqHs7Q5fpqr+TmpXZv6ym18/MZprZrJz/m59vVe402st1zAWZ2Z7uPqnI6+sCG7r7K62GD3T3QtfXBsel41wPrA3sDvya5MzWZ8voS1Acq1YMO5UiFcMKXZRvaXUrL3BRfom23+/u+9YyLoZ2mtnf3f3rZY5e9nwI7VstPxMz+193/0q54+fxIMltDGsWV8Xn8iWSy9pWABcCZwCHklzWNKrQD7Iq4jKdf7ZqlS4jqeC1LeSvRlXN9HJU9N1Sg7hVmNl4dy9Vh7rS78COwIkkJ3D9zd2fzHntAne/rEBoWSU3S2kXlb+KKVH1JrS6S1BcTvxMdx+Y87cL8Gd336tOcZVWDFvAqhflX0lyQgPuflOBuNAqUKFxWbez9Q8uA75CcjN0vHSt87LnQxV9C41bwudbEM1xawMfJ2G+boG4YoVJjq1DXGj/Mq2k1gbzL7SCV9D0cuKzqmJYqDyokZwA26vG0/s1yfL/LEmlsMnufnb6WlnVAa2KW/S2i8RcourN1919nQJxM4B9/fPqLjcD57v7ny29UUQt43Lim29G8TRJZZr3gBfcvW894iplZl1JLorfiKQi0Ftm9mqhlTsn7jMKV4Ha3t3XqnFc1u28h+TOXpcBn6TxjwM7Abh76y/GYFX0LTTu5yTnLpzr7u+kw17zEhWg0oT+feDTPC//j7t3r3FcaP9a1svWP9bNbIa7N9U4Luv5dw4BFbxCp5e1dJ19nZXXWU+f93T3L9R4ejObf2hbstv7OqA7cBTwdBnf8SeRnJ29gbtvYWZ9gevdfY9ypt9edmXvTFJ966NWw43krOVCOjbvinL3Z81sd2CimfWi+PGH0LhmE81sfZJ628+lMb+uV5wlxzdat2sxMJWkutBKBdbTX3VnmdkQ4Pdmdi/l3fBkDslJUHPztOHNWsdl3U53P9DMDgbGA1e6+z1mtqzchFzJfAjtWxVxZ6Qxt5rZ3cC1edqazz9IfhxOaf2CmV1c67gq5nnuODcXea0mcW0w/640s9uAq9Nl+CLKmH9VfJ5A5d8tVcS9Cuzh7m/kea9i3y2h02tJ9O6+HBhpyfk8f6d4+d1mp5HknmfS95hr6TH/crSXxPw08LEnN6leiZm9XCRuiZlt4elx4nQLeDfgbmCbOsSRjv/j9N87Lblpemd3X1yvOJKL3D8jKUsHSa1YSLb+JpBU/Mk3vWlm9nWS4zFPlDGdiym8Up9Rhzgg03bi7neZ2YPAj83sRHJW4DJUPB8C+hYcl8Z8AzidZI9C5zLCDiMpVZrv/YptrYXGNY9Taf/+YmZd3P0jd7+geaCZbQn8bx3iQtsZHOfuC4DDzeybwCSSXbF1m14q6LslIG4s8EVglcRMUmWt1u2camb7uPvfmge4+6Vm9jbwyzKm96m7/9ss2cC3Sk828zY6hT+GB+HVXaqqQkOywvwn6XWCJNfUlXP9c2hcNRXDMi0CUMW8zLyd6XLwvQrGD634FtS3aj4TktrA+7X1fI1tnjdK/wir4FXx9KpYpoO/kwI/y6yn9zOSG2G8RFJc5C7g8nLj293lUma2WfqrHzNbKz2Gkpe7P+/u81rHkexJKHTcOjgux40kx9Z2SJ8vIDluWa+4LmY2vPlJely8eXfM8kJB6XGSPwHNlwL1ItkrUDErfMJV1XFZt9MS3wEOcvfrzWzT9DMtpeL5ENq3KuLWNrP/BC519/vMrK+ZHVAqrsB7XVyvuEZZNgvEjS9jnKr65+6fuPsL5baziukFfbdUEbcKS86cr0s7m9cHM7shfV7u+jAGWERSEOpk4D7ggqIRuWr9S6EtH2RchaaKuKnp39yL3Z+vY1xoxbBaXpRft2pHWbeTZFfWL4A56fMvAv+ox3wI7VsVcbeT3ArwhfT5WlRQsajVewUVJiknrlGWzQJxWRfEqNs6VMV3S1Bcgfe6t47trNn6UMmjvRxjbhZ6wD3ruH9bcvq8A5jZFuQ/O7UmcZ5cyjXAzNYjORP/g5yX/1gktLrjJCu34aQ6xmXdzuHuPtjMpqcx75tZyWPNgfMhtG+hcVu4+wgzOypt8yfW/CYV8sDCJGXGNcqymS+u7IIYGfYvaHqh3y1VfCfle6/969VOAteH0JPimrW3xJz1l1ho3EUk9//dxMz+AOwIHFfHOCypGLYN0Lm5ve5+aYmwyVZlEQAz+y93P7+SmIC4rNu5zJICBM0/kDYkKUBRznQqnQ+hfQuNC/3R2MICC5NUGBf1smmBhUlyZN2/4OkFfrcEx+XEn+ru11Uwfsj0QteH0JPigPaXmLP+EguKc/dJZvYcsD3JJV2j3P1f9Yqz8IphY0guys89TlLw8ixru2pVmbQzxziSkzk2MrPLST7PksePAudDRX2rQVxFP/6sQGGS5uFeuFBIUFyO2JfNCXxemOQRksIk+5MUJrk+/VtM1st00PIS+t1SaZyZnd16EPBDM+sM4O5X1aOdhG8M7ejuO+Y8n2VmT7r7jpacn1JUuygw0swyqkJTbVwa25PkLMiWH0fu/lg94iywYlilLONqR1m3s9V79AP2SOMfdvc5ZcRkMh+qZWbd+PzH39PFfvxZeGGSoLhQWS+bFliYJFTW61DOdDOpYpj+YLsPmM3n/TuL5DIq3P2SerQzjS17fciJeZ7kjnPPpM+HkRzrH2TlFKFqT4m5UZjZT4ERJAtZ8y5Qd/cD6xQXWmmsouMklnG1o6zbmROfrzzgEndfViKu4vlQad9qEJfvDN7FwOueFFrIFzOEpOhNc2GSeeV8lqFxaWzUy6aZPe/ug9L/L/OVr4FuqSrV1v0LnV5OXCZVDM1sU+Aq4BWSO/h9XGH/QttZ8fqQxm0H/JbkzG8j2YV9Isl39/7uXvQ4ervald0GX2KhB/i/BWzl7hUdu6siLrTSWEXHSTzjakdZtzPHc8AmwPskK936wEIzexc4yQuf3BMyH7Iq4NDsOmAwyVmrBnw1/b+bmX3P3R9sHeBhhUmC41KxL5tVFSYh+2U6dHnJpIqhJxW/DjOzg4BJZnZ1GdOoRTsrXh/S9lZ3cpvX+bTvLB8kF3X/NzAgfVyePs4D/hpR3P1Al4D+BcW1eo81gfXKHLeawiRGctb67ytsX8VxWbeT5Bjh3jnP9yL5Nb898Ewt50No36qIuw3YJuf51iTXz/ehjMtECCxMUmlcoyyboY82WKarLsBRyXdLNXEkx4qvAB4L/Gwr+Q4MXh9Izin4AcnJfxcCF5bdxnovYFk+2uBLLDTuTmAeycX845ofZfQvNC60YtjzJJcGNT8fRnrdNDnXOxaI3YyMqlVl3U7S68nzDSu2sobMh9C+VRG3Svubh9W6b6vDslngfQbH1r8qlpdM53tOn7aq8DMPbWfo+nA9SW315rrls4DflN3ekAUr1kcbfImFxh2b71FG/0Ljgi6SJ/yi/KwLtmTdzgdJ9opslj5+QFKbuCPwXC3nQxV9q6agwi+BXdPHdSS73takSBGVKpaxdr1sFnivcgp+ZN2/TAtwVBF3IMltVl9LnzcB99RxeqHrw8xWf7sAD5a9jIQsWLE+qli42rx6TZ0/l6CKYTnjrgesX8H4mVaraoN2dgd+DkxP3+NaYMP0vVapoV6L+VBp36r4TNYiuR3jXSQnZZ1DsrXRgSKHUUL7trosm6GPrPpXxfSyrmI4LW1jbtzMOk4vdH14Jv37NPBlkkQ+t9zPtV2d/OUZV6GpNM7M/ujuRxQ4aQwvcKZmaFyO4KIRFnZRftYFWzJtpyeXSxS6C9W8IqFB8yGwb0Fx7v4J8D/po7XWt1XNlWk1u3Tchlg2c9pbUUGTjPuXdQGO0Ljl7r7YKi9GF1o1MXR9CD3ZDGhnZ2VDtl9iAXGj0r+V3hQgNK7ZRQRcJG/hF+VPtgwLtmTdTksqff2AdL43D3f3r5cIrXg+hPatiri+JCc0bs3KfSt1WUrQMhYaF/uyaVUW/GiDZTrrAhyhcS+Y2dFAx3RZPRNY5Z7etZpe6Prg4bfobXmDdvMg8IB71nFt9Nl0IzlL8ACge5kxQcdJSHbznATcQXLHmpNIr5mvU1zW7XyQpLDMHJLjTr8FflqP+VBF30LjniApnDKT5Pj5xSTXjdZlGWuvyybJnd9+D/wHn58PsojyzwvJun/Bx0Qznu9rk1z58g+Sy1IvJ0l69Zpe0PpAFSe3uXu7S8xZf4lVFAcsIbkusPVjCfBhreNavUdP4GvALs2PMmKqOk6S4XzPtJ3AtNz5nv4/uR7zIbRvVcQ1921WzrDHy4gbnOexBdCpTnFRL5tAV5KqVLcAPdNhr1YQn/UyHTy9kO+WauKq6GPId2Do+lDVXana267sT9K/H5vZl0mqu5RT4i+TOHcveG/oYkLjmlmBimFAqRKgQcdJChwLr2fBlkzbCTRX+FqYHsp4m+T+taWmFzIfMingkGOpJaVm55rZ6cBbJNWkSgkqxFBFXNTLpldf8CPrZTp0ekHfLVXE/ZXC/fuVuy+t5fQIXx+quktbuyrJackN3n9OsuvhF6QLl7v/Z0xxWTOzl4GBXnnFsNz3WJMyj5OY2c8oXEVoJ3fPW0UoNK4N2nkA8DhJ9a+fA+uS7N66p8T0qpoPlfQtNM6SUoJzSKqZ/Zikb1e4+9Ml4m4Dfuzus9PnWwPnpu/xZy9QGzo0ror+Zb5spl/IpwI7uHvJGxjkia97/6qYXtAyXUXcNSRXQNyaDhoB/JNki3Rddz+mxtMLXR+mkOSFJz25RewWwK3uPqysCddz10FbPsioCk21cRl9FqGVxkIvys+6YEum7cxyPlTRt+qOccE6FbYztBBD1gVNMl0209c3o/KCH1n3L3R6mVYxJE+1r+ZhwOxaTy8nvtL1YU+SMrOLSO4sNh/Yrdz4SnatRM/M1jaz/zSzGzz5ZbRRunUTVVwb+BiYYWa/MrNxzY8y4m4kuaRgh/T5AuCyMuK6mNnw5ieW3FmlS/q0YOH3KuIybaeZfcXMHjazF9LnA82s5G0fCZsPoX0LijOzHczsRZKtBMxskJmVc8/bl83sl2a2a/q4DvjfdOur2M09QuMaYtk0s5NITsL6VTqoF8n1sKVk3b/Q6YV+t4TGbWjJDS2AlptbdE+f/rvW0wtdH9x9EsnNMo4j2bof6u6Plopr1t6OMd9IcgF67sJ1BzAxsris3ZM+KhV6nORE4LeWXBrScmcVM1uH5NKDWsdl3c4bSHa3/iqd3kwzu4XSX2Qh8yG0b6FxY4G9m9vp7s+b2S5lxB1Hsrv2LJLP8gmSYgzLSC7BqXVcoyybp5EUHXombedcM6vnMcqs16HQ75bQuO8DT5jZKyT92xw4Ne1fsVtbhk4vdH2A5PKq90ny7NZmhpdxa19of4k56y+xqg7wZ8XD78UaelF+poVesm4nsLa7P9tqVhfbGmmeXsh8yLxwh7u/2apvn5URE1SIITSOxlk2Qwt+ZL1Mh04v6Lulirj7LLm2uB9JYn7JPz/ha2ytp5fGVrw+WPjJZkD7S8xZf4kFf/llwaqvGHYRYUUAsi70knU7/5XO6+b5fhiwsMg0qpkPoX0LjXvTzL4GuJl9gaSAw5xSQRZYiCE0jsZZNidbWEGTrPtX0fRCl+kafCdBcvx7K5LlZWC6JXpznaYXtD4QfoteoP2dlb0ncAHJSv4g6cJVat9+1nFZMbMe7r7QzDbL97q7v17Ge3QjuZ2hAU97Uo6yVEzeKkLu/t16xLVBO/sA40muiXyfpFb6d9x9foHxq5oPIX0LjTOz7sA1wDfSuAeBUV74MpvmuCdIvtyvJrl37/Ek3y8X1SMujY1+2bTkUpvvktwa1IAHSK7cKPnFm2X/Kp1e6DJdg3XhImA3ku/c+4B9gSfc/bA6TS90fbgfONzdi+31KRzfnhIzZPslVk1cozCzniRnlbbsXSl1nMTMZrr7wJy/XUgufdmrHnFZtzMnfh2ggyfXrNZNSN+qiQts4zR3H2Jms9x9QDrscXffuR5x6XgNsWyGynqZznJ5CZVu+Q4iuRnFIDPbmOSHTsnLwLJkZneStPNhcvaieolSrM3a265sCD/gnnVc3ZnZEvIfzzLA3X3dEvGhx0kyLfSSdTstOWP4UKA30KnU7sJq5kNo36qI25CkjGNvVv6CPqFYHOGFGILiGmjZDC1okvUyXdH0Qpfpar+TgE/cfYWZIaWV6AAAFUJJREFULTezdYF3gYKHPWrwHRi6PoSebAa0s8TcBl9iVR3grzevsmIY4cdJsq5WlXU7/0Ly5TqN8k6QqWY+hPYtNO4vJMVTHqKMk1xynEWyC/VMkkIMu5PUha5XXKMsm/dTuODHBJLd9/lk3b+Kphe6TNfgO2lq2r8bSNa/jyhys40aTC9ofajmZDNoZ7uyLfsqNFVX1IpZtcdJ0vfIolpVpu00sxfc/auh06pEaN+qiJvhZVTbKhK/jrv/v3rHNdCy+aS775hvWO7u+zxxWS/TVU8va2bWm6Ta18w6TqOi9cFqc3Jb+9piBl4F1qDyM6OzjmsUzRflV3ScxMzWJrnecFN3P8nMNjWznd296PXdoXFZtxOYYmYD3H1WifFqIahvVcRNNLP93P2+ShppZjsAvyEpZrGpmQ0CTnb3U+sRR+Msm13MbLi7P5O+T7kFP7JepkOXl0xZctzo20Afd7807d8wdy/nFpUhKl0fqr1FL9D+tpiDDrhnHdcozCzvLsVSu2nM7HaS3Uz/4e5fteSSsqdK/fKsIi7rdr4IbElyNvanfH68qqxfw5Woom+hcUuAdUj6tYzyj8U9Q3IG8D3uvm06rOSehSriGmXZ3I7ktqArFfwgOfy1v7vnvba4DfoXNL2smdkvSQ4bft3d+5vZF0nu6LddnaYXtD5Uq71tMWddhaaqA/yxq2KlzLRgSxu0c9/A6VUstG9VxAUfk/OAQgyhcQ20bIYWNMm0f7El4CKGe3JTiOkA7v6+JdcX10Wl64NVf3Ib0M4Scxt8iTXKwlyRGhwnyaRgS9btNLN13b35Pth1Fdq3KuL6uftLZjY43+vu/lyJJocWYqgorlGWzVZtLrvgRxss0zU5JpqhZWbWkc/7tyGfn3hbM6HrQzU/bFeafnvYld0GX2KNtjBXxKq/KD+Tgi1Zt9PMJrr7AWb2Gsl8z90ScS9drapsoX2rIu6G9JjkI/nD/Osl2htaiKGiuEZZNnPiKir40QbLdNVFiLJkZt8muRJmMElt7MOAC9z9jhpPp6r1oerpt5PEnPWXWEMtzG3BGqRgS9bTk7aX5bJpbVOYpF0v02bWj+RexwY87O7l7JlpKO0iMUtt1eI4iWVQrSrrdhbarZUTV2p3b9lC+1ZF3CHF2uPufy7R3qBCDJXGNcqymRPzjLsPN7OnSW4D+B7wgrv3LTB+1st0TY6JZindlb0xK/fvjRpPo6r1oVrt4hhzG3yJNdzCXIlqj5NYRgVb2qCd+e6A1NIcoGa7t0L7VsVnUqykoQOlvohCC5NUFNcoy2aOigp+ZN2/Wh0TzYqZnUFSW/0dkuXFSPpX68OH1a4PVdEWs9ScNUjBlqyn155ZYGGS0LhQbblsWmBBkwqn0a6XaTObR3JmdtFzFxpdh7ZugLRLzYVXsooLFTQ9M/tSq+c90i/dhtd6d72ZDU13jZYy0cz2C5hkaFyoTJdNM1vbzP4zPZnoU2AjMytZfMLMJrZ6/pCZ3V9GbNbrUNbeJCmHm4kq1ofqpqstZqk1a5CCLVW081533z/n+UPAFsCd7n5OPdqaleazUXOe30Sym/B/3X1EkbjQwiSZFnDIetm08IIfPdx9Yc7zLwM9gO3d/Re1bmejMLPfkNyL+V5W7t9VdZpe0PpQ9XSVmKXWLONqVaFqOT0zM2Brd59ddcMiZGZdvc63t8xC1summU1196FmNt0/r2z2vLsPKrfNlch6HcqaJfdjXoW7X5JxO+q6Pigxi9SAmXXxBroBQD7Nu+jd/Z/p2dI7Ay8X+7FhgYUYQuMajZlNIbm050lPKlZtAdzq7sMKjL8u8EOgF3C/u9+S89p1XrqGuNSJmf2Xu5+fybSUmKVWrEEKttRjemb2hrtvWpMGtgEzOxkYQ7Ir+afAcSRn9u4I/Mzdf1MgLqgQQ2hcqLZaNq3ygh93AnOBp4ETSHbvH+3un5rZc+6e94dM1utQ1sxsrLufZWZ/JX//Dqzx9Ma1HgQcA9ycTq+uhwaUmKVmrEEKtlTRzrMLvSXwI3ffoFZtzFr6hT4cWAt4Hdgy3XL+IvBIlmdO10NbLptWQcGP1mepm9mPgP2AA4FJRRJzuy56ZGZD3H2ame2a73V3n1zj6S0AHiX5MdVc4e9K4Jx0enU9NKDELFImM1tKcj1qvtv1jXb39TNuUs3kbo21Pgaae3w0T1xQIYbQuEZklRX8mANs4+4rcoYdC/wA6OLueROv1JaZdQV+DGwEnOvub5nZq17DsrvFtIsCIxIHa5CCLVVM7zngbnefluc9T6xhE9vCCjNbw92XAblnnHem+GWVoYUYMi3g0FbLplVemOSvJIVqHmoe4O43mdk7wM+LTKddFz0qtIu+Wa131acndp1lZkOA35vZvWR4ebG2mEXKZGZbAe/l2xVpZhu7+ztt0KyaMLNNgbfdfXmr4T2B/u7+UP5IKcbaecGPrBTaRd+snrvq0ysuTgV2cPfv1Gs6K01TiVlEqmVmg3PPpDazocBCd3+rHnGNwszuBw6vxRn7rT8ryUb6o6Cvuz+UXofeqd6XDqryl0gNmNnItm5DvZjZ+DJGO6XV8zNIqnrdXqe4RvExMMPMfmX2/9u7+5jJyruM498LMV3YpRAxKBQWKjbFKtsX6AKuNSAtoouFKJZqam0CKyKKf9TESEtBS0oA/5GmKaQrpVo0oUhK3Q20vgCNRF62QNlqs22CNGLbGEjZtgF528s/znnY52XmmTn3zJl55sz1SSbPPjPzm/vsvXfue8+Zc66jGxYehZ+1vK/mzvKxKOkzkj4p6edaam8bcDtwU/3U0cDn22hrSbvZY44YnaSLbd80+J2zZ+GM2MLaoiCGeQ80id6Wj0VJbwc2Aptt/2kL7T0GbAYeXBQQs9v2ieNua0m7WZgjYhQlwSSj1M2bSQZbxFLaf9vOR22/VdKBwCNtXxeeQ9kRDUg6QdKZqm54v/j5s6e1TeMg6SfrQ4KfkHS4pKsk7ZZ0m6QjV6m7GPh34AFJlwA7gHOAOyRdOO66WSHptvrnbkmPL3+sUnfDssfHgT8Y8RD4zFN184h7JH1W0jGS/knSXkkPS2rzGvv7JF0OHKQqLOZzVGfOtyp7zBFDknQZcCnwdeAtwB/bvrN+rW8q0yyQdDfVjQHWA78N3Ar8PXAu8E7b5/apKwomKa2bFaWBH9MOtlirJD1EdR/mw4DrqHIDbpd0JnC17dNaavcA4ELgLKp/jy8C293ywpmFOWJI9WJymu0fSjqO6qSQv7X9V6uFcMwCLb3JwpJ4Ua1yz2SVB5MU1XWdphxssVYNGJ+dGy8JGIkY3o8sXPZi+0lJpwO313tFWrVy7Vv8tdbfrPLacqXBJKV1M0GFgR/TDrZYw/5P0lnAoYAlnWf786oiOl9pq1H1DjbZC+yi2lN/po12szBHDO+7kt5i+zGAes/5HOBmoNWzNCfgTtV3yLL94YUnJf008I1V6n6deuKy/dSi5w8HPthC3UywfciI9V+R9EtUwRb/Np6tmmm/T3UIex/wy8Alkm4B/gdo81LFu6gW/oW7fL23/vl94BZWT7ArlkPZEUOSdDTwsu3v9nhti+37p7BZEdESSffb3tLruTYvm8ohkogh2X6q16Jcv9bZRVl97pk8RN0wwSRjq+u69EtvpeNzSBsknbKorc3AwhUZvW5mMxZZmCPGQNKOaW9Di0oTp0oDVzoZ1DIG6Zfe2kxEuwjYLum/JD0JbAe2SVoPXNNWozmUHTEGC5fHTHs7ImL8JB1KtV4+O4n2ssccMQZdXJQlfWyI95QGkxTVdV36ZXjDjM8xtbMVuBi4TNJHJH2k7TZzVnbEGEi6y/avTHs7SvVIlRLwO6oTzmxf1qf0FvYHk9xDFUyylSqY5Mb65zjruu4W0i8rjDA+R233RuBg4Ayqw9jnAw+10daSdnMoO2I4q5xkImCH7ZndoylNnBohmKSoruvSL71NKxFN0uO2Ny36uQG4w/ZZbbS3IHvMEcN7GLiP3mEih014W8btZ6gSp85mf+LUlUNMeKXBJKV1XZd+6a10fI7q+frnc5KOAp4BXt9ym1mYIxr4OnCx7W8uf0HSf09he8ZmhMSp0mCS0rquS7/0MMVEtB2SDgOuBx6hCsXZ3najOZQdMSRJ5wO7be/p8dp5tlu/gfokSBJV4tRptt837e2JWGxa41PSa4B1tve23lYW5ohYrM7+foPtf5Z0EHBgvcfS9HPeZvuRSdV1XfqlMq7xOWRbB1NFxG60vU3SG4A32m41t2Cev7OIGJuW04cmRtI2qrtmLYRZHA2UHgkoDX5oMzBils19v4x5fA7j08ALwMJtJZ8Crm6xPSALc8S4dGXSvBTYQhXST/19+hElH2R72yTrui79AoxxfA7peNvXAS/V7T3PBO4kl4U5Ygw6NGm+YPvFhV8kHUjv2xf2VRr8MKnAiFmTflli5PHZ0Iv14XLX7R1PtQfdqpyVHTECSR+zffm0t2OM7pN0OXCQpHdRnWTzj/3eXBr8MK3AiLUu/TJQo/E5BlcCdwPHSLqVam/9Ay22B+Tkr4ih9Zs0qa837cKkKekA4ELgLKq/3xeB7e4zUYwQTDKVwIi1Lv2yuqbjc0xtHg6cWrf3gO2n22rr1TazMEcMJ5PmSpIOoQp+OIL9wQ9P2P6pNuq6Lv2y9kh6HXAsi44w2/5yq21mYY4YzjxMmpJ2s/I7u73ALuBq28/0qTuJ6j8pO4E/tH3ckO0V1XVd+qW30vE5QnvXAhcA/wHsq5+27XePs50V7WZhjmimy5OmpOuAV4C/q596b/3z+8Av2P61VWqLgh8SaNJb+mWlUcZnYXt7gE22Wz/ha0m7WZgjmuvqpCnpfttbej0nabftE1epLQp+mGRgxCxJv6w0yvgsbO8u4Ddt/3CcnztILpeKKLMR2GP7fZIOqg9zd8EGSacs/CJpM7Ch/vXlfkWlwQ9TCIyYCemXvorG5wieAx6TdJOkGxYeLbSzRC6XimionjR/D/gx4HiqSfNG4MxpbteYXATcXF+eI6pDhBdJWg9cs0rdpcBm4EGogh8kDRP8UFrXdemX3krHZ6kv1I+JysIc0VxnJ03bDwMnSjqU6quuZxe9fNsqpS/YfrE6wt8o+KG0ruvSLz2MMD5L25vKlRZZmCOa6/SkKWkr8LPAuoW/o+2/GFBWGvww6cCIWZF+6aNwfDZt4zbb7+lzFji2N42zvRXt5+SviGbqM0OfBd4P/BHVpPmftj801Q0bA0k3AgcDZ1Ddd/Z84CHbFw6oKwp+mEZgxCxIv/RWOj4L2jnS9nfqE/BWsP2tcba3ov05/3eOaKzLk6akx21vWvRzA3CH7bOmvW0R8zI+cyg7oiHb+4BP1Y+ueb7++Zyko4BngNcPKhohmGSigRGzIv3SV9H4bErSD+j99ZSoAkZeO+42F8vCHNFQxyfNHZIOA64HHqH6e24fou4u+gc/3AL0C34oreu69EtvpeOzEdtTvfwxh7IjGpp0+tC0SHoNsM723iHeWxT8MOnAiFmRfhmsyficNQkYiWhui+0/s727fnwION32tcBxU962kUg6WNIVkj5VxxAeIemcIUpLgx8mHRgxK9IvPYwwPmdKDmVHNLdB0im2H4TOTZqfBr4CnFb//hTwOWDHgLrS4IdJB0bMivRLb6Xjc6bkUHZEQ5LeDtxMtRi/OmlS3YFmq+2xBx1MiqRdtk+W9Kjtt9bPfdX2m4es7xX80Fpd16Vflhp1fM6K7DFHNDTp9KEJe7G+YYIBJB0PDHVnndLgh0kERsyi9EtPxeNzlmRhjijQ4UnzSuBu4BhJtwJbgA8MKuoX/NBWXdelX/oqGp+zJoeyIxqaVPrQtEg6HDiV6jD9A7afHqKmKPhhXgIjmkq/9FcyPmdNzsqOaO7nbb8f+J7tP6c6EeWYKW/TOK0Dvkf13fmbJP3iEDXLgx9eYrjgh9K6rku/9FcyPmdKDmVHNDeR9KFpkHQtcAHViWz76qcNfHlAaWnww0QCI2ZQ+qWHEcbnTMmh7IiGJF0BfJzq/sufoJ40bV8x1Q0bA0l7gE31NaKln1EU/NDlwIhRpF/2G8f4nAU5lB3RkO2P2n7W9j8AxwIndGFRrj0B/GjTotLgh3kJjGgq/dJX0ficNTmUHdGQpIOBDwIbbW+TtFHSO2x3IeTgOeAxSf/CostQbF82oK40+GEuAiMKpF96Kx2fMyULc0RzXZ40v1A/mjre9gWSfgvA9vNauI6snbquS7/0Vjo+Z0oW5ojmOjtp2v5MYWlp8MNcBEYUSL/0MML4nClZmCOa69ykKek22+9R71taYnvTgI8oDX6Yi8CIAumXRcYwPmdKzsqOaEjSu4APA28CvkQ9adq+d5rbNQpJR9r+jqRje71u+1tDfEZR8MM8BEaUSL/sN47xOUuyMEcUyKS5kqTXUZ2l/uqRONsDry8treu69Mv8yqHsiDIL6UMHUqUPzfSkKekH9DhESPUfD9t+7YD6ouCHeQmMaCr9stSo43PWZI85oqF+k6btd09vq6arNPhhXgIjmkq/zLfsMUc0dx7wxkyaSywEPzTtk9K6rku/zLEszBHNZdJcqTT4YS4CIwqkX+ZYFuaI5jJprlQa/DAXgREF0i9zLN8xRzQk6Xd7PT8v4QcR0a4szBFRrDT4Yd4CI4aVfgnIwhwxtEyaK5UGP8xbYMSw0i8BWZgjhpZJMyImIQtzRBQrDX6Yt8CIYaVfArIwRwwtk2ZETEIW5oiIiDXkgGlvQEREROyXhTkiImINycIcMWckXSXpT3o8f5ykr9V/PlnSDZPfuohIJGdErGB7F7Br2tsRMY+yxxzRAZLWS9op6auSvibpAklPSvrx+vWTJd27qOTNkv5V0jclbevxeadL2lH/+SpJN0u6V9ITkuY5EzyiddljjuiGs4Fv294KIOlQ4NpV3r8JOBVYDzwqaeeAzz8BOAM4BNgj6ZO2Xxp9syNiuewxR3TDbuCdkq6V9A7bewe8/07bz9t+GrgH2Dzg/Tttv1C//3+BnxjDNkdED9ljjugA29+QdBLwq8A1kr4EvMz+/3yvW14y4PflFt97+hUyd0S0JnvMER0g6SjgOdufBf4SeBvwJHBS/ZbfWFZyrqR1kg4HTgcentCmRsQA+V9vRDecCFwvaR/wEnAJcBDw15IuBx5c9v6HgJ3ARuCjtr8t6bjJbW5E9JNIzoiIiDUkh7IjIiLWkCzMERERa0gW5oiIiDUkC3NERMQakoU5IiJiDcnCHBERsYZkYY6IiFhD/h88Q7hO9aF7vQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 576x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Stacked Column Chart: visual form of the two-way table\n",
    "\n",
    "drtype_sublin_table.plot(kind=\"bar\", \n",
    "                 figsize=(8,8),\n",
    "                 stacked=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}