{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plot somme basis plots with `python` and `matplot` library"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Put the correct file name of your csv file generated previously with the python script "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "\n",
    "csv_file=\"/ifb/data/mydatalocal/shortReads.statsMerged.csv\"\n",
    "\n",
    "# Load the csv file into a python/panda dataframe\n",
    "df1 = pd.read_csv(csv_file, index_col=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Plot the reads number by fastq file "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABKUAAAF+CAYAAABTfKkAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABDq0lEQVR4nO3dd7gkZZn38e89gRwkjIozwACCEiQOMIASRAEz5sEAriiKYFjDLrDugqCuGRcDiEpQVJIJF5GgIuqSZhAlL7iojPAKAiKiIOF+/6jnMDXHMzNnznRXd1d/P9fV1+l+urp+dXfNOdPnPk9VRWYiSZIkSZIkNWlSrzdAkiRJkiRJw8emlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqUkSVIrRMSjEXF1RFwbEd+LiCd0aL2/iYi1O7GuceadEhGv6OD6/tKpdXVTRLwhIj7b6+2QJEnNsSklSZLa4m+ZuXVmbgHcAxzS6w1qWkRM7vU2SJIkjZdNKUmS1EaXAtMBImKjiPhBRMyLiJ9GxNPL+Isi4vKI+EVEXBQRTyrja0XEBWX8C0CMFRARf4mID0XELyPistrrF5rpNDJTKSJ2j4ifRMSZEfG/EfGRiHhtRFwREddExEa11T+nbOv/RsQLy+snR8THI+LKiPhVRLyltt4fR8TXgWsWsa2fjIirIuKHETGtvCdX1Z7fOCLmjfG6d0TE9SXv9DK2Q0T8T3l//icinlbG3xAR3ymz1G6NiEMj4t1lucsiYs2y3MUR8eny2msjYocxcqdFxDdLrVdGxC5j72ZJkjTIbEpJkqRWKbOF9gTOKUMnAm/PzO2A9wKfL+M/A2Zn5jbA6cC/lPEjgZ+V8XOA9RYRtTJwWWZuBVwCvHkcm7cV8E7gGcDrgU0ycwfgS8Dba8vNBHYDXgCcEBErAAcC92Xm9sD2wJsjYoOy/A7Av2XmZovYzqsyc1vgJ8CRmflr4L6I2Los80/AKWO89jBgm8zcEnhrGbsR2LW8P/8BfLi2/BbAa8r2fAj4a1nuUmD/+jZl5s7A24CTxsj9L+DYUuvLy/sjSZJaZkqvN0CSJKlDVoyIq6kaOvOACyNiFWBn4KyIxyc8LV++zgDOiIh1gOWAW8v4rsDLADLz3Ii4dxF5fwf+u9yfBzx3HNt4ZWbeARARvwYuKOPXAHvUljszMx8Dbo6I/wOeDuwFbFmbhbU6sHHZjisy81bG9hhwRrl/GvCtcv9LwD9FxLuBV1M1kkb7FfC1iPgO8J1a7qkRsTGQwNTa8j/OzPuB+yPiPuB7tfq2rC33DYDMvCQiVhvj/F/PATar7bPVImLVsm5JktQSzpSSJElt8bfM3BpYn6rJdAjVZ50/lXNNjdw2Lct/BvhsZj4DeAuwQm1dOY68hzNzZLlHWfDHvkdKLlF1VZarveah2v3Hao8fY+E/Fo7OT6rDCN9eq2ODzBxpaj0wju0dve5vAs8DXgjMy8y7x1j2BcDngO2AeRExBTiGqvm0BfAiFn7flqW+uknATrVap9uQkiSpfWxKSZKkVsnM+4B3UB2q9zfg1oh4JVRNoojYqiy6OvD7cv+A2iouAV5bln8esMZSbsJvqJo4AC9h4ZlE4/XKiJhUzjO1IXATcD5wcERMLdu2SUSsPI51TQJGZle9huqwRTLzwbLO44GTR78oIiYB62bmj6kObXwCsAoLv29vmEBtUM3MIiKeSXVI4n2jnr8AOLS2LVtPMEeSJPUxm1KSJKl1MvMXwC+BOVQNpgMj4pfAdVSNIoCjqA7r+ynwx9rLPwDsWk4Evhfwu6WM/yKwW0RcAezI0s1iGnET1fmfzgPeWhpIXwKuB66KiGuBLzC+UzE8AGxeTmT+bODo2nNfo5qldMEYr5sMnBYR1wC/oDrH05+AjwH/GRE/L8tMxL0R8T/ACVTnyhrtHcCscoL161lwPitJktQisWDWuSRJkoZJRLwXWD0z/73BzIuB92bm3KYyJUlSf/JE55IkSUMoIr4NbEQ1e0qSJKlxzpSSJEmSJElS4zynlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxnmi82LttdfOmTNn9nozJEmSJEmSWmPevHl/zMxpYz1nU6qYOXMmc+d6ZWJJkiRJkqROiYjfLuo5D9+TJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4zym1GA8//DDz58/nwQcf7PWmNG6FFVZgxowZTJ06tdebIkmSJEmSWsim1GLMnz+fVVddlZkzZxIRvd6cxmQmd999N/Pnz2eDDTbo9eZIkiRJkqQW8vC9xXjwwQdZa621hqohBRARrLXWWkM5Q0ySJEmSJDXDptQSDFtDasSw1i1JkiRJkpphU2qInHLKKRx66KG93gxJkiRJkiTPKbU0Oj15KHNpl08yk0mT7CVKkiRJkqTBZnejz/3mN79h00035W1vexvbbrstxxxzDNtvvz1bbrklRx555OPL7bvvvmy33XZsvvnmnHjiiY+Pn3zyyWyyySbstttu/PznP398/KyzzmKLLbZgq622Ytddd220JkmSJEmSJGdKDYCbbrqJk08+mX333Zezzz6bK664gszkxS9+MZdccgm77rorJ510EmuuuSZ/+9vf2H777Xn5y1/O3//+d4488kjmzZvH6quvzh577ME222wDwNFHH83555/P9OnT+dOf/tTbAiVJkiRJ0tBxptQAWH/99Zk9ezYXXHABF1xwAdtssw3bbrstN954IzfffDMAxx13HFtttRWzZ8/mtttu4+abb+byyy9n9913Z9q0aSy33HK8+tWvfnydu+yyC294wxv44he/yKOPPtqr0iRJkiRJWiYRE7sNSl6bOVNqAKy88spAdU6pww8/nLe85S0LPX/xxRdz0UUXcemll7LSSiux++678+CDDwKLvoreCSecwOWXX865557L1ltvzdVXX81aa63V3UIkSZIkSZIKZ0oNkL333puTTjqJv/zlLwD8/ve/58477+S+++5jjTXWYKWVVuLGG2/ksssuA2DHHXfk4osv5u677+bhhx/mrLPOenxdv/71r9lxxx05+uijWXvttbntttt6UpMkSZIkSRpOzpQaIHvttRc33HADO+20EwCrrLIKp512Gvvssw8nnHACW265JU972tOYPXs2AOussw5HHXUUO+20E+ussw7bbrvt44fqve997+Pmm28mM9lzzz3ZaqutelaXJEmSJEkaPpGZvd6GvjBr1qycO3fuQmM33HADm266aY+2qPeGvX5JkiRJUv+b6PmaJtoOaTpv0EXEvMycNdZzHr4nSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqWWYFjPuTWsdUuSJEmSpGbYlFqMFVZYgbvvvnvoGjSZyd13380KK6zQ602RJEmSJEktNaXXG9DPZsyYwfz587nrrrt6vSmNW2GFFZgxY0avN0OSJEmSJLWUTanFmDp1KhtssEGvN0OSJEmSJKl1PHxPkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcV1rSkXEuhHx44i4ISKui4h3lvGjIuL3EXF1uT2/9prDI+KWiLgpIvaujW8XEdeU546LiCjjy0fEGWX88oiYWXvNARFxc7kd0K06JUmSJEmStPS6efW9R4D3ZOZVEbEqMC8iLizPHZuZn6gvHBGbAXOAzYGnABdFxCaZ+ShwPHAQcBnwfWAf4DzgQODezHxqRMwBPgq8OiLWBI4EZgFZss/JzHu7WK8kSZIkSZLGqWszpTLzjsy8qty/H7gBmL6Yl7wEOD0zH8rMW4FbgB0iYh1gtcy8NDMT+Aqwb+01p5b7ZwN7lllUewMXZuY9pRF1IVUjS5IkSZIkSX2gkXNKlcPqtgEuL0OHRsSvIuKkiFijjE0Hbqu9bH4Zm17ujx5f6DWZ+QhwH7DWYtYlSZIkSZKkPtD1plRErAJ8E3hXZv6Z6lC8jYCtgTuAT44sOsbLczHjE31NfdsOioi5ETH3rrvuWlwZkiRJkiRJ6qCuNqUiYipVQ+prmfktgMz8Q2Y+mpmPAV8EdiiLzwfWrb18BnB7GZ8xxvhCr4mIKcDqwD2LWddCMvPEzJyVmbOmTZu2LKVKkiRJkiRpKXTz6nsBfBm4ITM/VRtfp7bYS4Fry/1zgDnlinobABsDV2TmHcD9ETG7rHN/4Lu114xcWe8VwI/KeafOB/aKiDXK4YF7lTFJkiRJkiT1gW5efW8X4PXANRFxdRk7AtgvIramOpzuN8BbADLzuog4E7ie6sp9h5Qr7wEcDJwCrEh11b3zyviXga9GxC1UM6TmlHXdExHHAFeW5Y7OzHu6UqUkSZIkSZKWWlQTizRr1qycO3durzdDkiRJkiQthRjrrNLjMNF2SNN5gy4i5mXmrLGea+Tqe5IkSZIkSVKdTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqUkSZIkSZLUOJtSkiRJkiRJapxNKUmSJEmSJDXOppQkSZIkSZIaZ1NKkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcTalJEmSJEmS1DibUpIkSZIkSWqcTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxXWtKRcS6EfHjiLghIq6LiHeW8TUj4sKIuLl8XaP2msMj4paIuCki9q6NbxcR15TnjouIKOPLR8QZZfzyiJhZe80BJePmiDigW3VKkiRJktTPIiZ2G5Q8Da5uzpR6BHhPZm4KzAYOiYjNgMOAH2bmxsAPy2PKc3OAzYF9gM9HxOSyruOBg4CNy22fMn4gcG9mPhU4FvhoWdeawJHAjsAOwJH15pckSZIkSZJ6q2tNqcy8IzOvKvfvB24ApgMvAU4ti50K7FvuvwQ4PTMfysxbgVuAHSJiHWC1zLw0MxP4yqjXjKzrbGDPMotqb+DCzLwnM+8FLmRBI0uSJEmSJEk91sg5pcphddsAlwNPysw7oGpcAU8si00Hbqu9bH4Zm17ujx5f6DWZ+QhwH7DWYtYlSZIkSZKkPtD1plRErAJ8E3hXZv55cYuOMZaLGZ/oa+rbdlBEzI2IuXfddddiNk2SJEmSJEmd1NWmVERMpWpIfS0zv1WG/1AOyaN8vbOMzwfWrb18BnB7GZ8xxvhCr4mIKcDqwD2LWddCMvPEzJyVmbOmTZs20TIlSZIkSZK0lLp59b0AvgzckJmfqj11DjByNbwDgO/WxueUK+ptQHVC8yvKIX73R8Tsss79R71mZF2vAH5Uzjt1PrBXRKxRTnC+VxmTJEmSJElSH5jSxXXvArweuCYiri5jRwAfAc6MiAOB3wGvBMjM6yLiTOB6qiv3HZKZj5bXHQycAqwInFduUDW9vhoRt1DNkJpT1nVPRBwDXFmWOzoz7+lSnZIkSZIkSVpKUU0s0qxZs3Lu3Lm93gxJkiRJkjoqxjrr8jhMtF1gXmfzBl1EzMvMWWM918jV9yRJkiRJkqQ6m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqUkSZIkSZLUOJtSkiRJkiRJapxNKUmSJEmSJDXOppQkSZIkSZIaZ1NKkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcYttSkXE5Ig4ramNkSRJkiRJ0nBYbFMqMx8FpkXEcg1tjyRJkiRJkobAlHEs8xvg5xFxDvDAyGBmfqpbGyVJkiRJkqR2G09T6vZymwSs2t3NkSRJkiRJ0jBYYlMqMz8AEBErZ+YDS1pekiRJkiRJWpIlXn0vInaKiOuBG8rjrSLi813fMkmSJEmSJLXWEptSwKeBvYG7ATLzl8CuXdwmSZIkSZIktdx4mlJk5m2jhh7twrZIkiRJkiRpSIznROe3RcTOQEbEcsA7KIfySZK0KBETe13mYORJkiRJWjbjmSn1VuAQYDrwe2Dr8liSJEmSJEmakPFcfe+PwGsb2BZJkiRJkiQNifFcfW/DiPheRNwVEXdGxHcjYsMmNk6SJEmSJEntNJ7D974OnAmsAzwFOAv4Rjc3SpIkSZIkSe02nqZUZOZXM/ORcjsN8LSwkiRJkiRJmrBFnlMqItYsd38cEYcBp1M1o14NnNvAtkmSJEmSJKmlFnei83lUTaiRi2y/pfZcAsd0a6MkSZIkSZLUbotsSmXmBk1uiCRJkiRJkobH4mZKARARk4EXADPry2fmp7q3WZIkSZIkSWqzJTalgO8BDwLXAI91d3MkSZIkSZI0DMbTlJqRmVt2fUskachELHmZseQEr3/adJ4kSZKkZTcIvzdMNGvSOJY5LyL2mtjqJUmSJEmSpH80nplSlwHfjohJwMNUV+PLzFytq1smSZLUJW3+i6MkaekNwv8Ly5In9avxNKU+CewEXJPpt4AkSZIkSZKW3XgO37sZuHZpG1IRcVJE3BkR19bGjoqI30fE1eX2/Npzh0fELRFxU0TsXRvfLiKuKc8dF1H1lCNi+Yg4o4xfHhEza685ICJuLrcDlma7JUmSJEmS1H3jmSl1B3BxRJwHPDQymJmfWsLrTgE+C3xl1PixmfmJ+kBEbAbMATYHngJcFBGbZOajwPHAQVSHEX4f2Ac4DzgQuDcznxoRc4CPAq+OiDWBI4FZQALzIuKczLx3HLVKkiRJkiSpAeOZKXUr8ENgOWDV2m2xMvMS4J5xbsdLgNMz86HMvBW4BdghItYBVsvMS8tMra8A+9Zec2q5fzawZ5lFtTdwYWbeUxpRF1I1siRJkiRJktQnljhTKjM/0OHMQyNif2Au8J7SOJpONRNqxPwy9nC5P3qc8vW2so2PRMR9wFr18TFes5CIOIhqFhbrrbfeslUlSRoqnhBVkiRJWjZLnCkVET+OiB+Nvk0w73hgI2BrqsMCPzkSM8ayuZjxib5m4cHMEzNzVmbOmjZt2mI2W5IkSZKGV8TEboOSJ6k3xnNOqffW7q8AvBx4ZCJhmfmHkfsR8UXgv8vD+cC6tUVnALeX8RljjNdfMz8ipgCrUx0uOB/YfdRrLp7I9i7Y1om9zr+GS5IGlf/3SZIkqduWOFMqM+fVbj/PzHcDO04krJwjasRLgZEr850DzClX1NsA2Bi4IjPvAO6PiNnlfFH7A9+tvWbkynqvAH5Uzjt1PrBXRKwREWsAe5UxSQNoEP4qJ0mSJElaekucKVWuZjdiErAd8ORxvO4bVDOW1o6I+VRXxNs9IramOpzuN8BbADLzuog4E7ieahbWIeXKewAHU13Jb0Wqq+6dV8a/DHw1Im6hmiE1p6zrnog4BriyLHd0Zo73hOtDyb+GDzb3nyRJkiRpEEUu4TfFiLiVBedqeoTqanxHZ+bPur95zZk1a1bOnTt3zOfa/kt/2/Paru37bxDy2lybeeaZ17s8/98b2yDsu2XJk/pV27/3zDPPvGXLW1xWRMzLzFljPTeeq+9tsPSbI0mSJEmSJC3aeE50TkTsDMysL5+ZX+nSNkmSJEmSJKnlxnNOqa8CGwFXAyPneUrAppQkSZIkSZImZDwzpWYBm+WSTj4lSZIkSZIkjdOkcSxzLeO42p4kSZIkSZI0XuOZKbU2cH1EXAE8NDKYmS/u2lZJkiRJkiSp1cbTlDqq2xshSZIkSZKk4bLEplRm/qSJDZEkSZIkSdLwGM85pSRJkiRJkqSOsiklSZIkSZKkxo3nnFKPi4g1gHUz81dd2h5JkiQtg4iJvS5zMPLaru37z7ze5/m9J6mfLHGmVERcHBGrRcSawC+BkyPiU93fNEmSJEmSJLXVeA7fWz0z/wy8DDg5M7cDntPdzZI6J2Jit0HJkyRJkiRpEI2nKTUlItYBXgX8d5e3R5IkSZIkSUNgPE2po4HzgVsy88qI2BC4ububJUmSJEmSpDZb4onOM/Ms4Kza4/8DXt7NjZIkSZIkSVK7LbIpFRGfARZ5bYbMfEdXtkiSJEmSJEmtt7jD9+YC84AVgG2pDtm7GdgaeLTrWyZJkiRJkqTWWuRMqcw8FSAi3gDskZkPl8cnABc0snWSJEmSJElqpfGc6PwpwKq1x6uUMUmSJEmSJGlClniic+AjwC8i4sfl8W7AUV3bIkmSJEmSJLXeeK6+d3JEnAfsWIYOy8z/193NkiRJkiRJUpuN5/A9gIeAO4B7gU0iYtfubZIkSZIkSZLabokzpSLiTcA7gRnA1cBs4FLg2V3dMkmSJEmSJLXWeGZKvRPYHvhtZu4BbAPc1dWtkiRJkiRJUquNpyn1YGY+CBARy2fmjcDTurtZkiRJkiRJarPxXH1vfkQ8AfgOcGFE3Avc3s2NkiRJkiRJUruN5+p7Ly13j4qIHwOrAz/o6lZJkiRJkiSp1cYzU4qIeCawcWaeHBHTgOnArV3dMkmSJEmSJLXWEs8pFRFHAv8KHF6GpgKndXOjJEmSJEmS1G7jOdH5S4EXAw8AZObtwKrd3ChJkiRJkiS123iaUn/PzAQSICJW7u4mSZIkSZIkqe3G05Q6MyK+ADwhIt4MXAR8sbubJUmSJEmSpDZb7InOIyKAM4CnA38Gngb8R2Ze2MC2SZIkSZIkqaUW25TKzIyI72TmdoCNKEmSJEmSJHXEeA7fuywitu/6lkiSJEmSJGloLHamVLEH8JaI+C3VFfiCahLVll3dMkmSJEmSJLXWeJpSz+v6VkiSJEmSJGmoLLEplZm/bWJDJEmSJEmSNDzGc06pCYmIkyLizoi4tja2ZkRcGBE3l69r1J47PCJuiYibImLv2vh2EXFNee64ckVAImL5iDijjF8eETNrrzmgZNwcEQd0q0ZJkiRJkiRNTNeaUsApwD6jxg4DfpiZGwM/LI+JiM2AOcDm5TWfj4jJ5TXHAwcBG5fbyDoPBO7NzKcCxwIfLetaEzgS2BHYATiy3vySJEmSJElS73WtKZWZlwD3jBp+CXBquX8qsG9t/PTMfCgzbwVuAXaIiHWA1TLz0sxM4CujXjOyrrOBPcssqr2BCzPznsy8F7iQf2yOSZIkSZIkqYe6OVNqLE/KzDsAytcnlvHpwG215eaXsenl/ujxhV6TmY8A9wFrLWZd/yAiDoqIuREx96677lqGsiRJkiRJkrQ0mm5KLUqMMZaLGZ/oaxYezDwxM2dl5qxp06aNa0MlSZIkSZK07JpuSv2hHJJH+XpnGZ8PrFtbbgZwexmfMcb4Qq+JiCnA6lSHCy5qXZIkSZIkSeoTTTelzgFGroZ3APDd2vicckW9DahOaH5FOcTv/oiYXc4Xtf+o14ys6xXAj8p5p84H9oqINcoJzvcqY5IkSZIkSeoTU7q14oj4BrA7sHZEzKe6It5HgDMj4kDgd8ArATLzuog4E7geeAQ4JDMfLas6mOpKfisC55UbwJeBr0bELVQzpOaUdd0TEccAV5bljs7M0SdclyRJkiRJUg9FNblIs2bNyrlz5475XIx1lqpxmOhba5555g1WXptrM88883qX1+bazDPPvN7ltbk288wzr3d5i8uKiHmZOWus5/rlROeSJEmSJEkaIjalJEmSJEmS1DibUpIkSZIkSWqcTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqUkSZIkSZLUOJtSkiRJkiRJapxNKUmSJEmSJDXOppQkSZIkSZIaZ1NKkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcTalJEmSJEmS1DibUpIkSZIkSWqcTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmN60lTKiJ+ExHXRMTVETG3jK0ZERdGxM3l6xq15Q+PiFsi4qaI2Ls2vl1Zzy0RcVxERBlfPiLOKOOXR8TMxouUJEmSJEnSIvVyptQembl1Zs4qjw8DfpiZGwM/LI+JiM2AOcDmwD7A5yNicnnN8cBBwMbltk8ZPxC4NzOfChwLfLSBeiRJkiRJkjRO/XT43kuAU8v9U4F9a+OnZ+ZDmXkrcAuwQ0SsA6yWmZdmZgJfGfWakXWdDew5MotKkiRJkiRJvderplQCF0TEvIg4qIw9KTPvAChfn1jGpwO31V47v4xNL/dHjy/0msx8BLgPWKsLdUiSJEmSJGkCpvQod5fMvD0inghcGBE3LmbZsWY45WLGF/eahVdcNcQOAlhvvfUWv8WSJEmSJEnqmJ7MlMrM28vXO4FvAzsAfyiH5FG+3lkWnw+sW3v5DOD2Mj5jjPGFXhMRU4DVgXvG2I4TM3NWZs6aNm1aZ4qTJEmSJEnSEjXelIqIlSNi1ZH7wF7AtcA5wAFlsQOA75b75wBzyhX1NqA6ofkV5RC/+yNidjlf1P6jXjOyrlcAPyrnnZIkSZIkSVIf6MXhe08Cvl3OOz4F+Hpm/iAirgTOjIgDgd8BrwTIzOsi4kzgeuAR4JDMfLSs62DgFGBF4LxyA/gy8NWIuIVqhtScJgqTJEmSJEnS+IQTiCqzZs3KuXPnjvncRK/bN9G31jzzzBusvDbXZp555vUur821mWeeeb3La3Nt5plnXu/yFpcVEfMyc9ZYz/Xq6nuSJEmSJEkaYjalJEmSJEmS1DibUpIkSZIkSWqcTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxtmUkiRJkiRJUuNsSkmSJEmSJKlxNqUkSZIkSZLUOJtSkiRJkiRJapxNKUmSJEmSJDXOppQkSZIkSZIaZ1NKkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcTalJEmSJEmS1DibUpIkSZIkSWqcTSlJkiRJkiQ1zqaUJEmSJEmSGmdTSpIkSZIkSY2zKSVJkiRJkqTG2ZSSJEmSJElS42xKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmNsyklSZIkSZKkxrW6KRUR+0TETRFxS0Qc1uvtkSRJkiRJUqW1TamImAx8DngesBmwX0Rs1tutkiRJkiRJErS4KQXsANySmf+XmX8HTgde0uNtkiRJkiRJEu1uSk0Hbqs9nl/GJEmSJEmS1GNTer0BXRRjjOVCC0QcBBxUHv4lIm6aQM7awB/H3ICxtmDZmWeeeX2W1+bazDPPvN7ltbk288wzr3d5ba7NPPPM613eErLWX9QTbW5KzQfWrT2eAdxeXyAzTwROXJaQiJibmbOWZR3mmWeeeb3MMs8884Ynr821mWeeeb3La3Nt5plnXnfz2nz43pXAxhGxQUQsB8wBzunxNkmSJEmSJIkWz5TKzEci4lDgfGAycFJmXtfjzZIkSZIkSRItbkoBZOb3ge93OWaZDv8zzzzzzOuDLPPMM2948tpcm3nmmde7vDbXZp555nUxLzJzyUtJkiRJkiRJHdTmc0pJkiRJkiSpT9mUkiRJkiRJUuNsSo1TRERDOZMiYurI/SYyS5b1dSan8fraXFvJsr7O5Fhfd3JbW1+baytZra6vF5p6T5vMGpb91+S+6wXrG2xtrq/NtYH1DbrG/q/1nFKLFhFvB27KzAvK48guvmERcTzwV2At4H2ZeVe3skqe9XU2r7H62lxbybO+zuZZX2fzWltfm2srea2ur2SumZn3dDunlrc7cAfwm8x8qMtZrd5/bd53Jc/6Opu3O9bXqazW1lbyrK+zebtjfR3V2r/4LKuIOAJ4I/DMiHgRQGZmt7qFEXEksBJwOHAFMKcbObU86+tsXmP1tbm2kmd9nc2zvs7mtba+NtdW8lpdX8n8EnB4RGzV7ayS9yngvcCrgCMjYnIXs1q9/9q870qe9XU2z/o6l9Xa2kqe9XU2z/q6wKbUov0EeBFwLjA7Il4M1Qcg6MpUthuo/gr393J/3fqT0fkp49bXWU3W1+bawPqsb+lYX+e0uTZoeX0R8Wxgesl6ZkRsWXuu442biHgVsG5mvhA4HngUeKzTOTWt3X9t33fW1/E86+tcVmtrK3nW19k86+sSD98bJSKmZOYjtcdrANsBewDnAb8G/pyZD3Qjr4ytDnw0M98aEU8Dbs/M+7uRZ32dzStjXamvzbWNlWd9nc0rY9bXobw21dfm2sbKa1t9o3KWB9amat4sD/wkM6/udE4tb53MvKPc/zbwkcy8PKJzh9UNy/5r474blWd9nc2zvs5ltba2kmF9nc2zvi5wplRNREzOzEei8m8RsUpm3gvMA84C3gacA2zQhbwjImKV8tSKwN8j4oXAcV3Ks77O5nW1vjbXNkae9XU2z/o6m9eq+tpc2xh5rauvZNY/y/09M38PXAQ8CGwZEa+KiP0jOvNX1VF5/692/0bgrojYAHhjJ/Lavv/avO/GyLO+zuZZX+eyWlXbGHnW19k86+uGzPRWu1E16t4HXAh8HViljM8CfgE8r8t5q5bxM4ErgT2tz/raXpv1WZ/19a6+NtfW9vqAqeVrAO8HPgZMKmNrAh8H/g94VpfyPlrLez1wAvBDYHf333DvO+uzvn6tr821WZ/19Xt9i7o5UwqIiDfVOn/vBlajmiJ3HXBKRKwE/Ak4KDPPW9Yu4RLyToqI5YBrgQ9m5g+XJWscedbX2byO1tfm2saRZ32dzbO+zuYNdH1trm0ceQNf34jMfLjkfoVqNs8qwNHlr573AusDB2fmT7uUtypwTMmbBBwE/EdmXrwsOcOw/9q67xaTZ32dzbO+zmW1prZF5FlfZ/Osrxu62fEahBvwCeCrtceHAMfWHp8GnASsVB5P6XLe16k6kis0lGd9fVpfm2uzPuuzvt7V1+bahqG+so4PA68t958FnFXuj/wV9f3l8RPK1+hi3ieA95bHW3egtlbvvzbvO+uzvn6ur821WZ/19Xt9S7oN9UypiDgIeDZwYHm8BfAFYHJEHB0R2wB3AXcA34uIlXPUyS+7kPcH4G7g3IbyrK8P62tzbdZnfdbXu/raXNsw1FcyPgNMz8yvlaErgQci4uWZeQ/VjJ4tIuLfgfthwRXqupR3DbBDRBxR7lP+6jqRrFbvvzbvu3HmWV9n86zP2sabZ32dzbO+Dhvaq+9FdfWWbYGtgCvK/R2opqgl1fGUqwDfyMwLI+JQ4DuZOd8887qZ1+bazDPPvN7ltbm2YcgrmdsB/wm8JzOviYh9qabzPwS8EfgzsA1wJPAc4JHM/OAg5LV9/7V535lnXj/ntbk288zr97xxywamY/XbDdgH+BywLvBcqhOI3QA8bdRyK5lnXpN5ba7NPPPM611em2sbhryyrleUr/8EHEr1ofIiFpyIe03gDcAO5fETgLUHIa/t+6/N+8488/o5r821mWdev+ctzW2YZ0q9h+qN/gIwnWr69jzgssz8s3nm9SqvzbWZZ555vctrc21Dknca8LvMPCIi9gNeB3wuM7/f6awe5bV2/w3BvjPPvL7Ma3Nt5pnX73lLY5ibUp8HNgEeoOoUPhXYGfg1cGZmPmaeeb3Ia3Nt5plnXu/y2lxbm/MiYkqW8xhFxM+AX2bmIRGxP9VfNS8FrhzUvFpu6/Zf2/edeeb1a16bazPPvH7Pm4hJvQrupYh4HbBiZj4HuAz4ANWHkCuAy7vwwcc88/ouyzzzzBuevDbX1va8zHwkIqZExDeBbwMbR8SnMvMrVCc8fS6w8qDmQXv3X9v3nXnm9Wtem2szz7x+z5voRrb+Bkwe9XgH4Ija428Dn6Zc7tc885rKa3Nt5plnXu/y2lzbMOSVdU6q3V8BOKH2+AzgjHJ/xqDltX3/tXnfmWdeP+e1uTbzzOv3vGW5DcVMqcx8NCImRcSnI+JtwKrArIg4LCJeCUwBzs7MB80zr8m8Ntdmnnnm9S6vzbUNQ15ETM7Mx6KyBfAYMC0iDimLfATYMiI2zGW4ml+v8tq8/9q+78wzr1/z2lybeeb1e96yavU5pSIishQYEWcClwCTge2odsRs4MnAVZn5A/PMayqvzbWZZ555vctrc23DkDcqexLwdeBPwO+AnwHvpzoZ917AOzPzZ4OUNyz7r437zjzzBiGvzbWZZ16/5y2T7PFUrW7dKFO1qT58rAq8CdiQaor2S8rYk8wzr+m8Ntdmnnnm9S6vzbUNQ17Jqk+13w94L7ASMB/YgOqcDxsAWwxaXtv3X5v3nXnm9XNem2szz7x+z+vUrbWH72Xmo+XukVSXO9wQ+AbVh5GrgYuAdc0zr+m8Ntdmnnnm9S6vzbUNQ15tqv2kiHhzGZ4EfJnq6nBTgTdl5q2Zee2g5bV5/7V935lnXr/mtbk288zr97xOal1TKiLeFhGzyv1PA+tn5vHAh4BfUnUKvwz8e2bONc+8pvLaXJt55pnXu7w21zYMeSVno6zOexTAfwJ/BS4E9qa6Ms4vgc8BtwxaXtv3X5v3nXnm9XNem2szz7x+z+u47IPpWp28AWuUr9tSTU27CHhBGVsFWBFYzzzzms5rc23mmWde7/LaXNuQ5G0MzCn3NwTuBJ5fHm8EHAX8B/C8Ac1r7f4bgn1nnnl9mdfm2swzr9/zunHr+QZ0rJDa5X6BpwDXUXUGnwpc0OmdYJ55/ZhlnnnmDU9em2sbhrwx8t8LPAPYFbgU2L6MTxrEvGHaf23bd+aZNyh5ba7NPPP6Pa+j297rDejQDhg5mWUAG5b761BNWdsLeBrVFVfWMs+8JvPaXJt55pnXu7w21zYMefXM2uMXAqcCTwe2B66ifKActLy277827zvzzOvnvDbXZp55/Z7XzVuUAgZeVJc8PBO4h+rs8t+gOn7yW8AxwE8y86/mmdd0XptrM88883qX1+ba2p4XEZOynIwUOJxqJs/tVNPuDwQ+TnV1uMjMSwctbySTFu6/tu8788zr17w212aeef2e120D35SKiMjMjIiPAb8HfgCcRHVllS8B/4/q3AGXm2dek3ltrs0888zrXV6baxuGvFruJOA7wI1UF565CzgbmAG8BXhzZj4waHnDsP/auu/MM6/f89pcm3nm9XteNw1sUyqqSx4+Wnu8GvB3qilr3wOeSHVir3/p0Ics88zruyzzzDNvePLaXNsw5JWMSZn5WLm/BbB3Zn4yIi6lmt0zl+oD5p2Zee8g5bV9/7V535lnXj/ntbk288zr97ymTOn1BkzEyAeR0h08AfgFcE9mnhER91NdAnF74GOd/OBjnnn9lGWeeeYNT16baxuGvJI5MtU+gLWBGzLz2oj4BHAi1V86DwbO6OQH1yby2r7/2rzvzDOvn/PaXJt55vV7XpMGeaZUAJ+juuThucD7gfOAW4F/ofogcr555jWd1+bazDPPvN7ltbm2tudFPH6IWVCd6+hhqkPNzgNWAl4EbAUclpk/HbS8kUxauP/avu/MM69f89pcm3nm9Xte47IPzrY+kRvVX8BuBZ5ZHq8PHFPur2Seeb3Ka3Nt5plnXu/y2lzbMOSV9X4Y+FS5vytwHNXV4V4E7DHIeW3ff23ed+aZ1895ba7NPPP6Pa+p20AevgeQmVdGxBHAOyLib8AjwC5RnVPgL+aZ16u8Ntdmnnnm9S6vzbUNQ15ETAMeAJ4eEetm5iURsSdwW2beMOh5bd5/bd935pnXr3ltrs088/o9r1G97oot7Y1yyGHt8SuB66kug7jTIOYBkxbz3MDnLWFbWpvXptpGfx+0La/t+6+tecC6bc6rrbt1/+8NU96orCcCBwFfBF4KXAbsOMh5w7L/2rjvzDNvEPLaXJt55vV7XlO3nm/AON74NwDPAZ5cG5s0apnnA18DNhvAvJXHymhR3uvL+7l62/KoTiT3IuCJDdXWdN77gX1rj8dqFA1y3rHADktYZpDzPgc8d9TY6F/mBjnvs8BnGnw/G8sDjgBeDWzdUG1N570W2B1YqzbWzf9nG82rrXPKIsbXoTrn0feBV421Pf2cx5B8Lmtq3/Ugb/WW563X8rwN25q3uO/fLtXWdN7zWp43reV5K496PKlNeb2+TaKPRcS/A3sCewP7R8QcgCxnnS8n+iIzv091wq93RMSKA5R3LvCLiFizZEwa9fyg532R6r3cj+qXgtHPD2xeydqN6oPw80fndKm2JvOeRHW+jvdHxAvLujNqBjzvJOChzLwiIlYd9Vwb8t4IPJyZF0bEOhGxfFl/tiTv08B0YKuIeMkYzw9sXlRXUNkUeDLV/0f15yZ1obam8z4GPBt4HnBARLyorH/k8sYDnVfWeUhETM3MRyJi8ujnM/MOqgbKWcDWETFzZHv6PS/a/7nsSOATEfGKsZ7vwr5rOu8Y4PiIePWo8ZH3cdDz3gx8MSKeMdbzLcj7JHDAGOPdej8by4uILwBvH2N8UqezepR3HHBulM+4o7Ug78PAFyLiDRGxQm28W/82m857P/CRiHhZRKxTMh5rS14/6OumFLAqcEpmvg/4ObBRRLwMql92Mqu2YETsB6wA/HNm/m0Q8iLiycA1wJepvqnGbBQNcN7mwIOZ+TrgPcDeEbH7IpYdqLyI2AS4PzPnZObBwF4R8ewy3o1/K43mlXX+gerSomcD74mIV9WzBjkvIqYDWwLXR8TRVB8o3xcRs+uZg5pXc2X5QPlfwGEjvxQMel5E7FJW+1LgbVTf60+qPR+DmhdVQ2FtqhNY/hewY0TsFRFbUW3EYx2urdG8YlXgk8BhVFeOeWVEPKdkdHTf9SIvIj4CfBw4LapG0aNRaxTVfuH5PXAbcBXV+SEGIg9YjWY/lzWWFxEfAjaj+r/okIjYZ9Tznd53TecdDWwMfBR4VdQaKZmZg55X/I3qKlTvi4hnjtqeyYOcF1WDfbPMPDIitoyIbUee69L+aywvInYCHsjMgyPiJeUz7hPL+h/rwnvZdN4Hy91dqL4XZo56vtP7rum8dwGbA4dSnXh7h5Hnyr+VTr+fTee9GdgW+DSwHfDmiNi0ltfp97PRvL6RfTBda1E34FnAe6mmqE0C9gDeB6xYW+ZlwMXApoOUBwRlqj3wbuA7lGmIo/JeOqB5q1E7zAx4I/Cscn/5Qc4Dlgcml/vPozqXxWeAo4CndOHfStN5y5Wvz6I6pGd14HyqX+y2HvS8sq7NgEupDsl6EtWx2e+jdsWmAc/bpGQdVx7vBvw78NRBzwPWqN1fB/jPkXWOfJ+U+5362dJYHjAZeE353j4W+AVwOHAMMLsLtTWaV9b1eqrDsdYu/25OBQ4etcxA5lF9bphd7v8rcCYLfr5NonY4K9V5kH7CYg4/6Kc8YNXydVeqn11d/ZzUg7xVgFfV3r89gY93cd81nbcysE/t8ceommEHUP2fO9B5tXWtCLyJ6v+gk6l+oVurPBeDmgdMpToK4MNU54/5PHAJ8E/lve70/ms674lU/++cSHUlsS8B7wK26Ma+azIPmAa8tvb4MMq576j+D+70e9loXlnPXsAh5f5/AKeV7/VnjFpuUPN2Bg4q99cCfkTVEJvMwofUDWRev9x6vgFj7IhDgKnl/gZUH5BfB8woYycD25b7LwV+CmwyoHn1Zsm7yzfVXlS/sEYL8qbWxg8FXk7VzX7nIOaNyppc1rkNC87TdSwLfkHodG2N5o28n8Dx5eungPuB/VqUt27t/hSqX+i2bEMe1aFYRwDfZEED5RNQXSq2BXn1/5QPBs6jduz9oOWNylqL6lDBL1DOB0H1ofKlXaqt6bwtqf7483mqE1PvT/WzbPqg5tX/jbDw/7NHA2eV+yuw4A9DrwbOBZ4+CHnAkVQ/m/cFngD8G1Wzr1ufk5rOG3kvV6X62Typ/Lv5ZhlfhwV/GOrkvms0r/Z4c+DbwAuAM4CX1J4b6DyqZt+XgOdSff6bB3ywDXnle2E/4APl8UyqJvsLBzWvlrVK+T7/IlXDaFWqox5e26msXuaNGjuwrHf098nA5lH9/Hob1f+zvwZ2ovp9+sCW5G0NnED1f9L2wPeojuzYcZDz+u3W8w0YtTM+AvyV6j+tkQ+U25QfFB+iOvb3J5ROKNVfJJblg0g/5C1Xe/4LwL2UkwW3JG/kQ/GLqI57/SGw96DljZU16vmtqTrZ2y5rVj/llR96X6aamfUMYLu25I3K3aK8n5u2KO/pVDMKPkDVpP0ZC5qYbcirv58fBF5TezwwebWsM0dlvJGqoT6H6i/Uu3e4tl7kjcwKWavcZlA1NL/JgibRQOWVdYw0EIKqwVdvWP4z1SzJHwPrUv3x66cs2wfXxvKoPgudQXUowcVUh4OsS9Xo68bnpKbz6u/lU1i4yXcUVUPlJ8D6Xdh3TedNZ8HPs1XK1zlUh9ctB2w04HnrlvtPo5p5fVh5L0+nOsRu4wHPexK1mbtlfH+qkx1P7sL72dW8UVlPoLqYwZepzgk7ler8qSdTNdif2uHams6bzsIz4z8CvL32uBvf603mrUE1g24zqkOoAV5I1UxZpQvvZ9N5K5V/H2+haoTtTHWuyheXZTYctLx+vPV8A2o7Y/Q09LMp/1lTfYjcnOpDycghWct6Rvte59U/MI/kfpbyF4g25ZWxzYCbgb0GLW8J7+VkYD3gW8DzG9h3TefNoppVV39vBz3vrFreFKof7udR/nrbgrzRP8s2pZqWvnNL8ur/XkZ+4XkptUMFByVvCf9WtqL6sPwluvdzs+m8s0d9b68EnMKCWVn/cMXNfs4bI/vrVIeCHAbsVsbXB+6i/LwuY6sPQh7/eIjZc1hwiNmT6PznpEbzFvNe7lrGbgKuo/xhq4v7rsm8I4BdyviOwJWUBnRL8v4VeAfwXapZMFOBFVqQ941a3jPL+E5Up3PYdVDzStbpJeutVLOX3kN1nryfAHt2obam8+rf67uX8WdTneJgxRbknV7LeztVo+89VLMGn9OivHey4PPFRsCvKIdFDmpev916vgEjO6F8HT0N/ezR4y3LG5lmv2L5OjLLJ1iGD8v9mEf1gXLDkbxuv5+dyhtH1nIla5GXBB7QvLNrj//hEKYW5I28n1OpzjUzclhIt/9tNpU3bD87O9JcaDJvPN8LTe67BvMeP7ysfB2ZQdHt772O5o2RewTVh+T1gMuB2VR/RNiUBb+Ud/JnWlfzajmjDzH7dhl/cpf+vTSSt4T3cheqz0UfZcEvdZ38nu913myq83C+hwWHWLcpb2PKLPL6ci3Km0010+efKU2UBvZfR/PGyJpZsrYoWWsAT+tibU3nPf5elvGVR/K6/W+l4byNqf7AtQ3lULMGvtebyFu/5O1M9Tl+dxb8IWjg8vr11vsNGP809E79Em5eb/KeMmh5ba5tnHmXlbx1W5o38n5Ob3le0++necueNfK90NTP6abzBnbfjc6sPX421TkuzqW6IMX6VOenrB8iuSx/IGksb9T7udhDzLqw/xrLW8J7+TIWXAxmWX8B78e8lwNrNvh+NpX3wtr3QtCZRnA/5Y28n69g1KF1g5C3hKwXlKwDOlFXH+aNvJcHLOv3eJ/mjbyfr2tp3sj7+RpqM84GKa/fb5PosawuXzwJ+BrVNNG3R8Ru5elvUR0X+vHM/J15A513+6Dltbm2ceZtVPJua2neyPv5+5bnNf1+mrfsWSPfC039nG46b2D3HVSX0x7JjIgPRcSOVLN7dqI6/Ope4CvAnZn5cG07cxDyRr2fhwH/HBG7lrH9qA5p/XBm/nYi6+9l3jjfy9OAv2TmXWX7JvQ+9nHeV4H7M/Oeieb0cd4fWfC9kJn5WMvyRt7PP2fmvcuS1XTeOLLupnov71iWnD7OG/k5fceyfI/3cd7dVP9W7mxp3sj3wh8z82+DljcIogP7ceLh1Q55LCKOAB6jOg72LKrjKK+kunzzmpn584iIZf1HZ555/ZhlnnnmDU9em2sbhrxablCd92EK8CjVeeJWK0/vAnwnMy/qRFaTeYt5P98N/A/VCXPPy8yLlzWrF3kls5X7zjzz+j2vzbWZZ16/5/W7njSlImJyZj5ae/xsqis3vYDq5NvXU13J4oys/dXPPPO6ndfm2swzz7ze5bW5tmHIKxmPN7Ui4jVUVxI7gOpD5OuoLq5xPtWsgk785baxvHG+n9sBP83MuzrQnG06r7X7zjzz+jmvzbWZZ16/5w2Sxg/fiwlOQzfPvG7ntbk288wzr3d5ba5tGPJK5uTaB8nJwFVU5zl6LXA/1SFmW1KdSL0TH1wbyxvn+znIh7S1dt+ZZ14/57W5NvPM6/e8QdOrmVKtnh5n3uDmtbk288wzr3d5ba5tGPJK5iTgeKqr4XyX6mqoqwMrUTXBHsoOnvuhyby277827zvzzOvnvDbXZp55/Z43SKY0FRSx0NTr/VhwRYDR09VO71A30jzz+i7LPPPMG568Ntc2DHklc33gd2V9H6f6a+apwB5UJzn+JbArsFp24ETxTea1ff+1ed+ZZ14/57W5NvPM6/e8QdXI4Xttnx5n3uDmtbk288wzr3d5ba5tGPJKzhHACcCOZWgucElmXgOcQ3XJ5huBz3bog2tjeW3ff23ed+aZ1895ba7NPPP6PW+QNXb4Xtunx5k3uHltrs0888zrXV6ba2t7XkR8HFgR+C2wOfAZYDqwG9XV4ZYDjgQOysxfD1peyWzl/mv7vjPPvH7Na3Nt5pnX73mDrqszpSJi/YiI8nBkutqxVFO2p1JNV5tKNV2tEx9EzDOv77LMM8+84clrc23DkFcyDwSenJmHZubHge8AbwJuoZrhszPwr8AnOvTBtbG8tu+/Nu8788zr57w212aeef2e1wqZ2ZUbcATViStnl8f7AS8u92cCXwVWBZY3z7wm89pcm3nmmde7vDbXNgx5tdyDgf3K/Unl68uAU4Dty+OVBi1vGPZfW/edeeb1e16bazPPvH7Pa8OtKzOlynS1pwA/At4aEdsBDwC7RcRsYL3y/BMz8yHzzGsqr821mWeeeb3La3Ntw5A3yjXAMyLiiZn5GEBmfqvkvysiVgc6eXWcrucN0f5r3b4zz7wByWtzbeaZ1+95A6/jV9+rTVd7fXm8L9V0tc9Q/UVsZ2BPOj89zjzz+ibLPPPMG568Ntc2DHlj+BWwC7BzRMzNzPkRsT2wTsm8b5Dyhmz/tWrfmWfeAOW1uTbzzOv3vMGXHZ56Rcunx5k3uHltrs0888zrXV6baxuGvEVsw/rAu4DTgX+humLOcwcxb9j2X5v2nXnmDVJem2szz7x+zxv0W8evvhcRzwSeD3w6M++sjX8OeALwNuDP2aFg88zrxyzzzDNvePLaXNsw5C1mO6YCmwBrlrxfDmLeMO6/tuw788wbtLw212aeef2eN8i6cU6pXwH3UU1XmwFQm672ucy8r8MfRMwzrx+zzDPPvOHJa3Ntw5A3psx8ODOvy8yfNvFBsot5Q7f/WrTvzDNvoPLaXJt55vV73iDr+EwpgIhYH3gpMJvqsoevAg7PzAs7HmaeeX2aZZ555g1PXptrG4a8tnP/SZKkftWVphS0f3qceYOb1+bazDPPvN7ltbm2YchrO/efJEnqR11rSkmSJEmSJEmL0o1zSkmSJEmSJEmLZVNKkiRJkiRJjbMpJUmSJEmSpMbZlJIkSZIkSVLjbEpJkiRJkiSpcTalJEmS+lhEXBwRs3q9HZIkSZ1mU0qSJEmSJEmNsyklSZK0FCJi5Yg4NyJ+GRHXRsSry/h/RMSVZezEiIgyfnFEHBsRl0TEDRGxfUR8KyJujogPlmVmRsSNEXFqRPwqIs6OiJXGyN4rIi6NiKsi4qyIWGWMZS6OiI9GxBUR8b8R8axaxk/La6+KiJ3L+O4R8ZOIOLMs/5GIeG15/TURsVFZblpEfLPUeGVE7NK9d1mSJA0Dm1KSJElLZx/g9szcKjO3AH5Qxj+bmduXsRWBF9Ze8/fM3BU4AfgucAiwBfCGiFirLPM04MTM3BL4M/C2emhErA28H3hOZm4LzAXevYhtnJKZOwDvAo4sY3cCzy2vfTVwXG35rYB3As8AXg9sUl7/JeDtZZn/Ao7NzO2Bl5fnJEmSJmxKrzdAkiRpwFwDfCIiPgr8d2b+tIzvERH/AqwErAlcB3yvPHdO7bXXZeYdABHxf8C6wJ+A2zLz52W504B3AJ+o5c4GNgN+XiZhLQdcuoht/Fb5Og+YWe5PBT4bEVsDjwKb1Ja/srZNvwYuqG3vHuX+c4DNSjbAahGxambev4htkCRJWiybUpIkSUshM/83IrYDng/8Z0RcAHwM+DwwKzNvi4ijgBVqL3uofH2sdn/k8cjnsRwdNepxABdm5n7j2MyRjEdr6/9n4A9Us6ImAQ+Osfzobaxv3yRgp8z82zjyJUmSlsjD9yRJkpZCRDwF+GtmnkY1k2lbFjSg/ljO8/SKCax6vYjYqdzfD/jZqOcvA3aJiKeW7VgpIjZh/FYH7sjMx6gO0Zu8lNt3AXDoyIMy40qSJGnCbEpJkiQtnWcAV0TE1cC/AR/MzD8BX6Q63O07wJUTWO8NwAER8Suqw/+Orz+ZmXcBbwC+UZa5DHj6Uqz/82X9l1EduvfAUm7fO4BZ5UTs1wNvXcrXS5IkLSQyR88MlyRJUpMiYibV+am26PW2SJIkNcWZUpIkSZIkSWqcM6UkSZIkSZLUOGdKSZIkSZIkqXE2pSRJkiRJktQ4m1KSJEmSJElqnE0pSZIkSZIkNc6mlCRJkiRJkhpnU0qSJEmSJEmN+/8aYp+qUT3fZgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1440x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define the size of the figure\n",
    "plt.figure(figsize=(20, 5))\n",
    "ax = plt.gca()\n",
    "\n",
    "df1.plot(x='sample', y='reads',kind=\"bar\", color='blue', ax=ax, linewidth=3)\n",
    "plt.title('Read number by sample', fontsize=10)\n",
    "plt.xlabel('sample name')\n",
    "plt.ylabel('reads number')\n",
    "plt.xticks(\n",
    "        rotation=45,\n",
    "        horizontalalignment='right',\n",
    "        fontweight='light',\n",
    "        fontsize='8',\n",
    "    )\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Plot the length mean "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3sAAAF+CAYAAAA7nwN5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA8u0lEQVR4nO3dd7xkdX3/8ddnd0FYirRFkbaAgkpxxRUREBEFiYKBWFGjxoINNVFAwQLYC4qJBQVFVJTYiQYVS36KGkhYEFiwJUFEkMgiREQUET6/P873usNwy5RzZmbPfT0fj3ncmTPnvud9zq3fOS0yE0mSJElSuywYdwFJkiRJUv0c7EmSJElSCznYkyRJkqQWcrAnSZIkSS3kYE+SJEmSWsjBniRJkiS1kIM9SdKsIuKOiLgkIi6NiIsjYq8Bcz4SEQ8cosdGEfGSuubr8TW/ExHLa8iprVMPr/WciHh/r9Nrer37dDy+KiI2q/t1JEn9c7AnSZrLHzJzWWY+CDgWeNsgIZn5/Mz80RA9NgJ6GTD1Ot8obUQPnaKypv1tfg5wn7lmkiSN3pr2B0WSNF4bAjcBRMT6EfHtsrVvZUT8dZm+XkScU7YEXh4RTy3T/7KVLCIOKp93aUR8u/tFImLniPjPskXxsoi4H/B2YIcy7V0zvf40820REeeVx5dHxCOmeb03RMSF5flTIyI6nn5mRPx7eW6PMv8mEXF26XZBROxWpp8QEUd15F4eEUu7O3W99tKI+HFEfBC4GNg6Io4ufS6LiBM75j07Ii6KiCsi4oiO6X8XET+LiO8Ce8/1RYyIJRHxhfIaF0bE3h39Ty9fqysj4uUdn/P6iPhJRHwzIs6KiKMi4knAcuBTZdnWLbO/rOPrcv+5+kiSmrFo3AUkSRNv3Yi4BFgH2ALYv0z/I3BYZt5cdtu7ICK+DBwE/CozHw8QEffsDIuIJcBpwL6Z+fOI2GSa13wR8I+Z+amIWBtYCLwG2CUzl5WcRTO8fvd8rwLOzcy3RMRCYPE0r/f+zHxjmf+TwMHAV8pz62XmXhGxL3A6sAtwIvDDzDw0IvYHPgEsm2Ud3qXTNHYC/i4zXxIRBwL3A/YAAvhyROybmecBz83MG8ug6sKI+AKwdunzEOC3wP8DfjhLF4B/BE7OzO9HxDbAucADynP3Bx4FbAD8NCJOAR4EPBF4MNX/DhcDF2Xm5yPiSOCozFxR1h/ADZm5e1S7rh4FPH+OPpKkBjjYkyTN5Q8dA6eHA5+IiF2oBiJvLYOgO4EtgXsBK4GTIuIdwL9m5ve68vYEzsvMnwNk5o3TvOb5wGsjYivgi5n5X3fd2AazvH63C4HTI2It4OzMvGSaeR4VEcdQDQQ3Aa5g9WDvrNLzvIjYMCI2AvahGvyQmf8WEZt2D2r79IvMvKDcP7DcpgZs61MN/s4DXh4Rh5XpW5fp9wa+k5mrACLiM8COc7zeY4AHdqzTDSNig3L/nMy8DbgtIq6nWqf7AP+SmX8or/GV7sAuXywfLwL+Zo55JUkNcTdOSVLPMvN8YDNgCfCM8vEhZTD4a2CdzPwZ1VamlcDbIuINXTEB5Byv82ngCcAfgHPL1rNu077+NFnnAfsC1wKfjIhn3aVMxDrAB4EnZeauVFsdO3O6u2ZZhru9FPBn7vq39W59ZvD7zkrA28pxkssy876Z+dGI2I9qkPbwcvzkDzvyZ12f01hQcqZeY8vM/F157raO+e6gemN4uuWdzVTG1OdLksbAwZ4kqWfl+KuFwG+AewLXZ+btEfEoYNsyz32AWzPzTOAkYPeumPOBR0bEdmX+u+3GGRHbA1dm5j8BXwZ2A35HtWvhlGlfv3u+iNi2zHca8NFp+kwNmG6IiPWBJ3U9P3XM4T7AbzPzt1Rb2Z5Rpu9HtdvizcBVU/kRsTuw3XSd5nAu8NzShYjYMiI2L8t7U2beWr4Oe5b5/wPYr2xdXAt4cg+v8Q3gyKkHEbFsjvm/DxwSEeuUXo/veK6fZZMkjZDvtkmS5jJ1zB5UW3ienZl3RMSngK9ExArgEuAnZZ5dgXdFxJ3A7cCLO8Myc1U5ucgXozrz5PXAAV2v+VSqE6PcDvwv8MZyrNoPIuJy4GvAO6Z7/cz8Tdd8lwNHl6xbgLts2cvM/4uI06i2RF5Ftdtnp5si4t+pTk7z3DLtBOBjEXEZcCvw7DL9C8Czyvq6EPjZdJ0y8+hp1vNUn29ExAOA88tulrcAzwS+DryovOZPgQvK/NdFxAlUg+jrqI6nWzhTfvFy4AMlaxHV4PVFs3S6sBwPeSnwC2AF1fGBAGcAH4qIPwAPn+N1JUkjFJn97vkhSZLmm4hYPzNviYjFVIPDIzLz4nH3kiTNzC17kiSpF6dGxAOpdnv9uAM9SZp8btmTJEmSpBbyBC2SJEmS1EIO9iRJkiSphRzsSZIkSVILNXaClog4HTiY6tpGu5RpJwAvAFaV2Y7LzK+W6wJ9hOraRIuAT2Tm2+Z6jc022yyXLl3aQHtJkiRJmnwXXXTRDZm5ZLrnmjwb5xnA+4FPdE0/OTNP6pr2ZOAemblrOaXzjyLirMy8arYXWLp0KStWrKirryRJkiStUSLiFzM919hunJl5HnBjr7MD60XEImBd4E/AzU11kyRJkqS2G8cxe0dGxGURcXpEbFymfR74PXAdcDVwUmb2OlCUJEmSJHUZ9WDvFGAHYBnVwO7dZfoewB3AfYDtgFdFxPbTBUTEERGxIiJWrFq1arpZJEmSJGneG+lgLzN/nZl3ZOadwGlUgzyApwNfz8zbM/N64AfA8hkyTs3M5Zm5fMmSaY9DlCRJkqR5b6SDvYjYouPhYcDl5f7VwP5RWQ/YE/jJKLtJkiRJUps0eemFs4D9gM0i4hrgeGC/iFhGdUKWq4AXltk/AHyMavAXwMcy87KmukmSJElS2zU22MvMw6eZ/NEZ5r2F6vILkiRJkqQajONsnJIkSZKkhjnYkyRJkqQWcrAnSZIkSS3kYE+SJEmSWsjBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaiEHe5IkSZLUQg72JEmSJKmFHOxJkiRJUgs52JMkSZKkFnKwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUaG+xFxOkRcX1EXN4x7YSIuDYiLim3x3U8t1tEnB8RV0TEyohYp6lukiRJktR2TW7ZOwM4aJrpJ2fmsnL7KkBELALOBF6UmTsD+wG3N9hNkiRJklqtscFeZp4H3Njj7AcCl2XmpeVzf5OZdzTVTZIkSZLabhzH7B0ZEZeV3Tw3LtN2BDIizo2IiyPimDH0kiRJkqTWGPVg7xRgB2AZcB3w7jJ9EbAP8Izy8bCIePR0ARFxRESsiIgVq1atar6xJEmSJK2BRjrYy8xfZ+YdmXkncBqwR3nqGuC7mXlDZt4KfBXYfYaMUzNzeWYuX7JkyWiKS5IkSdIaZqSDvYjYouPhYcDUmTrPBXaLiMXlZC2PBH40ym6SJEmS1CaLmgqOiLOozqq5WURcAxwP7BcRy4AErgJeCJCZN0XEe4ALy3NfzcxzmuomSZIkSW3X2GAvMw+fZvJHZ5n/TKrLL0iSJEmShjSOs3FKkiRJkhrmYE+SJEmSWsjBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaiEHe5IkSZLUQg72JEmSJKmFHOxJkiRJUgs52JMkSZKkFnKwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWamywFxGnR8T1EXF5x7QTIuLaiLik3B7X9TnbRMQtEXFUU70kSZIkaT5ocsveGcBB00w/OTOXldtXu58DvtZgJ0mSJEmaFxY1FZyZ50XE0l7nj4hDgSuB3zfVSZIkSZLmi3Ecs3dkRFxWdvPcGCAi1gNeDZw4hj6SJEmS1DqjHuydAuwALAOuA95dpp9ItXvnLXMFRMQREbEiIlasWrWqsaKSJEmStCZrbDfO6WTmr6fuR8RpwL+Whw8DnhQR7wQ2Au6MiD9m5vunyTgVOBVg+fLl2XhpSZIkSVoDjXSwFxFbZOZ15eFhwOUAmfmIjnlOAG6ZbqAnSZIkSepNY4O9iDgL2A/YLCKuAY4H9ouIZUACVwEvbOr1JUmSJGk+a/JsnIdPM/mjPXzeCfW3kSRJkqT5ZRxn45QkSZIkNczBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaiEHe5IkSZLUQg72JEmSJKmFHOxJkiRJUgs52JMkSZKkFnKwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWcrAnSZIkSS3U2GAvIk6PiOsj4vKOaSdExLURcUm5Pa5MPyAiLoqIleXj/k31kiRJkqT5oMkte2cAB00z/eTMXFZuXy3TbgAOycxdgWcDn2ywlyRJkiS13qKmgjPzvIhY2uO8P+x4eAWwTkTcIzNva6ScJEmSJLXcOI7ZOzIiLiu7eW48zfNPBH7oQE+SJEmSBjfqwd4pwA7AMuA64N2dT0bEzsA7gBfOFBARR0TEiohYsWrVqgarSpIkSdKaa6SDvcz8dWbekZl3AqcBe0w9FxFbAV8CnpWZ/zNLxqmZuTwzly9ZsqT50pIkSZK0BhrpYC8ituh4eBhweZm+EXAOcGxm/mCUnSRJkiSpjeY8QUtE7AgcDWzbOX9mznp5hIg4C9gP2CwirgGOB/aLiGVAAlexenfNI4H7Aq+PiNeXaQdm5vV9LIskSZIkqYjMnH2GiEuBDwEXAXdMTc/Mi5qtNrfly5fnihUrxl1DkiRJksYiIi7KzOXTPdfLpRf+nJmn1NxJkiRJktSgGQd7EbFJufuViHgJ1clT/nI5hMy8seFukiRJkqQBzbZl7yKqY+uiPD6647kEtm+qlCRJkiRpODMO9jJzO4CIWCcz/9j5XESs03QxSZIkSdLgern0wr/3OE2SJEmSNCFmO2bv3sCWwLoR8WBW7865IbB4BN0kSZIkSQOa7Zi9xwLPAbYC3tMx/XfAcQ12kiRJkiQNabZj9j4OfDwinpiZXxhhJ0mSJEnSkHq5zt62EfHKrmm/BS7KzEvqryRJkiRJGlYvJ2hZDryI6vi9LYEjgP2A0yLimOaqSZIkSZIG1cuWvU2B3TPzFoCIOB74PLAv1bX43tlcPUmSJEnSIHrZsrcN8KeOx7cD22bmH4DbGmklSZIkSRpKL1v2Pg1cEBH/Uh4fApwVEesBP2qsmSRJkiRpYHMO9jLzTRHxNWBvqmvtvSgzV5Snn9FkOUmSJEnSYHrZsgfwQ+BXU/NHxDaZeXVjrSRJkiRJQ5lzsBcRLwOOB34N3EG1dS+B3ZqtJkmSJEkaVC9b9l4B7JSZv2m6jCRJkiSpHr2cjfOXVBdRlyRJkiStIXrZsncl8J2IOIeOSy1k5nsaayVJkiRJGkovg72ry23tcpMkSZIkTbheLr1wIkBErJeZv2++kiRJkiRpWHMesxcRD4+IHwE/Lo8fFBEfbLyZJEmSJGlgvZyg5b3AY4HfAGTmpcC+DXaSJEmSJA2pl8EemfnLrkl3NNBFkiRJklSTni69EBF7ARkRa0fEUZRdOmcTEadHxPURcXnHtBMi4tqIuKTcHtfx3LER8d8R8dOIeOxASzMBVl55PZ/+1qWsvPL6icuz2/iz6s6z2/izJj3PbuPPqjvPbuPPqjvPbuPPqjvPbuPPqjuv7m6jEpk5+wwRmwH/CDwGCOAbwMsz88Y5Pm9f4BbgE5m5S5l2AnBLZp7UNe8DgbOAPYD7AN8CdszMWbcgLl++PFesWDFr/1H5zW//wOPf8TpWbXL2X6YtufFQznn1m9n0nuuONc9udrPb5HerO89udrOb3exmN7tNXrcmRMRFmbl82ufmGuzNEHhSZh7Vw3xLgX/tYbB3LEBmvq08Phc4ITPPny1/kgZ7ex73qrt8E0xZcuOhXPDWd481z252s9vkd6s7z252s5vd7GY3u01etybMNtjr6Zi9aTxliD5HRsRlZTfPjcu0LYHO4wKvKdPWCCuvvH7abwKAVZuc3ffm3jrz7GY3u01+t7rz7GY3u9nNbnaz2+R1G4dBB3sx4OedAuwALAOuA6aGw9PlTbvJMSKOiIgVEbFi1apVA9ao18orrxvq+Sbz7DZYnt0Gy7PbZOTZbbA8uw2WZ7fB8uw2WJ7dBsuz22B5dXcbhxkHexGxyQy3TRlwsJeZv87MOzLzTuA0qmP0oNqSt3XHrFsBv5oh49TMXJ6Zy5csWTJIjdrtuv0WQz3fZJ7dBsuz22B5dpuMPLsNlme3wfLsNlie3QbLs9tgeXYbLK/ubuMw25a9i4AV5WPnbQXwp0FeLCI618hhwNSZOr8MPC0i7hER2wH3A/5zkNcYh12335wlNx467XNLbjyUXbfffGx5drOb3Sa/W915drOb3exmN7vZbfK6jcOMg73M3C4zty8fu2/bzxUcEWcB5wM7RcQ1EfE84J0RsTIiLgMeBfxDea0rgM8CPwK+Drx0rjNxTppzXv3mu30zTJ2pZ9x5drOb3Sa/W915drOb3exmN7vZbfK6jdpAZ+OcFJN0Ns4pK6+8npVXXseu229Ry2i/zjy7jT/LbpORN8nd6s6z2/iz7DYZeXYbf5bdJiPPbuPPqlvtl16YFJM42JMkSZKkUWni0guSJEmSpAm2aKYnImKT2T4xM2+sv44kSZIkqQ4zDvaozryZVJdZ2Aa4qdzfCLga2K7pcpIkSZKkwcx5Nk7gXOCQzNwsMzcFDga+OKqCkiRJkqT+9XLM3kMz86tTDzLza8Ajm6skSZIkSRrWbLtxTrkhIl4HnEm1W+czgd802kqSJEmSNJRetuwdDiwBvgScDWxepkmSJEmSJtScW/bKWTdfMYIukiRJkqSazDnYi4glwDHAzsA6U9Mzc/8Ge0mSJEmShtDLbpyfAn5CdamFE4GrgAsb7CRJkiRJGlIvg71NM/OjwO2Z+d3MfC6wZ8O9JEmSJElD6OVsnLeXj9dFxOOBXwFbNVdJkiRJkjSsXgZ7b46IewKvAt4HbAj8Q6OtJEmSJElD6eVsnP9a7v4WeFSzdSRJkiRJdZjzmL2I2DEivh0Rl5fHu5WLrEuSJEmSJlQvJ2g5DTiWcuxeZl4GPK3JUpIkSZKk4fQy2Fucmf/ZNe3PTZSRJEmSJNWjl8HeDRGxA5AAEfEk4LpGW0mSJEmShtLL2ThfCpwK3D8irgV+Djyj0VaSJEmSpKH0cjbOK4HHRMR6VFsC/wA8FfhFw90kSZIkSQOacTfOiNgwIo6NiPdHxAHArcCzgf8GnjKqgpIkSZKk/s22Ze+TwE3A+cALgGOAtYFDM/OS5qtJkiRJkgY122Bv+8zcFSAiPgLcAGyTmb8bSTNJkiRJ0sBmOxvn7VN3MvMO4Of9DPQi4vSIuH7qYuxdzx0VERkRm5XHa0XExyNiZUT8OCKO7WchJEmSJEl3Ndtg70ERcXO5/Q7Ybep+RNzcQ/YZwEHdEyNia+AA4OqOyU8G7lG2JD4EeGFELO11ISRJkiRJdzXjbpyZuXCY4Mw8b4YB28lUx//9S+fswHoRsQhYF/gT0MuAUpIkSZI0jV4uql6biHgCcG1mXtr11OeB31NdrP1q4KTMvHGU3SRJkiSpTXq5qHotImIx8FrgwGme3gO4A7gPsDHwvYj4VrnGX3fOEcARANtss01zhSVJkiRpDTbKLXs7ANsBl0bEVcBWwMURcW/g6cDXM/P2zLwe+AGwfLqQzDw1M5dn5vIlS5aMqLokSZIkrVlGNtjLzJWZuXlmLs3MpcA1wO6Z+b9Uu27uH5X1gD2Bn4yqmyRJkiS1TWODvYg4i+qC7DtFxDUR8bxZZv8AsD5wOXAh8LHMvKypbpIkSZLUdo0ds5eZh8/x/NKO+7dQXX5BkiRJklSDkZ6NU5IkSZI0Gg72JEmSJKmFHOxJkiRJUgs52JMkSZKkFnKwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWcrAnSZIkSS3kYE+SJEmSWsjBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaqHGBnsRcXpEXB8Rl0/z3FERkRGxWce03SLi/Ii4IiJWRsQ6TXWTJEmSpLZrcsveGcBB3RMjYmvgAODqjmmLgDOBF2XmzsB+wO0NdpMkSZKkVmtssJeZ5wE3TvPUycAxQHZMOxC4LDMvLZ/7m8y8o6lukiRJktR2Iz1mLyKeAFw7NajrsCOQEXFuRFwcEceMspckSZIktc2iUb1QRCwGXku1FW+6HvsADwVuBb4dERdl5renyTkCOAJgm222aa6wJEmSJK3BRrllbwdgO+DSiLgK2Aq4OCLuDVwDfDczb8jMW4GvArtPF5KZp2bm8sxcvmTJkhFVlyRJkqQ1y8gGe5m5MjM3z8ylmbmUaoC3e2b+L3AusFtELC4na3kk8KNRdZMkSZKktmny0gtnAecDO0XENRHxvJnmzcybgPcAFwKXABdn5jlNdZMkSZKktmvsmL3MPHyO55d2PT6T6vILkiRJkqQhjfRsnJIkSZKk0XCwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWcrAnSZIkSS3kYE+SJEmSWsjBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaiEHe5IkSZLUQg72JEmSJKmFHOxJkiRJUgs1NtiLiNMj4vqIuHya546KiIyIzbqmbxMRt0TEUU31kiRJkqT5oMkte2cAB3VPjIitgQOAq6f5nJOBrzXYSZIkSZLmhcYGe5l5HnDjNE+dDBwDZOfEiDgUuBK4oqlOkiRJkjRfjPSYvYh4AnBtZl7aNX094NXAiT1kHBERKyJixapVqxpqKkmSJElrtpEN9iJiMfBa4A3TPH0icHJm3jJXTmaempnLM3P5kiVL6q4pSZIkSa2waISvtQOwHXBpRABsBVwcEXsADwOeFBHvBDYC7oyIP2bm+0fYT5IkSZJaY2SDvcxcCWw+9TgirgKWZ+YNwCM6pp8A3OJAT5IkSZIG1+SlF84Czgd2iohrIuJ5Tb2WJEmSJOmuGtuyl5mHz/H80hmmn9BEH0mSJEmaT0Z6Nk5JkiRJ0mg42JMkSZKkFnKwJ0mSJEkt5GBPkiRJklrIwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWcrAnSZIkSS3kYE+SJEmSWsjBniRJkiS1kIM9SZIkSWohB3uSJEmS1EIO9iRJkiSphRzsSZIkSVILOdiTJEmSpBZysCdJkiRJLeRgT5IkSZJayMGeJEmSJLWQgz1JkiRJaiEHe5IkSZLUQg72JEmSJKmFGhvsRcTpEXF9RFw+zXNHRURGxGbl8QERcVFErCwf92+qlyRJkiTNB01u2TsDOKh7YkRsDRwAXN0x+QbgkMzcFXg28MkGe0mSJElS6zU22MvM84Abp3nqZOAYIDvm/WFm/qo8vAJYJyLu0VQ3SZIkSWq7kR6zFxFPAK7NzEtnme2JwA8z87YZMo6IiBURsWLVqlWN9JQkSZKkNd2iUb1QRCwGXgscOMs8OwPvmG2ezDwVOBVg+fLlOdN8kiRJkjSfjXLL3g7AdsClEXEVsBVwcUTcGyAitgK+BDwrM/9nhL0kSZIkqXVGtmUvM1cCm089LgO+5Zl5Q0RsBJwDHJuZPxhVJ0mSJElqqyYvvXAWcD6wU0RcExHPm2X2I4H7Aq+PiEvKbfNZ5pckSZIkzSIy19zD3iJiFfCLcfeYxmZUl5OYxDy7jT+r7jy7jT9r0vPsNv6suvPsNv6suvPsNv6suvPsNv6suvPq7laXbTNzyXRPrNGDvUkVESsyc/kk5tlt/Fl159lt/FmTnme38WfVnWe38WfVnWe38WfVnWe38WfVnVd3t1EY6aUXJEmSJEmj4WBPkiRJklrIwV4zTp3gPLuNP6vuPLuNP2vS8+w2/qy68+w2/qy68+w2/qy68+w2/qy68+ru1jiP2ZMkSZKkFnLLniRJkiS1kIM9SZIkSWohB3tDiIioIWNBRKw1dX/4Vn/JHbpbyam9n90GzrTbYJl2GyzTbhOgrmWtM2u+rbe62W0wdhuM3QY36f165TF7fYqIlwE/zcxvlMeRA67EiDgFuBXYFDg6M1dNSre6+9nNbnazW1u7lbxNMvPGYTI6svYDrgOuyszbhsxyvQ2eZ7fB8uw2WJ7dBsurrVvJ248a+02CiXxHblJFxHHAc4F9IuIQgMzMQUb+EXE8sBg4FvhP4GmT0q3ufnazm93s1tZuJe8jwLER8aBhckrWe4CjgKcAx0fEwiGyXG92s5vd7DamfpPCwV5/vgscApwD7BkRT4Dqjyf0vbn3x1Tviv6p3N+688nof9eYOrvV3c9udrOb3VrZLSL2B7YsOftExG4dz/XVKyKeAmydmQcDpwB3AHf2k9HF9TYAu9nNbvOrWxP9JkpmepvjBizqerwx8BjgLcA+wBbAeoNklWn3BD5U7u8EbDCObnX3s5vd7Ga3tnbryrgH1T8dLwJeASzrN6Mja4uO+18CHlbux6DL6XobKM9udrPbPOrWRL9Jubllbw4RsTAz/xyV10bE+pl5E3AR8DngJcCXge36zDouItYvT60L/CkiDgb+qZesurvV3c9udrOb3drareR1/v38U2ZeC3wL+COwW0Q8JSKeFTH3O8xdWf/bcf8nwKqI2A54bo9ZrrcB1pvd7Ga3+detiX4TadyjzTXhRrW769HAN4FPA+uX6cuBHwJ/NUTWBmX6Z4ELgUePq1vd/exmN7vZrY3dgLXKxwBeB7wTWFCmbQK8C7gSeMQAWe/oyPpb4EPAt4H9XG/NrTe72c1u86tbE/0m9eaWvRlExPM7RvGvBDakOvbhCuCMiFgM/B9wRGZ+bbYR/xxZp0fE2sDlwJsz89uj7FZ3P7vZzW52a2u3KZl5e8n8BNUWrvWBN5Z3iG8CtgVenJnfGyBrA+BNJWsBcATwhsz8zhDL6Xqzm93sZreG+02scY82J/EGnAR8suPxS4GTOx6fCZwOLC6P73Y8Qx9Zn6Z652CdubLq7lZ3P7vZzW52a2u38vxbgWeU+48APlfuT72j/LryeKPyccbjPObIOgk4qjxeNlsn19vg681udrPb/OvWRL9Jv7llr0tEHAHsDzyvPN4F+DCwMCLeGBEPBlZRXYPjKxGxXmb+eYisXwO/Ac6ZLavubnX3s5vd7Ga3tnYrn/8+YMvM/FSZdCHw+4h4YlbXeLoc2CUiXg/8Dlaf9XKArJXAHlFdOmFl+Zxpt8S53gZbb3azm93mX7cm+q0JvKh6h4jYGNgdeBDVNYR2B/ag2oybVPvyrg+clZnfjIgjgbMz85omsyY9z252s5vd2tqt5D0EeBvwqsxcGRGHUu0ydBvVdexuBh4MHE911ss/Z+abR5Dlehsgy252s9v869ZE3hojJ2Dz4iTcgIOAD1BdN+gAqoM+fwzs1DXf4lFmTXqe3exmN7u1tVuZ70nl498BR1L9o/AtVp+kZBPgOcAe5fFGwGYjyHK9DZBlN7vZbf51ayJvTbq5Za9DRLyK6ov7YarrduxPdbrqCzLz5nFlTXqe3exmN7u1uNuZwNWZeVxEHA48E/hAZn51gF61ZZU819sA7GY3u82vbk3krUkc7HWIiA8COwK/pxr13xfYC/gf4LOZeec4siY9z252s5vd2tYtIhZlORYtIr4PXJqZL42IZ1G9A3w+cOGos+pezrrzJnm92c1udptf3ZrIWxMtGHeBSRERzwTWzczHABcAJ1L9gftP4D/6/MeltqxJz7Ob3exmtzZ2y+pC4osi4gvAl4D7RcR7MvMTVCcBOABYb9RZdS9n3XmTvN7sZje7za9uTeStkXIC9iUdxw1Y2PV4D+C4jsdfAt5LObX0qLImPc9udrOb3drarcy/oOP+OsCHOh5/BvhMub/ViLNcbwNk2c1udpt/3ZrIW5Nv83bLXmbeERELIuK9EfESqgspLo+I10TEk4FFwOcz84+jzJr0PLvZzW52a2u3iFiYmXdGZRfgTmBJRLy0zPJ2YLeI2D5nOANlE1l1L2fdeZO83uxmN7vNr25N5K3p5t0xexERWRY6Ij4LnAcsBB5C9cXfE7g3cHFmfn1UWZOeZze72c1ube3WlbuA6qLh/wdcDXwfeB3ViUoOBF6Rmd8fVZbrbbgsu9nNbvOvWxN5a7ScgM2Lo7pRdluh+sO2AfB8YHuq3VX+uky716izJj3Pbnazm93a2q3kdO7uczhwFLAYuAbYjup4ju2AXUac5XobIMtudrPb/OvWRF5bbvNqN87MvKPcPZ7qlKvbA2dR/aG7hOp6G1uPOmvS8+xmN7vZra3dOnb3WRARLyiTFwAfpToj5VrA8zPz55l5+aiy6l7OuvMmeb3ZzW52m1/dmshrlXGPNkdxA14CLC/33wt8vNxfDzgVeBHVH7kDR5k16Xl2s5vd7NbWbuXzdigfA3gH8AxgM+DbJW874JvA40ec5XobIMtudrPb/OvWRF7bbmMvMJKFhI3Lx93LF/xbU19wYH1gXWCbUWdNep7d7GY3u7W42/2Ap5X72wPXA48rj3cATgDeAPzVKLNcb0OtN7vZzW7zqFsTeW28jb1AowvXcWpp4D7AFcBjqS4U+41+vvB1Zk16nt3sZje7tbXbNNlHAbsC+1JdXPehZfqCUWe53obPspvd7Db/ujWR16bb2As0tmCrD0QPYPtyfwuqzbgHAjtRnYVs01FmTXqe3exmN7u1tVtnXsfjg4GPA/cHHgpcTPknYRxZrrf+suxmN7vNv25N5LX51upLL5TTrn4WuJHqTDxnAb8Dvgi8CfhuZt466qxJz7Ob3exmtzZ2i4gFWQ7gB46l2rr1K6pdf54HvIvqjJSRmeePKqvu5aw7b5LXm93sZrf51a2JvLZr5WAvorqOUES8E7gW+DpwOtXZxj4C/C/V8Qn/McqsSc+zm93sZre2duvIXACcDfyE6kxtq4DPA1sBLwRekJm/H2WW623wLLvZzW7zr1sTeW3WqsFeVKddvaPj8YbAn6g2634F2JzqYM1j5voGqDNr0vPsZje72a2t3crnL8jMO8v9XYDHZua7I+J8qi1eK6j+abg+M28aYZbrbYAsu9nNbvOvWxN588WicReoy9QfuTLS/xDwQ+DGzPxMRPwOuJVqH9539vqPSx1Zk55nN7vZzW5t7Vbypnb3CapTcf84My+PiJOoTsn9E+DFwGd6/UeopizX2wBZdrOb3eZftyby5pO2bdkL4ANUp109B3gd8DXg58AxVH/kzh111qTn2c1udrNbG7tF/GV3xqA6Xu12qt0avwYsBg4BHgS8JjO/N6qsupez7rxJXm92s5vd5le3JvLmnZyAs8TUdaN61/LnwD7l8bbAm8r9xePKmvQ8u9nNbnZra7fyOW8F3lPu7wv8E9UZKQ8BHjXGLNfbAFl2s5vd5l+3JvLmy601u3ECZOaFEXEc8PKI+APwZ2DvqI5buGVcWZOeZze72c1ube0WEUuA3wP3j4itM/O8iHg08MvM/PG4ssD1NkiW3exmt/nXrYm8eWXco806bpTdUTsePxn4EdWpWB8+6ixmuYDjJOSNIqut3bq/P+r+fmvrepvEbsDWdXarO2+m75Fh19uk/b5sMq8jZ3PgCOA04DDgAuBh48pyvQ2XZTe72W3+dWsib77cxl5gwC/2c4DHAPfumLaga57HAZ8CHjiqrDLvetNlTEIe8LdlWe9ZU7fa8qgOqj0E2LymbrXlUR33cmjH4+kGaP10qy0POBnYY455+ulWWx7VsUMHdE3r/qe0n2615QHvB95X43qrLQ84DngqsKymbrXlAc8A9qPjQtzdv5v67FZrXpl/0QzTt6A6bu2rwFOme62msqj/70zdeevN8ly/y1pbVnn+njV2qy2rPL9Njd1qyyrPb19jt9qyyvMzfk8O0K22rPL8X9XYrbas8vySGrvVllWeX6/r8YJJypvvt7EX6LswvB74JNUFE48BntbxXNDxzx9wMNVZydZtOqvMcw7wM2CT8rj7j+/Y8qjeBTkT+Cjwkmme77dbbXkl65+BU4DndOcM2K2WPOBewL9Qnc734O6cAbrVlkd13ay3l/sb1LDeassDngu8t9zfArjHkN1qywPeC3wJOA/46xp+FmrLA06i+p30CuBVXc8tGKBbbXnAO6l+3t8BvBI4ZMj1VnfeS4G1yv2FM8yzJfB3VMd9LB1RVt1/Z+rOO57qd+WTZpmn12WtLavM+ybg08BTZ/ne6LVbbVll3hcA5wK71rDeassq874bOHGa6YOst9qyyrwfBk6ZZvqCAbrVllXm/SfgTjr+Ng/xNa0tq8z7VuCLVG/0rDPk17S2rDLv64D3AX8DbDFped6SBax5NgDOyMyjgR8AO0TE3wBkARARhwPrAP+QmX9oOisi7g2spPrn5cMRsUlWp4i92zoedV5E7Az8MTOfCbwKeGxE7DfD6/bSrba8iNgR+F1mPi0zXwwcGBH7l+mDfB1qzcvMX1Od0vfzwKsi4imdOX1+v9WWFxFbArsBP4qINwKnRcTREbHngOut1rziwoh4N/CPwGsi4qlDZNWSFxF7l085DHgJ1ffuvTqejz7XW215EbGQ6nTS78nMfwQeFhEHRsSDynLe2We3WvOofl++G3gN1ZnQnhwRjxlkvdWdFxFvpxr8nBkRa2V1CYKFHc8vKMt8LfBL4GKqYz8azSo2pL6/WbXmRcRbgAdS/U56aUQc1PV8P+uttqwy/xuB+1G9GfCUiNh16rnMzD671ZbV4Q9UZwE8OiL26Xq9hX3m1ZYVEe+k2tp1fETsFhG7D7qsdWaVvIcDv8/MF0fEX5e/zZuXjDv7WdY6s0rem8vdvam+R5Z2Pd/Peqstq8z/98DOwJFUJyXZY+q58nXoZ73VllXyXgDsTvWm50OAF0TEAzry+l3WWvNU5ASMOPu5AY8AjqJ6d38B8CjgaDreuaR6N+A7wANGmBWUXX6o3qU+m7KZvCvvsFHnUf1zsHnH4+cCjyj37zFAt9rygHtQ3jkH/orqOJP3AScA9xng61BbHrB2x/fJU4F7Ur37+m46donro1vdeQ8EzqfahfBeVPuxH03HWfd6zao7D9ix5PxTefxIqi0S9x2wWy15wMYd97cA3jY1Px1bcHr53q07D1gIPL18r55MdZ20Y6m2SOw5QLe68/6W6p3gzcrX4+PAi7vm6Smrzjyq3917lvuvBj7L6p+17q2XTwa+ywy7MNWctUH5uC/Vz9FQf2cayFsfeErH8j0aeNeAy1pbVplnPeCgjsfvpBpEPpvq92Y/3WrL6spdF3g+1e+ij1H9Y7ppeS76yasrC1gLOJxqa8cFwAep9jj4u7Ie+llvtWV1zLc51e+fU6m2fH0E+HtglwGWtc6sJcAzOh6/hnKsK9Xv0X7WW21ZHfMdCLy03H8D1V5Vz6ZrS3CPy1pbVplvL+CIcn9T4N+oBpILuesW1rHkeSvra9wFeip5111qtqP6h+WZwFZl2seA3cv9w4DvATs2nTVNXudA55Xlh+hAqn+aY9R5XVlrdUw/Engi1Ts6rxiw21B5XVkLy+c8mNXHKJ7M6n+6+u02VF5n1tSyUu2atBbwHuB3wOGDdBs2b5qsrTvuL6L653S3IboNnNf1Nbg31bFiX2D14OckyqmRB/iaDpXXldX5B+PFVNfpWa9jWr/dhsrrytqUaveUD1OO96D6Z+GwAbsNldeVtRvVm2MfpDrxx7Oofq62HLBbHXkLysfO35VvBD5X7q/D6jfNnkq1a/z9R5B1PNXP+KHARsBrqQa3g/6dqTtvalk3oPo5X1C+Hl8o07dg9Ztmva63obO6f57K452pdpN+PPAZOnaV7rVbHVnT/KyvTzXAOIDqb99FwJsH6TZsVtfXYSOqQdqJ5fFSqjdSDu63Wx1ZXXnrl+/d06gGaxtQ7R30jAG6DZ013fdImfa88nnd3z99fb8Nk9W1rLtR7TnyQeB/gIdT/e/6vAHW29BZXXnLqHYVP5Tq0i9fodpz6WHjzPPWtX7HXWDOgvB24FaqX8xT/yQ8uPxQvwV4GdXoftfy3EOY+R+X2rJmyVu74/kPAzdRTioxyrwZsqb+UTkE+BzwbeCxQ3QbKG+6rK7nl1G9m7P7oN0GzZspq/xy+SjVVsJdgYcM022QvBm+Bp2Zu5TlfMCg3QbNmyHr/lRbHE6kekPg+6wecA/SbaC8HpbzzcDTOx4Pu956zuvI+mxXxnOp3jh5GtW76fv12W3ovK6sqS02m5bbVlQD8C+wenDWT7c68qYGEEE1oO0cYP8D1Vbq/wdsTfXG3veY+R+hOrPeUr43dqfaurZ3+byjGOzvTN15nct6H+46uD2BakD0Xapr9PWz3obKmuHrMPXztX75+DSqXTHXBnbo82s6cNY0eVuX+ztR7anxmrKc/0y1O+b9+ug2VNY0efeiY4+DMv1ZVMd3LuxzvQ2VNU3eRlQnFvoo1cmZ1qI6kdDHqN5MuW8f3YbKmuF7pHMvlrcDL+t43O/vkIGzpsnbmGpr6gOpdsuG6njcr1ANevtZb0NlTZO3uKz3F1INIPcC9geeUObZftR53qZZx+MuMGu5u+9S83nKHxOqfwx2pvqDN7UL4WxnHqota4a8zn9gpnLfT3kHbJR5s2WVaQ8E/gs4cNhu/ebNsZwLgW2oDhx+XA3rra+8ObKWU21d7VzuYbr1lTdN1uc6shZR/QL8GuWd6gG6DZw3TVb3z9YDqHav2WvAbgPnzfE1mPrn7zA6dgcd4mvaV94cX4MHUf0D8xEG/zkdOG+Gr0Hn9+pi4AxWbyG82xllm8zryv001W5crwEeWaZvC6yi/NyXafdsOou77874GFbvzngv+v87U2veLMu6b5n2U+AKypt2A663gbOmyTsO2LtMfxhwIeWNigG6DZU1Td6rgZdTnXDrNKrBxjoDdhsqqyPvrI68fcr0h1MdzrBvn91qyerI++eS9yKqLXKvojo29rvAo/vsVkvWDN+/+5Xp+1MdKrBun91qyZpmWV9D9cbOR8vyXgQ8ZsD1NlTWNHmvYPXv9x2Ay+i45Ms48rx1rd9xF5jtC18+du9S8/nu6aPMmiNvanefdcvHqS1fwSz/vNSZ10sW1T8J209lDbusveb1kLV2yZrxFNRN5c31PVIe3223vVHk9bCca1Ed9zS1S9ewX9Oe8+ZaTur72arz577756rngUVdeb18f9S53mrK+svujOXj1BaSYb/f+srryjyO6h+XbYD/APakepPnAaz+x77Xn6+hsjpyundn/FKZfu+5lqvJvDmWdW+qvyvvYPU/qL1+XYfO6uHrsCHVP6WPGrJb31lz5N2PssdIL99rdWf1sKwbUW2ZfnQN662vrBnylpa8XUrexsBOA3YbOGuuZS3T15vKG2VWD98jD6LaI+1hQ35N+86aIW/bkrcX1f8O+7H6jbKR53mbYT2Pu8AMX/xed6np5Z/42rImPa+PrPvU3G3OvDW82wUla8aLZjeV18dyblnzepszr4+sutZbz3lreLep74+6fof0nDfJ660zr+Px/lTHn5xDdUKmbamOwe7chXWmN8Zqz6KH3RkHWG+15c2xrH/D6pOAzfYGT21ZfeQ9kXIJolFl9ZB3cMf3SDD3mwq1ZfWxrE+iazfMprN6yHt8yXv2qLP6WNZnz/V9W3dWH8v6zBq69ZXVx7I+nVku89JknrfZbwuYQFmd6noB1QVhXw28LCIeWZ7+ItU+xe/KzKtHmTXpeX1k/armbnPmreHddihZv6ypW895fSzntTV16zmvj6y61lvPeWt4t6nvj7p+h/ScN8nrLSIWTOVFxFsi4mFUW7weTrW74E3AJ4DrM/P2jg7ZZNY0y/ka4B8iYt8y7XCq3X3fmpm/mGs5687rcVnPBG7JzFWzLWedWX3kfZLqcjo31rCcPWX1mHcDq79HMjPvHEVWn8t6c2beNKqsHvN+U5b1ulFm9bGsnwCum+37tu6sPpb1k5TvkVFl9bGsnwRuyLkvoVR7nuYWPXydR6p8E9wZEcdRXZDy01THnryCah/7HanelftBxOprMTWdNel5drOb3ezW1m4dmUF1TMci4A6q40o3LE/vDZydmd+aK6fOrFmW85XAv1OdnOFrmfmdHnvVmlcyJ2692c1udpuf3ZrI0xxyAjYvlr/xfe9SM4qsSc+zm93sZre2diuf33ldqqcDX6Y6PmdbqpMzPIWu66WNKGscuzP2kzeR681udrPb/OvWRJ633m8TsRtnj5t077ZLTdNZk55nN7vZzW5t7VbyFubUfwkRC4GLqY5VewbV9Sk/S3XSkvWn5htR1jh2Z+wnbyLXm93sZrf5162JPPVnYnbjrHOTbt2bhyc5z252s5vd2tqt5C2guqD4WlSnpl+b6t3fxVQDx9uyx+M6as5yvQ2QZTe72W3+dWsiT71bNM4Xj7jLcRqHs/rMRRtS7ebzX8C5wD/38K5GbVmTnmc3u9nNbm3tVvK2Ba4u876L6p3fjwOPojqZxaVU13PbMOc48UzNWa63AbLsZje7zb9uTeRpQDnCfUY7b3Qcn0B1LaP7U1375khgE6prgLyZPk4BX0fWpOfZzW52s1tbu5WM46i2ak1do+pw4Anl/lKqs7RtQA/XXKw5y/U2QJbd7Ga3+detiTxvg9/Guhunm6/bt6x2s5vd7DZoXkS8C1gX+AXVteTeR3V9vkdSnZFybeB44IjM/J9RZdW9nHXnTfJ6s5vd7Da/ujWRp+EsGPULRsS25dgEWL1J92Sq3VfWotqkuxbVJt1Z/8jVmTXpeXazm93s1tZuJe95wL0z88jMfBdwNvB84L+pDubfi+pafSf18I9QnVmutwGy7GY3u82/bk3kqQb9bAYc9oabr1u3rHazm93sNmxe+bwXA4eX+wvKx78BzgAeWh4vHmWW623wLLvZzW7zr1sTed6Gv41sy17ZpHsf4N+AF0XEQ4DfA4+MiD2Bbcrzm2fmbaPKmvQ8u9nNbnZra7cuK4FdI2LzzLwTIDO/WLL/PiLuCfS6i+TQWa63obPsZje7zb9uTeRpSCM5G2fHJt2/LY8Ppdqk+z6qdzH3Ah5Nf5uvh86a9Dy72c1udmtrt2lcRnWJgb0iYkVmXhMRDwW2KHm/HVWW662WLLvZzW7zr1sTeRpWjmDzIW6+bt2y2s1udrPbsHnT5G8L/D3wz8AxwArggFFnud7qybKb3ew2/7o1kedtuNtIzsYZEfsAjwPem5nXd0z/ALAR8BLg5uyhTJ1Zk55nN7vZzW5t7TbLa6wF7Eh1+YGbM/PSUWe53urLspvd7Db/ujWRp8GN6pi9y4DfUm3S3QqgY5PuBzLzt338kasza9Lz7GY3u9mtrd2mlZm3Z+YVmfm9Yf85GCLL9VZTlt0mI89u48+aT92ayNPgRnadvYjYFjgM2JPq1KtPAY7NzG+OM2vS8+xmN7vZra3dJpnrTZLUBiO9qLqbr8efZTe72c1uk9RtkrneJElrupEO9iRJkiRJozGy6+xJkiRJkkbHwZ4kSZIktZCDPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmS5q2IeG1EXBERl0XEJRHxsAZf6zsRsbypfEmSui0adwFJksYhIh4OHAzsnpm3RcRmwNpjriVJUm3csidJmq+2AG7IzNsAMvOGzPxVRLwhIi6MiMsj4tSICPjLlrmTI+K8iPhxRDw0Ir4YEf8VEW8u8yyNiJ9ExMfL1sLPR8Ti7heOiAMj4vyIuDgiPhcR6490ySVJ84KDPUnSfPUNYOuI+FlEfDAiHlmmvz8zH5qZuwDrUm39m/KnzNwX+BDwL8BLgV2A50TEpmWenYBTM3M34GbgJZ0vWrYgvg54TGbuDqwAXtnMIkqS5jMHe5KkeSkzbwEeAhwBrAI+ExHPAR4VEf8RESuB/YGdOz7ty+XjSuCKzLyubBm8Eti6PPfLzPxBuX8msE/XS+8JPBD4QURcAjwb2LbOZZMkCTxmT5I0j2XmHcB3gO+Uwd0Lgd2A5Zn5y4g4AVin41NuKx/v7Lg/9Xjqb2p2v0zX4wC+mZmHD70AkiTNwi17kqR5KSJ2ioj7dUxaBvy03L+hHEf3pAGityknfwE4HPh+1/MXAHtHxH1Lj8URseMAryNJ0qzcsidJmq/WB94XERsBfwb+m2qXzv+j2k3zKuDCAXJ/DDw7Ij4M/BdwSueTmbmq7C56VkTco0x+HfCzAV5LkqQZRWb33iWSJGkQEbEU+NdychdJksbK3TglSZIkqYXcsidJkiRJLeSWPUmSJElqIQd7kiRJktRCDvYkSZIkqYUc7EmSJElSCznYkyRJkqQWcrAnSZIkSS30/wGsctti0/h+CgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define the size of the figure\n",
    "plt.figure(figsize=(15, 5))\n",
    "ax = plt.gca()\n",
    "\n",
    "#dot plot of the mean length by file\n",
    "df1.plot(x='sample', y='len',kind=\"scatter\", color='blue', ax=ax, linewidth=3)\n",
    "df1.plot(x='sample', y='len_mean',kind=\"scatter\", color='forestgreen', ax=ax, linewidth=3)\n",
    "plt.title('Basic stats about read length', fontsize=10)\n",
    "plt.xlabel('Sample')\n",
    "plt.ylabel('Read Length')\n",
    "plt.xticks(\n",
    "        rotation=45,\n",
    "        horizontalalignment='right',\n",
    "        fontweight='light',\n",
    "        fontsize='8',\n",
    "    )\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### EXERCICE\n",
    "\n",
    "* 1 - Plot the total_base column (barplot)\n",
    "* 2 - Plot the qual-mean (dot plot)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJkAAAF+CAYAAAAhs8u2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABD30lEQVR4nO3deZgcZbn38e+dBSIkgEJUJKyKC4QQMMQgiCiIgCCCqCCiuEUWd0UjR8XtuHIEESWionJEBBURj+AuCB62gAFkUVFfJciRCEIIixpyv39UjWnGmWSSqt6qv5/r6ivd1TX1e+6qzKRzz1NVkZlIkiRJkiRJVYzr9gAkSZIkSZLU/2wySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkqZKI2DAiFpaP/4uI21perzVs3TdHxDpj2OZFETFrhOX/LyI2qnP8gyAitoiIX42wfPeI+J9ujGl1RcSXI+Lgbo9DkiSNbkK3ByBJkvpbZt4JzASIiPcBSzPzhFFWfzPwVeD+ToxNkiRJneNMJkmSVLuI2CMifhkR10fE6RGxdkS8EXgc8LOI+Fm53qkRsSAiboiI949x88dGxJXl4wnldvaPiCvKzB9HxGPK5c9smVX1y4iYUi4/NiKuiojrRsqNiPHlzJlflTW8pVz++Ij4fkRcHRGXRMSTy+VbRsRl5TY/GBFLy+UPmykUEadExBHl86dGxMXltn4QERuXyy+KiI+V9f0mIp7RMqYTyvFcFxFvWMV2nhoR10bEZcAxK9mf60XEtyPixoiYHxHjIuLVEXFiy7hfGxGfHOM+em25H66NiG8NzVwr1z01In4WEb8vj83pEXFTRHy5ZbtLI+K/IuKaiPhJREwd4fiMWLMkSequxjeZyg8vd4w0RXyEdU9s+SD6m4i4uwNDlCSpaSYBXwZekpnbUcycPiozTwb+DDwrM59VrvsfmTkLmAE8MyJmjGH7SzJzNnAKcFK57FJgTmbuAHwdeEe5/O3AMZk5E3gG8EBE7AVsDcymmIH11IjYbVjGTGCTzJxe1vClcvlpwBsy86nltj9bLv8UcGpm7gT836oKiIiJwKeBg8ttnQ78Z8sqE8oa3wwcXy6bC2wJ7JCZM4AzV7GdLwFvzMydVzGc2cDbgO2AxwMHUezD55fbB3hlyz4YMpOR99G5mblTZm4P3AS8uuVrHgk8G3gL8F3gRGBbYLuImFmusy5wTWbuCFzcUj8wpn0nSZK6ZBBOl/syxYfQM1a1Yma+Zeh5+dvBHdo3LEmSGms88IfM/E35+isUM2lOGmHdF0fEXIrPJBsD2wDXrWL7Z7X8OTTbZhpwdjmjZS3gD+XyXwCfjIgzKZofi8om017AL8t1JlM0nX7ekvF7YKuI+DTwPeCHETEZeDrwjYgYWm/t8s9dgBeWz/8b+NgqangSMB34Ubmt8cDtLe+fW/55NbBF+XxPYH5mLgPIzLsiYvpI24mI9YENMvPiljHtM8pYrszM3wNExFnArpn5zYj4KbBfRNwETMzM64d93b/to3L59Ij4ELABxb79QcvXfDczMyKuB/4ytM2IuKGscyGwHDi7XP+rLftiyKr2nSRJ6pLGN5ky8+cRsUXrsoh4PPAZYCrFNSFem5k3D/vSQxn2mzNJkjQm941lpYjYkmI20E6Z+bfylKlJY/jSHOH5p4FPZub5EbE78D6AzPxoRHwP2Be4PCL2BAL4SGZ+btSAYjzbA8+laJC9mGJW0d3lrKhVjWvIMh4+c3yovgBuWMkso7+Xfz7Eis9rMULGiNuJiA1GGc9Ihq839PoLwHHAzfz7LKbR9tGrKH7B94LMvLY8NXD3li8bqmt5y/Oh16N9Lh1TzZIkqfsaf7rcKEab6g5ARGxOMR39p10YmyRJ/W4SsEWU10sCDqc47QngXmBK+Xw9iobUPVFcQ2m0mTbDvaTlz8vK5+sDt5XPXzG0YkQ8PjOvz8yPAQuAJ1PMrHlVOTOJiNgkIh7dGhDFHezGZea3gPcAO2bmEuAPEfGicp0omyxQzJg6pHx+WMum/ghsE8U1qdYH9iiX/xqYGhE7l9uaGBHbrqLuHwJHRsSE8mseNdp2MvNuiv266whjGm52FNeUGkexTy8FyMwrgE2Bl7Ji9thK91H51hSK2VQTV5E7mnHA0F3kXjo0nhZrsu8kSVIHNH4m03CrmOo+5BDgm5n5UCfHJklSQzxIcQ2fb5QNkauA+eV7pwEXRsTtmfmsiPglcAPFqVe/GOP2146IKyiaEYeWy95X5t0GXE7xyyKAN0fEsyhmBN0IXJiZf4+IpwCXlZ8FlgIvA+5oydgE+FLZeAF4V/nnYcCpEfFuYCLFtYuuBd4EfC0i3gR8a2gjmXlrRJxDcQrgbylP0cvMf0TEwcDJZfNpAsXphDespO4vAE8ErouIfwKfz8xTVrKdVwKnR8T9PPyUteEuAz5KcU2mnwPfbnnvHGBmZv5thK8bbR+9B7iCosF2PSuaimN1H7BtRFwN3MOKpiKwxvtOkiR1QGSOdSZ1/ypPl/ufzJweEesBv87MUe9CUn7gPSYz/7dTY5QkSc0REUszc3K3x1FVFHfGOzEzf9LBzEbsO0mSBtHAnS63iqnuRMSTKO58ctkom5AkSWq0iNggIn4DPNDJBpMkSepvjZ/JVN4lZXdgI+AvFBfz/ilwKsVdbCYCX8/MD5Trvw+YlJnzujFeSZIkSZKkftT4JpMkSZIkSZLab+BOl5MkSZIkSVL9bDJJkiRJkiSpsgndHkA7bbTRRrnFFlt0exiSJEmSJEmNcfXVV/81M6cOX97oJtMWW2zBggULuj0MSZIkSZKkxoiIP4603NPlJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmVNfqaTCP55z//yaJFi3jwwQe7PZSBMmnSJKZNm8bEiRO7PRRJkiRJktQGA9dkWrRoEVOmTGGLLbYgIro9nIGQmdx5550sWrSILbfcstvDkSRJkiRJbTBwp8s9+OCDbLjhhjaYOigi2HDDDZ09JkmSJElSgw1ckwmwwdQF7nNJkiRJkpptIJtMkiRJkiRJqtfAN5ki6n2srg9/+MOrXOfuu+/ms5/97Ji2N3ny5FHfu+iii9hvv/3GPDZJkiRJkqSxGvgmU7fV3WSSJEmSJEnqBptMXTRv3jweeOABZs6cyWGHHQbAJz/5SaZPn8706dM56aST/rXe7373O2bOnMmxxx7L0qVL2WOPPdhxxx3Zbrvt+M53vjPmzCVLlnDggQeyzTbbcOSRR7J8+XIAjjrqKGbNmsW2227L8ccf/7AxbrPNNsyYMYO3v/3tACxevJgXvvCF7LTTTuy000784he/AODiiy9m5syZzJw5kx122IF77723jt0kSZIkSZL6QGRmt8cAQESMBxYAt2XmfsPeC+BTwL7A/cARmXnNqrY5a9asXLBgwcOW3XTTTTzlKU9p2Xb1sbda3d05efJkli5dCsDVV1/NEUccweWXX05m8rSnPY2vfvWrPPKRj2S//fbjV7/6FQDLli3j/vvvZ7311uOvf/0rc+bM4be//S0R8bDtDXfRRRex9957c+ONN7L55puz995787rXvY6DDz6Yu+66i0c96lE89NBD7LHHHpx88slMmzaNnXfemZtvvpmI4O6772aDDTbgpS99KUcffTS77rorf/rTn3juc5/LTTfdxP7778+8efPYZZddWLp0KZMmTWLChAn/yh++7yVJkiRJGm5N/5++pu2NTuc1QURcnZmzhi/vpZlMbwJuGuW9fYCty8dc4NRODaqTLr30Ug488EDWXXddJk+ezEEHHcQll1zyb+tlJscddxwzZsxgzz335LbbbuMvf/nLmDJmz57NVlttxfjx4zn00EO59NJLATjnnHPYcccd2WGHHbjhhhu48cYbWW+99Zg0aRKvec1rOPfcc1lnnXUA+PGPf8zrX/96Zs6cyfOf/3yWLFnCvffeyy677MJb3/pWTj75ZO6+++6HNZgkSZIkSVKz9USTKSKmAc8DvjDKKgcAZ2ThcmCDiNi4YwPskLHOKjvzzDNZvHgxV199NQsXLuQxj3kMDz744Ji+Noa1aCOCP/zhD5xwwgn85Cc/4brrruN5z3seDz74IBMmTODKK6/khS98Ieeddx577703AMuXL+eyyy5j4cKFLFy4kNtuu40pU6Ywb948vvCFL/DAAw8wZ84cbr755tXbAZIkSZIkqW/1RJMJOAl4B7B8lPc3AW5teb2oXNb3Jk6cyD//+U8AdtttN8477zzuv/9+7rvvPr797W/zjGc8gylTpjzs+kb33HMPj370o5k4cSI/+9nP+OMf/zjmvCuvvJI//OEPLF++nLPPPptdd92VJUuWsO6667L++uvzl7/8hQsvvBCApUuXcs8997Dvvvty0kknsXDhQgD22msvTjnllH9tc2j57373O7bbbjve+c53MmvWLJtMkiRJkiQNkK6fzxQR+wF3ZObVEbH7aKuNsGzEaT8RMZfilDo222yzVeZ3+xzKuXPnMmPGDHbccUfOPPNMjjjiCGbPng3Aa17zGnbYYQcAdtllF6ZPn84+++zDO9/5Tvbff39mzZrFzJkzefKTnzzmvJ133pl58+Zx/fXXs9tuu3HggQcybtw4dthhB7bddlu22mordtllFwDuvfdeDjjgAB588EEykxNPPBGAk08+mWOOOYYZM2awbNkydtttN+bPn89JJ53Ez372M8aPH88222zDPvvsU/PekiRJkiRJvarrF/6OiI8AhwPLgEnAesC5mfmylnU+B1yUmWeVr38N7J6Zt69s22O58Lc6x30vSZIkSVoVL/zd+3r2wt+Z+a7MnJaZWwCHAD9tbTCVzgdeHoU5wD2rajBJkiRJkiSpc7p+utxoIuJIgMycD1wA7AvcAtwPvLKLQ+t5119/PYcffvjDlq299tpcccUVXRqRJEmSJElqup5qMmXmRcBF5fP5LcsTOKY7o+o/22233b8uxi1JkiRJktQJXT9drhu6fR2qQeQ+lyRJkiSp2QauyTRp0iTuvPNOmx4dlJnceeedTJo0qdtDkSRJkiRJbdJTp8t1wrRp01i0aBGLFy/u9lAGyqRJk5g2bVq3hyFJkiRJktpk4JpMEydOZMstt+z2MCRJkiRJkhpl4E6XkyRJkiRJUv1sMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqrOtNpoiYFBFXRsS1EXFDRLx/hHV2j4h7ImJh+XhvN8YqSZIkSZKkkU3o9gCAvwPPzsylETERuDQiLszMy4etd0lm7teF8UmSJEmSJGkVut5kyswElpYvJ5aP7N6IJEmSJEmStLq6frocQESMj4iFwB3AjzLzihFW27k8pe7CiNi2syOUJEmSJEnSyvREkykzH8rMmcA0YHZETB+2yjXA5pm5PfBp4LzRthURcyNiQUQsWLx4cbuGLEmSJEmSpBY90WQakpl3AxcBew9bviQzl5bPLwAmRsRGo2zjtMyclZmzpk6d2uYRS5IkSZIkCXqgyRQRUyNig/L5I4A9gZuHrfPYiIjy+WyKcd/Z4aFKkiRJkiRpFF2/8DewMfCViBhP0Tw6JzP/JyKOBMjM+cDBwFERsQx4ADikvGC4JEmSJEmSekDXm0yZeR2wwwjL57c8PwU4pZPjkiRJkiRJ0th1/XQ5SZIkSZIk9T+bTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIq63qTKSImRcSVEXFtRNwQEe8fYZ2IiJMj4paIuC4iduzGWCVJkiRJ6raINXv0S57614RuDwD4O/DszFwaEROBSyPiwsy8vGWdfYCty8fTgFPLPyVJkiRJktQDuj6TKQtLy5cTy0cOW+0A4Ixy3cuBDSJi406OU5IkSZIkSaPrepMJICLGR8RC4A7gR5l5xbBVNgFubXm9qFwmSZIkSZKkHtATTabMfCgzZwLTgNkRMX3YKiOdzTl8tlOxYsTciFgQEQsWL15c80glSZIkSZI0kp5oMg3JzLuBi4C9h721CNi05fU04M+jbOO0zJyVmbOmTp3ajmFKkiRJkiRpmK43mSJiakRsUD5/BLAncPOw1c4HXl7eZW4OcE9m3t7ZkUqSJEmSJGk0vXB3uY2Br0TEeIqm1zmZ+T8RcSRAZs4HLgD2BW4B7gde2a3BSpIkSZIk6d91vcmUmdcBO4ywfH7L8wSO6eS4JEmSJEmSNHZdP11OkiRJkiRJ/c8mkyRJkiRJkiqzySRJkiRJkqTKbDJJkiRJkiSpMptMkiRJkiRJqswmkyRJkiRJkiqzySRJkiRJkqTKbDJJkiRJkiSpMptMkiRJkiRJqswmkyRJkiRJkiqzySRJkiRJkqTKbDJJkiRJkiSpMptMkiRJkiRJqswmkyRJkiRJkiqzySRJkiRJkqTKbDJJkiRJkiSpMptMkiRJkiRJqswmkyRJkiRJkiqzySRJkiRJkqTKbDJJkiRJkiSpMptMkiRJkiRJqswmkyRJkiRJkiqzySRJkiRJkqTKut5kiohNI+JnEXFTRNwQEW8aYZ3dI+KeiFhYPt7bjbFKkiRJkiRpZBO6PQBgGfC2zLwmIqYAV0fEjzLzxmHrXZKZ+3VhfJIkIGLNvi6zP/IkSZIkVdP1mUyZeXtmXlM+vxe4Cdiku6OSJEmSJEnS6mhLkyki1l3Dr9sC2AG4YoS3d46IayPiwojYdiXbmBsRCyJiweLFi9dkGJIkSZIkSVpNtTaZIuLpEXEjxWwkImL7iPjsGL92MvAt4M2ZuWTY29cAm2fm9sCngfNG205mnpaZszJz1tSpU9ekDEmSJEmSJK2mumcynQg8F7gTIDOvBXZb1RdFxESKBtOZmXnu8Pczc0lmLi2fXwBMjIiN6hy4JEmSJEmS1lztp8tl5q3DFj20svUjIoAvAjdl5idHWeex5XpExGyKcd9Zw3AlSZIkSZJUg7rvLndrRDwdyIhYC3gj5alzK7ELcDhwfUQsLJcdB2wGkJnzgYOBoyJiGfAAcEim9w+SJEmSJEnqFXU3mY4EPkVxd7hFwA+BY1b2BZl5KbDSG1Vn5inAKTWNUZIkSZIkSTWrtcmUmX8FDqtzm5IkSZIkSep9dd9d7uMRsV5ETIyIn0TEXyPiZXVmSJIkSZIkqffUfeHvvTJzCbAfxelyTwSOrTlDkiRJkiRJPabuazJNLP/cFzgrM+8qbwonSY23pj/u1vQ2Bp3OkyRJklRdk//fUHeT6bsRcTPFHeCOjoipwIM1Z0iSJEmSJKnH1H3h73kR8TFgSWY+FBH3AQfUmSFJkjRck38jKElaff67IHVH3TOZADYBnhMRk1qWndGGHEmSJEmSJPWIWptMEXE8sDuwDXABsA9wKTaZJEmSJEmSGq3uu8sdDOwB/F9mvhLYHli75gxJkiRJkiT1mLqbTA9k5nJgWUSsB9wBbFVzhiRJkiRJknpM3ddkWhARGwCfB64GlgJX1pwhSVLtvECoJEmSVE3dd5c7unw6PyK+D6yXmdfVmSFJkiRJqsZfrkhqh9rvLhcRBwG7Aklx0e+eazL5A1WSNGj8t0+SJEntVvfd5T4LPAE4q1z0uojYMzOPqTNHUv/yP7qSJEmS1Ex1z2R6JjA9s/jvWUR8Bbi+5gytgv+J728eP0mSJElSP6q7yfRrYDPgj+XrTenB0+UkSZLUG/zliiRJzVFLkykivktxDab1gZsi4sry9dOA/60jQ5IkSZIkSb2rrplMJ9S0HUmSJEmSJPWhWppMmXlxHduRJEmSJElSfxrX7QFIkiRJkiSp/9lkkiRJkiRJUmU2mSRJkiRJklRZXXeXu57ibnL/9haQmTmjjhxJkiRJkiT1prruLrffmn5hRGwKnAE8FlgOnJaZnxq2TgCfAvYF7geOyMxr1ny4kiRJkiRJqlNdd5f7Y4UvXwa8LTOviYgpwNUR8aPMvLFlnX2ArcvH04BTyz8lSZIkSZLUA2q9JlNEzImIqyJiaUT8IyIeioglK/uazLx9aFZSZt4L3ARsMmy1A4AzsnA5sEFEbFzn2CVJkiRJkrTm6jpdbsgpwCHAN4BZwMuBJ4z1iyNiC2AH4Iphb20C3NryelG57PYRtjEXmAuw2WabjX3kkiRJPShizb4uR7paZg/mNV3Tj595/Z0nSXWr/e5ymXkLMD4zH8rMLwHPGsvXRcRk4FvAmzNz+OynkX7cjvijNDNPy8xZmTlr6tSpqzN0SZIkSZIkraG6ZzLdHxFrAQsj4uMUM43WXdUXRcREigbTmZl57girLAI2bXk9DfhzDeNVH/I3SpIkSZIk9Z66ZzIdXm7z9cB9FI2hg1b2BeWd474I3JSZnxxltfOBl0dhDnBPZv7bqXKSJEmSJEnqjrpnMr0gMz8FPAi8HyAi3gR8aiVfswtFc+r6iFhYLjsO2AwgM+cDFwD7ArcA9wOvrHnckiRJkiRJqqDuJtMr+PeG0hEjLPuXzLyUka+51LpOAsdUHZwkSZIkSZLao5YmU0QcCrwU2DIizm95az3gzjoyJEmSJEmS1Lvqmsn0vxQX+d4I+K+W5fcC19WUIUmSJEmSpB5VS5MpM/8I/BHYOSIeA+xUvnVTZi6rI0OSJEmSJEm9q9a7y0XEi4ArgRcBLwauiIiD68yQJEmSJElS76n7wt/vBnbKzDsAImIq8GPgmzXnSJIkSZIkqYfUOpMJGDfUYCrd2YYMSZIkSZIk9Zi6ZzJ9PyJ+AJxVvn4JcGHNGZIkSZIkSeoxtTaZMvPYiDgI2BUI4LTM/HadGZIkSZIkSeo9tTaZIuJjmflO4NwRlkmSJEmSJKmh6r5e0nNGWLZPzRmSJEmSJEnqMbXMZIqIo4Cjga0i4rqWt6YAv6gjQ5IkSZIkSb2rrtPlvkZxge+PAPNalt+bmXfVlCFJkiRJkqQeVUuTKTPvAe4BDq1je5IkSZIkSeovdV+TSZIkSZIkSQPIJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIqs8kkSZIkSZKkymwySZIkSZIkqTKbTJIkSZIkSarMJpMkSZIkSZIq64kmU0ScHhF3RMSvRnl/94i4JyIWlo/3dnqMkiRJkiRJGt2Ebg+g9GXgFOCMlaxzSWbu15nhSJIkSZIkaXX0xEymzPw5cFe3xyFJkiRJkqQ10xNNpjHaOSKujYgLI2Lb0VaKiLkRsSAiFixevLiT45MkSZIkSRpY/dJkugbYPDO3Bz4NnDfaipl5WmbOysxZU6dO7dT4JEmSJEmSBlpfNJkyc0lmLi2fXwBMjIiNujwsSZIkSZIklfqiyRQRj42IKJ/Pphj3nd0dlSRJkiRJkob0xN3lIuIsYHdgo4hYBBwPTATIzPnAwcBREbEMeAA4JDOzS8OVJEmSJEnSMD3RZMrMQ1fx/inAKR0ajiRJkiRJklZTX5wuJ0mSJEmSpN5mk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSJEmSJFVmk0mSJEmSJEmV9USTKSJOj4g7IuJXo7wfEXFyRNwSEddFxI6dHqMkSZIkSZJG1xNNJuDLwN4reX8fYOvyMRc4tQNjkiRJkiRJ0hj1RJMpM38O3LWSVQ4AzsjC5cAGEbFxZ0YnSZIkSZKkVemJJtMYbALc2vJ6UblMkiRJkiRJPaBfmkwxwrIcccWIuRGxICIWLF68uM3DkiRJkiRJEvRPk2kRsGnL62nAn0daMTNPy8xZmTlr6tSpHRmcJEmSJEnSoOuXJtP5wMvLu8zNAe7JzNu7PShJkiRJkiQVJnR7AAARcRawO7BRRCwCjgcmAmTmfOACYF/gFuB+4JXdGakkSZIkSZJG0hNNpsw8dBXvJ3BMh4YjSZIkSZKk1dQvp8tJkiRJkiSph9lkkiRJkiRJUmU2mSRJkiRJklSZTSZJkiRJkiRVZpNJkiRJkiRJldlkkiRJkiRJUmU2mSRJkiRJklSZTSZJkiRJkiRVZpNJkiRJkiRJldlkkiRJkiRJUmU2mSRJkiRJklSZTSZJkiRJkiRVZpNJkiRJkiRJldlkkiRJkiRJUmU2mSRJkiRJklSZTSZJkiRJkiRVZpNJkiRJkiRJldlkkiRJkiRJUmU2mSRJkiRJklSZTSZJkiRJkiRVZpNJkiRJkiRJldlkkiRJkiRJUmU2mSRJkiRJklRZTzSZImLviPh1RNwSEfNGeH/3iLgnIhaWj/d2Y5ySJEmSJEka2YRuDyAixgOfAZ4DLAKuiojzM/PGYatekpn7dXyAkiRJkiRJWqVemMk0G7glM3+fmf8Avg4c0OUxSZIkSZIkaTX0QpNpE+DWlteLymXD7RwR10bEhRGxbWeGJkmSJEmSpLHo+ulyQIywLIe9vgbYPDOXRsS+wHnA1iNuLGIuMBdgs802q3GYkiRJkiRJGk0vzGRaBGza8noa8OfWFTJzSWYuLZ9fAEyMiI1G2lhmnpaZszJz1tSpU9s1ZkmSJEmSJLXohSbTVcDWEbFlRKwFHAKc37pCRDw2IqJ8Ppti3Hd2fKSSJEmSJEkaUddPl8vMZRHxeuAHwHjg9My8ISKOLN+fDxwMHBURy4AHgEMyc/gpdZIkSZIkSeqSrjeZ4F+nwF0wbNn8luenAKd0elySJEmSJEkam144XU6SJEmSJEl9ziaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMp6oskUEXtHxK8j4paImDfC+xERJ5fvXxcRO3ZjnJIkSZIkSRpZ15tMETEe+AywD7ANcGhEbDNstX2ArcvHXODUjg5SkiRJkiRJK9X1JhMwG7glM3+fmf8Avg4cMGydA4AzsnA5sEFEbNzpgUqSJEmSJGlkvdBk2gS4teX1onLZ6q4jSZIkSZKkLpnQ7QEAMcKyXIN1ihUj5lKcUgewNCJ+vZrj2Qj468jbXs0tmWeeeeZ1J8s888wbnLwm12aeeeZ1L6/JtZlnnnn15G0+0sJeaDItAjZteT0N+PMarANAZp4GnLamg4mIBZk5a02/3jzzzDOv21nmmWfe4OQ1uTbzzDOve3lNrs0888xrb14vnC53FbB1RGwZEWsBhwDnD1vnfODl5V3m5gD3ZObtnR6oJEmSJEmSRtb1mUyZuSwiXg/8ABgPnJ6ZN0TEkeX784ELgH2BW4D7gVd2a7ySJEmSJEn6d11vMgFk5gUUjaTWZfNbnidwTIeGs8an2plnnnnm9UiWeeaZNzh5Ta7NPPPM615ek2szzzzz2pgXRf9GkiRJkiRJWnO9cE0mSZIkSZIk9TmbTJIkSZIkSapsYJtMEREdyhkXEROHnncis8yyvnpyOl5fk2srs6yvnhzra09uY+trcm1lVqPr64ZO7dNOZg3K8evksesG6+tvTa6vybWB9fW7jv1bO0jXZIqINwC/zswflq8j27gDIuJUirvhbQgcm5mL25VV5llfvXkdq6/JtZV51ldvnvXVm9fY+ppcW5nX6PrKzEdl5l3tzmnJ2x24Hfh/mfn3Nmc1+vg1+diVedZXb97uWF9dWY2trcyzvnrzdsf6atXY3+AMFxHHAa8Cdo2I/aG4a127unkRcTywDvAu4ErgkHbktORZX715HauvybWVedZXb5711ZvX2PqaXFuZ1+j6yswvAO+KiO3bnVXmfRJ4O/Bi4PiIGN/GrEYfvyYfuzLP+urNs776shpbW5lnffXmWV8bDEyTCbgY2B/4HjAnIp4PxQcaaMvUsZsofkv2j/L5pq1vRv1TtK2vXp2sr8m1gfVZ3+qxvvo0uTZoeH0R8WxgkzJr14iY0fJe7Y2YiHgxsGlm7gecCjwELK87p0Vjj1/Tj5311Z5nffVlNba2Ms/66s2zvjZp/OlyETEhM5e1vH4k8FTgWcCFwO+AJZl5XzvyymXrAx/LzCMj4knAnzPz3nbkWV+9eeWyttTX5NpGyrO+evPKZdZXU16T6mtybSPlNa2+YTlrAxtRNGPWBi7OzIV157TkbZyZt5fPvw18NDOviKjvNLZBOX5NPHbD8qyv3jzrqy+rsbWVGdZXb571tUGjZzJFxPjMXBaF/4iIyZn5N+Bq4BvA0cD5wJZtyDsuIiaXbz0C+EdE7Aec3KY866s3r631Nbm2EfKsr94866s3r1H1Nbm2EfIaV1+Z2frZ7B+ZeRvwY+BBYEZEvDgiXh5Rz289h+X9X8vzm4HFEbEl8Ko68pp+/Jp87EbIs75686yvvqxG1TZCnvXVm2d97ZCZjX5QNNKOBX4EfA2YXC6fBfwS2KfNeVPK5ecAVwF7WJ/1Nb0267M+6+tefU2uren1ARPLPwN4N/BxYFy57FHAJ4DfA89oU97HWvIOB+YDPwF29/gN9rGzPuvr1fqaXJv1WV+v1zfao5EzmSLiNS2dubcC61FMSbsB+HJErAPcDczNzAurdvFWkXd6RKwF/Ar4UGb+pErWGPKsr968Wutrcm1jyLO+evOsr968vq6vybWNIa/v6xuSmf8sc8+gmG0zGfhA+VvJvwGbA0dl5iVtypsCfLDMGwfMBd6bmRdVyRmE49fUY7eSPOurN8/66stqTG2j5FlfvXnW1w7t7GB14wGcAPx3y+tjgBNbXn8VOB1Yp3w9oc15X6PoGE7qUJ719Wh9Ta7N+qzP+rpXX5NrG4T6ym18GDisfP4M4Bvl86Hfcr67fL1B+We0Me8E4O3l65k11Nbo49fkY2d91tfL9TW5Nuuzvl6vb1WPRs1kioi5wLOBV5evpwOfA8ZHxAciYgdgMXA78N2IWDeHXQyyDXl/Ae4EvtehPOvrwfqaXJv1WZ/1da++Jtc2CPWVGZ8GNsnMM8tFVwH3RcQLM/Muihk30yPiPcC9sOIObG3Kux6YHRHHlc8pfyu6JlmNPn5NPnZjzLO+evOsz9rGmmd99eZZX80ac3e5KO5OsiOwPXBl+Xw2xZSwpDgfcTJwVmb+KCJeD5yXmYvMM6+deU2uzTzzzOteXpNrG4S8MvOpwEeAt2Xm9RHxAorp838HXgUsAXYAjgf2BJZl5of6Ia/px6/Jx84883o5r8m1mWder+eNWXZgulS7H8DewGeATYHnUFxQ6ybgScPWW8c88zqZ1+TazDPPvO7lNbm2Qcgrt3Vw+ecrgddTfEj8MSsuTP0o4Ahgdvl6A2Cjfshr+vFr8rEzz7xezmtybeaZ1+t5q/No0kymt1HsuM8Bm1BMl74auDwzl5hnXrfymlybeeaZ1728Jtc2IHlfBf6UmcdFxKHAy4DPZOYFdWd1Ka+xx28Ajp155vVkXpNrM8+8Xs9bHU1qMn0WeCJwH0Un7wnA04HfAedk5nLzzOtGXpNrM88887qX1+TampwXEROyvA5QRFwKXJuZx0TEyyl+63gZcFW/5rXkNu74Nf3YmWder+Y1uTbzzOv1vDUxrlvBdYqIlwGPyMw9gcuB91N8qLgSuKINH2TMM6/nsswzz7zByWtybU3Py8xlETEhIr4FfBvYOiI+mZlnUFwA9DnAuv2aB809fk0/duaZ16t5Ta7NPPN6PW9NB9l3D2D8sNezgeNaXn8bOIny9rTmmdepvCbXZp555nUvr8m1DUJeuc1xLc8nAfNbXp8NnF0+n9ZveU0/fk0+duaZ18t5Ta7NPPN6Pa/Koy9nMmXmQxExLiJOioijgSnArIiYFxEvAiYA38zMB80zr5N5Ta7NPPPM615ek2sbhLyIGJ+Zy6MwHVgOTI2IY8pVPgrMiIitssLd6rqV1+Tj1/RjZ555vZrX5NrMM6/X86rqq2syRURkOeCIOAf4OTAeeCrFjp0DPBa4JjO/b555ncprcm3mmWde9/KaXNsg5A3LHgd8Dbgb+BNwKfBuiotT7wW8KTMv7ae8QTl+TTx25pnXD3lNrs0883o9r5Ls8lSqsT4op0ZTfJiYArwG2IpiSvQB5bLHmGdep/OaXJt55pnXvbwm1zYIeWVW69T2Q4G3A+sAi4AtKa6ZsCUwvd/ymn78mnzszDOvl/OaXJt55vV6Xl2PvjldLjMfKp8eT3F7vq2Asyg+XCwEfgxsap55nc5rcm3mmWde9/KaXNsg5LVMbR8XEa8tF48Dvkhx97OJwGsy8w+Z+at+y2vy8Wv6sTPPvF7Na3Jt5pnX63l16vkmU0QcHRGzyucnAZtn5qnAfwLXUnTyvgi8JzMXmGdep/KaXJt55pnXvbwm1zYIeWXO47O4blAAHwHuB34EPJfizi/XAp8Bbum3vKYfvyYfO/PM6+W8Jtdmnnm9nle77IHpVCt7AI8s/9yRYirYj4HnlcsmA48ANjPPvE7nNbk288wzr3t5Ta5tQPK2Bg4pn28F3AHsW75+PPA+4L3APn2a19jjNwDHzjzzejKvybWZZ16v57Xj0fUBrGTnjm95/jjgBorO3ROAH9a9U80zrxezzDPPvMHJa3Jtg5A3Qv7bge2A3YDLgJ3K5eP6MW+Qjl/Tjp155vVLXpNrM8+8Xs+rdezdHsAoO3To4o4BbFU+35hiithewJMo7iiyoXnmdTKvybWZZ5553ctrcm2DkNea2fJ6P+ArwJOBnYBrKD8g9lte049fk4+deeb1cl6TazPPvF7Pa+cjygJ6ThS36DsHuIvi6ulnUZx/eC7wQeDizLzfPPM6ndfk2swzz7zu5TW5tqbnRcS4LC/OCbyLYqbNnymmub8a+ATF3c8iMy/rt7yhTBp4/Jp+7Mwzr1fzmlybeeb1el679VyTKSIiMzMiPg7cBnwfOJ3iziFfAP6P4tz7K8wzr5N5Ta7NPPPM615ek2sbhLyW3HHAecDNFDdWWQx8E5gGvA54bWbe1295g3D8mnrszDOv1/OaXJt55vV6Xjv1TJMpilv0PdTyej3gHxRTxL4LPJriQlfvqOlDk3nm9VyWeeaZNzh5Ta5tEPLKjHGZubx8Ph14bmb+V0RcRjH7ZgHFB8Y7MvNv/ZTX9OPX5GNnnnm9nNfk2swzr9fzOmVCtwcAKz5YlN27+cAvgbsy8+yIuJfiln07AR+v84OMeeb1UpZ55pk3OHlNrm0Q8srMoantAWwE3JSZv4qIE4DTKH4TeRRwdp0fRDuR1/Tj1+RjZ555vZzX5NrMM6/X8zqpl2YyBfAZilv0fQ94N3Ah8AfgHRQfLH5gnnmdzmtybeaZZ1738ppcW9PzIv51SldQXCvonxSndl0IrAPsD2wPzMvMS/otbyiTBh6/ph8788zr1bwm12aeeb2e13HZA1cfLxtdO1F8kNi1fL058MHy+TrmmdetvCbXZp555nUvr8m1DUJeud0PA58sn+8GnExx97P9gWf1c17Tj1+Tj5155vVyXpNrM8+8Xs/r1KMnTpcDyMyrIuI44I0R8QCwDNglinPyl5pnXrfymlybeeaZ1728Jtc2CHkRMRW4D3hyRGyamT+PiD2AWzPzpn7Pa/Lxa/qxM8+8Xs1rcm3mmdfreR3V7S4X5Sl7La9fBNxIcdu+nfsxDxi3kvf6Pm8VY2lsXpNqG/590LS8ph+/puYBmzY5r2Xbjft3b5DyhmU9GpgLfB44ELgceFo/5w3K8WvisTPPvH7Ia3Jt5pnX63mdenQ+EI4A9gQe27Js3LB19gXOBLbpw7x1R8poUN7h5f5cv2l5FBdW2x94dIdq63Teu4EXtLweqfHTz3knArNXsU4/530GeM6wZcP/c9bPeacAn+7g/uxYHnAc8BJgZodq63TeYcDuwIYty9r572xH81q2OWGU5RtTXDPoAuDFI42nl/MYkM9lnTp2Xchbv+F5mzU8b6um5q3s+7dNtXU6b5+G501teN66w16Pa1Jetx/j6KCIeA+wB/Bc4OURcQhAlldVLy98RWZeQHEBrDdGxCP6KO97wC8j4lFlxrhh7/d73ucp9uWhFB/yh7/ft3ll1jMpPtjuOzynTbV1Mu8xFNe7eHdE7FduO6NFn+edDvw9M6+MiCnD3mtC3quAf2bmjyJi44hYu9x+NiTvJGATYPuIOGCE9/s2L4o7hDwFeCzFv0et741rQ22dzvs48GxgH+AVEbF/uf2h2/H2dV65zWMiYmJmLouI8cPfz8zbKRoi3wBmRsQWQ+Pp9bxo/uey44ETIuLgkd5vw7HrdN4HgVMj4iXDlg/tx37Pey3w+YjYbqT3G5D3X8ArRljerv3ZsbyI+BzwhhGWj6s7q0t5JwPfi/Iz7nANyPsw8LmIOCIiJrUsb9ffzU7nvRv4aEQcFBEblxnLm5LXCzraZAKmAF/OzGOBXwCPj4iDoPjPS2bRtouIQ4FJwFsy84F+yIuIxwLXA1+k+CYZsfHTx3nbAg9m5suAtwHPjYjdR1m3r/Ii4onAvZl5SGYeBewVEc8ul7fj70pH88pt/oXiVpjfBN4WES9uzernvIjYBJgB3BgRH6D4gHhsRMxpzezXvBZXlR8QPwXMG/qQ3+95EbFLudkDgaMpvtcf0/J+9GteFA2CjSgu6Pgp4GkRsVdEbE8xiOU119bRvNIU4L+AeRR3RnlRROxZZtR67LqRFxEfBT4BfDWKxs9D0dL4afkPzG3ArcA1FNdX6Is8YD06+7msY3kR8Z/ANhT/Fh0TEXsPe7/uY9fpvA8AWwMfA14cLY2RzMx+zys9QHGXpWMjYtdh4xnfz3lRNMy3yczjI2JGROw49F6bjl/H8iJiZ+C+zDwqIg4oP+M+utz+8jbsy07nfah8ugvF98IWw96v+9h1Ou/NwLbA6ykuRD176L3y70rd+7PTea8FdgROAp4KvDYintKSV/f+7Ghez8gOTpsCngG8nWJK2DjgWcCxwCNa1jkIuAh4Sj/lAUE5tR14K3Ae5bS/YXkH9mneerSc1gW8CnhG+Xztfs4D1gbGl8/3obgWxKeB9wGPa8PflU7nrVX++QyKU2jWB35A8R+1mf2eV25rG+AyilOgHkNxbvOxtNyRqM/znlhmnVy+fibwHuAJ/Z4HPLLl+cbAR4a2OfR9Uj6v62dLx/KA8cBLy+/tE4FfAu8CPgjMaUNtHc0rt3U4xelPG5V/b74CHDVsnb7Mo/jcMKd8/k7gHFb8fBtHy+mjFNcRupiVTPfvpTxgSvnnbhQ/u9r6OakLeZOBF7fsvz2AT7Tx2HU6b11g75bXH6dobr2C4t/cvs5r2dYjgNdQ/Bv0JYr/oG1Yvhf9mgdMpJil/2GK6698Fvg58MpyX9d9/Dqd92iKf3dOo7hT1heANwPT23HsOpkHTAUOa3k9j/LacRT/Bte9LzuaV25nL+CY8vl7ga+W3+vbDVuvX/OeDswtn28I/JSiwTWeh5/C1pd5vfJofwAcA0wsn29J8YH3ZcC0ctmXgB3L5wcClwBP7NO81ubHW8tvkr0o/gMaDcib2LL89cALKbrNb+rHvGFZ48tt7sCK61ydyIoP/HXX1tG8of0JnFr++UngXuDQBuVt2vJ8AsV/0GY0IY/i1KfjgG+xoiFyAhS3Nm1AXus/skcBF9Jy7nq/5Q3L2pDi1LzPUV5PgeJD4oFtqq3TeTMofpnzWYoLNb+c4mfZJv2a1/p3hIf/O/sB4Bvl80ms+EXPS4DvAU/uhzzgeIqfzS8ANgD+g6J5167PSZ3OG9qXUyh+No8r/958q1y+MSt+0VPnsetoXsvrbYFvA88DzgYOaHmvr/MomndfAJ5D8fnvauBDTcgrvxcOBd5fvt6Comm+X7/mtWRNLr/PP0/RAJpCcVbCYXVldTNv2LJXl9sd/n3St3kUP7+Opvh39nfAzhT/n351Q/JmAvMp/k3aCfguxZkXT+vnvF57tHfj8FHgfop/hIY+IO5QfuP/J8W5sxdTdiopfmNQ5YNFL+St1fL+54C/UV48tyF5Qx9y96c4b/QnwHP7LW+krGHvz6ToNO9YNauX8sofYl+kmDm1HfDUpuQNy51e7s+nNCjvyRS/8X8/RdP1UlY0JZuQ17o/PwS8tOV13+S1ZJ0zLONVFA3yQyh+g7x7zbV1I29o1saG5WMaRYPyW6xo+vRVXrmNoYZAUDTsWhuQb6GYxfgzYFOKX2ZdQrUPoh3Lo/gsdDbF1P2LKE6/2JSicdeOz0mdzmvdl4/j4U2791E0SC4GNm/Dset03ias+Hk2ufzzEIrT2dYCHt/neZuWz59EMTN6Xrkvv05xStvWfZ73GFpm1pbLX05x8d/xbdifbc0blrUBxcX9v0hxTdWJFNcf/RJFw/wJNdfW6bxNePjM9Y8Cb2h53Y7v9U7mPZJihts2FKcsA+xH0RyZ3Ib92em8dcq/H6+jaGw9neJaj88v19mq3/J68dG+Df/7tO9vUv7jS/GhcFuKDxlDp0BVvWJ7t/NaPwAP5Z5C+RuCJuWVy7YBfgvs1W95q9iX44HNgHOBfTtw7DqdN4ti1lvrvu33vG+05E2g+GF9IeVvVxuQN/xn2VMopoE/vSF5rX9fhv4DcyAtp+b1S94q/q5sT/Hh9wu07+dmp/O+Oex7ex3gy6yYNfVvd5Ts5bwRsr9GcerFPOCZ5fLNgcWUP6/LZev3Qx7/fkrXnqw4pesx1P85qaN5K9mXu5XLfg3cQPmLqjYeu07mHQfsUi5/GnAVZUO5IXnvBN4IfIdilspEYFID8s5qydu1XL4zxeUTduvXvDLr62XWkRSzi95GcZ25i4E92lBbp/Nav9d3L5c/m+KSAo9oQN7XW/LeQNG4exvFrL49G5T3JlZ8vng8cB3laYj9mtdrj/ZsdPRp398cvrxheUPT2h9R/jk0Cyeo8OG3F/MoPiBuNZTX7v1ZV94YstYqs0a9hW2f5n2z5fW/nTLUgLyh/TmR4lotQ6dhtPvvZqfyBu1nZy3Ngk7mjeV7oZPHroN5/zqdq/xzaIZDu7/3as0bIfc4ig+9mwFXAHMofinwFFb8J7vOn2ltzWvJGX5K17fL5Y9t09+XjuStYl/uQvG56GOs+E9and/z3c6bQ3Edy7ex4pTmJuVtTTnLu3W9BuXNoZiJ8xbKpkgHjl+teSNkbVFmTS+zHgk8qY21dTrvX/uyXL7uUF67/650OG9ril9Y7UB5alcHvtc7kbd5mfd0is/xu7PiFzt9l9erj/o3OPZp33X9p9q87uQ9rt/ymlzbGPMuL/M2bWje0P7cpOF5nd6f5lXPGvpe6NTP6U7n9e2xG57Z8vrZFNeI+B7FDRo2p7i+Y+spiVV+4dGxvGH7c6WndLXh+HUsbxX78iBW3Byl6n+oezHvhcCjOrg/O5W3X8v3QlBPY7eX8ob258EMO5WtH/JWkfW8MusVddTVg3lD+/IVVb/HezRvaH++rKF5Q/vzpbTMCOunvF5/jKNmWdxudxxwJsW0zDdExDPLt8+lOK/yE5n5J/P6Ou/P/ZbX5NrGmPf4Mu/WhuYN7c/bGp7X6f1pXvWsoe+FTv2c7nRe3x47KG7/PJQZEf8ZEU+jmH2zM8XpTn8DzgDuyMx/towz+yFv2P6cB7wlInYrlx1KcQrphzPzj2uy/W7mjXFffhVYmpmLy/Gt0X7s4bz/Bu7NzLvWNKeH8/7Kiu+FzMzlDcsb2p9LMvNvVbI6nTeGrDsp9uXtVXJ6OG/o5/TtVb7HezjvToq/K3c0NG/oe+GvmflAv+X1g6jhOK7YWLGDl0fEccByivNIv0FxHuJVFLcbflRm/iIioupfIvPM68Us88wzb3DymlzbIOS15AbFdRMmAA9RXGdtvfLtXYDzMvPHdWR1Mm8l+/OtwP9SXED2wsy8qGpWN/LKzEYeO/PM6/W8Jtdmnnm9ntframkyRcT4zHyo5fWzKe5M9DyKi1HfSHGnhrOz5bdy5pnX7rwm12aeeeZ1L6/JtQ1CXpnxryZVRLyU4k5Zr6D4UPgyiptN/IDit/51/Ga1Y3lj3J9PBS7JzMU1NFs7ndfYY2eeeb2c1+TazDOv1/P6SeXT5WINp32bZ16785pcm3nmmde9vCbXNgh5Zeb4lg+G44FrKK4TdBhwL8UpXTMoLixexwfRjuWNcX/28ylkjT125pnXy3lNrs0883o9r9/UNZOp0dPRzOvfvCbXZp555nUvr8m1DUJemTkOOJXibi/fobjb5/rAOhRNrb9njddO6GRe049fk4+deeb1cl6TazPPvF7P6ycT1vQLIx421flQVlzxfvj0sK/X1C00z7yeyzLPPPMGJ6/JtQ1CXpm5OfCncnufoPht41eAZ1Fc9PdaYDdgvazhwumdzGv68WvysTPPvF7Oa3Jt5pnX63n9ao1Ol2v6dDTz+jevybWZZ5553ctrcm2DkFfmHAfMB55WLloA/DwzrwfOp7jF8M3AKTV9EO1YXtOPX5OPnXnm9XJek2szz7xez+tna3y6XNOno5nXv3lNrs0888zrXl6Ta2t6XkR8AngE8EdgW+DTwCbAMynufrYWcDwwNzN/1295ZWYjj1/Tj5155vVqXpNrM8+8Xs/rd6s1kykiNo+IKF8OTQ87kWKK9ESK6WETKaaH1fHBwjzzei7LPPPMG5y8Jtc2CHll5quBx2bm6zPzE8B5wGuAWyhm4DwdeCdwQk0fRDuW1/Tj1+RjZ555vZzX5NrMM6/X8xohM8f0AI6juJDjnPL1ocDzy+dbAP8NTAHWHus2zTOvjrwm12aeeeZ1L6/JtQ1CXkvuUcCh5fNx5Z8HAV8Gdipfr9NveYNw/Jp67Mwzr9fzmlybeeb1el4THmOayVROD3sc8FPgyIh4KnAf8MyImANsVr7/6Mz8+1i2aZ55deQ1uTbzzDOve3lNrm0Q8oa5HtguIh6dmcsBMvPcMv/NEbE+UOfdX9qeN0DHr3HHzjzz+iSvybWZZ16v5/W9Vd5drmV62OHl6xdQTA/7NMVvrJ4O7EH909HMM69nsswzz7zByWtybYOQN4LrgF2Ap0fEgsxcFBE7ARuXmff0U96AHb9GHTvzzOujvCbXZp55vZ7X/7LHpoeZZ14vZplnnnmDk9fk2gYhb5QxbA68Gfg68A6KO8I8px/zBu34NenYmWdeP+U1uTbzzOv1vH5/rPLuchGxK7AvcFJm3tGy/DPABsDRwJJc1YbGyDzzejHLPPPMG5y8Jtc2CHkrGcdE4InAo8q8a/sxbxCPX1OOnXnm9Vtek2szz7xez+tnY7km03XAPRTTw6YBtEwP+0xm3lPzBwvzzOvFLPPMM29w8ppc2yDkjSgz/5mZN2TmJZ34YNjGvIE7fg06duaZ11d5Ta7NPPN6Pa+frXImE0BEbA4cCMyhuE3fi4F3ZeaP2jIo88zrwSzzzDNvcPKaXNsg5DWdx0+SJPWqMTWZoPnT0czr37wm12aeeeZ1L6/JtQ1CXtN5/CRJUi8ac5NJkiRJkiRJGs1YrskkSZIkSZIkrZRNJkmSJEmSJFVmk0mSJEmSJEmV2WSSJEmSJElSZTaZJEmSJEmSVJlNJkmSpA6KiIsiYla3xyFJklQ3m0ySJEmSJEmqzCaTJEkaaBGxbkR8LyKujYhfRcRLyuXvjYirymWnRUSUyy+KiBMj4ucRcVNE7BQR50bEbyPiQ+U6W0TEzRHxlYi4LiK+GRHrjJC9V0RcFhHXRMQ3ImLyCOtcFBEfi4grI+I3EfGMloxLyq+9JiKeXi7fPSIujohzyvU/GhGHlV9/fUQ8vlxvakR8q6zxqojYpX17WZIkDQKbTJIkadDtDfw5M7fPzOnA98vlp2TmTuWyRwD7tXzNPzJzN2A+8B3gGGA6cEREbFiu8yTgtMycASwBjm4NjYiNgHcDe2bmjsAC4K2jjHFCZs4G3gwcXy67A3hO+bUvAU5uWX974E3AdsDhwBPLr/8C8IZynU8BJ2bmTsALy/ckSZLW2IRuD0CSJKnLrgdOiIiPAf+TmZeUy58VEe8A1gEeBdwAfLd87/yWr70hM28HiIjfA5sCdwO3ZuYvyvW+CrwROKEldw6wDfCLcpLUWsBlo4zx3PLPq4EtyucTgVMiYibwEPDElvWvahnT74Aftoz3WeXzPYFtymyA9SJiSmbeO8oYJEmSVsomkyRJGmiZ+ZuIeCqwL/CRiPgh8HHgs8CszLw1It4HTGr5sr+Xfy5veT70eujzVQ6PGvY6gB9l5qFjGOZQxkMt238L8BeKWUvjgAdHWH/4GFvHNw7YOTMfGEO+JEnSKnm6nCRJGmgR8Tjg/sz8KsVMox1Z0VD6a3mdpIPXYNObRcTO5fNDgUuHvX85sEtEPKEcxzoR8UTGbn3g9sxcTnFK3PjVHN8PgdcPvShnREmSJK0xm0ySJGnQbQdcGRELgf8APpSZdwOfpzi97DzgqjXY7k3AKyLiOorT7U5tfTMzFwNHAGeV61wOPHk1tv/ZcvuXU5wqd99qju+NwKzywuQ3Akeu5tdLkiQ9TGQOn7ktSZKkKiJiC4rrO03v9lgkSZI6xZlMkiRJkiRJqsyZTJIkSZIkSarMmUySJEmSJEmqzCaTJEmSJEmSKrPJJEmSJEmSpMpsMkmSJEmSJKkym0ySJEmSJEmqzCaTJEmSJEmSKvv/Z9OUY7t+o1QAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1440x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define the size of the figure\n",
    "plt.figure(figsize=(20, 5))\n",
    "ax = plt.gca()\n",
    "\n",
    "df1.plot(x='sample', y=' total_bases',kind=\"bar\", color='blue', ax=ax, linewidth=3)\n",
    "plt.title('Total base sequenced by sample', fontsize=10)\n",
    "plt.xlabel('sample name')\n",
    "plt.ylabel('total base')\n",
    "plt.xticks(\n",
    "        rotation=45,\n",
    "        horizontalalignment='right',\n",
    "        fontweight='light',\n",
    "        fontsize='8',\n",
    "    )\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4UAAAF+CAYAAAA1Gr9TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA9eElEQVR4nO3deZwkdX3/8ddnD+UW0AUXWVgPPBFXGAEvxBUJEi+I4YiLGA/wSkw8wRhBo4lBUBM1CJ4EIoK3IWLEgyj+ZGFWASGQ8EskAXcD6y8gl6LC5/dH1bjNMEd3T3X3d6Zez8ejHtNdVf2Zd1Wfn66jIzORJEmSJLXTolEHkCRJkiSNjk2hJEmSJLWYTaEkSZIktZhNoSRJkiS1mE2hJEmSJLWYTaEkSZIktZhNoSSptSJi54j4SkRcGxH/GREfjoj7z6HehREx1mRGSZIGzaZQktRKERHAF4EvZ+ZuwG7A5sBJIw0mSdKQ2RRKktpqNfDLzPwUQGbeDfwp8JKI2CoiXhoRH56YOSLOi4j968unRsR4RFwVEe+c7R9FxHUR8ZcR8YP6dntGxD9HxH9ExKs65ntzRFwaEVd01o2IL0fEuvr/HdMx/vaIeE9EXB4RF0fEjg2sF0lSy9gUSpLa6nHAus4RmXkrcB3wiFlu+2eZOQbsATwjIvbo4v9dn5lPBr4HfBp4EbAv8C6AiDiQamvl3sAqYK+I2K++7csycy9gDPjjiHhgPX5L4OLMfALwXeCVXeSQJOlebAolSW0VQE4zfjaHRcQPgR9RNZeP7eI2X63//hhYm5m3ZeZG4JcRsS1wYD38CPgh8GiqJhGqRvBy4GJgRcf4XwHn1ZfXASu7yCFJ0r0sGXUASZJG5Crg9zpHRMQ2wI7AvwG7c+8vTzer53ko8CbgSZl5c0R8emLaLO6q/97TcXni+hKqZvSvMvO0SZn2Bw4AnpyZd0bEhR3/79eZOdHY3o3v65KkPrilUJLUVt8CtoiIlwBExGLgFODDmfkLqt1IV0XEoohYQbVbJ8A2wB3Az+tj+J7TUJ5/Bl4WEVvVeR4SETsADwBurhvCR1PtcipJUmP8RlGS1EqZmRFxCPCRiPhzYBlwTma+p57l+8BPqHb3vJJql04y8/KI+BHVlsb/rOdrIs83IuIxwA+qE6NyO7AG+Drwqoi4gmoL5sVN/D9JkibEpr1OJElqr4h4CnA2cGhmrpttfkmSFgqbQkmSJElqMY8plCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWa8VPUjzoQQ/KlStXjjqGJEmSJI3EunXrfpaZy6aa1oqmcOXKlYyPj486hiRJkiSNRET813TT3H1UkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQktSoDRvgkkuqv5IkqXw2hboPP9BpmHy89afE9XbnnXDUUbDTTrDPPtXfo46qxpeixPUmSdKo2RTqt+bDBzr1p8QPwj7e+lPyejv2WDjrrHuPO+usavyolbze5oMSX0OkUfC5oIXKpnBESnxRKfkD3XxQ4n1a8gfh+fB4K/E+LXW9bdhw31wTzjpr9Ouw1PXWqcTHW8mvIRNKXG8TzNafErP5XJgbs80Dmbngh7322itLcccdmWvWZMKmYc2aavworV9/70yTh/XrR5tvwvr1mWvXlpMns9z7NPO+uTrzjVLpj7dS79OS19vatTNnW7t2dNlKXm+Z5T7eMst9Dckse72ZbeFl87lgtoUAGM9p+qWRN2zDGEpqCkt9USn5A11m2U/cUu/Tkj8Il/54K/U+LXm9+XjrX6mPt5Lv08xy15vZFl42nwtmWyhG0hQCmwGXAJcDVwHvrMefA1xWD9cBl01z++uAH9fzjXeMPxH4aUeNg2fLUkpTWPKLSsnZMst94pa83kr+IFzyejNb/3yeLqxsvoaYzWwVnwtmWyhmagoHeUzhXcDqzHwCsAo4KCL2zczDM3NVZq4CvgB8cYYaz6znHZs0/gMTNTLzawNJPwDXXz+36YO0fDmsWTP1tDVrqumjUvKxSiXfpytWzG36IJX8eCv5Pi15vQGcdtp9861ZU40fpZLXW8mPt5JfQ0peb2brT8nZfC70x2zzy8Cawrohvb2+urQecmJ6RARwGHD2oDKUpuQXFSj3A13JT9yS79OSPwhDuY+3ku9TKHe9AWyxBZx5JqxfD2vXVn/PPLMaP2qlrreSH28lv4aUvN7M1p+Ss/lc6I/Z5pnpNiE2MQCLqXbxvB3460nT9mOm/VrhJ8APgXXAMR3jT6TatfQK4JPAdtPc/hhgHBjfZZddmtnm2oBSd6/qVNrJXErfxF/yfVrysZgTSnu8ZZZ9n04ocb3NByWut5IfbyW/hpS83sy28LL5XDDbQjBj7zXdhCYHYFvgO8DuHeNOBd44w212qv/uQHVc4n719R3rZnMR8B7gk7P9/1KOKcws+0WlZCU/cefDfVriB+GSzYf7VAvHfHi8lfgaUvJ6M9vCyzbB54LZ5rOZmsKopg9eRJwA3JGZJ0fEEqqTxeyVmTd0cdsTgdsz8+RJ41cC52Xm7jPdfmxsLMfHx/vOPggbNlS7Pa5YMfrd+OaDO++87++MTez6VcKuaeB9uhB5n2qYfLz1p+T1Zrb+lJytZCWvN7OVISLW5X3P1VJNG1RTGBHLgF9n5i0RsTnwDapdSM+LiIOA4zPzGdPcdktgUWbeVl++AHhXZn49IpZn5oZ6vj8F9snMI2bKUmJTqP606YkrSZIkNWWmpnDJAP/vcuCMiJjY1fPczDyvnnYEk04wExE7AR/PzIOpdhH9UnUuGpYAn8nMr9eznhQRq6hOWnMdcOwAl0GFWb7cZlCSJElq0tB2Hx0ltxRKkiRJarOZthQO8ncKJUmSJEmFsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFrMplCRJkqQWsymUJEmSpBazKZQkSZKkFlsyqMIRsRnwXeD+9f/5fGaeEBHnAI+qZ9sWuCUzV01x++uA24C7gd9k5lg9fnvgHGAlcB1wWGbePKjlkCRJkqSFbGBNIXAXsDozb4+IpcBFEXF+Zh4+MUNEnAL8fIYaz8zMn00adxzwrcx8b0QcV19/a9PhJUmSJKkNBrb7aFZur68urYecmB4RARwGnN1j6RcAZ9SXzwBeOLekkiRJktReAz2mMCIWR8RlwE3ABZm5tmPy04EbM/PaaW6ewDciYl1EHNMxfsfM3ABQ/91hmv99TESMR8T4xo0b57wskiRJkrQQDbQpzMy76+MFdwb2jojdOyYfycxbCZ+amXsCzwFeGxH79fi/T8/MscwcW7ZsWa/RJUmSJKkVhnL20cy8BbgQOAggIpYAh1KdMGa626yv/94EfAnYu550Y0Qsr+ssp9oKKUmSJEnqw8CawohYFhHb1pc3Bw4ArqknHwBck5k3THPbLSNi64nLwIHAlfXkrwJH15ePBr4ykAWQJEmSpBYY5NlHlwNnRMRiqubz3Mw8r552BJN2HY2InYCPZ+bBwI7Al6pz0bAE+Exmfr2e9b3AuRHxcuC/gd8f4DJIkiRJ0oIWmTn7XPPc2NhYjo+PjzqGJEmSJI1ERKyb+O33yYZyTKEkSZIkqUw2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiA2sKI2KziLgkIi6PiKsi4p31+HMi4rJ6uC4iLpuhxuKI+FFEnNcx7sSI+GlHjYMHtQySJEmStNAtGWDtu4DVmXl7RCwFLoqI8zPz8IkZIuIU4Ocz1Hg9cDWwzaTxH8jMkxtPLEmSJEktM7AthVm5vb66tB5yYnpEBHAYcPZUt4+InYHfBT4+qIySJEmS1HYDPaaw3v3zMuAm4ILMXNsx+enAjZl57TQ3/yDwFuCeKaa9LiKuiIhPRsR20/zvYyJiPCLGN27c2P9CSJIkSdICNtCmMDPvzsxVwM7A3hGxe8fkI5l+K+FzgZsyc90Uk08FHg6sAjYAp0zzv0/PzLHMHFu2bFn/CyFJkiRJC9hQzj6ambcAFwIHAUTEEuBQ4JxpbvJU4PkRcR3wWWB1RJxV17qxbjbvAT4G7D3Q8JIkSZK0gA3y7KPLImLb+vLmwAHANfXkA4BrMvOGqW6bmcdn5s6ZuRI4Avh2Zq6pay3vmPUQ4MrBLIEkSZIkLXyDPPvocuCMiFhM1Xyem5kTPy1xBJN2HY2InYCPZ+ZsPzFxUkSsojppzXXAsU2GliRJkqQ2icycfa55bmxsLMfHx0cdQ5IkSZJGIiLWZebYVNOGckyhJEmSJKlMNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYjaFkiRJktRiNoWSJEmS1GI2hZIkSZLUYktmmyEilgGvBFZ2zp+ZLxtcLEmSJEnSMMzaFAJfAb4HfBO4e7BxJEmSJEnD1E1TuEVmvnXgSSRJkiRJQ9fNMYXnRcTBA08iSZIkSRq6bprC11M1hr+IiFsj4raIuHXQwSRJkiRJgzfr7qOZufUwgkiSJEmShq+bYwqJiO2A3YDNJsZl5ncHFUqSJEmSNBzd/CTFK6h2Id0ZuAzYF/gBsHqgySRJkiRJA9ftMYVPAv4rM58JPBHYONBUkiRJkqSh6KYp/GVm/hIgIu6fmdcAjxpsLEmSJEnSMHRzTOENEbEt8GXggoi4GVg/yFCSJEmSpOHo5uyjh9QXT4yI7wAPAL4+0FSSJEmSpKHo9uyjTwN2y8xPRcQy4CHATwaaTJIkSZI0cLMeUxgRJwBvBY6vRy0FzhpkKEmSJEnScHRzoplDgOcDdwBk5nrAH7SXJEmSpAWgm6bwV5mZQAJExJaDjSRJkiRJGpZumsJzI+I0YNuIeCXwTeBjg40lSZIkSRqGbs4+enJEPBu4ler3Cd+RmRcMPJkkSZIkaeC6OvtoZl4QEWsn5o+I7TPzfweaTJIkSZI0cLM2hRFxLPAu4BfAPUBQHV/4sMFGkyRJkiQNWjdbCt8EPC4zfzboMJIkSZKk4ermRDP/AdzZa+GI2CwiLomIyyPiqoh4Zz3+nIi4rB6ui4jLZqixOCJ+FBHndYzbPiIuiIhr67/b9ZpNkiRJklTpZkvh8cD/qY8pvGtiZGb+8Sy3uwtYnZm3R8RS4KKIOD8zD5+YISJOAX4+Q43XA1cD23SMOw74Vma+NyKOq6+/tYvlkCRJkiRN0s2WwtOAbwMXA+s6hhll5fb66tJ6yInpERHAYcDZU90+InYGfhf4+KRJLwDOqC+fAbywi2WQJEmSJE2hmy2Fv8nMN/RTPCIWUzWQjwA+kplrOyY/HbgxM6+d5uYfBN4CbD1p/I6ZuQEgMzdExA7T/O9jgGMAdtlll37iS5IkSdKC182Wwu9ExDERsbw+nm/7iNi+m+KZeXdmrgJ2BvaOiN07Jh/J9FsJnwvclJmzbpGc4X+fnpljmTm2bNmyfstIkiRJ0oLWzZbCP6j/Ht8xrqefpMjMWyLiQuAg4MqIWAIcCuw1zU2eCjw/Ig4GNgO2iYizMnMNcGNELK+3Ei4Hbuo2hyRJkiTp3mbdUpiZD51imLUhjIhlEbFtfXlz4ADgmnryAcA1mXnDNP/z+MzcOTNXAkcA364bQoCvAkfXl48GvjJbFkmSJEnS1LrZfbRfy6l2Pb0CuBS4IDMnflriCCbtOhoRO0XE17qo+17g2RFxLfDs+rokSZIkqQ+RmbPPNc+NjY3l+Pj4qGNIkiRJ0khExLrMHJtq2iC3FEqSJEmSCjftiWYiYs+ZbpiZP2w+jiRJkiRpmGY6++gpM0xLYHXDWSRJkiRJQzZtU5iZzxxmEEmSJEnS8HXzO4XUPzr/WKrfDAQgM/9+UKEkSZIkScMxa1MYEScA+1M1hV8DngNcBNgUSpIkSdI8183ZR18EPAv4n8z8Q+AJwP0HmkqSJEmSNBTdNIW/yMx7gN9ExDbATcDDBhtLkiRJkjQM3RxTOB4R2wIfA9YBtwOXDDKUJEmSJGk4Zm0KM/M19cWPRsTXgW0y84rBxpIkSZIkDUM3J5rZb6pxmfndwUSSJEmSJA1LN7uPvrnj8mbA3lS7kfrj9ZIkSZI0z3Wz++jzOq9HxArgpIElkiRJkiQNTTdnH53sBmD3poNIkiRJkoavm2MKPwRkfXURsAq4fICZJEmSJElD0tVPUnRc/g1wdmZ+f0B5JEmSJElD1E1T+DngEfXlf8vMuwaYR5IkSZI0RNMeUxgRSyPig8D1wKeAM4D/jIjj6ulPHEpCSZIkSdLAzLSl8BRgC2BlZt4GEBHbACdHxKnAQcBDBx9RkiRJkjQoMzWFBwO7ZebESWbIzFsj4tXAz4DnDDqcJEmSJGmwZvpJins6G8IJmXk3sDEzLx5cLEmSJEnSMMzUFP5rRLxk8siIWANcPbhIkiRJkqRhmWn30dcCX4yIlwHrqH6r8EnA5sAhQ8gmSZIkSRqwaZvCzPwpsE9ErAYeBwRwfmZ+a1jhJEmSJEmDNevvFGbmt4FvDyGLJEmSJGnIZjqmUJIkSZK0wNkUSpIkSVKL2RRKkiRJUovZFEqSJElSi9kUSpIkSVKLDawpjIjNIuKSiLg8Iq6KiHfW48+JiMvq4bqIuKzb29bTToyIn3bUOHhQyyBJkiRJC92sP0kxB3cBqzPz9ohYClwUEedn5uETM0TEKcDPe7jtxfX0D2TmyQPMLkmSJEmtMLCmMDMTuL2+urQecmJ6RARwGLC619tKkiRJkpox0GMKI2JxvXvoTcAFmbm2Y/LTgRsz89o+bvu6iLgiIj4ZEdtNc/tjImI8IsY3btzYyPJIkiRJ0kIz0KYwM+/OzFXAzsDeEbF7x+QjgbP7uO2pwMOBVcAG4JRpbn96Zo5l5tiyZcvmuiiSJEmStCAN5eyjmXkLcCFwEEBELAEOBc7p9baZeWPdMN4DfAzYexCZJUmSJKkNBnn20WURsW19eXPgAOCaevIBwDWZeUOvt42I5R2zHgJcOYj8kiRJktQGgzz76HLgjIhYTNV8npuZ59XTjmDSrqMRsRPw8cw8eJbbnhQRq6hOPHMdcOwAl0GSJEmSFrSoTvS5sI2NjeX4+PioY0iSJEnSSETEuswcm2raUI4plCRJkiSVyaZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklrMplCSJEmSWsymUJIkSZJazKZQkiRJklpsyaAKR8RmwHeB+9f/5/OZeUJEnAM8qp5tW+CWzFzVzW3radsD5wArgeuAwzLz5kEthyRJkiQtZIPcUngXsDoznwCsAg6KiH0z8/DMXFU3gl8AvtjtbetpxwHfyszdgG/V1yVJkiRJfRhYU5iV2+urS+shJ6ZHRACHAWf3eNsXAGfUl88AXth4eEmSJElqiYEeUxgRiyPiMuAm4ILMXNsx+enAjZl5bY+33TEzNwDUf3cYVH5JkiRJWugG2hRm5t31bqI7A3tHxO4dk49kiq2EXd52VhFxTESMR8T4xo0bew8vSZIkSS0wlLOPZuYtwIXAQQARsQQ4lOqEMT3dFrgxIpbXdZZTbUmc6nanZ+ZYZo4tW7ZsbgsgSZIkSQvUwJrCiFgWEdvWlzcHDgCuqScfAFyTmTf0cduvAkfXl48GvjKI/JIkSZLUBoPcUrgc+E5EXAFcSnVc4Hn1tCOYtOtoROwUEV/r4rbvBZ4dEdcCz66vS5IkSZL6EJk5+1zz3NjYWI6Pj486hiRJkiSNRESsy8yxqaYN5ZhCSZIkSVKZbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxQbWFEbEZhFxSURcHhFXRcQ76/HnRMRl9XBdRFw2xW1XRMR3IuLq+rav75h2YkT8tKPGwYNaBkmSJEla6JYMsPZdwOrMvD0ilgIXRcT5mXn4xAwRcQrw8ylu+xvgjZn5w4jYGlgXERdk5r/W0z+QmScPMLskSZIktcLAmsLMTOD2+urSesiJ6RERwGHA6iluuwHYUF++LSKuBh4C/OvkeVW+DRvg+uthxQpYvnzUaSRJkiR1GugxhRGxuN499Cbggsxc2zH56cCNmXntLDVWAk8EOm/7uoi4IiI+GRHbNRxbDbnzTjjqKNhpJ9hnn+rvUUdV4yVJg7VhA1xySfVXknrla0i7DLQpzMy7M3MVsDOwd0Ts3jH5SODsmW4fEVsBXwD+JDNvrUefCjwcWEW1NfGUaW57TESMR8T4xo0b57Qc6s+xx8JZZ9173FlnVeNLUfILXsnZJJXLL+Q0Cr5nLRy+hrTTUM4+mpm3ABcCBwFExBLgUOCc6W5TH4f4BeAfMvOLHbVurJvNe4CPAXtP8z9Pz8yxzBxbtmxZU4tSrNJejDdsuG9DOOGss0afs+QXvJKzzQelPRc0dyXfpyVmmw9fyGnh8D1r4ZkPryElvvZOKDnbTAZ59tFlEbFtfXlz4ADgmnryAcA1mXnDNLcN4BPA1Zn5/knTOo9KOwS4suHo80qpL8bXXz+36YNW8gteydlKVupzoVPJbxQlZiv5Pi01W+lfyE0o8fGm/syH96ySH2+lZSv9NaTU197Ss3UlMwcyAHsAPwKuoGrc3tEx7dPAqybNvxPwtfry06hOSnMFcFk9HFxPOxP4cT3tq8Dy2bLstddeuVCtWZMJ9x3WrBltrvXrp841Maxfb7b5lq3T+vWZa9eWkyez3OdCZuYdd9w335o11fhRKzlbyfdpqdnWrp35NWTt2tHmK/nxpt6V/p5V8uOt1Gylv4aU+tpberYJwHhO17tNN2EhDQu1KSz9xbjUJ0fJL3glZ8ss903M58LCy1byfWq2/pX6eOtU4pdeE0rLVvp7VsmPt1KzlfwaYra5m6kpHMoxhRqM0nfRPO00WLPm3uPWrKnGj9KKFXObPkglZ4NydxMq+blQ8q44JWcr+T4tOdvy5fd93Z2wZs1ofxao5McblL3rV6nZSn7PKvnxVnK2kl9DSn7tLTlbt2wK57GSX4wBttgCzjwT1q+HtWurv2eeWY0fpZJf8ErOVvKbWMnPhZLfKErOVvJ9WnI2KPcLuZIfb1Dul15QbraS37NKfryVnA3KfQ0p+bW35Gzdsimcx0p+Me60fDnsvXc5eaDcFzwoN1vJb2IlPxdKfqMoOVvJ92nJ2aDcL+RKfryV/KVXydmg3Peskh9vJWeDcl9DSn7tLTlb16bbr3QhDQv1mMLMco/xmi9KOz6jU2nZSt9fvuTnQqnHjpSereT7tORsJSv18VbysXElZ+tU2ntWZrmPt9Kzlazk196Ss01ghmMKo5q+sI2NjeX4+PioYwzUhg3VlpoVK+bJtxGal446aupvrNesqb5FLEGJz4U777zv7l8T36SP+pvXkrNNKPE+nVBythKV+njbsKE6Tm8669eP7v4tOVvpSn28QdnZ5oOSX3tLzhYR6zJzbMppNoWSuuWb2NyU/EZRcjYtPCU+3kr+0qvkbPNBiY+3CSVn08JjU2hTKDXKNzFJC03JX3qVnE3S/GFTaFMoSZK6UPKXXiVnk1S+mZrCJcMOI0mSVKrly8ttuErOJml+8ycpJEmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFbAolSZIkqcVsCiVJkiSpxWwKJUmSJKnFWvHj9RGxEfivUeeYwoOAnxVaz2yjr9V0PbONvlbT9cw2+lpN1zPb6Gs1Xc9so6/VdD2zjb5W0/VKztakXTNz2VQTWtEUlioixjNzrMR6Zht9rabrmW30tZquZ7bR12q6ntlGX6vpemYbfa2m65lt9LWarldytmFx91FJkiRJajGbQkmSJElqMZvC0Tq94HpmG32tpuuZbfS1mq5nttHXarqe2UZfq+l6Zht9rabrmW30tZquV3K2ofCYQkmSJElqMbcUSpIkSVKL2RRKkiRJUovZFA5BREQDNRZFxNKJy3NP9du6Zuuvrtn6q2u2/uqarb+6xWZrWhPL2mStNq63ppmtP2brj9n6U3K2XnlM4YBExB8B/5aZ36ivR/a5siPiVOBO4IHAmzNzo9nMZjazmW3eZds+M/93LjUm1dsf2ABcl5l3zbGW663/eo3lM1vftfbHbP3U2h+z9VNrfxrMVopiv72bzyLibcDLgKdFxPMAMjP7+TYhIk4AtgCOBy4BjjCb2cxmNrPNu2wfB46PiCfMpU5HvfcDbwIOA06IiMVzqOV6KyCf2cxmtvZlK4lN4WD8C/A84J+AfSPi+VC9yULPm5qvpvqm9Vf15RWdE6P33XLMZjazmc1sQ8wWEauBh9R1nhYRe3RM66fxOgxYkZnPBU4F7gbu6bVOB9dbH5rMZzazma192YqTmQ4NDcCSSde3Aw4A3gM8DVgObNlPrXrcA4CP1pcfBWxtNrOZzWxmKzdbR437U30weRXwemBVrzUm1VvecflLwD715XC9DXa9DSqf2cxmtvZlK2lwS2FDImJxZv4mKn8WEVtl5s3AOuBzwGuArwIP7bHW2yJiq3rS5sCvIuK5wN92U8tsZjOb2cw2smyd77G/ysyfAt8EfgnsERGHRcRLIrr7tnpSvf/puHwNsDEiHgq8rJt6rjegj/XWdD6zmc1s7ctWrFF3pQtpoNod983ABcBngK3q8WPAj4DnzKHW1vX4c4FLgWeZzWxmM5vZyswGLK3/BvB24CRgUT1ue+B9wH8CT++z3l931DsK+CjwLWB/19vg1luT+cxmNrO1L1vJg1sK5ygiXtHxzcAbgG2ojs24Cvh0RGwB3AIck5nnz/Qtwiy1PhkR9wOuBN6dmd8ym9nMZjazlZcNIDN/Xdf8e6otZlsB76q/cb4Z2BV4dWZ+r896WwN/UddbBBwDvCMzL5ypjuutv/U2iHxmM5vZ2petaKPuSufzAJwMnNlx/bXABzqunwV8Etiivn6fYy56qPUZqm8jNputltnMZjazmW1k2f4SeHF9+enA5+rLE99Ov72+vm39d8bjUGapdzLwpvr6qpnquN76X29N5zOb2czWvmzzYXBLYZ8i4hhgNfDy+vruwGnA4oh4V0Q8EdhI9Tsm/xgRW2bmb+ZQ60bg/wH/NFMts5nNbGYz28iyfQh4SGb+Qz3qUuCOiPi9rH4f60pg94j4c+A22HSGzz7r/RjYO6qflPhxfZspt+y53vpbb03nM5vZzNa+bPOFP17fh4jYDtgTeALV7zDtCexNtQk5qfY33go4OzMviIjXAV/OzBsGWctsZjOb2cw2smx7AX8FvDEzfxwRL6TaTekuqt8AvBV4InAC1Rk+f5OZ756qVtP1XG9l1DOb2czWvmzzShawuXI+DcBBwEeofnvp2VQHr14NPGrSfFsMs5bZzGY2s5ltZNleVP/9Q+B1VB8ovsmmk61sD7wU2Lu+vi3woGHUc72VUc9sZjNb+7LNt8EthX2IiDdSPRBOo/rdk9VUp/G+ODNvHVUts5nNbGYz20hqnQX8d2a+LSKOBNYAH8nMr/VSZ4D1XG8jrmc2s5mtfdnmG5vCPkTE3wGPBO6g+ibhEcBTgP8Azs3Me0ZRy2xmM5vZzDa8WhGxJOtj5SLiIuDyzHxtRLyE6hvlHwCXdpur6XpNLmuTtUpfb03WM5vZzNa+bPPVolEHmG8iYg2weWYeAFwMvJPqzfASYG2Pb66N1TKb2cxmNrMNt1ZWP9i+JCK+AHwJ2C0i3p+Zf091IoNnA1t2m6vpeuB6G3U9s5nNbO3LNm9lAfuwljwAiydd3xt4W8f1LwEfpD7t9rBqmc1sZjOb2UaWbVHH5c2Aj3ZcPwc4p76887Drud7KqGc2s5mtfdnm++CWwllk5t0RsSgiPhgRr6H60cqxiDguIn4fWAJ8PjN/OcxaZjOb2cxmtuFni4jFmXlPVHYH7gGWRcRr61neC+wREQ/Lac64Och6rrfR1zOb2czWvmwLgccUTiMiIuuVExHnAt8FFgN7UT1Q9gUeDPwwM78+rFpmM5vZzGa20WTrqLuI6sfZbwH+G7gIeDvVCVcOBF6fmRcNs57rrax6ZjOb2dqXbd7LAjZXljZQ7zJD9Sa4NfAK4GFUu8q8oB6347Brmc1sZjOb2UaWrXM3oyOBNwFbADcAD6U63uShwO7Drud6K6Oe2cxmtvZlW0iDu49OITPvri+eQHU62ocBZ1O9KV5G9ZslK4Zdy2xmM5vZzDb8bB27GS2KiFfWoxcBn6A6A+dS4BWZ+ZPMvHLY9Vxvo69nNrOZrX3ZFpxRd6UlDcBrgLH68geBM+rLWwKnA6+iekM8cJi1zGY2s5nNbCPL9vD6bwB/DbwYeBDwrbreQ4ELgN8ddj3XWxn1zGY2s7Uv20IcRh6gpAHYrv67Z/3g+ObEgwPYCtgc2GXYtcxmNrOZzWzDzwbsBhxRX34YcBNwcH394cCJwDuA54yonuttxPXMZjaztS/bQh1GHqCEgY7TbgM7AVcBv0P1o7zf6OVB0mQts5nNbGYz22iyTar9JuDxwH5UP2L8pHr8olHUc72VV89sZjNb+7IttGHkAUY9sOmg+gAeVl9eTrUJ+UDgUVRnXnvgMGuZzWxmM5vZRput4/pzgTOARwNPAn5I/WFi2PVcb2XUM5vZzNa+bAt98Ccp+O0pac8F/pfq7ENnA7cBXwT+AviXzLxz2LXMZjazmc1sw80WEYuyPhEBcDzV1rL1VLscvRx4H9UZOCMzfzDsek0ua5O1Sl9vTdYzm9nM1r5sbdDqpjCi+i2miDgJ+CnwdeCTVGdY+zjwP1THUKwdZi2zmc1sZjPbaLLVNRcBXwauoToz3Ubg88DOwLHAKzPzjmHXc72VU89sZjNb+7ItdK1sCqM6Je3dHde3AX5FtUn5H4EdqA48fctsD5Yma5nNbGYzm9lGlm1RZt5TX94d+J3MPCUifkC1BW2c6sPFTZl58zDrud7KqGc2s5mtfdnaZMmoAwzbxBti/e3BR4EfAf+bmedExG3AnVT7GZ/U7ZtrE7XMZjazmc1sI8s2sZtRUJ2i/OrMvDIiTqY6Vfk1wKuBc3r5gNNEPddbGfXMZjaztS9b27R1S2EAH6E6Je0/AW8Hzgd+AryF6g3xn4ddy2xmM5vZzDbcbBG/3Y0yqI6n+zXV7pTnA1sAzwOeAByXmd8bdr0ml7XJWqWvtybrmc1sZmtftlbKAs52M+yB6pvQnwBPq6/vCvxFfXmLUdUym9nMZjazjSzbXwLvry/vB/wt1Rk4nwc8c5T1XG9l1DOb2czWvmxtGlq3+yhAZl4aEW8D/jgifgH8BnhqVMdW3D6qWmYzm9nMZrbh14qIZcAdwKMjYkVmfjcingVcn5lX91JrEPVcb6OvZzazma192Vpn1F3pMAfq3WU7rv8+8K9Up6l98rBrMcOPZfZar8laXfyvxuot1GyTHx9zqddkrdLXW4nZgBVN1Wuy1myPk1G/vs2XbB11dgCOAT4GHAJcDOwzynqut7Lqmc1sZmtftjYNIw8w0IWDlwIHAA/uGLdo0jwHA/8APHZYtep5t5yqRp/ZGqtVz3tUvawPaCBbY7XqeV9NtQvADg1ka6xWPe/bgRd2XJ+qmes2W2O16nk/AOw9yzzdZmusVj3vR4BnTxo3+QNst9kaq1XP+2HgQw2tt8Zq1fO+DTgcWNVAtsZq1fO+GNifjh89n/z61EO2xmp1zL9kmvHLqY6r+xpw2FT/a5D1mAfvWU2ttwHUe0BT9ZqsVU/fpcFsjdWqpz+swWyN1aqnT/u47CNbY7Xq6c9pMFtjterpyxrM1litevqWk64vmkO2xmo51Ott1AEGtmDw58CZVD9O+RbgiI5pQccHReC5VGdi23zQtep5/gn4d2D7+vrkN+pesjVWq57nY8BZwCeA10wxvZdsjdXqqPdZ4FTgpZNr9ZGtkVr1PDsCX6E61fFzJ9fqMVtjtep5Pgm8t7689RzXW2O16nleBnywvrwcuP8csjVWq57ng8CXgO8CL5jjc6GxWvU8J1O9Jr0eeOOkaYt6zNZYrXqek6ie838NvAF43hzWW2O16nleCyytLy+eZp6HAH9IdVzKyulqNV2Pst+zTqB6rXzRDPP0st6arvcXwGeAw2d4fHR7PzRWq573lcA/A4+f67I2Waue9xTgnVOM72e9NVarnvc04NQpxi/qI1tjtep5/xa4h4735zncp43Vquf9S+CLVF8KbTbH+7SxWvW8bwc+BBwKLJ9jtsZqOWwaFrFwbQ18OjPfDHwfeHhEHAqQNYCIOBLYDPjTzPzFoGtFxIOBH1N90DktIrbP6vS597kvZqvXZK16nscBv8zMNcAbgd+JiP2nmXe2bI3Vqud5JHBbZh6Rma8GDoyI1fX4nu6HJmtNyMwbqU53/HngjRFxWGetXuo1WSsiHgLsAfxrRLwL+FhEvDki9u1jvTVWa5JLI+IU4G+A4yLi8DnUa6RWRDy1vtkhwGuoHr87dkyPHtZbY7XqeRZTnWr7/Zn5N8A+EXFgRDyhXtZ7esjWWK0OW1N9QDyO6uxvvx8RB/SzrE3Wioj3UjVJZ0XE0qx+mmFxx/RF9TL/FLge+CHVsSlTaroesA3NvWc1Visi3gM8luo16bURcdCk6b2ut6brvQvYjeqLg8Mi4vET0zIze6nXZK0Ov6A68+GbI+Jpk/7f4h7rNVYrIk6i2np2QkTsERF79rusTdaq6z0ZuCMzXx0RL6jfn3eoa9zTy7I2Wauu9+764lOpHiMrJ03vZb01Vque/0+AxwGvozrByt4T0+r7oZf11litut4rgT2pviDdC3hlRDymo14v662xWpokC+hMBzEATwfeRLXFYBHwTODNdHwbSvUNw4XAY4ZYK6h3NaL65vvL1JvnJ9U7ZLZ6Tdaq59uGjl0pqba8PL2+fP8eszVWa+I21N/EA8+hOhbmQ8CJwE693A9N1qrnu1/H4+Rw4AFU3+aeQsfueF1ma6xWx7yPBX5AtfvijlT72r+ZjjMN9rCsjdWq531kXetv6+vPoNrK8Yg+sjVZa7uOy8uBv5q4DR1bhLp5/DZZa+I2wB/Uj9cPUP3W3PFUWzj27TFbY7U65j2K6pvlB9X3yRnAqyfN0+2yNlKL6nV73/ryW4Fz2fRcm7w19PeBf2Hm3aYaq0e9xZ3qg9ebmcP7TJO16nm2Ag7rWLZnAe+bw3prut6WwEEd10+iajaPpnrt7OV+aKzWpLqbA6+gej36FNWH2AfW06KXek3VApYCR1JtQbkY+DuqvRj+sF4Pvay3xmp1zLcD1WvQ6VRb0j4O/Amwex/L2mStZcCLO64fR308LtVraS/rrbFaHfMdCLy2vvwOqj21jmbSluUul7WxWvV8TwGOqS8/EPg2VcO5mHtvse0mW2O1HCat21EHaHRh7r0rz0OpPtysAXaux30K2LO+fAjwPeCRg641Rb3OpugN9ZPtQKoP2NFjtjnVmqLe0o7xrwN+j+obotf3kW1Otaaot7i+3RPZdBzlB9j04ayXbHOqNbnexPJS7RK1FHg/cBtwZK/Z5lprmnorOi4vofoQu0ef2fquNcX98GCq49m+wKZG6WTq00b3eJ/OqdYU9TrfXF5N9VtHW3aM6yXbnGpNUe+BVLvGnEZ9PArVh4pD+sg2p1pT1NuD6ku0v6M6iclLqJ5bD+kj25xqda577v1a+S7gc/Xlzdj05drhVLvkP3oY9di0G+ULgW2BP6NqhPt5z2qs1qTl3Jrqeb6ovj++UI9fzqYv13pZb43W67j+OKpdtH8XOIeO3bS7uB8aqzXF830rqkbk2VTvf+uAd/eTba61Jt0P21I1c++sr6+k+tLlub1ma6LWpHpb1Y/fj1E1dVtT7XH04j6yzbnWVI+RetzL69tNfvz09HibS61Jy7oH1d4ofwf8B/Bkqs+uL+9jvc251qR6q6h2U38h1c/i/CPV3lD79JFtzrUcprm/Rh2gsQWB9wJ3Ur2AT3ygeGL95H8P8EdU3xg8vp62F9O/uTZWa4Z69+uYfhpwM/UJMvrI1letGepNfKB5HvA54FvA7/SZra9a09WbNH0V1TdEe/aTrd9aM9WrX4g+QbXl8fHAXv1m66fWDPdDZ93d62V9TD/Z+q01Q71HU23FeCfVlwcXsak57zVbX7W6XNZ3A3/QcX0u663rWpPqnTupzsuovmQ5gurb+f17yDbnWlPUm9gC9MB62JmqWf8Cmxq5brPNqVY9faLJCKrGt7MR/1Oqrd7fAVZQffn3PWb+gNNYPar3knOodoO6kGoXshVUTXCv71mN1ZpiOXfi3g3wiVRN079Q/cZhr+ut6XoPYdNzbKv67xFUu4DeD3j4LPdDY7WmqLeivvwoqr0/jquX87NUu4Hu1kO2OdWaot6OdOzFUI9/CdUxqIt7XG9zqjVFvW2pTpL0CaoTTS2lOinSp6i+dHlED9nmVGuax0jnnjHvBf6o4/psz/vGak1RbzuqrbOPpdolHKpjhv+RqjnuZb3NqdYU9bao1/uxVI3mU4DVwPPreR7WQ7Y51XKYfhh5gEYW4r678nye+k2H6kPE46jeHCd2X5zpbEaN1ZqmXueHnYm6H6b+Rq3HbH3Xmq1ePe6xwLXAgXPJ1mutLpZ1MbAL1QHQB89xvfVUq4t6Y1RbbDuXvd9sPdWapt7nOuotoXqxPJ/6m+8es/Vdq8vn1mOodut5SgPP065rdXE/THxIPISOXVGHUauL++EJVB90Pk5/z9O+a81wP3Q+XrcAPs2mrY73OYvuIGpNUfczVLuOHQc8ox6/K7CR+nlfj3vAMOpx390oD2DTbpQ70tt7VmO1ZlnO/epx/wZcRf3l3hzWW5P13gY8tR6/D3Ap9Zca3dRrstYU9d4K/DHVycM+RtWUbNZntjnV6qh3dke9p9Xjn0x1KMV+PWZrpFZHvc/W9V5FtYXvjVTH7v4L8KweszVSa5rH7/71+NVUhyls3mO2RmpNsazHUX0J9Il6edcBB/S53uZUa4p6r2fTa/zDgSvo+DmcHrPNqZbDNOt41AHmvADT78rz+cnjh1lrlnoTuxltXv+d2JoWTPNBp8la3daj+kDxsIl6w6jVZb371fWmPTX3IGp18xipr99nl8FB1+pyWZdSHZs1sTvZXO7Trmt1s6zMj+fprE1Ik7W6fYw0td4arPfb3SjrvxNbXObyeOu61hQ130b14WYXYC2wL9WXQY9h04f/Xp5fc6rXUWfybpRfqsc/eLYsg6jVxXI+lep95a/Z9CF2LvdD0/X2pTqW/Y1s2mW829elOdeapd5u1HuhNPB467lWF8u6LdWW7mc1sN56qjVNvZV1vd3retsBj+ozW9+1ZlvWevyWE/WGWauLx8gTqPZw22eO92nPtaapt2td7ylUnx/2Z9MXar1m67uWwyz326gDzCl897vydNM8NFZrAWXbaZi1FkC2i+t60/5A+SBq9bisDxlmrR7rNbHeuq61ALJNPEaaeA3putZ8WW8d11dTHRvzT1QnltqV6vjwzl1nZ2pYG6s3aVln3I2yx/U2p1o9LOehbDqZ2Wwf5EZR7/eof55pWLW6qPfcjsdIMPuXhY3V6mFZX8Sk3T8HXauLer9b1zt62LV6WNajZ3vcNl2rh2Vd00C2nmr1sKx/wAw/gTOIWg7dDYuYx7I6Bfgiqh/ffSvwRxHxjHryF6n2eX5fZv73MGstoGzrh1lrAWR7eF3v+mHW6nFZfzrMWj3Wa2K9dV1rAWSbeIw08RrSda2ml7XJWhGxaKJeRLwnIvah2oL2ZKpdFG8G/h64KTN/3ZEhh1Fv0rIeB/xpROxXjzuSalfjv8zM/5ptWZus1eVyngXcnpkbZ1rGEdY7k+qnhv63gWXtqlaX9X7GpsdIZuY9w6jV47Lempk3D6tWl/X+X72sG4ZZq4dl/Xtgw0yP26Zr9bCsZ1I/RoZVq4dlPRP4Wc7yU1VN1lL3oov7uUj1A+aeiHgb1Q9/fobq+JjXUx0D8Eiqb/m+H7Hp96wGXctsZjOb2cw2mmx1zaA65mQJcDfVca/b1JOfCnw5M785W52m682wrG8A/g/VSSbOz8wLh1mr6eWcD/XMZjaztS+bupAFbK7sZaCPXXmGUctsZjOb2cw2smydv+n1B8BXqY4d2pXqBBOHMem35oZVr8tl7Wo3yiZrlb7eBnA/mM1sZmtZNofehnm1+2iXm5PvsyvPoGuZzWxmM5vZRpZtcdafICJiMfBDqmPpXkz1+57nUp18ZauJ+YZVr8tl7Wo3yiZrNb2cpdczm9nM1r5s6t2823205E3TZjOb2cxmtpFkW0T1w+1LqU7Xfz+qb5O3oGow78oejjtpsp7rrYx6ZjOb2dqXTb1ZMuoA3Yi417EkR7LpTE3bUO1idC3wz8Bnu/gmorFaZjOb2cxmtpFl2xX473re91F9k3wG8EyqE3JcTvVbeNtkFyfQabKe662MemYzm9nal01zkEPcV7WfgY5jKKh+C+rRVL8d9Dpge6rfUHk3PZxqv4laZjOb2cxmtpFlexvVVrKJ3/c6Enh+fXkl1VnptqbL36xssp7rrYx6ZjOb2dqXzWFuw7zYfbTkTdNmM5vZzGa24dWKiPcBmwP/RfVbfB+i+n3DZ1CdgfN+wAnAMZn5H8Ou1+SyNlmr9PXWZD2zmc1s7cumuVs06gDTiYhd6+MnYNPm5A9Q7TqzlGpz8lKqzckzviE2WctsZjOb2cw2smwvBx6cma/LzPcBXwZeAfxfqpMSPIXqtw5P7vIDTmP1XG9l1DOb2czWvmxqSC+bFYc1UPCmabOZzWxmM9todgsCXg0cWV9eVP89FPg08KT6+hbDrud6K6ee2cxmtvZlc2hmKG5LYb05eSfg28CrImIv4A7gGRGxL7BLPX2HzLxrWLXMZjazmc1so8nW4cfA4yNih8y8ByAzv1jX/pOIeADQy5np5lzP9VZcPbOZzWzty6YGFHX20Y7NyUfV119ItTn5Q1TfjD4FeBa9bZqecy2zmc1sZjPbaLJNcgXVTy88JSLGM/OGiHgSsLyu9/Nh1nO9FVnPbGYzW/uyqQlZwObKiYGCN02bzWxmM5vZhp9tivq7An8CfBZ4CzAOPHsU9VxvZdYzm9nM1r5sDnMfijr7aEQ8DTgY+GBm3tQx/iPAtsBrgFuzi9BN1jKb2cxmNrONJts0/2Mp8Eiqn2W4NTMv77fWXOq53sqtZ7bR1zKb2YadTXNT2jGFVwA/p9qcvDNAx+bkj2Tmz3t4Q2yyltnMZjazmW002e4jM3+dmVdl5vea+BAxh3qut0LrmW30tcxWRr02ZdPcFLWlECAidgUOAfalOi3tYcDxmXnBKGuZzWxmM5vZRpOtZK43SdJCUFxTCGVvmjbb6OuZbfS1zGa2YWcrmetNkjTfFdkUSpIkSZKGo7RjCiVJkiRJQ2RTKEmSJEktZlMoSZIkSS1mUyhJkiRJLWZTKEmSJEktZlMoSdIMIuLPIuKqiLgiIi6LiH0G+L8ujIixQdWXJGkqS0YdQJKkUkXEk4HnAntm5l0R8SDgfiOOJUlSo9xSKEnS9JYDP8vMuwAy82eZuT4i3hERl0bElRFxekQE/HZL3wci4rsRcXVEPCkivhgR10bEu+t5VkbENRFxRr318fMRscXkfxwRB0bEDyLihxHxuYjYaqhLLklqDZtCSZKm9w1gRUT8e0T8XUQ8ox7/4cx8UmbuDmxOtTVxwq8ycz/go8BXgNcCuwMvjYgH1vM8Cjg9M/cAbgVe0/lP6y2SbwcOyMw9gXHgDYNZRElS29kUSpI0jcy8HdgLOAbYCJwTES8FnhkRayPix8Bq4HEdN/tq/ffHwFWZuaHe0vifwIp62vWZ+f368lnA0yb9632BxwLfj4jLgKOBXZtcNkmSJnhMoSRJM8jMu4ELgQvrJvBYYA9gLDOvj4gTgc06bnJX/feejssT1yfed3Pyv5l0PYALMvPIOS+AJEmzcEuhJEnTiIhHRcRuHaNWAf9WX/5ZfZzfi/oovUt9EhuAI4GLJk2/GHhqRDyizrFFRDyyj/8jSdKs3FIoSdL0tgI+FBHbAr8B/i/VrqS3UO0eeh1waR91rwaOjojTgGuBUzsnZubGejfVsyPi/vXotwP/3sf/kiRpRpE5eY8VSZI0KBGxEjivPkmNJEkj5+6jkiRJktRibimUJEmSpBZzS6EkSZIktZhNoSRJkiS1mE2hJEmSJLWYTaEkSZIktZhNoSRJkiS1mE2hJEmSJLXY/we77d5zE8QAeQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define the size of the figure\n",
    "plt.figure(figsize=(15, 5))\n",
    "ax = plt.gca()\n",
    "\n",
    "#dot plot of the mean length by file\n",
    "df1.plot(x='sample', y='qual_mean',kind=\"scatter\", color='blue', ax=ax, linewidth=3)\n",
    "plt.title('Qual mean', fontsize=10)\n",
    "plt.xlabel('Sample')\n",
    "plt.ylabel('Qual mean')\n",
    "plt.xticks(\n",
    "        rotation=45,\n",
    "        horizontalalignment='right',\n",
    "        fontweight='light',\n",
    "        fontsize='8',\n",
    "    )\n",
    "plt.show()"
   ]
  },
  {
   "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.8.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
