notebooks/Step1_Taxonomic_analysis_extented_DFA_set.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "6338a08d",
   "metadata": {},
   "source": [
    "## Extented Divergent Fungal Actin (DFA) set - Taxonomic analysis - Kingdom level distribution\n",
    "\n",
    "**Notebook summary**\n",
    "\n",
    "This Jupyter Notebook:\n",
    "- Sorts taxonomic information of the species from the extented cluster into kingdom, phylum and order taxonomic ranks\n",
    "- Is associated with Step 1 of the Pub approach: Expanding the initial set of fungal species that possess DFA\n",
    "- Provides the code for the Pub Figures 4A\n",
    "\n",
    "\n",
    "**Context/Goal reminder**\n",
    "Divergent Fungal Actins (DFAs) have been identified in another Arcadia's project as potentially new actin isoform mainly found in Fungi. In this work, we have identified an extented set of 436 DFA sequences by combining the 2 new clusters identified in the DFA ProteinCartography run and the original DFA cluster. While the original cluster mainly contain fungal species, we want to characterize the kingdom-level distribution of the organisms that possess DFA in the extended cluster. \n",
    "\n",
    "**Notebook purpose**\n",
    "The purpose of this notebook is to investigate the extended set of DFA sequences and the taxonomy of the associated species. In this notebook, we:\n",
    "- characterize the different kingdom where DFA is found\n",
    "- visualize the distribution of DFA at the kingdom level"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "41cb4f96",
   "metadata": {},
   "source": [
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f7917312",
   "metadata": {},
   "source": [
    "### Setup path and environment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "8c792da8",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message:\n",
      "“package ‘ape’ was built under R version 4.2.2”\n",
      "\n",
      "Attaching package: ‘ape’\n",
      "\n",
      "\n",
      "The following object is masked from ‘package:ggtree’:\n",
      "\n",
      "    rotate\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "setwd('..')\n",
    "\n",
    "library(tidyverse)\n",
    "library(ggplot2)\n",
    "library(ggtree)\n",
    "library(ape)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d74b026b",
   "metadata": {},
   "source": [
    "### Function definition"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "21231778",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [],
   "source": [
    "# For each protein ID, the function extracts the associated organism information as well as its lineage information -\n",
    "# We further parse lineage information associated with 'Kingdom', 'phylum' ,'class' and 'order' when available\n",
    "\n",
    "\n",
    "lin_inf_str=function(data){\n",
    "  dat=data[,c('protid','Organism','Taxonomic.lineage')]\n",
    "  \n",
    "  data_taxo=data.frame()\n",
    "  \n",
    "  for (i in 1:dim(dat)[1]){\n",
    "    line=dat[i,'Taxonomic.lineage']\n",
    "    vec_lin=c(strsplit(line, \",\")[[1]])\n",
    "    \n",
    "    #extraction of kingdom info\n",
    "    loc_king=grep('(kingdom)',vec_lin,fixed=TRUE)\n",
    "    \n",
    "    king_inf=vec_lin[loc_king]\n",
    "    \n",
    "    if (length(loc_king)>=1){\n",
    "      vec_kin=c(strsplit(king_inf, \" \")[[1]])\n",
    "      \n",
    "      kingdom=vec_kin[2]\n",
    "      \n",
    "    }\n",
    "    \n",
    "    if (length(loc_king)==0){\n",
    "      kingdom='no_kingdom_rank'\n",
    "    }\n",
    "    \n",
    "    #extraction of phylum info\n",
    "    loc_phy=grep('(phylum)',vec_lin,fixed=TRUE)\n",
    "    \n",
    "    phy_inf=vec_lin[loc_phy]\n",
    "    \n",
    "    if (length(loc_phy)>=1){\n",
    "      vec_phy=c(strsplit(phy_inf, \" \")[[1]])\n",
    "      \n",
    "      phylum=vec_phy[2]\n",
    "      \n",
    "    }\n",
    "    \n",
    "    if (length(loc_phy)==0){\n",
    "      phylum='no_phylum_rank'\n",
    "    }\n",
    "    \n",
    "    #extraction of order info\n",
    "    loc_ord=grep('(order)',vec_lin,fixed=TRUE)\n",
    "    \n",
    "    ord_inf=vec_lin[loc_ord]\n",
    "    \n",
    "    if (length(loc_ord)>=1){\n",
    "      vec_ord=c(strsplit(ord_inf, \" \")[[1]])\n",
    "      \n",
    "      orderO=vec_ord[2]\n",
    "      \n",
    "    }\n",
    "    \n",
    "    if (length(loc_ord)==0){\n",
    "      orderO='no_order_rank'\n",
    "    }\n",
    "    \n",
    "    dat_temp=data.frame(dat[i,],\n",
    "                        'Kingdom'=kingdom,'Phylum'=phylum,'Order'=orderO)\n",
    "    \n",
    "    data_taxo=rbind(data_taxo,dat_temp)\n",
    "\n",
    "  }\n",
    "  \n",
    "  data_taxo\n",
    "  \n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7db4e869",
   "metadata": {},
   "source": [
    "### Analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "b7635560",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [],
   "source": [
    "## Data import\n",
    "\n",
    "data_ext=read.csv('data/step1/Extended_cluster_hits_DFA.csv')  # file with extended cluster information (resulting from the combnation of the new ProteinCartography DFA clusters and original DFA cluster)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "a7f30740",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table class=\"dataframe\">\n",
       "<caption>A data.frame: 10 × 6</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>protid</th><th scope=col>Organism</th><th scope=col>Taxonomic.lineage</th><th scope=col>Kingdom</th><th scope=col>Phylum</th><th scope=col>Order</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>A0A015ICG5</td><td>Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)                                      </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                        </td><td>Fungi</td><td>Mucoromycota   </td><td>Glomerales     </td></tr>\n",
       "\t<tr><th scope=row>2</th><td>A0A0B7FRH7</td><td>Thanatephorus cucumeris (strain AG1-IB / isolate 7/3/14) (Lettuce bottom rot fungus) (Rhizoctonia solani)</td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Ceratobasidiaceae (family), Rhizoctonia (genus), Rhizoctonia solani (species), Rhizoctonia solani AG-1 (no rank)</td><td>Fungi</td><td>Basidiomycota  </td><td>Cantharellales </td></tr>\n",
       "\t<tr><th scope=row>3</th><td>A0A1E3I214</td><td>Cryptococcus depauperatus CBS 7841                                                                       </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Tremellomycetes (class), Tremellales (order), Cryptococcaceae (family), Cryptococcus (genus), Cryptococcus depauperatus (species)                                                                        </td><td>Fungi</td><td>Basidiomycota  </td><td>Tremellales    </td></tr>\n",
       "\t<tr><th scope=row>4</th><td>A0A4P9W5S9</td><td>Blyttiomyces helicus                                                                                     </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiomycetes incertae sedis (no rank), Blyttiomyces (genus)                                                                                       </td><td>Fungi</td><td>Chytridiomycota</td><td>no_order_rank  </td></tr>\n",
       "\t<tr><th scope=row>5</th><td>C6H2X6    </td><td>Ajellomyces capsulatus (strain H143) (Darling's disease fungus) (Histoplasma capsulatum)                 </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Eurotiomycetidae (subclass), Onygenales (order), Ajellomycetaceae (family), Histoplasma (genus), Ajellomyces capsulatus (species)  </td><td>Fungi</td><td>Ascomycota     </td><td>Onygenales     </td></tr>\n",
       "\t<tr><th scope=row>6</th><td>A0A507EZQ4</td><td>Chytriomyces confervae                                                                                   </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiales (order), Chytriomycetaceae (family), Chytriomyces (genus)                                                                                </td><td>Fungi</td><td>Chytridiomycota</td><td>Chytridiales   </td></tr>\n",
       "\t<tr><th scope=row>7</th><td>A0A1M8A8G9</td><td>Malassezia sympodialis (strain ATCC 42132) (Atopic eczema-associated yeast)                              </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Malasseziomycetes (class), Malasseziales (order), Malasseziaceae (family), Malassezia (genus), Malassezia sympodialis (species)                                                                       </td><td>Fungi</td><td>Basidiomycota  </td><td>Malasseziales  </td></tr>\n",
       "\t<tr><th scope=row>8</th><td>U5HBM6    </td><td>Microbotryum lychnidis-dioicae (strain p1A1 Lamole / MvSl-1064) (Anther smut fungus)                     </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Microbotryales (order), Microbotryaceae (family), Microbotryum (genus), Microbotryum lychnidis-dioicae (species)                                                            </td><td>Fungi</td><td>Basidiomycota  </td><td>Microbotryales </td></tr>\n",
       "\t<tr><th scope=row>9</th><td>A0A8A3P372</td><td>Monilinia vaccinii-corymbosi                                                                             </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), sordariomyceta (no rank), Leotiomycetes (class), Helotiales (order), Sclerotiniaceae (family), Monilinia (genus)                                           </td><td>Fungi</td><td>Ascomycota     </td><td>Helotiales     </td></tr>\n",
       "\t<tr><th scope=row>10</th><td>A0A2T0AG86</td><td>Rhodosporidium toruloides (Yeast) (Rhodotorula gracilis)                                                 </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Sporidiobolales (order), Sporidiobolaceae (family), Rhodotorula (genus)                                                                                                     </td><td>Fungi</td><td>Basidiomycota  </td><td>Sporidiobolales</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 10 × 6\n",
       "\\begin{tabular}{r|llllll}\n",
       "  & protid & Organism & Taxonomic.lineage & Kingdom & Phylum & Order\\\\\n",
       "  & <chr> & <chr> & <chr> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & A0A015ICG5 & Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)                                       & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                         & Fungi & Mucoromycota    & Glomerales     \\\\\n",
       "\t2 & A0A0B7FRH7 & Thanatephorus cucumeris (strain AG1-IB / isolate 7/3/14) (Lettuce bottom rot fungus) (Rhizoctonia solani) & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Ceratobasidiaceae (family), Rhizoctonia (genus), Rhizoctonia solani (species), Rhizoctonia solani AG-1 (no rank) & Fungi & Basidiomycota   & Cantharellales \\\\\n",
       "\t3 & A0A1E3I214 & Cryptococcus depauperatus CBS 7841                                                                        & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Tremellomycetes (class), Tremellales (order), Cryptococcaceae (family), Cryptococcus (genus), Cryptococcus depauperatus (species)                                                                         & Fungi & Basidiomycota   & Tremellales    \\\\\n",
       "\t4 & A0A4P9W5S9 & Blyttiomyces helicus                                                                                      & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiomycetes incertae sedis (no rank), Blyttiomyces (genus)                                                                                        & Fungi & Chytridiomycota & no\\_order\\_rank  \\\\\n",
       "\t5 & C6H2X6     & Ajellomyces capsulatus (strain H143) (Darling's disease fungus) (Histoplasma capsulatum)                  & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Eurotiomycetidae (subclass), Onygenales (order), Ajellomycetaceae (family), Histoplasma (genus), Ajellomyces capsulatus (species)   & Fungi & Ascomycota      & Onygenales     \\\\\n",
       "\t6 & A0A507EZQ4 & Chytriomyces confervae                                                                                    & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiales (order), Chytriomycetaceae (family), Chytriomyces (genus)                                                                                 & Fungi & Chytridiomycota & Chytridiales   \\\\\n",
       "\t7 & A0A1M8A8G9 & Malassezia sympodialis (strain ATCC 42132) (Atopic eczema-associated yeast)                               & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Malasseziomycetes (class), Malasseziales (order), Malasseziaceae (family), Malassezia (genus), Malassezia sympodialis (species)                                                                        & Fungi & Basidiomycota   & Malasseziales  \\\\\n",
       "\t8 & U5HBM6     & Microbotryum lychnidis-dioicae (strain p1A1 Lamole / MvSl-1064) (Anther smut fungus)                      & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Microbotryales (order), Microbotryaceae (family), Microbotryum (genus), Microbotryum lychnidis-dioicae (species)                                                             & Fungi & Basidiomycota   & Microbotryales \\\\\n",
       "\t9 & A0A8A3P372 & Monilinia vaccinii-corymbosi                                                                              & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), sordariomyceta (no rank), Leotiomycetes (class), Helotiales (order), Sclerotiniaceae (family), Monilinia (genus)                                            & Fungi & Ascomycota      & Helotiales     \\\\\n",
       "\t10 & A0A2T0AG86 & Rhodosporidium toruloides (Yeast) (Rhodotorula gracilis)                                                  & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Sporidiobolales (order), Sporidiobolaceae (family), Rhodotorula (genus)                                                                                                      & Fungi & Basidiomycota   & Sporidiobolales\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 10 × 6\n",
       "\n",
       "| <!--/--> | protid &lt;chr&gt; | Organism &lt;chr&gt; | Taxonomic.lineage &lt;chr&gt; | Kingdom &lt;chr&gt; | Phylum &lt;chr&gt; | Order &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|\n",
       "| 1 | A0A015ICG5 | Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)                                       | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                         | Fungi | Mucoromycota    | Glomerales      |\n",
       "| 2 | A0A0B7FRH7 | Thanatephorus cucumeris (strain AG1-IB / isolate 7/3/14) (Lettuce bottom rot fungus) (Rhizoctonia solani) | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Ceratobasidiaceae (family), Rhizoctonia (genus), Rhizoctonia solani (species), Rhizoctonia solani AG-1 (no rank) | Fungi | Basidiomycota   | Cantharellales  |\n",
       "| 3 | A0A1E3I214 | Cryptococcus depauperatus CBS 7841                                                                        | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Tremellomycetes (class), Tremellales (order), Cryptococcaceae (family), Cryptococcus (genus), Cryptococcus depauperatus (species)                                                                         | Fungi | Basidiomycota   | Tremellales     |\n",
       "| 4 | A0A4P9W5S9 | Blyttiomyces helicus                                                                                      | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiomycetes incertae sedis (no rank), Blyttiomyces (genus)                                                                                        | Fungi | Chytridiomycota | no_order_rank   |\n",
       "| 5 | C6H2X6     | Ajellomyces capsulatus (strain H143) (Darling's disease fungus) (Histoplasma capsulatum)                  | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Eurotiomycetidae (subclass), Onygenales (order), Ajellomycetaceae (family), Histoplasma (genus), Ajellomyces capsulatus (species)   | Fungi | Ascomycota      | Onygenales      |\n",
       "| 6 | A0A507EZQ4 | Chytriomyces confervae                                                                                    | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiales (order), Chytriomycetaceae (family), Chytriomyces (genus)                                                                                 | Fungi | Chytridiomycota | Chytridiales    |\n",
       "| 7 | A0A1M8A8G9 | Malassezia sympodialis (strain ATCC 42132) (Atopic eczema-associated yeast)                               | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Malasseziomycetes (class), Malasseziales (order), Malasseziaceae (family), Malassezia (genus), Malassezia sympodialis (species)                                                                        | Fungi | Basidiomycota   | Malasseziales   |\n",
       "| 8 | U5HBM6     | Microbotryum lychnidis-dioicae (strain p1A1 Lamole / MvSl-1064) (Anther smut fungus)                      | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Microbotryales (order), Microbotryaceae (family), Microbotryum (genus), Microbotryum lychnidis-dioicae (species)                                                             | Fungi | Basidiomycota   | Microbotryales  |\n",
       "| 9 | A0A8A3P372 | Monilinia vaccinii-corymbosi                                                                              | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), sordariomyceta (no rank), Leotiomycetes (class), Helotiales (order), Sclerotiniaceae (family), Monilinia (genus)                                            | Fungi | Ascomycota      | Helotiales      |\n",
       "| 10 | A0A2T0AG86 | Rhodosporidium toruloides (Yeast) (Rhodotorula gracilis)                                                  | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Sporidiobolales (order), Sporidiobolaceae (family), Rhodotorula (genus)                                                                                                      | Fungi | Basidiomycota   | Sporidiobolales |\n",
       "\n"
      ],
      "text/plain": [
       "   protid    \n",
       "1  A0A015ICG5\n",
       "2  A0A0B7FRH7\n",
       "3  A0A1E3I214\n",
       "4  A0A4P9W5S9\n",
       "5  C6H2X6    \n",
       "6  A0A507EZQ4\n",
       "7  A0A1M8A8G9\n",
       "8  U5HBM6    \n",
       "9  A0A8A3P372\n",
       "10 A0A2T0AG86\n",
       "   Organism                                                                                                 \n",
       "1  Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)                                      \n",
       "2  Thanatephorus cucumeris (strain AG1-IB / isolate 7/3/14) (Lettuce bottom rot fungus) (Rhizoctonia solani)\n",
       "3  Cryptococcus depauperatus CBS 7841                                                                       \n",
       "4  Blyttiomyces helicus                                                                                     \n",
       "5  Ajellomyces capsulatus (strain H143) (Darling's disease fungus) (Histoplasma capsulatum)                 \n",
       "6  Chytriomyces confervae                                                                                   \n",
       "7  Malassezia sympodialis (strain ATCC 42132) (Atopic eczema-associated yeast)                              \n",
       "8  Microbotryum lychnidis-dioicae (strain p1A1 Lamole / MvSl-1064) (Anther smut fungus)                     \n",
       "9  Monilinia vaccinii-corymbosi                                                                             \n",
       "10 Rhodosporidium toruloides (Yeast) (Rhodotorula gracilis)                                                 \n",
       "   Taxonomic.lineage                                                                                                                                                                                                                                                                                                                                                                    \n",
       "1  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                        \n",
       "2  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Ceratobasidiaceae (family), Rhizoctonia (genus), Rhizoctonia solani (species), Rhizoctonia solani AG-1 (no rank)\n",
       "3  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Tremellomycetes (class), Tremellales (order), Cryptococcaceae (family), Cryptococcus (genus), Cryptococcus depauperatus (species)                                                                        \n",
       "4  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiomycetes incertae sedis (no rank), Blyttiomyces (genus)                                                                                       \n",
       "5  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Eurotiomycetidae (subclass), Onygenales (order), Ajellomycetaceae (family), Histoplasma (genus), Ajellomyces capsulatus (species)  \n",
       "6  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Chytridiomycota (phylum), Chytridiomycota incertae sedis (no rank), Chytridiomycetes (class), Chytridiales (order), Chytriomycetaceae (family), Chytriomyces (genus)                                                                                \n",
       "7  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Malasseziomycetes (class), Malasseziales (order), Malasseziaceae (family), Malassezia (genus), Malassezia sympodialis (species)                                                                       \n",
       "8  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Microbotryales (order), Microbotryaceae (family), Microbotryum (genus), Microbotryum lychnidis-dioicae (species)                                                            \n",
       "9  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), sordariomyceta (no rank), Leotiomycetes (class), Helotiales (order), Sclerotiniaceae (family), Monilinia (genus)                                           \n",
       "10 cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Pucciniomycotina (subphylum), Microbotryomycetes (class), Sporidiobolales (order), Sporidiobolaceae (family), Rhodotorula (genus)                                                                                                     \n",
       "   Kingdom Phylum          Order          \n",
       "1  Fungi   Mucoromycota    Glomerales     \n",
       "2  Fungi   Basidiomycota   Cantharellales \n",
       "3  Fungi   Basidiomycota   Tremellales    \n",
       "4  Fungi   Chytridiomycota no_order_rank  \n",
       "5  Fungi   Ascomycota      Onygenales     \n",
       "6  Fungi   Chytridiomycota Chytridiales   \n",
       "7  Fungi   Basidiomycota   Malasseziales  \n",
       "8  Fungi   Basidiomycota   Microbotryales \n",
       "9  Fungi   Ascomycota      Helotiales     \n",
       "10 Fungi   Basidiomycota   Sporidiobolales"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "## Sorting taxonomic information into kingdom/phylum, and order\n",
    "DFA_clust_taxo=lin_inf_str(data_ext)\n",
    "head(DFA_clust_taxo,10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "81a103ee",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table class=\"dataframe\">\n",
       "<caption>A data.frame: 10 × 6</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>protid</th><th scope=col>Organism</th><th scope=col>Taxonomic.lineage</th><th scope=col>Kingdom</th><th scope=col>Phylum</th><th scope=col>Order</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>A0A015ICG5</td><td>Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)   </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                               </td><td>Fungi</td><td>Mucoromycota </td><td>Glomerales        </td></tr>\n",
       "\t<tr><th scope=row>2</th><td>A0A060SBE2</td><td>Pycnoporus cinnabarinus (Cinnabar-red polypore) (Trametes cinnabarina)</td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Polyporales (order), Polyporaceae (family), Trametes (genus)                                                                                   </td><td>Fungi</td><td>Basidiomycota</td><td>Polyporales       </td></tr>\n",
       "\t<tr><th scope=row>3</th><td>A0A061AJ95</td><td>Cyberlindnera fabianii (Yeast) (Hansenula fabianii)                   </td><td>Cyberlindnera (genus), Phaffomycetaceae (family), Saccharomycetales (order), Saccharomycetes (class), Saccharomycotina (subphylum), saccharomyceta (no rank), Ascomycota (phylum), Dikarya (subkingdom), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                                    </td><td>Fungi</td><td>Ascomycota   </td><td>Saccharomycetales </td></tr>\n",
       "\t<tr><th scope=row>4</th><td>A0A061HA90</td><td>Pseudozyma flocculosa PF-1                                            </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Ustilaginomycetes (class), Ustilaginales (order), Ustilaginaceae (family), Pseudozyma (genus), Pseudozyma flocculosa (species)                                                                               </td><td>Fungi</td><td>Basidiomycota</td><td>Ustilaginales     </td></tr>\n",
       "\t<tr><th scope=row>5</th><td>A0A066VBJ4</td><td>Tilletiaria anomala (strain ATCC 24038 / CBS 436.72 / UBC 951)        </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Exobasidiomycetes (class), Georgefischeriales (order), Tilletiariaceae (family), Tilletiaria (genus), Tilletiaria anomala (species)                                                                          </td><td>Fungi</td><td>Basidiomycota</td><td>Georgefischeriales</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>A0A067LYG7</td><td>Botryobasidium botryosum FD-172 SS1                                   </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Botryobasidiaceae (family), Botryobasidium (genus), Botryobasidium botryosum (species)                                 </td><td>Fungi</td><td>Basidiomycota</td><td>Cantharellales    </td></tr>\n",
       "\t<tr><th scope=row>7</th><td>A0A067QC48</td><td>Jaapia argillacea MUCL 33604                                          </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Jaapiales (order), Jaapiaceae (family), Jaapia (genus), Jaapia argillacea (species)                                                                        </td><td>Fungi</td><td>Basidiomycota</td><td>Jaapiales         </td></tr>\n",
       "\t<tr><th scope=row>8</th><td>A0A067TP80</td><td>Galerina marginata (strain CBS 339.88)                                </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Agaricales (order), Agaricineae (suborder), Strophariaceae (family), Galerina (genus), Galerina marginata (species)                                        </td><td>Fungi</td><td>Basidiomycota</td><td>Agaricales        </td></tr>\n",
       "\t<tr><th scope=row>9</th><td>A0A068RS63</td><td>Lichtheimia corymbifera JMRC:FSU:9682                                 </td><td>Lichtheimia corymbifera (species), Lichtheimia (genus), Lichtheimiaceae (family), Mucorales (order), Mucoromycetes (class), Mucoromycotina (subphylum), Mucoromycota (phylum), Fungi incertae sedis (no rank), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                              </td><td>Fungi</td><td>Mucoromycota </td><td>Mucorales         </td></tr>\n",
       "\t<tr><th scope=row>10</th><td>A0A072PQI0</td><td>Exophiala aquamarina CBS 119918                                       </td><td>cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Chaetothyriomycetidae (subclass), Chaetothyriales (order), Herpotrichiellaceae (family), Exophiala (genus), Exophiala aquamarina (species)</td><td>Fungi</td><td>Ascomycota   </td><td>Chaetothyriales   </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 10 × 6\n",
       "\\begin{tabular}{r|llllll}\n",
       "  & protid & Organism & Taxonomic.lineage & Kingdom & Phylum & Order\\\\\n",
       "  & <chr> & <chr> & <chr> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & A0A015ICG5 & Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)    & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                                & Fungi & Mucoromycota  & Glomerales        \\\\\n",
       "\t2 & A0A060SBE2 & Pycnoporus cinnabarinus (Cinnabar-red polypore) (Trametes cinnabarina) & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Polyporales (order), Polyporaceae (family), Trametes (genus)                                                                                    & Fungi & Basidiomycota & Polyporales       \\\\\n",
       "\t3 & A0A061AJ95 & Cyberlindnera fabianii (Yeast) (Hansenula fabianii)                    & Cyberlindnera (genus), Phaffomycetaceae (family), Saccharomycetales (order), Saccharomycetes (class), Saccharomycotina (subphylum), saccharomyceta (no rank), Ascomycota (phylum), Dikarya (subkingdom), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                                     & Fungi & Ascomycota    & Saccharomycetales \\\\\n",
       "\t4 & A0A061HA90 & Pseudozyma flocculosa PF-1                                             & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Ustilaginomycetes (class), Ustilaginales (order), Ustilaginaceae (family), Pseudozyma (genus), Pseudozyma flocculosa (species)                                                                                & Fungi & Basidiomycota & Ustilaginales     \\\\\n",
       "\t5 & A0A066VBJ4 & Tilletiaria anomala (strain ATCC 24038 / CBS 436.72 / UBC 951)         & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Exobasidiomycetes (class), Georgefischeriales (order), Tilletiariaceae (family), Tilletiaria (genus), Tilletiaria anomala (species)                                                                           & Fungi & Basidiomycota & Georgefischeriales\\\\\n",
       "\t6 & A0A067LYG7 & Botryobasidium botryosum FD-172 SS1                                    & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Botryobasidiaceae (family), Botryobasidium (genus), Botryobasidium botryosum (species)                                  & Fungi & Basidiomycota & Cantharellales    \\\\\n",
       "\t7 & A0A067QC48 & Jaapia argillacea MUCL 33604                                           & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Jaapiales (order), Jaapiaceae (family), Jaapia (genus), Jaapia argillacea (species)                                                                         & Fungi & Basidiomycota & Jaapiales         \\\\\n",
       "\t8 & A0A067TP80 & Galerina marginata (strain CBS 339.88)                                 & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Agaricales (order), Agaricineae (suborder), Strophariaceae (family), Galerina (genus), Galerina marginata (species)                                         & Fungi & Basidiomycota & Agaricales        \\\\\n",
       "\t9 & A0A068RS63 & Lichtheimia corymbifera JMRC:FSU:9682                                  & Lichtheimia corymbifera (species), Lichtheimia (genus), Lichtheimiaceae (family), Mucorales (order), Mucoromycetes (class), Mucoromycotina (subphylum), Mucoromycota (phylum), Fungi incertae sedis (no rank), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                               & Fungi & Mucoromycota  & Mucorales         \\\\\n",
       "\t10 & A0A072PQI0 & Exophiala aquamarina CBS 119918                                        & cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Chaetothyriomycetidae (subclass), Chaetothyriales (order), Herpotrichiellaceae (family), Exophiala (genus), Exophiala aquamarina (species) & Fungi & Ascomycota    & Chaetothyriales   \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 10 × 6\n",
       "\n",
       "| <!--/--> | protid &lt;chr&gt; | Organism &lt;chr&gt; | Taxonomic.lineage &lt;chr&gt; | Kingdom &lt;chr&gt; | Phylum &lt;chr&gt; | Order &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|\n",
       "| 1 | A0A015ICG5 | Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)    | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                                | Fungi | Mucoromycota  | Glomerales         |\n",
       "| 2 | A0A060SBE2 | Pycnoporus cinnabarinus (Cinnabar-red polypore) (Trametes cinnabarina) | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Polyporales (order), Polyporaceae (family), Trametes (genus)                                                                                    | Fungi | Basidiomycota | Polyporales        |\n",
       "| 3 | A0A061AJ95 | Cyberlindnera fabianii (Yeast) (Hansenula fabianii)                    | Cyberlindnera (genus), Phaffomycetaceae (family), Saccharomycetales (order), Saccharomycetes (class), Saccharomycotina (subphylum), saccharomyceta (no rank), Ascomycota (phylum), Dikarya (subkingdom), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                                     | Fungi | Ascomycota    | Saccharomycetales  |\n",
       "| 4 | A0A061HA90 | Pseudozyma flocculosa PF-1                                             | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Ustilaginomycetes (class), Ustilaginales (order), Ustilaginaceae (family), Pseudozyma (genus), Pseudozyma flocculosa (species)                                                                                | Fungi | Basidiomycota | Ustilaginales      |\n",
       "| 5 | A0A066VBJ4 | Tilletiaria anomala (strain ATCC 24038 / CBS 436.72 / UBC 951)         | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Exobasidiomycetes (class), Georgefischeriales (order), Tilletiariaceae (family), Tilletiaria (genus), Tilletiaria anomala (species)                                                                           | Fungi | Basidiomycota | Georgefischeriales |\n",
       "| 6 | A0A067LYG7 | Botryobasidium botryosum FD-172 SS1                                    | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Botryobasidiaceae (family), Botryobasidium (genus), Botryobasidium botryosum (species)                                  | Fungi | Basidiomycota | Cantharellales     |\n",
       "| 7 | A0A067QC48 | Jaapia argillacea MUCL 33604                                           | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Jaapiales (order), Jaapiaceae (family), Jaapia (genus), Jaapia argillacea (species)                                                                         | Fungi | Basidiomycota | Jaapiales          |\n",
       "| 8 | A0A067TP80 | Galerina marginata (strain CBS 339.88)                                 | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Agaricales (order), Agaricineae (suborder), Strophariaceae (family), Galerina (genus), Galerina marginata (species)                                         | Fungi | Basidiomycota | Agaricales         |\n",
       "| 9 | A0A068RS63 | Lichtheimia corymbifera JMRC:FSU:9682                                  | Lichtheimia corymbifera (species), Lichtheimia (genus), Lichtheimiaceae (family), Mucorales (order), Mucoromycetes (class), Mucoromycotina (subphylum), Mucoromycota (phylum), Fungi incertae sedis (no rank), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                               | Fungi | Mucoromycota  | Mucorales          |\n",
       "| 10 | A0A072PQI0 | Exophiala aquamarina CBS 119918                                        | cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Chaetothyriomycetidae (subclass), Chaetothyriales (order), Herpotrichiellaceae (family), Exophiala (genus), Exophiala aquamarina (species) | Fungi | Ascomycota    | Chaetothyriales    |\n",
       "\n"
      ],
      "text/plain": [
       "   protid    \n",
       "1  A0A015ICG5\n",
       "2  A0A060SBE2\n",
       "3  A0A061AJ95\n",
       "4  A0A061HA90\n",
       "5  A0A066VBJ4\n",
       "6  A0A067LYG7\n",
       "7  A0A067QC48\n",
       "8  A0A067TP80\n",
       "9  A0A068RS63\n",
       "10 A0A072PQI0\n",
       "   Organism                                                              \n",
       "1  Rhizophagus irregularis (strain DAOM 197198w) (Glomus intraradices)   \n",
       "2  Pycnoporus cinnabarinus (Cinnabar-red polypore) (Trametes cinnabarina)\n",
       "3  Cyberlindnera fabianii (Yeast) (Hansenula fabianii)                   \n",
       "4  Pseudozyma flocculosa PF-1                                            \n",
       "5  Tilletiaria anomala (strain ATCC 24038 / CBS 436.72 / UBC 951)        \n",
       "6  Botryobasidium botryosum FD-172 SS1                                   \n",
       "7  Jaapia argillacea MUCL 33604                                          \n",
       "8  Galerina marginata (strain CBS 339.88)                                \n",
       "9  Lichtheimia corymbifera JMRC:FSU:9682                                 \n",
       "10 Exophiala aquamarina CBS 119918                                       \n",
       "   Taxonomic.lineage                                                                                                                                                                                                                                                                                                                                                                           \n",
       "1  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Fungi incertae sedis (no rank), Mucoromycota (phylum), Glomeromycotina (subphylum), Glomeromycetes (class), Glomerales (order), Glomeraceae (family), Rhizophagus (genus), Rhizophagus irregularis (species)                                                                               \n",
       "2  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Polyporales (order), Polyporaceae (family), Trametes (genus)                                                                                   \n",
       "3  Cyberlindnera (genus), Phaffomycetaceae (family), Saccharomycetales (order), Saccharomycetes (class), Saccharomycotina (subphylum), saccharomyceta (no rank), Ascomycota (phylum), Dikarya (subkingdom), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                                    \n",
       "4  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Ustilaginomycetes (class), Ustilaginales (order), Ustilaginaceae (family), Pseudozyma (genus), Pseudozyma flocculosa (species)                                                                               \n",
       "5  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Ustilaginomycotina (subphylum), Exobasidiomycetes (class), Georgefischeriales (order), Tilletiariaceae (family), Tilletiaria (genus), Tilletiaria anomala (species)                                                                          \n",
       "6  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetes incertae sedis (no rank), Cantharellales (order), Botryobasidiaceae (family), Botryobasidium (genus), Botryobasidium botryosum (species)                                 \n",
       "7  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Jaapiales (order), Jaapiaceae (family), Jaapia (genus), Jaapia argillacea (species)                                                                        \n",
       "8  cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Basidiomycota (phylum), Agaricomycotina (subphylum), Agaricomycetes (class), Agaricomycetidae (subclass), Agaricales (order), Agaricineae (suborder), Strophariaceae (family), Galerina (genus), Galerina marginata (species)                                        \n",
       "9  Lichtheimia corymbifera (species), Lichtheimia (genus), Lichtheimiaceae (family), Mucorales (order), Mucoromycetes (class), Mucoromycotina (subphylum), Mucoromycota (phylum), Fungi incertae sedis (no rank), Fungi (kingdom), Opisthokonta (no rank), Eukaryota (superkingdom), cellular organisms (no rank)                                                                              \n",
       "10 cellular organisms (no rank), Eukaryota (superkingdom), Opisthokonta (no rank), Fungi (kingdom), Dikarya (subkingdom), Ascomycota (phylum), saccharomyceta (no rank), Pezizomycotina (subphylum), leotiomyceta (no rank), Eurotiomycetes (class), Chaetothyriomycetidae (subclass), Chaetothyriales (order), Herpotrichiellaceae (family), Exophiala (genus), Exophiala aquamarina (species)\n",
       "   Kingdom Phylum        Order             \n",
       "1  Fungi   Mucoromycota  Glomerales        \n",
       "2  Fungi   Basidiomycota Polyporales       \n",
       "3  Fungi   Ascomycota    Saccharomycetales \n",
       "4  Fungi   Basidiomycota Ustilaginales     \n",
       "5  Fungi   Basidiomycota Georgefischeriales\n",
       "6  Fungi   Basidiomycota Cantharellales    \n",
       "7  Fungi   Basidiomycota Jaapiales         \n",
       "8  Fungi   Basidiomycota Agaricales        \n",
       "9  Fungi   Mucoromycota  Mucorales         \n",
       "10 Fungi   Ascomycota    Chaetothyriales   "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "## Data export for manually curation of no kingdom rank entries\n",
    "    # 16 entries have no kingdom rank\n",
    "    # we have to curate them manually and add clade information instead: This includes: Discoba, Sar, Amoebozoa and Opisthokonta\n",
    "\n",
    "write.csv(DFA_clust_taxo,\"results/step1/DFA_clust_taxo_final_intermed.csv\")\n",
    "\n",
    "## Manually corrected data import\n",
    "data_taxo_DFA=read.csv('data/step1/DFA_clust_taxo_manually_corrected.csv')[-1]\n",
    "head(data_taxo_DFA,10)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "7db1c191",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table class=\"dataframe\">\n",
       "<caption>A data.frame: 7 × 3</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>tip.label</th><th scope=col>Taxo</th><th scope=col>nb_hits</th></tr>\n",
       "\t<tr><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;fct&gt;</th><th scope=col>&lt;int&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td>Acanthamoeba_castellanii</td><td>Amoebozoa    </td><td>  1</td></tr>\n",
       "\t<tr><td>Trypanosoma_cruzi       </td><td>Discoba      </td><td>  5</td></tr>\n",
       "\t<tr><td>Aureobasidium_pullulans </td><td>Fungi        </td><td>406</td></tr>\n",
       "\t<tr><td>Rattus_norvegicus       </td><td>Metazoa      </td><td> 10</td></tr>\n",
       "\t<tr><td>Sphaeroforma_arctica    </td><td>Opisthokonta </td><td>  1</td></tr>\n",
       "\t<tr><td>Perkinsus_olseni        </td><td>Sar          </td><td>  9</td></tr>\n",
       "\t<tr><td>Oryza_sativa            </td><td>Viridiplantae</td><td>  4</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 7 × 3\n",
       "\\begin{tabular}{lll}\n",
       " tip.label & Taxo & nb\\_hits\\\\\n",
       " <chr> & <fct> & <int>\\\\\n",
       "\\hline\n",
       "\t Acanthamoeba\\_castellanii & Amoebozoa     &   1\\\\\n",
       "\t Trypanosoma\\_cruzi        & Discoba       &   5\\\\\n",
       "\t Aureobasidium\\_pullulans  & Fungi         & 406\\\\\n",
       "\t Rattus\\_norvegicus        & Metazoa       &  10\\\\\n",
       "\t Sphaeroforma\\_arctica     & Opisthokonta  &   1\\\\\n",
       "\t Perkinsus\\_olseni         & Sar           &   9\\\\\n",
       "\t Oryza\\_sativa             & Viridiplantae &   4\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 7 × 3\n",
       "\n",
       "| tip.label &lt;chr&gt; | Taxo &lt;fct&gt; | nb_hits &lt;int&gt; |\n",
       "|---|---|---|\n",
       "| Acanthamoeba_castellanii | Amoebozoa     |   1 |\n",
       "| Trypanosoma_cruzi        | Discoba       |   5 |\n",
       "| Aureobasidium_pullulans  | Fungi         | 406 |\n",
       "| Rattus_norvegicus        | Metazoa       |  10 |\n",
       "| Sphaeroforma_arctica     | Opisthokonta  |   1 |\n",
       "| Perkinsus_olseni         | Sar           |   9 |\n",
       "| Oryza_sativa             | Viridiplantae |   4 |\n",
       "\n"
      ],
      "text/plain": [
       "  tip.label                Taxo          nb_hits\n",
       "1 Acanthamoeba_castellanii Amoebozoa       1    \n",
       "2 Trypanosoma_cruzi        Discoba         5    \n",
       "3 Aureobasidium_pullulans  Fungi         406    \n",
       "4 Rattus_norvegicus        Metazoa        10    \n",
       "5 Sphaeroforma_arctica     Opisthokonta    1    \n",
       "6 Perkinsus_olseni         Sar             9    \n",
       "7 Oryza_sativa             Viridiplantae   4    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "## Kingdom / Clade species count\n",
    "    # To investigate the distribution of DFA at the kingdom level, we now count how many protein hits were found for each kingdom or clade\n",
    "\n",
    "taxo_table=as.data.frame(t(table(data_taxo_DFA$Kingdom)))\n",
    "\n",
    "colnames(taxo_table)=c('tip.label','Taxo','nb_hits')\n",
    "\n",
    "# we add representative species name for each kingdom/clade - they are the same ones used to build the phylogenetic tree (tree tip.label)\n",
    "rep_species=c('Acanthamoeba_castellanii',\n",
    "              'Trypanosoma_cruzi','Aureobasidium_pullulans',\n",
    "              'Rattus_norvegicus',\n",
    "              'Sphaeroforma_arctica','Perkinsus_olseni',\n",
    "              'Oryza_sativa')\n",
    "\n",
    "taxo_table$tip.label=rep_species\n",
    "\n",
    "taxo_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "11ac50e3",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAYAAAD958/bAAAEDmlDQ1BrQ0dDb2xvclNwYWNlR2VuZXJpY1JHQgAAOI2NVV1oHFUUPpu5syskzoPUpqaSDv41lLRsUtGE2uj+ZbNt3CyTbLRBkMns3Z1pJjPj/KRpKT4UQRDBqOCT4P9bwSchaqvtiy2itFCiBIMo+ND6R6HSFwnruTOzu5O4a73L3PnmnO9+595z7t4LkLgsW5beJQIsGq4t5dPis8fmxMQ6dMF90A190C0rjpUqlSYBG+PCv9rt7yDG3tf2t/f/Z+uuUEcBiN2F2Kw4yiLiZQD+FcWyXYAEQfvICddi+AnEO2ycIOISw7UAVxieD/Cyz5mRMohfRSwoqoz+xNuIB+cj9loEB3Pw2448NaitKSLLRck2q5pOI9O9g/t/tkXda8Tbg0+PszB9FN8DuPaXKnKW4YcQn1Xk3HSIry5ps8UQ/2W5aQnxIwBdu7yFcgrxPsRjVXu8HOh0qao30cArp9SZZxDfg3h1wTzKxu5E/LUxX5wKdX5SnAzmDx4A4OIqLbB69yMesE1pKojLjVdoNsfyiPi45hZmAn3uLWdpOtfQOaVmikEs7ovj8hFWpz7EV6mel0L9Xy23FMYlPYZenAx0yDB1/PX6dledmQjikjkXCxqMJS9WtfFCyH9XtSekEF+2dH+P4tzITduTygGfv58a5VCTH5PtXD7EFZiNyUDBhHnsFTBgE0SQIA9pfFtgo6cKGuhooeilaKH41eDs38Ip+f4At1Rq/sjr6NEwQqb/I/DQqsLvaFUjvAx+eWirddAJZnAj1DFJL0mSg/gcIpPkMBkhoyCSJ8lTZIxk0TpKDjXHliJzZPO50dR5ASNSnzeLvIvod0HG/mdkmOC0z8VKnzcQ2M/Yz2vKldduXjp9bleLu0ZWn7vWc+l0JGcaai10yNrUnXLP/8Jf59ewX+c3Wgz+B34Df+vbVrc16zTMVgp9um9bxEfzPU5kPqUtVWxhs6OiWTVW+gIfywB9uXi7CGcGW/zk98k/kmvJ95IfJn/j3uQ+4c5zn3Kfcd+AyF3gLnJfcl9xH3OfR2rUee80a+6vo7EK5mmXUdyfQlrYLTwoZIU9wsPCZEtP6BWGhAlhL3p2N6sTjRdduwbHsG9kq32sgBepc+xurLPW4T9URpYGJ3ym4+8zA05u44QjST8ZIoVtu3qE7fWmdn5LPdqvgcZz8Ww8BWJ8X3w0PhQ/wnCDGd+LvlHs8dRy6bLLDuKMaZ20tZrqisPJ5ONiCq8yKhYM5cCgKOu66Lsc0aYOtZdo5QCwezI4wm9J/v0X23mlZXOfBjj8Jzv3WrY5D+CsA9D7aMs2gGfjve8ArD6mePZSeCfEYt8CONWDw8FXTxrPqx/r9Vt4biXeANh8vV7/+/16ffMD1N8AuKD/A/8leAvFY9bLAAAAOGVYSWZNTQAqAAAACAABh2kABAAAAAEAAAAaAAAAAAACoAIABAAAAAEAAANIoAMABAAAAAEAAANIAAAAAN/ryxkAAEAASURBVHgB7N0HnBRF2sfxZzazZJAlCUhGASWInqDIKXh4KGdEVIQzgQrmdGA6ODxERVFET0FBRUEx60uSpKhEEURQEAQkB4mysGnm7aeYGXp2F3Znd2Z3ZvpXnw9sT4earm/3efunqqtdHqsIBQEEEEAAAQQQQAABBBBAQOIwQAABBBBAAAEEEEAAAQQQOCpAQOJOQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCtAQOJWQAABBBBAAAEEEEAAAQS8AgQkbgUEEEAAAQQQQAABBBBAwCuQgAQCCCCAAAIIIIAAAtEikJ6eLl9++aWsXbtWsrKypH79+tK5c2epWrVqtDSB84xwAQJShF8gTg8BBBBAAAEEEEDgqMD06dPllltukc2bNweQlC1bVh566CF5/PHHA9bzAYGiCDDErihqHIMAAggggAACCCBQogIbNmyQHj16SFJSkkydOlV2794t27dvl2nTpskpp5wiTzzxhLz88sslek58WWwKuDxWic2m0SoEEEAAAQQQiEQBT+ZByVg1RrI3zRDJSpe48vUk6bRbJKH2+ZF4upxThAgMGzZMBg0aJO+8845cd911AWe1c+dOqVWrlrRq1UqWLFkSsI0PCAQrQA9SsGLsjwACCCCAAAJFFvBkHZL0Wf+U7N8+FsnYJ+LOFPf+X+XIgkGSuebdItfLgbEv8Ouvv5pGduzYMU9j09LS5OGHH5amTZsGbMvIyJAPP/xQbrvtNjnrrLPM8Lzx48fLnj17/PutWrVKbrzxRtGfc+bMkauuukomT57s386C8wQISM675rQYAQQQQACBUhPI/PkN8RzeaQWjrMBz8ORI5oqXxJ1ubaMgkI/Aqaeeatbec889smnTpjx7PPnkk6Z3yb5Be5o08MyaNcsEpHXr1smtt94qnTp1kuzsbLPr1q1bRUOTDs+78MILZcqUKXLkyBF7NSw7TIAhdg674DQXAQQQiFSB1157TUaMGCFnnHGGxMXx73eRep2Ke17D/75Jqqbm5F9NfLIkn3GvJDa4LP/trHW0wP79++X888+X5cuXi8vlkubNm5ug0759e7nooovyzGK3bds2M+xOA9Xzzz/vt7vjjjvklVdeke+//17atGkjM2fOlC5dukhiYqI89dRTpjepcuXK/v1ZcJ4As9g575rTYgQQQCAiBcaMGSNr1qwxfyLyBDmpkAg807WtVc9xfv2wnorWIXgUBPITqFixonz33Xfy1ltvyaeffirz5s2Tn376SV566SVJTk42PUXai1SvXj1zuE7m8PHHH8s555wTUF3jxo3NZ/swO13RtWtXue+++wL25YMzBY7zXyhnYtBqBBBAAIHSE2jdurV5uFqn8W3WrFnpnQjfHFaBcj8PFs/BleLK71s82RJX6egvr/ltZh0Cqamp5nkifaYoJyfH9CZpCJowYYIZXrds2TLz35GUlBTTo3TZZZeJrnvhhRdEh9f99ttvx53E4R//+AfACBgBAhI3AgIIIIBARAjokBktNWvWlLp160bEOXESoRfIKTNADn/d3+oqcgdW7kqQuAr1JT6tXeB6PiFgCegLYbXHSGeq8/0DSnx8vBkip8PkhgwZYpY1DC1cuNAMxdOZ7TQgzZ8/Xxo2bChnnnmmXHrppdKtWzcZPHhwHlcNXxQEVIBB3twHCCCAAAIIIFBiAvHVWkvymY9Zv4EkHf0TnyLiipe4io0k5byR5tmSEjsZvihqBPS5RA02vXv3zvec9R9Y+vTpY7bppAtadLidhqOnn35a1q5dK5MmTTIvkq1du7bZzl8IHE+AHqTjybAeAQQQQAABBMIikFjvYkmo2V6ydywSsd6JFFexgcSf1Cos30WlsSGgvUX6jiPtHVqxYoW0bNkyT8MWLbLuJ6vocF0t+kyjBiudtc5e9MWyFAROJEAP0ol02IYAAggggAACYRFwJVWUxDrWzGENryAchUU49ip98cUXRYNShw4d5Nlnn5WvvvrKTNLw0UcfiT4/NHHiROnevbv/XUgaqNxut+k10meP9B1HvXr1MlN+q472Kh06xKQgsXenFL9F9CAV35AaEEAAAQQQQAABBMIs0K5dO9FJXHTa7gcffDDg23SGO12v03T7nmccNGiQfPPNNzJq1CjzR8PV5ZdfbkKVTut9++23S5UqVcyfgMr44HgBApLjbwEAEEAAAQQQQACB6BC44IILzKx0v/76q+if9PR0M6nL6aefLrknWShfvrx8/fXXZv99+/aJBqxy5cqZhq5cuVJWrVolOuW3ThHu8VhzzFMQ8AoQkLgVEEAAAQQQQAABBKJGQJ8ratq0qX8o3YlOXHuTfM8k2ffTOlq0aGFfxTICfgGeQfJTsIAAAggggAACCCCAAAJOFyAgOf0OoP0IIIAAAggggAACCCDgFyAg+SlYQAABBBBAAAEEEEAAAacLEJCcfgfQfgQQQAABBBBAAAEEEPALEJD8FCwggAACCCCAAAIIIICA0wUISE6/A2g/AggggAACCCCAAAII+AUISH4KFhBAAAEEEEAAAQQQQMDpAgQkp98BtB8BBBBAAAEEEEAAAQT8AgQkPwULCCCAAAIIIIAAAggg4HQBApLT7wDajwACCCCAAAIIIIAAAn4BApKfggUEEEAAAQQQQAABBBBwugAByel3AO1HAAEEEEAAAQQQQAABvwAByU/BAgIIIIAAAggggAACCDhdgIDk9DuA9iOAAAIIIIAAAggggIBfgIDkp2ABAQQQQAABBBBAAAEEnC5AQHL6HUD7EUAAAQQQQAABBBBAwC9AQPJTsIAAAggggAACCCCAAAJOFyAgOf0OoP0IIIAAAggggAACCCDgFyAg+SlYQAABBBBAAAEEEEAAAacLEJCcfgfQfgQQQAABBBBAAAEEEPALEJD8FCwggAACCCCAAAIIIICA0wUISE6/A2g/AggggAACCCCAAAII+AUISH4KFhBAAAEEEEAAAQQQQMDpAgQkp98BtB8BBBBAAAEEEEAAAQT8AgQkPwULCCCAAAIIIIAAAggg4HQBApLT7wDajwACCCCAAAIIIIAAAn4BApKfggUEEEAAAQQQQAABBBBwugAByel3AO1HAAEEEEAAAQQQQAABvwAByU/BAgIIIIAAAggggAACCDhdgIDk9DuA9iOAAAIIIIAAAggggIBfgIDkp2ABAQQQQAABBBBAAAEEnC5AQHL6HUD7EUAAAQQQQAABBBBAwC9AQPJTsIAAAggggAACCCCAAAJOFyAgOf0OoP0IIIAAAggggAACCCDgFyAg+SlYQAABBBBAAAEEEEAAAacLEJCcfgfQfgQQQAABBBBAAAEEEPALEJD8FCwggAACCCCAAAIIIICA0wUISE6/A2g/AggggAACCCCAAAII+AUISH4KFhBAAAEEEEAAAQQQQMDpAgQkp98BtB8BBBBAAAEEEEAAAQT8AgQkPwULCCCAAAIIIIAAAggg4HQBApLT7wDajwACCCCAAAIIIIAAAn4BApKfggUEEEAAAQQQQAABBBBwugAByel3AO1HAAEEEEAAAQQQQAABvwAByU/BAgIIIIAAAggggAACCDhdgIDk9DuA9iOAAAIIIIAAAggggIBfgIDkp2ABAQQQQAABBBBAAAEEnC6Q4HQA2o8AApEpcPvtt8vYsWOlUqVKEhfHv+VE5lUK7Vnt3bvXVLht2zZp2bJlaCunNgQQQAABBAopQEAqJBS7IYBAyQqsW7dOsrOzJTk52fwp2W/n20pD4MCBA5KVlSXx8fGl8fV8JwIIIIAAAkaAgMSNgAACESnQsGFD+fLLL+WLL76QVq1aReQ5clKhFejXr5+89tprkpaWFtqKqQ0BBBBAAIEgBBi3EgQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAggggEBsCxCQYvv60joEEEAAAQQQQAABBBAIQoCAFAQWuyKAAAIIIIAAAk4Q2LZtmyxdulS2b99+3Oa63W6zz/Lly80+P//8s6xateq4+/s27N+/3xy3b98+36p8f3o8Hvn6669lw4YNZnthj8uvsjVr1hTq3PI7lnXOEyAgOe+a02IEEEAAAQQQQOCEAj/++KO0bdtW7rjjjuPup+FF93nwwQfNPtdee6306NHjuPv7NsyZM8ccN3PmTN+qfH9mZWXJ+eefL6+++qrZXtjj8qvspptukssvvzy/TSFb99lnn5ngF7IKqajUBAhIpUbPFyOAAAIIIIAAApEp0LlzZ6lZs6ZMnTpVDh48mO9JTp482ay/4YYbzM9LLrlEunfvnu++9pV16tQRPaZevXr21QUuF/W4AisO0Q4awJ577rkQ1UY1pSmQUJpfzncjgAACCCCAgPMEPO4syVz5mmSt/0wk60+Jq9hIks+4R+KrtXYeRoS2OD4+Xq677joZMWKEfP7552bZfqo6vO7DDz+UsmXLyhVXXGE2DR061L7LcZe11+mtt9467vbjbSjqccerj/UIHE+AHqTjybAeAQQQQAABBMIicGT+QMn6daJI5n4RT464962Ww1/3l5xdP4Tl+6i0aAK+nqH3338/TwU6vG7Hjh0mHGlI0vLUU0/Jf//7X/++uqx/9uzZY7Z17drVbNPhe3feeacsW7bMv68uvPHGGyaIaS+UBrMDBw4EbM/vuIceekjefvtt0eefdPnss8+Wq666SgoavqcV//LLL/LMM89Ily5d5KKLLpIhQ4bIt99+G/Cd+hzWjTfeKOvWrZMvvvjCnN9f/vIX6du3r6xevdrsq+el+2honD9/vlnevHmz2ZaRkSEbN24MqJMPkS9AD1LkXyPOEAEEEHCEgP5yoeW3336ThAT+7ylWL3pS+jpJ2/aduOTo9fa30+OWjOUjJbXzm/5VLJSuwBlnnCEtW7aUadOmmbBSoUIF/wnlHl6nGyZNmiTZ2dkyaNAgs9+UKVMkMzNT9NmchQsXSvPmzc16/d/4Sy+9ZJ4vatWqlVmnzwiNGzdOmjVrZvZ74YUX5N133zXbfH/ld9xrr70mDRo0MBM5NG7c2CzPnj1bPvroI9Ftt9xyi+/wgJ+LFi0y36//3bn00kulTJkyMnbsWHniiSdMO6655hqz/++//y7jx483www1IGnIq1Spkglz+lnDkbbxjz/+MPtrINJlddCivXDdunUTbR8legT4f6DouVacKQIIIBDTAj/8cLT34LLLLovpdjq9cbf+LU2G3lBXUpPj81C496/Ns44VpSugvUjaM/Ppp5+a54b0bHzD62rVqiUXXnjhCU9w8eLF0rFjRzMb3XnnnZfvvhrANBxpr8wrr7wicXFxcujQIfn73/+e7/65V+p/O4YNGyb/+te/zCadHU+Dl04eoc8FVa1aNfch8vrrr5t12rtTo0YNs7xp0yZp1KiRvPPOO+ILSL4DtcdM26KBUYt+nwZBDX4agDQE6rDETp06yYQJE8w+OhOgBrWHH37YfOav6BEgIEXPteJMEUAAgZgWuOeee+TZZ5+V1q1bm1+QYrqxDm7c6TX2S2LCFkvAk1chsVzedawpVYHrr7/eBA8dZucbcucbXqcBRMNMQeX555+XNm3aHHc37U3SYXrDhw/316efdXjeueeee9zjfBuqVatmQpzvs/bwaHjp16+fzJ07V6688krfJv/P3r17i8665wtHuqFy5cpy0kknmSGB/h29C3fffbc/HOkq7UnS7/jpp59MQMq9v35OSkqSjz/+WE4//fT8NrMuggUISBF8cTg1BBBAwEkCvXr1Ev1DiW0BT9YhOfR/1kxn2YcCGxqXKImnWOspESXg6yWaMWOGaM+Mho/8htcd76R15rkThSM9TkNGkyZNTN32enzD7+zr8ltu2rSpP1j5trdo0cIsrl2bf69khw4dTC+V9iTpcDt915I+E7Vz506pX7++rxr/T+1ZshcNUlpO9C4n7bnSHnF9BislJcV+OMsRLlBw7I/wBnB6CCCAAAIIIBA9Aq7EslLm3BEiCdaD/XFJIvHJIq44ia9+liS16Bs9DXHQmWrPkT5n88knn/iH12l48Q03OxGFPttTUNm1a1e+AUJDRWF6qE4++eQ8X6FBTsvxnmccPXq01K5dWwYMGCAaojQA6RTd+YUjrSc1NVV/BFU0bLVv3964BXUgO5e6AD1IpX4JOAEEEEAAAQScJRB/Uisp2+0zyd7+nXgy9kt85WYSX/Xov/g7SyI6WqvP8eiQNx1md8opp5jZ6/S5pFAVrfP777+X9PT0gCCyZcsWE8gK+p7169fn2UUnV9CivUu5iwYyHTKXlpYmup998olQtuvJJ580s9rpZBGU6BKgBym6rhdniwACCCCAQEwIaE9SYp0uktToKsJRhF/RcuXKmckOdOrsV1991UxGoLOzharo1NzaQ5V7iu3cs9gd7/tWrFghW7duDdis71nS3iffUDv7Rp0NLycnxzyDZA9HOiNd7nrsxwW7vGbNGnMOOnyQEl0CBKToul6cLQIIIIAAAgggUOICOswuKyvLTIGt7w2yT25Q3JN5/PHHJTk52UzJPXXqVPM80JgxY8w7ihITEwusXnue9FkfnQ1P3200cOBAmThxoujEL9o7lbtor5IO/dNpyTX06RA7fQeTvgupfPnysnv3btHeq2BL3bp1Zd68eTJ9+nT5888/zUx6OuOfvkuJEl0CBKToul6cLQIIIIAAAgggUOICOp23TtigxTebXahOQkOMBqPDhw+bqb31OaAHHnjATJGtgaWgotNsax06Lfipp55qXkqrgUmHuOVX9PkkDWA6wYKGPX1/ks5Ip71j+rJb7WHS9cGWu+66y9SpM9xpHVqnzsK3atWqYKti/1IWcHmsUsrnwNcjgAACeQRuv/12+d///if6fovCzmSUpxJWIIAAAghEjcCRI0dkyZIlZmIFne5fe5UKKhp2NLx9+OGHpudH/z9DA1buWefyq0d7inTmOp1pz/6skr6/SIfg5Tf5Q3712NdpyNMXxfqO1V+z9+7dK1WqVLHvxnKECzBJQ4RfIE4PAQQQQAABBBBwgoDOWleY9x4dz0Kn3g6m50f379y5c57qatasmWddYVfo0D1fONJjXC4X4aiweBG0H0PsIuhicCoIIIAAAggggAACCCBQugIEpNL159sRQAABBBBAAAEEiiigQ+N0iBwFgVAKMMQulJrUhQACCCCAAAIIIFBiAgsXLiyx7+KLnCNAD5JzrjUtRQABBBBAAAEEEEAAgQIECEgFALEZAQQQQAABBBBAAAEEnCNAQHLOtaalCCCAAAIIIIAAAgggUIAAAakAIDYjgAACCCCAAAIIIICAcwQISM651rQUAQQQQAABBBBAAAEEChAgIBUAxGYEEEAAAQQQQAABBBBwjgAByTnXmpYigAACCCCAAAIIIIBAAQIEpAKA2IwAAggggAACCCCAAALOESAgOeda01IEEEAAAQQQQAABBBAoQICAVAAQmxFAAAEEEEAAAQQQQMA5AgQk51xrWooAAggggAACCCCAAAIFCBCQCgBiMwIIIIAAAggggAACCDhHgIDknGtNSxFAAAEEEEAAAQQQQKAAAQJSAUBsRgABBBBAAAEEEEAAAecIEJCcc61pKQIIIIAAAggggAACCBQgkFDAdjYjEBECQ4cOlREjRkiDBg0kPj4+Is6JkwivwJo1a8wX7NmzJ7xfRO0IIIAAAggggIBNgIBkw2AxcgVmzZol+/btkx9++EFcLlfknihnFjIBt9tt6jp48GDI6qQiBBBAAAEEEECgIAECUkFCbI8IgebNm8vcuXNlyZIl0qZNm4g4J04ivAK33367/O9//5N69eqF94uoHQEEEEAAAQQQsAnwDJINg0UEEEAAAQQQQAABBBBwtgABydnXn9YjgAACCCCAAAIIIICATYCAZMNgEQEEEEAAAQQQQAABBJwtQEBy9vWn9QgggAACCCCAAAIIIGATICDZMFhEAAEEEEAAAQQQQAABZwsQkJx9/Wk9AggggAACCCCAAAII2AQISDYMFhFAAAEEEEAAAQQQQMDZAgQkZ19/Wo8AAggggAACCCCAAAI2AQKSDYNFBBBAAAEEEEAAAQQQcLYAAcnZ15/WI4AAAggggAACCCCAgE2AgGTDYBEBBBBAAAEEEEAAAQScLUBAcvb1p/UIIIAAAggggAACCCBgEyAg2TBYRAABBBBAAAEEEEAAAWcLEJCcff1pPQIIIIAAAggggAACCNgECEg2DBYRQAABBBBAAAEEEEDA2QIEJGdff1qPAAIIIIAAAggggAACNgECkg2DRQQQQAABBBBAAAEEEHC2AAHJ2def1iOAAAIIIIAAAggggIBNgIBkw2ARAQQQQAABBBBAAAEEnC1AQHL29af1CCCAAAIIIIAAAgggYBMgINkwWEQAAQQQQAABBBBAAAFnCxCQnH39aT0CCCCAAAIIIIAAAgjYBAhINgwWEUAAAQQQQAABBBBAwNkCBCRnX39ajwACCCCAAAIIIIAAAjYBApINg0UEEEAAAQQQQAABBBBwtgABydnXn9YjgAACCCCAAAIIIICATYCAZMNgEQEEEEAAAQQQQAABBJwtQEBy9vWn9QgggAACCCCAAAIIIGATICDZMFhEAAEEEEAAAQQQQAABZwsQkJx9/Wk9AggggAACCCCAAAII2AQISDYMFhFAAAEEEEAAAQQQQMDZAgQkZ19/Wo8AAggggAACCCCAAAI2AQKSDYNFBBBAAAEEEEAAAQQQcLYAAcnZ15/WI4AAAggggAACCCCAgE2AgGTDYBEBBBBAAAEEEEAAAQScLUBAcvb1p/UIIIAAAggggAACCCBgEyAg2TBYRAABBBBAAAEEEEAAAWcLEJCcff1pPQIIIIAAAggggAACCNgECEg2DBYRQAABBBBAAAEEEEDA2QIEJGdff1qPAAIIIIAAAggggAACNgECkg2DRQQQQAABBBBAAAEEEHC2AAHJ2def1iOAAAIIIIAAAggggIBNgIBkw2ARAQQQQAABBBBAAAEEnC1AQHL29af1CCCAAAIIIIAAAgggYBMgINkwWEQAAQQQQAABBBBAAAFnCxCQnH39aT0CCCCAAAIIIIAAAgjYBAhINgwWEUAAAQQQQAABBBBAwNkCBCRnX39ajwACCCCAAAIIIIAAAjYBApINg0UEEEAAAQQQQAABBBBwtgABydnXn9YjgAACCCCAAAIIIICATYCAZMNgEQEEEEAAAQQQQAABBJwtQEBy9vWn9QgggAACCCCAAAIIIGATICDZMFhEAAEEEEAAAQQQQAABZwsQkJx9/Wk9AggggAACCCCAAAII2AQISDYMFhFAAAEEEEAAAQQQQMDZAgQkZ19/Wo8AAggggAACCCCAAAI2AQKSDYNFBBBAAAEEEEAAAQQQcLYAAcnZ15/WI4AAAggggAACCCCAgE2AgGTDYBEBBBBAAAEEEEAAAQScLUBAcvb1p/UIIIAAAggggAACCCBgEyAg2TBYRAABBBBAAAEEEEAAAWcLEJCcff1pPQIIIIAAAggggAACCNgECEg2DBYRQAABBBBAAAEEEEDA2QIEJGdff1qPAAIIIIAAAggggAACNgECkg2DRQQQQAABBBBAAAEEEHC2AAHJ2def1iOAAAIIIIAAAggggIBNgIBkw2ARAQQQQAABBBBAAAEEnC1AQHL29af1CCCAAAIIIIAAAgggYBMgINkwWEQAAQQQQAABBBBAAAFnCxCQnH39aT0CCCCAAAIIIIAAAgjYBAhINgwWEUAAAQQQQAABBBBAwNkCBCRnX39ajwACCCCAAAIIIIAAAjYBApINg0UEEEAAAQQQQAABBBBwtgABydnXn9YjgAACCCCAAAIIIICATYCAZMNgEQEEEEAAAQQQQAABBJwtQEBy9vWn9QgggAACCCCAAAIIIGATSLAts4gAAggggAACCCCAQEgFli5dmqe+smXLSuPGjSUuLu+/1e/fv1/WrVsnDRo0kEqVKuU5tiRWbNmyxZzDWWedJSkpKSXxlXxHBAnkvSsj6OQ4FQQQQAABBBBAAIHoFcjOzpa2bdvm+dOsWTOpUKGCXHDBBfLzzz8HNHDOnDlm/5kzZwasL8kPkydPlvPPP1+2bdtWkl/Ld0WIAD1IEXIhOA0EEEAAAQQQQCBWBbS3aOjQof7mbd26VX744QeZOHGinH322fLee+/JxRdfbLbXqVNHbrjhBqlXr55/fxYQKEkBAlJJavNdCCCAAAIIIGAEcvavk+xNX4onY5/EV24mCfW6iiueoUyxentUrVpVevTokad5N998s1xyySWiP3VYXZkyZUzv0VtvvZVnX1YgUFICDLErKWm+BwEEEEAAAQSMQOa6j+Xwl70ka807kr3+E8lYNkLSp18r7sO7EXKYQMeOHaVv375mKNvYsWNN63/88Ue58847ZdmyZX6N3377TW655RYTnrTHSZd1Xe7y+++/m21nnHGGdO3aVUaPHi36PJG96LC5O+64Q/T5otNPP1369Okjq1evtu/iX87KypInn3xS9DzPOeccuf/++0WfkbKXjIwM+fDDD+W2224zdeq5jR8/Xvbs2WPfjeUoEiAgRdHF4lQRQAABBBCIdgH3wd8l84enrWZ4RNxZR5tj/fQc3ikZ3z8Z7c3j/IsgoAFFy7fffmt+avB56aWXZO3atebzxo0b5cwzz5Tp06dLmzZtpHXr1jJlyhSz7NtHd1yyZInZ9sEHH4g+43To0CETtPQZqO3bt5u6NHRpeNJ9TjvtNP+zTq1atTJ1mp1sf1111VUyatQoqVu3rpQvX16ef/558x379u3z73XdddeJ7jdr1iwTkLQn7NZbb5VOnTqJPoNFiT4BhthF3zXjjBFwhIDv/8wefPDBUpvFyBHQNBKBEhboduoB6X6aR5Jy/wbiyZGc7QvEk31EXAkMtSvhy1KqX1e/fn3z/Ros8ivjxo2TvXv3yldffSUtW7Y0uzzxxBPmGaV3331XHn/8cbPurrvukszMTFm4cKEJP7ry448/liuuuELGjBkjjz32mAwYMEBycnJk8eLFZpY83Wfnzp3Svn176devn5kwoly5crralPj4eLOucuXK5rMGN+3dGjZsmAwfPtz0fH300Udyzz33mPDkPcz0UL3yyiuivWEa6ijRJZD7P0/RdfacLQIIxKxAenq6aVtpzmIUs7g0DIFSFGh5fR2R02pYZ5DfIBaPFZDSCUileH1K46s1kKSmpppwk9/3//nnn2b1p59+Ki1atBCXyyU1a9aUX3/9VZKSksw2nfBh/vz5Znid9gz5ymWXXSYanDTg6D7aS9WzZ09/ONL90tLSpHfv3qKhS+vo0qWL73AzsYQvHOnKf/7zn3LvvffKF198YQKSfr+GMB1+Zy86KYUWhtnZVaJnmYAUPdeKM0XAUQI6lEL/z08f2KUggEDsCCT+sUAS1460RtjlHXrkSq4qcSlVYqextKRQAkeOHBH9RzEdCpdf0Wd6dNIG7QHSnqDOnTvL3/72N/OMkU4VrmXNmjXmpz5TZC8apl544QWzSmfK06LTd+cuvl4efRbJHpB8PVa+/TXM6bNLGra06OQTGsJ06J5+j/aC6RBBHe5HiV4BAlL0XjvOHIGYF/D9C1zMN5QGIuAgAU/tK+Twrs/Fbc1iFxCSXHGS1OpeB0nQVJ+Ab6RAu3btfKsCfurzRBqAJk2aJO+//76888478sYbb5j3KL399tvSvXt3/zNG+qzQ8cqOHTvMpiZNmuTZRYOUFo/HejbOVnSIXe6izyL5ig7P04CkPU8NGzY0z0pdeuml0q1bNxk8eLBvN35GmUB+/dtR1gROFwEEEEAAAQSiRcBlBaEy578sCXW7WqPsjg6PcqXWkJSzn5TEOp2jpRmcZwgFdMY3Ldozk1/RiQ40lOgscbNnzzbD1jQs6fNG/fv3N4f43pmUX8+NhqipU6f636vkmwzC/l0acLQ0aNDAvjrP7HZut9uEId9zUzrDnR779NNPm0kl9Lz0majatWsH1MOH6BIgIEXX9eJsEUAAAQQQiHoBV2JZSWn3qJS9/Kujf/7+iSSc/NeobxcNCE5AJ0t4+OGHzRTZOgucb5hb7lp0Wm/78Dt9Xumaa66Rv/71r/LHH3+YXh/tfdLngebMmRNw+Lx588zzRRpidPa7hIQEM9tcwE7Wh6+//tqEsNzPEn3yyScBuy5dulQOHDhghvjpBu3ZiouLM7PW2XecNm2a/SPLUSbAELsou2CcLgIIIIAAArEiYIY1xSfHSnNoxwkE9P1EAwcO9O+hs9ItWrTIPMujU3i/+eabZvIF/w62hX/84x9mAoWnnnrKvLNIQ45O863PqurwOr2PtMdG3200cuRIeeihh6RXr17y/fffm4kU9LkhnYRBh9/pO5defvllcy76cloNaa+99ppokHrkkUekSpXAZ+B0W6VKlUSHze3atcvMVle2bFl54IEHzBnq9OAahrTXSGey0ynJX3/9dX8I02nIdbpvPWdKFAlYYy0pCES8gNWFroOCPdZ/7CL+XDlBBBBAAAEEEDgqYL1o1fz/t/5/uP1PxYoVPVbw8VgTG3is53gCuKxZ4cy+kydPNuutSRw8V199dcDxWpc12YLHeumr/1j9rrvvvttj9ej497WGzHlmzJjh38eaDMJjTeft3671WAHL8+ijj3qssOTfz3rfkf8cqlev7t/fCmIea4pw/35Wb5Ln3HPP9W+3nlnyWL1hnk2bNnmsZ6fMemvSBv/+LESHgEtP07o5KAhEtIC+t0Dfhq3/GnS8LviIbgAnhwACCCCAAALFEtDeGB3SZgUhMyGCTvmdX9FpwVesWGGm9m7UqFG+vTf6rr3ly5eL9i7pzHf2iRdy16m/KuusdSkpKeYFtDqkzl50u85ipy+P1aF+WqcWfV5p1apV5p1MuY+xH89y5AkQkCLvmnBG+QgQkPJBYRUCCCCAAAIIIIBAyAUCI3DIq6dCBBBAAAEEEEAAAQQQQCB6BAhI0XOtOFMEEEAAAQQQQAABBBAIswABKczAVI8AAggggAACCCCAAALRI0BAip5rxZkigAACCCCAAAIIIIBAmAWiclL2hQsXmnntdTYz8w6FMCNRfekL+N56rTPCUBBAAAEEEEAAAQQQCJdAVM5ip285XrBgQbhMqDeCBfSFbbfeemsEnyGnhgACCCCAAAIIIBDNAlHZg9S6dWsTkEaMGCFNmjSJZn/OvZACL730knljdtu2bQt5BLshgAACCCCAAAIIIBC8QFQGpMTERNPS8847z7yQK/hmc0S0CUybNs0EpGg7b84XAQQQQAABBBBAILoEmKQhuq4XZ4sAAggggAACCCCAAAJhFCAghRGXqhFAAAEEEEAAAQQQQCC6BAhI0XW9OFsEEEAAAQQQQAABBBAIowABKYy4VI0AAggggAACCCCAAALRJUBAiq7rxdkigAACCCCAAAIIIIBAGAUISGHEpWoEEEAAAQQQQAABBBCILgECUnRdL84WAQQQQAABBBBAAAEEwihAQAojLlUjgAACCCCAAAIIIIBAdAkQkKLrenG2CCCAAAIIIIAAAgggEEYBAlIYcakaAQQQQAABBBBAAAEEokuAgBRd14uzRQABBBBAAAEEEEAAgTAKEJDCiEvVCCCAAAIIIIAAAgggEF0CBKToul6cLQIIIIAAAggggAACCIRRgIAURlyqRgABBBBAAAEEEEAAgegSICBF1/XibBFAAAEEEEAAAQQQQCCMAgSkMOJSNQIIIIAAAggggAACCESXAAEpuq4XZ4sAAggggAACCCCAAAJhFCAghRGXqhFAAAEEEEAAAQQQQCC6BAhI0XW9OFsEEEAAAQQQQAABBBAIowABKYy4VI0AAggggAACCCCAAALRJUBAiq7rxdkigAACCCCAAAIIIIBAGAUISGHEpWoEEEAAAQQQQAABBBCILgECUnRdL84WAQQQQAABBBBAAAEEwihAQAojLlUjgAACCCCAAAIIIIBAdAkQkKLrenG2CCCAAAIIIIAAAgggEEYBAlIYcakaAQQQQAABBBBAAAEEokuAgBRd14uzRQABBBBAAAEEEEAAgTAKEJDCiEvVCCCAAAIIIIAAAgggEF0CBKToul6cLQIIIIAAAggggAACxRY4fPiwbNy4UfRnqMr27dtl27ZtBVZ35MgR893p6ekF7lsaOxCQSkOd70QAAQQQQAABBBBAoBQFZs2aJaeccorMmDEjZGfRvXt36dq1a4H1zZ0713z3lClT/PuuWLFCxo4d6/9cmgsEpNLU57sRQAABBBBAAAEEEHCYwEknnSRdunSR6tWr+1vetm1bWbhwof9zaS4klOaX890IIIAAAggggAACCCDgLIEzzzwzT89VdnZ2xCAQkCLmUnAiCCCAAAIIIIAAAggUT+DXX3+VCRMmyJ133inff/+96DC2nTt3Srt27eT222+XMmXK5PmC9957T6ZOnSopKSly4YUXytVXX51nn2BW7N27V/73v//JsmXLpH79+qJD79q3b++vYu3atfL222/LNddcI1WrVpWXX35ZPB6POd8nnnhCbrnlFqlTp45s2rTJDLtbtWqVVKxYUVq2bCk333yzlCtXzl9XOBYYYhcOVepEAAEEEEAAAQQQQKAUBDQgDRkyRB5++GHzPNDXX38tP/zwg9x///1y7rnnitvtDjirf//739K3b185ePCg2a9Hjx7Sq1evgH2C+bBjxw7R4XKjR4+WrVu3yjPPPCOdOnWS//u///NXowFJz1GDj07UoOeoRSd50OUDBw6I7tO6dWt58cUXzT5r1qyRBx54QFq1aiV//PGHv65wLBCQwqFKnQgggAACCCCAAAIIlKKABhKd+EDD0S+//GJC0NKlS2XmzJkBZ6UhRtd/+OGH5hmggQMHyjvvvCOffPJJwH6F/aAB6ZJLLpHNmzfLvHnzZOXKlaLD50aNGpVvFdrDNGfOHHG5XNKtWzez3Lx5c3nttddMaNM2aFs0OL377ruybt06mThxYr51hWolASlUktSDAAIIIIAAAggggECECOhwuhYtWvjPxjdsTqf2tpf77rtPGjZs6F/12GOPSaVKleT111/3rwtmIS4uzvQa+Y5p1qyZtGnTxoQ037rC/NR6MjMz5ZtvvpGcnBxziLZBpxEfMGBAYaoo8j4EpCLTcSACCCCAAAIIIIAAApEp0KhRo4ATS0tLM59zv/dIn02yF31GqUmTJrJ69Wr76kIv16tXT5KTkwP213MJdlicPmtUu3Ztufbaa0XPvWfPnua5pdTU1IC6w/GBgBQOVepEAAEEEEAAAQQQQKAUBXIHCR3CpkUnQ7CXChUq2D+aZZ0EIXeQyrPTcVaUL1/+OFuCW924cWMzacPQoUNFlydPniy9e/c2y99++21wlQW5NwEpSDB2RwABBBBAAAEEEEAgVgR0YoTcRYfh5e6Byr1PSXyuXLmyPPLII7JgwQLZtWuXmRlv9+7d8q9//SusX09ACisvlSOAAAIIIIAAAgggELkC2jNjL/rMj06EoDPPlWSJj483zxz5vlOnANepvg8dOmRWValSRfr16yf6TFN+oc53XCh+EpBCoUgdCCCAAAIIIIAAAghEoYDOCPfQQw/J8uXLzUx2Os13rVq1zLTgJdkc7S3S2ez0nUj6/iOdllzf36RTjussdgsXLjRTl+vU4L4JJ8J1fgSkcMlSLwIIIIAAAggggAACES7wxhtvyEcffWTeL3TVVVdJgwYNRJ/xCffLWHOzPProo2Yih/79+8usWbPkrLPOEn3+SKcl12nD//KXv5ghdoMGDTLrcx8fys8u60GtwCe1Qll7mOq6++67zUujFi1aZN4KHKavodoIEtDpHPWFY/pGaJ0qkoIAAggggAACCCAQOoH169eLTthQtWrV0FUaZE06nfeePXvkpJNOMu9F0sN1qm/tUdKi4c032YRZEaa/EsJUL9UigAACCCCAAAIIIIBAlAjoC1tLu+hzSNWqVQs4jaSkpID3NAVsDNMHAlKYYKkWAQQQQAABBBBAAIFoFli7dq2MGzeuUE247LLLYmZkFwGpUJecnRBAAAEEEEAAAQQQcJaAvgtJp/wuTDl48GBhdouKfQhIUXGZOEkEEEAAAQQQQAABBEpWoGXLljJhwoSS/dII+DZmsYuAi8ApIIAAAggggAACCCCAQGQIEJAi4zpwFggggAACCCCAAAIIIBABAgSkCLgInAICCCCAAAIIIIAAAghEhgABKTKuA2eBAAIIIIAAAggggAACESBAQIqAi8ApIIAAAggggAACCCCAQGQIEJAi4zpwFggggAACCCCAAAIIIBABAgSkCLgInAICCCCAAAIIIIAAAghEhgABKTKuA2eBAAIIIIAAAggggAACESBAQIqAi8ApIIAAAggggAACCCCAQGQIEJAi4zpwFggggAACCCCAAAIIIBABAgSkCLgInAICCCCAAAIIIIAAAghEhgABKTKuA2eBAAIIIIAAAggggAACESCQEAHnwCkggAACCCCAAAIIIIBACQu4D/4u2dvmifvABvFkH5a4srUkPq2d9aeNuFzxJXw2kfN1BKTIuRacCQIIIIAAAggggAACYRfwZOyTI0uHS86WuSLxSSI5GeY7c1wJkrXmHXGlVpfkto9IQlrbsJ9LJH4BQ+wi8apwTggggAACCCCAAAIIhEHAfWibpH/ZS3KsniMRjz8cma/yZFurcsRzaKscmXenZG34vzCcgfWVOTkydOhQ2bNnT576V69eLSNGjJB33nlH9u/fn2f7vn375K233pIXXnhBdN9wFAJSOFSpEwEEEEAAAQQQQACBCBPwuLPl8Ly7xZOxV8RaPmHxuCXj+ycl54+fTrhbUTY++OCD8thjj+UJQMOGDZPmzZvLokWLZOTIkdKhQwfZuXOn/ytWrlwpNWrUkFGjRsn8+fOlTZs2Mm3aNP/2UC0QkEIlST0IIIAAAggggAACCESwQNZvn4gnfZvpJSrUaXo8kmENxQtV2bRpk3Tr1k1eeeWVPFWuWbNGBg8eLLNnz5b33ntPvvvuOylTpow899xz/n1vuukm6du3rwlQkyZNkkGDBsmAAQPEY51nKAsBKZSa1IUAAggggAACCCCAQIQKZK2bbPUcZQVxdh5x718rOft/C+KY4+968803i9vtls8//zzPTtOnT5cGDRpIx44dzbbExETp3bu3TJw40Xzevn27CUb9+vWzJpBwmXVa37p168z6PBUWYwUBqRh4HIoAAggggAACCCCAQDQIeLIOiefgxuBP1ZrEwb17afDH5XPE2LFjZerUqXLyySfn2bp+/Xpp2LBhwHoNTFu2bDGhasOGDWabfR8dbqe9TNozFcpCQAqlJnUhgAACCCCAAAIIIBCBAp6MvBMiFOo03TniPlLEY3N9Qd26dXOtOfZx48aNUrVq1WMrrKXKlSubCR12794tGpBSU1MlJSUlzz47duwIWFfcDwSk4gpyPAIIIIAAAggggAACkS6QkFq0M4yLF1dC2aIdG8RRSUlJkpUVOPwvMzPT1FC+fHnJb7tu1GN0eygLASmUmtSFAAIIIIAAAggggEAECriSq4gklgv+zKwepLiKDYI/LsgjatWqlWfab50GvHr16mYYnW7XMHTw4MGAmnWf+vXrB6wr7gcCUnEFOR4BBBBAAAEEEEAAgQgX0IkNEk7uLOJKDO5MrWeQ4tPODO6YIuzdokULWbJkiWRnH5t+fMGCBf7nkho3bmx6kXSdryxevNgMwdNnlUJZCEih1KQuBBBAAAEEEEAAAQQiVCDp1JusgBTElNhxiZLUvK+4rJ/hLj179jRfMXz4cDMpw08//STjxo0zU3nrBn0+6frrr5chQ4aY9yelp6ebdyn16dNHateuHdLTIyCFlJPKEEAAAQQQQAABBBCITIG41DRJbvdvKyQVIgJYoUh7jhIbXVMijdHZ6CZPniyjR4+WSpUqSZcuXaR///7mvUm+E9DwlJycLGlpaaJD7nQqcPt7knz7FfdnQnEr4HgEEEAAAQQQQAABBBCIDoHEOp2tfJQgRxY9Yb0w1p3Pe5E0PHkkoe7fJLnNw/53DoWydc2aNcv35a6dOnWSrVu3mmm7tVcoLi4wyFWrVk1mzpxpnlVKSEiQChUqhPK0/HURkPwULCCAAAIIIIAAAgggEPsCCbU7Sdm/fyaZa9+X7E0zxHNoq5WJckSSKkhC9XOsXqOrJb5qi1KDqFOnzgm/u0oVa8KJMBYCUhhxqRoBBBBAAAEEEEAAgUgUcCVXlOTmt5o/en4eKyC5XPGReKolfk6B/VYl/vV8IQIIIIAAAggggAACCJS2AOHo2BUgIB2zYAkBBBBAAAEEEEAAAQQcLkBAcvgNQPMRQAABBBBAAAEEEEDgmAAB6ZgFSwgggAACCCCAAAIIIOBwAQKSw28Amo8AAggggAACCCCAAALHBAhIxyxYQgABBBBAAAEEEEAAAYcLEJAcfgPQfAQQQAABBBBAAAEEEDgmQEA6ZsESAggggAACCCCAAAIIOFyAgOTwG4DmI4AAAggggAACCCCAwDEBAtIxC5YQQAABBBBAAAEEEEDA4QIEJIffADQfAQQQQAABBBBAAAEEjgkkHFtkCQEEEEAAAQQQQAABBJwg4HFnS86upZKzc4m4D20RyTkirtSaEl+5qSTUPE9cyZWcwJBvGwlI+bKwEgEEEEAAAQQQQACB2BPweDySveFzyVgxWiTrT6uBLhFPtrehLsmOS5QM938loX53SW5xmxWUKsceQgEtIiAVAMRmBBBAAAEEEEAAAQRiQcCTnS5HFjxq9RotFnFn5dMkj7U+06zP3jhFsrfMlTLnjpT4Kqfms2/sruIZpNi9trQMAQQQQAABBBBAAAEj4PHkyOFv7jtBOMoFpQEq84Ac/uo2cR9Yn2tjbH+kBym2ry+tQwABBBBAAAEEQiawf/9+WbdunaSmpkpn2DlwAABAAElEQVSzZs1CVm+4KtqyZYs537POOktSUlLC9TVRUW/myjHi3rPyOD1Hx2uC1aOUkyWH590jqRd/IC5r+F1xy/r16+WTTz6RhIQEueSSS6R+/foBVa5evVq++OILqVGjhtlesWLFgO379u2Tzz77TPbu3Stdu3aVpk2bBmwPxQd6kEKhSB0IIIAAAggggIADBIYPHy5t27Y1fw4ePBjxLZ48ebKcf/75sm3btog/13CeoPvwLslaPSHIcOQ7I7d4MvZI1rqPfCuK/PPqq6+WM844Q5YuXSrjxo2T0047TaZMmeKvb9iwYdK8eXNZtGiRjBw5Ujp06CA7d+70b1+5cqUJTqNGjZL58+dLmzZtZNq0af7toVogIIVKknoQQAABBBBAAIEYFnC73fLWW29JYmKipKeny4cffhjDrY2tpunzRFb3T9EbZQ23y/p1UtGPt47UUPTxxx/LqlWr5O233zafL7/8crnnnntMvWvWrJHBgwfL7Nmz5b333pPvvvtOypQpI88995z/e2+66Sbp27evCVCTJk2SQYMGyYABA0QnnghlKYZUKE+DuhBAAAEEEEDAaQLu9B2Ss2eVeHIynNb0qGzvl19+KTpk7bHHHpP4+HiZMMHqkaBEhUD25ln+yReKesKe9G3i/tOaDryIZceOHTJkyBA5+eST/TX89a9/lY0bN5qAM336dGnQoIF07NjRbNcg3rt3b5k4caL5vH37dhOM+vXrJy6XNfOeVW6++WYzhFJ7nEJZeAYplJrUhQACCCCAAAIFCng8bslYMlTMv2rrFMOJZSXlnOGSkNa2wGPZofQExo8fL8nJyXLnnXfKvHnzZNasWbJ58+aAX3i1l0CHP913332iz5q8//77snv3bvMsif7LvwasZ5991vQe6C/Dt912m5x33nkBjdLhcP/5z39kyZIlcuTIEWndurXpKcj9rElGRoY8//zz8u2338qBAwfk3HPPlRtvvFEaNWoUUJ9+yMrKkieffFL0l3Bdbt++vTz++OOS+/kW7dn44IMP5Oeff5Z69erJpZdeKv379zeBUOvR5160bccr2p6zzz7bbPa14+KLLzb1HO+YklhfnGDjPz/r+SP3n5slrlxt/6pgFtRB/9iL9gK1a9fOBB69Xxo2bGjfbAKT3jPae7lhwwazzb6PPqekvUybNm3yuwdUUMQPBKQiwnFYyQro/zC0rF27VpKSkkr2y/k2BBBAAIGQCpTb97VU3TFDo5FVrKEx1rtYjsx/WMpeMkVc8fw3PqTYIapMH4zXB+s1MFSqVEmuueYa0R6ld999Vx566CH/t/z++++iQUrDwYwZM8zzPz/88IMJJl9//bUJVtnZ2Sb06HCrTz/9VBYvXiynnnp0Gully5bJRRddZOr7+9//bn5x1no0aOmQPl2nZdeuXfK3v/1NfvrpJ7NOf1EeO3asvPLKK/LRRx9Jp06dzH6+v6666irzLEvnzp3NTw1W+v0a6LQ9WjTAjRkzxvRg6C/yWrcO/9L9Zs6caUJSZmamCW6+en0/f/nlF9F26ZAxLfZ2XHvttb7dSu+nNb13sYs1RM+Tub/Y1fgq0Gug98SCBQvMKu1Jqlq1qm+z+Vm5cmXJyckxIVsDkk4OknuyDd1He6dCWQhIodSkrrAJ6H/AtOh/kCkIIIAAAtEtMP7eRnLFOVWsZyKORiTTGiskufevc9z7VqLlSuq/9Gtvzg033GBO+YorrpA77rjDPEtiD0i+9ixfvlx0NrLGjRuLznzXokUL0QkT7rrrLnnmmWfMP3ZqnRoeNHz4ApI+T6K/EGto0h4mLfqQvvb46NAq7dkpV66cPProo6LBa+7cuSaE6X46BEt7kfQ5FZ1pzzcMS7fpkEA9Vn+Z1vLSSy+ZnjCdFEAnntBz0HCkgUh/cfeVESNGyAMPPGD2v/vuu6V69eqyYsUK32bz8/XXX5dbbrlF1EQDpJb82mE2lNZfVi/t0ZfCFuMErJ5fV3LgjHJFrU2fNVJ3DbM66YcW/Qdw7d2zFw2kWsqXL5/vdt2mx+j2UBYCUig1qStsAtqd/dRTT5l/cYqL49G5sEFTMQIIIFACAlVrbZMczx+S4Ap8sNqVEvivxyVwKnxFIQW0V+ikk07yD5HSoKE9PTods/aWtGrVKqCm22+/3YQjXanD2Fq2bGmG4w0cONA/EkSnaNaiPTVaNPDocLmePXv6w5GuT0tLM8+iPPHEE2bmsjPPPNP0FunU3TpDna9oL5KGo0ceecTsp6HKV4YOHeoPR7run//8p9x7773m/PUXdR0WqEWHyNmLBjqtT3uwNCDlLt98840Jitp+ncBCQ9nx2pH72JL8HFf2ZHHv+6V4X2lN1BBXrk6x6tARQXpvaDieOnVqwPWrVauWmcDB/gV79uwxoVSH0el2DUM6e6I9EOk+uacKt9dRlGUCUlHUOKbEBXr06CH6h4IAAgggEP0C7kPbJP3LXkcfGteXUVrPNiTU7SpxqWnR37gYbIH2vCxcuFA0mGhvia/ocCct+txO7oCkPUf2UrZsWTP7nYYYX9H34GjRniktOouZFnvoMSusv3Q6Zy3aK6W/HOsv2jpzmfbU2IsO09KiP+0BSQOavWgvlAYsDTNa9Lv13HI/56QTBWjvl35v7qLfob1GOkRPhwpqG7WcqB256yipzwl1OlvvfP2tWBM1uMrWlriytYp1ytoDqcPqNAirq73oZw2ZOlTRd2/o8DvfM0d6T2kvk67r0qWLOVR7GrXH0dfbaK+vOMsEpOLocSwCCCCAAAIIBC0QV7ampHaZIFlrJ4v78E5JqH62JJzSLeh6OKBkBLT3SItOmqB/chedZezpp5/2T2Sg2zWA2IuGmYJGgPieI2nSpIn9ULPsGy6n9ehQOi061bjvwX2zwvtXt27d/M8V+dbrELvcxd4Lod+d+xd23/763bmnkf7zzz+le/fuZvjgnDlzpG7dur7d/c/D5NcO/04lvJBQ72LJXPlq0b/V+keMpCbXFf1468g333zTzEj36quvmskudKIPX9Ewqz2HOlxTe/S0p1GnA9d3JfnuP30+6frrrzcz4Wm41fCqMyr26dNHatcu2sQRvu/P/ZOAlFuEzwgggAACCCAQdgENScln3BX27+ELiieg/zqvPUT6L/T6XFHuosPUdHIEfYZHJ00oTtFZ47Ro78IFF1wQUJW+FFSLnodvmmidLc7eo6XbNcjoOecORNoDZP8lWnugtE7f0Cz9bg1/OjOeztTnK4cPHzbt1peb+op+h/aE/Pjjj+YXeHtPle5zonb46ijpn3HW8NXEZv+UrF/GB/+yWFe8uMqkSUL9fxTrtF988UVzfXQyjNxFh81pqNbn1K677joTkrRHTmcQ1MDrKxqe9Lk1HXapw+50BkT7e5J8+xX3Jw9zFFeQ4xFAAAEEEEAAgRgV0BnkdEa6Xr16mV9g9ZdY+x/9ZVZLKN6JpNN569AqnT48d9FhWdrjc84554j2zOiwtilTpuTp2dHnhHQYlk77bC86A5+96ORPOjW4L9TpVNMahnxBzLevftbnXnz76XqdIELru//++82zTL59fT9P1A7fPqXxM+nUGyX+pNZmSGvhv9+KCnFJUubc5633zBavX+X7778310sDZu4/vh5HnX1w69atsnLlSjMl/L///e+AU61WrZoJ43pP6oyJn3/+uVSpYk34EuJCQAoxKNUhgAACCCCAAAKxIqBDnLRoQMqv6PNCOrObToWtw86KU3SYmvYufPXVV2aIlb7aQ3t+NIjocCydMEF/GdaeA31Pkg6105nt9JdpnbVO33Ok03zrkCv7kDc9p9dee82890ifWdFgpcFOeyh0hjotGnp0+midjU4Dmtat4VDfq6STTGhPmRZt53//+18T0LTd2nth/6O9X7nbob0jkVBc1jTdKR2soZA1z7VmkCxE2LGG1blSKkuZC8ZIXPljQwhLoi116tQ54ZBMvQ8qVKgQtlMphE7YvpuKEUAAAQQQQAABBCJUQGcH++yzz8wLOHNPuuA7ZX2uSN8xNHr0aBMe7M/1+PYJ5qe+RFaHyOnMtfpHiz4DpAFGp4b2FR16pc8g6TMoOj23Fg1OOuutznaXu2gPl07ooMFKiw63m2tNEa4zo2nRUKOBSIOTvivJV7S3Smdb872fRyer0KLvhcpvenOd8a5Dhw7mZbi+dqidzq4XCcUVnyJlzhkmWb/PkMwfrSFvGXuOhiV35tHTs4bTWS8js8Yq5khCgysl+bSbxJUUviASCSb5nYPL6uLy5Lchktdp96mOY1y0aJF5+24knyvnhgACCCCAAAIIIBCcgPbg6DNPOvTq9NNPD5jW2V6TvmNJZ6LTZ4d0Jj3t1Tle0V95dV/tKWrWrFm+PRT63h2dHGDz5s3m3Uz6zJNvgojj1Xui9doODUr2559OtH9JbvNY7zVy//GT5OxcLO707eLJPmzNJFld4io3k4Qa7cWl705yaCEgOfTC02wEEEAAAQQQQAABBBDIK8AzSHlNWIMAAggggAACCCCAAAIOFSAgOfTC02wEEEAAAQQQQAABBBDIK0BAymvCGgQQQAABBBBAAAEEEHCoAAHJoReeZiOAAAIIIIAAAggggEBeAQJSXhPWIIAAAggggAACCCCAgEMFCEgOvfA0GwEEEEAAAQQQQAABBPIKEJDymrAGAQQQQAABBBBAAAEEHCqQ4NB202wEEEAAAQQQQAABBBwt4HFni/vg7+I5vEPEnSOu5IoSV6G+9ZLYco52ISA5+vLTeAQQQAABBBBAAAGnCeTsWyOZv7wpOVvnWcEoU8RlRQKXy1rOMhRxlZpKYqMeklDvYmu18wacEZCc9r8I2osAAggggAACCCDgSAFP1iHJWPqUZG+aeTQQedxHHTzZIp5jJO59q81+mT+Pk5SzB0t8lebHNjpgyXmR0AEXlSYigAACCCCAAAIIIGAXcB/aKukze0v2lrnWaisN+cKRfSf7stWb5Dm0RQ7P6StZv0+3b4n5ZQJSzF9iGogAAggggAACCCDgZAFP1p9y+Kv+4knf7h9GVzgPDVI5krF4sGRvX1C4Q2JgLwJSDFxEmoAAAggggAACCCCAwPEEjlgBx3Nktwk7x9vnhOut3qYj8/8l7iN/nHC3WNnIM0ixciVpBwIIIIAAAggggAACuQRydv9oTcbwrbXW+7xRru2F/mj1JGWufFVS2g4q9CH57bhixQqZOnWqNGzYULp27Sply5YN2G316tXyxRdfSI0aNeSSSy6RihUrBmzft2+ffPbZZ7J3715zfNOmTQO2h+IDPUihUKQOBBBAAAEEEEAAAQQiUCBz1RjrrGwzMBT1HK1nkrLXf16sXqRrr71WOnfuLBqCBg4cKB07dpQ//jjWKzVs2DBp3ry5LFq0SEaOHCkdOnSQnTt3+s945cqVJjiNGjVK5s+fL23atJFp06b5t4dqgYAUKknqQQABBBBAAAEEEEAgggR01rqcnUusMwpBQNJ2xSUenRq8CG387rvv5IMPPpCFCxfK66+/LtqTtHHjRnnzzTdNbWvWrJHBgwfL7Nmz5b333hPdv0yZMvLcc8/5v+2mm26Svn37mgA1adIkGTRokAwYMEA8nhC1z/tNBCQ/OQsIIIAAAggggAACCMSOQM7eX6zpvEP46771zqSc3cuLBFSrVi0zdO6UU04xxycmJkqVKlVk+3Zr4girTJ8+XRo0aGB6lfSzbu/du7dMnDhRP5r9tGepX79+1ruZrHc2WeXmm2+WdevWmcBkVoToL55BChEk1SCAAAIIIIAAAgggEEkCnvRtVkCKL/rkDPk0xv3npnzWFrxKg5EvHC1fvlzGjx9vhtfdcMMN5uD169eb55LsNWlg2rJli7jdbtmwYYPZpM8u+Yo+p6S9TJs2bZKzzz7bt7rYPwlIxSakAgQQQAABBBBAAAEEIlDAbb0A1tvbErKzy8koVlUaePTZowMHDsgtt9wiTZo0MfXpcLuqVasG1F25cmXJycmR3bt3m4CUmpoqKSkpefbZsWNHwLrifghhn1txT4XjEUAAAQQQQAABBBBAIFQCrqTyVu9RMWevy30yyZVyrwnqc+3atWX//v2ybNky8zzSlVdeaY5PSkqSrKysgLoyMzPN5/Lly0t+23WjHqPbQ1kISKHUpC4EEEAAAQQQQAABBCJEwFW+vjW7t9WLFKriSpD4SqGZVvuMM86Qe++9V2bMmGF6k/QZpT179gScqX6uXr26GUan2zUMHTx4MM8+9etb7QxhISCFEJOqEEAAAQQQQAABBBCIFIG4ClZw0F6kkBW3xFdvV6Tann/+ebnooosCjtWeJB1CFxcXJy1atJAlS5ZIdvaxQLdgwQL/c0mNGzc2vUi6zlcWL15sjtdnlUJZCEih1KQuBBBAAAEEEEAAAQQiREBne0tscLk1PXdSSM7IlVxZ4tPOLFJdOpRu7ty5MnbsWBOC5s2bJ/o+I11frlw56dmzp6l3+PDhZlKGn376ScaNG2em8tYN+nzS9ddfL0OGDDFD9NLT0+Wxxx6TPn36iA7bC2UhIIVSk7oQQAABBBBAAAEEEIgggaQmvayAZM1kV9xiDa9LOv1Oa86HotVVt25deeGFF+See+6RChUqyPnnny/nnHOOjBmjL7IVM4xu8uTJMnr0aKlUqZJ06dJF+vfvL926dfOfuYan5ORkSUtLEx1yp1OB29+T5N+xmAsu68VKoX2zUjFPqDCH33333fLiiy+aOc/btStaN19hvod9EEAAAQQQQAABBBCIdoGsjVMlY8nQok/3bb0gNq7q6VKm40v+dxAV1USfI9IZ62rWrClly5bNtxqdtlt7hXToXX5Fn01KSEgwQSu/7cVdxzTfxRXkeAQQQAABBBBAAAEEIlggsd7F4t6/TrJ+nWSFpGPP+BTqlK1w5CqTJmXaDy92ONLv016fRo0anfCr69Spc8Lt+oLZcJb8Y1k4v5G6EUAAAQQQQAABBBBAoEQFkk8fIEkt7rDei6S//hcyAmjPUZXmknrhOHEllivR8y3NL6MHqTT1+W4EEEAAAQQQQAABBEpIIKnpdRJfrbVkLHtO3Ht+snJSojUN+NF3DflPQQOU/olPleSWd0hC/UuL/NyRv84oWyAgRdkF43QRQAABBBBAAAEEECiqQHyVUyX1gjGSs2eVZG+ZKzm7fhDP4R1Hn09KqijxlZtJfI32klDrPHHFJxf1a6L6OAJSVF8+Th4BBBBAAAEEEEAAgeAF4qucJvqHklegkAMQ8x7IGgQQQAABBBBAAAEEEEAg1gQISLF2RWkPAggggAACCCCAAAIIFFmAgFRkOg5EAAEEEEAAAQQQQACBWBMgIMXaFaU9CCCAAAIIIIAAAgggUGQBAlKR6TgQAQQQQAABBBBAAAEEYk2AWexi7YrSHgQQQAABBBBAAAEEChBwp++UnB0LJWffGvGkb7Peh5QtruTKElexgcSf1Nq8INblchVQS2xuJiDF5nWlVQgggAACCCCAAAII5BHQQJSx/AVx7/reelGs9Z4jd5a1j/vYfnFJ1juRNCxVkaQW/SShXjfrRbHOGnTmrNYeu/QsIYAAAggggAACCCDgGAGPFYSOLPmvHJ7ZR9y7fzjabneG9dMWjnStO9MKSG7xHNktGUuflvTpPcV9cOPR/R3yNwHJIReaZiKAAAIIIIAAAgg4U8CTuV8Oz+kr2b9PtQA8JgAVSsIKVZ5DmyXdClXZOxYV6pBY2ImAFAtXkTYggAACCCCAAAIIIJCPgMd6tujwN/eJe9+v3uF0+ex0olVWb5LkHJEj395vPa9k1eGAQkBywEWmiQgggAACCCCAAALOFMhYPtIKR2vMc0XFEnDnyJF5d1uPJ6UXq5poOJiAFA1XiXNEAAEEEEAAAQQQQCBIAX12KHvdR0XrOcrzXdZzSVl/SubqCXm2xNoKZrGLtStKexBAAAEEEEAAAQQQsAQyVrwsojPQ6TC5UBRrAoes1W9LUuOe4kqqUKwaZ82aJbt27ZKePXsG1LN69Wr54osvpEaNGnLJJZdIxYoVA7bv27dPPvvsM9m7d6907dpVmjZtGrA9FB/oQQqFInUggAACCCCAAAIIIBBBAtrbk7NtXvGH1uVpk0uyN8/OszaYFb///rtceeWV8u677wYcNmzYMGnevLksWrRIRo4cKR06dJCdO3f691m5cqUJTqNGjZL58+dLmzZtZNq0af7toVogIIVKknoQQAABBBBAAAEEEIgQgZxdS63eo/jQn43Vi5S97Zsi1+t2u+WGG26w3q0U+BLaNWvWyODBg2X27Nny3nvvyXfffSdlypSR5557zv9dN910k/Tt29cEqEmTJsmgQYNkwIAB4vFYM/OFsBCQQohJVQgggAACCCCAAAIIRIKAe/866zQCQ0iozstM+lDEyp555hkTjnr06BFQw/Tp06VBgwbSsWNHsz4xMVF69+4tEydONJ+3b99uglG/fv384ermm2+WdevWmfUBlRXzAwGpmIAcjgACCCCAAAIIIIBApAl4juwJ0eQMeVum71UqSlm6dKk8++yz8uabb0pcXGAMWb9+vTRs2DCgWg1MW7ZsEe112rBhg9lm30efU9Jepk2bNgUcV9wPgWdW3No4HgEEEEAAAQQQQAABBCJAQCdmCO3QM3+jijCk7fDhw9KrVy8TkOrVq+evyrewceNGqVq1qu+j+Vm5cmXJycmR3bt3m4CUmpoqKSkpefbZsWNHwLrifiAgFVeQ4xFAAAEEEEAAAQQQiDSBJGv2t7jEsJyVK7Fc0PU++OCDcuqpp0qfPn3yPTYpKUmysrICtmVmZprP5cuXl/y260Y9RreHshCQQqlJXQgggAACCCCAAAIIRIBAXIX61lmE5xkkV8UGQbVQZ60bPXq0mZGue/fuon+mTp1qnh3SZe0hqlWrluzZYw0LtBX9XL16dTOMTrdrGDp48KBtDzHH1K+vbQ1d4T1IobOkJgQQQAABBBBAAAEEIkIgvlob6xmkoz0wIT2huCRJqHFOUFWWK1dO/vOf/wQco5MuHDhwQM466yzTO9SiRQt56623JDs7WxISjkaUBQsW+J9Laty4sdlP13Xp0sXUtXjxYjMET59VCmUhIIVSk7oQQAABBBBAAAEEEIgAgbiUqhJ3Uitx7/7ROpsQvShW2+XOkoSTjwaUwjazSpUq8uijjwbsrpMv6B/fen1h7EMPPSTDhw+XgQMHyqpVq2TcuHEyfvx4c5w+n3T99dfLkCFDTKjSWe4ee+wxM2Svdu3aAXUX9wND7IoryPEIIIAAAggggAACCESgQHKLO0J7VtYzTQn1L5W41LTQ1mvVprPRTZ482QzFq1Spkukl6t+/v3Tr1s3/XRqekpOTJS0tzQzJ05Bkf0+Sf8diLtCDVExADkcAAQQQQAABBBBAIBIF4k86XeJrdZCc7QtCM+W3K06SmvcLSVNfeeWVPPV06tRJtm7daqbt1l6h3FOBV6tWTWbOnGmeO9JheBUqVMhTRyhW0IMUCkXqQAABBBBAAAEEEEAgAgVS2j0hLmu4nVjhpnglTlLOeUp06F64S506dfKEI/t36pC9cIUj/Z7iStnPlWUEEEAAAQQQQAABBBCIIAGdkrtMx1HiSq5shaQiDh6zwlVy24eDnpwhghiCOhUCUlBc7IwAAggggAACCCCAQHQJxJWrI6ld3pG4ys2CC0nWjHWSYL2c9dyRklj/H9HV6GKcbRFjZDG+kUMRQAABBBBAAAEEEECgRAVcyZWkzF/HSPbmLyVz+Yviydhrfb/1niRPdq7zsNbpC2Y9OZLY8EpJOvVGcSWF51mfXF8cMR8JSBFzKTgRBBBAAAEEEEAAAQTCJ+ByuSSxzkXWNN2dJWfnEjN5Q86en8RzeJcJRJJYXuIrNpL46mdKQq3zHReMfPIEJJ8EPxFAAAEEEEAAAQQQcICAy3qmKKH6WeaPA5obdBN5BiloMg5AAAEEEEAAAQQQQACBWBUgIMXqlaVdCCCAAAIIIIAAAgggELQAASloMg5AAAEEEEAAAQQQQACBWBUgIMXqlaVdCCCAAAIIIIAAAgggELQAASloMg5AAAEEEEAAAQQQQACBWBUgIMXqlaVdCCCAAAIIIIAAAgggELQAASloMg5AAAEEEEAAAQQQQACBWBWIyoDkdrtj9XrQLgQQQAABBBBAAAEEEChFgagMSBMmTDBks2fPLkU6vhoBBBBAAAEE/r+9OwG3dS77B37vMw8ZO6aEOmYZjkyZ8q8MSRGVKUO6KkQKKZFSvVzlJbO8xUslV9KbpEEoEhlDSDKEYziGU4hznHGv/7of797vPvsMlu051vOs/Xmu6zhrr/Ws37p/n1s63/N7nt8iQIAAgU4TqGVAWnPNNYs+TJgwodP6YT4ECBAgQIAAAQIECLRRoJYBaaONNirIllxyyTbS+WgCBAgQIECAAAECBDpNoJYBqdOaYD4ECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQQEpAo0QQkECBAgQIAAAQIECFRDQECqRh9UQYAAAQIECBAgQIBABQSGVaAGJRAgQIAAAQIECHSQwMMPPxz/+te/YsiQITFhwoT5zuyxxx6Lp59+unj97W9/+3zP8wKB11PACtLrqe2zCBAgQIAAAQKDQODLX/5ybLDBBrH++uvH3//+9/nOeI899ijOy3MbjcZ8z5vfC/fdd1/84he/iClTpszvFM8TeNUCAtKrJvMGAgQIECBAoAyBRqM7GrOmlTGUMSoscNFFF82zuscffzyuv/76eb7W6pM59k477RSTJk1q9S3OI/CKAgLSKxI5gQABAgQIEChToHvqU/HSn74YU362ZUz5+f+LKZe9L2Y88NMyP8JYFRHYeOONY34B6eKLLy5WjTbccMOKVKsMAi8LuAfJvwkECBAgQIDA6ybQ/dLTMfWqvSNmNi+JaswuPrcx/V8x485TovuFR2LU+oe/brX4oIUv8KEPfSi++MUvxt133x1rr732HB/4k5/8JDbaaKNYaaWV4tZbb53jtfzhxhtvjP/+7/+Ou+66K1ZeeeXYdtttY++9946urq7i3K9//evx85//vHicn/Gud70rDj744OLne++9Ny677LK44oorivO32GKLeM973hObb7558fqzzz4bhx12WPF4Xv844IADYpNNNul96Yc//GH89Kc/jb/97W9FvR/4wAfioIMOiqFDh/aeM3369PjlL38ZV155Zdx2222x7rrrRn7unnvuGSNGjOg9z4PqC1hBqn6PVEiAAAECBDpGYMadZzbD0dTecNQ7se5ZMevBn8bs5+7vfcqD+gvsvPPOxUYN/VeRJk6cWASgXXfdNWbPfjko953t2WefHVtuuWX84Q9/iNVXXz3+8Y9/xL777luEjRkzZhSnPv/88zF1avPfpebx3HPPxQsvvFA8vvnmm4t7n/I+qMUWWyyWWWaZOOecc4qw0lNHjpGhrP+vCy64IM4///x46qmnirHyH5/61Kdin332KT5j++23j+7u7vjc5z4XW2+99Ry1ZxD68Ic/HL/73e8iV84efPDB+OQnPxm//vWve8fyoB4CVpDq0SdVEiBAoOMFzjrrrDjppJOKv3XNna8cnSlw9oefiNHD53Mz/pDhMXvS9TF08VU7c/KDcFZLLbVUbLXVVpGrRd/4xjd6BfLyujwyIF133XW9z+eDhx56qAgg7373u4sVmeHDhxevn3jiiXHEEUfEdtttFx/72MeK/14svvji8ZWvfCX+67/+K1ZZZZXivHPPPbf4/ZFHHolll122ePzoo48Wr//oRz+K3XbbrQhNuTLV98j3feITn4hddtklcoUoj6uuuiq+973vFfWcfPLJvafnf6s+//nPxxlnnBGf/exni3ugfvazn8113qc//eliN7/eN3pQCwEBqRZtUiQBAgQ6XyD/1jb/ljh/OTpX4NxdN25O7uVLpOaaZbFpw8srAnO95onaCmQgyUvWbr/99mJlJyeSgSkvYVtxxRXn2r3u9NNPj7xcLVdpesJRvueQQw6JY489Ni688MIiIOVz8zpytSd3x+sJR3nOEkssEePGjZtvWMmQlmEmtyT/wQ9+0HsZX9aSR9bf98hajj766GIeGZDyErpLLrkkNt10076nxaqrCvtzgNTkBwGpJo1SJgECBDpdILcDvuWWW4q/sV1zzTU7fbqDdn5dd3wuYurEec+/eW/JEKtH87ap8bN5H1LeG5SXt+X/zvM7kvIyuL4rMn2nd8899xQBJe8h+tWvftX3pXjDG94QuTK0oCPvM8ptv3NFKD8nP++OO+4ovm/prW9961xvzfFy1ShXoy699NIYO3Zs7zm5jXgGrbzMr++RwS3vqerZwvyNb3xjfPCDHyw+59RTTy0ur8u/7MlL+HpWtPq+3+NqCwhI1e6P6ggQIDDoBJZeeul405veNOjmPVgmPKtxUEy78ai570HqGhJdo8bFsOX/32ChGDTzzJWbvFwuV42++c1vFr/nRgsf+chH5mnw5JNPxrBhwyLvU+p/5I53GWQWdJx55pnF6k6uQm222WaxxhprFPcQHXPMMXO97cUXX4wdd9wx8n6mq6++uljR6ntS3ovUf3OJntdzDj3f3ZRfdpsB6YYbbig2lMg68zK9HXbYoed0v9dIQECqUbOUSoAAAQIE6i4wbPmtYsQ6B8eMu86IGNL8Y0jx5aDd0TV6mRi95anR1bwPydF5Arvvvnt8/OMfL1Z0MijlKs/yyy8/z4nmKs9f/vKX+M53vhMrrLDCHOfMa0OHvic888wzxT1B+RctGbAWXXTR3pe/8IUv9D7OBxlucle8O++8M84777wiTM1xQvOHnh32MmyNHDmy9+WXXnqpqHG99dYrnjvuuOOKcHTCCScU90n1nJj3LznqJ+Au2Pr1TMUECBAgQKDWAiNW2yPGbH9JjFzv0Bix1idi1KbfijHbXRRD3vDmWs9L8fMXyN3s8j6d448/Pv785z8XGyXM7+ye7bX7X1732GOPxSKLLFLsKtf/vT0rOXlZW4aovAepbzjKEPTEE0/M8bbc5S63CT/88MPne09TbkOeYShXhvoe+fPMmTOLDSPy+bwULzeXyV3r+h6XX3553x89romAgFSTRimTAAECBAh0ksCQMUvH8PEfjBFr7BPDltu8uXLkopZO6m//ueRlcfk9RnmPTwaJ3A57fkdugJCbN+T3HOX5jz/+eFxzzTXFpXB5WVvuHtdz5Hl55Dbe+V1Lea/Q6NGj48c//nFxP+MDDzxQfJdSfnaGq8mTJxfj5YYKGdayrtwG/Nvf/vYcv66//vpi3AxRo0aNKna3y+278/K//G6l/fbbr9hC/NBDDy3Oy80dcvvv3FEvQ1perrfXXnsVW34XJ/hHrQT816hW7VIsAQIECBAgQKCeArmbXX6R6jvf+c45dpjrP5sxY8YUX7b60Y9+tLivp+f1DTbYIH7zm99E301ccsvvvN8nL227//77I7fazsva8ruLttlmm+KtGYByG/AMWrkzXj6f9x3lkd+f1P/Su3w+d6jLywAzgGUgyu84yu896jlWW221opbcnCGPo446qtiuPHe9y1/5BbK5apahLTeKcNRLoKu5JDmfLyOo7kRyO8XTTjut+Bculz4dBAgQIFB/gf333z+++93vFvcDrLPOOvWfkBkQIPCaBHJFJkNP7hSXO8nln/lyBWleR4afXA3q2YEuV4py57q8h6nvDnSTJk0qLsF785tf3eWc+cWyubteXuaXAW38+PFz1ZJ/pM7PzNCVteaOe456ClhBqmffVE2AAAECBAgQ6GiBvBQvw03fgDO/Cfff8CF3zuu74tPzvuWWW67n4av6Pe+fysvo8tf8jgxvuY25o/4C7kGqfw/NgAABAgQIECBAgACBkgQEpJIgDUOAAAECBAgQIECAQP0FBKT699AMCBAgQIAAAQIECBAoSUBAKgnSMAQIECBAgAABAgQI1F9AQKp/D82AAAECBAgQIECAAIGSBASkkiANQ4AAAQIECBAgQIBA/QUEpPr30AwIECBAgAABAgQIEChJQEAqCdIwBAgQIECAAAECBAjUX0BAqn8PzYAAAQIECBAgQIAAgZIEBKSSIA1DgAABAgQIECBAgED9BQSk+vfQDAgQIECAAAECBAgQKElAQCoJ0jAECBAgQIAAAQIECNRfQECqfw/NgAABAgQIECBAgACBkgQEpJIgDUOAAAECBAgQIECAQP0FBKT699AMCBAgQIAAAQIECBAoSUBAKgnSMAQIECBAgAABAgQI1F9AQKp/D82AAAECBAgQIECAAIGSBASkkiANQ4AAAQIECBAgQIBA/QUEpPr30AwIECBAgAABAgQIEChJQEAqCdIwBAgQIECAAAECBAjUX0BAqn8PzYAAAQIECBAgQIAAgZIEBKSSIA1DgAABAgQIECBAgED9BQSk+vfQDAgQIECAAAECBAgQKElAQCoJ0jAECBAgQIAAAQIECNRfQECqfw/NgAABAgQIECBAgACBkgQEpJIgDUOAAAECBAgQIECAQP0FBKT699AMCBAgQIAAAQIECBAoSUBAKgnSMAQIECBAgAABAgQI1F9AQKp/D82AAAECBAgQIECAAIGSBASkkiANQ4AAAQIECBAgQIBA/QUEpPr30AwIECBAgAABAgQIEChJQEAqCdIwBAgQIECAAAECBAjUX0BAqn8PzYAAAQIECBAgQIAAgZIEBKSSIA1DgAABAgQIECBAgED9BQSk+vfQDAgQIECAAAECBAgQKElAQCoJ0jAECBAgQIAAAQIECNRfQECqfw/NgAABAgQIECBAgACBkgQEpJIgDUOAAAECBAgQIECAQP0FBKT699AMCBAgQIAAAQIECBAoSUBAKgnSMAQIECBAgAABAgQI1F9AQKp/D82AAAECBAgQIECAAIGSBASkkiANQ4AAAQIECBAgQIBA/QUEpPr30AwIECBAgAABAgQIEChJQEAqCdIwBAgQIECAAAECBAjUX0BAqn8PzYAAAQIECBAgQIAAgZIEBKSSIA1DgAABAgQIECBAgED9BQSk+vfQDAgQIECAAAECBAgQKElgWEnjGIYAAQKlCnzyk5+M8847LxZZZJEYMsTf5ZSKW9HB/v3vfxeVPfHEE7HOOutUtEplESBAgECnCwhInd5h8yNQU4GJEyfG7NmzY9FFF41Ro0bVdBbKfjUC06ZNi1mzZsWIESNezducS4AAAQIEShUQkErlNBgBAmUJjB8/vhjq0ksvjQkTJpQ1rHEqLLD//vvHd7/73Rg3blyFq1QaAQIECHS6gOtWOr3D5keAAAECBAgQIECAQMsCAlLLVE4kQIAAAQIECBAgQKDTBQSkTu+w+REgQIAAAQIECBAg0LKAgNQylRMJECBAgAABAgQIEOh0AQGp0ztsfgQIECBAgAABAgQItCwgILVM5UQCBAgQIECAAAECBDpdQEDq9A6bHwECBAgQIECAAAECLQsISC1TOZEAAQIECBAgQIAAgU4XEJA6vcPmR4AAAQIECBAgQIBAywICUstUTiRAgAABAgQIECBAoNMFBKRO77D5ESBAgAABAgQIECDQsoCA1DKVEwkQIECAAAECBAgQ6HQBAanTO2x+BAgQIECAAAECBAi0LCAgtUzlRAIECBAgQIAAAQIEOl1AQOr0DpsfAQIECBAgQIAAAQItCwhILVM5kQABAgQIECBAgACBThcQkDq9w+ZHgAABAgQIECBAgEDLAgJSy1ROJECAAAECBAgQIECg0wUEpE7vsPkRIECAAAECBAgQINCygIDUMpUTCRAgQIAAAQIECBDodAEBqdM7bH4ECBAgQIAAAQIECLQsICC1TOVEAgQIECBAgAABAgQ6XUBA6vQOmx8BAgQIECBAgAABAi0LCEgtUzmRAAECBAgQIECAAIFOFxCQOr3D5keAAAECBAgQIECAQMsCAlLLVE4kQIAAAQIECBAgQKDTBQSkTu+w+REgQIAAAQIECBAg0LKAgNQylRMJECBAgAABAgQIEOh0AQGp0ztsfgQIECBAgAABAgQItCwgILVM5UQCBAgQIECAAAECBDpdQEDq9A6bHwECBAgQIECAAAECLQsISC1TOZEAAQIECBAgQIAAgU4XEJA6vcPmR4AAAQIECBAgQIBAywICUstUTiRAgAABAgQIECBAoNMFBKRO77D5ESBAgAABAgQIECDQsoCA1DKVEwkQIECAAAECBAgQ6HQBAanTO2x+BAgQIECAAAECBAi0LCAgtUzlRAIECBAgQIAAAQIEOl1AQOr0DpsfAQIECBAgQIAAAQItCwhILVM5kQABAgQIECBAgACBThcQkDq9w+ZHgAABAgQIECBAgEDLAgJSy1ROJECAAAECBAgQIECg0wUEpE7vsPkRIECAAAECBAgQINCygIDUMpUTCRAgQIAAAQIECBDodAEBqdM7bH4ECBAgQIAAAQIECLQsICC1TOVEAgQIECBAgAABAgQ6XUBA6vQOmx8BAgQIECBAgAABAi0LCEgtUzmRAAECBAgQIECAAIFOFxjW6RM0PwIECBAgQIAAgVcvMHHixLj55pvjrrvuihVXXDHWW2+9eNvb3hajR49+9YP97zv+9re/RaPRiLXWWmvAY+Qbb7vttujq6or111//NY3T8+Znn302HnrooVhjjTVizJgxPU/7fZAKCEiDtPGmTYAAAQIECBCYl8DkyZPjgx/8YFx//fVzvbzIIovE97///dh5553neq2VJ/bYY4+YNWtW3H333a2cXpxz3333xb333hvvec97YuzYscVzH//4x2Po0KHx5z//ueVxFnTi5ZdfHnvuuWfccsstseGGGy7oVK8NAgGX2A2CJpsiAQIECBAgQKAVgQceeCA23XTTIhzttdde8bvf/S6eeuqp+Mtf/hJnnnlmEUo+9KEPxde+9rVWhpvrnPe///2x4447zvX8gp646KKLYqeddopJkyYt6DSvEShNwApSaZQGIkCAAAECBF6NQPeUSdGY9s8Ysvgq0TV01Kt5q3MXksCnPvWpyJB01llnxYEHHtj7KUsvvXSsu+66scMOO8Rmm21WBKRddtkl1llnnd5zWnnwH//xH62c5hwCbRWwgtRWfh9OgAABAgQGn0Cj0R3Tbv5aTP3NzvHS1Z+MKZftELOevnXwQVRsxjfeeGNcffXVsfnmm88RjvqWudJKK8U3vvGN4j6ib37zm8VL+Z7PfOYz8cwzz8S5554bGZxyjC996Uvx4osv9n175HuOP/743ue6u7uL92yzzTbF/U0f+MAH4jvf+U7v61//+tfjkksuKX7+4he/GGeccUbva/kg75P68pe/HFtssUXk6lTf9/Y98Yc//GGxCrXaaqtFftZpp50Ws2fP7nvKXI9feumlOPTQQ2P//fePp59+uvf1XMn69Kc/HRtvvHERGvfdd9/4+9//3vt6PvjqV78aJ510Ujz33HNx663+3Z4DpwY/WEGqQZOUSIAAgcEgkH9QyiPvN8ibuB2dK7DIc3+McU9fGV3FFJu9njUlpt1wZIx9/6+aK0kjO3fiFZ/ZL3/5y6LCvitH8yp5n332iQMOOCAyUOWRl99lcPnrX/8ad955ZxFA8n/DJ5xwQlx66aXFRg9veMMbinN//OMfF/cgHXXUUcXPRxxxRJxyyinF/UXvfe9745577inCR46VYz7//PMxderU4twMGy+88ELxOP+RgWzLLbeMKVOmFJcF/va3v41f/epXxc+f//zne8/LVbHvfe978c53vjO233774v6nz33uc0Xwuuqqq4rLBntP/t8H+ZkZuP70pz/FxRdfHLmClscdd9wR2267bfH4fe97X7FRxBVXXBE/+clP4n/+538in8sj550BLD83A6H7mgqW+vyj+S9w7Y5DDjkk/5+z0dxZpXa1K5gAgdYEmv/nW/zv/Pbbb2/tDc6qvcDb3/72ouf533e/Otvg+4eu0vj3TzZuvHDxJnP8mvXPe2r/73GdJ9DcQKH4314z+LziNJorMY1hw4Y1mhsuNE4++eTifW9605sa//znP3vf2/P8Mccc0/tccye8RnMnvOLn6dOnN0aMGNForjj1vp4PDjrooEZzM4ZGM/gUzzdXkYrx77///t7zcpz870QzYDVynDyaq0mNZZddttFc2ek978orryzOawai3ufywYknnlg83wxnxfMXXnhh8XNzk4ZGc9WrsdVWWzVGjRrV+M1vfjPH+5orY40ll1yy8eCDD/Y+37xHq7Hyyis33vzmNzeaAa54vqe+gw8+uNH8y5/ecz2oh4AVpOb/uhwECBAg0H6BL3zhC9H8Q0uxbe+QIa4Ab39HFl4FSy7/eMxuTI5hXfln0v87uka/8f9+8Oh1F3jyySeLz1xhhRVe8bMXW2yxYiVoxowZvefmJXXN8ND7c67S5CpSXiKXl8r1P/ISttzR7qabbor87Ga4KU7Jy/A++9nPRjOA9X/LHD+PHDky8p6mZsgqns+6m3/REtdcc02xCp3bgJ9++unFa7ni1fdo/mV7HH300cXKT35Wz5ErVLmSlbvjXXbZZbH11lv3vBTNv7ArNq/YfffdY/z48b3P5+pSrqrlZXU33HBDsYKWL+aOfzn/rMNRL4EF/5tXr7molgABAgRqLLDbbrtF/nJ0vkD31Cdj6hUfjehu/uG6e2bEkOExbKXtY8joly9j6nyBas6wuQJSFJa71jVXgxZY5MMPPxzLLbfcHN+JtOaaa871nrXXXrsIFc11g7mCQoas5mpREWIy3OQlcO9+97uLjSAmTJgw11j9n8jvZuoJRz2vjRs3rrgkL4NbBqi8ZDeD1+qrr95zSvH78OHDI2vrf+9QBp3HHnss8pLAVVdddY735Fh5NFeX5ng+f8hglkeOl/c45bH88svP4VM86R+1EPBXdLVokyIJECBAgEDnCAwZs2yM2fbCGL7KbjFshW1j5NuPbP76UudMsKYz6QkEef/Pgo5HHnmkuP9nlVVWmeO0noDV98nFF188ckV4fqsouVlCrrp84hOfiOZla8WGC/nlr7mKM23atL5DzfU4V2j6HxnE+h4Z9nrm1ff5fJw19T8/N5XIzRXyvqb+q045Vh650UP/o2d+fcfzhbP9lerzs4BUn16plAABAgQIdIzAkDHLxMh1D45Rm3w9hr9lh/n+AbpjJlyDifRsPpChZUFHz05yudrS93jooYf6/lg8zl3m+q/e9JyUYSIvsXvHO95R7D6Xq1IZzvJLaHPDhR/96Ec9pw7499x1L3eRa96nNMcYeXlfbi7R91K5PKF5L1IcdthhRWDLL4/N3e96jhwrj3l9gW6GvDz6j1c86R+1ExCQatcyBRMgQIAAAQIEyhfYZJNNiq2wmxsVFPf29F0N6fm05qYFxRfGZhD42Mc+1vN08XsGir7HvffeGzlWfn/SvI7cNS8vdfv5z3/e+/Jaa60VRx55ZPHz5MmTe5/PB/OqZ44T5vHDRhttFBmGegJMzyn588yZM2O77bbrear4famllip+/9a3vhX5OLf57tniO1e28r6o/PLc/se1115b3HOUX7LrqL+AgFT/HpoBAQIECBAgQKAUgby8LC9Ja+48V4SHU089NZo7wcUFF1xQ3C+00047FVte55bW/TdRyA0RcoOF3Ko737PjjjsWoeFrX/vaPGvLLbpzU4evfOUrxUYNeVlbbo6QmycMHTq0916evNcoj3POOafYonueg83nyfyOpOZudMWKUAab3Awit+Xeb7/9Iu+BygA0r2OJJZYoLrVr7soXuaFDHllHbhn+hz/8ofiOp/xC3bzn6PDDD48//vGPxXl9N6mY17ieq4eATRrq0SdVEiBAgAABAgQWukBzu+q47bbbork9dXG5WQadniM3RMgNCM4///xidaXn+Z7fcxfK4447rggP+VyuwOQlavPbFS/vT8rglWElL7PrOfLeorPOOqt344Nc5cnvEcod4ZpbfcfPfvaznlNf8fcMNRmI9txzzzl2pMv7iHI17I1vnP/OiXvvvXecd955cdFFFxXvz8CXc8zvN8ogmL/yyPuPMojNLwi+YpFOqJxAV3O5cs672SpX4twF5XaMeX1s83uQIpdOHQQIdJ5AflHh2WefXWyr2spuRp0nYEYECBBor0DuBJerQbl7W27IkLu+9d81LivML3rNlZj8ktjmdxwV/93Oy9rykrTm9xm94iRyY4S77rqr2D1umWWWiXXWWSdyBaf/8fjjj0eGqlbG7P/enrnkDnW5215eItizsUL/c1v5OVei8h6m3O0uLyGc14YRrYzjnGoKWEGqZl9URYAAAQIECBBoq0CGofwLqlfzl1S5Y90GG2zwqurOkNHKvTu5bfZAj4HMZUGflVuH93xv04LO81o9BdyDVM++qZoAAQIECBAgQIAAgYUgICAtBFRDEiBAgAABAgQGi0BuTJC7z+VmCA4CnSAgIHVCF82BAAECBAgQINAmgfw+pPz+ovl9IWubyvKxBAYsICANmM4bCRAgQIAAAQIECBDoNAEBqdM6aj4ECBAgQIAAAQIECAxYQEAaMJ03EiBAgAABAgQIECDQaQICUqd11HwIECBAgAABAgQIEBiwgIA0YDpvJECAAAECBAgQIECg0wQEpE7rqPkQIECAAAECBAgQIDBgAQFpwHTeSIAAAQIECBAgQIBApwkISJ3WUfMhQIAAAQIECBAgQGDAAgLSgOm8kQABAgQIECBAgACBThMQkDqto+ZDgAABAgQIECBAgMCABQSkAdN5IwECBAgQIECAAAECnSYgIHVaR82HAAECBAgQIECAAIEBCwhIA6bzRgIECBAgQIAAAQIEOk1AQOq0jpoPAQIECBAgQIAAAQIDFhCQBkznjQQIECBAgAABAgQIdJrAsDpO6OGHHy7K3muvvWLs2LF1nIKaCRB4BYEHHnigOGPy5MmvcKaXCRAgQIAAAQLlCdQyIC211FKFQAalrq6u8jSMRIBAZQRmzJhR1DJlypTK1KQQAgQIECBAoPMFuhrNo/OnaYYECNRN4MADD4yzzz47br/99pgwYULdylcvAQIECBAgUFMB9yDVtHHKJkCAAAECBAgQIECgfAEBqXxTIxIgQIAAAQIECBAgUFMBAammjVM2AQIECBAgQIAAAQLlCwhI5ZsakQABAgQIECBAgACBmgoISDVtnLIJECBAgAABAgQIEChfQEAq39SIBAgQIECAAAECBAjUVEBAqmnjlE2AAAECBAgQIECAQPkCAlL5pkYkQIAAAQIECBAgQKCmAgJSTRunbAIECBAgQIAAAQIEyhcQkMo3NSIBAgQIcstKEAAAEpxJREFUECBAgAABAjUVEJBq2jhlEyBAgAABAgQIECBQvoCAVL6pEQkQIECAAAECBAgQqKmAgFTTximbAAECBAgQIECAAIHyBQSk8k2NSIAAAQIECBAgQIBATQUEpJo2TtkECBAgQIAAAQIECJQvICCVb2pEAgQIECBAgAABAgRqKiAg1bRxyiZAgAABAgQIECBAoHwBAal8UyMSIECAAAECBAgQIFBTAQGppo1TNgECBAgQIECAAAEC5QsISOWbGpEAAQIECBAgQIAAgZoKCEg1bZyyCRAgQIAAAQIECBAoX0BAKt/UiAQIECBAgAABAgQI1FRAQKpp45RNgAABAgQIECBAgED5AgJS+aZGJECAAAECBAgQIECgpgICUk0bp2wCBAgQIECAAAECBMoXEJDKNzUiAQIECBAgQIAAAQI1FRCQato4ZRMgQIAAAQIECBAgUL6AgFS+qREJECBAgAABAgQIEKipgIBU08YpmwABAgQIECBAgACB8gUEpPJNjUiAAAECBAgQIECAQE0FBKSaNk7ZBAgQIECAAAECBAiULyAglW9qRAIECBAgQIAAAQIEaiogINW0ccomQIAAAQIECBAgQKB8AQGpfFMjEiBAgAABAgQIECBQUwEBqaaNUzYBAgQIECBAgAABAuULCEjlmxqRAAECBAgQIECAAIGaCghINW2csgkQIECAAAECBAgQKF9AQCrf1IgECBAgQIAAAQIECNRUQECqaeOUTYAAAQIECBAgQIBA+QICUvmmRiRAgAABAgQIECBAoKYCAlJNG6dsAgQIECBAgAABAgTKFxCQyjc1IgECBAgQIECAAAECNRUQkGraOGUTIECAAAECBAgQIFC+gIBUvqkRCRAgQIAAAQIECBCoqYCAVNPGKZsAAQIECBAgQIAAgfIFBKTyTY1IgAABAgQIECBAgEBNBQSkmjZO2QQIECBAgAABAgQIlC8gIJVvakQCBAgQIECAAAECBGoqICDVtHHKJkCAAAECBAgQIECgfAEBqXxTIxIgQIAAAQIECBAgUFMBAammjVM2AQIECBAgQIAAAQLlCwhI5ZsakQABAgQIECBAgACBmgoISDVtnLIJECBAgAABAgQIEChfQEAq39SIBAgQIECAAAECBAjUVEBAqmnjlE2AAAECBAgQIECAQPkCAlL5pkYkQIAAAQIECBAgQKCmAgJSTRunbAIECBAgQIAAAQIEyhcQkMo3NSIBAgQIECBAgAABAjUVEJBq2jhlEyBAgAABAgQIECBQvoCAVL6pEQkQIECAAAECBAgQqKmAgFTTximbAAECBAgQIECAAIHyBQSk8k2NSIAAAQIECBAgQIBATQUEpJo2TtkECBAgQIAAAQIECJQvICCVb2pEAgQIECBAgAABAgRqKiAg1bRxyiZAgAABAgQIECBAoHwBAal8UyMSIECAAAECBAgQIFBTAQGppo1TNgECBAgQIECAAAEC5QsISOWbGpEAAQIECBAgQIAAgZoKCEg1bZyyCRAgQIAAAQIECBAoX0BAKt/UiAQIECBAgAABAgQI1FRAQKpp45RNgAABAgQIECBAgED5AsPKH9KIBAgQIECAAIF5CzRmTonZz/w5up+7L7qnPhkxe0Z0jVgsuhZZMYaOWy+GLr7avN/oWQIECLxOAgLS6wTtYwgQIECAwGAW6H7hkZhxz7kx67HfR3R1vUzRPfN/SZo/Dx3xclgavVQMX32fGD5+p+gaMnwwk5k7AQJtEnCJXZvgfSwBAgQIEBgMAo1GI6Y3g9HU3+4Rsx6/OqIxKyKDUW84SoVGMxxNL35vvPR0zLjz1Jh6+a4x+/kHBgORORIgUDEBAaliDVEOAQIECBDoFIHG7Gkx7U9fiJn3ntecUne/ULSAWTbDU2PqU/HS7/aLWU9cu4ATvUSAAIHyBQSk8k2NSIAAAQIECDQFpt10bMx+6sZmMGquGr3q4+VANe2GL8XsyXe86nd7AwECBAYq4B6kgcp5HwECBAgQIDBfgen3fj9mT7ru5Uvq5ntWCy80ZsdL1x0eY7b7cQxp3p/k6FyBO++8M2bNWnCYXn311WPs2LFtQ3j++efjwQcfjPHjx8fiiy/etjp88MIVEJAWrq/RCRAgQIDAoBPonvJEzPzrd5vhaHY5c++eEdPvOCVGb3pcOeMZpZIC73rXu+Jf//rXAmu7/vrrY7PNNlvgOQvzxauvvjp23nnnuPjii+PDH/7wwvwoY7dRQEBqI76PJkCAAAECnSgwo7l61NyCrsSANDNmP/776H5hYgxpbgfu6FyBcePGxbHHHjvfCebKTTuPFVZYIfbee+9YaaWV2lmGz17IAgLSQgY2PAECBAgQGEwCjeb3Gs2aeHnrGzK0ijNkWMx8+Jcxcp1Pt/oO59VQYLHFFouDDjqospVvsMEG8YMf/KCy9SmsHAGbNJTjaBQCBAgQIECgKTD7n3cNcFOGV+Br7mw364k/vMJJXh4MAjfddFPst99+cccdc27ekfcG5fPXXnttL8NXv/rVOOmkk+K5556LI488Mrbaaqt4//vfH6effnrvOT0PJk6cGIcddlhss802xbm33nprXHnllXHUUUf1nBJ5n9RnPvOZuT679wQPOkJAQOqINpoEAQIECBCohkD38/dHNFd7FsbRePHRaJR1X9PCKNCYr4vAP/7xjzj//PPj0UcfnePznn766eL5++67r/f5Sy+9tHjuHe94R/ziF7+I5ZZbLu6+++445JBD4vjjj+8977bbbotcHTrzzDNj9uzZcd1118Xmm28exxxzTPznf/5n73n52WeccUY88IDv6OpF6cAHC+e/YB0IZUoECLy+ApMmTSo+8NBDD7VT0OtL79MIvCaB3Sb8O9675vQYtjD+CrbRHY3pz0XXqDe+phq9uboCkydPLsLLvCrM1aH1119/Xi8t8LkMRAcffHCcfPLJMWzYsGI1aY011ohLLrmkd3UoL+ubPn16EZ5WXXXVYryLLroodt999+I9C/wAL3acgIDUcS01IQKdITBjxoxiItdcc01nTMgsCAwSgU0Wa26isMayzdl2LZwZd5e0M97Cqc6or1Egt9Ge1+VvOewWW2wxoIC0yCKLxAknnNAbdHJ77k022SSuuuqq5opko7hc7sYbb4yjjz46esJRft5uu+0Wp512Wtx88835o2MQCQhIg6jZpkqgTgK//vWv45FHHolRo0bVqWy1Ehj0AsMnXhRDJ17c3MFu5kKx6BqxyEIZ16DVEFh55ZVLv3xt+eWXj9GjR88xwdwtb+rUqZF/GXf//c3LQpvHWmutNcc5+cOECRMEpLlUOv8JAanze2yGBGorYBvV2rZO4YNYYOaMtWL6o83Vo8ZCQBixaHQNm/MPugvhUwxZU4FcfZrXMWbMmHk93ftcXtaXxzLLLNP7XM+DIUMWxrWiPaP7vaoCul7VzqiLAAECBAjUUGDouAnNXexevkS23PK7YuhSG5Q7pNFqLZCXx/U9+m7O0Pf5V3q84oovf7dWXmbX/+i/U17/1/3cmQICUmf21awIECBAgEBbBIaMHhdDlpj7UqXXXEzX0Bi+4vaveRgD1F+gJ9D03Ukuw1Jemj2QY+21145cKcotvfsed911V+RW347BJ+ASu8HXczMmQIAAAQILVWDEGvvGtBu/1LzMrru0z+lqBq+hb9q8tPEMVF+B3Mlu+PDhceKJJ8Zb3/rWyE0YLrzwwrjlllsGNKm3vOUtccABB8RZZ51V7KC30047xV//+tfi+5NGjhwZs2bNGtC43lRfAStI9e2dygkQIECAQCUFhi2/VQxZ8m3NjexK+nvYriExcsLh0dVcRXIQyHuKMhDll7/usssu8b73vS+eeuqpuOCCCwaM8+1vf7vYxS7H2HrrrYvvSNp3331j1113jUUXXXTA43pjPQW6mkuSc17AWc95qJoAAQIECBCokED3S8/E1Cv3ipjx72ZVr+GPGkOGx7Dxu8SoCYdWaHZKqYLAzJkz45577om85G6JJZYoraRnn322d7wMYBnEfv/735c2voGqL2AFqfo9UiEBAgQIEKidwJDRS8XoLU6JGNrcqn+gKz/NcDR06Y1i5Hqfrd38FbzwBfIyu/XWW683zAz0E6dMmRIbbrhhHHHEEcUQPWHrySefjGuvvbb4jIGO7X31FChp7buek1c1AQIECBAgsPAEhi65ZozZ+vvx0h8Pica0fzZ3t2v1u5Ga24R3dcXwlT8SI9Y9uPnQ3+cuvC4ZeezYsbHsssvGKaecEhmKtt9+++KSvXPOOae41+nAAw+ENMgEXGI3yBpuugQIECBA4PUWaMxqfiHnPefGzAcuan50M+zMbxvwIgh1RdfY5Yt7joYtu8nrXarPG6QCL7zwQuR9SFdddVXcdNNNsfTSS8fGG28cxx13XKy55pqDVGXwTltAGry9N3MCBAgQIPC6CjSmPxczH70iZj1+TXQ/d1/EzBdf/vxmMOoavUxxOd2wFbZp/r5hc9WouYrkINAGgby3KS/fcwxeAQFp8PbezAkQIECAQFsFGt3N7ZObv7qGNe9TchAgQKAiAgJSRRqhDAIECBAgQIAAAQIE2i/grsf290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsFBKT290AFBAgQIECAAAECBAhUREBAqkgjlEGAAAECBAgQIECAQPsF/j9itZSN0H2EPQAAAABJRU5ErkJggg==",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "## Tree representation of DFA kingdom/clade distribution (Figure 4A)\n",
    "\n",
    "tree_king=read.tree('data/step1/DFA_king_tree.nwk')  # import the newick phylogenetic tree file\n",
    "\n",
    "tree_tips=c(intersect(tree_king$tip.label,taxo_table$tip.label)) # make sure that data are consistent between table and tree tips\n",
    "\n",
    "tree_pruned=keep.tip(tree_king,tree_tips) # prune the tree to only keep tips of representative sequences\n",
    "\n",
    "p=ggtree(tree_pruned, branch.length=\"none\") %<+% taxo_table \n",
    "\n",
    "p_tree <- p + geom_tiplab(aes(label=Taxo),offset = 1.2, hjust = .5) +\n",
    "  geom_tippoint(aes(size = nb_hits),color ='#F7B846') + \n",
    "  theme(legend.position = \"right\")\n",
    "\n",
    "p_tree\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "14a09e82",
   "metadata": {
    "vscode": {
     "languageId": "r"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "R version 4.2.1 (2022-06-23)\n",
       "Platform: x86_64-apple-darwin13.4.0 (64-bit)\n",
       "Running under: macOS Big Sur ... 10.16\n",
       "\n",
       "Matrix products: default\n",
       "BLAS/LAPACK: /Users/manonmorin/miniconda3/envs/R_analysis_env/lib/libopenblasp-r0.3.21.dylib\n",
       "\n",
       "locale:\n",
       "[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8\n",
       "\n",
       "attached base packages:\n",
       "[1] stats     graphics  grDevices utils     datasets  methods   base     \n",
       "\n",
       "other attached packages:\n",
       " [1] ape_5.7-1       ggtree_3.6.0    forcats_0.5.2   stringr_1.5.0  \n",
       " [5] dplyr_1.0.10    purrr_0.3.5     readr_2.1.3     tidyr_1.2.1    \n",
       " [9] tibble_3.1.8    ggplot2_3.3.6   tidyverse_1.3.2\n",
       "\n",
       "loaded via a namespace (and not attached):\n",
       " [1] Rcpp_1.0.9          lubridate_1.8.0     lattice_0.20-45    \n",
       " [4] assertthat_0.2.1    digest_0.6.30       utf8_1.2.2         \n",
       " [7] IRdisplay_1.1       R6_2.5.1            cellranger_1.1.0   \n",
       "[10] repr_1.1.4          backports_1.4.1     reprex_2.0.2       \n",
       "[13] evaluate_0.17       ggfun_0.1.3         httr_1.4.4         \n",
       "[16] pillar_1.8.1        yulab.utils_0.0.4   rlang_1.0.6        \n",
       "[19] lazyeval_0.2.2      googlesheets4_1.0.1 uuid_1.1-0         \n",
       "[22] readxl_1.4.1        labeling_0.4.2      googledrive_2.0.0  \n",
       "[25] munsell_0.5.0       broom_1.0.1         compiler_4.2.1     \n",
       "[28] modelr_0.1.9        gridGraphics_0.5-1  pkgconfig_2.0.3    \n",
       "[31] base64enc_0.1-3     htmltools_0.5.3     tidyselect_1.2.0   \n",
       "[34] fansi_1.0.3         crayon_1.5.2        tzdb_0.3.0         \n",
       "[37] dbplyr_2.2.1        withr_2.5.0         grid_4.2.1         \n",
       "[40] nlme_3.1-160        jsonlite_1.8.3      gtable_0.3.1       \n",
       "[43] lifecycle_1.0.3     DBI_1.1.3           magrittr_2.0.3     \n",
       "[46] scales_1.2.1        tidytree_0.4.5      cli_3.4.1          \n",
       "[49] stringi_1.7.8       farver_2.1.1        fs_1.5.2           \n",
       "[52] xml2_1.3.3          ellipsis_0.3.2      generics_0.1.3     \n",
       "[55] vctrs_0.5.0         IRkernel_1.3        treeio_1.22.0      \n",
       "[58] tools_4.2.1         ggplotify_0.1.2     glue_1.6.2         \n",
       "[61] hms_1.1.2           parallel_4.2.1      fastmap_1.1.0      \n",
       "[64] colorspace_2.0-3    gargle_1.2.1        aplot_0.2.2        \n",
       "[67] rvest_1.0.3         pbdZMQ_0.3-8        haven_2.5.1        \n",
       "[70] patchwork_1.1.1    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sessionInfo()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}