ETE_methods.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import re\n",
    "import numpy as np\n",
    "import matplotlib.mlab as mlab\n",
    "import matplotlib.pyplot as plt\n",
    "import pickle\n",
    "import pandas as pd \n",
    "import csv\n",
    "from ete3 import Tree\n",
    "from ete3 import NCBITaxa \n",
    "from Bio import SeqIO"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "NCBI database not present yet (first time used?)\n",
      "Downloading taxdump.tar.gz from NCBI FTP site (via HTTP)...\n",
      "Done. Parsing...\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Loading node names...\n",
      "2110795 names loaded.\n",
      "208594 synonyms loaded.\n",
      "Loading nodes...\n",
      "2110795 nodes loaded.\n",
      "Linking nodes...\n",
      "Tree is loaded.\n",
      "Updating database: /home/adam/.etetoolkit/taxa.sqlite ...\n",
      " 2110000 generating entries... nerating entries... 166000 generating entries...  generating entries...  generating entries... "
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Inserting synonyms:       5000 "
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Uploading to /home/adam/.etetoolkit/taxa.sqlite\n",
      "\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Inserting taxid merges:  25000  "
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Inserting taxids:       10000  "
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Inserting taxids:       2110000 "
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n"
     ]
    }
   ],
   "source": [
    "## Install a local copy of the NCBI taxonomy database\n",
    "ncbi = NCBITaxa()\n",
    "\n",
    "## The first time you attempt to use NCBITaxa, ETE will detect that \n",
    "## your local database is empty and it will attempt to download the \n",
    "## latest NCBI taxonomy database (~300MB) and will store a parsed \n",
    "## version of it in your home directory: ~/.etetoolkit/taxa.sqlite. \n",
    "\n",
    "\n",
    "## To upgrade the database when previously installed:\n",
    "ncbi.update_taxonomy_database() \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "## parse a GenBank file, store the taxon identifiers in a dictionary\n",
    "# taxDict = {}\n",
    "\n",
    "# recs = [rec for rec in SeqIO.parse(\"viral.genomic.gbff\", \"genbank\")]\n",
    "## print(len(recs))\n",
    "## 10544\n",
    "\n",
    "## get taxids of RNA viruses\n",
    "for rec in recs:\n",
    "#     print(rec)\n",
    "    for feature in rec.features:\n",
    "        if 'source' in feature.type:\n",
    "            taxid = ''.join(feature.qualifiers['db_xref'])\n",
    "            taxid = re.sub(r'.*taxon:','',taxid)\n",
    "    taxonomy = rec.annotations['taxonomy']  \n",
    "    if (\"ssRNA viruses\" in taxonomy) or (\"dsRNA viruses\" in taxonomy):\n",
    "        taxDict[rec.id] = taxid\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "taxID_set = set(taxDict.values())\n",
    "\n",
    "## make a dictionary of all the taxon identifiers for each RNA viral family\n",
    "## (keys = families; values = list of taxon identifiers)\n",
    "rank2tax_dict = {}\n",
    "\n",
    "for taxid in taxID_set:\n",
    "    ## get lineage --> hierarchically sorted list of parent taxids.\n",
    "    try:\n",
    "        lineage = ncbi.get_lineage(taxid)   \n",
    "    except:\n",
    "        continue\n",
    "    names = ncbi.get_taxid_translator(lineage)\n",
    "    lineage2ranks = ncbi.get_rank(names)\n",
    "    ranks2lineage = dict((rank, taxid) for (taxid, rank) in lineage2ranks.items())\n",
    "\n",
    "    ranks = {'family_id': ranks2lineage.get('family')}\n",
    "    if ranks is not None:\n",
    "        for rank, id in ranks.items():\n",
    "            if rank.split('_id')[0] == 'family' and str(id).isdigit():\n",
    "                if id not in rank2tax_dict.keys():\n",
    "                    rank2tax_dict[id] = []\n",
    "                rank2tax_dict[id].append(taxid)\n",
    "\n",
    "# print(len(rank2tax_dict))\n",
    "\n",
    "## Dict of scientific names (instead of taxon identifier numbers)\n",
    "RNAfamDict = {}\n",
    "for family in rank2tax_dict.keys():\n",
    "    taxList = rank2tax_dict[family]\n",
    "    famname_scientific = (ncbi.get_taxid_translator([str(family)]))[family]\n",
    "    taxDict_scientific = ncbi.get_taxid_translator(taxList)\n",
    "    taxList_scientific = list(taxDict_scientific.values())\n",
    "    if family not in RNAfamDict:\n",
    "        RNAfamDict[famname_scientific] = taxList_scientific\n",
    "#         print(f\"{famname_scientific}\")\n",
    "\n",
    "\n",
    "## save a pickled version of the dict to file\n",
    "with open(\"RNA.family.pickle\", \"wb\") as handle:\n",
    "    pickle.dump(RNAfamDict, handle, protocol = pickle.HIGHEST_PROTOCOL)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/YAAAFuCAYAAAArsoNTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xe8bFV5+P/PA6iggvWqRKXYWwhRNPZONKIYuyRGgwVNrAkpaDREjT9b1KhRCfYWxUqwRCEo2KLSmx28JkQULFG/Fozw/P7Y+4Tx3Dkz+1lnZu6cez/v12teZ2aftfZefa89s0tkJpIkSZIkaWPaYWsnQJIkSZIktfPAXpIkSZKkDcwDe0mSJEmSNjAP7CVJkiRJ2sA8sJckSZIkaQPzwF6SJEmSpA3MA3tJkiRJkjYwD+wlSZIkSdrAPLCXJEmSJGkD22lrJ2A9rnnNa+Zee+21tZMhSZIkSdLMnXLKKd/LzE3Twm3oA/u99tqLk08+eWsnQ5IkSZKkmYuIbw0J56n4kiRJkiRtYB7YS5IkSZK0gXlgL0mSJEnSBuaBvSRJkiRJG5gH9pIkSZIkbWAe2EuSJEmStIF5YC9JkiRJ0gbmgb0kSZIkSRuYB/aSJEmSJG1gHthLkiRJkrSBze3APiLeFBEXRsTZI8uOiojT+9fmiDi9X75XRPx85H9HzCtdkiRJkiRtS3aa47rfAvwT8LaVBZn5iJX3EfEy4Ecj4c/NzH3nmB5JkiRJkrY5czuwz8xPRcRe4/4XEQE8HLjnvLYvSZIkSdL2YGtdY38X4LuZ+fWRZXtHxGkRcWJE3GUrpUuSJEmSpA1lnqfiT3IQ8K6RzxcAe2Tm9yPiNsDREXHLzPzx6ogRcQhwCMAee+yxkMRKkiRJkrSsFv6LfUTsBDwYOGplWWZenJnf79+fApwL3GRc/Mw8MjP3y8z9Nm3atIgkS5IkSZK0tLbGqfj3Br6SmeevLIiITRGxY//+BsCNgfO2QtokSZIkSdpQ5vm4u3cB/wHcNCLOj4jH9f96JL9+Gj7AXYEzI+IM4H3AkzLzB/NKmyRJkiRJ24p53hX/oDWW//GYZe8H3j+vtEiSJEmStK3aWnfFlyRJkiRJM+CBvSRJkiRJG5gH9pIkSZIkbWBb6zn2C7XXYR+ZGmbziw5YQEokSZIkSZotf7GXJEmSJGkD88BekiRJkqQNzAN7SZIkSZI2MA/sJUmSJEnawDywlyRJkiRpA/PAXpIkSZKkDcwDe0mSJEmSNjAP7CVJkiRJ2sA8sJckSZIkaQPzwF6SJEmSpA3MA3tJkiRJkjYwD+wlSZIkSdrAPLCXJEmSJGkD88BekiRJkqQNzAN7SZIkSZI2sJ22dgKW0V6HfWRQuM0vOmDOKZEkSZIkaTJ/sZckSZIkaQPzF/sZ8Bd+SZIkSdLWMrdf7CPiTRFxYUScPbLs7yLivyPi9P51v5H/PTMivhERX42I+8wrXZIkSZIkbUvmeSr+W4D7jln+iszct399FCAibgE8ErhlH+e1EbHjHNMmSZIkSdI2YW4H9pn5KeAHA4M/EHh3Zl6cmd8EvgHcbl5pkyRJkiRpW7E1bp73lIg4sz9V/2r9susC/zUS5vx+2RYi4pCIODkiTr7ooovmnVZJkiRJkpbaog/sXwfcENgXuAB4Wb88xoTNcSvIzCMzc7/M3G/Tpk3zSaUkSZIkSRvEQg/sM/O7mXlJZl4KvJ7LTrc/H7j+SNDrAd9eZNokSZIkSdqIFnpgHxG7j3x8ELByx/xjgEdGxBUiYm/gxsAXF5k2SZIkSZI2ork9xz4i3gXcHbhmRJwPHA7cPSL2pTvNfjPwRIDMPCci3gN8CfgV8OTMvGReaZMkSZIkaVsxtwP7zDxozOI3Tgj/AuAF80qPJEmSJEnboq1xV3xJkiRJkjQjHthLkiRJkrSBeWAvSZIkSdIG5oG9JEmSJEkbmAf2kiRJkiRtYB7YS5IkSZK0gQ163F1EHADcEth5ZVlmPm9eiZIkSZIkScNMPbCPiCOAKwL3AN4APBT44pzTtU3b67CPDAq3+UUHzDklkiRJkqSNbsip+HfMzEcDP8zM5wJ3AK4/32RJkiRJkqQhhhzY/7z/+7OI+A3gf4G955ckSZIkSZI01JBr7D8cEVcFXgqcCiTdKfmSJEmSJGkrm3pgn5nP79++PyI+DOycmT+ab7IkSZIkSdIQU0/Fj4jnR8ROAJl5MZAR8ea5p0ySJEmSJE015Br7nYAvRMQ+EfG7wEnAKfNNliRJkiRJGmLIqfjPjIjjgS8APwTumpnfmHvKJEmSJEnSVENOxb8r8ErgecAJwD/1d8eXJEmSJElb2ZC74v8D8LDM/BJARDwY+ARws3kmTJIkSZIkTTfkwP4OmXnJyofM/EBEnDjHNEmSJEmSpIGGXGN/SUQcANwS2HnkX8+bW6okSZIkSdIgQ66xPwJ4BPBUIICHAXvOOV2SJEmSJGmAIY+7u2NmPhr4YWY+F7gDcP35JkuSJEmSJA0x5MD+5/3fn/V3w/9fYO9pkSLiTRFxYUScPbLspRHxlYg4MyI+GBFX7ZfvFRE/j4jT+9cRLZmRJEmSJGl7M+TA/sP9AfhLgVOBzcC7B8R7C3DfVcuOA26VmfsAXwOeOfK/czNz3/71pAHrlyRJkiRpuzfkrvgvycyLgfdHxIfpbqD3i2mRMvNTEbHXqmXHjnz8PPDQ4UmVJEmSJEmrDfnF/j9W3mTmxZn5o9Fl6/BY4N9GPu8dEadFxIkRcZe1IkXEIRFxckScfNFFF80gGZIkSZIkbVxr/mIfEdcBrgvsEhG/TXdHfIDdgCuuZ6MR8TfAr4B39osuAPbIzO9HxG2AoyPilpn549VxM/NI4EiA/fbbL9eTDkmSJEmSNrpJp+LfB/hj4HrAy7jswP7HwLNaNxgRjwHuD9wrMxO6MwGAi/v3p0TEucBNgJNbtyNJkiRJ0vZgzQP7zHwr8NaIeEhmvn8WG4uI+wJ/DdwtM382snwT8IPMvCQibgDcGDhvFtuUJEmSJGlbNvXmea0H9RHxLuDuwDUj4nzgcLq74F8BOC4iAD7f3wH/rsDzIuJXwCXAkzLzBy3blSRJkiRpezLkrvhNMvOgMYvfuEbY9wMzOStAkiRJkqTtyZC74kuSJEmSpCU19Rf7iHjwpP9n5gdmlxxJkiRJklQx5FT8xwF3BD7Rf74HcALwIyABD+wlSZIkSdpKhhzYJ3CLzLwAICJ2B16TmQfPNWWSJEmSJGmqIdfY77VyUN/7Lt0z5iVJkiRJ0lY25Bf7EyLi48C76H69fyTwybmmSpIkSZIkDTLkOfZPiYgH0T1rHuDIzPzgfJMlSZIkSZKGGPoc+1OBn2Tmv0fEFSNi18z8yTwTJkmSJEmSppt6jX1EPAF4H/DP/aLrAkfPM1GSJEmSJGmYITfPezJwJ+DHAJn5deBa80yUJEmSJEkaZsiB/cWZ+cuVDxGxE91N9CRJkiRJ0lY25MD+xIh4FrBLROwPvBf40HyTJUmSJEmShhhyYH8YcBFwFvBE4KPAs+eZKEmSJEmSNMyQx91dCrw+It4K3BL478z0VHxJkiRJkpbAmr/YR8QREXHL/v1VgNOBtwGnRcRBC0qfJEmSJEmaYNKp+HfJzHP69wcDX8vM3wRuA/zV3FMmSZIkSZKmmnRg/8uR9/vTP7s+M78z1xRJkiRJkqTBJh3Y/09E3D8ibk33HPuPwf897m6XRSROkiRJkiRNNunmeU8EXgVcB3jGyC/19wI+Mu+ESZIkSZKk6SYd2D8uM+8bEQ/PzPesLMzMjwMfn3/SJEmSJEnSNJNOxb9fRFyO7jn2kiRJkiRpCU06sP8Y8D1gn4j4cUT8ZPTvkJVHxJsi4sKIOHtk2dUj4riI+Hr/92r98oiIV0XENyLizP7afkmSJEmSNMGaB/aZ+ZeZeRXgI5m5W2buOvp34PrfAtx31bLDgOMz88bA8Vx2RsDvATfuX4cAryvkQ5IkSZKk7dKkX+wByMwHRsSeEXFvgIjYJSJ2HbLyzPwU8INVix8IvLV//1bg90eWvy07nweuGhG7D9mOJEmSJEnbq0k3zwMgIp5A9wv61YEbAtcDjqC7O36La2fmBQCZeUFEXKtffl3gv0bCnd8vu2BVeg7p08Mee+zRmISNZ6/Dpj+IYPOLDlhASiRJkiRJy2TqL/bAk+meY/9jgMz8OnCtiTHaxJhlucWCzCMzc7/M3G/Tpk1zSIYkSZIkSRvHkAP7izPzlysfImInxhxwF3x35RT7/u+F/fLzgeuPhLse8O11bEeSJEmSpG3ekAP7EyPiWcAuEbE/8F7gQ+vY5jHAY/r3jwH+dWT5o/u7498e+NHKKfuSJEmSJGm8IQf2hwEXAWcBTwQ+Cjx7yMoj4l3AfwA3jYjzI+JxwIuA/SPi68D+/Wf69Z4HfAN4PfCnhXxIkiRJkrRdmnrzvMy8lO5A+/UryyLiTsBnB8Q9aI1/bXHjvcxMuuv5JUmSJEnSQGse2EfEjsDD6e5M/7HMPDsi7g88C9gF+O3FJFGSJEmSJK1l0i/2b6S7md0XgVdFxLeAOwCHZebRi0icJEmSJEmabNKB/X7APpl5aUTsDHwPuFFmfmcxSZMkSZIkSdNMunneL/vr68nMXwBf86BekiRJkqTlMukX+5tFxJn9+wBu2H8Ounvd7TP31EmSJEmSpIkmHdjffGGpkCRJkiRJTdY8sM/Mby0yIZIkSZIkqW7SNfaSJEmSJGnJeWAvSZIkSdIGtuaBfUQc3/998eKSI0mSJEmSKibdPG/3iLgbcGBEvJvubvj/JzNPnWvKtC57HfaRQeE2v+iAOadEkiRJkjRPkw7s/xY4DLge8PJV/0vgnvNKlCRJkiRJGmbSXfHfB7wvIp6Tmc9fYJokSZIkSdJAk36xByAznx8RBwJ37RedkJkfnm+yJEmSJEnSEFPvih8RLwSeDnypfz29XyZJkiRJkrayqb/YAwcA+2bmpQAR8VbgNOCZ80yYJEmSJEmabuhz7K868v4q80iIJEmSJEmqG/KL/QuB0yLik3SPvLsr/lovSZIkSdJSGHLzvHdFxAnAbekO7P86M78z74RJkiRJkqTphvxiT2ZeABwz57RIkiRJkqSiodfYS5IkSZKkJTToF/tZioibAkeNLLoB8Ld0N+h7AnBRv/xZmfnRBSdPkiRJkqQNZeIv9hGxQ0ScPcsNZuZXM3PfzNwXuA3wM+CD/b9fsfI/D+olSZIkSZpu4oF9/+z6MyJijzlt/17AuZn5rTmtX5IkSZKkbdqQU/F3B86JiC8CP11ZmJkHzmD7jwTeNfL5KRHxaOBk4NDM/OEMtiFJkiRJ0jZryIH9c+ex4Yi4PHAg8Mx+0euA5wPZ/30Z8Ngx8Q4BDgHYY495nUiw/dnrsI8MCrf5RQfMOSWSJEmSpIqpd8XPzBOBzcDl+vcnAafOYNu/B5yamd/tt/PdzLykP/3/9cDt1kjPkZm5X2but2nTphkkQ5IkSZKkjWvqgX1EPAF4H/DP/aLrAkfPYNsHMXIafkTsPvK/BwEzvWmfJEmSJEnboiGn4j+Z7tfzLwBk5tcj4lrr2WhEXBHYH3jiyOKXRMS+dKfib171Py0ZT92XJEmSpOUw5MD+4sz8ZUQAEBE70R18N8vMnwHXWLXsj9azTkmSJEmStkdTT8UHToyIZwG7RMT+wHuBD803WZIkSZIkaYghB/aHARcBZ9GdHv9R4NnzTJQkSZIkSRpm6qn4mXlpRLyV7hr7BL6ames6FV+SJEmSJM3G1AP7iDgAOAI4Fwhg74h4Ymb+27wTp22HN9uTJEmSpPkYcvO8lwH3yMxvAETEDYGPAB7YS5IkSZK0lQ25xv7ClYP63nnAhXNKjyRJkiRJKljzF/uIeHD/9pyI+CjwHrpr7B8GnLSAtEmSJEmSpCkmnYr/gJH33wXu1r+/CLja3FIkSZIkSZIGW/PAPjMPXmRCpNWG3HDPm+1JkiRJ2t4NuSv+3sBTgb1Gw2fmgfNLllTnnfclSZIkbY+G3BX/aOCNwIeAS+ebHEmSJEmSVDHkwP4XmfmquadEkiRJkiSVDTmwf2VEHA4cC1y8sjAzT51bqiRJkiRJ0iBDDux/E/gj4J5cdip+9p8lSZIkSdJWNOTA/kHADTLzl/NOjCRJkiRJqtlhQJgzgKvOOyGSJEmSJKluyC/21wa+EhEn8evX2Pu4O0mSJEmStrIhB/aHzz0V0lbgc+8lSZIkbQumHthn5omLSIgkSZIkSaqbemAfET+huws+wOWBywE/zczd5pkwSZIkSZI03ZBf7Hcd/RwRvw/cbm4pkpaUp+5LkiRJWkZDrrH/NZl5dEQcNo/ESNuSli8ChsTxiwNJkiRJo4aciv/gkY87APtx2an5zSJiM/AT4BLgV5m5X0RcHTgK2AvYDDw8M3+43m1J2yrPIpAkSZI05Bf7B4y8/xXdAfcDZ7T9e2Tm90Y+HwYcn5kv6s8KOAz46xltS5IkSZKkbc6Qa+wPXkRCeg8E7t6/fytwAh7YS5IkSZK0pjUP7CPibyfEy8x8/jq3ncCxEZHAP2fmkcC1M/OCfgMXRMS1xqTrEOAQgD322GOdSZAkSZIkaWOb9Iv9T8csuxLwOOAawHoP7O+Umd/uD96Pi4ivDInUfwFwJMB+++237mv9JUmSJEnayNY8sM/Ml628j4hdgacDBwPvBl62VryhMvPb/d8LI+KDdI/Q+25E7N7/Wr87cOF6tyPpMt5sT5IkSdr27DDpnxFx9Yj4e+BMui8Bbp2Zf52Z6zrgjogr9V8WEBFXAn4XOBs4BnhMH+wxwL+uZzuSJEmSJG3rJl1j/1LgwXSnvf9mZv6/GW732sAHI2IlDf+SmR+LiJOA90TE44D/BB42w21KkiRJkrTNmXSN/aHAxcCzgb/pD8IBgu7mebu1bjQzzwN+a8zy7wP3al2vJEmSJEnbm0nX2E88TV+SJEmSJG19U59jL2n75c32JEmSpOXngb2kmRryZYBfBEiSJEmz44G9pK2qelaAZxFIkiRJv87r6CVJkiRJ2sA8sJckSZIkaQPzwF6SJEmSpA3MA3tJkiRJkjYwb54naZvmzfYkSZK0rfMXe0mSJEmSNjAP7CVJkiRJ2sA8FV+SRrScuj8kjqf6S5IkaV48sJekBat+eeB9AiRJkjSJp+JLkiRJkrSBeWAvSZIkSdIG5oG9JEmSJEkbmAf2kiRJkiRtYB7YS5IkSZK0gXlXfEnaxngXfUmSpO2LB/aStJ3ziwBJkqSNzVPxJUmSJEnawBZ+YB8R14+IT0bElyPinIh4er/87yLivyPi9P51v0WnTZIkSZKkjWZrnIr/K+DQzDw1InYFTomI4/r/vSIz/2ErpEmSJEmSpA1p4Qf2mXkBcEH//icR8WXguotOhySp3ZDr8r0mX5IkaTG26jX2EbEX8NvAF/pFT4mIMyPiTRFxtTXiHBIRJ0fEyRdddNGCUipJkiRJ0nLaagf2EXFl4P3AMzLzx8DrgBsC+9L9ov+ycfEy88jM3C8z99u0adPC0itJkiRJ0jLaKgf2EXE5uoP6d2bmBwAy87uZeUlmXgq8Hrjd1kibJEmSJEkbyda4K34AbwS+nJkvH1m++0iwBwFnLzptkiRJkiRtNFvjrvh3Av4IOCsiTu+XPQs4KCL2BRLYDDxxK6RNkiRJkqQNZWvcFf8zQIz510cXnRZJkiRJkja6rXpXfEmSJEmStD5b41R8SdJ2Zshz7wE2v+iAOadEkiRp2+OBvSRp6fhFgCRJ0nCeii9JkiRJ0gbmgb0kSZIkSRuYp+JLkja8llP3h8TxVH9JkrQReGAvSdIAXvcvSZKWlQf2kiTNQfWLAL84kCRJrbzGXpIkSZKkDcwDe0mSJEmSNjAP7CVJkiRJ2sC8xl6SpA3IJwFIkqQVHthLkqSxvKGfJEkbgwf2kiRpJnwSgCRJW4fX2EuSJEmStIH5i70kSdoQFnFfAc8ikCRtRB7YS5IkNVrE5Qfz/nLCLzMkaePzwF6SJEmD+UWAJC0fD+wlSZI0N34RIEnz54G9JEmSloqXH0hSjXfFlyRJkiRpA1u6X+wj4r7AK4EdgTdk5ou2cpIkSZK0HdsWbpIoadu2VAf2EbEj8Bpgf+B84KSIOCYzv7R1UyZJkiRtXMv45YSk2VmqA3vgdsA3MvM8gIh4N/BAwAN7SZIkaRviWQfS7ERmbu00/J+IeChw38x8fP/5j4DfycynjIQ5BDik/3hT4KuNm7sm8L3tKPwitmEeZh9+EdtYtvCL2IZ5mH34RWxj2cIvYhvLFn4R2zAPsw+/iG0sW/hFbMM8zD78IraxbOEXsQ3zMPvwi9oGwJ6ZuWlqqMxcmhfwMLrr6lc+/xHw6jlt6+TtKfwypml7zIN5Xo5tmAfzbJ633zwvY5rMs3neKGkyz+Z5o6RpW8hz9bVsd8U/H7j+yOfrAd/eSmmRJEmSJGnpLduB/UnAjSNi74i4PPBI4JitnCZJkiRJkpbWUt08LzN/FRFPAT5O97i7N2XmOXPa3JHbWfhFbMM8zD78IraxbOEXsQ3zMPvwi9jGsoVfxDaWLfwitmEeZh9+EdtYtvCL2IZ5mH34RWxj2cIvYhvmYfbhF7WNwZbq5nmSJEmSJKlm2U7FlyRJkiRJBR7YS5IkSZK0gXlgL0mSJEnSBuaBvX5NRDwyIv6mf3/9iLjNlPDXi4h79O+vEBFXWkQ6tdxsFxuT9bYxWW9aS0TcPiIe3b+/RkTssbXTJEmaj+3qwD4irh0Rb4yIf+s/3yIiHjchfETEoyLib/vPe0TE7Wa1/j7MnhFx7/79LhGx65Twu0TETSeFaU1TRPwTcA/gUf2inwJHTFj3Y+keR/iGftGewL9OSU+pTKtxIuKKEfGciHh9//nGEXH/Cesuhe/DlOqgj1Ot55tExPERcXb/eZ+IePYs8lFddzVOS7vo4905Ig7u32+KiL2nxenDDjqIaSjTcjlV87GgMWNQelrrbST+1Hpo6J/l8unDldpStUyr25hnPa+jv1X3JS1lNDQPTfXchy19idE4fle3MTh8Q38eHL4frw4HVsatnYF/GZCmubWNhnG4Om9r2ae39M/BZdSYpnnvD1vyXGnXc9/fNmwjoj7/rIzz5Xruw1XHl+r+begcoKWdtpRptd7KY3Yfr9Jey/2tUK7N+7cmmbndvIB/Ax4OnNF/3gk4a0L41wGvAb7cf74acNIM1/8E4CTg3P7zjYHjJ4R/APBV4Jv9532BY2aVZ+DU/u9pI8vOmLDu04HLrwp/5pT0lMq0Ggc4Cvgr4Oz+8y7A6RPWXQ3fUgeleu7DnAjcblXZnj2LfFTXXY3T2C4OBz4EfK3//BvAZ6fEuSPwJeA/+8+/Bbx2hmXaUk6lfFT6Z0tbqqSnpd6q9dDQ30rl01gHLf1zaeq5sb+VxrHGMqrkoaWeS/2/Md/VMaYavtqfq+FPB2KZ2gb1cbjad0pjTOM2qmVUHfcWsT+szA1b+toi9rfVbVTn9NVxvlrPLeVaTVNlDtDSd1rm9JX5ZMucu9oXWvpbpVzL+7f1vOay0mV9rTS2VY1pUqerHuhW1z9uQjZpZ3IKcBVqO+nBaQK+QHcWx0q+rzEab0z4z4+uG9hxWmOtlmk1DnBysc6q4VvqoFTPjW1pcD6q625oRy3tomUC+gXg+gzfqVfLtKWcSvloSFN1zBicnpZ6q9ZDQ39bRB209M+lqefG/lYax9ZRRkPz0FLPpf7fmO/qGFMN39KfK+G/0P9d2YdecVJ+F9E2GvrC3PaF69hGtYyq495S7Q8b+9oi9rfzntNXx/lqPbeUa8v+begcoKXvtMzpK22vZc7dMg5X+1ulXMttez2v7epUfOCnEXENIAEi4vbAjyaE/9+I2HEk/Cbg0hmu/+LM/OXKh4jYaSXuGn6VmZPWt940vQZ4P7ApIp4LfAZ48YR1fzYi/grYObrrO48CPjwlPdUyrcb5ZUTsMhL2hsDFE9ZdDd9SB9V6Bvhen5aVdD0UuGBC+Eo+quuuxmlpF7/MbrRbWf+gU6gy879WLbpkQvBqvlvKqZqPeY8ZlfS01BtQqodqf6uWD9TroKV/LlM9t9RbdRxrKaNKnJZ6rvZ/aBi/q9sohq+WazX8ByLiNcBV+lNGjwXeNCE8zL9tVMfVatuojjEt26iWUXlesoD9YSnPDX1tEfvbapzq/LNaD+W211Cu5TQVwrf0nZY5faXeWubc1XJt6W+VOE37t2bz+sZgGV/ArYHP9gX6WeBrwD4Twv8h3bWL5wMvoDsd5GEzXP9LgGcBXwH2Bz4IvGBC+DcCfwCcSXe626uBI2ac51sCTweeAdxqyrp3BP6kT/fR/fsdpsQplWk1Tl+OJwIXAe8ENgN3n7DuaviWOijVcx/nBsC/Az8D/pvuS5a9ZpGP6rqrcRrbxV8A/wycR3dq538AT50S5310p1ydSveL0V8A755hmbaUUykfzH/MGJyelnqr1kNDfyuVT2MdtPTPpannxv5WGscay6iSh5Z6LvX/xnxXx5hq+Gp/bqmH3wNeAfwj8HsD+vNc2wb1cbjad0pjTOM2qmVUHfcWsT8cnOfGvraI/W11G9U5fXWcr9ZzS7lW01SZA7T0nZY5fWU+2TLnrvaFlv5WKdfy/m09r+g3ut3ov02+Kd0pFF/NzP+dEv5mwL368Mdn5pdntf6I2AF4HPC7ffiPA2/INSolIq4I/M2q8M/PzF/MKk19+KvT3WQHgMz89qTwVdUyrcbpvxm7fR/285n5vSnrHhy+pQ6q9bwq7pXoJuk/GRC2mu/B615PnMK692ekjDLzuCnhrwm8Erh3H+dY4OmZ+f0p8Up5aAhfzcfcxoyW9FRV66GhnZbGrz7O4Dy39s9lq+eK6jjW2O6q+7fqfqrc/xvyXW3b1fDVMppru+i3Mfe20cer7NuqbaM0xlS30TgPqI57c98fDs1z6762jzvX/W01TsOcvloPlflk6xymmqbK/rCl75Tn9H28qfXW2Nda+kJ5nlQs1/I8ptX2eGB/K+AW/PqB69umxLnWqvD/Ocv1z9vQNEXEAXTf7F8P+B5wXeDrmXmzNdZ7Q7pv6Fav+yYD0jRiFqWKAAAgAElEQVS4TKtxIuJqdN/sjYb91IT1lsIvSl8ft+TX0/W8CeEH56O67kqc9bSLeWso03I5NaRpKcaMRdVbQ/9civJZr3nlY5n7W8W2Us/LJCJuS/cL182BK9BNKi/OzN22crqq43CpbbTs0+fd/pZxnrGAPM99f9sYpzz/HGoZ67miNf3VMl3E3GqZLHL/ttM8VrqsIuJw4O50hftRulPUPgOMLdyIOBB4Gd3dDi+ke4zQl+ka4yzWf2PghWxZ2TdYI/wmujtWru4M9xwXviFNLwDuBBybmb/dfxv1kLXWDbwF+HvgH/r1HsyUa2uqZVqNExGPp7uU4Hp0N7e4Pd0pMmPLqCF8Sx2U6rmPcwTdjY7uQfcYq4cCX5wQfnA+qutuiPMW6u3i9lw2Ab083enFP500AY2Inel+KVpdF4+dQR5ay6mUjwWMGZX0vIVivfXbGFwPDf2tVD59nGodtPTPZarnt1Dvb6VxrLGMBsdprOdS/+/jVPNdHWOq4av9uVoPr6V7fO276e5C/cd0N5Va07zbRsM4XO07pTGmcRvVMqqOe4vYHw7Oc2NfW8T+trqN6py+Os5X67mlXKtpGhy+se+0zOkH11vjnLvaF1r6W6Vcy/u3dck5neO/jC/gLLq7vq88cuDawIcmhD+DkTvD0zXCI2e4/s/Qnb5yJl1n+DvguRPCH0vXWL8M3I3uJjgvnlWeueyOmGdw2dkcX5yw7lNWtjGy7NNT0lMq02qcPr87099xErgZcNSU8qmEb6mDUj33cc5c9ffKdF+4rDsf1XVX4zS2i5OBGwGn0Q2QBzP92tH3As8HzgUe09fNK2dYpi3lVMoH8x8zBqenpd6q9dDQ30rl01gHLf1zaeq5sb+VxrHGMqrkoaWeS/2/Md/VMaYavtqfq+HHtY3Pbc22Qdu+rdJ3SmNM4zaqZVQd9xaxP6zMDVv62iL2t9VtVOf0LeN8pZ5byrWapsocoKXvtMzpK/PJljl3tS+09LdquZb2b+t5zWWly/qiP0ile3zCbnSnpZ0zqeL6v2fQ34yIyQe61fWXJmQj4c8cWXbirPIMHA9cie7u+G+n+xbu8xPW/bl+fUcDT6J/3uSU9JTKtBqHyx4rcTpwhZX3E9ZdDd9SBy0T75XHFH2e7pvQK9BdFrHufFTXXY2zznYxWq7TJqCnjcYBLgd8YoZl2lJOpXww/zFjcHpa6q1aDw39rVQ+jXXQ9EXUstRzY38rjWONZVTJQ0s9l/p/Y76rY0w1fOscYGj4T9H9mvQO4P8Dnjqa963RNqiPw9W+UxpjGrdRLaPquLeI/WFlbtjS1xaxv61uozqnr47z1XpuKddqmipzgJa+0zKnr8wnW+bc1b7Q0t8q5Vrev63ntV2dig+cHBFXBV5PV8D/j8mn+vxPRFyZbuf4zoi4EPjVDNf/i+huPPP1iHgK3d0hrzUh/MrNFi7or0/5Nt0pM5NU0vT7dI+2eAbwaLpnRz5gwrr/jO6btqfRncZ/FWDNU4h61TKtxjm/z+/RwHER8UO6clpLNXxLHVTrGeDDfbpeSndnz6Q7ZWkW+aiuuxqnpV38LCIuD5weES+he/TJtEeOrNTF//TXL30H2GtGeWgJ35KPeY8ZlfS01BvU6qHa36rlA/U6aOmfy1TPLfVWHcdayqgSp6Weq/1/NM7QfFe3UQ1fLddq+D+m+6XoKcChdNfOPnRC+NE8zKttVMfVatuojjEt26iWUTVNi9gfVvLc0tcWsb+txqnOP6v10DqfrJRrNU2V8C19p2VOX6m3ljl3tVxb+lslTsv+rd28vjFY9hddJU97nM6V6E6x2InudI6nAdeY4fpvSzchux7wZuADwO0nhL8/3aTtVsAn+wZy4CzzvIByL5dpaz3QnbZzIHD5gWmbGr6lDqr1PCb+FYCrFMIPznd13a1xBqxzT2AXum8zDwdeDtxoSpzHA1fr83se3fVdT5pTmQ4K35KPkbjzGDOa01Mom6Z6aOifg8avap5b+uey1XNDnZXGscYyaspDoZ7L7a4h36VtNISv9ue5totFtY2RuNVxeFDbGAlfGmOGbqNaRtU0tYwvLf1haJ7Xs+7Gep7LvITiXLKlHor13DKGVfdvTXkY2neqZVqtt5a+Vi3Xxv7WWq4T+9osXtvFXfEj4taT/p+Zpy7z+ltU0tR/K7dmQ8jMq69a9wenhH/w8JTORnSP6FtTZv5gPeEXJSImll1mfmBV+MH5qK67GmcZ2wU0lWm5nBrStDRjxiLqraF/Lk35rMc887Gs/a1iW6nnZRMRpzG5bUws93loGIdLbaNln76AueHSzTMWkOe5728XsY+uWMZ6rlhU+pet3uZta+3ftpdT8V/W/90Z2I/+5nDAPsAXgDuPBo6InzB5p7j6rofV9X9oyvoPXBX+1VPCP23M4kqartn/73DgIrrr6wP4Q7o7V672T/3fB9JdH/PO/vNBdDer2EJDmVbjnNKHDWAP4If9+6sC/wnsvSp6KXxLHVTrubdy6cO1gDsCn+g/3wM4ge7XkNZ8VNddjdPSLs5ichntMybOn68Vvo/z8lWLqvkul1NDPuY9ZlTSU663fhuVeqj2z1L59Okp1UFL/1yyem7pb6VxrLGMKnFa6rna/1vyXdpGQ/hqf67Ww8rp9k+i+zXt7f3nPwTGPjN6AW2jOq5W20Z1jClvo2EeUJ1nLGJ/ODjPLX2NBexvq3Gq88+GeqjWc8sYVt2/VcKX+07LnJ5CvTXOuavjcEt/q8Qp799mYbs4sM/MewBExLuBQzLzrP7zrYC/GBN+1/7/z6O7NmP0QHfX9a6f7rFEAA8GrkN3YxvoJmSbx4Q/uf97J7rHJRzVf34YXYdcV54z85L+f7+bmb8z8q9XR8TngRevCn98H/7wzLzryvKIOBo4cY30lMq0Gicz9+7DHgEck5kf7T//HnDvMesuhaehDqjXM5l5cJ+ODwO3yMwL+s+7093UsDkf1XVX47S0C7rTrACe3P8dnYD+bI04K3V/U7rTQI/pPz+A7jqv5jy0hG/JxwLGjMHpaaw3KNRDQ/+slg/U21K5f1a3Mc96bqy36jjWUkaVPLTUc6n/96r5rm6jGr5arqXwmXkuQETcMTPvNPKv0yLis8Bzx2xjrm2jYRyuztuq+/SW9lcqo4Y0LWJ/WMlzua8tYn/bsI3q/LM6zlfruWUMq7aNyhygpe+0zOkr9dYy566Wa0t/q5Rry/5t/XKO5/kv24sxd3cct2zkf18Ysmwd6//UkGUj//skcLmRz5cDPjmrPNPdofIR8H+XaDyCyXfF/wqw18jnPYEvT0lPqUyrcejvoLlq2ckT1l0N31IHpXru/3/2qs87rF7Wmo/quqtxGtvFZ4csW/X/Y4FdRz7vCnxshmXaUk6lfCxgzBicnpZ6q9ZDQ38rlU9jHbT0z6Wp58b+VhrHGsuokoeWei71/8Z8V8eYavhqf66GP4OR692B36F/5NLWahvUx+Fq3ymNMY3bqJZRddxbxP6wMjds6WuL2N9Wt1Gd01fH+Wo9t5RrNU2VOUBL32mZ01fmky1z7mpfaOlvlXIt79/W89oufrEf8eWIeAPdN8sJPIru2YhruSQi/hB4dx/+IOCSGa5/U0TcIDPPA4iIvYFNE8L/Bl0DXbne5cr9skkqafoD4NXA6yLiUroD/T+csO5DgU9HxFf7zzcG/mRKeqplWo3zvYh4Nr+e3+9PWHc1fEsdVOsZ4ISI+Djwrj5dj6Qb4GaRj+q6q3Fa2sWVIuLOmfkZ6H5pYvpdSfcAfjny+ZdMvvNpNd8t5VTNx7zHjEp6WuoNavVQ7W/V8oF6HbT0z2Wq55Z6q45jLWVUidNSz9X+D/V8V7dRDV8t12r4xwNvjoid+88/Z/oTE+bdNqrjarVtVMeYlm1Uy6iapkXsDyt5bulri9jfVuNU55/VeqjWc0u5VtNUCd/Sd1rm9JV6a5lzV8u1pb9V4rTs39rN6xuDZXzRXefwZ8AH+9efATtPCL8X8K/A9+iuPT+akV9GZrD++9Jdv3JC/9oM3GdC+IOBbwFv6V/fBB4zyzw3lOkuwG361y4DwpfKtBoHuDrwSuC0/vVK4OoT1l0N31IHpXoeifdg4BX960FTwlbzMXjdjemptovb0P26tLl/nQ7cekqcv+nj/B3d/SFOB541qzw0hi/lo9o/q22pIT2leqvWQ0M7LY9fDXku988lrOdqfyuNY41lNDhOYz239P9qvkvbaAhfrefWfck1GP5En0W0jcq+pNp3SmNM4zaqZVQd9+a+P6zkubrulnpuCd/QlvaiNqevjvPVem4Zw6ppGhy+se+UyrRab7TNuat9oaW/Vcp1rsdhq1/bxV3xl1lEXAG4Wf/xK5l58ZTw16E7nQ66012+M4M0HJqZL4uIVzDmphCZ+eerwt8tM0+MiHE3fyMzjxm3fFvRUgfVet6IZtEuImI3uktBfjRwm7cG7tJ//FRmnjY4wXNUzUdx3eW2NCk9M6q3pauHSh209s+tWc/rrbfqONbY7uY67rW0u4Z8l7bREL46B5gaPiIOysx3RcS4G+uSma+aso25t42Nbh5zsTHbWJr94TKO8Ysy53G+qVwb2sbc8jBvjXPuln1DuYyWsVy3iwP7iHhPZj481ribYW55N8m/ysyXxBp3Zcwt7w5bXf89M/MTscajH3LLx4HcLDO/Ems8OiHHPDKhkqaI+P3MPDoiHrfG+t+4at1/n5nPjoi3jw+ejx6TnlKZVuNExD9m5jNijTv15pZ3Ga6Gb6mDUj33cT6TmXeOLe84Gl2ULe7eOjgf1XVX4zS2i0dl5jtijbuZ5vi7w+6WmT+ONR7Rkls+Oq1api3lVMrHAsaMwelpqbc+3uB6aOhvpfLp41TroKV/Lk09N/a30jjWWEaVPLTUc6n/93Gq+a6OMdXw1f5cDf+nmfnaiHj+GuGfs3rZvNtGwzhc7TulMaZxG9Uyqo57i9gfVuaGLX1tEfvb6jaqc/rqOF+t55ZyraapMgdo6Tstc/rKfLJlzl3tCy39rVKu5f3bLGwv19g/vf97/4mhLrNy7cPJE0O1r/9udI95eMCY/yVbPt7jz4FDuOzRCavD33M9acruoH5H4MaZediA8M/uwx+dme+fFr5XLdNqnJXJ7T9MDNUevqUOqvVMZt65/zv2rqJjDM5Hw7pLcRrbxco1SYPTBPwLXbteeUTLiug/32BVukr5bikn6vmY95gxOD2N9Qa1eqj2t2r5QL0Oyv2zYRtzq+fGequOYy1lVInTUs+l/t+r5ru6jWr4armWwvcH9TsCF+WUX+dHzLVtNIyr1bZRHWNatlEto2qa5r4/pJbncl9bxP62IU51/lmth2o9t4xh1TRVwrf0nfKcvlhvLXPuarm29LdKnJb92/rlnM7xX7YX3bNc/70Y/qVzXP8OwMOL4e805zx/orj+TzekZ3CZVuP0Yd9RXHclfEsdlOp5JM7Eu8G25qO67sb0tLSLPyvGCWCPOeahpZxK+VjAmFFNT6neqvXQ2D8Hl09jnlv657LVc7W/lcaxdYxhg+I01nOp/zfmuzrGtIxJ1TlAtR5OWKa20bhvq87bBo8xjduollHLuDfv/eHgPK+jr817f9vSlqpz+uo4X6nnlnJtSdOg8OvoO9U5/eB6q/a1lnJt7G/Vci3t39b72oHtRHbPav9ZRFylEP42c1z/pcBTCuu/lNq3aeU0AadGxAci4qCIOHDlNSH8xyPiGRGxe0TstvKakp7BZVqN04fdFBGXn1P4ljoo1fNInDMiYo+B4Qfno7ruxjgt7WJSOxsXJ+luQjI0fLVMW8qplI8FjBnVci3VW7+NwfXQ2D8r41dLHbT0z6WqZ+r9rTSOrWMMGxSnsZ5L/X8kTZV8V8eYljGpOgco1QPwmYh4ZUTcISL2WXlN2cbc2kbjvq06bxs8xjRuo1pGLePevPeHg/O8jr427/1tS1uqzumr43ylnlvKtSVNg8Kvo+9U5/SD661xzt3SF6r9rVqupf3bem0vp+Kv+AVwVkQcB/x0ZWGOuRakd1pEHAO8d1X4cacftqz/uIj4C+CoVeG3uL6md2xEPAT4QN94h6ik6dp9mPuNLEtgrZswPbH/e+iq8JM6bLVMq3E2A5/tw4+G3eJamcbwLXVQrWeA3YFzIuKLq+KsNZhsZng+quuuxmlpF5+LiH9iyzLa4jqqEZ+PiNtm5kkTwoyq5rulnKr5mPeYUUlPS71BrR42U+tv1fKBeh209M9lqueWequOYy1lVInTUs/V/g/1fFe3UQ1fLddq+Lv1f0evU03grhPSNO+2UR1Xq21jM7UxpmUb1TKqpmkR+8NKnlv62iL2t9U41flntR42U6vnlnKtpqkSvpp+aJvTV+qtZc5dLdeW/laJ07J/a7Zd3DxvRUQ8ZtzyzHzrGuHfPD54jn0ObMP6v7nG+sddX0N0N5y4EvAruoay5s1FWtM0b9UyrcaJiMPHrSMzn7vGuqvhW+qgVM99nLuNW56ZJ64RfnA+qutujVMREZ8cv/ocdx3VSpwvATehexTKT7msLsb+ItVQpi3lVMrHAsaMcrlWVeqhob+Vx6+GOmjpn0tVz1XVcayxjAbHaaznUv/v41TzXR1jquGr/Xmu7aLfxlzbRsM4XO07pTGmcRvVMqqOe4vYHw7Oc2NfW8T+trqN6py+Os5X67mlXKtpGhy+se+0zOkH11vjnLvaF1r6W6VcF3octl0d2Gu88HF3moFFt4uI2HON7XxrltvZ1q233qyHrWN7H4cX0e6q21iWvhDrfNydNp55tr1ladfbGst1Prb3ct0uTsUPH3c3LU3n9n/PHrf+MfYHTgQeNi45jDl1v1qm1Tjh4+62+uPuaGsXzY/3AX4yLs4689BaTkvzGLSG9JTrrd/G4Hpo6G8+7m56mlr62zbzuDsG9v8+TtPj7oZuoyH8XB93B1yt/7tpSHr6bfi4u+nbWNrH3TG87ZUfdzd03X0cH3e3ZfjV9dxSrj7u7tfXP27OXe0LPu5uA/NxdxNk5tH92y9m5lkDwj+7f/uY7G5uMYSPu9syzhaT4tzgj7vr31baxSwe7xOjycDH3ZXS01hvUKsHH3c34zQ11tu2+Li7if2/t97H3U3bRjV8tVxL4TPztf3bV+Tk+yCM8nF3022Ex91Na3vredzd1L62iP1tQ5xlf9zdkDHMx92NrIphj7ubVq4+7m6jv4A/A65bCP/bc17/jsuUZ+DTwDnA4cDNBoT/FvBa4G7zKtNqHOBBwBXmFb6xDsr1DLwSuMM88lFddzVOY7vYNM86aCzTlnIq5WPeY0YlPS311lAH1f5ZKp/GOmjpn0tTzwuqt5YyquShXM/bwquhP1fDnwt8FHgMsNuS5KE6Dlf7TnmfPu/21zDuLWJ/OO88L2J/W91GdU5fHecXMZ+spqkyB2jpOy1z+nJdL1OZNpTrQvdv283j7nq70d1h8dMR8eSIuPaU8C+PiK9ExPMj4pZzWP83I+LIiLhXRMSUsETEGRHxzIi44YC0lNOUmXcBfpfuNJa3RsRpEXHYhHXfAvgMcGhEnBcR/xgRd5iSnmqZVuMcCHwtIt4eEQdExLSzUkrhG+ugVM+9U4HnRMQ3IuKlEbHflPCVfFTXXY3T0i4+FxHHRsTjIuJqU8ICEBH/Gt2jGa84JHwxDy3hoZ6PuY4ZxfS01Fu1Hqr9s1o+UK+Dlv65TPVcrreGcayljCpxyvXc0P/L+a5uoyFN1XIthc/MGwJ/T/dIqjMj4uiIeOSkOAtoG9Vxtdo2qmNMeRsNZVRN0yL2h4Pz3NLXWMz+thqnOv+s1kN1PtlSrtU0VcK39J2WOf3gemvoay3lWu5vxTgt85h2W/ubkq307cw+wAuArwD/PiXsdYCnAZ8FzgKePav1A7sAD6c7ZW0z8E/AnSeE3xP4K7rTTE4C/gLYY9Z57sPfHHgb8L8D13/1PvwlA8K2lOngOMDl6Aaod9L9mvWGKeseHL6lDqr1PKZcnwAcD3x9Vvmorrs1TrFd3A54OXAe8GHgUVPC343ul8pv0T1q5aHAznPKQyV8KR99nLmMGetIT6XeSvVQbaeV8mnJc2v/XLZ6rtQbxXGssd21xKnsm8v9vyHf1bZdDV+dA6x3X7IUbWMkPZVxtdI2ymNMZRvVMmpJE4vbH07Nc+u6G+t5rvMSivPPhnqozCdb66yapsr+sGX/XJ7TD623xr7Wsm9o2Z9X66E0j2l9zWWly/7qG+FT+0Z45sA4v0l3Hcov57T+qzFwIt2Hv3Ex/NQ09et8NnA63S9ATwV2n7LeOwGv6hv2B4BHFOphcJlW4/SD0wP6NF00YL2l8C110FjPt6O7xuhc4EOzzEd13ZU462wX1yyW0Y50NxJ7D/DjOZRpuZyq+VjQmDE1Peust8H10NA/y+XT2JZKZbos9bzOeqvuS1rKaFCcxvIp9f+WfDeMMeU0NfTnIe3iysAfAh8CvtGPY7+zDG2D+jhcahvVMWYd7a/SjlrSNO/94eA8N7brue9vG+OU5p+VeqjW8zrGsGrbGBS+pZ02lmm1bVTHo5b2WirTYrk2zWOqr7msdFlfwJ8AJ9BdR/5c4BZTwt8c+Du6u8Wf2Me/1qzW38dZ+Wbpm33je8iU8Htx2bdXXwQOnVWe6b4NO5ThZwGcRzdh+CNg14FxSmW6Ks450+IA9wXeQvdN3VuB+wE7TVh3KXxLHTTW84uBrwMfAx4LXHVK+MH5qK67GqexXexGdx3ovwFf67d3mwHxVn4ten9ftq+eYZm2lFMpH5X+2dKWKulpqbdqPTT0z5byKbelhv65NPXcWm/U9yWlMirmoVw+lXa3znyXttEQvtr2Kv1/M/Bq4C5D28W82wb1cbg6b2vZp7f0z8FlVE0Ti9kfVsu12q4Xsb+tbqM6p6+O8y1tr1qu1TRV5gAt6R88P19HvQ3ua419oWXOUCnXpv1b62tuK17GF/AiYN9C+M/T3dXwN+a0/m8CHwQOAq40IPwX6K5NeSZwg3mkqVieV2uIUyrTahzg3cDvM/xGctXwLXVQquc+zpOAaxbKaHA+quuuxmlsF98EXkHtRjhH0U1cj6C7Q+oOMy7TlnIq5WMBY8bg9LTUW7UeGvpbefxqqIOW/rk09dzY30rj2DrKaGgeWuq51P8b810dY6rhW/pzJfzUMll022gYh6t9pzTGNG6jWkbVcW8R+8PBeW7sa4vY31a3UZ3TV8f5aj23lGvL/m3oHKCl77TM6SvzyZY5d8s4XO1vlXKd23HYuFf0G91uRMRvAXfpP346M8+YEPYBwEey9ggoIuJawM4rnzPzP9cIt/LMxaHrvVlmfqWSlj7exDxHxLsy86CIOI3xz5Yc+xzJiHgh8ELgZ8BHgH2BP8vMf6mmcUAeLg/cpP/41cz83wlhrw3ctv/4xcy8cIbpKNdBtZ5H4l2N7tSj0bb0qQnhB+c7Ig4E7tp/PDEzPzQgPYPiRMRv0N319M79ok/RtYtvT1h3ZHEwioj7Asdl5iVTwlWfPVx+fupI3MjMjIhdu6D5/6bnZK5jxuByjYgrAH8M3HJVWg6ZEm9QPYyEn9pOV/IZEVcft46c8Aivah209M9lqufG/lYaxxrLqCXOoPLpw5baXR+nmu9q266Gr/bnQeEj4mWZeWhEfJDxz08e++z5Pu4i2sbUfds6x4DyHCAi7gzcODPfHBGbgCtn5jfXCNs0FxuqcX94H7rrdgf3hz7e1D5XadeL3N/28UvzpIqWeiiuv2UMK6WpIXxL3xk0P2+p68Y5d3Ucbulv1XLdEbg2I4+Zn7R/W4/t6sA+Ip5G91zEleerPgg4MjNfvUb4dwB3oDuV482Z+eVx4UbCP4DuRgq/AVxId9OHL2fmLVeF+6vMfElEvJrxO92nrQr/qMx8R0T8+bjtZubLJ6Rpap4j4nqZef5ad53MzHPXWPfpmblvRPw+8BC6Rzocn5m/NSbsezLz4RFxFuO/PNhnQh7uRnf9yuY+/PXpnt28xeAdEQ+jexbnCX3YuwB/mZnvWxXuM5l554j4yRrp2W1Ceg5gy4Of500IvzPwuDFxHjshzuPpvgW9Ht09D24P/Edmjnt25+B892FfSHdt0zv7RQcBJ2fmMyekZ3CciPg48D66OoPuFOGHZeZ9Jqx/E92pVqvLaGx+R+Ldiu6u4KNx3rYqzJGZeUhEfHLMKnL1Nqrhx6Tn7XQ3hQngIrq2evYa4QeNGSPhS22pUq4RcRTdad2PoLvByx8A56wej8bEuxzdqWb/96UPcMS4nXuhf344M+8fEd9kzPNoM3Pcc35X4lbroKV/zruer0d3CvWd6PL/GeDpmXn+mLCD+1vrvqSYntL+rY9TKp8+zuB2tyre1PE7Iu6ZmZ+IiLEHwJn5gfWEH4k3uFwr4SPidpn5xYi41xrpOX7c8pH4g/dxDWPSoH1b6xhQ2ReOxDkc2A+4aWbepP+y7L2ZeacJcSpldHu6ers5cHm6a39/utY8o2V/GBEnA28C3pWZP1wr3Ej4qX2upV0veH87tC2V5p8R8Y+Z+YyI+BDjx7AD10hPqZ77OHekO9V89IDvbRPCD2obLXlo7DuV+XlTXVfn3H2cqXPDkbCfZHwZTWp7lbnVU+guV/gucOllQdc+7lmXXNCpAcvwAs5k5FQx4EpMv1nIbsAT6U43+Q+6g+Sx1zECZwDXAE7rP9+D7iB6dbgH9H8fM+41JvwT+7+Hj3vNKs/A/zdk2cj/zu7/Hgncr39/+hphd+//7jnuNSUPp9DtcFc+3wQ4ZUIdXGvk8ybgjBm2oSPoBrH/6sv/LOCNU+K8F3g+3Q1CHgMcC7xySpyz6AaL0/vPNwOOmhB+cL77NrHDyOcdB/SDwXHGtYG12sXI/4+lmxx+me56zTcBL54S53Dgk3SD5ZuB7wDvm1VdN7aPzwH3GPl8d+BzU+pt6pjR2pYq5TqShtM7TYkAACAASURBVDP7v5cDPjEgz2+gux7vnv3rzaxxJ91iOw0G3u9jnXXQ0j/nXc/HAQfTTfZ2ojuT4rg1wg7ubzTuS4rpKe3fWsqn2u5G4gwav4Hn9n/fPOb1pvWGbynXlvAtr6Fl1Np/KO7bGtJfngPQHRTGSvvrl625T2woo5OBGwGn0e07DwZeMCF8y/7wRnRfyH6D7pTq+9D/eDehnCb2udZ2vajX0LZEcf5Jf610X/ZbvGZYz2+n25e8lu4LgVcDr5qS50FtoyUPjX1n8Py8sY5b5tyluSHd40BXXnei+8LrJbOohz7sN4BrLKxfLGpDy/BaGQRGPu8MnDUg3jWBZ9B9I/VvdDd9eOqYcCf3f8+gPwiiO5VlrfX+9gzydPlZ5Rk4dcyyNTs18FK6m5CcQXcQcE3gC1PSM24AmrbD2mIHO27ZSn5Xfd5hUh0D9x6z7DHT0jLy98rAsVPSXz5gAk7q/55Of70TEw6OK/mmO0i/+sjnq69Vni1xgE8Aj6SbKAXdL8DT8nvK6nqlO91/WtveYaWN0p3mNOlO/U9m5CYtdHdw/tMJ4R897jUlTVv0lyl9qDpmlNpSpVxXtkt3KvfN+zo+d1J+q3lu6J/lCUJDHbT0z3nXc+VgvaW/bSqWaenLOrpJ7UsL6y+VT0sdrKrjweP3PF8N5VoNf3u6OcuX6G7w9HXga7Mso4YxqbRv6/+/D90juB688poQtjTGjLY1+jkQU370aSijk0fD9+8nfRFY3h+uyu+BwH/THQw9l5F995g0De5zlRfdlz07jnzeje7M17XCn0F3HfUNC9toaUvX6cvnAcB1poTd4iZwjBzEzqCev8yEL1/W2zboxuF3FNbd0ncGz89Xhbkj3VmBE+dW1b62kg8Kc8M11jFt/lmph08y5SaEs3z936kf24k3A1+I7roz6G4S8ca1Akd3TfHBwA3pvlm7XWZeGBFXpOuQq0/h/5+IuDLdxPidEXEh8KsJ6Xl5ROxO9433uzPznEmJj4gTgD/OzM3959vS/WqxxanvI6bmOSKeSHczi5tExOg1LrvSfQM5Vmb+ZUS8FPhBZv4qIn5Bt9OdZH/gr1ct+70xy0adHBFvpKsD6B7fc8oaYT/Wn5r6rv7zI+gmNmv524h4CN2zMa9MV54X0/0SNM7P+78/60/X+z6w94T1A6ycGvo//elB36E79WqS8yPiqsDRwHER8UNgzWtmqeX7hcBp/elHQXcq65qn4TfEeSzdN9CvoTu96fN032xOslJGF/SnXX2b7vS6SX6emZdGxK8iYje60wnXPEUbeEJmvmblQ2b+MCKe0Kd1nNuOvN8ZuBfdTVzWPE0OOC8insNlbfVRdDdZWUt1zKi2pUq5vjG66xUPBz4OXLF/P80lEXHD7C/ZiYgbAGtd21btn5+PiNtm5kkD0rGiWgct/XPe9fy9iHgUl5XTQXRjzTij/Q26M8um9bfP9ac4HwV8IKefultJD5l5SUTcZso6R1XLB2rtbkVp/I7uvhMPYctTZdc65foadH3mzlx2qvzzMnOtsiqVa0P4N3PZ3aSHXsdb3cdV+09p3xYRb6I7sD+HkdNYuezywtWqYwzAeyLin4Gr9vuExwKvnxC+WkY/i+465NMj4iXABXRfHqylZX9IROxDN2e9H90lpO+ka4ufoLsH0qjBfa6hXUPXX74YEQfTHUyv/CK9lgPp6uo9EXEp3dj0npx8HXK1LT0e+Fu68gjg1RHxvMx80xpRPh0Rz8nM9/TxD6UbW2+xRvhqPZ9NVzYXTAiz2uC20Y/DmyLi8pn5ywHrbuk7lfk5ABHxdrrjqtO5bFxKxs+tWubcpblh/Pp9PHag++X+OlO2Uemj5wEnRMRH6I4vgMmXUa/Lor5BWJYXcGvgaXTX5Uz8xZyukd11jf/da8yyK9F9Q7YT3SlpT2PK6Rd0jedpdM81PAt49oSw9wG+Avwp3SlXpwK3Xm+e6X65vBHdFww3HHlNe2TFH4x7rRH2T/r8/ZTu19+V1zeZ8o0icAXgz+l25B+ku5Z/zbt20n258HK6O1Y+aMq6g+6g/uv966Ap4Z8DXJVusvcdugH5+VPiPL4v47vRdfALgScV2uzd6HZ6087OmJpvLrsGavd+nQ9k+rfW5TjVF3B/4CrArei+3TwFOHBKnNf2dfGkvu5OY/IvAmcy8u1431fPKaTxKsAxU8Jcje554qf26XklE+5aTjdm7MCvjxlb/LrS2pZayrWh7u4F/CfddXkn0p3ZdI8J4R9S6J9fotvxn9vX31lMP7tktA5OBf5xSh2U+2djPQ/eNwB7AMfQXbt/Id3Edc8Z19vt+no4D/gw8KiB6bloSHronk98DN01/xN/aa2WT0u76+OUxm+6xzEdRXdwfOjKa0L44/pt7N2/nk13Q7OZ1HND+Iln0M2ojJr3bwzYtwFfasjD4DnASJz96c5CfCmw/4zLaE+6x2/tRneA/HLgRhPCt+wPTwGOp5uHXWHV/z4wJvzgPldt1yPx7k13YPbtSfkdE6/0vPJCW/rqaB7pLkX46oTwu9M9SvS9dF+A/DPdTRUn1fPOhXr+JPBDui/Sj1l5TclnqW30aT6pr78/X3nNqu9QnJ/3cQafqVDta32c6tzwm3Rj1zf78McCd55VPdBwGfV6XtvFzfNijbuqrsg17q4aES/OzL+etmwWIuI36SYPj8jMy08Id3e6QfZ7dAfp3xmw7srdXm/FZXdX/nROOIsgIl438nFnuuscT8kxd9yNiKvQ7fxfCBw28q+frFX+I3GvlavuzBkRN83Mr44JuzdwQWb+ov+8C3Dt7M9yGBP+6nQD365037a9g+7SgKkdo/81Z+fM/NG0sFUR8Q9017B9qRBnT7p6/vf+rJIdM/MnY8KdkpmVX9NKcSJiD+ApbPkr17SzOZpFxF7Abpl55oQwL+3TdATdt8NPAv4rMw8duI3L0R1U3ny96R1Z58LGmAFpGXeTvB/R9em1bgq3A93pvqcAN6X7EugrmXnxuPAj8Xbj19vGWmPwnuOWZ+a31gi/I/CizPzLSduflX5cu3RcP5vDtsb+6tK3/VfQ3egVui+JD/3/2zvz+Num8o+/P9c8l0gKXWPmKRJJg1SGROYhkkJJhuRXKKKEVCRDwnUzEwoZ4xqu2TXcK0PkkkKh0k0I9/n98ax9zz777L3P3vuc73S/6/N6ndf3e85Za+91ztl7rfU8z+f5PEVzXk7/hfCN3E5mNktfBuvHHZfzslmJKGGDc8xBjesup2/p/C3pITNbucYxO+ZJSfea2VpVj9FPyEVPwTfd6UhR4TyZ6T8ga5xqKESHaOCPq66FTedUSe/CnV2GU7wL91aS5kiuteQ7Al6rc/2VHHsW4Otm9tOa/ZYysyd7PX/BsWtf15I2AE7B91Sr4KldX7Tyah1j8drj2+EO3QvN7MddxlbnWroB2DiZR0N0/Soz+0TJ8ffG2YnT8cDPbSXjGG9mO5eNN9PnI3mvm9nNVY9R4Ry5rDsz+15O20HZj0i6GL/G6zAVqs7ZAhYzs2fC87F02RvObBgthv1UWqqqS+AeMuEenT+bWS6tQ9J9lin1JmmyFSgZypVDjwHeSSvf0axY+XQFfALbGqeXXABckjViU+2/g096e+DUtP3xDdzvSj77YVRUew0T2N54FAA8MnuSmRVRlbP93w6cZWafzXmvl9I1jwEddCgz66BDyZVh18tM3LeZ2drZtuH9P+KGwJnBCXAMsJaZrVfQPs84fRnPQyr63WrROUOfL9ESSRqHK92WTWZfxq+LBc1saUnL4grRHarIkk7Cf6fK9OY6fSQ9gHvbp9CiTmIlSszB4fRlOr+jUiNA0ntwL3m6T26pm2CE7olH+oR7ZU+3gpIoaleTHYPT7y4ys2/ltQ991gIOzvkcRXNGpTlGBSrmqeMXqZlX/l4lXYCnH1wZXtoEuBvPtz+3aIMl6Q4zWzfvvZy2ewJH4FGc6bTmyDKaXGXHZGh/o3WpppBp3+T+XBsXy5kvvPQyvmmdlGmXrbrRhpK14SZy0q4sv+LIHbiAaVKxYkdcJK/wNwmOlS3x3Pyl8UjLRdnxp9ofC3wf/92uwdO/9jOzc4rOUQcF39PLeCrYN/IMlrCR3pTO362U3qgaKtSSTgNONLMpFT/HcWHMF4WXtgZWMrPczXVwRO+TM57NM+1yKwyk2udWrpB0a35z2yDn9XS/Ot/R2/Ac2Wz7ojHtg0esKilEBwPxCjxa9zqtOaOnOTXzfpai/RGcap5L0S44R95rRWrsQOm6MMHMPlY03oI+ddNGKu9X617Xoc/d+Bz2cOp8R5nZ8gXt78L1GS7GDfquTooG19KvcCfDb/Hf47P4GvfH0PEnmfbX4xHir+OBnzOBW8zswILjX4uLh1ahvSd9apWXa7pXqjiWJvfOZrieQrIPq1JZagKeGnI37Q7HomoDdSsHNAleNalOUHVvNYGaqvu9YFQY9gkknYrTXK4KzzfGxdO+kWn3FZzuvjSuZphgPtxIzPXISXoCv6lLy+Kl2t+J57JcXObFTLU/AfiWmb0anr8X3+xtVNLnAWANXBRmjfBa7o0qaTJuFP8nPJ8XF/6oVJJB0qy4gdsR0VRv5asWxTeur+Ge2UfwzV5H7WiFEnyZ1x7M2xCH95awjHdX0gYlxuHv8MjYhPDSR/Ec8uXwjcDZOX2uIUQ+SeU5FhlKmb7vww38HfBI3C/NbEJOuwfwaMNdqd95ipmtktP24TDep/HUiColBx/GI2NPdesj6W4z+0C3z5bpcztwK53f0SUlfY7BHWMPp/pY0eJQFxlv+pvA01ZQhirV5zHgm3Q6NZ7OtEvmmKVwmnmC3DlGBV731PE7vO+hX+XvNWxKtrYQfZbXaL8I3yzem+dIC+2+h9PkL7UuC4qkx4F1zezFsnap9pUdk6k+P8apnBfj1ypQWnKs9v0Z5sq9zezW8Hx94OSSDeURuGFyNn7v7IRXVzm2oP2ncHr/z4D34E6W3S2/zu9dZrZOt9cy70/FHbgXmdkdRe1S7ZPSplviOi37AxOK5tXQp3IZtHANPQuch38/2+Npao8BXzGzj+b0uQpfE7L3Wu69EPrk5naWGKEP42lqU6lmVE7DKc5vhbZjaF2DHZtdSQ/imjfZz3Bzpt2uRZ8ptC/ShKmNBt/R7fgamP0MuWMK+6R1rDw/O9v+gJzj9zSnZvo+hu99XgrP34Hvfd6Xafcu/H48B3egJfuY+XFH+vKZ9oua2XOqzzz6AU7zvZD2OaywpnvdeazOfrXudR36zGIZp7mkdxT97mpWr7zutVRrHZW0hZn9JvV8VuDbZnZkwfF/gae+Xk7771bkeN8WT/24CSqXl6u0pqugzF2q/eapto1sntD3CZy+P6Xb+p/qU5mpUHc+Cn1qBa8anqPO3irtZJgT31O9aWYHVRlfXYw28by1zWyv5ImZXS0p7wY9DxeMqEsb/1tVoz6c/4NV24b2+wJImsfMXgmLQqFRH/A/MzNJlvQtaStaghCE/1XQFrkgXzqiuRKtaH927JuFv91EL/L6PhcWrYQO9e08oz7gBUmbm9nlYYyfxdMWivCMXIxoKTM7Qk4jf62k/XRgBTP7Wzj+IjjdbB08B6vDsMdpQZ8uOWYu5BGp5cPjRVw19gBJe5rZ9pnmr5vZ/yQlfWeleFLfuO5YavY5UdKheN5YVfrn3Faf7rUFbvB1o303iprkLTIV8EJy7XVB3TlmbjP7P0nbmNnFNcZT53tdgpZQDfhvN9bM/iup7Ds+AN/0JQKaZR77PwH/rTge8KjyGnguO2b2bHA4lGFBnAGV9oYbxWJbTe7PaYlRH8Y1MWx+i/CpjKF9ijxClWvYm9m1kvaiWtrVjZIOxBlfhju7rpBH5TGzf+f0WSqsC2XrQRqzhb+b4OyhfyRzTQnOxjVhPoWzNHbCnbJ5+HTm+zlN0p1hTj64oM9iZc7IAqwFrFh1A0rNudLMul2bWbxmZj+rcNzGhntwEmWdK0eVdKn7Hc1pZqWMogyewQ3QqvjzAM2pafwFSN+/08I4s/gUXmJwMTx9Jd2+4zq1QDVODHhlUpBKkDAG09F2o31Oy6LuPFZ5v9rgugZYSNJRwHvM7NOSVsSDIkWC1VMl7UgN5hQ1r6W04S5n8M1bMD8m7X+jVIojPg8eX3KKZ8NjDC02V9l9dAhul/w9jGlh4PdAoWFP9TX9uPD3c7iTNGFX7YAHaNLo5d55Bi99XTlKbGY3qzpToe58BF66cU9JVYNXTc5ReW9lnUy42yT1Ld0ii9Fm2L8YDI5z8JttZ3IUZc0pzy8DO0haDfeigXtnyi7yeyVdiBu3aYOmKFK0LH4jrUj7opsbvZaUTIrzAkuEse1pZl8tGVOe2uvpmePOamZv4huxOyUlHqctKVaHB/h56v8kovlUSXsk3WAZenjea5n3EzrUygQ6lKQiOtReuMJrosj+F5wmWISTcWP94/giOg1Xk82l7uOGzt9Sz/8OLBc2um8U9Lld0ipWkc4JIOknuBDMDTh97e7w1jEhupDFzWEDPJekjXDv6xUFh/++mX0+c76zcZGrXJjZ08qhRBc0Xw4XVNqYdhXjMvrnlZI2scCmqYgn8YW2W17jvuHvZjWOXYuqmMJhkk7Hf7fCOSAzx6RzBOeVNK915ghuEuaub+OR6Kqo871eBNwhKXHObY7PH/PgkdNc1Nz0fRu/H+6i/fsp8ozXcUwmx9qtxnigwf2Jqz3/AmdcJcb0TZLWDGPIRtbekrQTLeN7B0pUytVKu9oAT7u6SVJR2lUSTdk38/qe4VxL5PT5oDx3uepacoWkR3HHz1fD/V/mAAUXjdpG0mfNbLyk83BnXx6mh+hVsqHdOvVe0WbrakmfNLPruowjjVoq1GHea9sDmNmD2XZJtDH5/XOOUxRpPSFEEa+j/X5oay/peDPbrygKZ8UU1kREagM8pWsrPLpehrpK3WeHvcWVmc9QtFeqqxD9aLh2rqDCnBrmyefN7HW5JtGqkn5lZv8q+Qx/xasHtVG0FVKgkrEFB8t4SVtZCZssC7WnICW/n1Gs1L27Zajo8qoPZag7j1Xer0r6EF5G7pUQCFkTOD5nnUrjLPyaOyQ8/yPOQCgy7H9Li3FQVaug1rUUrqO98Ll3ErCApJ+Y2Y8K2s9IccQjuovhGj1F+9WHs453SduUjH9MxqB9CXcKlKHSmp4EJiQdae2pN1dIuiXT9uXgmF7FClgkJTgIuCoYqpUU39XJVDhRUhFToUnlgLrBqybnqLy3Uqfq/lp0V91vjNFGxV8Qz8dJLvJbgO8VLUByMak9aEV6tgROM7Pckh2qKRYkaWIYz0/xmpq74b9JUT7eXfiG53Jr0a27ivsEQ++T+A10rZldn3l/Rl6NPJfzw6HtLdaFyhI2eGvji9S9ZvZCQbs58fJZE3D6eprCdrWVCJKpJh0qtJkX/y5LRa2Szy7p/tR3WkbdPxnfKCeT91a48+CbwJWWkxenmnTO0OeLeAnEjuimpAUsk28fvM+70/4755brUSaPKhiWU6yAah3a1NFqeAxY1eoJWSVUv/+FR5U8rUvwXN+sEV1In6oD1UytCX3OwRkWbaWZSuaArwGH0yVHUC78twf+HaWviW46Hsn3+jotBk5Z+3Vw8UwBE82smxFQy1knz7ucSHXa7oE4rX4j3An6RTxi3BHllHSQmR2rgnzkouui4f3ZkQ7Tfqr23Dm5gM8JQHK/TMRz1J8qOH7ttKs6aLKWyDVU/m1eQmkePJWgTGDsbjP7QNhEfhVPRbg7z3EdjJYT8Iie4cbn/rjB9X4zm5jTZ0vcST+GCtd26FM3t3NfPI+ydA8g6TQz26Pguui4HlL9fog7VP9E+/2fvX7eb2aTVFNsSyHtLlnT5GyXS8zsk3ntQ5+639HeeJWef5EyWvN+59C+sphXaF93X/UAvlaNpaU2/j4z2ySvfdmYsmOTtLOZnSPX+cmbY4oMyropSHm5zqV5w3XnsTrfqzz1aDXcyXg2bpx/zsxyr8fQ5x4zWzuzt+pIlUy1ryVUGfrUvZaSlKKd8JJm/4eLwxZ9R5VTHMN7lbQXUu/9CP9O0+XlJltJJDizV0qCSWVr+iPApomjSK7rcVXenlvSufjeusxhk+1zHfAf6qVEPYhXnmhjKuTtu+vOR6HP2ZYTvMp5LXGUztfgHJX3VmqlIIMHQZ/CU3c71rV+YFRF7IMBv6+qKxl/Cc/feQVAntN7Bzm1OINxNNnqKZnOZWY3SJK5l+xwudhN4SJjZs+onQLZtTZtMOSvT8YpaSczOzfVRKm29+ClMbpCXp/0CLzUkIBTJX23YJO+J7Af8G7cU5qc89+06i8XYTFJb7PgcTdnFxTlOC0CHAW828w2VqB/mVmRl/iN8NslEcGFSU1OOdgbpzYlxs+v8I2S4fSfPNSmvpuL+b1H0upkhOGyRn3APmZ2Aqnau5L2Da8lz7+N0wXnkpTQz4QvEKd1GVIdSvRkfKKsbNhbM6pfUh6mFCoWMOtmCNRKrQlYrWjRL8B++KazNEfQXOX9m5J+aznilHmQTxQr1Vmkcdraf/Hvq5tTLHHWLRSMvrSz7t0F3d60GrRdMzsuOCb/jWs8fNcyjskUkt/q3qrHD2iSmvIJKxBdzEMw4Cv9bgG/TIz60L8w7UrOLLnezKZJ+hYeTfuB5USWM2OqvJbIRUnPxDeg/wxr4itF7QNOC9fFofh9Oi8uUJY3lidx53YeijY/P8YdAZVzO3EnWh3sToU9gJntEf7WEjzD59WlrLvY1gvh+HXpm8k19Jo8P/wlymvMQ/3v6ACcndHVaA1r7bxWo2qF1WfgTDezN8N9cbyZnSjp/i7nSAz3eZLfugAJYyiPrVZ2DVZKQZK0PJ42sYDahXrnJ8XqLEDleazBfvVNMzN5auMJZnaGuug+AK/ItQqSvdUHKafN12ZOpX63+fxpYYpmgtnk1W22AH5uZm8osMEKUCnFUa7XtQnwHklpp/P8uCFXNP5vht852U+eZmaXlX2ABnul/XFWQ8IAGYvvx/OwKPCH4IBPayiUrV0LljkKC1CHqXB4zWOD30MzEK73PKfY5ThbMisy+hHcqZyLqnsreaD0GQspyOGe2Qo37CtXvKqLUWXYK6NkLClXyTjdhfbNTiIc0gHzKMbmePS9Kl6TR1ofl0fu/orTfovwjFy50eRq71+nIGdRnsu1Ny70cjlu2O+NR5YfoKWgDLCwSlS3i7zQeC7Omhai9PKySbeRQ983sxMk/Rw42Eoi7QV4F04buw///a4t2cidRT36189wReh3ygVrtsY3oh0Ik8O15qVRKtPwrCKdM3Ouo3EBqTZhOJxlkodd8YhXGl9Iv2ZmPwR+KOmHZvbtquMPqEOJfgdOn8zSrQvL3YWJcidgSTM7UtLiwKLWSkHItp8F9/h2LS3T0GkANVNrAu6UtKJVL1NYN0ewsnEYfq/LyF/QOhDmoK/i94NwGn5ZVYwmzroJkvagk1bbreTo9TmvtcHMrgh/K+cjh/n3d1YzUgQ8IenXeEnKKuJTi+HG4Ifw+3gisK8VizGeGub4s4DzrJxGfLiZXRrWhs/gub+/wMsQFqHyWhKwPc4ouycY+eOA68oMajNLUr5uoYByrIZq7wGPUyO3M8wZ37GS0lZ53ai4Bwjn2BuvIPGv8PzteImsonvoQZwqX6qEjc9BCavuEjPbqtrwuVquWn8cLWGowvuj4Xf0ByrqZoR9Um70smRMlUUYA96QtAOegpc4i2YraJuco2qa4+/CuTuikZKKHFNQPQXpfXjK2Ntod3RNw5kjhbAa6XIN9qvTQmDg88CHw3XSzX44AN97Li3pNmBh2lNsslgf+EKIblZlTq2MMwgWDM9fBHax4jLNv8CNqgeBW+RMqMIce6qnOD6LO5Q3x9fCBNNww7oMt+ERX8Mjxl0RfruEeXyTmV1Z1NbMrpGn/SbCjmVlQdPXtfDfZIcuw/m96qdEXSMX600zFa7ONqo7H+UEr5K5uih49VncJmnTf5L0Ch5gzbUbauytfgF8IhxzA5x1uA/ODjiN8vuhMUYbFb+ukvEBuMGUeNC2wJUWc8UzVFPJNDgaHsEn8SND32OtgP4aDOcT8AtFeF7evnnRPnmu2D/x6MKGeA352UP7BzJtn8MF4IqcFkW0phtxUag3wvPZgGusPF++cmmsTD/hNPPdcJrdRcAZZvanTLta9K/w/vK0SqDdULZRl3Q58HmrUddXFemcmT6V6OxhA7MjPgGnvY7zAW/lTYjKz5c7wUpyq1SPEp37+1t5ubtTCFoHZrZC2BBfZwVlCkOf2qVlQr930r5BLKp5W4sCGvo8gufiVaVDnoFv5kpzBNXOOkju06S6hFkxDa9OmcJGVTEk7VN2LWfa5pWpMyum7TYpv7MwTrHMapcU0aGb0A/no2XsjsEdjhdYgRCTXCfkPFrimjvjdePLKposi99n2+AbvrPyNk/JXCcXqvqDmZ2bnv8Kjp23lny9yMGS6jcGNzyS+/VMfO7o6BfmvXH45vaX+DzzrfRnUA9q75LOwh0GV1M9t7PW/J3ZAwjfCJbtAfKqshT+FvKyhqviLLlCCmhmPSv9bUs+y1w4S7Dbb1z3O7oMN7onUCElSvWrVlyMizDuSEqE0YKYcE77FfE86jvM7Hw59Xg7Mzu65DNUSk0J6/KnLJNCI2cvHmpmSxccv24K0rpWoVpFpk+tCiJ19qtytseOeCrNxGCkjCv6vKl+s+Lrm4DHkr1iQdtalQNCn9uBQyxUCpJrKhxlBeWKi8ZozgLNe68jxRFPico1nMIeeFZgCTMr1KVJtW+iin80nv6aBOd2wNMJysrwrkznelhUvnJ1/LfeFt/HXNplr1or3S/VL81UuMUKmAoN99yVgld593jqvcKUi/B+172VUmm9of0LZnZ4eF5ql/SCURWxp6aSsZn9JCy8ycW3m5mVdfX+UwAAIABJREFUUbpqKZmmLoj/4BvEUphT3Xbq1i5gqeSilIt5vYhPNnmf9zkrVx5tg1x7AODPtMS2DHd8dDMgrpO0FRVKY6URPGTP43mab+KOil9Lut7aS0ZUon+pXczi77Q8h0hasGTz8xowJWzU04thWWSpEp0zg6rCcLfjgh8L4dTUBNNwSnweTgFWk0clDsK9kr/C6Ue5sHqU6GVxo7+O8vE6FrQOwvn+KY8kluEpXF20ammZzfHv6N34b/5e3LG2Ul57q08BBairrv7n8Jg9PHJhzVkHddRha1XFSOF5SfOZU8EPxQ247+dtEK1iVQylyu8Eh0OC+fAIRxnOxTerm+Ib/F0JVOYC5NEPS/Prwjz6S+CXYZN7PvBTeRT/SDN7ItNlYTNLO4rOkrRf2Ycws8fD93kvzixaIzg4D84YQc+FTcPGwPvDfVMqwJS3loTxFCo+S1oVX6c2wRlL5+Jr4414BCKLL5oztT6FM9F2I0T6U+NorPaObzqn0uXeyaDW/N1gDzBGkpK1TR5xKhtbaW53eigF/5dCXtt8T3z8BkyU6wGUrSt117jfUFANpwB1q1bUEWEEd/bsZ2bTw7inAoVG/YwBVEtN2R+4Xi6a9TjMiBLuSMn6Sc0UJOAlSTcAi5jZyuHe29zMvl/Sp24Fkcr7VTN7Xh7I2VGuIzOVgrlC7SkEaSwnqcOBI2n+4BDtlhqbh3ksVf7XzG5SDptQLW2Eot8gd88QrqFfkkpx7IJP4+yY2YElg5F8RMla0kQVfxNg9eT6ljQeuJ92NfsZCA6fj+KG/VX4OjER3/MlbZbDHdU74PfmhYCsQmpR071JuA4uDefPSxFO0GTPfYg8aNWNAVqW3jJXl49QZW81S8pxtCGuk5RgwOzvUWHYq0X9ylUyzmk/J74hXAb3sJ5c5NFLo8pNEI5fub5kaN+ErvhG6v23JE0tMOqh2uY9jYXD32fCY47w/JoKfeuUxvLBuSNhV9w5cTru0XwjeFMfx43T9PGr0L8m0Yp4Qns0tEyt9nfhUQe16JwB/wUeCIt7YRQkeLOfxnNNq6J2vpxqUKLxHK77QhTkTPMyMd1QV+sA8kvLlOFInJ78e/MI58cooZmFxe4UamyurGbahZUIzPQJXfMu1bwqRoLvmNnFcgbUp/CNTVICMnuuufF7dAlzsbFl8QhTlkrYS/mdd4Rrel/znOSbVV5apvZvEK7VTXFjdSzuMDoX/92vwitDpPFi2GgkDsRkA1V0/MSI3hS/5z5jZvfJo3B30G4EbYtv9k4MDrF3U7DJ64IDKN6sT8LF0c7Ao+7JnHSXnAGU2y383QSP7j0o5dfIU02WRXivyb1Ta/6WtDTOgrhPHg38cFhLi1IjrsNTWE7F57K9KFkXrXrO/GpqUUuzGill6+d4fP1IDJMdwmvZcqlp1PqO0s4ZOdNqcSspbdrAYZrsZf4ljzw+T7lOwPZ4tYFL8Ouuik5KpdQUM7tKXvrzaklb4DpMawMbmNk/S45fKwUJ/72+iVN5MbPJwaFRZtjXqiBSZb/a0OArS0nIc+Cch7OAsnuypH1ZNYAn5RVE0kyoPFZY8l3k7RPycuZzS+PO6FDMGDscF9u7KbR7QNLYouPQTBUfnOWbXDsLdGm7NS58eL+Z7SbXoTo90+ZRnPH5mcQpLalbCsEMSHoPHihp04PKaVcnRThBkz33SbSqXR2JB09PorPa1T2SvmwZoWlJu9OeUpGHKpoW5+P7jxdxvZOELb4M9Up+1sKooOKrvoLxhfhiciv+4z1lZqXRldBvAdpV92/GvXVZBfMyz27HYt/N6MqLekh6i5Z3S7j36b/kbARUHqEeckg6Aqfdd1CyJK1gZo+oJVLxvJz+tScuUvEwHl3+R6bf+uaMjTnNrFvZpl7HXyulI/TJ/c2LIlxqp2rPjkf7X8nb8AUj5xqc5vthPJp5f8liVZsSHZwuG+PGyWr4BHemFauA74Q72tbEN55b49TGOqXdSiHpXjNbS67IuoaZTVdQ7i5ofzNhc2XVlcNrpV0Eg+YgOvNHy+oV10I3R4N6qIoR+iRU8B/iQmbnqYAuHObWSXge5MpyevAdlkNJC9fQ5LLvu2A8d5rZB+WpGj/DnT+/thLaqKrX1E3aP4lTj88ws9sz7/0s64CTtAReHjRRfb8dT4vKpZnKleRPx6m0r2be+7yZnZ15bWU8Kgv+GxflmJZ9pmfMbPGc18fgxnxZ7fO8443DN3BL4nPALHg+aEdeolxZ+ULgQFIsizzHoRqWfmsCtSusX4MbZoUK6+G72pNWatd1OHU3V5iwzrzdcPwdc7RKqr40PMdNeG7xrPjm/AXgZstEqNW8asWXcIbIqjjjY17cmfiLkjHNjxuju4VzjcNZZLnBDdVIcwzt18dZCrcD23bbQ6h+ClKTlMK8dLnzsmuPukSvLcV4kzQd3wfvnjL4niwad3h/X3Omzvo2QKrfqXO9HXfMJnPfLbjmSJkmSfYY+2X3YvL0zFcLuhSmB0i6y8zWyfxuZfukJqr4O+AMlAn4tboBnkp2QUH7pDrJJDzKPA3XJlkp1WZL3IGzHj7PXYDPW10ZdnIG6nZk9KDy5mHVSBHO9JuLiukNoX2laldh3b8Mz8FPDPm1wpi2tJKqL6ljlKZ2ylnDi+KppQlrdzlcRDQ3TbtnmFl8ZB745jT5f1bgvor9LsEnmaXC4zCccp5tt0TN8Zwd/u471N9NZlxr4rnud+P0r/uqfFf4zfwBfELaAPd2V+mzajjnmrhoX/r9+3B1TsIxn8UN+yPxTX32eJOSfhU/60XJtYFT3NseFb+rr+O1ptcYhN9mCzzXLO+9d+HRufVT39efCtp+JXzmVzKfeSpwTpcxrIRHbx/DIxAPAj8sab887r39GrBChc+4MJ6fdhVOB74RuLGk/e/xTeGJ+EJ6Ap5DXtT+nvD3/tRrD3QZ02ScHpg8n6fs+sA3kLvj0aGP4DnLx/TxOtgXr9F6RHhMwSsopNvc3+M5rsQjS3/CIwlzAA8WtL035zvNbRveO5f68+VmeBRjZXzzMwlnWhS13xZnvYzH6YlTga27nGPefv1GffiNv4Zvqo4Kjz8AX21wnD+XvHdLg+ONCfPe28Lzd+C6IXltk/l4cuq1mwvavj/8/Ujeo6BPo/mbsD7gzrd9stduwbU3poffsnDebni8X+FU3xnfHXBqlz7L4lTgh/GUsCeBJ0va3x/+fgkvH9z2O6bafSb83TXvUXL8WRp+9oVwcc+ncPbP49m5r8Exp+HpaNNwY+CV1PN/F/QZA3yo5nmuxrVakutva7wscLd+G+Fr4nG4uGxemz3D38PyHpm2W+IOt2fwNXxDYGqXMTwQ/lbdW61Z9ujSd5sqr3U5Rse8l/rez655rDPwtIzJ4T46scL9thWeCvBT3Jiscp5FcWfaZ4F3dWl7Mr4u7xXugftxJkte23nwFK0r8UDgKcAnuxz/MWCOiuNO21az4Eb+fF36fCacY2p4vjquhVHW565w/OR3XJjyeftjuKjdPrjOU5XPsnn4Pl/B9wzTcXZX5etloB6jImKfhqRN6YyOHZFpk63zXViHMtMvTzgn77V0hKyrwq28PunGOHXlo2Ro3DZE0XZJj+IKlFlBmD+V9PkSbmwshnv3P4hH7AojlCFivxu+wcit9auaIhWS7sSNqU1x72QbrDPidgZudP2V/GhDkQe3adRxasF5ymhp2WPcaWa5ytiqKJASWChvpwIlWoHSLemruCL/v/GF7lIzez18F09kP0MP31GlKF+gPS2CX2+v4hutnXDq2O+soCqGpKtxo+lic+/v1njkopCCJafwrW0hgiNP67nHiuveTjKz96e9+pJutpL6wHUgz09f11qe4nnw+23VVJu/UJBjCOViZKH/3Hhu4RTzvPBFgVUsX+jtdnxzeFv4TpfGI2lFrIkb8Uh6nfI7taB6NXVrpUVJyi3v1mqeXyFE0ma4U3IsvkEppFyrhuihyks/zmVmuel5crrrq3QKbRWuPXLtgQ5YPkWzMstC0hJWr4QjkhY1s+dUU6BLnk50PJ4L+xkzm6pywaVzcFZGHRp49hiF83aNY9yP/86z4+kNSamrJfFoXWHEXtJE3Mj7Kb6p3g2QmR1W0H4KLi42Hhcyu6csSpnqNz9+TXcrqzkVdzSMswrVRuTq9F/EDeOzgfFm9vcwTz1iZu9Nte2lKkNlqKZosKSlcNXs9XDjZyoutlkoJJfpvxDwkvVpgx/WjS1wFsTH8d/6soI5/nz8HlgYd/bOeIscfRfVZNRm+taqG19wjA6mkqSHcAfJd3HWXnZQRUKPc+NzRVL+7Vpcc6Vy6d+KY14VXxvS1Peyaj1Jv7HA/FaSKpNquyAu3Lpdl9/gatyZ0q3UYCPbKjANPo6zvZLoezdhuzwG6HfM7KJuY6yKsG/4OJnUTgulT4cSoyLHPoE8721u3DtzOv5j55WXSPLZgLactm75bK+m6Ufy3MM8Ok/aMK9iqJ2K02OWor2sFHTPQRpIvFRlMslgX3yjfqeZfSxQnrrlSm4HLG3l6ud1RSo2w2l3H6d7Lg24B/Y43FN6IW6MFFKHEpjTvR9ssBldK/X/nPgEu2BBW9QuWDMm9LdMm9r5cuZpJC/LRbyeDwb6R4FVJf3K2ilvd+MT6WLA9ua1qdPHmi4XsMueo+l3VDWX+nhccCwxSKYD4yWthefEFeUE7o1vrpaX9Fd8c9WtvN44PO84nXZRVGoRWvmjzwWn47P499cvVNF3mAVnMtTS2lBL9GhOQk5h2Ay8TnEt+cPwuWxxuRr9h3AnUBHSc4MoKb/T1IimXp5j+nN9j+7iZ3k1sefBWRrvwI33PBwPfI5qNdorix5acxHGpBLE3unDUb72pDfEc+IsrWSTlsX3gxPxG3iUa36Ky0TVLv1mZs+FvzMMo4rGz2640/AHwahfEjin5Dw7q0UDHyfPdx5HAQ28yrzdEL2UUZrLzG6QpPB9HS7pVoqv9SNwA2ZiMOqXwiNZuQjz7jg831mS/kV52eFV8XXr9OAEzq1AERy478LXyp8mDiRJH5aLe/5JUraiSd37uSnqigabmX0iGNRjzIVJc2nRcqrv0XjO9ZG4M2MhXMhxFzPL1XgIv9MJeHDFcHr0/tl1OwzmFZw9dW7K4PsWKSHMVNsd5Cr61+JRzW4ftJI2VWbsjerGFw0h57W9cOd/tuxg0r5o37upmR1Cq9wykrbBK0B0INz/x+DioqK7jYGkM/F74g+kAl1FY8pzsEraIM/BmkZw2v4iPMpQSQ8qoIlt9aaZvax2eZbSe8i8MswkWilRWzRxsnbBG2b2kqQxksaY2QR5WsLQw4YBbWCwHgR6WOrvvHjeQ7+OvzpON34qPO4HVstpd1/e/xWOf8pQf4eZ8XwSv+m3wSfwzSmhvIY+Cb35AQJ9h+705kuAd3Zpcwiulv3b8L0nbJRl8OhgUb+O36fLed6Lizzdj0f8vwss16XPjThV7wacdXE5XahEBceZWPLeuNTjl+H7WDjTZjqu+7BM6rVCimWm7wO4g2QZ3Av/U+CqTJvGlO6C7+i3XfrcGf5eizMv1iAnpQCPUBUdY0qFsc1DF7pYpn3ltAtq0sYbfK8HhDnp8PB4AFeMTrepPAdl+l0Z/k7FI4JTU48y6u47wu+1GbBQhfOsDhyLz6kTKKDT4kZh9vFdnGb/n5Lj/yhcQ18Ij6upkA5R93rHDZlDw/dzDCVzWvicpZRuYNbw9yB8Pjo0PCYBB/brGurjtbg4buD2epz78/7v0ueDuPPp0jBPPISLsP0d+HSXvnPhefV1xliJBk7+vF261jX4vlbGDZW9gJUqtL8NdzJcijOWtsRLlfVrPJOBD6eer0+FdLbQdgOcNfcKHolLr2VXkpPqgTtLrqhzXfX7ga9t03H6fil1P7TvmJMJ6So5r9+L78W2waP7HwyvL1/2mYA78br0s4bHzsBdA/UdVPiOZsPXzV+Hx9eA2QraroYz9J6mPaXjc8DbC77/f+c8puGGY9GYdq/5GfJ+t8L1FXiCCmmHmT4P12x/RepxPS7aVpiy2OB32zXv0cfjN0lv6EihyHutx3HVSu0czMeoouKrJWxxJz4BvIRv+pft83nmB7DimsaJsF1a1A4qeOtC/0p1uAca8jIbq+K5eGmK/C4lfS7DoyD74ZGbf+KTd64YUeizFm6wP0R5rd/aIhVyqvTudKZnFNYqT/VdA48erGpms5S0+0je61aiiKxWJQdoRXK+YjVEj5QRhVFvAimJGMlBwKtmdqIyAmnqgdKd+Y5mRGYtJfCS02czXNhncVpRvu+Z2eWZdk+Y2TIFxyh7bw48/20s7ZS33NKQaphSMNAI11K6Xuz9mfcb1cVuMIZCZO/PAnbJgZai0XY533y4Y2V3XAfkx1YiiKeKNXUzfaqmaC2IO1h2wo2RE6xcQRu5iOGRuCMut0a7ehQ9rAu5avhY2u+F3FrIBf2F3x8dFMoQjdwn5/h5Akzpz131N7gXTxtbAGfhbGxmdwbG2PlF17+c1n0cMLuZLaku5atUgwY+GJD0NbxsZFKO7rPASWZ2ckmftXGn9dvwa3AB4Fgzu7Ogfa01VNJtlqmtnvda6r1sBYqzaVWgOMrMlgvtGtekDm1qUbgHAuF6XAl3ZKYZL/Pj1YA61kOlUg0lPWJmK6TeK5zbk/1w5rWe00BSx1oWT+HLVrooEgw8HTfuE4HgzwNvmdmXSs4xm5m9UfR+P6AKNeBTDIJt8bUqwfzAilacalZ43ZeM5wx8PeuallLQf3H8fi6sCDScoPz0hu9biWBl9l4Oc8gUM1uxj+Oah/bUzgWAc61AcHMwMaqo+MCVkt6GR2juw+kc2bIPjSHpKPyG+Vd4/nbgG2Z2aLpdmRHY5fifwQ2nSnW4BwHvr2vEmNmW4d/D5blVC9C9TN54PMLVlsufc+yOjYeZ/bHLsc/GS318CqcU7kROmZsEkmbD84m3x2k+N9MllcDMblZN1W3aa9K/iUd+tu3SJ4u28lXBWLlMrXy5/YFFJJ1CQb5cCm/I1Vh3oUVNmy3TphGlO4ztZnXm/Z/apU9SIu1lPL2mCE1Lmvw2HHsSKeOqZDyVUwo0wLmdwZhM8FR4zHjP2nOjN+zxXFviEYCXw/O3AR81s3Rt6x/ndnYYnfTsRuV3cozoNbsZ0TAjP/FSBXp2t/ZVIVc9/hxuTK5iFfIQA36Al+iZk+I66DPus2DI992Yn3Ei6WzcWH2AlPIxqVrIOX3S1/gYWoy2PPwGj8xcQfcyl01Kv82azG+SjkjWCjN7VPkV+BIcTmf5qjJHaBsNPIGZ/VedNPDEgVWrpGZN7AF8wFr6C0fhSu6Fhn3KKfQf3JjuhkprqGqWHU7hcZzB8iNrr0Dx6wzNuJea1AOKPDo05OpNvA9nMmUp4NPwait5SN8v2dTPjjUmtTZMkPQt3MGf/A51y4qVYRwtrYaPEbQaStqvnQlc3CjPYy7DWHk1lkrOg7pQhRrwAc/izInNad9TTKM4pQjgXnmlmN/Q7sAtS3EdD9wh6fnQJ1e7oAR/wVk8PUHSRWa2rTpLA9YdTze8zzLpDSVj+jbuwE1T/MGZMqf1aTzAjPQUgOmSfkcfNS16xaiK2KcRonFzWqYUXY/H7PCO9tMLrGEm1hA8h8dahRIUwau/F07lnoKXiaqUC6U+ionlHDsp1TXZzFYNhvu11lkCcSM8ergpnkt+AfCb1M1ddo5tcWfSTfhE82Hc+/7r/n6ajvPmlq/KtKkqkLIi/vvdYWbnh83tdmZ2dKpN7Wu9SWRWNXOp1bCkSVkEqGRseWJveeySXVNPO3I7raCsYY1xTIUZNYGXwJkxwjeMf7YKLI0a58oTp+yJBaAG7JKMEX1SNyNaJbmpeDm+Doej2sXn5qYL20peLup13DmXt/nJZWcplGbsMv6eRA/rQNIjeOSp8oYhc42/iZeNva2gbUf0sJ8oi/KXzVuqWb6qwbhql9SsefwpwFoWxLvCvufevOi1pMuzr6WRncNS/aquoRPKD5+//kiat4pDTC7admOBA/eTZrZdTp9a93NTyEszJpihN1Hymdc1szsqHruMATqnmc2WaZ9eG7KwPhrFiTDsDLaEpFvN7MMF7e/DRdj+FJ4vhQtoFu4pVFPoscFnmEKrBvxqYS9xupnl6vKEa1/AcuGlx6yEUSAvCZqFWQljVNITuPM6K1pdJACadbCugSvMd9MLKoWkxc3sGeULkr43x2nV9DwTcCbuxbi2RtdSrpJ+aGbf7sf5c45de98w2BhtEfsOOqGkWnTCLphF0hypRXQuvPRTvzDcxBo+AEwOE03ac5g3EY/HRZ1uxb2eK+JU2SqYFLyyl9Pu1exHDchk0v2XnHL1PH59ZHEwcB5udNatQnAI7o1uU93G88hyIekd+IK1Pj4pT8QpoHWiiV034VZRIMXMHpbXyV0ufE+PpY36ZNg1xpagSWS2liCZmf0NWC84wpIN8+/M7MYu57ld0ipmNqVLuzS6CUEmY5phuMtTJnoy5HOOv2Q49qm4nsNV4fnGuGhkP5EnNFekrj4bXkIxiWDdhBs1bZsfa8Yu+QY+PxwKHJKKxBZt0n9Oi559Ixl6NjlMIqspPmdmRSJ83fB7SZ8s+JwJGjNkGuAhXJjsuW4NFRgrNa/pE0J07Dr6P8dDeZS/LNL7kKQd8bV9WTwH+PZsI3VWG2gTuC0xEOc2s7vVzhqoK/7VAbWEZM8G7pR0SXhrS1pU5yzWxcuanY+Xi6p6XVVaQ62mSFragatOVkWHAxdP77tMrojd4cDNO0fd+7kpsoagAh26pMv9kvamQnqD1WSA9tOp2wWvydPTHpenhPwVF4krwjdxFkEi3jeW7oyRukKPdfGqORPvTXmK7d8pFwxdD4/mP4XfP4tL2rXIyDWzKoyYLP5smXTDLngUXyvAgyfnFzlYa+LmsL/4SZhrkiDKj3HmydplnavCXGT7XTiT87TwO1xoJawmM/u2XKh5xj7DWgzPXlF73zDYGFURexXQCa1/JU0Owqk44/BF/ov4prpsAq9z/N/jm9wf4h6iv+MG43r9OH6D8XSUIgKwnHJ3Ga/trDgdvVJ0t8DTX+jhrwN5+b1LcK2AcfhG+btmVkoDr3mOtvy+sNg9mBc1SbW5HriFlgLzTji9+ROZdo3KV9WFXAl/PKkFCxdIuSXVJkvxrnLcxnn/oX+tXOqKx0yoZbPiYi1P0oDypoolh8oihr0iiZpkXusaDa55jjOBfwEn4d/bPriA0Rdy2tbOo0z1rcQuqTHuRrmpg4FwX8+Ds0wSw6nNQBzI6yZ1jivw33Q+nEp/NyU6J9lxqaJqfXDcfh4X5swtaToUUHt+p2iVryrM76x5/NolNSset7b+gjwPNWGmrYrTss/vFiGru4bKqx8cRmvTfTPutH450+4bOd3nBr6EV0SZt+D4aQfuHyo4cAcdUrHeRHj/Ytwg25FUeoOZVQ2GVB1H1/zxHo6d1WqYH2d43pXT7hkze17OKNkTdz4/D3yrbE8h6Tb82v41bmT9FTjazN7Xp89wMm7EbY87jv+Diz3nGuRyJfYdLbBY5YzE87NrcKp9ZW2RzJjehqctFdL3gxP9R3j65FP4/f9O4EQzO1rSGpbR26kDearx0fjebV9gFZxJcCwu9N0tparJOVfBxWK3M7OiFLVkPfkArsUBPqfd248o/nDeN8wYxygz7GvTCRucY2NaJRauM7Nr+3jsPLGGcxpEkPuGsDCsH57eWrQJyG5CB2NTOlwgpwivinvzwHPZJlum3nqmz4AbZHVQd8FqcPzKdXJD+9qCZDXGUipwZTmUN/VAzxpgw/5anBFxDm6g7QxsYGaf6uM55gG+g2/GhEddv285aSqSHrSMAGTea4MBNaRnDxcMxiZC0peBRfBrKI2PAH81s45SjmqnrVcao6RHcRHSspKmIwKSVsONDXBDurBmtPLrle9sZk/1OIZeU2HmwOfiH+FG94m9jCdz7EtwBkjaubeamX2upE/fHbiDCeXrTTxlBXRoVUxv6HFMufnjZtZLqcT08bcxs4srvHYf8Akz+4dci+AC3NhdHVeMLxxPjvOgVOixF0gaS5ca8MpJ08l7LfXeg7i2SJZWXyasXIm+Ly8DODdewnBaeG1+XAz0LbwaSM/sDUn74qkQz+IVGf7S6zEzx18B3zNvjTMOLgAuKbv/JU0GVk+cC8FpeX/VoEyX8Qz7fcNoM+wvBr5uoabtSIOkY7LGYN5rgzieyoq7auWBAW25YFUrAWxKJy0tV5285meopXxe89jL4KJIt6lddfufuHpmB7Mh1fc4XIzlovDS1nipon5RzGqh7oLV47lKI7OqmUvd4Py19SBUU3Vbg5fbuSDt0bFb8MoBQ+IMVIM8ygEcS63c1MGGulAJ1YAh02AMVwIHZzez8kolh1lOrmnZxqfkPBfi5eCGhbGm5vnm++IiZ0n0bEvgtG6GsVL1yhsMN+94jfQXwnq4KW7Uj8XT3840s7+WnKtu9ZA8TY6O18LrA+bAHUyoht5EaH+3mX1A0i34Hut5nOXYl/z3cI5a+eMNjt9x7xe8NsOxK+kk4AUzOzw8z70uBhNyQcuxtF/bRTXjz8TX9bPDSzvhwp1FEf7a2iJV5315iuyy2UBmMHJfJOxR6pw7c5y34aLW6+BR9E3woOa+1keWjLyK2fk4q+nZin0m4wzXf4TnC+JraD8M+2G9b4BRYtirAZ2w5vEnmtn6ys+z6+cmPW9SHBDjquJ4JgPrWUtxd168jmNfxyPP45kbV1Y9HTdy7zaz3ftw7GtoKZ8n6RmYWZmKd9VjN9kUJ9eQcCpuMqZZ8FrcfbmW6qLugjXAY2kkSFbj+BfSrgfxtHWhQGoE0LP6DUnHm9l+qfm1DXnzqqQNcbpuWx6lmZUJa406SDoaz1FMUwknmdm3BnkctUuIddn45N6fkm5KdFE4AAAXDUlEQVTCWU330Fqbzcw+2/OHaABJL1CSb14UUQtr4rrWKrc6Dy44WhSxGxDHsqTncLX93Dx5M+vQApGXr10ZuBoXqXqo4rlqraGS7sDFYyeG5x8CjjOzdTPtBtSBO9iQNDvVRdUGI0UwcR5MwvdW0/Dyzz1VWVLN0m+SHsKjq2/KmTt7WEjxK5p/mjre6iLse1YF/kB7ilBRKcc5gL1JlU4FTragu5XTfkc83a+ytoikx/F04nHA1VnDPdXujxZKQdZ5ryrkWggnA8dbK8d+9fDa0zaE5fTkFZyOxqtpCHeQH2xm55d2nEkwWgz7UkX1okV6uEDSV3Cv7VJ4DmKC+YDbrEd1yx7GVVlxt8fzJHS05O+8wKVm9smunbsfu28KxHWOXbQpHq6ou2CNZKiBHkRZlLJq1HIgIE+ZOJBOw6EfOervN7NJRfNrel5VD3mUoxEaQCphzXE8YWbL1H2vwXnS15DweWaHXo2MHsbTNN98Cq5781p4PidwT9FcP1CO5SZzTnCYJqy6OhUcaq2hYfM/Hmc3gTPYvmBmD2baDagDdzChCho1QzCmWvnjNY67Gq68/j0gXcFmGjAhy7iQdAjuCHgRr+CyppmZnPE43nLqvDd1vDX4LA9bH2uf5xy/traIJOFr5xfxPPILgbMsU9pZ0m/wffKvMq/vjLPmenKaSlrMCmj3yikv3OD42TJ6M94CpluX9D1Ji+LOcQF3WUH1o5kRo8KwL0JYvLc3s3O7Ni4/Tpq6OxmnrvWsbJs6/gLA23HRvHTEZtpQbIgVFHflYoE74J5lcNrh+WZ2XJ/Pl3iW78Q9+C/hnuVl+3Ds03AxkTrK51WPXXtTLGl589rKuZuyMk9uRH/QxDAfrvQseQ7fqXQaDpMKOzU7z8LhuC8UvN84j3I0YiCphDXHUbuEWA/nWh0XC9sWzzW/1PqY290UqpFvLukAYFe8vCa4bshZZnZ8QfsBcSwPJkuo6Roqz/fFzP7dre1Ih+qLqg1YimDB+cbSJX+8wTGTygxV2n4QL2l2XYrtshwwb96ep6njrS7kJZ1/bGYPV2y/GZ7r/178d+vmFOtJW0QuFHkOzu58EHeS3xHeew+eEvQqvv4bbujOhZf5LUyvGQ5QvtaRgMXw6PsmJX1vMLMNu702s2JUlLsLC8jewHvwnLHrw/Nv4pSWngx72ku5bYLngvdNvdRcLfZlYAdJ6+N5M+MkLSRpSTOb2q9zVcTduFf1WLlifaK4u5cVKO72iCvk+Tw/Au7DJ6h+eQNnBXYLtKLayuddcE+e5zJsiosMqwOAPfCSIVkYLiw3aCjxmgIw2IbGICEpjQW0lccqXKStZsmhQcSbZnbKQBw4RA4Ow1W9BYyR9Ca+yc9uQGdJOSG3w/OOLwEukfTAQIxvhOOHeMmrNJVwQOrydkHtEmJ1EDbv2+Mb9Jfw6JOsZmm0gYA6881/Rit3Phdm9hN5WkHCbNrNypWnm5TUrIIB38A2XUMlHYULnP0rPH878A0zO3SgxzyEmC0x6gHM7I9yQbwi/JYWk2NAWHF5wQN5paOnewlMpfcM6ixTmLtnsJxc72wEOvPeW3gVnWtSjrebJPVV6BHf298h6Xmq7Q+Px4NPU6xa1PRBXPivsraIvBTyznik/2+4g/xy3El+MbAkPsi/AutI+jhukwin7t9Q9VxDCUuJFOc4fS/J6xOCrHMDC4V5JbkA5wfePaADHkYYFRF7Sb/F6V534Ave2/GNyb5m1vOmUj2Ucqt5nsPwTdX7zGw5Se/GBSU6qEoDiUGOBozBlTZvD8/nwCOgL5f37Hrc2srnDc6xCB65+R85m+IialD4zOtaf2qN9gR57eZFcNpbGu8FnrVQez5ieELS4fim4TLac/h6ZvpI2h93ZO6ROBflYninANeY2U9TbWvnUY52DCcqoQaohFigW9+Kl3h7Irz2pPVRJKzhuGrlm6uB4Gbo9zA9lNQcSjRdQ/P2D1VYUSMZqi+qNuBzYmBAromzTIVf75OBd+BBmtyKNBWOO+B7q3Ce2kKPDc7xBB5syarWF13bE4ANrWKpNzXQFpH0R/w6GmcZKryk/zOzY6qce7ijwOl7oJkVXl9y8dL9cCM+LbT3b+CXZvbzgRvx8MFoMezThneiCLmE9U+BdlByakNkaw3gPmuVExp08Tw1VNzt4Xx3WEZYp8/HXxOPsBiuWdBXunuTTfFAf+aqUAMBwIjhA0l5bB7rh+Ek6X5gIzN7MfP6wjilco3Ua7XzKEc7ApUyoXQCYEOYkzsQkLQlvnlbD4/AXYArc/dchqnHcdXKN1en4OZTZrZfhfPkblL7ZfwMJHpwZkzGdQgSbZ65cG2eIdFTGAwoX6PmpCIKdtP0hppjugA4MqGvS1oRZ7EeiafB9KxGH67vZc3s9+F3nrUf++66jrceznOj1dCjkWvJHAncTLsjvagKRW1tEUmqyAYY0ejF6Stpnz4zN0YURgUVH19wAafwSJraL6M+oDZ1tyH+FzbDCc1pnj4dty5mwVVacxV3BwDXSdoKX2z6OqFJ+i5eWi2hWI6TdLGZfb9f5zBX/K6r+j1gn7kmxmaNegAzu1eelxcxjDHABtJsWaM+nPOFLM3UzH4g6QZaeZTpms77DOAYRyQkHYOnLLSpMeMGwUwDM7sMp/rPg+ej7w8sIukU4LKmUcM+jGtMzS4rpoIHZ+DpaoVoahQPM6RTEDfG66FXSUE8B7hBXo/bcBGw8eVdRjz2CsbdDAMvRBdPSDdqmt7QEMtbKifdzB6WtIaZPZlHoa8LSV/G0woXBJbGc6NPpT9pIp/HHW/LAV9Pjbffe+5HJZ0HXEG7oV6UkvMDXIRwTpyZWQozuzmHZp5b+UCpSgAFKQ59qQQwjLAV7vSdIBcZvYDqNsfLknbJvmgZIcGZFaMlYt9TDfXhAkkH4rS9jfAczC8C5w22Z2qwaXPyEnDz4Oq4r9HH303SI8Aa1lIxngtnRKxQ3nNgkfrMb+HiJ0NyrWqQVLEj+gtJB5nZseH/bczs4tR7R5nZwX04R+E8MLNTawcakh7DRZVmuqoT3SAXCtwG2K5OtGwoUZe1lxPh71pSc7hBPaQgykuibYiva9eZ2bUDN9KhR971UJCSMCg09nCuC4F/4AYTuCNxIdxonmhma/d4/Adw1fa7UgzTkVYNaFzOy2bF5e7uNbO1Khy3Cc18UCoBDDeknL474BpT4+ni9JWUtonmxOea+2yUiPSOCsN+ZoKkjYBP4jf1tWZ2/RCMYaapyy3papz6lAj5vA04x8w2G9qRDQ9oEFWxI/oHDUL5vYzDtO0thrASwMyAMC9tYyO8dvdoQd3gQS9G8XDBYKUgjmTI62nviFOsb029NT8ubPqJkr4DnSI4F15GOUkPmIjXIH8NmLvXuUfSXWa2TrJfDNf5fX1mHQwrSDoa3y+VMo2a0Mw1SJUAhjOaOn3llcXOnglZDbmIhv0IhaSFgJeGgqYtaUEb5DJ7coXLZXHvG9CffFN5rc+18UoJhk+cEwkqpWb29V7P0XBcwgV2ljSzIyUtDixqZqUUzwEYRyMBwIihRdr5lnXEzUyOuZkNIdJgeAWX1YAbaKeADsl8FNFfzAxGcQNnxkQzWz+w0UZ0XfqqCBH4JckpVQxMLkq/yEkR3AIXSu5bimA4z1y43tRjXRvXP/axwL+AXfB0q68CD5vZIf0+10BB0pL42MfSrnWSayCmmJav44yconuhJ20R1SjBGQEhNXDyUDNxBwvRsB8BkNf4PBqnTR2JK2IuhOen7mJm1wzh8AYckr6E5+4thpcn/CBwRz9ompJ2LXvfzIYk9y/kmE4HPm5mKwTHxnW90uN6GM+AqGJHDAwGI2If0X90mY9stOQIzuyYWdIDI6ojOMmT9ftuMysscTYYKYKSNscNw9nNbMmQ631Ev6Ka8uo+u5NimOLG64gxOiQ9CJxBpyp+B+09BGMWN7M/1zh+LZq5BqESwMwASVfQch6OwfU/LjKzbxX3mnkQDfsRAEn3AgcDCwCnARub2Z2SlsepODN19C0IyqwN3Glmq4fP/b2ZmQaeGF+ZyOuDZrbaUI8tYvgjZTikjQaINPkRAUn7mllWWKvjtYiIkQK1CwZOxg2SkSYY2AiStgGOA27C5+APA980s18XtB/wFEFJk3Bj8iYboCpL8gopmNkL/TrmYCJJJ6jRfpKZvb/huUpp5hqkSgAzA9RebeBNXMPkL0XtZzZEw34EQNIDFkqPSHok7bUdDbRaSfeY2dpBjGUdM3s9/Z00PGaiPpsHG2oDWtJdOE3rnmDgd5QQi4iImDlRVWwrImKkYGYQDGyKEPndKInSh/X890X7jMFIEczmwIfXejbsQ+T6MOBruBNDuAjwiWZ2RK/jHkxI2hFPAb2O9pSoXL0DSScBZ5nZPQMwllolOCNA0vy0p1AMagrxUGG0lLsb6Zie+v/VzHujwTPzl+Cx/g1wvaR/As/2eMw8z7dwun/PiuF9wM/w3PZ3SvoBsDVw6NAOKSIiYiCREttaSqnyRsB8uHpyRMRIRa2SgDMZxmSo9y/hFOEiXBYeCW4agDE9FAzXWSQtC3wduL0Px90P+BCwtplNBZC0FHCKpP3N7Kd9OMdgYRW8SsDHaS87WpQG+jFgT0lP02LMWT9YEFa/BOeohaQ98LTlV/HfTfjvVihOODMhRuxHACKttoVAsVkAuMbM/tenY2briF5iZj/vx7F7QUg5SEoC3WBmjwzxkCIiIgYQTcW2IiKGO0az1oekH+FK5ueHl7bD7+f/G8IxzQ0cQnsO/JFJXn8Px70fZye8mHl9xLEOJT2Klx2ttNdUQblC62OZwojukPQ4sG72GhwtiIZ9xLBFJidvCnBGvza2alBHdDAwkJ85IiJi+COUNbrWSkphRUSMNIx2wUBJn6NVWu4WM7ssp82wThGsAkkPmdnKdd8bjgjpI/uUCR0W9Hsn7RWcKgvqRfQOSdcAnzOz/3ZtPBMiUvEjhjPG056TtyKujt8PPBqO+xlr1RHdv0/H7gXZz7wCTm2LiIgYBTCztyT9V9ICZvbyUI8nIqIfMLNZhnoMQwkzuxS4VKFUcUGzAU8RzKT4dKAPqvhl0e2+sCwHEYsAj0q6h/Yc+6Jyd5sDPwbejeshvBd4BFhp4IcakcK3gduDVtWoKxcbDfuI4YyBzMnbCo/YTwjevQvwBXSoMZrzECMiIhyvAVMkXU8ryjlqNiYRETMDykoVS+ooVZymbOelCPZpWOsCz+BpAXfR/33PapL+nfO6SEWxRwgOq9n+SLwc8+/NbI1QJniH/g8rogt+AdxIpkzhaEE07COGM95I/jGzN11stT8INLjLUnVE9wcWCfXjC+uIDgIG7DNHRESMGEzAWTvTcUXprGhqRETE8MfPaZUqvpFMqWKgzbAvSBGUmX2sj2N6F66ynwh1/g4vm/yHfhx8ZmJm5NWr74I3zOwlSWMkjTGzCZKOGZDBRZThTTM7YKgHMVSIOfYRwxaDnZPXrY7oYGC05yFGRIxmSJoVOAr4IvA0rpy9ODAOONjM3ijpHhERMYxQt1RxKGl2K7B7KkXwSTMbEDVvSXPgBv6PgCPM7MSBOM9IRWBcnIinRM4OzAK8UrQPk/R7PFD0Q5yZ8Xe8OsB6gzPiCIBQSepp4AraqfijotxdNOwjIiIiIiKGAST9FC9tt7+ZTQuvzQ8cB/zXzKLeRkTECEFa+b9KVQBJW+IR+/XwaP4FwOlmtmSfxzUHsClu1I8FLgfONLO/9vM8Ix2S7sV/j4uBtYBdgGXN7OBMu2XwfPwHcHbVGGAnPMf+d2Y2aTDHPdohaWrOyzZQDrLhhmjYR0REREREDAOEMj3LWWZhDkr5j5rZskMzsoiIiLpoWqo4lSK4A14zfTx9ShGUNB5YGbgauMDMHur1mDMrJN1rZmtJmpzUopd0ezYCL+lKnFE1OfP6WsBhZvaZwRv16IakMXipu9uGeixDhWjYR0REREREDANI+qOZLVf3vYiIiJkT/U4RDHT/JN0vbQDEdL8MJN0CfAI4A3guPL6QLTvYpcTflEQQOWJwIOkOM1t3qMcxVBgz1AOIiIiIiIiIAOBhSbtkX5S0M16iMyIiYhTBzP5hZr/ol+6PmY0xs/nCY/7UY75o1Hfg87idtDfuDFkMr6iURZna/1wDMK6IclwnaSuNUvXpGLGPiIiIiIgYBpD0HuBSPE9zEh5RWxvfHG4Zc2AjIiIiBhaSPgssZmYnhed3Ae/E5+ODzOzXmfbnAzea2S8zr+8OfNLMthuckUcASJoGzEOrosyoYqNEwz4iIiIiImIYQdLHgZXwDckfzOyGIR5SRERExKiApNuA7c3smfD8AVzrYF5gnJltmGm/CHAZ8D/cIQsutjc77pB9frDGHhER69hHREREREQMI5jZjXjd64iIiIiIwcXsiVEfMDGUSvtHEDZsg5n9DVhP0sdwYUJwNfw4hw8BAgV/J2BJMztS0uLAomZ29xAPbVAQI/YREREREREREREREaMekp4ws2UK3vuTmS092GOKqA5JpwDTgY+b2QqS3g5cZ2ZrD/HQBgVRPC8iIiIiIiIiIiIiIgLukvTl7IuS9gRGRdR3hGMdM9sbeA3AzP6Jp0WMCkQqfkREREREREREREREBOwP/EbSjsB94bX3A3MAWwzZqCKq4g1JsxDKOUpaGI/gjwpEKn5ERERERERERERERERASsQUXMQ05syPAEjaCdgOWBMYD2wNHGpmFw/pwAYJ0bCPiIiIiIiIiIiIiIiIGPGQtDywIV5Z5gYze2SIhzRoiIZ9RERERERERERERERExIiEpDmBvYBlgCnAGWb25tCOavARDfuIiIiIiIiIiIiIiIiIEQlJFwJvALcCGwNPmdl+QzuqwUc07CMiIiIiIiIiIiIiIiJGJCRNMbNVwv+zAneb2ZpDPKxBRyx3FxERERERERERERERETFS8Ubyz2ik4CeIEfuIiIiIiIiIiIiIiIiIEQlJbwGvJE+BuYD/hv/NzOYfqrENJqJhHxERERERERERERERERExghGp+BERERERERERERERERERIxjRsI+IiIiIiIiIiIiIiIiIGMGIhn1ERERERMQogqS3JD2QeoztwzH3krRL+P8sSVv3esyIiIiIiIiI6ph1qAcQERERERERMah41cxW7+cBzezUfh4vIiIiIiIioh5ixD4iIiIiImKUQ9JYSbdKui881guvf1TSzZIukvRHSUdL2knS3ZKmSFo6tDtc0oGZY24o6bLU840kXTq4nywiIiIiImJ0IBr2ERERERERowtzpWj4ieH9d2AjM1sT2A74War9asC+wCrA54HlzOwDwOnAPiXnuRFYQdLC4fluwLg+fo6IiIiIiIiIgEjFj4iIiIiIGF3Io+LPBvxc0urAW8ByqffuMbPnACT9CbguvD4F+FjRSczMJJ0N7CxpHLAusEufPkNERERERERECtGwj4iIiIiIiNgf+BsenR8DvJZ67/XU/9NTz6fTfR8xDrgiHO9iM3uzL6ONiIiIiIiIaEM07CMiIiIiIiIWAP5iZtMl7QrM0o+Dmtmzkp4FDgU26scxIyIiIiIiIjoRc+wjIiIiIiIiTgZ2lXQnTsN/pY/HPhd4xswe7uMxIyIiIiIiIlKQmQ31GCIiIiIiIiJmUkj6OXC/mZ0x1GOJiIiIiIiYWREN+4iIiIiIiIgBgaRJePR/IzN7vVv7iIiIiIiIiGaIhn1ERERERERERERERERExAhGzLGPiIiIiIiIiIiIiIiIiBjBiIZ9RERERERERERERERERMQIRjTsIyIiIiIiIiIiIiIiIiJGMKJhHxERERERERERERERERExghEN+4iIiIiIiIiIiIiIiIiIEYz/B/w+MWHUA74bAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1008x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "## load the dict from file\n",
    "with open(\"RNA.family.pickle\", \"rb\") as handle:\n",
    "    RNAfamDict = pickle.load(handle)\n",
    "\n",
    "## get taxon counts per family \n",
    "taxCounts = {k:sum(1 for x in v if x) for k,v in RNAfamDict.items()}\n",
    "with open(\"RNA_virus_family_counts.tsv\", \"w\") as fo:\n",
    "    writer = csv.writer(fo, dialect=\"excel-tab\")\n",
    "    for key, value in taxCounts.items():\n",
    "        writer.writerow((key, value))\n",
    "\n",
    "## bar plot of the counts per family  \n",
    "def barPlotFromDict(inDict, xlab, ylab, fname):\n",
    "    sortedX = sorted(inDict, key=lambda k: inDict[k], reverse = True)\n",
    "    sortedCounts = sorted(list(inDict.values()), reverse = True)\n",
    "    x = np.arange(len(inDict))\n",
    "    xlabs = [elem[:20] for elem in sortedX]\n",
    "    fig, ax = plt.subplots(figsize = (14, 5))\n",
    "    plt.bar(x, sortedCounts)\n",
    "    plt.xticks(x, xlabs, rotation = 90)\n",
    "    ymax = max(sortedCounts) + 20\n",
    "    plt.ylim(0, ymax)\n",
    "    plt.tight_layout()\n",
    "    plt.margins(x = 0.01)\n",
    "    ax.set_ylabel(ylab)\n",
    "    ax.set_xlabel(xlab)\n",
    "    plt.savefig(fname, bbox_inches='tight', dpi = 300)\n",
    "\n",
    "## bar plot of counts by taxon\n",
    "barPlotFromDict(\n",
    "    taxCounts, \n",
    "    xlab = 'Family', ylab = 'Number of RefSeq taxa', \n",
    "    fname = 'RNA_family_counts.png'\n",
    ")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4627\n",
      "\n",
      "            /-YP_009329842.1|Hubei_narna_like_virus_24\n",
      "         /-|\n",
      "        |  |   /-YP_009336767.1|Hubei_narna_like_virus_23\n",
      "        |   \\-|\n",
      "        |     |   /-YP_009333146.1|Beihai_narna_like_virus_26\n",
      "        |      \\-|\n",
      "        |         \\-YP_009336749.1|Hubei_narna_like_virus_22\n",
      "        |\n",
      "        |                     /-ALD89118.1|Rhizoctonia_solani_mitovirus_13\n",
      "        |                  /-|\n",
      "        |                 |  |   /-YP_009182161.1|Botrytis_cinerea_mitovirus_3\n",
      "        |                 |   \\-|\n",
      "        |                 |     |   /-ALM62240.1|Soybean_leaf_associated_mitovirus_5\n",
      "        |                 |      \\-|\n",
      "        |                 |         \\-YP_009272900.1|Fusarium_poae_mitovirus_3\n",
      "        |                 |\n",
      "        |               /-|      /-CAJ32466.1|Ophiostoma_mitovirus_1a\n",
      "        |              |  |   /-|\n",
      "        |              |  |  |  |   /-ALD89122.1|Rhizoctonia_solani_mitovirus_3\n",
      "        |              |  |  |   \\-|\n",
      "        |              |  |  |     |   /-YP_004564622.1|Tuber_aestivum_mitovirus\n",
      "        |              |  |  |      \\-|\n",
      "        |              |  |  |        |   /-AGT55876.1|Ophiostoma_mitovirus_1c\n",
      "        |              |   \\-|         \\-|\n",
      "        |              |     |            \\-CAJ32467.1|Ophiostoma_mitovirus_1b\n",
      "        |              |     |\n",
      "        |              |     |   /-ALD89117.1|Rhizoctonia_solani_mitovirus_12\n",
      "        |              |     |  |\n",
      "        |              |     |  |      /-YP_005352912.1|Clitocybe_odora_virus\n",
      "        |              |      \\-|   /-|\n",
      "        |            /-|        |  |  |   /-AIF33766.2|Heterobasidion_mitovirus_1\n",
      "        |           |  |        |  |   \\-|\n",
      "        |           |  |        |  |      \\-AIS37555.1|Cytospora_ribis_mitovirus_1\n",
      "        |           |  |         \\-|\n",
      "        |           |  |           |      /-ALD89120.1|Rhizoctonia_solani_mitovirus_15\n",
      "        |           |  |           |   /-|\n",
      "        |           |  |           |  |   \\-AMM45292.1|Macrophomina_phaseolina_mitovirus_3\n",
      "        |           |  |            \\-|\n",
      "        |           |  |              |   /-ALD89119.1|Rhizoctonia_solani_mitovirus_14\n",
      "        |           |  |               \\-|\n",
      "        |           |  |                 |   /-ALD89126.1|Rhizoctonia_solani_mitovirus_7\n",
      "        |         /-|  |                  \\-|\n",
      "        |        |  |  |                     \\-ALD89128.1|Rhizoctonia_solani_mitovirus_9\n",
      "        |        |  |  |\n",
      "        |        |  |  |   /-AIS39016.1|Leucostoma_persoonii_mitovirus_1\n",
      "        |        |  |   \\-|\n",
      "        |        |  |     |   /-AAO21337.1|Sclerotinia_homoeocarpa_mitovirus\n",
      "        |        |  |      \\-|\n",
      "        |        |  |        |   /-ANJ77670.1|Sclerotinia_nivalis_mitovirus_2\n",
      "        |        |  |         \\-|\n",
      "        |        |  |            \\-AHX84134.1|Sclerotinia_sclerotiorum_mitovirus_6\n",
      "        |      /-|  |\n",
      "        |     |  |   \\-YP_009272901.1|Fusarium_poae_mitovirus_4\n",
      "        |     |  |\n",
      "        |     |  |      /-AHF48627.1|Sclerotinia_sclerotiorum_mitovirus_11\n",
      "        |     |  |   /-|\n",
      "        |     |  |  |  |   /-ALD89101.1|Macrophomina_phaseolina_mitovirus_2\n",
      "        |     |  |  |   \\-|\n",
      "        |     |  |  |     |   /-ALD89127.1|Rhizoctonia_solani_mitovirus_8\n",
      "        |     |  |  |      \\-|\n",
      "        |     |  |  |         \\-AIT71973.1|Rhizoctonia_cerealis_mitovirus\n",
      "        |     |  |  |\n",
      "        |     |   \\-|      /-AHY03257.1|Buergenerula_spartinae_mitovirus_1\n",
      "        |     |     |     |\n",
      "        |     |     |   /-|      /-AHF48628.1|Sclerotinia_sclerotiorum_mitovirus_12\n",
      "        |     |     |  |  |   /-|\n",
      "        |     |     |  |  |  |  |   /-AHE13866.1|Sclerotinia_sclerotiorum_mitovirus_7\n",
      "        |     |     |  |   \\-|   \\-|\n",
      "      /-|     |     |  |     |      \\-AHF48623.1|Sclerotinia_sclerotiorum_mitovirus_7\n",
      "     |  |     |      \\-|     |\n",
      "     |  |     |        |      \\-CEZ26304.1|Grapevine_associated_narnavirus_1\n",
      "     |  |     |        |\n",
      "     |  |     |        |   /-ALD89121.1|Rhizoctonia_solani_mitovirus_2\n",
      "     |  |     |        |  |\n",
      "     |  |     |         \\-|      /-ALD89116.1|Rhizoctonia_solani_mitovirus_11\n",
      "     |  |     |           |   /-|\n",
      "     |  |   /-|           |  |  |   /-NP_660174.1|Cryphonectria_parasitica_mitovirus_1_NB631\n",
      "     |  |  |  |            \\-|   \\-|\n",
      "     |  |  |  |              |      \\-AAR01973.1|Cryphonectria_cubensis_mitovirus_2a\n",
      "     |  |  |  |              |\n",
      "     |  |  |  |               \\-ALM62241.1|Soybean_leaf_associated_mitovirus_1\n",
      "     |  |  |  |\n",
      "     |  |  |  |   /-AGT55877.1|Ophiostoma_mitovirus_7\n",
      "     |  |  |  |  |\n",
      "     |  |  |  |  |      /-YP_009337200.1|Wenling_narna_like_virus_9\n",
      "     |  |  |  |  |     |\n",
      "     |  |  |  |  |     |         /-AHL25281.1|Rhizoctonia_mitovirus_1_RS002\n",
      "     |  |  |  |  |     |      /-|\n",
      "     |  |  |  |  |     |     |  |   /-YP_009165597.1|Binucleate_Rhizoctonia_mitovirus_K1\n",
      "     |  |  |  |  |   /-|     |   \\-|\n",
      "     |  |  |  |  |  |  |   /-|      \\-YP_009249807.1|Rhizoctonia_oryzae_sativae_mitovirus_1\n",
      "     |  |  |  |  |  |  |  |  |\n",
      "     |  |  |  |  |  |  |  |  |   /-AEP83726.1|Tuber_excavatum_mitovirus\n",
      "     |  |  |  |  |  |  |  |   \\-|\n",
      "     |  |  |  |  |  |   \\-|     |   /-BAN85985.1|Rhizophagus_sp_HR1_mitovirus_like_ssRNA\n",
      "     |  |  |  |  |  |     |      \\-|\n",
      "     |  |  |  |  |  |     |         \\-BAJ23143.2|Rhizophagus_sp_RF1_mitovirus\n",
      "     |  |  |  |  |  |     |\n",
      "     |  |  |  |  |  |     |   /-APG77166.1|Shahe_narna_like_virus_6\n",
      "     |  |  |  |  |  |      \\-|\n",
      "     |  |  |  |  |  |         \\-YP_009336548.1|Shahe_narna_like_virus_5\n",
      "     |  |  |   \\-|  |\n",
      "     |  |  |     |  |            /-NP_660181.1|Ophiostoma_mitovirus_6\n",
      "     |  |  |     |  |         /-|\n",
      "     |  |  |     |  |        |  |   /-AEY76153.1|Gremmeniella_abietina_non_host_specific_mitochondrial_RNA_virus_S1\n",
      "     |  |  |     |  |        |   \\-|\n",
      "     |  |  |     |  |        |      \\-YP_009009144.1|Sclerotinia_sclerotiorum_mitovirus_6\n",
      "     |  |  |     |  |        |\n",
      "     |  |  |     |  |        |         /-ALM62249.1|Soybean_leaf_associated_mitovirus_4\n",
      "     |  |  |     |  |        |      /-|\n",
      "     |  |  |     |  |        |     |  |   /-AHI43534.1|Fusarium_circinatum_mitovirus_2_1\n",
      "     |  |  |     |  |        |     |   \\-|\n",
      "     |  |  |     |  |      /-|     |     |   /-AHI43535.1|Fusarium_circinatum_mitovirus_2_2\n",
      "     |  |  |     |  |     |  |   /-|      \\-|\n",
      "     |  |  |     |  |     |  |  |  |         \\-YP_009272898.1|Fusarium_poae_mitovirus_1\n",
      "     |  |  |     |  |     |  |  |  |\n",
      "     |  |  |     |  |     |  |  |  |   /-YP_077184.1|Gremmeniella_abietina_mitochondrial_RNA_virus_S2\n",
      "     |  |  |     |  |     |  |  |   \\-|\n",
      "     |  |  |     |  |     |  |  |     |   /-ALM62242.1|Soybean_leaf_associated_mitovirus_2\n",
      "     |  |  |     |  |     |  |  |      \\-|\n",
      "     |  |  |     |  |     |  |  |         \\-AKN79252.1|Alternaria_brassicicola_mitovirus\n",
      "     |  |  |      \\-|     |  |  |\n",
      "     |  |  |        |     |   \\-|      /-NP_660180.1|Ophiostoma_mitovirus_5\n",
      "     |  |  |        |     |     |     |\n",
      "     |  |  |        |     |     |   /-|   /-YP_009182160.1|Botrytis_cinerea_mitovirus_2\n",
      "     |  |  |        |     |     |  |  |  |\n",
      "     |  |  |        |     |     |  |   \\-|      /-YP_009121785.1|Sclerotinia_sclerotiorum_mitovirus_1_HC025\n",
      "     |  |  |        |     |     |  |     |   /-|\n",
      "     |  |  |        |     |     |  |     |  |   \\-AEX91878.1|Sclerotinia_sclerotiorum_mitovirus_1\n",
      "     |  |  |        |     |     |  |      \\-|\n",
      "     |  |  |        |     |     |  |        |   /-AHX84132.1|Sclerotinia_sclerotiorum_mitovirus_5\n",
      "     |  |  |        |     |     |  |         \\-|\n",
      "     |  |  |        |   /-|      \\-|            \\-ALD89135.1|Sclerotinia_sclerotiorum_mitovirus_18\n",
      "     |  |  |        |  |  |        |\n",
      "     |  |  |        |  |  |        |      /-YP_002822229.1|Thielaviopsis_basicola_mitovirus\n",
      "     |   \\-|        |  |  |        |   /-|\n",
      "     |     |        |  |  |        |  |  |   /-YP_009272899.1|Fusarium_poae_mitovirus_2\n",
      "     |     |        |  |  |        |  |   \\-|\n",
      "     |     |        |  |  |        |  |     |   /-YP_009126873.1|Fusarium_coeruleum_mitovirus_1\n",
      "     |     |        |  |  |        |  |      \\-|\n",
      "     |     |        |  |  |         \\-|        |   /-YP_009126872.1|Fusarium_globosum_mitovirus_1\n",
      "     |     |        |  |  |           |         \\-|\n",
      "     |     |        |  |  |           |            \\-AHI43533.1|Fusarium_circinatum_mitovirus_1\n",
      "     |     |        |  |  |           |\n",
      "     |     |        |  |  |           |   /-NP_660179.1|Ophiostoma_mitovirus_4\n",
      "     |     |        |  |  |            \\-|\n",
      "     |     |        |  |  |              |   /-AGC24233.1|Sclerotinia_sclerotiorum_mitovirus_4\n",
      "     |     |        |  |  |               \\-|\n",
      "     |     |        |  |  |                  \\-ANJ77669.1|Sclerotinia_nivalis_mitovirus_1\n",
      "     |     |        |  |  |\n",
      "     |     |        |  |  |      /-YP_009259483.1|Cronartium_ribicola_mitovirus_4\n",
      "     |     |        |  |  |   /-|\n",
      "     |     |        |  |  |  |   \\-YP_009259482.1|Cronartium_ribicola_mitovirus_3\n",
      "     |     |        |  |   \\-|\n",
      "     |     |         \\-|     |   /-BAD72871.1|Helicobasidium_mompa_mitovirus_1_18\n",
      "     |     |           |      \\-|\n",
      "     |     |           |        |   /-YP_009259369.1|Cronartium_ribicola_mitovirus_1\n",
      "     |     |           |         \\-|\n",
      "     |     |           |            \\-YP_009259487.1|Cronartium_ribicola_mitovirus_5\n",
      "     |     |           |\n",
      "     |     |           |         /-YP_009259481.1|Cronartium_ribicola_mitovirus_2\n",
      "     |     |           |      /-|\n",
      "     |     |           |     |   \\-ALM62243.1|Soybean_leaf_associated_mitovirus_3\n",
      "     |     |           |     |\n",
      "     |     |           |     |         /-AHF48631.1|Sclerotinia_sclerotiorum_mitovirus_15\n",
      "     |     |           |     |      /-|\n",
      "     |     |           |     |     |  |   /-AAR01970.1|Cryphonectria_cubensis_mitovirus_1a\n",
      "     |     |           |   /-|     |   \\-|\n",
      "     |     |           |  |  |   /-|      \\-AMS38556.1|Hymenoscyphus_fraxineus_mitovirus_1\n",
      "     |     |           |  |  |  |  |\n",
      "     |     |           |  |  |  |  |   /-AHX84129.1|Sclerotinia_sclerotiorum_mitovirus_2\n",
      "     |     |           |  |  |  |   \\-|\n",
      "     |     |           |  |  |  |     |   /-AHF48625.1|Sclerotinia_sclerotiorum_mitovirus_9\n",
      "     |     |           |  |  |  |      \\-|\n",
      "     |     |           |  |   \\-|        |   /-AHF48632.1|Sclerotinia_sclerotiorum_mitovirus_16\n",
      "     |     |           |  |     |         \\-|\n",
      "     |     |           |  |     |           |   /-ALD89136.1|Sclerotinia_sclerotiorum_mitovirus_19\n",
      "     |     |            \\-|     |            \\-|\n",
      "     |     |              |     |               \\-AHX84135.1|Sclerotinia_sclerotiorum_mitovirus_7\n",
      "     |     |              |     |\n",
      "     |     |              |     |   /-YP_009336622.1|Shahe_narna_like_virus_7\n",
      "     |     |              |      \\-|\n",
      "     |     |              |        |   /-ALD89125.1|Rhizoctonia_solani_mitovirus_6\n",
      "     |     |              |         \\-|\n",
      "     |     |              |            \\-AQM32767.1|Agaricus_bisporus_mitovirus_1\n",
      "     |     |              |\n",
      "     |     |              |   /-AGW51760.1|Mitovirus_AEF_2013\n",
      "     |     |               \\-|\n",
      "     |     |                 |   /-AHF48622.1|Sclerotinia_sclerotiorum_mitovirus_6\n",
      "     |     |                  \\-|\n",
      "     |     |                    |   /-ALD89102.1|Rhizoctonia_solani_mitovirus_10\n",
      "     |     |                     \\-|\n",
      "     |     |                        \\-YP_009336494.1|Hubei_narna_like_virus_25\n",
      "     |     |\n",
      "     |      \\-BAU79521.1|Diatom_colony_associated_ssRNA_virus_2\n",
      "     |\n",
      "     |            /-YP_009337087.1|Hubei_narna_like_virus_4\n",
      "   /-|           |\n",
      "  |  |           |         /-YP_009336672.1|Hubei_narna_like_virus_5\n",
      "  |  |           |      /-|\n",
      "  |  |           |     |   \\-YP_009333214.1|narna_like_virus_6\n",
      "  |  |           |   /-|\n",
      "  |  |           |  |  |   /-YP_009333144.1|Beihai_narna_like_virus_3\n",
      "  |  |           |  |   \\-|\n",
      "  |  |         /-|  |     |   /-YP_009337264.1|Wenling_narna_like_virus_1\n",
      "  |  |        |  |  |      \\-|\n",
      "  |  |        |  |  |        |   /-YP_009337193.1|Wenling_narna_like_virus_2\n",
      "  |  |        |  |  |         \\-|\n",
      "  |  |        |  |  |           |   /-YP_009333266.1|Beihai_narna_like_virus_4\n",
      "  |  |        |  |  |            \\-|\n",
      "  |  |        |  |  |              |   /-YP_009336518.1|Wenzhou_narna_like_virus_2\n",
      "  |  |        |  |  |               \\-|\n",
      "  |  |        |  |  |                  \\-YP_009336520.1|Wenzhou_narna_like_virus_3\n",
      "  |  |        |   \\-|\n",
      "  |  |        |     |         /-YP_002019754.1|Epirus_cherry_virus\n",
      "  |  |        |     |      /-|\n",
      "  |  |        |     |     |  |   /-YP_002019757.1|Ourmia_melon_virus\n",
      "  |  |        |     |     |   \\-|\n",
      "  |  |        |     |   /-|      \\-YP_003104770.1|Cassava_virus_C\n",
      "  |  |        |     |  |  |\n",
      "  |  |        |     |  |  |   /-YP_009337787.1|Hubei_narna_like_virus_3\n",
      "  |  |        |     |  |   \\-|\n",
      "  |  |        |     |  |      \\-YP_009344979.1|Wuhan_spider_virus_7\n",
      "  |  |        |     |  |\n",
      "  |  |        |     |  |            /-YP_009333283.1|Beihai_narna_like_virus_1\n",
      "  |  |        |     |  |         /-|\n",
      "  |  |        |      \\-|        |   \\-YP_009336796.1|Hubei_narna_like_virus_2\n",
      "  |  |        |        |      /-|\n",
      "  |  |        |        |     |  |   /-YP_009336615.1|Shahe_narna_like_virus_2\n",
      "  |  |        |        |     |   \\-|\n",
      "  |  |        |        |     |     |   /-YP_009336536.1|Shahe_narna_like_virus_1\n",
      "  |  |        |        |   /-|      \\-|\n",
      "  |  |        |        |  |  |        |   /-YP_009337331.1|Changjiang_narna_like_virus_1\n",
      "  |  |        |        |  |  |         \\-|\n",
      "  |  |      /-|         \\-|  |            \\-APG77214.1|Hubei_narna_like_virus_1\n",
      "  |  |     |  |           |  |\n",
      "  |  |     |  |           |   \\-YP_009333277.1|Beihai_narna_like_virus_2\n",
      "  |  |     |  |           |\n",
      "  |  |     |  |            \\-APG77283.1|Wenzhou_narna_like_virus_1\n",
      "  |  |     |  |\n",
      "  |  |     |  |   /-YP_009337307.1|Wenzhou_narna_like_virus_5\n",
      "  |  |     |  |  |\n",
      "  |  |     |  |  |      /-YP_009336843.1|Hubei_narna_like_virus_7\n",
      "  |  |     |  |  |     |\n",
      "  |  |     |  |  |     |               /-ALD89138.1|Sclerotinia_sclerotiorum_ourmia_like_virus_1_RNA_1\n",
      "  |  |     |  |  |     |            /-|\n",
      "  |  |     |  |  |     |           |  |   /-ALM62238.1|Soybean_leaf_associated_ourmiavirus_1\n",
      "  |  |     |  |  |     |           |   \\-|\n",
      "  |  |     |  |  |     |         /-|      \\-ALM62250.1|Soybean_leaf_associated_ourmiavirus_2\n",
      "  |  |     |  |  |     |        |  |\n",
      "  |  |     |  |  |     |        |  |   /-ALD89131.1|Rhizoctonia_solani_ourmia_like_virus_1_RNA_1\n",
      "  |  |     |  |  |     |      /-|   \\-|\n",
      "  |  |     |  |  |     |     |  |      \\-SBQ28480.1|Magnaporthe_oryzae_ourmia_like_virus\n",
      "  |  |     |  |  |     |     |  |\n",
      "  |  |     |  |  |     |     |  |   /-YP_009345044.1|Phomopsis_longicolla_RNA_virus_1\n",
      "  |  |     |  |  |   /-|     |   \\-|\n",
      "  |  |     |  |  |  |  |   /-|     |   /-ALD89139.1|Sclerotinia_sclerotiorum_ourmia_like_virus_2_RNA_1\n",
      "  |  |     |  |  |  |  |  |  |      \\-|\n",
      "  |  |     |  |  |  |  |  |  |         \\-YP_009182165.1|Botrytis_ourmia_like_virus\n",
      "  |  |     |  |  |  |  |  |  |\n",
      "  |  |     |  |  |  |  |  |  |   /-YP_009333280.1|Beihai_narna_like_virus_20\n",
      "  |  |     |  |  |  |  |  |   \\-|\n",
      "  |  |     |  |  |  |  |  |     |   /-ALD89132.1|Rhizoctonia_solani_ourmia_like_virus_1_RNA_1\n",
      "  |  |     |  |  |  |  |  |      \\-|\n",
      "  |  |     |   \\-|  |  |  |        |   /-AQM49944.1|Agaricus_bisporus_virus_14\n",
      "  |  |     |     |  |  |  |         \\-|\n",
      "  |  |     |     |  |  |  |            \\-AQM49945.1|Agaricus_bisporus_virus_15\n",
      "  |  |     |     |  |  |  |\n",
      "  |  |     |     |  |  |  |      /-APG77281.1|Wenzhou_narna_like_virus_6\n",
      "  |  |     |     |  |   \\-|   /-|\n",
      "  |  |     |     |  |     |  |  |   /-YP_009336702.1|Sanxia_narna_like_virus_1\n",
      "  |  |     |     |  |     |  |   \\-|\n",
      "  |  |     |     |  |     |  |      \\-YP_009337343.1|Changjiang_narna_like_virus_3\n",
      "  |  |     |     |  |     |  |\n",
      "  |  |     |     |  |     |  |      /-YP_009336532.1|Shahe_narna_like_virus_3\n",
      "  |  |     |     |  |     |  |     |\n",
      "  |  |     |     |  |     |  |     |            /-YP_009333152.1|Beihai_narna_like_virus_16\n",
      "  |  |     |     |  |     |  |     |         /-|\n",
      "  |  |     |     |  |     |  |     |        |  |   /-YP_009337385.1|Changjiang_narna_like_virus_4\n",
      "  |  |     |     |  |     |  |     |        |   \\-|\n",
      "  |  |     |     |  |     |  |     |        |      \\-YP_009336759.1|Hubei_narna_like_virus_12\n",
      "  |  |     |     |  |     |  |     |      /-|\n",
      "  |  |     |     |  |     |  |     |     |  |      /-YP_009337619.1|Wenzhou_narna_like_virus_9\n",
      "  |  |     |     |  |      \\-|     |     |  |   /-|\n",
      "  |  |     |     |  |        |     |     |  |  |   \\-APG77303.1|Wenzhou_narna_like_virus_10\n",
      "  |  |     |     |  |        |   /-|     |  |  |\n",
      "  |  |     |     |  |        |  |  |     |   \\-|      /-YP_009333164.1|Beihai_narna_like_virus_19\n",
      "  |  |     |     |  |        |  |  |     |     |   /-|\n",
      "  |  |     |     |  |        |  |  |     |     |  |   \\-APG76981.1|Wenzhou_narna_like_virus_11\n",
      "  |  |   /-|      \\-|        |  |  |   /-|      \\-|\n",
      "  |  |  |  |        |        |  |  |  |  |        |   /-YP_009333347.1|Beihai_narna_like_virus_17\n",
      "  |  |  |  |        |        |  |  |  |  |         \\-|\n",
      "  |  |  |  |        |        |  |  |  |  |            \\-YP_009333315.1|Beihai_narna_like_virus_18\n",
      "  |  |  |  |        |        |  |  |  |  |\n",
      "  |  |  |  |        |        |  |  |  |  |      /-YP_009337436.1|Wenzhou_narna_like_virus_8\n",
      "  |  |  |  |        |        |  |  |  |  |   /-|\n",
      "  |  |  |  |        |        |  |  |  |  |  |  |   /-YP_009333147.1|Beihai_narna_like_virus_15\n",
      "  |  |  |  |        |        |  |   \\-|  |  |   \\-|\n",
      "  |  |  |  |        |         \\-|     |   \\-|      \\-APG77163.1|Wenzhou_shrimp_virus_10\n",
      "  |  |  |  |        |           |     |     |\n",
      "  |  |  |  |        |           |     |     |   /-YP_009333153.1|Beihai_narna_like_virus_14\n",
      "  |  |  |  |        |           |     |      \\-|\n",
      "  |  |  |  |        |           |     |         \\-YP_009333241.1|Beihai_narna_like_virus_13\n",
      "  |  |  |  |        |           |     |\n",
      "  |  |  |  |        |           |     |   /-YP_009337386.1|Wenzhou_narna_like_virus_7\n",
      "  |  |  |  |        |           |      \\-|\n",
      "  |  |  |  |        |           |         \\-YP_009333255.1|Beihai_narna_like_virus_12\n",
      "  |  |  |  |        |           |\n",
      "  |  |  |  |        |           |      /-YP_009337133.1|Wenling_narna_like_virus_4\n",
      "  |  |  |  |        |           |   /-|\n",
      "  |  |  |  |        |           |  |   \\-YP_009337146.1|Wenling_narna_like_virus_5\n",
      "  |  |  |  |        |            \\-|\n",
      "  |  |  |  |        |              |      /-YP_009337789.1|Hubei_mosquito_virus_3\n",
      "  |  |  |  |        |              |   /-|\n",
      "  |  |  |  |        |               \\-|   \\-APG77221.1|Hubei_narna_like_virus_14\n",
      "  |  |  |  |        |                 |\n",
      "  |  |  |  |        |                  \\-YP_009337805.1|Hubei_narna_like_virus_13\n",
      "  |  |  |  |        |\n",
      "  |  |  |  |        |   /-APG77202.1|Hubei_narna_like_virus_8\n",
      "  |  |  |  |        |  |\n",
      "  |  |  |  |        |  |      /-YP_009336802.1|Hubei_narna_like_virus_9\n",
      "  |  |  |  |        |  |   /-|\n",
      "  |  |  |  |        |  |  |  |   /-YP_009336851.1|Hubei_narna_like_virus_11\n",
      "  |  |  |  |         \\-|  |   \\-|\n",
      "  |  |  |  |           |  |     |   /-YP_009333271.1|Beihai_narna_like_virus_10\n",
      "  |  |  |  |           |  |      \\-|\n",
      "  |  |  |  |           |  |        |   /-YP_009333278.1|Beihai_narna_like_virus_11\n",
      "  |   \\-|  |           |  |         \\-|\n",
      "  |     |  |            \\-|            \\-YP_009333318.1|Wenzhou_narna_like_virus_4\n",
      "  |     |  |              |\n",
      "  |     |  |              |      /-YP_009333316.1|Beihai_narna_like_virus_9\n",
      "  |     |  |              |   /-|\n",
      "  |     |  |              |  |  |   /-YP_009333267.1|Beihai_narna_like_virus_8\n",
      "  |     |  |              |  |   \\-|\n",
      "  |     |  |              |  |     |   /-YP_009337159.1|Wenling_narna_like_virus_3\n",
      "  |     |  |              |  |      \\-|\n",
      "  |     |  |               \\-|        |   /-APG77081.1|Beihai_narna_like_virus_7\n",
      "  |     |  |                 |         \\-|\n",
      "  |     |  |                 |            \\-YP_009333251.1|Beihai_narna_like_virus_6\n",
      "  |     |  |                 |\n",
      "  |     |  |                 |   /-YP_009337311.1|Changjiang_narna_like_virus_2\n",
      "  |     |  |                  \\-|\n",
      "  |     |  |                     \\-APG77208.1|Hubei_narna_like_virus_10\n",
      "  |     |  |\n",
      "  |     |   \\-APG77011.1|Beihai_narna_like_virus_5\n",
      "  |     |\n",
      "  |     |            /-ANI86034.1|Leptomonas_seymouri_Narna_like_virus_1\n",
      "  |     |         /-|\n",
      "  |     |        |   \\-YP_009256557.1|Phytomonas_serpens_narnavirus_1\n",
      "  |     |      /-|\n",
      "  |     |     |  |   /-YP_009333140.1|Beihai_narna_like_virus_21\n",
      "  |     |     |   \\-|\n",
      "  |     |     |     |   /-YP_009241365.1|Phytophthora_infestans_RNA_virus_4\n",
      "  |     |     |      \\-|\n",
      "  |     |   /-|         \\-YP_009333179.1|Beihai_barnacle_virus_10\n",
      "  |     |  |  |\n",
      "  |     |  |  |      /-YP_009333261.1|Beihai_narna_like_virus_23\n",
      "  |     |  |  |   /-|\n",
      "  |     |  |  |  |   \\-YP_009333139.1|Beihai_narna_like_virus_22\n",
      "  |     |  |   \\-|\n",
      "  |     |  |     |   /-SCW25778.1|Beauveria_bassiana_narnavirus\n",
      "  |     |  |      \\-|\n",
      "  |     |  |         \\-APG77272.1|Wenling_narna_like_virus_6\n",
      "  |     |  |\n",
      "--|     |  |      /-NP_660177.1|Saccharomyces_23S_RNA_narnavirus\n",
      "  |     |  |     |\n",
      "  |      \\-|     |         /-YP_009342440.1|Wuhan_insect_virus_18\n",
      "  |        |     |      /-|\n",
      "  |        |   /-|     |  |   /-NP_660178.1|Saccharomyces_20S_RNA_narnavirus\n",
      "  |        |  |  |   /-|   \\-|\n",
      "  |        |  |  |  |  |      \\-APG77168.1|Shahe_narna_like_virus_4\n",
      "  |        |  |  |  |  |\n",
      "  |        |  |  |  |   \\-APG77120.1|Hubei_narna_like_virus_16\n",
      "  |        |  |   \\-|\n",
      "  |        |  |     |      /-YP_009342297.1|Wuchang_cockroach_Virus_4\n",
      "  |        |  |     |   /-|\n",
      "  |        |  |     |  |   \\-YP_009337783.1|Hubei_narna_like_virus_15\n",
      "  |        |  |      \\-|\n",
      "  |        |  |        |   /-YP_009337166.1|Wenling_narna_like_virus_7\n",
      "  |         \\-|         \\-|\n",
      "  |           |           |   /-APG77263.1|Wenling_narna_like_virus_8\n",
      "  |           |            \\-|\n",
      "  |           |               \\-YP_009333245.1|Beihai_narna_like_virus_24\n",
      "  |           |\n",
      "  |           |      /-YP_009272902.1|Fusarium_poae_narnavirus_1\n",
      "  |           |   /-|\n",
      "  |           |  |  |   /-APG76998.1|Beihai_narna_like_virus_25\n",
      "  |           |  |   \\-|\n",
      "  |           |  |      \\-APG77107.1|Changjiang_narna_like_virus_5\n",
      "  |           |  |\n",
      "  |           |  |         /-APG77159.1|Hubei_narna_like_virus_20\n",
      "  |            \\-|      /-|\n",
      "  |              |     |  |   /-YP_009337792.1|Hubei_narna_like_virus_19\n",
      "  |              |     |   \\-|\n",
      "  |              |   /-|      \\-YP_009330065.1|Hubei_narna_like_virus_18\n",
      "  |              |  |  |\n",
      "  |              |  |  |   /-YP_009342444.1|Wuhan_horsefly_Virus_3\n",
      "  |              |  |   \\-|\n",
      "  |               \\-|     |   /-YP_009336688.1|Sanxia_water_strider_virus_13\n",
      "  |                 |      \\-|\n",
      "  |                 |         \\-YP_009337744.1|Hubei_narna_like_virus_17\n",
      "  |                 |\n",
      "  |                 |   /-YP_009337775.1|Hubei_narna_like_virus_21\n",
      "  |                  \\-|\n",
      "  |                    |   /-AJT39597.1|Narnaviridae_environmental_sample\n",
      "  |                     \\-|\n",
      "  |                        \\-YP_009333331.1|Zhejiang_mosquito_virus_3\n",
      "  |\n",
      "  |            /-AMQ23551.1|unidentified_ssRNA_virus\n",
      "  |         /-|\n",
      "  |        |   \\-APG77058.1|Beihai_levi_like_virus_13\n",
      "  |        |\n",
      "  |        |      /-APG77079.1|Beihai_levi_like_virus_14\n",
      "  |        |   /-|\n",
      "  |        |  |  |   /-NP_085473.1|Acinetobacter_phage_AP205\n",
      "  |        |  |   \\-|\n",
      "  |        |  |     |   /-APG77171.1|Shahe_levi_like_virus_1\n",
      "  |        |  |      \\-|\n",
      "  |        |  |         \\-APG76993.1|Beihai_levi_like_virus_20\n",
      "  |        |  |\n",
      "  |        |  |                     /-YP_717671.1|Pseudomonas_phage_PRR1\n",
      "  |        |  |                  /-|\n",
      "  |        |  |                 |  |   /-APG77123.1|Hubei_levi_like_virus_6\n",
      "  |        |  |                 |   \\-|\n",
      "  |        |  |                 |     |   /-APG77018.1|Beihai_levi_like_virus_25\n",
      "  |        |  |                 |      \\-|\n",
      "  |        |  |               /-|         \\-YP_007111575.1|Enterobacteria_phage_M\n",
      "  |        |  |              |  |\n",
      "  |        |  |              |  |   /-YP_007237175.1|Enterobacteria_phage_Hgal1\n",
      "  |        |  |              |  |  |\n",
      "  |        |  |              |  |  |      /-APG77179.1|Shahe_levi_like_virus_2\n",
      "  |        |  |              |   \\-|   /-|\n",
      "  |        |  |              |     |  |  |   /-APG77225.1|Hubei_levi_like_virus_5\n",
      "  |        |  |              |     |  |   \\-|\n",
      "  |        |  |            /-|      \\-|      \\-YP_007237129.1|Enterobacteria_phage_C_1_INW_2012\n",
      "  |        |  |           |  |        |\n",
      "  |        |  |           |  |        |   /-APG77021.1|Beihai_levi_like_virus_22\n",
      "  |      /-|  |           |  |         \\-|\n",
      "  |     |  |  |           |  |           |   /-APG77037.1|Beihai_levi_like_virus_24\n",
      "  |     |  |  |           |  |            \\-|\n",
      "  |     |  |  |           |  |               \\-APG77003.1|Beihai_levi_like_virus_23\n",
      "  |     |  |  |         /-|  |\n",
      "  |     |  |  |        |  |  |   /-APG77136.1|Beihai_levi_like_virus_26\n",
      "  |     |  |  |        |  |   \\-|\n",
      "  |     |  |  |        |  |      \\-APG77025.1|Beihai_levi_like_virus_27\n",
      "  |     |  |  |        |  |\n",
      "  |     |  |  |        |  |      /-ACT66728.1|Escherichia_virus_BZ13\n",
      "  |     |  |  |        |  |   /-|\n",
      "  |     |  |  |        |  |  |   \\-ACT66724.1|Escherichia_virus_BZ13\n",
      "  |     |  |  |        |   \\-|\n",
      "  |     |  |  |        |     |   /-AFN39796.1|Enterobacterio_phage_MS2\n",
      "  |     |  |  |        |      \\-|\n",
      "  |     |  |  |        |        |   /-P15965.2|Enterobacteria_phage_fr\n",
      "  |     |  |  |        |         \\-|\n",
      "  |     |  |  |        |           |   /-ACN90651.1|Enterobacterio_phage_MS2\n",
      "  |     |  |  |        |            \\-|\n",
      "  |     |  |  |        |               \\-AEQ25571.1|Escherichia_virus_MS2\n",
      "  |     |  |  |        |\n",
      "  |     |  |  |        |            /-APG77106.1|Changjiang_levi_like_virus_2\n",
      "  |     |  |  |        |         /-|\n",
      "  |     |  |  |        |        |  |   /-APG77245.1|Hubei_levi_like_virus_3\n",
      "  |     |  |  |        |        |   \\-|\n",
      "  |     |  |  |        |      /-|      \\-APG77228.1|Hubei_levi_like_virus_2\n",
      "  |     |  |  |        |     |  |\n",
      "  |     |  |  |        |     |  |   /-APG77231.1|Hubei_levi_like_virus_7\n",
      "  |     |  |  |        |     |   \\-|\n",
      "  |     |   \\-|        |     |      \\-APG77234.1|Hubei_levi_like_virus_8\n",
      "  |     |     |        |     |\n",
      "  |     |     |        |     |            /-APG77195.1|Hubei_levi_like_virus_14\n",
      "  |     |     |      /-|     |         /-|\n",
      "  |     |     |     |  |     |        |   \\-APG77217.1|Hubei_levi_like_virus_13\n",
      "  |     |     |     |  |     |      /-|\n",
      "  |     |     |     |  |     |     |  |   /-APG76989.1|Beihai_levi_like_virus_34\n",
      "  |     |     |     |  |   /-|     |   \\-|\n",
      "  |     |     |     |  |  |  |     |      \\-APG77242.1|Hubei_levi_like_virus_10\n",
      "  |     |     |     |  |  |  |   /-|\n",
      "  |     |     |     |  |  |  |  |  |      /-APG77076.1|Beihai_levi_like_virus_35\n",
      "  |     |     |     |  |  |  |  |  |   /-|\n",
      "  |     |     |     |  |  |  |  |  |  |   \\-APG77174.1|Shahe_levi_like_virus_3\n",
      "  |     |     |     |  |  |  |  |   \\-|\n",
      "  |     |     |     |  |  |  |  |     |   /-APG77220.1|Hubei_levi_like_virus_12\n",
      "  |     |     |     |  |  |  |  |      \\-|\n",
      "  |     |     |     |  |  |  |  |        |   /-APG77271.1|Wenling_levi_like_virus_4\n",
      "  |     |     |     |  |  |  |  |         \\-|\n",
      "  |     |     |     |  |  |   \\-|            \\-APG77286.1|Wenzhou_levi_like_virus_7\n",
      "  |     |     |     |  |  |     |\n",
      "  |     |     |     |  |  |     |         /-APG77252.1|Wenling_levi_like_virus_2\n",
      "  |     |     |     |  |  |     |      /-|\n",
      "  |     |     |     |  |  |     |     |   \\-APG77292.1|Wenzhou_levi_like_virus_6\n",
      "  |     |     |     |  |  |     |   /-|\n",
      "  |   /-|     |     |  |  |     |  |  |   /-APG77255.1|Wenling_levi_like_virus_3\n",
      "  |  |  |     |     |  |  |     |  |   \\-|\n",
      "  |  |  |     |     |  |  |     |  |      \\-APG77127.1|Beihai_levi_like_virus_33\n",
      "  |  |  |     |     |  |  |     |  |\n",
      "  |  |  |     |     |  |  |      \\-|      /-APG77295.1|Wenzhou_levi_like_virus_4\n",
      "  |  |  |     |     |  |  |        |     |\n",
      "  |  |  |     |     |  |  |        |   /-|   /-APG77015.1|Beihai_levi_like_virus_32\n",
      "  |  |  |     |     |   \\-|        |  |  |  |\n",
      "  |  |  |     |     |     |        |  |   \\-|      /-APG77095.1|Beihai_levi_like_virus_29\n",
      "  |  |  |     |     |     |        |  |     |   /-|\n",
      "  |  |  |     |     |     |         \\-|      \\-|   \\-APG77071.1|Beihai_levi_like_virus_28\n",
      "  |  |  |     |     |     |           |        |\n",
      "  |  |  |     |   /-|     |           |         \\-APG77130.1|Beihai_levi_like_virus_30\n",
      "  |  |  |     |  |  |     |           |\n",
      "  |  |  |     |  |  |     |            \\-APG77237.1|Hubei_levi_like_virus_11\n",
      "  |  |  |     |  |  |     |\n",
      "  |  |  |     |  |  |     |            /-APG77119.1|Changjiang_levi_like_virus_3\n",
      "  |  |  |     |  |  |     |         /-|\n",
      "  |  |  |     |  |  |     |        |   \\-APG77248.1|Hubei_levi_like_virus_4\n",
      "  |  |  |     |  |  |     |      /-|\n",
      "  |  |  |     |  |  |     |     |  |   /-APG77067.1|Beihai_levi_like_virus_18\n",
      "  |  |  |     |  |  |     |     |  |  |\n",
      "  |  |  |     |  |  |     |     |   \\-|      /-APG77133.1|Beihai_levi_like_virus_15\n",
      "  |  |  |     |  |  |     |     |     |   /-|\n",
      "  |  |  |     |  |  |     |   /-|      \\-|   \\-APG77009.1|Beihai_levi_like_virus_16\n",
      "  |  |  |     |  |  |     |  |  |        |\n",
      "  |  |  |     |  |  |     |  |  |         \\-APG77031.1|Beihai_levi_like_virus_17\n",
      "  |  |  |     |  |  |     |  |  |\n",
      "  |  |  |     |  |  |     |  |  |   /-APG77201.1|Hubei_levi_like_virus_9\n",
      "  |  |  |     |  |  |     |  |   \\-|\n",
      "  |  |  |      \\-|  |     |  |      \\-APG76985.1|Beihai_levi_like_virus_19\n",
      "  |  |  |        |  |     |  |\n",
      "  |  |  |        |  |      \\-|         /-YP_009208148.1|Escherichia_virus_FI\n",
      "  |  |  |        |  |        |      /-|\n",
      "  |  |  |        |  |        |     |  |   /-ACY07245.1|Enterobacteria_phage_SP\n",
      "  |  |  |        |  |        |     |   \\-|\n",
      "  |  |  |        |  |        |   /-|      \\-ACT66752.1|Escherichia_virus_FI\n",
      "  |  |  |        |  |        |  |  |\n",
      "  |  |  |        |  |        |  |  |   /-4R71|Escherichia_virus_Qbeta\n",
      "  |  |  |        |  |        |  |   \\-|\n",
      "  |  |  |        |  |         \\-|     |   /-AAC06252.1|Enterobacteria_phage_M11\n",
      "  |  |  |        |  |           |      \\-|\n",
      "   \\-|  |        |  |           |         \\-NP_046752.1|Escherichia_virus_Qbeta\n",
      "     |  |        |  |           |\n",
      "     |  |        |  |            \\-APG77085.1|Beihai_levi_like_virus_21\n",
      "     |  |        |  |\n",
      "     |  |        |   \\-APG77299.1|Wenzhou_levi_like_virus_3\n",
      "     |  |        |\n",
      "     |  |        |   /-NP_042307.1|Pseudomonas_phage_PP7\n",
      "     |  |         \\-|\n",
      "     |  |            \\-SLL96484.1|LeviOr01_phage\n",
      "     |  |\n",
      "     |  |   /-AMQ23527.1|unidentified_ssRNA_virus\n",
      "     |   \\-|\n",
      "     |      \\-AMQ23519.1|unidentified_ssRNA_virus\n",
      "     |\n",
      "     |               /-APG77040.1|Beihai_levi_like_virus_3\n",
      "     |            /-|\n",
      "     |           |   \\-APG77043.1|Beihai_levi_like_virus_4\n",
      "     |         /-|\n",
      "     |        |  |   /-APG77049.1|Beihai_levi_like_virus_5\n",
      "     |        |   \\-|\n",
      "     |        |      \\-APG77064.1|Beihai_levi_like_virus_6\n",
      "     |        |\n",
      "     |        |         /-APG77046.1|Beihai_levi_like_virus_8\n",
      "     |      /-|      /-|\n",
      "     |     |  |     |   \\-APG77006.1|Beihai_levi_like_virus_7\n",
      "     |     |  |   /-|\n",
      "     |     |  |  |  |   /-APG77267.1|Wenling_levi_like_virus_1\n",
      "     |     |  |  |   \\-|\n",
      "     |     |  |  |      \\-APG77055.1|Beihai_levi_like_virus_9\n",
      "     |     |   \\-|\n",
      "     |     |     |      /-AHA86936.2|Marine_RNA_phage_MB\n",
      "     |     |     |   /-|\n",
      "     |   /-|     |  |   \\-YP_007002852.1|Caulobacter_phage_phiCb5\n",
      "     |  |  |      \\-|\n",
      "     |  |  |        |   /-APG77052.1|Beihai_levi_like_virus_10\n",
      "     |  |  |         \\-|\n",
      "     |  |  |            \\-APG77028.1|Beihai_levi_like_virus_11\n",
      "     |  |  |\n",
      "     |  |  |      /-APG77115.1|Changjiang_levi_like_virus_1\n",
      "     |  |  |   /-|\n",
      "     |  |  |  |  |   /-APG77280.1|Wenzhou_levi_like_virus_1\n",
      "     |  |  |  |   \\-|\n",
      "     |  |   \\-|      \\-APG77192.1|Sanxia_levi_like_virus_1\n",
      "      \\-|     |\n",
      "        |     |   /-APG77277.1|Wenzhou_levi_like_virus_2\n",
      "        |      \\-|\n",
      "        |        |   /-APG77239.1|Hubei_levi_like_virus_1\n",
      "        |         \\-|\n",
      "        |           |   /-APG77061.1|Beihai_levi_like_virus_1\n",
      "        |            \\-|\n",
      "        |               \\-APG77034.1|Beihai_levi_like_virus_2\n",
      "        |\n",
      "        |   /-AMQ23523.1|unidentified_ssRNA_virus\n",
      "        |  |\n",
      "        |  |      /-APG77142.1|Beihai_levi_like_virus_12\n",
      "         \\-|   /-|\n",
      "           |  |  |   /-AMQ23539.1|unidentified_ssRNA_virus\n",
      "           |  |   \\-|\n",
      "            \\-|      \\-AOG30806.1|ssRNA_phage_DC\n",
      "              |\n",
      "              |   /-AMQ23516.1|unidentified_ssRNA_virus\n",
      "               \\-|\n",
      "                 |   /-AMQ23536.1|unidentified_ssRNA_virus\n",
      "                  \\-|\n",
      "                     \\-AMQ23530.1|unidentified_ssRNA_virus\n"
     ]
    }
   ],
   "source": [
    "## Parsing and viewing trees\n",
    "\n",
    "## Koonin supplementary Data Set S2 = full tree of RdRps + RTs \n",
    "kooninTree = Tree(\"Data_Set_S2.nwk\")\n",
    "print(len(kooninTree))\n",
    "# 4627\n",
    "\n",
    "kooninNodes = kooninTree.get_leaves()\n",
    "# print(kooninNodes)\n",
    "\n",
    "## Koonin supplementary Data Set S4 = \n",
    "## representatives from Branch 1\n",
    "kooninBranch1 = Tree(\"Data_Set_S4_Branch_1.nwk\")\n",
    "# print(len(kooninBranch1))\n",
    "# 45\n",
    "\n",
    "branch1Leaves = kooninBranch1.get_leaf_names()\n",
    "\n",
    "## get the Branch 1 clade on the full tree\n",
    "branch1_clade = kooninTree.get_common_ancestor(branch1Leaves)\n",
    "# print(len(branch1_clade))\n",
    "# 299\n",
    "\n",
    "print(branch1_clade)"
   ]
  }
 ],
 "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.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}