16S_2023_solutions.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "colab_type": "text",
        "id": "view-in-github"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/Gibbons-Lab/isb_course_2023/blob/main/16S_2023_solutions.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "nQRlv0ypvthh"
      },
      "source": [
        "# 🦠 Amplicon Sequencing Data Analysis with Qiime 2\n",
        "\n",
        "This notebook will accompany the first session of the 2023 ISB Virtual Microbiome Series. The presentation slides can be [found here](https://gibbons-lab.github.io/isb_course_2023/16S). \n",
        "\n",
        "Save your own local copy of this notebook by using `File > Save a copy in Drive`. At some point you may be prompted to trust the notebook. We promise that it is safe 🤞\n",
        "\n",
        "**Disclaimer:**\n",
        "\n",
        "The Google Colab notebook environment will interpret any command as Python code by default. If we want to run bash commands we will have to prefix them by `!`. So any command you see with a leading `!` is a bash command and if you wanted to run it in your terminal you would omit the `!`. For example, if in the Colab notebook you ran `!wget` you would just run `wget` in your terminal. \n",
        "\n",
        "## Setup\n",
        "\n",
        "QIIME 2 is usually installed by following the [official installation instructions](https://docs.qiime2.org/2023.7/install/). However, because we are using Google Colab and there are some caveats to using conda here, we will have to hack around the installation a little bit. But no worries, we provide a setup script below which does all this work for us. 😌 \n",
        "\n",
        "So...let's start by pulling a local copy of the project repository down from GitHub."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 1,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "QFsZdvOuvZmq",
        "outputId": "13f9143b-51a7-4507-ae11-fc2436b27d80"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "Cloning into 'materials'...\n",
            "remote: Enumerating objects: 1124, done.\u001b[K\n",
            "remote: Counting objects: 100% (313/313), done.\u001b[K\n",
            "remote: Compressing objects: 100% (162/162), done.\u001b[K\n",
            "remote: Total 1124 (delta 168), reused 283 (delta 151), pack-reused 811\u001b[K\n",
            "Receiving objects: 100% (1124/1124), 263.50 MiB | 27.26 MiB/s, done.\n",
            "Resolving deltas: 100% (470/470), done.\n",
            "Checking out files: 100% (578/578), done.\n"
          ]
        }
      ],
      "source": [
        "!git clone https://github.com/gibbons-lab/isb_course_2023 materials"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "2eenIEY76mcv"
      },
      "source": [
        "We will switch to working within the `materials` directory for the rest of the notebook."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 2,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "h3NALJ7u6mBP",
        "outputId": "cd899a8e-5a12-4f4f-a347-d636f2678f5a"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "/content/materials\n"
          ]
        }
      ],
      "source": [
        "%cd materials"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "fQVnDjFauZ-n"
      },
      "source": [
        "Now we are ready to set up our environment. This will take about 10-15 minutes. \n",
        "\n",
        "**Note**: This setup is only relevant for Google Colaboratory and will not work on your local machine. Please follow the [official installation instructions](https://docs.qiime2.org/2021.8/install/) for that."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 3,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 272
        },
        "id": "BAyTazMCvrpf",
        "outputId": "1aa2d30f-ab4a-4f2e-c1e3-b75dfe331b38"
      },
      "outputs": [
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:37:15] </span>🐍 Downloading miniconda<span style=\"color: #808000; text-decoration-color: #808000\">...</span>                                                              <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:37:15]\u001b[0m\u001b[2;36m \u001b[0m🐍 Downloading miniconda\u001b[33m...\u001b[0m                                                              \u001b]8;id=581195;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=978348;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:37:16] </span><span style=\"color: #000080; text-decoration-color: #000080\">🐍 Done.</span>                                                                                 <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:37:16]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m🐍 Done.\u001b[0m                                                                                 \u001b]8;id=661823;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=942666;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>🐍 Installing miniconda<span style=\"color: #808000; text-decoration-color: #808000\">...</span>                                                               <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m🐍 Installing miniconda\u001b[33m...\u001b[0m                                                               \u001b]8;id=188324;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=779753;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:37:36] </span><span style=\"color: #000080; text-decoration-color: #000080\">🐍 Installed miniconda to `/usr/local`.</span>                                                  <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:37:36]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m🐍 Installed miniconda to `\u001b[0m\u001b[34m/usr/\u001b[0m\u001b[34mlocal\u001b[0m\u001b[34m`.\u001b[0m                                                  \u001b]8;id=744886;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=513412;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>🔍 Downloading Qiime <span style=\"color: #008080; text-decoration-color: #008080; font-weight: bold\">2</span> package list<span style=\"color: #808000; text-decoration-color: #808000\">...</span>                                                   <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m🔍 Downloading Qiime \u001b[1;36m2\u001b[0m package list\u001b[33m...\u001b[0m                                                   \u001b]8;id=948143;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=275161;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:37:37] </span><span style=\"color: #000080; text-decoration-color: #000080\">🔍 Done.</span>                                                                                 <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:37:37]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m🔍 Done.\u001b[0m                                                                                 \u001b]8;id=131093;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=438406;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>🔍 Installing Qiime <span style=\"color: #008080; text-decoration-color: #008080; font-weight: bold\">2</span>. This may take a little bit.                                       <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "<span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span> 🕐                                                                                      <span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">               </span>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m🔍 Installing Qiime \u001b[1;36m2\u001b[0m. This may take a little bit.                                       \u001b]8;id=225727;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=651063;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n",
              "\u001b[2;36m           \u001b[0m 🕐                                                                                      \u001b[2m               \u001b[0m\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:46:23] </span><span style=\"color: #000080; text-decoration-color: #000080\">🔍 Done.</span>                                                                                 <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:46:23]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m🔍 Done.\u001b[0m                                                                                 \u001b]8;id=659860;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=293696;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>🌲 Installing Empress<span style=\"color: #808000; text-decoration-color: #808000\">...</span>                                                                 <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m🌲 Installing Empress\u001b[33m...\u001b[0m                                                                 \u001b]8;id=739409;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=194963;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:46:27] </span><span style=\"color: #000080; text-decoration-color: #000080\">🌲 Done.</span>                                                                                 <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:46:27]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m🌲 Done.\u001b[0m                                                                                 \u001b]8;id=165586;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=5683;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>📊 Checking that Qiime <span style=\"color: #008080; text-decoration-color: #008080; font-weight: bold\">2</span> command line works<span style=\"color: #808000; text-decoration-color: #808000\">...</span>                                           <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#56\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">56</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m📊 Checking that Qiime \u001b[1;36m2\u001b[0m command line works\u001b[33m...\u001b[0m                                           \u001b]8;id=794084;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=587137;file:///content/materials/setup_qiime2#56\u001b\\\u001b[2m56\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">[15:46:37] </span><span style=\"color: #000080; text-decoration-color: #000080\">📊 Qiime </span><span style=\"color: #000080; text-decoration-color: #000080; font-weight: bold\">2</span><span style=\"color: #000080; text-decoration-color: #000080\"> command line looks good 🎉</span>                                                    <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#62\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">62</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m[15:46:37]\u001b[0m\u001b[2;36m \u001b[0m\u001b[34m📊 Qiime \u001b[0m\u001b[1;34m2\u001b[0m\u001b[34m command line looks good 🎉\u001b[0m                                                    \u001b]8;id=915305;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=567084;file:///content/materials/setup_qiime2#62\u001b\\\u001b[2m62\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span>🧹 Cleaned up unneeded files.                                                            <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#51\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">51</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m🧹 Cleaned up unneeded files.                                                            \u001b]8;id=71343;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=535573;file:///content/materials/setup_qiime2#51\u001b\\\u001b[2m51\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        },
        {
          "data": {
            "text/html": [
              "<pre style=\"white-space:pre;overflow-x:auto;line-height:normal;font-family:Menlo,'DejaVu Sans Mono',consolas,'Courier New',monospace\"><span style=\"color: #7fbfbf; text-decoration-color: #7fbfbf\">           </span><span style=\"color: #008000; text-decoration-color: #008000\">Everything is A-OK. You can start using Qiime </span><span style=\"color: #008000; text-decoration-color: #008000; font-weight: bold\">2</span><span style=\"color: #008000; text-decoration-color: #008000\"> now 👍</span>                                  <a href=\"file:///content/materials/setup_qiime2\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">setup_qiime2</span></a><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">:</span><a href=\"file:///content/materials/setup_qiime2#159\" target=\"_blank\"><span style=\"color: #7f7f7f; text-decoration-color: #7f7f7f\">159</span></a>\n",
              "</pre>\n"
            ],
            "text/plain": [
              "\u001b[2;36m          \u001b[0m\u001b[2;36m \u001b[0m\u001b[32mEverything is A-OK. You can start using Qiime \u001b[0m\u001b[1;32m2\u001b[0m\u001b[32m now 👍\u001b[0m                                  \u001b]8;id=49598;file:///content/materials/setup_qiime2\u001b\\\u001b[2msetup_qiime2\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=146594;file:///content/materials/setup_qiime2#159\u001b\\\u001b[2m159\u001b[0m\u001b]8;;\u001b\\\n"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        }
      ],
      "source": [
        "%run setup_qiime2"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "BGWsqUfq5oSi"
      },
      "source": [
        "---\n",
        "\n",
        "While this is running let's get back to the presentation. 📚\n",
        "\n",
        "---\n",
        "\n",
        "<br>\n",
        "\n",
        "## Our first QIIME 2 command\n",
        "\n",
        "Let's remember our workflow for today.\n",
        "\n",
        "![our workflow](https://github.com/Gibbons-Lab/isb_course_2023/raw/main/docs/16S/assets/steps.png)\n",
        "\n",
        "The first thing we have to do is to get the data into an 'artifact'.\n",
        "We can import the data with the `import` action from the tools. For that we have to give\n",
        "QIIME 2 a *manifest* (list of raw files) and tell it what *type of data* we\n",
        "are importing and what *type of artifact* we want. \n",
        "\n",
        "**QoL Tip:** QIIME 2 commands can get very long. To split them up over several lines we can use `\\` which means \"continue on the next line\"."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 4,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "6mGgUrRf5pyb",
        "outputId": "a61cefdd-0f79-44fe-85e1-76a9bcca45bf"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mImported data/manifest.tsv as PairedEndFastqManifestPhred33V2 to sequences.qza\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime tools import \\\n",
        "  --type 'SampleData[PairedEndSequencesWithQuality]' \\\n",
        "  --input-path data/manifest.tsv \\\n",
        "  --output-path sequences.qza \\\n",
        "  --input-format PairedEndFastqManifestPhred33V2"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "lw7TmzbiAoLj"
      },
      "source": [
        "Since we have quality information for the sequencing reads, let's also generate\n",
        "our first visualization to inspect sequence quality. \n",
        "\n",
        "---\n",
        "\n",
        "Qiime 2 commands can become pretty long. Here are some pointers to remember the\n",
        "structure of a command:\n",
        "\n",
        "```\n",
        "qiime plugin action --i-argument1 ... --o-argument2 ...\n",
        "```\n",
        "\n",
        "Argument types usually begin with a letter denoting their meaning:\n",
        "\n",
        "- `--i-...` = input files\n",
        "- `--o-...` = output files\n",
        "- `--p-...` = parameters\n",
        "- `--m-...` = metadata\n",
        "\n",
        "---\n",
        "\n",
        "In this case we will use the `summarize` action from the `demux` plugin with the previously generated artifact as input and output the resulting visualization to the `qualities.qzv` file."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 5,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "feGa41cEApGs",
        "outputId": "251a9eda-e135-459a-87bc-161f96328bd0"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: qualities.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime demux summarize --i-data sequences.qza --o-visualization qualities.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "O2GenFy3C8gu"
      },
      "source": [
        "You can view the plot by downloading the .qzv file and opening it using http://view.qiime2.org. To download the file click on the folder symbol to the left, open the `materials` folder, and choose download from the dot menu next to the `qualities.qzv` file.\n",
        "\n",
        "🤔 What do you observe across the read? Where would you truncate the reads?"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "os5MeRvLIq0v"
      },
      "source": [
        "# Denoising amplicon sequence variants\n",
        "\n",
        "We will now run the DADA2 plugin which will do 3 things:\n",
        "\n",
        "1. filter and trim the reads\n",
        "2. find the most likely set of unique sequences in the sample (ASVs)\n",
        "3. remove chimeras\n",
        "4. count the abundances of each ASV\n",
        "\n",
        "\n",
        "Since this step takes a bit, let's start the process and use the time to\n",
        "understand what is happening:"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 6,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "RdLSMp9-Dl-d",
        "outputId": "9ea7f19f-2d6a-4d04-c57e-5615ad321d3e"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "Running external command line application(s). This may print messages to stdout and/or stderr.\n",
            "The command(s) being run are below. These commands cannot be manually re-run as they will depend on temporary files that no longer exist.\n",
            "\n",
            "Command: run_dada.R --input_directory /tmp/tmpr1ug5vyu/forward --input_directory_reverse /tmp/tmpr1ug5vyu/reverse --output_path /tmp/tmpr1ug5vyu/output.tsv.biom --output_track /tmp/tmpr1ug5vyu/track.tsv --filtered_directory /tmp/tmpr1ug5vyu/filt_f --filtered_directory_reverse /tmp/tmpr1ug5vyu/filt_r --truncation_length 150 --truncation_length_reverse 150 --trim_left 0 --trim_left_reverse 0 --max_expected_errors 2.0 --max_expected_errors_reverse 2.0 --truncation_quality_score 2 --min_overlap 12 --pooling_method independent --chimera_method consensus --min_parental_fold 1.0 --allow_one_off False --num_threads 2 --learn_min_reads 1000000\n",
            "\n",
            "R version 4.1.3 (2022-03-10) \n",
            "Loading required package: Rcpp\n",
            "\u001b[?25hDADA2: 1.22.0 / Rcpp: 1.0.9 / RcppParallel: 5.1.5 \n",
            "\u001b[?25h\u001b[?25h2) Filtering \u001b[?25h\u001b[?25h.........\u001b[?25h\u001b[?25h\u001b[?25h\n",
            "\u001b[?25h\u001b[?25h3) Learning Error Rates\n",
            "\u001b[?25h26887650 total bases in 179251 reads from 9 samples will be used for learning the error rates.\n",
            "26887650 total bases in 179251 reads from 9 samples will be used for learning the error rates.\n",
            "\u001b[?25h3) Denoise samples .........\n",
            ".........\n",
            "\u001b[?25h5) Remove chimeras (method = consensus)\n",
            "\u001b[?25h\u001b[?25h6) Report read numbers through the pipeline\n",
            "\u001b[?25h\u001b[?25h7) Write output\n",
            "\u001b[?25h\u001b[?25h\u001b[?25h\u001b[?25h\u001b[?25h\u001b[?25h\u001b[?25h\u001b[32mSaved FeatureTable[Frequency] to: dada2/table.qza\u001b[0m\n",
            "\u001b[32mSaved FeatureData[Sequence] to: dada2/representative_sequences.qza\u001b[0m\n",
            "\u001b[32mSaved SampleData[DADA2Stats] to: dada2/denoising_stats.qza\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime dada2 denoise-paired \\\n",
        "    --i-demultiplexed-seqs sequences.qza \\\n",
        "    --p-trunc-len-f 150 \\\n",
        "    --p-trunc-len-r 150 \\\n",
        "    --p-n-threads 2 \\\n",
        "    --output-dir dada2 --verbose"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "aFpQQ5GHdbzX"
      },
      "source": [
        "---\n",
        "\n",
        "While this is running let's get back to the presentation. 📚\n",
        "\n",
        "---\n",
        "\n",
        "<br>\n",
        "\n",
        "\n",
        "If this step takes too long or fails, you can also copy the results from the treasure chest with the following command."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 7,
      "metadata": {
        "id": "PFrfo-uCdoRz"
      },
      "outputs": [],
      "source": [
        "# obscure magic that will only copy if the previous command failed\n",
        "![ -d dada2 ] || cp -r treasure_chest/dada2 ."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "4YzrovxVj7U_"
      },
      "source": [
        "Ok, this step ran, but we should also make sure it kind of worked. One good way to tell if the identified ASVs are representative of the sample is to see how many reads were maintained throughout the pipeline. Here, the most common issues and solutions are:\n",
        "\n",
        "**Large fraction of reads is lost during merging (only paired-end)**\n",
        "\n",
        "![read overlap](https://gibbons-lab.github.io/isb_course_2023/16S/assets/read_overlap.png)\n",
        "\n",
        "In order to merge ASVs DADA2 uses an overlap of 12 bases between forward and reverse reads by default. Thus, your reads must allow for sufficient overlap *after* trimming. So if your amplified region is 450bp long and you have 2x250bp reads and you trim the last 30 bases of each read, truncating the length to 220bp, the total length of covered sequence is 2x220 = 440 which is shorter than 450bp so there will be no overlap. To solve this issue trim less of the reads or adjust the `--p-min-overlap` parameters to something lower (but not too low).\n",
        "\n",
        "<br>\n",
        "\n",
        "**Most of the reads are lost as chimeric**\n",
        "\n",
        "![read overlap](https://gibbons-lab.github.io/isb_course_2023/16S/assets/chimera.png)\n",
        "\n",
        "This is usually an experimental issue as chimeras are introduced during amplification. If you can adjust your PCR, try to run fewer cycles. Chimeras can also be introduced by incorrect merging. If your minimum overlap is too small ASVs may be merged randomly. Possible fixes are to increase the `--p-min-overlap` parameter or run the analysis on the forward reads only (in our empirical observations, chimeras are more likely to be introduced in the joined reads). *However, losing between 5-25% of your reads to chimeras is normal and does not require any adjustments.*\n",
        "\n",
        "Our denoising stats are contained in an artifact. To convert it to a visualization we can use `qiime metadata tabulate`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 8,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "n_Tr3Ujcj61w",
        "outputId": "1bd3b467-67b0-4635-fe1b-3175739b1946"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: dada2/denoising-stats.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime metadata tabulate \\\n",
        "    --m-input-file dada2/denoising_stats.qza \\\n",
        "    --o-visualization dada2/denoising-stats.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "o75g7DSO2tu5"
      },
      "source": [
        "What proprotion of reads was retained throughout the entire pipeline? Look at the final number of used reads (non-chimeric). What do you observe when comparing those values between samples and how might that affect diversity metrics?"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "c0ZV3_sGQ20u"
      },
      "source": [
        "# Phylogeny and ecological diversity metrics"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "aEit8p7V9Ht9"
      },
      "source": [
        "## Building a tree\n",
        "\n",
        "We can build a phylogenetic tree for our sequences using the following command:"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 9,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "2PVNlD0g9MsX",
        "outputId": "f3921866-1c1e-4407-d9f7-9ea78055f1ae"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved FeatureData[AlignedSequence] to: tree/alignment.qza\u001b[0m\n",
            "\u001b[32mSaved FeatureData[AlignedSequence] to: tree/masked_alignment.qza\u001b[0m\n",
            "\u001b[32mSaved Phylogeny[Unrooted] to: tree/tree.qza\u001b[0m\n",
            "\u001b[32mSaved Phylogeny[Rooted] to: tree/rooted_tree.qza\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime phylogeny align-to-tree-mafft-fasttree \\\n",
        "    --i-sequences dada2/representative_sequences.qza \\\n",
        "    --output-dir tree"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "4lB8NI0IiJgk"
      },
      "source": [
        "We can create a visualization for the tree using the [empress](https://github.com/biocore/empress) Qiime 2 plugin."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 10,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "cSS_y3Kligrc",
        "outputId": "e28721ef-90cb-4c45-b407-867e1173034e"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: tree/empress.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime empress tree-plot \\\n",
        "    --i-tree tree/rooted_tree.qza \\\n",
        "    --o-visualization tree/empress.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "UmxTRou99a_n"
      },
      "source": [
        "This looks nice but is not particularly informative as is. The main utility of our tree will be in complementing our diversity analyses. It will tell us which ASVs are more or less related to one another, which will allow us to calculate different kinds of ecological diversity metrics. But what are those?\n",
        "\n",
        "---\n",
        "\n",
        "Let's get back to the presentation. 📚\n",
        "\n",
        "---\n",
        "\n",
        "<br>\n",
        "\n",
        "## Alpha and Beta Diversity\n",
        "\n",
        "![sample sources](https://github.com/Gibbons-Lab/isb_course_2023/raw/main/docs/16S/assets/sample_sources.png)\n",
        "\n",
        "One of our main goals will be to compare the microbial composition across different populations. QIIME 2 has \"run-all\" command for diversity analyses. This will\n",
        "\n",
        "1. Subsample our samples to the same total number of reads (Why?) \n",
        "2. Calculate alpha and beta diversity measures\n",
        "3. Visualize PCoA projections"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "p-uZBNVX-y2L"
      },
      "source": [
        "## Statistical analyses\n",
        "\n",
        "Let's first have a look at alpha diversity. Can we see a difference in the per-sample diversity across populations? "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 11,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "0js6xHuw-NqF",
        "outputId": "5071d344-76ac-4b07-b452-a0ffa4be2d47"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved FeatureTable[Frequency] to: diversity/rarefied_table.qza\u001b[0m\n",
            "\u001b[32mSaved SampleData[AlphaDiversity] to: diversity/faith_pd_vector.qza\u001b[0m\n",
            "\u001b[32mSaved SampleData[AlphaDiversity] to: diversity/observed_features_vector.qza\u001b[0m\n",
            "\u001b[32mSaved SampleData[AlphaDiversity] to: diversity/shannon_vector.qza\u001b[0m\n",
            "\u001b[32mSaved SampleData[AlphaDiversity] to: diversity/evenness_vector.qza\u001b[0m\n",
            "\u001b[32mSaved DistanceMatrix to: diversity/unweighted_unifrac_distance_matrix.qza\u001b[0m\n",
            "\u001b[32mSaved DistanceMatrix to: diversity/weighted_unifrac_distance_matrix.qza\u001b[0m\n",
            "\u001b[32mSaved DistanceMatrix to: diversity/jaccard_distance_matrix.qza\u001b[0m\n",
            "\u001b[32mSaved DistanceMatrix to: diversity/bray_curtis_distance_matrix.qza\u001b[0m\n",
            "\u001b[32mSaved PCoAResults to: diversity/unweighted_unifrac_pcoa_results.qza\u001b[0m\n",
            "\u001b[32mSaved PCoAResults to: diversity/weighted_unifrac_pcoa_results.qza\u001b[0m\n",
            "\u001b[32mSaved PCoAResults to: diversity/jaccard_pcoa_results.qza\u001b[0m\n",
            "\u001b[32mSaved PCoAResults to: diversity/bray_curtis_pcoa_results.qza\u001b[0m\n",
            "\u001b[32mSaved Visualization to: diversity/unweighted_unifrac_emperor.qzv\u001b[0m\n",
            "\u001b[32mSaved Visualization to: diversity/weighted_unifrac_emperor.qzv\u001b[0m\n",
            "\u001b[32mSaved Visualization to: diversity/jaccard_emperor.qzv\u001b[0m\n",
            "\u001b[32mSaved Visualization to: diversity/bray_curtis_emperor.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime diversity core-metrics-phylogenetic \\\n",
        "    --i-table dada2/table.qza \\\n",
        "    --i-phylogeny tree/rooted_tree.qza \\\n",
        "    --p-sampling-depth 8000 \\\n",
        "    --m-metadata-file data/metadata.tsv \\\n",
        "    --output-dir diversity"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 12,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "GnkeX6iY--g-",
        "outputId": "b0fa8b45-3052-4bc5-cb51-69be0bfede7d"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: diversity/alpha_groups.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime diversity alpha-group-significance \\\n",
        "    --i-alpha-diversity diversity/shannon_vector.qza \\\n",
        "    --m-metadata-file data/metadata.tsv \\\n",
        "    --o-visualization diversity/alpha_groups.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ngH5fUeDQMO0"
      },
      "source": [
        "Now, let's use beta diversity to see how different the samples are from one another. First download `diversity/weighted_unifrac_emperor.qzv` and take a look. Do samples separate based on the population?\n",
        "\n",
        "We can check whether that separation is 'significant' by using a PERMANOVA test."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 13,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "KzNW92riQguX",
        "outputId": "e002325d-8ed6-41f0-e6c4-85ed1551a175"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: diversity/permanova.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime diversity adonis \\\n",
        "    --i-distance-matrix diversity/weighted_unifrac_distance_matrix.qza \\\n",
        "    --m-metadata-file data/metadata.tsv \\\n",
        "    --p-formula \"ethnic_group\" \\\n",
        "    --p-n-jobs 2 \\\n",
        "    --o-visualization diversity/permanova.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "5OgGG9yabL4d"
      },
      "source": [
        "---\n",
        "\n",
        "Let's get back to the presentation. 📚\n",
        "\n",
        "---\n",
        "\n",
        "<br>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "FB7WYPZ1DEoR"
      },
      "source": [
        "# Taxonomy\n",
        "\n",
        "We will use a Bayes classifier trained on the NCBI Refseq database version 212. However, there are many more databases which can be downloaded from https://docs.qiime2.org/2023.8/data-resources/. For instance, the SILVA database is very popular and probably the largest 16S database available right now."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 14,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "EEBbgvnhDGTs",
        "outputId": "33ed815a-4262-4b61-e0bb-ec412743db5c"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved FeatureData[Taxonomy] to: taxa.qza\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime feature-classifier classify-sklearn \\\n",
        "    --i-reads dada2/representative_sequences.qza \\\n",
        "    --i-classifier ncbi-refseq-genus-515f-806r.qza \\\n",
        "    --p-n-jobs 2 \\\n",
        "    --o-classification taxa.qza"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "xISxzBo6Mpka"
      },
      "source": [
        "Now let's have a look at the relative abundances of the different bacterial taxa we have in each sample:"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 15,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "4VTL2yEBMqfO",
        "outputId": "aa94b2f2-d005-4078-d67e-ca702a316c18"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: taxa_barplot.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime taxa barplot \\\n",
        "    --i-table dada2/table.qza \\\n",
        "    --i-taxonomy taxa.qza \\\n",
        "    --m-metadata-file data/metadata.tsv \\\n",
        "    --o-visualization taxa_barplot.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "4_FVM8VfX_yQ"
      },
      "source": [
        "We can also collapse data on a particular taxonomic rank using the QIIME 2 [taxa plugin](https://docs.qiime2.org/2021.4/plugins/available/taxa/). Why might we want to look at different taxonomic ranks, rather than just looking at ASVs?"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 16,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "gl72G8KaYmBm",
        "outputId": "d660beca-57fe-4892-e71e-3bc3c5443b27"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved FeatureTable[Frequency] to: genus.qza\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime taxa collapse \\\n",
        "    --i-table treasure_chest/dada2/table.qza \\\n",
        "    --i-taxonomy treasure_chest/taxa.qza \\\n",
        "    --p-level 6 \\\n",
        "    --o-collapsed-table genus.qza"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "dlVdzgfbZcjt"
      },
      "source": [
        "We can export the table and convert it to a .csv file so that we can analyze these data using tools outside of the QIIME 2 environment."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 17,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "ZsxKf3yAZiHk",
        "outputId": "2f9e5c0c-23e9-4ae6-e7a0-b1484125a6e1"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mExported genus.qza as BIOMV210DirFmt to directory exported\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "!qiime tools export \\\n",
        "    --input-path genus.qza \\\n",
        "    --output-path exported\n",
        "!biom convert -i exported/feature-table.biom -o genus.tsv --to-tsv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "TezGLSEGZrcn"
      },
      "source": [
        "Now the data are in a common format and we can use them, for instance, to draw a heatmap using Pandas and Seaborn. No need to understand every bit of code here. This just serves to illustrate that you can get data out of QIIME 2 for custom visualizations."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 18,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 460
        },
        "id": "p5zUjU8JZqL4",
        "outputId": "bc407aad-601c-4aaa-cb16-e220d650a34e"
      },
      "outputs": [
        {
          "data": {
            "text/plain": [
              "<seaborn.matrix.ClusterGrid at 0x7f04bda1d850>"
            ]
          },
          "execution_count": 18,
          "metadata": {},
          "output_type": "execute_result"
        },
        {
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAABQgAAAGqCAYAAABH+6mYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzde7hkVX3n//eHbqC5CWKCIoiNSiSI1xyJKD8DjjHxig5Ro6AhIDCjxiTGMMaItIAjOsSJgUjSGAXByUQEBU2Ui+KNIHAQ5CZeAqgYiMNFkUA30Hx/f9Q+dNF2c7rrrOqq0/V+PU89p2rXrk+tXbXrcr611tqpKiRJkiRJkiRNpo1G3QBJkiRJkiRJo2OBUJIkSZIkSZpgFgglSZIkSZKkCWaBUJIkSZIkSZpgFgglSZIkSZKkCbZwWMErHvjKxB8eecFGe2fUbZAkSZIkSZIeztAKhNx//9Ci541NRt0ASZIkSZIk6eENr0B4331Di543LBBKkiRJkiRpzA2vQLhs2dCi540tRt0ASZIkSZIk6eENrUAYexBKkiRJkiRJY2+IQ4zvHVq0JEmSNixTU1PvBBaNuh2SJG2glk1PTx876kZofA2vQHivPQglSZK01hZNT08vGXUjJEnaEE1NTS0ZdRs03oY3xHj58mFFS5IkSZIkSWrEHoSSJEmSJEnSBBviHIQWCCVJkiRJkqRxN7wC4TIPUiJJkiRJkiSNO3sQSpIkSZIkSRNsiHMQ3j+0aEmSJEmSJEltDK9AuNwhxpIkSZIkSdK4swehJEmSJEmSNMGGOAfhiqFFS5IkSZIkSWrDHoSSJEmSJEnSBBtagbCWWyDMqBsgSZIkSZIkzWKIPQgdYixJkiRJkiSNu+H1ILRAKEmSJEmSJI29IQ4xtkAoSZIkSZIkjbsh9iB8YFjRkiRJkiRJkhpJVQ0l+N53vXE4wfPIJv/zEx6nRNLQTE1NvRNYNOp2SFIji4EbR9wGSWpp2fT09LGjboTmjyF/v1/M8D5n3dc3AMPrQXjfxNcHJWnYFk1PTy8ZdSMkSZL0y6amppaMug2ad+bl93v39Q3D0AqEDywbVrIkSZIkSZKkVoZXILxvWMmSJEmSpPlqgqZJWTwhPascXiptAIZYIHT6vbWR5HeBDwMLgI9WlW+skiRJkjZk83IYpVZvQoqg0gZvaAXC+5dbIJxNkgXA3wK/DdwEXJrk7Kq6drQtkyRJkiRJ0qQYWoFwxf0bDSt6Q7IH8IOquh4gyf8F9gUsEEqSJEmSJGm9sEA4REkOBQ7tW7S0qpb2Xd4B+HHf5ZuA31wfbZMkSZKkVa2n+QHX19x8zo3HenlOh/18+jwOUaP9o8U+4PM8YsMbYmyBkK4YuHTWFSVJkiRpPGww8wM6N96D5vVz6vM4dGOxf/g8j97QCoT33Te06A3JT4DH9V3esVsmjbUJOvLcuJuUI+ONM3/plCRJkjTvDa8H4Qp7EK6FS4FdkuxMrzD4+8DrR9skaa2Mxa9M0qhZoJUkSZK0IRheD8IHLBDOpqruT/JW4BxgAfCxqrpmxM2SJEmSJGlgjUcctR414wgQaTWGViC89/4Fw4reoFTVvwD/Mup2SJq/HPI9Ug7zHh2/3EvSBmKI32WG9TntZ9DsxnbEkd/dpNUb3hBjexBK0voytl/ApGHxy70kbVDm1XcZP4PG0zoUmtemcGwRWBPHIcaSJGlWY9hTdRx7b/rPhCRJo9Os0DyG3zGkoRtagXD5CocYSxo/Y1jkaGEcCyVzZaFl/Ixt744xel3vOiavRV8/kiRJWifDG2Jc9iCUNJbGtsihlcakyKL5Y72/rtdQlLxxfbbhYYxLoXJ1LF5KGgtz/HFp0B9n5/17oMN4pQ3XEIcYZ1jRkiRJo+aPDQMY48KlpCFoXEyCtgWlUfy4tF7vb0gcxqs1GlHhHSw2N+EchJIkSZKkYWhahLOgJI29kfyA6ntDG8Obg9AehJIaazTPWIv5+ibiF6oRz+s2qnkVJ+K5lSRJkqR+DjGWNJ+MxZC+CfqFaiwe7/Vpgp5bSZIkSXrQ8AqEZYFQkrT+2MN0ME42LkmSpHHk99T1a2gFwnvtQShpjHnkug3SWPR4nIe9EJ1sfB4b8VQAgxjV9AFz4Xu3JEmj4ffU9cghxpIeYi3/2RzFkeZa88h1kjYEY1EY35D53i1JkiaBQ4wlrcpfaTYwc+hhNJeePuNcHJY05kbUM3IvVv/deJupqam913Nb1uTG7jTOfP/XRJqgH9k1AIfKaj4YYg/CYSVLktaRvSUlzTcjed+a770xx2TI+a5z+AxYU5F2Xd1Iu0Kq/4hrbfkjux6O+4fG3hDnIBxWsiRJw2fPS0nz0Lwecj6ORVr/EZckTQp7EEqS1lrj4REw3gUxe16OmMNxJEmSpPXDg5RIUgMTNO9M06KZBTHNwuE4kiRJ0nowxCHGNaxoSRpHFjIkSZIkSfPS0AqE9zvEWJIkwKGykiRJksabcxBKkjR89jCVpPWs0VGd53LgqRn+sCNJGnvDKxA6wliSJGlWcyxizKujZtubVuvZWBzV2R92JEnzwfDmIFxhhVCSJGktjKSIMaKihb1p+zTq4davRW+3GRZgJUmaIEMcYmyBUJIkaX2zl968st6Lw+uwf+y6lsVG9xFJkjYAQzxIiQVCSZKkEbCX3og1LtJC2yJc06LkpOwjI5oKwOKrJGm9SdVwCnn7/coRE18hPOPWozPqNkgzNt1kh4l/TUoaD3+5+LBmWW999vXNsgC2P/28pnmT4n1PavecfuvWFc2yAM742Ueb5klq6/hdD22a95Vb2h0t8ymPXNAsC2DnLdq9v/3+m29rlgVw3kcf0Szr+v/ctFkWwJ98d2nTvElx97JPNsu690//oVkWwILN25VKbr2u7f722M+cNLF1HHsQSpIkSZIkSRNsiEcxbvvrryRJkiRJkqT2hlYgXG6BUJIkSZIkSRp7w+tByP3DipYkSZIkSZLUyPAKhLlvWNGSJEmSJEmSGhliD0ILhJIkSZIkSdK4G94chFk2rGhJkiRJkiRJjQxxiPHyYUVPhCRLgEOA/9cteldV/cvoWiRJkiRJkqQN0fAKhGUPwgb+d1UdN+pGSJIkSZIk6ZclORn4fFV9etRt6Zfki8D29Gp/XwfeUlUr1rT+EOcgvGdY0ZIkSZIkSZLW7DVVdWeSAJ8GXg383zWtPLQC4f32ICTJocChfYuWVtXSdYh4a5I3AtPAn1XVHU0bKEmSJEmSpLXW1WneARRwJbACeH6StwOPAQ6f6U2Y5M+B1wCbAp+pqiOTLAa+CFwGPAu4BnhjVd2d5D3Ay4HNgH8FDquqSvIV4GJgH2Ab4OCq+nqSzYGTgd2B7wKPpddTcLqq7uyavBDYpGvvGg2vQPiAcxD+xm/8xmOBz/cteuzU1NSSmQtXX331G1esWLHlqrfbbrvtvrTbbrvdtOmmm56ShOuvv36f+++//8tTU1NnrYdmS5IkSZIkTZzZOnoleQrwbuC5VXVrkm2BD9EbyrsXsCtwNvDpJC8CdgH2AAKcneT5wI+AJ9Mr8l2Y5GPAm4HjgBOq6qjuvk4FXgZ8rrv7hVW1R5KXAEcCL+xud0dV7ZZkd+CKVbbnnO7+v0CvF+EaDa1AuKIsEE5Pf2vJLKvMdj0AST4OfH56enqt1pdWZ9NNdjhy1G2QJEmSJGlcdcXAhxv5+QLg9Kq6tVv/9t4IXj5bVQ8A1yZ5dLfui7rT5d3lLekVDH8E/LiqLuyWnwa8jV6BcJ8khwObA9vS6104UyA8s/t7GbC4O78X8OGuLVcnuXKV7fmdJIuAT3ZtP29NGza8OQhXOMR4LpJsX1U3dxdfBVw9yvZIkiRJkiRptfp7yaXv7/ur6u/7V+yGGK863Le6Qt5HgKmq+nGSJcCi1dzHCtahnldVy5KcBezLKAqEKxxiPFcfTPIMejvNjcBho22OJEmSJEnSRPsy8JkkH6qq27ohxmtyDnB0kk9W1V1JdgDu667bKcmeVXUR8HrgG6wsBt6aZEvg95hlWDBwIb05Di9IshvwVIDu9ltV1c1JFgIvpXck4zUaYoHw3mFFT4SqesOo2yBJkiRJkqSeqromyfuAryZZwcrhw6tb99wkvw5c1A1Dvgs4gF4PwO8Cb+nmH7wWOLE7SMlJ9EaQ3gJcuhZN+ghwSpJrgevoDUn+ObAFvTkPNwU2Ai4A/u7hgoZYIHSIsSRJkiRJkjYcVXUKcMrDXL9l3/kP080ROKMbYnx/VR2wmtu+m95BUFZdvnff+VtZOQfhMuCAbhjxE4HzgR9W1b3As9d2m2CIBcIHHGIsSZIkSZIkDcvm9IYXb0xvzsM3d8XBdTa0AmHVfbOvJEmSJEmSJE2IqroR2L1R1i+AqRZZQ+xBeE9mX2t8TE1NLZmenl4y6nZIkiRJkiRJ69NGo26AJEmSJEmSpNGxQChJkiRJkiRNMAuEkiRJkiRJ0gQb2hyEksbL1/Z8XbOsU294RLOsZ2/7QLMsgP923ceb5klq7/0/+oeGWc2iANhoI78aDeKI69s9p621fE7fs/jgZlkAy1a0m7L73+9u+3l627L7m2V94T8/0SxLg/nSs1/TNO8ZT7+lWVbd/4NmWQC/+b1HNcu6d8WCZlkAv7LF3c2ytnjrvzTLAjjp1w9slvX273+sWRbAax55SLOsDzz3J82yAP7nZY9rmtdSbmn3Ol102POaZQFc8KYbm2X91qtva5Y16exBKEmSJOlBLYuDkiRpfrBAKEmSJEmSJE0wC4SSJEmSJEnSBLNAKEmSJEmSJE0wC4SSJEmSJEnSBLNAKEmSJEmSJE0wC4SSJEmSJEnSBLNAKEmSJEmSJE0wC4SSJEmSJEnSBLNAKEmSJEmSJI1AksVJrm6UtWuSi5IsT/KOdbntwhYNkCRJkiRJkjRStwNvA165rjdsWiCcmpp6J7CoZeZ6tHjUDZAkSZIkSdJ4SrIY+CLwTeC5wKXAx4H3AtsB+wPXAMcDuwMbA0uq6qwkBwKvArYGdgBOq6r3dtELkpzUZf4E2Leq7klyCHAosAnwA+ANVXV3kpcD7+6W3wbsX1X/UVU/BX6a5KXrum2texAump6eXtI4c72Ymppasobl87noKT3ob7fYe9RNkCRJkiRpbCU5lF5BbsbSqlq6ympPAl4NHESvQPh6YC/gFcC7gGuBL1fVQUm2AS5Jcn532z3oFQ7vBi5N8s/ArcAuwOuq6pAknwL2A04Dzqyqk7q2HQMcTK/4+A3gOVVVSd4EHA782Vy23SHGs5u3RU+p38W/9Y4jR90GSZIkSZLGVVcMXLUguKobquoqgCTXAF/qCnVX0RuduiPwir45ABcBO3Xnz6uq27rbnkmvsPjZLvOKbp3LWDnKdfeuMLgNsCVwTrd8R+CfkmxPrxfhDYNt8UoepESSJEmSJElaO8v7zj/Qd/kBeh3xAuxXVc/oTjtV1Xe6dWqVrJnL/ZkrWNmh72TgrVX1VHrDmGdGuB4PnNAtP4wGI18tEEqSJEmSJEltnAP8UZIAJHlm33W/nWTbJJvRO5DIhbNkbQXcnGRjevMbztia3lyFAH/QotEWCCVJkiRJkqQ2jqZ3cJIruyHIR/dddwlwBnAlcEZVTc+SdQRwMb1C4nV9y5cApye5jN4chgAkeUySm4C3A+9OclOSR6xNo52DUJIkSZIkSZpFVd1I7yAjM5cPXMN1h60h4qaqeuUsmcf1nT8ROHE17TgLOGs1y2+hNz/hOrMHoSRJkiRJkjTB7EEoSZIkSZIkDVFVnUzvoCNjyR6EkiRJkiRJ0gSzQChJkiRJkiRNMAuEkiRJkiRJ0gSzQChJkiRJkiRNMA9SIk2Iny3fpFnW726/vFnW/tf80pHZ52SThVs1zZsUr9361c2y/unnpzfL0obpE7u9sllW6186D7j2s40TtSE59qZPjboJ85KfzaP34su/0Dbw8rZxk+Ltj31ts6zTn35AsyyAn93bLuufn/W6dmHAq678dLOsI+7eu1kWwPHnPqZZ1n2nfaNZFkBt+tRmWT89YrpZFsCv/Wq7b3A3fXlBsyyAx/9p07h5xR6EkiRJkiRJ0gSzQChJkiRJkiRNMAuEkiRJkiRJ0gSzQChJkiRJkiRNMAuEkiRJkiRJ0gSzQChJkiRJkiRNMAuEkiRJkiRJ0gSzQChJkiRJkiRNMAuEkiRJkiRJ0gSzQChJkiRJkiQNIMnJSX5v1O3ol2TzJP+c5Lok1yQ5drbbLFwfDZsnlk1NTS1ZzfLF67kdkiRJkiRJ0lwcV1UXJNkE+FKSF1fVF9a0sgXCzvT09GqrqWsoGkqSJEmSJGnCJHkj8A6ggCuBFcDzk7wdeAxweFV9ulv3z4HXAJsCn6mqI5MsBr4IXAY8C7gGeGNV3Z3kPcDLgc2AfwUOq6pK8hXgYmAfYBvg4Kr6epLNgZOB3YHvAo8F3lJV08AFAFV1b5JvATs+3HZZIFxPpqam3gksGnU7NLnet/BFo26CJEmSJEljK8mhwKF9i5ZW1dK+658CvBt4blXdmmRb4EPA9sBewK7A2cCnk7wI2AXYAwhwdpLnAz8CnkyvyHdhko8BbwaOA06oqqO6+zoVeBnwue7uF1bVHkleAhwJvLC73R1VtVuS3YErVrNN29ArOn744bbdAuH6s2h6enrJqBuhyXXOc9515KjbIEmSJEnSuOqKgUsfZpUXAKdX1a3d+rcnAfhsVT0AXJvk0d26L+pOl3eXt6RXMPwR8OOqurBbfhrwNnoFwn2SHA5sDmxLr3fhTIHwzO7vZaycDm8vusJfVV2d5Mr+xiZZCPwj8DdVdf3DbbsFQkmSJEmSJGlwy/vOp+/v+6vq7/tX7IYY1yq3rySLgI8AU1X14yRLeOhI1Jn7WMHa1/OWAt+vqr+ebUWPYixJkiRJkiTN7svAq5M8CqAbYrwm5wAHJdmyW3eHJNt11+2UZM/u/OuBb7CyGHhrd5u1OTLyhfTmOCTJbsBTZ65IcgywNfAna7Nh9iCUJEmSJEmSZlFV1yR5H/DVJCtYOXx4deuem+TXgYu6Ych3AQfQ6wH4XeAt3fyD1wIndgcpOQm4GrgFuHQtmvQR4JQk1wLX0RuS/PMkOwJ/2S37Vnf/J1TVR9cUZIFQkiRJkiRJWgtVdQpwysNcv2Xf+Q+zysFBuiHG91fVAau57bvpHQRl1eV7952/lZVzEC4DDqiqZUmeCJwP/LCq7mXlUOe1YoFQkiRJkiRJmn82By5IsjG9guCbu+LgOrNAKEmSJEmSJK0HVXUjsHujrF8AUy2yPEiJJEmSJEmSNMEsEEqSJEmSJEkTzAKhJEmSJEmSNMGcg1CaEPtefuqom6Ax9olb13i0+w3K23Y4pGne3/zkpKZ5LX3uWb90ULSBvfxbpzXLAnjtt9d40DepuZav+3F+zUsP5707t/38236zFc2yXvLkHzXLAtjxM+c3zWvp/T9c2izrY1se1CwLYPOFDzTL+o9lmzTLAjjnWb/bLOuGn7dt2y5Zp4PEPqxaXs2yWvvuLY9qmrfNpsubZd11X9vn9PFN0+YXexBKkiRJkiRJE8wCoSRJkiRJkjTBLBBKkiRJkiRJE8wCoSRJkiRJkjTBLBBKkiRJkiRJE8wCoSRJkiRJkjTBLBBKkiRJkiRJE8wCoSRJkiRJkjTBLBBKkiRJkiRJE2zhqBswDyybmppa0iBncYMMSZIkSZIkjYkkJwOfr6pPj7ot/ZK8D3gj8Miq2nK29S0QzmJ6evrYFjmNioySJEmSJEnSbD4HnAB8f21WdoixJEmSJEmStBaSvDHJlUm+neTUbvHzk/xrkuuT/F7fun+e5NJu/fd2yxYnuS7JJ5N8J8mnk2zeXfeebv2rkyxNkm75V5J8IMklSb6X5P/rlm+e5FNJrk3ymSQXJ5kCqKpvVtXNa7td9iCc56ampt4JLBp1OyRJkiRJkuazJIcCh/YtWlpVS/uufwrwbuC5VXVrkm2BDwHbA3sBuwJnA59O8iJgF2APIMDZSZ4P/Ah4MnBwVV2Y5GPAm4HjgBOq6qjuvk4FXkavJyDAwqraI8lLgCOBF3a3u6OqdkuyO3DFoNtugXD+WzQ9Pb1k1I3Q+Fu06eOOHHUbJEmSJEkaV10xcOnDrPIC4PSqurVb//auk99nq+oB4Nokj+7WfVF3ury7vCW9guGPgB9X1YXd8tOAt9ErEO6T5HBgc2Bb4BpWFgjP7P5exsrjXOwFfLhry9VJrhxgswELhJIkSZIkSdJcLO87n76/76+qv+9fMclioFa5fSVZBHwEmKqqHydZwkNHjM7cxwqGUM9zDkJJkiRJkiRpdl8GXp3kUQDdEOM1OQc4KMmW3bo7JNmuu26nJHt2518PfIOVxcBbu9v8HrO7EHhNl78b8NR12Zh+9iCUJEmSJEmSZlFV1yR5H/DVJCtYOXx4deuem+TXgYu6Ych3AQfQ6wH4XeAt3fyD1wInVtXdSU4CrgZuAS5diyZ9BDglybXAdfSGJP8cIMkH6RUfN09yE/DRqlqypiALhJIkSZIkSdJaqKpTgFMe5vot+85/mG6OwBndEOP7q+qA1dz23fQOgrLq8r37zt/KyjkIlwEHVNWyJE8Ezgd+2K13OHD42m2VBUJJkiRJkiRpPtocuCDJxvTmPHxzVd07SJAFQkmSJEmSJGk9qKobgd0bZf0CmGqR5UFKJEmSJEmSpAlmgVCSJEmSJEmaYBYIJUmSJEmSpAlmgVCSJEmSJEmaYB6kRJoQy++9ZdRNUEN7bv3mpnkX/fwjTfPG1f+64ehRN2G9edE3jxt1E+alD+x6RLOs/3Hd5Oxv42ySXvfSmnzpP+5umnfEUx5olvXX009olgVw5d77Nsu67P9t2ywLYKO0y7rktgXtwoCtN2mX9+St7m+WBfC9X2zRLOs7d7Z93H7zfd9olrX1Cx7ZLAvgga22apb1/EPuaZYFcPNZ7fJ+/O/t9o9JZw9CSZIkSZIkaYJZIJQkSZIkSZImmAVCSZIkSZIkaYJZIJQkSZIkSZImmAVCSZIkSZIkaYJZIJQkSZIkSZIm2MJRN2CCLJuamloyhNzFQ8iUJEmSJEnShLBAuJ5MT08fO4zcIRUdJUmSJEmSNCEcYixJkiRJkiRNMAuEkiRJkiRJ0ggkWZzk6kZZ+ye5MslVSf41ydPX9rYOMZYkSZIkSZLmvxuA36qqO5K8GFgK/Oba3NAehJIkSZIkSdIsut5+1yU5Ocn3knwyyQuTXJjk+0n2SLJFko8luSTJ5Un27W57YJKzknylW/fIvugFSU5Kck2Sc5Ns1t3mkCSXJvl2kjOSbN4tf3mSi7v885M8GqCq/rWq7ugyvwnsuLbbZg9CaYSmpqbeCSwadTskSZIkSZp0SQ4FDu1btLSqlq6y2pOAVwMHAZcCrwf2Al4BvAu4FvhyVR2UZBvgkiTnd7fdA9gduBu4NMk/A7cCuwCvq6pDknwK2A84DTizqk7q2nYMcDBwPPAN4DlVVUneBBwO/Nkq7TwY+MLabrsFQmm0Fk1PTy9ZH3eUbHzk7GtJkiRJkjSZumLgqgXBVd1QVVcBJLkG+FJXqLsKWEyv194rkryjW38RsFN3/ryquq277Zn0Couf7TKv6Na5rMsB2L0rDG4DbAmc0y3fEfinJNsDm9AbWvygJPvQKxDutbbb7hBjSZIkSZIkae0s7zv/QN/lB+h1xAuwX1U9ozvtVFXf6dapVbJmLvdnrmBlh76TgbdW1VOB97JyBOLxwAnd8sP6lpPkacBHgX1nipFrwwKhJEmSJEmS1MY5wB8lCUCSZ/Zd99tJtu3mGHwlcOEsWVsBNyfZGNi/b/nWwE+6838wszDJTsCZwBuq6nvr0mgLhJIkSZIkSVIbRwMbA1d2Q5CP7rvuEuAM4ErgjKqaniXrCOBieoXE6/qWLwFOT3IZvTkMZ7wHeBTwkSRXJJkt/0HOQShJkiRJkiTNoqpupHeQkZnLB67husPWEHFTVb1ylszj+s6fCJy4mnacBZy1muVvAt4023asjj0IJUmSJEmSpAlmD0JJkiRJkiRpiKrqZHoHHRlL9iCUJEmSJEmSJpgFQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkmSJEmSJGmCeZCS+W/Z1NTUklE3QgNbPOoGaH563MZbN827qGmaNH9tlFG3QJLaO+9nf9U078j8abOsTTaqZlkAd927SbOsu1e07U/zyE1WNMu6fXnbx223rdvmtbRpw33kZ/c2iwJgwaYPtAuboC8hDzzQblt/cd+CZlmTzgLhPDc9PX3sqNugwVnclSRJkiRJo+YQY0mSJEmSJGmCWSCUJEmSJEmSJpgFQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkmSJEmSJGkEkixOcnWjrH2TXJnkiiTTSfZa29subNEASZIkSZIkSSP1JeDsqqokTwM+Bey6Nje0B6EkSZIkSZI0i66333VJTk7yvSSfTPLCJBcm+X6SPZJskeRjSS5JcnmSfbvbHpjkrCRf6dY9si96QZKTklyT5Nwkm3W3OSTJpUm+neSMJJt3y1+e5OIu//wkjwaoqruqqrrMLYBiLdmDUBqtZVNTU0tG3QhJkiRJkiZdkkOBQ/sWLa2qpaus9iTg1cBBwKXA64G9gFcA7wKuBb5cVQcl2Qa4JMn53W33AHYH7gYuTfLPwK3ALsDrquqQJJ8C9gNOA86sqpO6th0DHAwcD3wDeE7XU/BNwOHAn3XrvQp4P7Ad8NK13XYLhNIITU9PH7u+7ivZ+MjZ15IkSZIkaTJ1xcBVC4KruqGqrgJIcg3wpa5QdxWwGNgReEWSd3TrLwJ26s6fV1W3dbc9k15h8bNd5hXdOpd1OQC7d4XBbYAtgXO65TsC/5Rke2AT4Ia+bfgM8JkkzweOBl64NtvuEGNJkiRJkiRp7SzvO/9A3+UH6HXEC7BfVT2jO+1UVd/p1ll1yO/M5f7MFazs0Hcy8NaqeirwXnrFRuj1IjyhW35Y3/KVwVVfA56Q5FfWZqMsEEqSJEmSJEltnAP8UZIAJHlm33W/nWTbbo7BVwIXzv4okiIAACAASURBVJK1FXBzko2B/fuWbw38pDv/BzMLkzyp736fBWwK3LY2jbZAKEmSJEmSJLVxNLAxcGU3BPnovusuAc4ArgTOqKrpWbKOAC6mV0i8rm/5EuD0JJfRm8Nwxn7A1UmuAP4WeG3fQUselnMQSpIkSZIkSbOoqhvpHWRk5vKBa7jusDVE3FRVr5wl87i+8ycCJ66mHWcBZ61m+QeAD8y2HatjD0JJkiRJkiRpgtmDUJIkSZIkSRqiqjqZ3kFHxpI9CCVJkiRJkqQJZoFQkiRJkiRJmmAWCCVJkiRJkqQJ5hyE0oT42vP+tFnWzo/6WbOsn965ZbMsgOdf/LlmWYds94ZmWQC/+5hlzbL2u/Ifm2UBbLHZ4qZ5LV2050ubZf30ns2aZQHcfu/GzbKmb2/7kfzeV32/WdatP2j7uO163gVN81pacuOpzbLG+XUlzVdLFrf9bJ6+dUWzrDvvu7dZVmtfW/7ppnm/+eL/1yzr389c1CwLYNn9C5pl/foj/rNZFsCKSrOs03/+f5tlAez5qwc0y9p16zubZQG8/+qtm2V9/q52n/MAB35v32ZZj/2PXzTLAthxi/Pahe2yQ7ssYIcPbtss679+9GvNsiadPQglSZIkSZKkCWaBUJIkSZIkSZpgFgglSZIkSZKkCWaBUJIkSZIkSZpgFgglSZIkSZKkCWaBUJIkSZIkSZpgFgglSZIkSZKkCWaBUJIkSZIkSZpgFgglSZIkSZKkCWaBUJIkSZIkSZpFksVJrp7D7e8a8HZfTPLtJNck+bskCwZtw5pYIJQkSZIkSZLG12uq6unA7sCvAq9ufQcWCCVJkiRJkqS1syDJSV1vvnOTbJbkkCSXdr38zkiyOUCSnZNclOSqJMfMBCQ5KskV3eknST7eLf9sksu67ENn1q+qO7uzC4FNgGq9UQtbB65PU1NT7wQWjbod0nzwvxe9YNRNkCRJkiRpbHVFuUP7Fi2tqqWrrLYL8LqqOiTJp4D9gDOr6qQu4xjgYOB44MPAiVX1iSRvmQmoqvcA70myDfB14ITuqoOq6vYkmwGXJjmjqm7rcs8B9gC+AHy67ZbP8wIhsGh6enrJqBshzQdf3+vwI0fdBkmSJEmSxlVXDFy1ILiqG6rqiu78ZcBiYPeuMLgNsCVwTnf98+gVEAFOBT4wE5IkwGnAh6rqsm7x25K8qjv/OHrFyNu6tv1OkkXAJ4EXAOcNso1r4hBjSZIkSZIkae0s7zu/gl7nu5OBt1bVU4H38tDRrmsaDrwEuKmqZoYX7w28ENizm2/w8lVyqKplwFnAvnPdiFVZIJQkSZIkSZIGtxVwc5KNgf37ll8I/H53/sHlSV5Orxj4tr51twbuqKq7k+wKPKdbd8sk23fnFwIvBa5rvQEWCCVJkiRJkqTBHQFcTK8g2F+8+2PgLUmuAnboW/727vIl3YFKjgK+CCxM8h3gWOCb3bpbAGcnuRK4Avgp8HetN2C+z0EoSZIkSZIkDV1V3Qjs3nf5uL6rT1zN+jcAe/Ytene3fJ813MWL17D82evU0AHYg1CSJEmSJEmaYBYIJUmSJEmSpAlmgVCSJEmSJEmaYBYIJUmSJEmSpAlmgVCSJEmSJEmaYBYIJUmSJEmSpAm2cNQNkLR+POXxP22WtdGCapb14q/f2ywLYOdFz2uWdf6d1zfL6uW1y/rpQc9pFwYsfESaZX3nm9s2ywJI2u1vz37yT5plAdx00zbNsl4ydUezLIBnf3xBw7TxfZ2Os7fuuHPTvCdssaxZ1nab3dMsC+CAK37WNE9ak1NuafvZ3FS7j9LmWn9vWPDUxzbL2uOidt9RAV52Ybv3ysvf3iwKgMvO+5VmWUfs9IZmWQDH39TutfVfdti8WRbANff/sFlW6+8g//LvWzfL2u2uto/bq//9P5plZcWKZlkANx/3g2ZZ2x/+tGZZk84ehJIkSZIkSdIEs0AoSZIkSZIkTTALhJIkSZIkSdIEs0AoSZIkSZIkTTALhJIkSZIkSdIEs0AoSZIkSZIkTbCFjfOWTU1NLWmc+XAWr8f7kiRJkiRJkjY4TQuE09PTx7bMm816LkZKkiRJkiRJGxyHGEuSJEmSJEkTzAKhJEmSJEmSNIski5NcPYfb3zXg7d6X5MeD3n5tWCCUJEmSJEmSxtfngD2GeQcWCCVJkiRJkqS1syDJSUmuSXJuks2SHJLk0iTfTnJGks0Bkuyc5KIkVyU5ZiYgyVFJruhOP0ny8W75Z5Nc1mUfOrN+VX2zqm4e5ka1PorxWJmamnonsGjU7ZDGwblP3n3UTZAkSZIkaWx1RblD+xYtraqlq6y2C/C6qjokyaeA/YAzq+qkLuMY4GDgeODDwIlV9Ykkb5kJqKr3AO9Jsg3wdeCE7qqDqur2JJsBlyY5o6puG8Km/pINukAILJqenl4y6kZI4+D2/Q88ctRtkCRJkiRpXHXFwFULgqu6oaqu6M5fBiwGdu8Kg9sAWwLndNc/j14BEeBU4AMzIUkCnAZ8qKou6xa/LcmruvOPo1eMtEAoSZIkSZIkjZHlfedXAJsBJwOvrKpvJzkQ2LtvnVpDzhLgpqqaGV68N/BCYM+qujvJV1iPo2Kdg1CSJEmSJEka3FbAzUk2BvbvW34h8Pvd+QeXJ3k5vWLg2/rW3Rq4oysO7go8Z7hNfigLhJIkSZIkSdLgjgAuplcQvK5v+R8Db0lyFbBD3/K3d5cv6Q5UchTwRWBhku8AxwLfnFk5yQeT3ARsnuSmJEtab4BDjCVJkiRJkqRZVNWNwO59l4/ru/rE1ax/A7Bn36J3d8v3WcNdvHgN93s4cPg6Nned2INQkiRJkiRJmmAWCCVJkiRJkqQJZoFQkiRJkiRJmmAWCCVJkiRJkqQJZoFQkiRJkiRJmmAWCCVJkiRJkqQJtnDUDZijZVNTU0se5vrF66kd0th71P/5ZLOsw3c+olnWm3eoZlkAB+79g2ZZf/ulXZplAbxkh1ubZb3z7F9plgWw81btslr/8vTWQ29qlrXk+LbP6T0r2u2/293U9jm9+o6jm+Zp3f23O0bdAmk87Lvt/2iWddbtH2iWNUnuvfP1TfM2uebmZlmP/Y0VzbIA9vn2E5tl3fmd65tlASxacH+zrO03a/u4XX1Hu/8VnvTFZlEAHLDdXzbLOu2n72uWBfAnj2vXtm02brd/APC7e7TL+sVd7bKATRfd1izr2+/692ZZAM+4oGncvDKvC4TT09PHPtz1sxQPJUmSJEmSpInnEGNJkiRJkiRpglkglCRJkiRJkiaYBUJJkiRJkiRpglkglCRJkiRJkiaYBUJJkiRJkiRpglkglCRJkiRJkiaYBUJJkiRJkiRpglkglCRJkiRJkmaRZHGSq+dw+7sGuM3mSf45yXVJrkly7KD3/3AsEEqSJEmSJEnj67iq2hV4JvC8JC9ufQcWCCVJkiRJkqS1syDJSV1vvnOTbJbkkCSXJvl2kjOSbA6QZOckFyW5KskxMwFJjkpyRXf6SZKPd8s/m+SyLvtQgKq6u6ou6M7fC3wL2LH1Ri1sHThmlk1NTS0ZdSMkSZIkSZI03rqi3KF9i5ZW1dJVVtsFeF1VHZLkU8B+wJlVdVKXcQxwMHA88GHgxKr6RJK3zARU1XuA9yTZBvg6cEJ31UFVdXuSzYBLk5xRVbf1tW8b4OVdblMbdIFwenp6KOOypfko2fjIUbdBkiRJkqRx1RUDVy0IruqGqrqiO38ZsBjYvSsMbgNsCZzTXf88egVEgFOBD8yEJAlwGvChqrqsW/y2JK/qzj+OXjHytm79hcA/An9TVdcPtIEPY4MuEEqSJEmSJEkNLe87vwLYDDgZeGVVfTvJgcDefevUGnKWADdV1czw4r2BFwJ7VtXdSb4CLOpbfynw/ar66zlvwWo4B6EkSZIkSZI0uK2Am5NsDOzft/xC4Pe78w8uT/JyesXAt/WtuzVwR1cc3BV4Tt/6x3TX/8lwmm+BUJIkSZIkSZqLI4CL6RUEr+tb/sfAW5JcBezQt/zt3eVLugOVHAV8EViY5DvAscA3AZLsCPwlsBvwrW79N7XeAIcYS5IkSZIkSbOoqhuB3fsuH9d39YmrWf8GYM++Re/ulu+zhrt48RqWZ50aOgB7EEqSJEmSJEkTzAKhJEmSJEmSNMEsEEqSJEmSJEkTzAKhJEmSJEmSNMEsEEqSJEmSJEkTzAKhJEmSJEmSNMEsEEqSJEmSJEmTrKo8efLk6cETcOg4Ztm28cizbRtWlm0bjzzbNvos2zYeebZt9Fm2bTzyJqVtk7Kdtm088lq3bUM82YNQ0qoOHdOs1nm2bfRZrfMmpW2Tsp2t82zb6LNa59m20We1zrNto89qnWfbRp/VOm9cs1rn2bbRZ7XOa922DY4FQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkmrWjqmWa3zbNvos1rnTUrbJmU7W+fZttFntc6zbaPPap1n20af1TrPto0+q3XeuGa1zrNto89qnde6bRucdJM1SpIkSZIkSZpA9iCUJEmSJEmSJpgFQkmSJEmSJGmCWSCUJEmSJEmSJpgFQkkasSQbJXnNqNuh8ZRkQZLjRt0OtdO95p876naoHV+nG570PG7U7ZAkaX3xICWSmkqyMfDfged3i74K/F1V3Te6Vq2U5KXAU4BFM8uq6qgGuY8EHldVVw54++mqmpprO1pLsgC4pqp2bZj5x8DHgV8AHwWeCbyzqs4dZVaX9zRgMbBwZllVnTlIVpfXZH9L8s2qes6g7VhN3tbAEuD/6xZ9FTiqqn4+YN4mwK91F787l9f7OO8fXeYjgV146HP6tQFyLq+qZw7ajtXkPRG4qaqWJ9kbeBrwiar62QBZOwN/xC+/Fl7RoJ1zfa/cArinqh5I8mvArsAXxuEzZgiv0+cBV1TVfyY5AHgW8OGq+uE65jR/Hx9nSRYBB/PL770HDZB1VVU9tVG7FtB7Te7fIq/LbPaZ1epxG9L3hqG8J831/ai1JB8EjgHuAb5I7338T6vqtJE2rDPzePHQ5+Bbc8jbjofuaz+aQ9buwG6r5H1igJzm+1qSV9D3f1FVfW7QrC6v2eM2jLxxM+6vq3FjD0JpgiW5c5bTL5J8bx1jTwR+A/hId3pWt2zQNu6S5NNJrk1y/cxpwKy/A15L74M/wKuBx8+hbV9J8ogk2wLfAk5K8qEB485P8o4kj0uy7cxpDm17TpJLk9yV5N4kK5Lcua45VbUC+G6SnQZty2ocVFV3Ai8CHgm8ATh21FlJPgZ8DNgPeHl3etmA7Wq9v12e5Owkb0jyX2dOg7aN3nbeCbymO91Jr5C2zrpi1PeBv6X3mv9ekuc/7I0e3ljuHwBJ3gR8DTgHeG/3d8mAcV9Ksl+SDNqeVZwBrEjyJGApvX/i/s+AWZ8FbgSOB/6q7zSQxu+VXwMWJdkBOJfec3ryHNrW5L2y0/p1eiJwd5KnA38G/Buwzv/wDuN9vOXj1vg5ADgVeAzwO/R+/NiR3o8Eg/hWkmfPoS0P6p6Hx3c/qMxZ688sGj1uQ/re0Ow9qfH7Uev990XdZ9bL6G3vk4A/H4N2keRo4Ergb1j5+A/UazrJK5J8H7iB3r52I/CFObTtSHr7xvHAPsAHgUELeq0//94P/DFwbXd6W5L/OWBW68etWV6SlyW5PMntff9DzmV/G8vX1SRYOPsqkjZg/zZbD5Ykl69j5rOr6ul9l7+c5Nvr3rQHfRw4Evjf9D70/5DBf9x4blU9LcmVVfXeJH/FHD5Yga2r6s6uaPCJqjoyyaC/Qr+2+/uWvmUFPGHAvBOA3wdOB6aAN7Kyh9e6eiRwTZJLgP98sHGD/5o6UxB5CXBqVV0zhyJJy6znVNVuA952dVrub4uA24AX9C0rYNDejU+sqv36Lr83yRUDZv0VvS9f3wVIr1fXP9L7oWAQ47p/QO9L/rOBb1bVPkl2BQb6og8cBrydXlHvnq6tVVWPGDDvgaq6P8mrgOOr6vgB3r9nLKuqvxnwtqvT8r0yVXV3koOBj1TVB+ew70Lb98rWr9P7q6qS7AucUFX/0G33IFq/j7d83FpmATypql6dZN+qOiXJ/wG+PmDWbwL7J/khvcdt5nX6tAHzrgcuTHI2D30eBilQtf7Mavm4td7fWr4ntXw/grb778z/5i8FTq+qn8/hI6v16+o19L473DuHjBlHA88Bzq+qZybZBzhgDnm/BzwduLyq/jDJo4FBe4e1/vx7KfCMqnoAIMkpwOXAuwbIav24tcz7a+C/AldVmyGq4/q62uBZIJQm236zr7JW6/RbkeSJVfVvAEmeAKxY55attFlVfSlJumFVS5JcBrxngKx7ur93J3ksvX/ktp9D2xYm2Z7el6a/nEMOVbXzXG6/hswfJFnQ/Zr/8a5Y8BcDRB3RuGmXJTkX2Bn4iyRbAQ+MQdZFSXarqmsHvP2qmu1vVfWHjdo0454ke1XVN+DBoYz3zHKbNdl4pjgIUFXfS2+qgUGN6/4BvX8cliUhyaZVdV2SJw8SVFVbzaEdq3NfktcBf0CvJxHAoM/Dh7seGecCy2cWzmEoWbP3SnpTw+0J7E9vOCTAgrkEtnqvHMLr9BdJ/oLeP2zPT7IRgz+nrd/HW37GNM0CZoab/yy9YYe3ANsNmPU7A95uTf6tO20EzPU9oPVnVsvHrfX+1vI9qeX70Uw7Wu2/n09yHb3P4/+e5FeBZWPQLoCrgW2Anw7anj73VdVt6c3Hu1FVXZDkr+eQNzPtxP1JHtG1cdD5Q1t//kHvcbu9O7/1HHJaP24t834MXN2oOAiM7+tqQ2eBUJpgVfWQobrdh2r/fBu3r7rOWvhz4IL0hgGH3pDKufzTtLz7p+j7Sd4K/ATYcsCszyfZBvhf9IaVFL15yQZ1FL0hht+oqku7Yuj3BwlKsjm93kQ7VdWhSXYBnlxVnx+wbXenN4zpivTm3riZAXteVtVXkzwe2KWqzu/aOpd/yA8GngFc3/UEehSD7yMtsz5B7x+uW+h9KZxrT5Fm+1uSHekNd3let+jrwB9X1U0Dtu2/A6ekNxdh6H1xPXDArOkkH2Xlr/X7A9MDZsH47h8AN3XP6WeB85LcAazTfHAzup6M+wM7V9XR6R0MYfuqumTAtv0h8N+A91XVDenNo3TqgFlPpTd09wWsLKgWD+0Zty6avVfS68X5F8Bnuh6hTwAuGDALGr5Xptd79kTg0VW1e3rzw72iqo4ZsG2vBV4PHFxVt6Q3ZPN/DRI0hPfxZo9b4yyApenNlXYEcDa97wyD/KhIVf0wyV70HrePd/9YDvodhKp6L/Q+86vq7kFzOq0/s1o+bq33t5bvSS3fj6Dt9613dhk/r6oVSf4T2HfU7eq8n940Clfz0MLZIL1Cf5ZkS3pTRnwyyU/p62k6gOnus/kk4DLgLuCiAbNaf/7NPG4X0HuNPh9454BZrR+3lnmHA/+S5Ks8dP8YdPj+uL6uNngepEQSSQ6jN5/WMnofgtD7kjnQ8NYkmwIzvWq+W1XLH279WbKeDXyH3q9vR9P75e2DVfXNQTP72rioBjwoQ2tJ/onel5o3dv9Ybg78a1U9Y8C8x9P7BXVj4E/pPW4fqaofDJB1CHAosG1VPbErXv5dVf2XAdu22vnparADPcwUWZ5QVUd1/0A/ZpAiS5If0CvSXkVfL7NaxwMCrCF7TvtbkvPozSc3U/A5ANi/qn57ju16BED15mYZNGNTekPj9+oWfZ3evjbQ635c94/VZP8WvdfVFwcZcpXkRHr72Quq6te7f8zPraom853NRfda2K3RULKx1vi98qv0fiT7++qm70hydVXt3rDJAxnC+3jLx61ZVmtdT6Ipej/Y/Vp6vcFPr6rnzXLTNeXtCfwDsGVV7ZTe/JKHVdWbB8ga2mfWXA1hfxvb96Ru//0PYBPm/lp44//P3nmHS1KVW/+3ZgQGkAFURL1KcCSISM7gFVGMgKIIEhUVRbgEA2YuKiYQDKASFImDRBFERHIOA0MUkAumq4IX5QIO6CWu749313SdPt3ndO3aLfNBr+eZ50zXObV6d1fVDu9+37V6HXee2UbR50rSbcARjL/fLs3gWphYc1Rj9KLATNv357Sti3sZYLrzjbCK32uKjNVqbJ9l+y+ZPEW/t5J8imqNhxl/f3wxs20lx5hiz9WzAaMA4QgjjIBCoHZ9239rwbGJ7YvUR5DdLZxg26Jfmyrktk1lnRKvt72Was6mkm72WD3HpwUKfa91gGtrbct2dpRUd2+blrhn2268O1syyCLpatvrNz2vB0/xZ0HSTd3B4l7HBuDZwfYJkj7Wp23ZQu2lMC/eH5rEMMj2/070+z6cN9heo9QzL+l3dDZ46m1rvNEj6afAh2yXKCVD0tF92pbTVy5BZCp097u52R3FIOk622t3XdPGz2mNbw6d721+YqH0sO3GJWql+/F5GX36t4eIfqSRXmX63lYHbqh9b7fkZulJupbQSzurbRC54JjVczyokDMuDGHeUKxPKtwfFXWmlnRo7eU04PXEvbdVCf42qPq3p7sddUhaY6LfO6MsuNS9JmlFhwxJzzbmtG1exryyGdYL8/JzNS9iVGI8wggjQOjhtC11eS1wER3tqzoai7SnIEHfHYyGJQ292jSXinwB+eOBXxMaRV8idt/uyOR6TNKCqT1ImkEtRX9QSDrF9taSbqX3BDhnUfOo7ceUBH0lPacX96CwPeZ6KEorczVP1q2CLIn7AeU7RN6oEGX/GWPLI5reH0WfhYT7Je1AmH8AbEtoGjbFwulna/27Id1r8+r9MZv4jL1UrXPNhB5Pi8vqmV+CdvqIa9X+P41wzc51Ql8M+LWk62hfSgZQl0qYBmwJ3JPJNRM4mXAj3JXQXPxrU5J+922FzPv3b6nvrq7pVkRZVBZc06lM2bBvJwTlc1CkHy/53A/pGkA8C2sRfTnEvXILsKukU20f2IDrMduWVF3ThSc7YTLY/qPGCuTn6jSXGrNK66FC4XkDZfukYv2Ro1xxaUnzl8g4s71H/bWibPakJhxDfK4uVzjynkWmNl/XpseYX5Fn0jWRu3BuWXCpe+1jRBZtrzY2apukK2xv1OP7y/rehnAdIMqL32j7vIxz620rPrcs8Vw9mzDKIBxhhBGQtDrhFnwtYwfDPRvyTAG2sn1KgTa9dqLf55Q0lEaVJVJlEyiMGS633XgBJ+mNhFj2SoQw8obAzrYbaWtJerHte1Nq/jjklB4pdDseJBzE9gB2A263XUTcOy18b3OGG2PKxtgAuC4FgpYgMsQmdOfuw3V0j8POySoojXQ9DwXWJyZMVwF72P5jBtdUYE/b32rZpuL3Wp/3mSfuj9KQtD2hMbcGcCyRVbRvif6z9h6zbTd2k+7X/5bqd9NYcYXtDTLOnW17zXoWV05mS7/7tkJmX/ly4EjinnsA+B0hBVCs5LOendjwvCL9eMnnfhjXIPFeBrzV9sPp9XOBnwNvJrIIB+5LJH0CWA7YlNASez/wY2e6nEo6Dfgm4dC5LqGpuZbt92RwzctjVtF5wzD7pDb9UTr/OOCVROCsrTN1N/d8hPHDwGZYQ3yues1HPS9kb5dEyXst3Vvr276ydcPmcaSg48LEOvJx8oOXQ59b5jxXzyaMAoQjjDACkmYBVzBeN+LYDK7rba81+V/+6zCM8pnEO8v2Omkxshvh+DfL+dqNzyeyQwRc43Yl32+x/YuuY7vaPjyDawpRSv3G1LZf2v5Bi7YdSmdXcAphIvF72ztkcFVBljWBY4ggy+dtn5rbvlJQaPO9C1iGseY/X8rg2rB7gtnrWAO+WbbXyTm3i2cqcIHt17XlqnEO4/6oB+Ea3x/9SoQq5JYKSVqRKHURcKHt3Azk7jZOITKoPuJ5QKagGwrn55/bfkXGudfYXk/SL4FDiMyf02zPKN3OppC0rMMgZmFgiu051bFMvrpMQXVNX+uMstLS/fi8DIVb5attP55eLwDcbHvFnACrpE0Z+72d36JtLwC+A7wh8Z1HbNg0likoBUkTBjubbhYnzv9v7rc2/VE6f79ex52hvaax1TNTiE3jU2xnmVpIWpKx2ndFZCNyocJyHRqShFBJ5G7qdHGU/t6Ky6aUROm5ZddzNZUI6Gc/V890jEqMRxhhBID5bE8YRGuAC9KO+8mM3UltOnj1K5HIcekbRvkMFHT8k3ShQ7z75z2O5WBfSY/avihxfRJ4HdA4QEhkqn2HcIer2rZXOpaDusPtE0Q2Rlagy/ZMSbOJIAvAO3KDLCqoS5RwJkn3ioxy8S4cSgS5Jjs2KK6U9F3GP6eNAl2O8qqnJC3qcoY/w7o/RP79UbyMSdLxtnckZAq6j+Wg3sYngN8DW+cQqbf23SNNMwF68Cn9/AvwqRwu4MsK9+2PE8/AdELAPLdNpHZRb2PmZz0dWMN23QXyNGIDIwf1cvvqmuY6Lxbtx9PC/ADghcR31vh7G9I1gChDv1bSmen15sCJKXB7exMiSQfY/hRwfo9jOVjBXXp1kjYEGvdxBces2U3fewCUvt+K9UmF+6NsE4Y+OKj2/yeAP9j+Uw6RpK0J1/NLiM96qKR9bJ+Wyddzfttw07O0XEdxCaHS4x9woaR3AT9xfmZW6e+tuGyKChrMDWFuWey5ejZglEE4wggjIOmrxMKjW8cmR3i/V6aEm2bVDatEYl6DwuhkIeBiYGM6g/V0wh11xUzeFxA6O/sQZVUrAts6z231BttrdB3L3hFNi7T/s/1kej0VWMB2lg5mypzaiJjUXNkim+tdtZdzdYlysicSX2vBZoXj5QbA3kC9JHg6sGVudljJUqG0CF+dWEDXg41Z31sJzOu74zD+uUrPwa3OKKUeJqSO9t1ot703Uiboq4ADiT63wnRgH9uveloaVsMQ+vG7gc3bZL0OA+l+fSmwJCHVATEuXN//rAn5en1vbUxKevGNOzYgV9Exq8a7UO54XOMomRWaTAAAIABJREFUer918cwTfZKkb9veW300s52hjyhpWeBe2/+XXi8ILGn79xlcNwObVlmDComNC1rMGz5eezmN0Pa8o8Um6jyPEveaOqW3TwL/pP0GyDwJFTSYS3xF55aSXpTaZEJ2JstJ+tmAUYBwhBFGKBbUGxZSsHA52xekydJzbM/J4FkeOIyYbK0saRVgC9tfzmxX6xJSSXsRwZ+XAH+mEyD8O/AD29/NaVvifiFwAbFT+P6mO5eStgW2I4Jvl9d+tQjwVG52o6RrgDd4rD7Uec7TI/tPwozhdFKGGHBq7jXt4m6rS3QkcKjtW1u04bVE4HhXxmZ/zgF+ZvuuXO5SkPTeXsfdUKKgZNawOo6+ApYi9OBEiI//t+1lm7StxrsQITy+lO0PSVqOyAg6e5JT6xyfAT4LLEjHHErAY8CRtj/TsE1DkVDo8T45pZnDcJhcAtiF8f1u9iJV0kbEGHN02lxZxA3KgiW9neh7tiCyySvMAU6yfVXD9tTL7MehyQJpiP34lbY3nPwvB+ZrdQ26uFq7M0v6CCEdMgO4u/arRYCrurMAB+AbymZP13u0HbPWB44Cnmt7KUmrAh+2vVsDjqHcb33eq1GfVLo/krSm7dkqq1l3PbBBtZmrMNW60hnuwd3PQbo/bm77bNT4FiBKxzfOOFeEsd+ytveXtBTwItuzMtuyJPBV4CW23yJpJUL776gcvh78RYLbLd5/KI7Ipa9DF/fLgG/bftekf9z7/CJzy8T1QaLC6yJivvVa4Eu2f5TTtmc6RgHCEUYYoShSRtxudDK6LgcOr3ZDM/h2IVzAnmd7RlqQH54zyZR0KZHdcUQ10LfJ8pJ0Lp0S0rkuhLYnKkfsx7WH7UNz2tHF0106Mz+RTm+al38tDSxLCLPXd07nALfYfiKzjTfZXm2yYwNy3Qms2rXbfpMLCA+rvS7R7cArCKOCR8kIdNW4lnZZo4PnA/vReU6vICZLOc7I1SJmxcR1p/MyVYdhGvED4Azb56TXbyHKjD/clCudfzLxvO+UNhkWIoIFOffu15oGA/vw9NS/quA8Hawi2nd9MlVrTcvKWL2KGFe6+93Tm3Ilvv2Iz7eC7eUlvYTYZGgc/JK0vu2rc9rRxdNzYVShyQJpiP34d4AXAT+lnYNu0WuQ+I4Fvmv7upzzE8eiwOL0+N6cV10x9M2eAmPWtYRO61m5c6Qh3m+t+6Rh9Ec17tZjYOLpNT+6OSeALOkbwCrAj9OhbYgs9U/mtK0H/+JEJlaOluxhhOb5JrZfmbjOywmEJr5fEGaLn7O9qsI1+8acYGip8a8H59x1ke2fNjz/yLQpWdQopvR16OLONpgrjbRW2KCa46Y58FUl1grPRIw0CEcYYQQk7dTruO3jMuiOIyaCVbBrO+B4IssrB7sTKeHXpjbdlTLjcrCQ7VkxZs1F1mQ14aW239zi/DqekrSY7Qdh7sRrW9vfb0Jiu5jeYgrI/EFh9HBPVxDupURZeg4ekbRGteMpaU2i7CIH9xClDFUAegEiE7MxegRXW+kSAW9pcW43fijp3V33x0m235TJdxJwGZEBC7GDfDIhnN8Ikt4KHAH8hvjulpX0YXeZ5EyGkgHQGtazvUvtPX6hcNfMxQzb26QsGWz/Q10dSgPMUk1fR9JiwMZNFw45AcABUET7zgXNa2pYyPn6b72wJVHGdAOA7Xsk5faju0q6o+s5PdgNsxtzMiQm4BpWPz6dyIB9Y/3tyND8ouw1gHAH3l7SH4jStMabM+m5fCgFQv/XqWpB0nRJ69q+tkmDUjbZpZKOKdXXDWHMwvYfu7q0J/v9bZ/zh3W/te6ThtQfIeltRNC31RiY8FdJW9g+K3G/HcgyrLO9Ty0wBZGhfkYOV2pLPct/KrAE0Nh0LWFd22tIujG19YEUZM3FC2yfosjQx/YTkhrduzWU1H5F0veJzeIqULurpE1t7z4oh+0PpZ+l7+Fi10G9DeayMhsTX1UJMgbOq267n1ibVpiTjo3QA6MA4QgjjAAdhzOIYMvriU49J0C4ctdu0cUpkyoXj9p+rJqwpl3B3NTnv0maUZ0vaSvg3hZtu0rSq92ihLSGXWx/r3qRBuldgEYBwgqStgQuaht8SDiFKI2q8CRwKmPvmybYGzhV0j3EZPpFxM52Dh4CbpN0PnFdNyUCL4dAs1K8ksHVxPcHgBTQntaS7gVV0CFxP9AiUA7wYtv7115/WVLuNfgm8DrbdwOkZ+znQKPFkcYKg4/5Ffl6PfdI+jxwQnq9PRFUzsVjaaFb9SEzyDeg2a++WLP9YMqkappZUNyB1PbOTc/pBQ3HYfJsSW+tskIL4DHbllRd04VbcK3S4zltXJamPppmNd7G2mYU7sdL3SMJJa8BQO7GSS8cxlgzqId7HGuCYps9pccs4I+SNgAsaT5gLyBXY3Keu98kbWL7on79UmZ/BGEQ1XoMTNgVmKkwEQP4E9BzE38yKPQMz6k+l6QFJS3jDD3DhM1q/38C+J/cjFDgcYXubvXML0FksuXikZQVVvGtR8wPG6Nw3wZhYvZK21XbjgVuyyGS9G5Cn3xOmtusAexv+8bMtpW8DsUM5hLWqv1/GpFoMqG+dDfUkWC5m45xlYmA7y0t2vaMxihAOMIII2B7j/rrFEw6KZPuBknr2b4mca3L2EGjKS6V9FlgQUmbEuXLP5vknH7YHTgSWFHSn4nSzx1atG0j4H1pl6tVCSkwVZJqE4ipRHlwLooEHxKe41q5TArYZrfN9nUKUf8qtf9O249n0p2R/lW4pCmBhqftsgWxcHgJcB+wNLHYyjEseErSUrb/O3EvTX6gHOA8Se8hFnEQZWW/zOSaUy2MEn7L2J3agTCExS7AtkQpdXWPXJaO5WI/4FzgZZJmEiYI78vkmtLjWM68rJgDqQpq3yUUd5gkghaflfQYUPUbuQFkgFMkHQEsljZl3k/NebUhpkha3PYDAAqznJxretDkf9IYRftxSS8lKgWqMuDLgb2c5wxZ5BpImm7772T0PxPRVuMygO2n0kZlLlpv9gxrzCKCU98B/o3IxD+PmDfloMj9JumTtg/s1zc17JNeS2iQ9eqXcvsjKDQGAtj+DbCeQpsZJ63mTJxKgSDtBM/VdEm5pl+HEOPyCyV9hZiD7JvBU+FjhP7rDElXEtmNjSqXhjD+Vbib0EKuModfxlhd0ybY1/apCs3WNxAu1YcTWdM5KHYdSma+J77uDL9vS5pNaAkOimpe+Zv0r8KZPf52hIRRgHCEEUbohUcIDZmBUSs9mI/IrPvv9Hpp4Nct2vJp4APArcCHgXOAH+YQ2f4t8IaUnTDFGUYnXShZQnoucHJaJEF81nNb8JUKPkDBkpd0/u7ATNu/Sq8Xl9S4nBqKTUg+Ruhc9tKONLH7m4P9gfUI18DVJb2O/ID054ArFDqaAl5DtLkRNLYkbW+i/B+iXOhh4BMNuKosjOslnUMEG01Myhtrf1WLEPVxIM5ZhKRz9mp63gR850u6gbiuIgIiuc/C9ZK+CVSZw7uTEewrPClvs5kzDkPIxBhGpu9BafPp78SmxX/aPj+T7mDgakmnEvfHVsBXMtrU2NxgABTtxwm9rxPpLMJ3SMc2bUpU8BqcSGQ5zabTz819GyCnNO23kvYksgYhNil/m8FTocRmz1DGrNSXNTJfmQCl7rcqg7F132R7v/SzdL/UcwysxsgmmYmSvgoc2JVh+nHbn89oV6lNgeLPle2ZKdjz+sT3DrdzRL+NCACvkPjupPc8eCIUHf9qmeCLAHdImpVerwvkmoBUZdNvI0rGfy4p25SvxHWQdIrtrTXeaK5N0gRdGyCVFmSjdYy7JFgKBd6f8RiZlIwwwgjd5UxTgJWAU2x/uv9Z4ziKGwyURsqM3Inx7pe5u4IoXP5ek15ebvvmTJ4pRFCwMl85H/ih7SwNFUk/Ah5kbPDhebbfl8E1A5hJZMIJ+CNh0pC1A6reItxN3Qj7TUgAaDohSd//+i3LIbo5r7e9lqSbgdVT5kmW2HjiewERmAK4pkVgqggkHT3Br+2GumuSzra9mcY6ENf5Bl6ESPq27b3Vp1TTeSWaSNqQMMF5RNIORHnPd3L6t7RRsS8d3cfzgS/bfqQhz1A+awlI2sH2CerjtOxMh+WUnfvv6eUlbuAiPWxIehVQ6URdZLuxxMYwFlz/on48y2yqdv50xo7NOZlJRZGy+w4hgm4GLgT2tn1fJt+biUqGMZs9thtlcA9pzDoQ+DKhCXwuYXDxUdsnTHhib65i91uqqDjA9sAbWJPwLUDo7y7D2PstS0+v5FjYay4k6QbbjUvaFdIrh3YFafd0QSfpXEg63vaOkx1rwDfuO8r93kpBfdytK+RsBEk6m8ju3ZSYf/wTmNViXtn6Okh6se17+60Dc9d/GmvIUmlBHmT7zgyulYkN8WoD+m9Ef5RV6v1MxyhAOMIII3QPYk8Af3DDMqFhZP8k3s2ITKyliYlcth6Zwv3yGiIbca7GRm4GjqS9gF3olKVsSezotXYjbotSwYcuziI7b2nRu4o9ppz6FtsDl94OY0LSNEg5AN8FwDsIN8cXEGXGa9veYMIT+/MtDixHTc/Q9mUNOYZVmjZPQdKatmf3m6DnZmhJugVYlVg4Hw0cBWxte8KFwDAxjM+aJua9go2NMpMUQv1HqI/Tcvfu/oCcXydK5GamQ9sC1zvTETpl+RwAvJAYX9poXlacY3RHq2yxBucPZcGVuEv14xcSz0AlvL8tsHNO8EHSh4EvEmZTT9G5BjkZf0g6K7XrTNv/yOEYJkpt9gxhzLrJ9moKDePNiEzFy3KDD4mz1P12tVu4yHZxnUvo03U7offKyPyXIo0xa9t+NL1ekOjfGkuTlArS9psvVMiZN3QH79I88FY3dLyV9CKiJP4EwhSx2licDhxue8WMthUZ/4YBSQsBbya+q7skvRh4te3zMvlKXYepRLXMUIyA2iKt/z5n++L0emPgq7nz8Wc6RgHCEUYYoQhKZv908d4NvJMYsFp1WKV3E9NEbv0q6JaCcldnZncsRwSSVmLswjLreysNhUvfqxjbttzd9m8QAd96OfUfbX88g+sAdzma9jo2INdBwNXAT9rea4lvYWLBK6Jsa1GitLqxc5qkDxKlsi8FbiIWl1dnBGyOtP2hrp3ZCm7Cp7LaUHXekgY7KEqqlk8v2+hdzu1DJP0n8GfbR7XI7lgC+CTjn6usRUi63/5p+6n0eiqwQE6AROEsXmEakW3zhO1P5rStJFK/u1rX57wxp99N598NbO525W0VV0/d0ZzFfeJbFrjXY51gl3SmyUDhfnxpQoOwCtpcSWQnNQqGJq67iLG0SFZ0CpZvQ5ThXUdoKp9dfY8NuaYRMifd31ujDOkuztabPYmn9Jj1K9srS/ohcJrtc9Uu673k/XYYEQg6lZDBqfga6wZWnzOnHX34SmZeforQSKyyEncGzrJ9YIv2tQrS1uYL04gyz5uJec0qRPBy4MCtwmX4s8CChAt6tVZ4jNhgb7TRI+m9hA7wWowtEZ4DHJN5fxQZ/yRdYXsjjTdga5PoMAP4k+1HU5BrFeA413RNB+Qpeh0S54XAO6u5W1soTGf2I/TeDVwBfClzDj2uH2vTtz3TMQoQjjDCsxjDGLxKI01MXl8tBltyfZTQWTubmvOo87MbbyV2eqvF2zTgOtuvzuC6ghgIv0VMDncmdBKbiPHW+YoFHyQdDixElM39kNDVmmX7A5ltK1ZO3Ss4I+mWzCDtHGBhIqPgnxR6DlSgbK6614hsk9UUJi9ftT2hS+wwIWlz2z9LE/RxcH5mbusS9Np5GwPHEqUpIsTB35uzGE98lxILwJ2JEtf7gJszn/nzgJMJ3cddgfcCf80Jbie+a4A3VIvAtCg8r9QOuaRZttfJPPdYQq+xrqt1cE6QJQUIN66eI0XW+iUtAoRX2t5w8r8ciOtmohR1jO5oi77yemADJx2xFOy+0nZjJ9jS/XhJpIyud+YEsyfhnUpcj12AN2cuyE8ldJS3A75EbPbcYTtL27TUZk/iKjpmKbJz35G41gEWIwKrjQ0QhjBv6FXG68w+5Eii9PbWnLb04CuaeakoQ59b/eHm5efDknb4CWGAd2t6vTLwBdtbZXB9LScINQHfu2yfXoqvB3/2+Fe4HTcRwdBlCE32M4FX2X5rJl+x66BwCF6dmNPXg/i5m8XnE8ZyVaB9e2Lsf0P/s/pynQHcQEd3ewdgTdtb5rTtmY6RSckIIzyLYXuj9LOY6LvKl/Z8EjgnLczrQb2cCc5jhOPX5+gERHOFyyF2eK9NAw/ExPqoTK4FbV8oSY7ysS+ouVtXHTOJ4MNm1IIPmVwb2F4lBd6+KOlg4BeZXJUL5FHEbqCJrK5GwUFJHyHE4mekgEGFRYhMlpx2FTU/UJ+yOfLut/+z/X+SkLSAo0x4hclP69u2W4jn9BSHa2Jj2K7cxG912dLkkgY7BwNvdNKskbQ88bnXnPCs/tiGCBR8wPZfJC1F9Ck5eH7KQNzLUQZ8qaTG5i41TKtniNh+WFGO1BgaKxUxhfi+Fm3RtlU83rk1tzTya8CNafNIRKB2YL3cHrhe0smEw3t9jMlxNH3c9v2SpkiaYvtiSd9u0baSzsNF+3FJLyccb9cj+rWriaypHAOPzxDmZtcy9hq00QdekNhs24bQ6so183mF7XdLervtYyWdSDg252IvOps9r6s2e3KISo9Ztj+dsuEesv2kpEeAt2fSlZ43tDYWUUfT8znAzpJ+S9xvrcwU6IxPbwNOtf2QpIn+fjLcSBj+Of2/KRZOP4veH8AK9aCq7V9JemUOke3PqFAmbcKFCtOvSpv2UiLbrHE2W+nxT9IhwI9tX53LUcNTtp9QSGMcavtQSTn3CDD3OpTS9P0J+U7gvfBi2/vXXn9Z0jaZXO8n5uNV+y5Px0bogVGAcIQRnuVIO+y3OUOnow8OJibkX0+L3ezSnoSvEFl/04DchVGFjxOT/SJlTLa/KekSIv0dQn8pd6B+NGXW3SXpPwgR4ue2aF7J4MM/089/SHoJcD/w4tyG9crqktQ0q+tEYrHxNcYGB+a0yAitSoGXtb2/pJcRE5Rct7lPACsXut/+pCi1/SlwvqQHgDbGP9XC+RRJTxHB5FOcUR4IHKzQAToNONnJnboFirj7JsznmqC17f+SNF+Lts0hTEmeTMHGFelosDVFVep8r6IU7x46Ato5eETSGlWwVlEm9c9JzumHulvlE8DviDLLXEyRtLjtB1LbnkfmHNT2j1O/W2XRfcr2X1q0bTpRYvXG+tuQt9B5MGVuXg7MlHQftUyKDJR0Hi7ajxN98PcI7V2A9xDPQuNsM0Ju4iK69IFzIekUIgPuXOC7wKXOr0KontMHU8bUXwi9ylwU2+wZwpgF0actI6n+fB6XwVN63rA84SS9pKMMehVgC9tNHFw3y33/SXC2pF8Tn/kjigqOrDmvpK2JTadLiP73UEn72D5tUA7bR6SfjTVeJ8EtivLzekbXLRP8fV+oTyYtmQ7cxOb8r4Ct0+sdiU38nCqL0uPfbGDf9JyfAZxkO9cx+XFJ2xKGi5unY9lzGklfI/rKStN3L0kb2P5sU660gbIgsJQzjER64DxJ7yHcwSGykBtl09ba9gCQveH0bMOoxHiEEUao0sL3yAwO9ONsXdqTeIppxShK+t5Rsowp7YC+jLElpDmCzWsDdxAlPfsTu5UH2r4ms13X2F5P0i8JB8Z7CE2hGRlc+xJaU68nFoQGfuD88ufZwHbdWV22G2d1qZAeS+I6jFicbmL7lenanueMcr7EN6yyudcS98e59eyiFnzLEYY229uemsnxImJivg0RcDm54cKtzlXMYEdRlvYkYxc0U52pH5bu3dcAixOZqtcBj9nePoNrMyKQ9DLi+ZoOfLEKBmXwrU1syNxDLGxeBGxjOze4WgySdiL0jk4l2rYV8BXbx094Yn++VRjvQFoycyELKqg7mvhKOsGW7sfHSTkoU9NJ5c023kSUeTeWrejB9UHgdDrGRM8F9q2CMBl8ZxASBXsTc6QHiI2MxiWCQxizjgdmEAGb6rtzTibnEO63S4F9gCOqe6XN/FBhvFFpm13ZNgs+bXpUmZcLAdNzNi4UMgWbOrlkp2DjBZnPVYmgap1vGvAROtlmlwGH5SQAqLBsiobgql4a6R55F7GZspTt5TI4ViKqgq5Om2XLEkZpB2S2qZimr6TNgYOA+W0vK2k1Iotzi4Y8leyViGzYp9LrqcDDTdaTkn5GD43sCk3b9mzBKEA4wggjIOkyQjdiFmN1I7I6To0v7Tnb9h6ZXAcSk6Msh64urjMITb6LKVDGJGl/Qhz5N9RKlt3C6UyhV2fbc3I5Ek/R4EONdwGilDFbhLjPwjJXN7CYHos6BhRzF6u5C9507uqkMnRa3G8qn+Vb8S5NPKPbEIvBk93SxVHSqwlZgG1st834bY10v+5OJ8v3cuD7Tg6RGXzVPbIHIQtwYJt7pDRSdmSVjZRtyCJpdyKwVdcM3Nb291u07VWEHhmECc3tmTw/IoI1t9HJNnOLoG9RE4oULF+HGBOua5ndWHEWcYKt8WX34+qU332KCG6dRHzWbYDFnSdu/1Uio/xnlNEHno+xgYxLCUfTbIOiYaDtZs8Qxqw7gJVceHFYaN5wne21uz5rVgBIYTL1bjpZwu8gSoNzA2fF7jdJt7qmaauoLMnVuS0aVC2J2vW8CVg3bfLe5nxDp6uBfWxfkV5vCBzkDOfrYYx/iWcdop98O6FjuvkkpwwdKqjpmzZQN0nnzxP3W+pj+8JRYTVCF0YlxiOMMAJEtk4RqGxpD8Sk6xOSHiVKfdqIcP80/SuFrYEZhbK41iKCSYuk1w8B78/N/nFHQ+QhOovy3LZNI/T+5jqJScraNU64XuPLVHLLLUrqsTyegnGGuTv3be7dImVzKSPhTklLuVCWr0Lraz4io+vdztMNq7heSUx630WUkZ1MlPPn8hUx2EnX8uYUWM0SZe9Nq/WJe7YqOeqlmTgIUaXhtj5xf7TRcKuwAh0n9DUkYTunPHAX21WJd6UZuAuQvUCyfZukv6a20eJ+Xs/2Srnt6IHjCROKN1EzocghStlm/0k891V54Jds/yi3cao5wSrpmjnDCbZgP14vv4MwnKpgQk+wKbZNP+vnttEHPozo36r7dcd07INNiRROml8ANkxtuhzYPzcrNHF2Z6/lziFKj1m/IjKP723BQWpL6XnD3xQZtdVn3apFO7cHVnXHYO7rRNZkVoCQgvcbcK6i8qOSrtiGfO3GhWzP0lg9xCeakkg6xfbW6mg4jkFOMInysikfAY6VVGkFPkBob+eg6PiXEh22JJIJTib6j8ZVLolrOUJapxrnqzbm9pUlNX0f93j9zez+SGovoTAKAOZhFCAcYYQRsH1pyiZazvYFivKIrFJDQgdkWxco7UltKyay7ExX1QnwK6Ik+L4CXD8CdrN9OYCkjYiAYaOJl6RDmTidPidb8jhCe+3Q9Ho7YlH97gwuiInc7nT0QC4nP/BQUo/lEEIf5oWSvkKUQbYJns9nu6eLYAYWB26TVCTLlyhTLKERA3Hvngy8yfY9BfiKGOwMI7BKlAV+BjgjBbxeTmQk56CkhhuS9gM2JhYO5wBvIYyAcgKEUyWpyiRKQYjsjFCFCPrBRKnsfcDSRBAuJ1vkakkr5WYg9kBJE4p9gNWr4FEKLl1FPCONoT5OsJltK9KP21428/3/lZxrd2XRXaQo3czBSUQp5bvS6+2J/qmxkyb0zF47WlJu9lrpMesFwO1pnKlncuaMM6XnDbsDRwIrSvozoQvXWNoh4R4iuFIFKxcgdJ9zUex+s71P2vCsst6PtH3GROdMgFJB1cqxu5iGozvusV9IwalFicSCXNxqe1VFFQ62/96Cq+j4RwQG13cZPeqjgf2AbxHjws5kblJCcU3f2yRtR3x/yxFz/Kty20asC54ishL3J/Tov1dr68DoE9x+iEhO+HKbDZ9nIkYlxiOMMAJpZ+xDwPNsz0gd++G2X5/BVaTUQtKKDvHuNXr93nk6f0V33lLW35lEoLDVZFo9NJiUyoca8ky4Y5oTJJV0e3fGTq9jDfgWJoTan0yvpwILOEOrT+X1WFYkNJMEXGg7K5MocbUum1OI2D/ar0yi6e6opB1snyCpZ+DSee7gRSFptu01VSs7VypHyuAqKp9QEiqo4ZbOvRVYldAPWlXSksAJtjfN4PoGEcSrdNY+DPzRdlZmaFosb0LIRawu6XXADrYbC7+nZ+EswiyitQOppFm210n3ym6Jd1bOuCDpKqJc67H0en6i3GqDzLbd4o4T7CqKUuNf2H5NBlfRfjydvzLjx9PGAem0KfkxQpfrQ2msXsGZbpqSbiAyo3+TXr+c0OBtNJ6mc8eVyKmrDLQh352MzV5bELjJdq5RSckxq8g4k7hKzxuWtf27NH+YYntOdSyD66dEkOF8ImiwKTFG/AmyZEBK3m/LAvd23R9L2v59BtfLiaDqBkRG3e+IfjeHayrRf7eqSKnx9TLkmtN0rVDj+28iwHgyIWGRHeAoNf4NaR1TzY/m9kHVsaZc6dyfEZuVZzlD57mLayHgc3QMv35JZEu2lXRpLaGQsjifJD4rxIbsQsR4v5HngXLveQmjDMIRRhgBYmd2HUIrDdt3Scp16StVavExImjZSxPN5DmdFd15I5x4D6BFCWlt4nCppCOIDKJKz+mSpny9AoAK/ZQHW0yYbpC0npNhiqR1yS8JBriQyL6o9LQWBM4jJrGNYPt2SZ8Alk+L1TtbBAePt70jUW7YfSwHJcrmriZ0PD/Yoh11LJx+ts7MVf+yo1YBG8q6+xaRT5D0bdt7q4/gdZOAY21h9AtJn2ashts5LZr5T9tPSXoiZVHcR2iQ5uBTxKLoI+n1+UQGWy4et32/pCmSpti+WNK3M7mOIsaVIo5RduF5AAAgAElEQVS3wJGpj9yXCDw+lygTHhi1gPvdwLUK4y8TWlNZLp8JJZ1gi/bjKpuxejRRulyNAX8m5A+yAoSEg/zFkn5L9EdLE+N9Doo5aSYUy14rPWblBAInQOl5w+nAGl1BjNOAnMDIGelfhUtatAsic7jU/XYqY+dCT6ZjjTfIHHIVb6gHVTPbVGXkPyVpUbfQkqzhBmJ8eoD4zhYD/iLpf4gS36byOisSGY67A0dJOptwC74io22lxr9hrGMeVehS3iXpP4i+47kZPBUOIuYeX5d0HTEfOdt5UgBvs/05IkgIgKR3E/dvDkpKKLyhK2B/ay0AuUMm5zMWowDhCCOMAPCo7ceUdCMkPYcJylQnQZFSC9sfSj+L7FYmLGj7Qkmy/QeitGE2DReDNfzD9iEt29Q9cdiv9v/G10BRvnRK2rVcgNCuWQ14QtJ2ti9owFUFfuYDrko7tCYmv7+e6NxJMM01sX3bD6edx8ZQOBcfS2TqCXiZpPfaviyDbky5Y5qYZO3KQrGyufkVJRsbKMqOut+jkXOr7SPS5/q77W+1bFvxsqOELyt0hD5Ox2DnozlELiefULntHpTTji4MQ8MNQttzMeAH6T0eJgLMjeHQjT0s/SuBB1P22+XATEn3UcvobIi/uqXZUh22q4XfpeRr3lUB99+kfxXOzG1Xwtnpmn6DWFCbuL4DY4j9+FZ0MlZ3VspYzeSaYXsbhVwEtv8hjRWyGhSpf1sVWI6xhj2Nslg01klzbzqfbQrxbH2iIV8l//EQUYrXnb2WgyJjlqQrbG9U+8xzf0VDzefS91vKkHwVsGjXGDidWuZqEzikBOYHlk+Hsg2dEt+FSlmvNb6srCngOa5pUqa5eRt5hyIapgkPE4GV8xmbkZ8jXXM+kWX5y9TONxJl/EcTCQaNpDYcFSinAKekDZ/vEH164/G+1Pg3pHXMXkTm255E2e0m5GstVpsCl6a+YxNgF0ISI0fn/TOMDwb2OjYoekkofD6Ta6qkdZz0CyWtTefeaKzL+UzHKEA4wggjQAwOnwUWlLQpUWb1s0yuJyXN6Cq1yNYjVNnSo9I7b5dL+hqRdVIvIR24bKDwxAFiJ3D/9P/3EhP8JYiJ8LHAwAFCygd+KjwiaY3qe5K0Jp1MmaY4GHijk56epOWJLMyBF0mSPgNU93+lWyPgMaI8JwuF7t1dCZ2lxehoLFYwHR2rgZEyAbYlMmmzYfveNKk8puR97LIGO3PlE4AZwL8BhxMleU3aNDv9vFRR8rWUMzUcCwWOe/Hulv57uKRzgem2s7LXVF4I/e1ExtTexP28KGEIkoMbFTqB3aX7jZ8FGJP9V8dDwGzbNw3CYfuLOe89AG/Vl5+eMmJynGCH1Y+XzFh9LD1XVabIDGrXtgmq/i1tgGRnb7qg/nFClT03m5bZa6XHLNsbpZ8lPnPp+22FxNk9Bs4hghmNUXhjscKawDLE2no15RtE/VXSFtUmiKS3A1nadSqrYQox3+juZ3OTCdazPff62T5P0kG2P5w2txtDUSK/DfBm4nnbOpOn9PiHpA3o3B8VX+P7w/Z16b8Pk5+l2t22BYlnaxuiYqWRFJGktwBvBf5NUj1pYjotgm+2Z6YkjkpC4R3Ol1D4IPCjtFEp4O/AB1J27ddy2/hMxUiDcIQRRiAFzT5A6EYI+KXtRlkKNa5NgGOAMaUWtrOE/CWdTEyod7K9cgq6XGV7tQyutQlh/MWIINqiwIFVCUwGX6/PZDdwW1VhTTiN1eo4HTjP9hHpdWNNw3TeUn3almX8kK7DSUSplQjXxG0ySkr66biNOzYg19ds52Zv9eIree9+wPZRBdv2LSLD42TGZgLkaOJcCLyzUNlRtalQxN1X0k0k+YTac9FGP2xzIotwftvLSloN+JKblRhvYvuiXhmhkB/oStyVuL2BK5wpbi/pCjpyDJuT5Bhs52ZbI+lFxLUwcJ0zhdAlHd3jsG2/P5PvRGAtOptimxGBpWWAU20fOABHsRL0Lt5xTrBAIydYSdNt/129Nb8aaaJ28X6fCFK9h8j2fZjQ0mu8aE0bk58nFuTnEY7B77N9SWbbivVviW8LOrrKl2RuUNb5imSvlRqz+t0bFZrcI0O839a3nZUR3YNrNrBd98ai83Xcjic2oG6isyHunMy6FByfSWxmAfwR2LHadG/IVUzDNPHtZfs7kx0bkOs8Qm7mpHRoGyKb9s3E+NBUf/v3wI1EFmErPb3S41/h+2MtooR3acYGG3M1eE8hxuVKv/HSlEHZhGNVolLpS4ytyJoDXGz7gcy2zQD+5NDh3pgwbTzOmQ7QiXNRgFLz1WcqRgHCEUYYodign7KJ9iTKA0qUWiDpettrqYBI7bwISR92lH3u1+v3TTNTJF1D7JT9D3AnsKaTiLekX9teMaONVcmQiN3UZYnrmuNAWnHOx9h7JHeBdDQx4arKv7YHpuYECyRtSYhbP5ReL0YYDvw0s22t791hBZNKBLdrXGcSRiAlyo6qe/h7RCYoRABiD9uN3X0lXWt73eoaKOQTbmgxmZ5NlOFckhtwlPRF2/sNIdD1feAVdL63bYDf2N49g6u0EPoHiYXDRUQ/8loisJrl7lsSCnOStzrJHqRF9M+JhepsD2CqIGlN27NV0OQh8Z5CLLKq/m07YDHbAzvBSjrb9maSfsfY0vbUtPysmNp7LEOLjNXE8XxgPaJ917iF42fh/u3rhP7bzHRoW+D63MCcemSvAe91RvZaqTGrz71RodE9Mqz7TdKxwF5VcEBRRnpw5lhfbGMxnXsHsJILLqpTP4RrciwZHNX4dw3wTkLD9Dbbr8jkG7fJrB4GewNyvYAIwlVuzVcCXySyt5eyfXdDvulu51xc5yo9/hW7PxQmR/vQpcHrkE3K4XsTYT6TXelV45oOPOICBoTp/JuIzbtliDH5LOBVtt+awbUocb/VDTS/NAoU9saoxHiEEUaAKEXtDga+r8exCeFCpT1daF16NMTsjtYDjstqwkHok5xGlBV/qxYcfCuxu9oY3cEPhbHKbn3+fFJovNP1JZKOyAwS7kqIUlfBqMvpGOQ0xX6uZVzZfjAFbrMChJQpm3stEVTp5bCWVWKc8IYSE8KEXmVHbbCQ7eNrr0+QtE8m16UqJ58AYbbxkMbKozWa9NveL/0sUh5UwybAK6tFSFpQ35bJVVqOYR9gddv3p7Y9H7iK0DpqBEkvJbQpN0yHLicCB3/KbNsLGftcPk64hv5T0kDPq2sl6Jlt6IeVuwKUF0u6vQlBCtYIeK0zs77rSOPVgrWA6nrA/On/i7iBEUKdy2Fic3/iWikFHrJMFVxWuuOtwGpVVk16rm4kXyu0tSxGDUXGLBeUPbC9WWnOhFVcyxyy/YCkxoGphOsl/ZCxG4ttDFR+RVRC3NuCAxg/r5TUJpDRS8O0sdmGQpJkO2BZSXX91+lAVkZo2gDYo8+vGwUHE14k6Qyi715Z0irAFra/nMFVevwrdn9QWIOXGD8/I6mEjNN5FDIgTHjK9hNpg/y7tg+VlLWOIeYav6JTdr4joXfZc/P92Y5RgHCEEZ7FmGDQX4TMQR+4UtJ3KVTaQ0yUziU0YmaSSo8acpQ0GKijyIDjQppwietaws2t+/g51BxSFXo7jXRGalw3KBwJc1HE6TotLm92ZEU2KsXug16O1m3Gydb37hCDSXcpStB/5HxNlwr/C/zcDctSuqHhuPt+mpBPuJUwAznHmfIJCbcpTGOmpon0nkSgqzEkfZWQOKhnxHzcdq4I993AUkCVSfAy8hZaUFgInchcqQd75qRjOTgaOBGosuh2SMc2zeSbScd5GCIYf6JCm2igYJzGO3mPQW5mEoWcYG1b0s+BrNL6LhxA6A1Wpdc/JsbBaUQQ4lMtuE5MXAtmcM3FEJ6txejMiRbN5Kgwn2v6pbb/K22a5aDImJU2/foiZ/6WgtLbA8va3l8hVfIiJ5OADEyRtLhTuWIaL3LH549QbmMR4AXA7ZJmMVYXNWfzuVggw2U0TCHGuHuJz1k31ptDZjJACox/gvG6fDnOvhDmTfsARySeWxTyETkBwtLjX8n7Y78U3L6QAhq8lHWQL2ZAmPB4WhvtRGeTPLevnGH7XbXXX0wZiiP0wKjEeIQRnsVQuHsuSwi0frr2qznALbYbi8uWLO2pcRYrPSoJSTe5S0+u17EBuYpqJg3wfgPrEWqsPuIUQsT4+bbflPne48psex0bkOtMovy0RGbMj4AHifJWiAXE82y/rwVnkXtXIdr9LsZPprOMHiQtQpTu7kxc0x8BJ+WU6Eg6gdALrAKOWc6oKljmVuMsppmUzl2I0P+Zq9cK7O8GmnA1rnGlWU2ey9o5VWb0okQpZLUAXweYZXvjpm2rcU8nvvusTK5a37EaEZw6M7X17cQY874MzpL9roCXAkvSyUi80najIFwaS/vCDcu/NNYJdgVgjBOsByh77sF5LJGFcd2kfzwxz43A2tX8QJ3yfQGXOxle/Ku5unlLPFvpvG2BrwMXE8/8vwOftn1yZtt+RJQGlpDFKDJm9Zm3Vciav0k6jPicm9h+ZQrSnmd77aZciW8nQvPyVOI6bAV8pSvb/GmBCkoLFO7fdgdmdgXKt7WdFQxNmyaVOdHyxGb0L5xR/SHpZsIwbDY1I0NnaFEnvutsr62xki5Z31uNs9X4V+MpeX+cQHzvt9EpMXZO/5H4isk4SbqSmI/XDQi/a3v9zLatRFQJXW37x5KWBba2fUAG19XAPravSK83BA7KbdszHaMMwhFGeBYjLVr+IGl74J5qkasoi3wpoZHTlLO0Ky9EZsIDRJ+1ksIdbmC9niFmd/xT0kZdA06uG281iakHfEzsXA4DvYIw/VB3NnyC0AI5vcV7l3S6XpzI6prF2MBqzs7sHsC+RJAWQlMvR7+texFalZUspSjjyAn6nklyViXT3bOONOH9AfCDNHk9EfiWpNOIgNfAmWe2d0gT6W2BYySZ2JX+cZOJtYfj7ltEPiEtZq5M/35o+3MF2jZV0gJOGq2p381xbyydGY1CCP1o0rMv6SHg/RmLt6rv+E36V+HMHn87KO6XtAMdrcVtycxGTJl15zhkFLJLDJsGAAfAMJyH1wW2l/QHoq8U8RU0Hf+mdG0efgrmfpdNy/BKctVR6tkiLUwvIYLvAJ9ypsFOQsnstSJj1pDmbevaXiMFgauS4PlzyWwfp9CArdr6TtuNyu0rSNqMyAyrjB6qZ2F6ZttKSguUnFfuYrsKHlfXYBfy77fLgNdUwV7gOiLDf/sMridsH5bZjl74m0LGpZLY2IrMkt6C4x9Q/P5Y2/YKk//ZwCjmIA/sDZwqaYwBYW7DbN8u6RPA8pJWJvTKGwcHE3YFjlOU8IvICH9fbtue6RhlEI4wwghIuh7YwPZj6fX8RBbFwDu9KuzGW+M9gBhgunfLmriGFs3uqPGuRoiNjxlwbN+cwTXV5TThBnm/rGyKQu/9emLy1drpuuTObCnUsjGmEQLLNxOfcxVC3L7xjqWkX9leuWAbpwJvIzIIlyHK8GcCrwG+anv5/mf35Xw+UQ61N+EW/grgENuHDnh+MUMWdeQTNiIW4BUWIXRtXj8oV+JbmSjBqf4tTLgrX0k4U1/bhC9xfooom6nMSnYm3Bcndc2dgHNJOoGMWbbvy+S5Bdjd9uXp9UbA91tsphRD6s8PJbJWTZS/7el8V/UimXWJa73UtlcSWnpTCdH2RoEHDcEJtt84mJHdeAewTnfwPy28rnUDI6ySXF3nt362JK1o+9c9NnyA4WX3P51QD31gIEsfWNK1RF95XQoULkFkEObqBla8LyTGVgBynntJdxMlu7e6wEI4jVkHEHqmokXAUeEIexydUvYHCBObxqW8aXN8leozpnH/FmcazFXzRkl7ENqhB7bIbvwCIS9wBmNLZXNdrl8OHEnccw8AvwO2z5nflxr/JF1heyNJcxibpNDm/jga+EZucLwHX2kH+SIGhIlrYwoZOtU4pwO4kKHNMxWjDMIRRhgB4DlVcBDA9mMZO70Lp5+LTPhXzfEOQjA3O2tqCNkdFe9NwKqFBpySmnBFocJaMbYvVBJCToeyna5tX5oWvsvZvkBRAjo1hystYD4JvIqxC5BGn7PKxpD0E2AN27em1ysDX8hpG3CVpFdXXAVwF1Ey9w3bdQ290yT9e59zekLS24md2FcQC5t1bN+XrsXtRMBkEJQ0ZCmqmWT7V4Qu1JEACgfG9xDB0IPIuOdsH5AWIlWwcn/bv2zKU0HS1oQY/SXEZPpQSfvYPi2D7slqcZTaeoWkHMmJogZRaYH71abnTYJSmXUA3yXui1OJzYGdgMbBdiKjdzMiY3icEyzQuNy+Gge7AywZ+AFwsqRdq+BM6oMPo7kBQkmuuSj0bH0M+BBj+4+5b0HD7H4NoZKh1JhVQxF94IRDiODPCyV9hSgJztWARNIWxLV4CRFYWprYiMoJdv0R+FWJ4GDCgcDmbeduqX/b0XapeeW5xPN1RHr94XSsRRO1PpEx+IF0LGu+RUfTr25AltW3pe9tN9tvUJRBT+nedGiIIuOfk0SC7ZLrovWAmxSSLI/SYrxSGLEsTgTLKymcvdxOxmkFItg4DVhDUfF1XCZXMUMndcn0KBnNOVOm55mOUQbhCCOMgKTzgUOdnLHSgn/Pplk2w4CkXwDvdk34NoOje/du7q/I2MXrlyVZISdbUgU14QZ8vzmDTlpUWCsmcW7A+IBj40lEKpf5EKG7NCMFHg/PuXclnUeUan2CKEd4L+EYlyuUf1v3Tn2vYwNy3U4E4FpPChPf3BKm2rENbV+ZwXUscFSvXV1Jr7d9YU4bSyBlFXTLJyxp+/cNeaYCqxPZCRsCMwgx76sJfZynLWO1QnpON62yBlPw4AI30BKqZUvtRBhF/JiOUcz/2Z6w7+vBt6bt2SUzfSVdQeiaPTbpHw/GVySzLnFVek63VM+meujhPR3oF2DJ7I92JfTgqo3Bh4GvO6NksCTXvIzafVaVAFfaeTsQffmnx581KWfpMauYPnA6d0UiSCvgwjYBtNS/bUL0aatLeh2wg+0PTHJqL661iRLjSxmbvZZb6XKl7Q0n/8uBuK6xvV4hrilEULCaD51PSGRkVaukfvzjRIXRAWl83dv2npOcOnSU+N5Kj3813vWA26qgZZrvr+S8yoNi41Xiu972Wjnn9uDaD9iYCBCeA7wFuML2Vpl8c8fRiY4NyHUuHZme+jqm1wbQsx6jAOEII4xQaU7MJBYOInZXd3IDHbIa18sJfa/1iIH1auCjtn/bkOfQdP6/Aasy3rHraZuQpEEQYqdsbaBygN6cKOvboSV/pQm3GNBYE27A9xh40Spptu3GO3YT8B1PBFhuojNQO+eaKlzI1iHK0SqB5VsdmmJNuWbbXrNrcX+d80XVf0xkJNXF6J9re9sMrtKTwnEl5r2ODch1QPeCtNexBnzFDFlUQD4hnfcPIhvye8Altn/XtC01ruJlR4l3zH2fFoc3N3kWNASzgsRbzCxG0nFECe9ZjNUdzXYyV5T1vSa9vNwZMhGJ5zLgDUT221+ILNb3tQiwnEUsUs+0/Y8cjhpXsQBLjXMRmKtp2v2799o+9l/JNcRnq8iGVuIqaaBSesy6gdiQresDn5bTtnT+VMIAqP695UoBVMH3m4HVHUYZuWYK5xGB6FvpSNdg+4sNeSo5jNcSems/paWzrMLc5d+ILOR6/5blUpvGvBWI56FVyWeN87mpTW027nfqdbzFc9X6exvi+HcjUU1SlXpPIeRmsmV+VKDUPvF8Hfgb4w0Sc2QsbiXWazc6smCXBE6wvWlm27oNnXYgskNzDJ2KyvQ80zEqMR5hhBFIk8H1Sgz6RGDre8CW6fV7iAXOug15KsH42XQCcFlQYT2nahKZFoNr1HYFv0AYeOS0sVsT7mA6mnDnkFemNhEm3R2qfV8/k7QbhbRiiNK7larJUks86iiJB0DScxjgs/VBNXG+V9LbgHuAnvfMgNiZ0HPaK72+jCjXagzbfygRxFCUB20ALKGxmbDTyS8V2pRkLlDDW3ocGxQlDVlKyCdAlFOtT5Ta7SzpOjrZg39uQuThlB0BnCvpl3TMO7YBftGwbcMwK4BCZjEJleHJFApIWkjaC9iFTgn7CZKO9IDamV3YMbXrP4CPEppJPTU1B8TBxHX8errnTgLOdoZrNvC47fslTZE0xfbFkr7dom09g3k17EXoR/3LuIbxbPXb0CLkFDIpO9naKfg4JZOr9Ji1D3CxpDH6wDlECp26/YD/Ib43Ed9bro7pg2mOejkwU9J91AIaDfGSQsGCuhzGPwh3+wpNZTEqTCNMl+rBqCwu9dBwS8H2LA03Sa8m7vvnxUv9lUgmuC2Drh7EnkZkOd5A/nPV+nsb4vin+nw3BbezYjAqW2oPHRORurlRVqk3HYfrJxQl8vcRY2Auug2dLiPfYKe0TM8zGqMMwhFGGAGANLns1rHJydjplRKeXaJSApLOtr2ZQrNjnJ6T7ZyBEEl3EgLQlVPiAoQAdGOHsTQhv5go1byq63eHlM6YHCRjoc/3VaHN93YqUcKe5TDXxXUg8CBRErIHsBtwuzNcZhXOhpcTE5pDiaDZF51K759O9AhibAk0DmKk7NSNiXK0w2u/mgP8zPZdDbg+QnzfM4B6husiRJZeViZtyZ1eDUE+QaGtuA4RaN0ZmN/2hEZIE3AtTtxv9QybbPODlNGyUXp5ue0zGp5f1GxKhc1iEmfRSb5Cq25924+k1wsTgd+cMqa32P5F17FdbR/e75wBeacSC99dgDfnZMJJuoDQ9P0aoc15H+GIuUGbtk3wfsVKqwfl6rcJWCEzK+YOym1oIWlNQkKkMjd7gHBIbfzcD2PMSvOY1vrACiOQdW1nOYz34FsY+D/iO9ue+P5m5vCnecMFts8r0bZ5FQrX5+3cpeHmzIoQSVcBn3MylEsByK+W6EMkLUbI6ry5LVeLNgzLbPEnhDZwtUG8G/A62+/I4CqeCV4Kkr5PyEW8hyhFfxi4yXbWJkMP/jcC++RkJKqwTM8zHaMMwhFGGAFJhwMLAa8jSqO2AmZl0v1C0qeJTIdKu+OcauLedIIuaUPC2GFpos+qOvWBg1O2N0s/l23y3gPgOGCWpGoR/g4aZEx0YSf30YQrHRys6Cf7g9LflzpGBYsAt0uaxdiMxBzjgU8T2V23Elo755Apbm/77PTfh4hnIQsaghg98RnXrQUxDiAy2BoFCB2ab5dKOsbtzXtOJDLUvkZchwpzWmSXQtmd3l2JbJPvQkc+IYcoLVDXpaNDuHbia6zbmPj2J7LofkvNoZ2G5gc1vqqs+yc9jg2K0mZTRc1iEr6fghjHEAGCh1q1MO6LuiZXle2Ug30lPWr7IgBJ+xDXMztAqNDN3JwYS9cgf4x5OxFg+SidAMswBdpLZiAMyvU34E9AZSjQ2tyFMCd6EXEft4ZDu3dVhVMzbe7fgmPWDkTSyPEpIHhLOr6jpCdtn5hB+8fUriKw/YikFxGbM/8L/LJF8PEjwCckPQY8RvsS9GMJc4cH0+vFgYOdVwrZS6Znb+fJWsxXBQcBbP+XwmU2FwtXwcHEd0kaF0vgESB7zlnoexuW2eKuhGnP51PbLiS0s3NQJBNc0rqE6doMYv78frc02bG9W/rv4QrNv+nOc9+uxsyXEGX7BxCO9AK+ktm8t2Se96zEKINwhBFGmJv1V/v5XOAXtl8z6cnjuSYajBtnnUn6NbGg6RaWHXhimLIvFnQqnVYIBldlhjdOUto0GfcadMo+L7N9YyZPMU24dG6vHdCHgNm2b5L0Xdv/MQnHhKVxbqiJoz5GBTW+p9XooTbBXJ8I2OTqZ06YTZYTmEtBx7XdMduYBlznDK3FdP7F9HaVbRyckrRUr+PO18QpvtOrlvIJCg2hlxHSB1elf9fk8iXOO4FXu5zZRq8+JEvQe15HyoTZGXg3sZl1TG42UOor30tIKEBs9BxjO2fR9QLgbKJU883AisC2uddY0ilEQORcQiPqUttPTXzWvIE241cuV1oov44I2v+YEMhvtdBJfeVqxH2WvaE1jOykgmPWtcDru/uzFPy5LCfjTNJRRCbizyljBPJB4D8Jp3sRun9fsv2jHL6S6JXhmptBK+kaQqankop4D7CH7aYyPWi8htv2wNScwGXiO4MoA64b7Kxpe8v+Z/XlqjvbTyGMLU5xhllP4ivyvaX1wp62v5XTjmGjVCa4Qp/5M0TZ7hbAB22/qUD7qioGE/1voyqGxHEjse67mgjsnQB82vZ3W7atiNbwswGjAOEII4yApGttr5sG2HcSOh632X7F09y0uW1ryXEQcJ/tA9Pr3xFZAdOAGxpm2NR5WwdG1NGE2xuoT0imA1s6X9z+RELr72fp0GZEVsAywKnVdzEJx9ET/Nq5k8zEXWUCmAh0/SWTZzliorQSY8vjG2eKlJyY1ziXpKO1M8vJYTaDp1gQI/HVF3zTCFOQJ2x/MoOryphU4lqWKE3L0sTpF2DNzXhUAfkESasAt04WbFADYwZJpwMfyb0najwTlXpfZXv7BlyHTPT73GzmtGg4AHghcZ+0ythJnFOJ5+AQ4O+J87NNNy4S1xqMLc3O2uhJXC8ELiA2td7fJkAl6U1EKVmW62ji6DbsmPsrWl6Dyd7XhbQAm3BJEiGjsC0xxpwHHJaZgdV3Y6vphpakD9s+Qh2Ts26+RgYZibNUUKRvADZ3k6Hk50x8dxKGU/en188n+rccSZeqTHlZ2/tLehnwYttZlTOKss+NbT+QXj+PCObnmKUVk+lJmda7U+vbgO87v2x8ceCLXXxfqD53Q676c/UE8Afbf8ppV+Ir+b3Nsr1Oblt68JXMMC1Sat/9zJfY0FGUGL+CsTrIv7G9e/+zBmrbnTnPeRdnEZmeZwtGAcIRRhgBSfsSpYqvJyabBn5g+z8zuCqzjWUYq+wc9JEAACAASURBVKuVu2v8dcI84SeM3YUeWK8n7UatbfuJ6rVDu0PEYnCjiRn68tZLSRckIzCigppwXbyXAW91J2vyucRO/puJLMKVcnj7vNfAQZH098UyASRdQQihf4sow9uZcDnLuXeL6mdK2hr4BqE9I2Lnch/bp2XyFQti9OEvMilO7dzN9gdbcGwELGf7aElLEO7PjRf46iOf4CHp9TSZZEtaizBk+RXtMpMWBRanQKm3pPdO9Psmz3kX793A5m5ZwpS4ViGe87cB5xO6rTdIegmhHTiQHqQK6tXVgnCVEcP8xKLXZAThJG1i+yL1yeLOCYL+q5GbQVWKS6Fp9h5gfyJw/IMW719ko6c0So1ZCp3FtZwkLGrHFyE28FZs39p2UOjfbeyxjvSXNM2cSuceRmTWbWL7lSlgc57z3Z93IrTXTiX6gK2Ar9g+fsITe3MdQOhS1mV6FifmEm3M4Z7RKPm9SfoWMB/jnX2z9IFLZpiWgkL3/BO1QwfVX2dutP0aeGW1KaZwa77N9itbtu0bRFZ+m7YV0xp+NmAUIBxhhBHGIO04TnOmLo6kc4jdrVvp6Gq12TW+uMdhu0EpZPeEWdIbncrRJN1ke7WctvV4n+zAiKSlczOk+vD9mihffDy9XgC42faKpScmTXceC2cCzLa9pqRbq9366lgDjipQ8Cl6TDBtf6ZpuxLvzcCm1WIyBbouyNzRPoQQ8L5q0j8ejK8eHJkCrAkc0naXtsY/93pknLsfkf26gu3lU+DnVNsbZnAVk08Y8P0GfrYk3QYcwfi+MqvUXiGdcJs7rurTicn6tTl8iaOEsz2Srsy5fn24LgWOIu6Jf3b9bsdBF+WSnmICvTpnmjCVgKQv2t5PvbO4nZl5UlQKYID3ezpKjBcmtBa3AZYgNhZPafMZh7DR81JiQ7Z6Hi4nsosGzp4qPWZJ+gSxQbxrNQ+RtAyxYXyJ7W804UvnLwF8kvHZ241kLNQpyV4NeDWxqWLiOt9i+30ZbbvB9hr1/rrNZmA6/1V0dCAvsn17Jk8xmR6N1/CuSJpK/dTLgXs1auBNLQ0po7nw99Z63dHF1zrDtPT31mdsqZA7xpwN7F7rQ5YGvmt784nP/Je0rahMzzMdI5OSEUYYoeood6OmGyHpsKojbYiXltyRsZ0tvF3D/JIWqRbOteDgotQmrm2Rslhyy1GPkVREEy5hJnCtpDPT682BE9MCKmviOgGaCvrfT2RIVpiTjuXg0bRLeZek/wD+DDy3Icdsxro1f7j2OxM6LTmY0pVpcj8RjMvBbODzklYgyoxPsn19JlfFV33mJwi9v6ysOo3V1JpCGCnc06JtWwKrE1pH2L4nZbLkoAog/SMFGu8HXtyibZOhya7rP2xPWNLbEIcR332Fh3scGwiSViZ0pp4XL/VXwkjptsy2XS/pZEJwvJ4t2TgTwHZfLdOGGTuHUEivLm0O9UXTzBPbVXnml9yVOSspV8j/57X/z5UCIAI4zxTcB9xFBMzuIp7HtVK2bm7m5eeIheWYjR4gK0BIiO2fSOhnQui4HQ00ceYsOmbZPkjSw8Bl1aYA0X983fZhE5w6EWYSGVibERUS7wX+msFT9f2/Sf8qnNnjbwfF44pqlyrTaQlqmzQ5sH1b6ienJc6lcgLTLmsOdxQ9NLwzcFCZ5oBrUgElN6tLfm+F1h11HAxcLWlMhmnDNhU1TvGAzsIaoEJIYw0I71AYEJowdWtctl+ybTUcTayJ6jI9RzVt27MFowzCEUYYAYUQ+hw6QsbbAYvZfnf/s/pyHQBc6EzB+B58SwJfBV5i+y2SViLSxAfu2FMA4w3E7vh/p2NLE4vni2xnTX76BEae7wyhXxXUhKtxrkUnS+HKlgGlid6naQbhcfTIBEj/GpWjS1obuANYjCglWxQ40PY1A3+AIUHSN4BVGKvHcoszNS8T5/OIe+M9wFK2l2vd0JbQWK2pJ4DfA6dnbjDMLXWuZXlkl4KooHzCgO/XpMT4m0Sw7Cwy5RO6+MZlQytfP+wq4HNOjpWSNga+6oxyvnR+yUy41s72Na4ienV9Mk4qtMk86WU80yhDegLu1lIAk/D/y0uMJR1D/yB9duZJPcskbUjdnJt50uc5LVbJ0BbVZox7mLc1WZCrk90/tw+SdJ0zy3hLQtL2jHUF3wr4vO1TM/m2IIJALyGC1EsDdzhDh1fhMvwR4N/ToUuAI5yqQRpytdbw7uJbGPink1FSCrIuYPsfmXzFsowT38qM16M+LoOn5/zADbWLuzhXIhztISPDVAUlMRq+76TXSE+TAWHG2mOoMj3PJIwyCEcYYQSAlT1Wk+5ihZNoDq4BzkiT6MdpWTYAHEPs/Hwuvf4vYld64ACh7W9K+geRGblwOtx2dxw6O9sQgZGfA6fnENme3XXoyrQLl40UEBxKULAl+mUCNN4htX1d+u/DhC5ZY2hIel+295H0LjpB2iOd4ejWhVcQzqhLE4HRxkjB8Uds/01RlroRcLftn+bwOckHKEpa3Wth2RCnSDoCWEzSLsD7gSz9MNv7p/+enspfsuUTBkSTgG0V8Fivdsx0FhFN8VtJexIbHxBZ4Y3cTGtYuAoOAti+pNZ3NsagGQEDolRWDClj8GKFTm2lV3cXDe+3iTJO0oK/ESStSGT2LdrVL02nUNZ7y4z37g2yCg8RGrc3EZmZ/1IuD1hu2jDz5FxJv2TsRs85A57bC/dL2qHGty2ZGfQlg0kVJum/9yICaoOgasO9CqOoe4iM5EaQ9G3be6tPmasbaramc2ZKmk1sHAl4h9vpo+5P9OMXOPStX0dkhubgMEL/7vvp9Y7p2MCB/FpG88VpozJbw7sLFxKb7ZXkxILEpkrWxlFJpI3KjYkA4TmE++0VQOMAITXdQaK/3YzM+VZq21LEd3ZW/VjDDNPurOE6DAxLEmPSCqFhBQAHwKRt6wqs/j79m/u7YQVW/3/HKINwhBFGQNIJhE7ENen1uoSOxE4ZXL8jMsImdfwckO8622trrFZM9m57qd3xrvMWyt1BrXEMVRNumFCmW6UK6JtJWp4QL+7W2GmiUVlc76s0JB1IlN7+hgiQn+HkiNeQZ1/gfcSE8iRisn8JUQpys+29MzjXIoL41T3wEPCBNhmrkjYF3khMAH9p+/wMjqKB0AHf82kTHle45x5CBBhNLOb2doahQirDuYEoM4ZY7K5pe8vMtk0jSti79chyMrqKZMVoCHp1NW4R12E7YDPbSzY8/+1ECdQW1BaVRKZ/lhZpyYz3xHcioRX6s3RoMyILfBlCH/LAp4NrwPdrmnlS3+i5vM1GT+qXDgXWJ57Tq4A9c+47ST8kgknVvGVH4Ml5IStU0maEvuLL/l97Zx5uSVXd7ffXgNDQNqCiOITB/hAEZFYEiYqiGBRQQVGGSNBEIgaicYxxQBQnEAQMk8osQ4NGJEAQgW5B0EA3iAbUT1AcQYjMM6z8sfe5t+7pO9Su2nWqzj3rfZ773FN1Tu3adU4Ne6/1W2sRjnc+cLCZnT/thsu2s6WZXTeVSqmOcSLeM4v3o0rXvqRrzWwrhVxzm5vZk6peQXeZ7VLbalDRXFv92ufwmFAYI3aukkNWIcfcpsBSM9tUIfrodDNLCd2fqu0VCeOQV9boW62Chm1RUkF4pZltp2XzJNYViOToW2dzDXcZVxA6zghTeGitAPxQ0m1xeW3g5orN/hb4aQ7jYOQBhSIWvVwxLyUYICqR0TuOpG0IapZ5wFqSNgXebWbvqdC1bDnhWuD/p3xYE/ObIelOquc3W0io/nwiFdVEFvN95VI5TTJIGnuL6oOlXxFC6++s1bmgVnkhobLvbcCaZvagpOWB6yu2+Q1CqOIPABQqEH+DEF5diWgQTDYK9igaQiUVDaGvl/TKKobQkiTd96Kypt9oVimMKRoC31Zl20nYDziYYDSDMNGvYyg/jfBM2RH4NLAXiYqMBlQx2fPVxefTngTj3tOAA+ibAJfBzL4DfEfSNmZ2der2U5BN8R55HrBFz8ETFTz/SVCzXQekGPVytlWGpLy5ZnYe9b6rsNMQknloFdXbFLy4z3B0WTRSNUWp+1s8zvXM7ALCeK1yTjeL0RU5VUqaIiSY6vk4744Oz8XAGZLuYKIKLYUnJC0ws1/Fvj6fxLHNdIrmmjwgaYvevVYhNc5DM2zTT7FgxaK+ZWP8mZPKQ9Ew+3iMZriDYJzOwcqEe1QlrC8dQXyWVZkn9JxPewHrmtkhUZ24ppnVijiqg5ltF/9nzZOYiWy5hkcJNxA6zmjzhgbavAW4QtJFTJy8lc4r18f7CQqKBZKuIig9dq/dy8lJLbZxJGHCez6Amd0g6eXTbzI5ljcx9aBJfdieALzfJuY3O5FqYSqPW70w8TFKhLmVwhpIwm1mx0taXdJLmGhMWpzY1MNm9ijwqKRf9ZSvZva4pEcrdu+JnnEwtnWlpMen22AypjGs9tpNMaw2YQjNiqTjCP3bHvga4b5WeZCfU6VnodrigbHd1YG7aw6q/5+ZvUXSrmZ2SlSM/WDGrSZyeN/yVoXXVUKzF8bt1o9/RZImqpIOJRSduI0wCTkYuLaKIr2PpZIOIM9venDNvvTzTArPeEJI6bPM7CFJj0yxzSDaKkPpczkqnr5A6KOo4egxsyckrS3pKfE+XJfaxqQmiMf5duCIXG0qY+5R8oYEQ1AiP0RIfbAXIRdy1evtgwQnyC2EY1ybxPQpClFAJwALgBuB/axeCHWPfwYWSvpD7NuaBAV2aco6YitE9FwraTXCWPI6QkhvJedKn+JvOcK8o3L+wX6sXnqHfycU1HkV4Ty+j+C8aCq354xpU9RSfkRK9M1CeoJeruF9gKMlVco1PEq4gdBxRhgbL0W/VsZmb41/T4l/tYgP0lcQJnAiyPIr59eZaXfJG5j9Njx7xkgenLcRCtkyOfObfVfSewiVfYsG6SoDkq2YPMxtf0lVw9yyeColvYugcH0ewcD1UsLgN9Uoslqc8AqYXwj3EWFSk9KnnqJrkULOwDMJx7sHQa2XRM+wKukQ4I8E1VnPW55aebgJQ2gZUpwM25rZJgpJ/A+WdDhwUY1951DpfYIQZntzDKu6iBC29YSkPc3s0op9692z744K4j8RDC6lya2Ksbz56t5FyI97LPBdM3tEk1Slr0Dt37SHJs/hdg8hT+3xll5U6AxCVcheDtmdgW/Ge3lqDuOcbeXmi8DOmQwsEJyoV0k6n4LKrKITdTJjUpMpMVJyrF4l6RhCSozicVbNf5ct9yjwmJndJWmOpDlmdrmkI6s2Zma943uSGIUSHdovm3Kjqdv6vqT1GHda/NzMUo3kXyUolxcT0hT0nNm1MLP/VsiPWuxbU+Px0hE90QD0OQtpV46TdDEw38x+UnHfRfHE48DtZpbs9Cz0b7L0Dn+o2NzWFoq3LYXgzJNUe641DWUihNrKj1gqeik6N2vnGh4l3EDoOA6EUJ7ezX0lYn4MKoRb5FYpaNnCES+QdA8hx2Fybq3M/FbStoQwxhUIA5rUCXlboZA5SVVe3hKPu5jfrGoxhXfE/x8srKs6IBl0mFsKBxE8xNeY2fZxkH5ohXaKIT2LmRjek6pG7Fd0FasZ1zGO7NIXNndsDJtLqTyczRBapITKtHRhBsZDsx6U9BxCoYJUQ2iRHCq9PQiDZwjX1hyCIe8FhMlaVQPhCVGJ+HGC4noeab/nGAopJz5JcKQYIRH9p82sUqGHEpSZqD4beA1BuXqkQg6wuZKWrzOxJM9v2uMWghKmWGzjPsJveyJBXVGaGN52EeNGkP1tPO/oXm21VZIUQ9ftGY2DMF6kaw4VCnP1cSXhWMYMNjXbm4mUdCK9vHRF5VWdIkz3mFkdB0qRnCHBU5HkeFcoXCMzOy0aBH8S1+8j6Qkz+2ZCc3NsPG/vQkkfTenLNH1cpiiOpFpFcabbXdkPmplJuhB4UVz+daUdjivh+lMRzZdURwmXM73DYwoh/L20S2sQDNNNMeNYrsUIqBn7pslzDW9pGXINz2bcQOg4Tu78GGsAH2LZkKiqg8J3EpJ59xRnryQYa9aV9GkzO22qDSuQMmkA2B/4CvBc4PeEam4HJLYxDKGQOY0iMDG/mVEjv9lkA5Ma3tQsYW59Ru3V+parJuF+2MweloSkFaPCK7mATc7wntyKrgIPSNqLkBvOCNdI6uQtpyG0yEwq0/cmtHVBDIn6EqEgiBFCjatSW6UHPFoIJd4RONPMngBuivekSphZ77gWUV9NcBbhN9wtLu9FUCntULPdqShTxfEJ4GJCxdsVCefFXOD3kr5vZntW3HeO37THtmZWDEP7rsaLgFXJ/0o04lUuRtRUWyWY0dBVuG9fK+ls4D+YqFKvWt2+V/W9dnEz4GoLSfrHlFKSlhAUSk2QkoPwfDPLFmJM3oq8uwIPMzEkOFsIaSTVSfZPhKrK/XyLcL9LMRD2jzsmLFc9d8lQYTmB1O9viaQXm9l/19jnnUxT0IKKz66icEL1U3YcRYiYeaakzxJSk3y8YlvZice3HvVS4eQke67hUcANhI7jLEPN/BhnECZrbyAY0N4B/LlGd5YHXmhmtwMoVCY7lVB1dTHjKrQcJBXbsFAwoq66oa1QyBSyGUXixOFbuY1LMcRkrGookFQ1NFIMc1Nsp0qYWxNJuH8XjUn/AXxP0l+A31Ropywp4T2rEhRdPVXBIoKiq2oxoT0JhvevxOUr47rSNJjnKKfK9ItRKXKepAsIA+rUMM8iOVR6j0RD1O2E3IjFAhsrp3ZoCufCGBXDKp9tZocUlj8jKSkPViJJE7neb0r4XecTCpYAlc63bMpLYJ6ktXrKCYXUIvPie1151gyKMr9p8b79IKGqenH7qtVWaxc3k7QmwTE5V9LmjBsy5lPhOs2NNZCDkDDmg/q5R4shwZBQmK6fSSJcxt4iOAhSWKH3XCliZg9E5V4K/eOO4nKdQiCDLoqTwtbAXpJ+Q3Ao9nJUphRLy1rQQpOn7NgMeFwVU3aY2RmSriMYkwW8MbPCuZ/SSk7lS4VTepclPpMt1/Ao4QZCx3Fy58d4upl9XdJBFqrOLZJUx6P3Vz3jYOSOuO5/JeUOa0itQPpF4DOEUMGLCVVb32dmpyc000goZGayGUXixOFJSavWMCCNoUxVQ2PfsoS5NWGcMrM3xZefiuGLqxLOuaZICRv/BvBT4K1xeR/gJGCqydO0xPCgXatsW4GkyuXkLaZwNVHpE41Kj9RR/2RS6R0EnEsIxTnCYhJvSTsBSyu010RVw0skvQ04Jy7vDvxXA/upjZndS3Bo9Ug63zIrL/8FuFLSrwjX97rAe6IDpG4xlVlH7z4u6WVmNkElr1Awoyo5ipvtSEhN8jygaGS/F/jXGn2biZTnQtYchDmcipq6EFbVwjM7T/PeBYltzZW0Sp/xEklPJTGfd4MOskEWxUmN6MmRYzF3QYvJUnasQY2UHZJOM7N9CLlp+9clkzlCKFcqnGx9s7y5hkcGNxA6jgOZ82PE/3+U9HqCoXHaClczcEVU1yyMy7vFdasAd9doNwevNbMPSXoT8GuCQWQxkGIgbCoUMie5K0zeD9wo6XtMnDgcWLYBNVc19DFCPhdj/FxuihmNBZq8OtyN8f88oKnqcCnG8gVmtlth+WBJlcPjJT0POJpxQ+0PgIPM7HdV25xud4mfr11MIbf6J6dKz8x+BGwwyfoLgQsL+yw1mLb8lXMB/p5QTbOnHl+OEJb+bipWl52B1InqdCSdbyUmSKUxswsVih/0ft+f23hhksoFGoaUlN/haJY12k+2rjRWs7hZvPZOkbSbmVUdqy1DZmNB1hyEMXrkUOA5ZvY3kjYEtjGzr5dtw2IhrNjeUjPbvEpfCu3lNMR9HThX0v42XkBwHULBkdLHmEiqg6x2heUEUiN6et/ZMymEt6YSFYO5Clo0kbJjQm74GJWzZcW2IG/alCypcBrq20ykXguzGjcQOo4z6SQu9eElaQULiYo/E0MO/4UwiJ5PyPNSlQMIRsGeseBU4Lz40M2dAy3VWND7jl4PLDSze/oG/TPSoKc3J7krTH6LZWX9qWEc2auGSjqIYHw4j3AunC7pBDM7uk670+2yxGcmqw7XW26yOlwKD0nazsyuhDF1zUMzbDMdJxHyLb0lLu8d172mVi8nJzV8NIfKdCr1z31UU/80odKbiaTBdFSafIUQcmQE9eT7zCypOFFUd2xkg00wnjRRnYHUe1S2yuqSVgbeD6xtZn8vaT1J65tZqtKp8+QwdMVQ4G2BNfram08wSleldnGzAldJ+jo1jGZ9ZJuQ51D89XEy4Tnwsbj8C4I6seqx5qgyXpYZ75dmdpik+4HFCgVUIDhTP29mxzbUr6QBq+WpsFx6dykflrQLoXDacwiRRmsTrqvSxRaVv6BFtpQdCoVm/pXgWLyX8d/uUeCEiv2DvGlTcqfCGWThwNT536xGNULrHccZciRdaWbbxdcTJOqSllhIfl22rTsIITNnApfVydvRFDNNGiQdkzIAlvR5QmjrQ8BLgNWAC8ysav7G6faV9Hs0sP+tGDeKXFUwilRp6yAz+8pM62ZoYznGq4a+mlDEZgdC+HmlqqGSfkKYXD0Ql1chJIFPyWGTsr9Wf9PpkHRfUW0xw2c3I0x+ViUMsv4X2NfMKuUmknS9mW0207octPkb5Fb/DJJU9Y2kawhKmF4F3bcB/1TlXinpRusrrNUkOc+RCs/VxcBOhQnSPMIE6XWE59aGCW2dTZhU/a2ZbRwNhj9s4rpqG4Vqz5MZutYhOPNmnFhKegUh1HB/4LjCW/cRnFK/rNi3ZxCM5TsQ7peXAAdahQqp0VlxEvAxM9s0OnaXVr0+Mp9vtRV/fe31CuqM3XvqPBcGee+vcL98KoCZ9VfTzeosrvIdROP2OhRERmZ26pQbVKTCvfIGgjr1UjPbXNL2wN5m9s6ENh5g2YIWY1hiQQuFXO6nEIyNR1rMnauQsmMfM3t7Sntx28+ZWZaq1LG9m4EXRYEHCrkSbzCzDeqobOP9c1XgYgt51jvTtyn21dnxeBu4gtBxRptVCq837nsv1ZvyQkIuqH8jhL6cR5DTX1Ojf1PljLmHUO3wXxIVKFnl6mb2EYU8hPdYyK33IM3lTmvVu2V5K0y+g/ECFD32nWTddP1pomqomBjq9QQtf++S3mtmx8TXG5lZpWqjFSitmoqKnE0VCjL08q7V4S5JezNuTHo7cFfNNqciZ/hoKhtLWkbdYGZJ1TQlfcjMvijpaCZRXVhC6H4CqQ6glW1ixfnTJX2w4r5zVKtsi9TzLWd6hwVmtodC8QjM7EGlSt6Hh9rKExvPoXyyxfDFTKxvZhPUxlF1nRK+2+MZZnZOVBf1ipvVyQmX83w7mbyKvwckPZ1471HIP5yUx1jTVPaFRiuapirVlzEMFmgtFFLSacACQgGK3nlmTMyzmm13iZ9/zMzukjRH0hwzu1xSauqErAUtLHPKjrjtR6Naspe39IqaKvBsEULxmvyZmd1nZovimHBz4Edt981Jww2EjjPa2BSvJ1ueviGzu4DjgeMlPYcQHniEQj6Qs8zsY9M2MDVHAr8jhByKoDxZACwhFEd4ZUJbWeTq0St4QuzHjcB+wE1RefbAdNvWoHOKzFTixHRPYF1J5xfeeio1culZvqqhJxEGI9+Oy2+kudw/UM5YsB9wTHx9GjXyXiUy4/k2hSKXns3BqlWohXDMRxMqYBrwQzqS5ygzxYqVKxGM3FVCDXvb5DLgN8FFkj5CUGYYIYTrQsUcm4nqqRzVKlPIaURLPd9yTpAelTSXcQPLAiYag2YTOQ1dK0o6gWVVU1Urc+bMaVjbaNZHzvMtt/Hy/YQolQWSriKost4y/SbLMFVlXxieiqY570epDoutgA1zRAhNMX6omu8S4O6oeF0MnKEQ1ZQ0Hrf2ClqUNvpK+hwhYumM3raStjWzSsWJLFNxvsixTLyP3T/Jurb6NhNtOos7hxsIHWe0WU2hwMYcJnpTa1XQNbM/KOTF+QthUPcuxr3IqexiZpsWlk+IYSUflpT6QMw1afgqIZfIYmAXxqsSOtPzQ+CPwDMIuWJ63EdQctbGalQNNbMvS7oC2C6u+jszq1K5tSypxoKuqX16IcjrEyrX9Yy+OwM/rtpoVOvsUq9r5Xc3oP0su2Oz4jWApMOoUI3XzHqK6AfNbGHxPUmpE+iypA6mexWu3923/m2k59LMeq9tYKI6HW3kvOzxSYLq+q8knRHb3DexjWEhp6FrISHE+GvUqNiqZnIaTmY0271qHzOfb7mNlz8DXkF43gj4OWHsWhprL+dzTuNDzmdW6hjkp8CahHFcXXIXoNgFeJiQ83wvwhwmSY2fQG4VZ8rY7vXAZmb2JICkU4Cl1KhenjFCSEXjsZk9qerFWHpt5Ixemo42ncWdww2EjjPaLGJ8It7vTU2uoCtppdjG2wkD4YuBjwDfq9HHByW9FTg3Lu9OGARA+kAp16Rhjpn1jmlhz0PeMEPv3YqGn98A2wxwtzMOvDSxUvCv49/Ye1VyQ5WkzPlbNOLPH2BI1Izfm8XiRgp5q7bohUVJ+hRBmVttx9JJTB4qu1/VNoeElQkq56p8lPFq79Oty0Fqhcl1c+3YzH4jaVPgr+OqH1jFfJeRQVZKrMJKwL1mdpKkNSSta2a3pjQgaQ6wOvBmQqEYESqD35m/u+2T2dD1uOUpEvEUQuX55ZlYWOheKhr1zGyJQq6vMaNZL19XVTJOyHMo/opcHXOEjaXZkLSEZlT1uQ1AXTU+pI6hnwH8j6QfU3C2m1kVh15TET3vzGzcnXS3mdtL/R1WYzzqprKYowFukXQgQTUI8B4gqRBZiwx9lFZO3EDoOCNMWW9qGRSSgu9AMDSeAexpZg9Pv1Up9iLkpvt3wg38GmDvGCqVmjMw16ShP3fNhOWGjDZdHWAmE7+rLxAUnWI8RHB+A7sr89C/3ge6IwAAGeJJREFUkxDG3itu0l8xuM1KwUUj/mIyhURlVk09i1BJr8ejcV1Vivl0VgLeBPyhRnvT0ZoqU9KNjJ+fyxEm0MmKB0l/A+wEPFfSUYW35jN+Tuem1GBa0qvM7LJ+w/ZYIxXulRqvNt7btm618c5WSox92YpgADoJWAE4nfFnWCmikuNDZnYONYz3w0RdQ1fBcfRdSe8Bvs1Eo0iS42iynIaSVgfuTg3ZVMjRKjM7zUJRrp/F9ftIesLMvpnSXkPUVvwBSFoTeC6heuvmjF9D80msBJuy28zt5TQ+tOks/lTGtpqK6DmC5iN62jQmfQ5YKulywnn6coIQowvsDxxFyEVvwPeBf2i1R04l3EDoOCOMpDfMlNy2zGciFwPvniG5cnLohoUiJDtP8faVZdsptJfDO96vtiwuN5XHZjZ5t74I7GxmVfKtNcFRwPYEo9iZwJU5cuyUoIxKr6mQqJyqqVOBH/flbqzswbe+yr6SzqTCtR63HWT4aCpvKLx+HLjdqlXg/gPhnrYLwaDV4z5CuFWbvAK4jMnv4VXvle8EtrbxauNfAK4m5HGrQrZ8dQ2cb28iJHlfAmPpO0pVF5+ESyV9gFAoYiw3V4MK6WHnOsI52rtPF4vqJDuOJH0COMfMblYorHURsBnwuKQ9zezShOb+CXj1JOu/RTCUdMFAmEvxtyMhFP55QDGv7X3UCKucgS6Pt3I6i5MModHInYthjujJTWmjr5mdGVPhvDiu+rCZ/amRXiViZncQ0oYMI11L4dMqbiB0nNHmS5J+z/Q3xkOZqOiZFDMrW8UsKXSji+GGLeaxmS3cPkDj4IwDLzP7Z0kiFLzZBzha0iXAsamhfEUGbJxKDYnKppoys89GZW4v5DN37sb1CEacKnQufLSgTOp3psyXVEWZdANwQ1RxLw+sZWY/z9DV6Sg1mDazT8aXn+6/liRVDTvOXW08Z7663Ofbo2Zmknp53FZJ3L7IHvH/AYV1bSukO0vOsPjIHsAh8fU7CGq6NYAXEO7dKQbCFXr37iJm9oCkFep2tA65FX9xDHWKpN36nUdDRE7jQ2njZe4xSM7ojyGP6Mmt4ixt9JX0XYID4Pyek6wrSHoBIbz4WWa2saRNCHnkP9Ny17ruLO4cbiB0nNHmdiZ6ZCfjl5n3mTpQGmS4YW5y5rGZTd6tayWdDfwHE0O2WgvNjorByyUtJXhADyGc+yfW2PcgjVOp50dO1dRahDDtbxfXmdltiX3qbXsf48odA/4EfLhKWww2fLQs/cqkInUMNq8DDiPkOltX0mYEw1xyfqgGBtPnsaxy6Fxgy9S+kbnaeOZ8dbnPt3MkHU+Y9P49ocJ38j0p5iD8iJmdnbrtqFIM4+1bvw9QJYz30YIyfUfgTDN7ArhJ6Yn850papd9AENWlT0lsKzeNKP7M7DxJrwc2IowFe+ubKESRZADqsPEh9xgka/THEEf05E75k6JYPYzgbPi8pP8GzgIuyJTSqS4nEpTWxwOY2U+i47J1AyEddBZ3GTcQOs4IY2avLPO5zEq41CqO2cINW6C00abDA8wmmA88CLy2sK610OyoytmVMOhaI/Zjy6oGrgKDNE61VbAHwjH19j8XWJeQb2qjxHYAMLOqIZSTkc0QmosGlEk9PgW8BLgi7uf6Giq9LINpSRsQzoNV+1Qe8ylM8lOwBqqNZ5qoQubzzcwOk/QaQiGL9YFPFMLpUtp5UtIHCeHFTjlyh/E+ImljgmN2e0LetB6pyrqvA+dK2r+Qz3AdQj62ysbyHDSl+JN0HOF72p5QUXp34Me52u8j1QDUVWdg7jHIIKM/StFSRE9rIeg2nst0OeBVhHy83yA8U9tmZTP7cQjIGaOpPMipdNFZ3FncQOg4ThlyV3SrQ51ww0GTMogYJe/WB8zsrrY7UeAOglrwrPjfgK0kbQW1lI2dM071yKmaMrMXFZclbUGoXpdMVNL8DbBBXPU/wH9VzM0HeQ2hWYgD+7mFgepLGVf9LLUZ8rhOw2Nmdk/f4LzqRCbXYHp9wr1sNSaqPO4jTGxKo8mrVbbyG05D9vMtGgSTjYKT4DkI08gdxnsQQTW7BnBEL+Re0k5AkoE7Go7vBxZLmhdX3w983vJUXK6MpL3N7HRgnckcn2Y2U9TKVGxrZptI+omZHSzpcEIexyZIvW9mMz5kdhZnGYMUnDuDjP7ITZfmMf2kFq+aS3i27EFQ5nfluO6UtIB4/UjaHfhju10ao7Pj8S7iBkLHccqQM7w1NXQjZ7hhlxkl79Y1kq4nhApeVAi7aoIy5+5Cwrm1fvwrUkfZOEjjVHJOnIyqqf52l0RjThKSnksoaPFHwoRZBOPSlyVtb2bJqQUyh4/m4gsEo3Tvmj4T+ClBUbeE6ve3n0naE1hO0nrAgcAPK7aVZTBtZt8BviNpGzO7umJferRRrTKJXOebpFuZ2khhZragQvc8B2EaWcN4zexHjDs+iusvBC4stF9K6WRmxwHHxf4wmWMhs2qqLL08mfOm/VQ6D8X/D0p6DnAX8OzM+6hKTuNDTmdxrjFI0bkzqOiP3OScx6Qa9LIZfSWdQ4gUuBg4BlhkZk+m9KdBDiA48TZQyG9/K7B3u10ao3PO4i6jZudljuPMBiQtsVCNLkdbS81s8xxtdR1J95UNl5R0M/AiM3ssLq8I3GBmG8y270xB4rQDIZ/Wi4FzgJPN7BcV2pp24CXpmFzqyyqTrahC7BkLrioYC7LS5jnS9xvMIXi0n25mScYbSScD15vZkX3rDySEfL+jbl+7gEKeyxf3VJG93y5eFz8ws+2mb2HKdlcGPsb45O2/gM9UyU0k6eOEfK/FwfT5wOHACWaWZFyVdApwkJndHZdXBw63hGJT/c+hnM+lriHp6X2r5gBvJRhIl5jZboPv1WgR1ZavJhh5+8N4rzCzLzW035zjrVlzjcR70tGE3+SrBMPU18zs4w3sK+l5mvN+KWkxsFPBWTyP4Cx+HWFcs2HZtuL2AxmDdJ2UayH3uFIhD99kRt91gIVmlqIw3RG41EL+0k4SjW5zakRDNIJfC+VxA6HjODPSxoC1wTC8gZEyyMw9IR8WJG0PnE5QHdxASKRfWmmUc+BVYl+dnWy12beodu3xOPBr4LxUw5Skm81sGYVNfO/nZtav7hxKJN1gZpsWll9rZpfE19eb2Wbt9W6cnIPpye6FFSbhtzAxb9thxeUhCXNLQqG4yD6ExO/XA4emhlVrYu7HZZiN31suJO0PfJRxNVzjYbw5nT1tOI4kHTXd+2Z2YIZ9rAisZGb3VNw+u2Mx1/2yy87iHI6etkg0EGYdV+Y0+kZH4PuBtczsH2K0wPpmdsEMmw4EDa6QkNMgHmLsOE4ZksMXM9BUGN4gKe2B6WgoZCNEdczehInv7YRk8OcDmxHCfVMKKwwyNHs2VZLOhpkdDGHgamYP1mjqoWneq9Nu13iKpKf2nBwF4+CqVCzcEbf/HvCWvsnbWalKzh6ZQ9DnSFrdzP4S+/Y00segbVSrbAWFHHf7Ae8jFOV6o5lVrZy58zTvzarvLTcthfHmVG60oQK5Lv5/GbAh44Vx3kKFUD5NUk3azB6R9FZJVapJQwM5nzPeL7scCrlJ7/kCYGZ/kTQs0S0p85jc48qcIegnxT5sG5d/Txg3t24g1GALCTkN4gZCx3HKUHViMhlljSyvJoSf9rjbzHbuheFl7E9nyDwh7zJXA6cRJr2/K6y/Ng4wUhhk4uEuS+5bM15K2oZQOXMesJakTYF3m1lqoZL+Srdju6AbFfpycSJwtkIF0tsAJK0NHEsYVFflGZNM3rpS0Olw4GpJCwm/5+7AZ1MasHaqVbbFrQQ17pHAbcAmkjbpvZmi+hux760RZohY6HLxg4HTO4ck/SOwXSGVwnFUG7vlriYNHc753HFncQ5HT1ukzGNyjytzGn0XmNkekt4OYGYPxnlRFxhkISGnQYblonYcp11KG0ZmCt2gfDLeOTaxcumHIWRn13jVvq7TlYd211jfYn6LqHK6u7dsZl9IbKvL3vZsZLyumuBIQrGI8wHM7AZJL6/QTr9CrMjiin3rHGb2ZUkPAlfG8xTyhC4+KWmtPqNjJ4zaZnaqpOsIygKAN6eGyiYwGww2lxJ+u03jX5GmVH+z4Xtrg9zP+ZwRG21Ef/RYneDY6VXJnhfXpZK7mjR0vKJph53FtR09LZLyLMw6rsxs9H1UoYpxbwy9gInncpt0uZCQk4AbCB3HyU2u0I1GwvBy0nGjTeeQ9AngHDO7OebVuZgw+X1c0p5mdmlqmwP2trc52coeEpUTM/ttnxM7OYH2KCmdGgpd/BjB6LiIMHn7a+DdOfqbAzP7maQ/E+/fRWNmZobeMWNm+5b5XOZrYei/t5bIbYTPGbGRs61UPg8slXQ54dx6OfCpCu1krSYdGQnHYm4G7OhpjSbGlRmNvp8kjJ3/StIZhD7um6HdHFwgaTXgS4RUUEaImHCGDC9S4jjOjKQkRs6VjDca33YgPJj7w/AuM7PD0o8kL7kTGc92JP0M2DiqQP8BeDvhN34BcIqZvaTVDs5AmwnCcya5bqBv5wJfBo4BtiYokbYys7c1tL/OFovJSZXjlPQM4KVx8RrgHovJ7ttE0i4E9clzCLll1wZuMrONGtjXSJwf4BVvu0Du7202/aaS1iQ8EwB+ZGZ/qtBGI9Wk5RVNKxNTVxSLUDTh6MlKm+O3XCgUrdod+D7hOS/gGjO7s+V+vdfMjomvN4rOwFqFhJx2cQWh4zi5lXBZQjcaDMPLSWfz2HSUR3uhxISQ1LPM7AngJknD8Dxq06PW5ZCo/YGvAM8lJMy+BDigwf2NitIp+TjN7E5J/wm8ilDo6Q3As3J3rAKHECY0l5rZ5grVy/duuU+zgVG5FrpMm8ryzhLzou0APN/MPi1pLUkvMbOkogVmdpik+4HFhfQytceCHQ7j7SxTOXoIVWtbZRQieszsSUkfMrNzCHONrrAfwUEMIb/4Fmb2CN0JfXYSGYYJmeM4zZMzfDFb6EZDYXg56bLRpos8ImljQuXi7YEPFN5buZ0uDQ2dDYmK3utBJk8fldCHpOOU9FJgT+CNwNMIRtoPTLvR4HjMzO6SNEfSHDO7XNKRDe1rlAw2Oa+FUfrecpI7jDen0bdNA/K/A08SnBWfBu4DzmNi8blSDMFYcFTosqOn02lYMnJpVNWeDYyF3ZvZ/069yUBxp9UswA2EjuNARiVcQ7k7ulpBsLNGm45yEHAusAZwhJndCiBpJ2Bpmx0rSWsDnwHnWkxC0heBzxASVF8MbAK8z8xOb7Nfo4KkQ4G3EKrdngkcDFzbscny3VH9sxg4Q9IdFCY3mWkz79ow499bNVIN+dmUTh1XTW1tZltIWgpjVdWr5gwkttHVseCoMEhHTyqjEtGzR/xfjNIw4Pkt9KXHapLeBMwB5kt6c/FNM2uiqJbTIG4gdBwHMivhBhy64UabIcHMfgRsMMn6C4ELe8ttKQE6PtnqckjUa83sQ3GA+GvgzQRDUFMGwlFROpU9zncBvyDkZ/2umT0iqWsqy10JBuT3Ee6NqxIMmU3QtWNvkpzXwih9b22SU+nUZdXUY5KWY7za6hoERWFTuHKpeQbp6EllJCJ6zGzdtvswCYuAXeLrxQSxRA8D3EA4ZLiB0HEcGG4lXKuTmg4bbYaZtpQAXZ5sdZneWOL1hOI89/RVNM7NqCidyh7ns4HXEIr+HBmrhs6VtLyZPd5Y7xIoVCB9knhtS7qKceeKU41RuRa6TOrNLqfSqcuqqaOAbwPPlPRZQnGFf2twf27gbp5BOnpSGeZ5zIxI2ho4AVgA3AjsZ2Y3tdurgJn9XZnPeRqA4cENhI7juBLO6RptKQG6PNnqMhdIupkwcfjHqBR5uMH9jcpEsNRxxkI/FwMXx8qBbwDmAr+X9H0z27PBPtZhrYbaHSUlUc5rYZS+t9I0oCzPqXTqrGrKzM6QdB2hArGAN3bFoOFUo8uOnhGYx3yVkFN4MUGtdySh2N8w4WkAhgQ3EDqOAwy1Em5Uwg1HibYMQJ2dbHUZM/tIzEN4j5k9EauP79p2v0aRWDnwPOA8SfMJBUuATnrvK13nXU8F0FX8e6tMbmV5TqVT51RTkp5WWLyDkBd17L0Giyn4WLAdmnL0JDPE85gyzDGz78XXCyV9tNXeVMOdUEOCGwgdxxl2PMTKyUXnJltdZrKQF+CmqDJoMi/RqAwyax2nmd0LnFpYNXDvfX+y8uJbBJVjFTwVwDgp54h/b9XIqizPqXTqqGrqTuB3QC+9QfEcbbKYgo8F22FUFP1ts1rf83TC8pAUAvFzZUhwA6HjOMOOP3BmH60oATo62eoyjYS8jIrSqYXjbMOwuvM0711Qsc2RSQWQ+RwZme8tM9mV5TmVTh1UTR0FbE84N88ErjSzQYzTfCzYEA05epw0FjHxeVpc9kIgTlbcQOg4juN0jdaUAB2cbHWZpkJeRkXpNOjjHPgEuqHk5aOUCiDnOTJK31tOXFmegJn9s0KVqlcC+wBHS7oEONbMbm21c05VmnD0OAnMkkIgngZgSHADoeM4w86ohBuOEq4EGA6aCnkZFaXTqBxnGVLCn0fJYJPzHBml7y0brixPJyoGL5e0FHgbcAjwS+DEBnfrY8GGmCXGqVGhy4VAPA3AkOAGQsdxOs2ohBs6zhDSVMjLqCidBn2cXfbel57cj5jBJts5MmLfW1ZcWV6eaHDeFdgDWIPwHNjSzG6r2a6PBbtPl41To0KXDeXu/B8S3EDoOE7XGZVwQ2ecLg9wnEiDqoJRUToN+ji77L1PmjiMkMEm6zkyQt+b0x53ENSCZ8X/BmwlaSuopSz3sWD38bFb+7gRzqmNBpM31nEcpxqSFgM7FUKs5hFCrF5H8Bxv2Gb/nHRmUgJIOsYH+7MHSUvMbIvEbbZiXOl0VUHpNKsY5HFW+R0GRZf71jajci04swNJJzO1kcLMbL+K7fpYsOP4fbx9uvwbSFpqZpu33Q9nZlxB6DhO1xmVcMNRwpUAo0WyqmBUlE6jcpwl6HL4c6v4OeIME2a2b5nPVVCW+1jQcWamtWeppwGYPbiB0HGcrjMq4YajhBdnGC08VKEbdDn8q8vhz47j5Cc1X52PBbuPO3rap81nqTv/ZwluIHQcp9N4YvVZiSsBHCczQ+69dyOy44wWSQ4LHwsOBe7oaZ82n6Xu/J8luIHQcZzO4yFWsw5XAowWrioYDO69dxxnWEg2ZPhYsPO4o2e0cef/LMENhI7jOM5AcSXAyOGqgsEwzN77Loc/O47jOM4w0Oaz1J3/swQ3EDqO4zgDx5UAI4WrCgZDZ733Qx7+7DhOflxZPvtwR0/DdPlZ6s7/2YMbCB3HcRzHcYafLnvvPfzZcZwiriwfMrpsnBohOv0sdef/7EBm7th3HMdxHKcZJC01s83b7scoIGkrxr33VxW8960iaTGwUyH8eR4h/Pl1hMnlhm32z3GcwSJpiZlt0XY/nPJI+iaTG6fWARaaWZfTWMwK/FnqDAJXEDqO4ziOUxlXFXSHDnvvOxv+7DiO45RimPPczhb8Weo0jhsIHcdxHMepQ6dDXpxO0OXwZ8dxBo/nqxs+3DjVPv4sdRrHQ4wdx3Ecx6mMh7w4Zehq+LPjOPmZSVku6Rh3Hg0Xkj4OvAkoGqfOBw4HTjAzL0QxAPxZ6jSNGwgdx3Ecx6mMpJuBF5nZY3F5ReAGM9vA8w86juOMHp6vbnbixinHmf14iLHjOI7jOHXwkBfHcRyniOerm4V0OM+t4ziZcAWh4ziO4zi1cFWB4ziO08OV5Y7jOMOJKwgdx3Ecx6mFqwocx3GcAq4sdxzHGUJcQeg4juM4juM4juNkw5XljuM4w4cbCB3HcRzHcRzHcRzHcRxnhJnTdgccx3Ecx3Ecx3Ecx3Ecx2kPNxA6juM4juM4juM4juM4zgjjBkLHcRzHcRzHKYmkz0naXtIbJX20sF6S/k3SLyX9QtLlkjaK7/1I0vWSbpP05/j6eknrSLq/r/19JR0zyX7H1kv6lKTfxzZ+KelbkjZs+tgdx3Ecx5m9uIHQcRzHcRzHccqzNXAN8ApgcWH9AcC2wKZm9gLgc8D5klYys63NbDPgE8DZZrZZ/Pt1jX4cEdtYDzgbuEzSGjXacxzHcRxnhHEDoeM4juM4juPMgKQvSfoJ8GLgauBdwLGSPhE/8mHgvWb2IICZXQL8ENir6b6Z2dnAJcCeTe/LcRzHcZzZyfJtd8BxHMdxHMdxuo6ZfVDSOcDfAu8HrjCzlwFImg+sYma39G12LbDRgLq4BNhgQPtyHMdxHGeW4QpCx3Ecx3EcxynHFsANBEPcTQ3uxypso+y9cBzHcRxnZHAFoeM4juM4juNMg6TNgJOB5wF3AiuH1boe2MbM7pX0gKTn96kItwQWzdD8Q5KeYmaPxuWnxX2ksjlBseg4juM4jpOMKwgdx3Ecx3EcZxrM7PpYZOQXwIbAZcCOsUjIQ/FjXwKOkjQXQNIOwHbAN2dofhGwd9xmLvBW4PKU/knaDXgtcGbKdo7jOI7jOD1cQeg4juM4juM4MxArBP/FzJ6UtIGZ/U/fR44GVgdulPQE8Cdg14IBcSoOAo6XdCAhTPhUM1s8wzYA75O0N7AK8FPgVWb255RjchzHcRzH6SGzKilOHMdxHMdxHMdxHMdxHMeZDXiIseM4juM4juM4juM4juOMMG4gdBzHcRzHcRzHcRzHcZwRxg2EjuM4juM4juM4juM4jjPCuIHQcRzHcRzHcRzHcRzHcUYYNxA6juM4juM4juM4juM4zgjjBkLHcRzHcRzHcRzHcRzHGWHcQOg4juM4juM4juM4juM4I8z/AcXC/iEjVJp5AAAAAElFTkSuQmCC",
            "text/plain": [
              "<Figure size 1296x432 with 4 Axes>"
            ]
          },
          "metadata": {
            "needs_background": "light"
          },
          "output_type": "display_data"
        }
      ],
      "source": [
        "import numpy as np\n",
        "import pandas as pd\n",
        "import seaborn as sns\n",
        "\n",
        "abundances = pd.read_table(\"genus.tsv\", skiprows=1, index_col=0)\n",
        "abundances.index = abundances.index.str.split(\";\").str[5]       # Use only the genus name\n",
        "abundances = abundances[~abundances.index.isin([\"g__\", \"__\"])]  # remove unclassified genera\n",
        "abundances = abundances.sample(50)                              # use 50 random genera\n",
        "\n",
        "# Let's do a centered log-ratio transform: log x_i - log mean(x)\n",
        "transformed = abundances.apply(\n",
        "    lambda xs: np.log(xs + 0.5) - np.log(xs.mean() + 0.5),\n",
        "    axis=0)\n",
        "\n",
        "sns.clustermap(transformed.T, cmap=\"magma\", xticklabels=True, figsize=(18, 6))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "QbRoCh1BR0GZ"
      },
      "source": [
        "# Exercises\n",
        "\n",
        "Okay, that's enough time in the back seat. \n",
        "\n",
        "It's time to take the wheel 🚗 \n",
        "\n",
        "Now you can dive into the data 🏊"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "hy5HN4ZoS2-9"
      },
      "source": [
        "## Exercise - Plant a Tree\n",
        "\n",
        "One visualization that we did not spend a lot of time on was the phylogentic tree of our ASVs. Let's change that! We have seen that there are genera that appear in multiple populations in the previous step. But are the organisms in that genus actually the same?\n",
        "\n",
        "Let's annotate the tree with our taxonomic classifications and abundances. We will use the empress plugin again but this time with the `community-plot` option. I filled in a template of the command for you. Can you figure out what has to go in the empty spaces?\n",
        "\n",
        "**QUESTIONS:**\n",
        "\n",
        "1) Are some of the branch lengths on the tree longer than you would expect? Do you notice anything interesting or suspicious about the taxonomic identities of these branches?\n",
        "\n",
        "2) Can you find examples of phyla that are polyphyletic (i.e. where clusters of ASVs from the same phylum are found in different locations on the tree, showing different commmon ancestors)? What about polyphyletic taxa at lower taxonomic levels, like at the family or genus levels? Why do you think these patterns exist?\n",
        "\n",
        "3) Use the \"Shear Tree\" functionality to only visualize the phylogenetic tree of the genus *Prevotella*. What do you observe? Can you figure out which of the unique V4 sequences are found only in one population? What impact might this have on studies trying to link *Prevotella* to human health?"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 19,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "tq9TWRMjkER5",
        "outputId": "ba14ac27-6b9e-4207-9347-27e373e6a328"
      },
      "outputs": [
        {
          "name": "stdout",
          "output_type": "stream",
          "text": [
            "\u001b[32mSaved Visualization to: community-tree-viz.qzv\u001b[0m\n",
            "\u001b[0m"
          ]
        }
      ],
      "source": [
        "# This won't run until you fill in the [EMPTY] spots with the right files ;)\n",
        "\n",
        "!qiime empress community-plot \\\n",
        "    --i-tree tree/rooted_tree.qza \\\n",
        "    --i-feature-table dada2/table.qza \\\n",
        "    --m-sample-metadata-file data/metadata.tsv \\\n",
        "    --m-feature-metadata-file taxa.qza \\\n",
        "    --o-visualization community-tree-viz.qzv"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EaH5k-SV_XiW"
      },
      "source": [
        "This is what the command looks like with everything filled in. This command will create the same tree as before, but you can now annotate it with a lot more information, and hovering over a node will now also show its taxonomic classification. Another thing that can be helpful is to collapse nodes on a particular taxonomic rank. For instance, you could collapse on the genus level to reduce the large number of external nodes (and the unclassified ones). \n",
        "\n",
        "![Example visualization](https://github.com/Gibbons-Lab/isb_course_2021/raw/main/docs/16S/assets/community_tree.png)\n",
        "\n",
        "Let's get to the questions:\n",
        "\n",
        "**1)** There are definitely some very long branches, what could be happening? Hovering over those very separate nodes you see that some of them have no classification, sometimes not even on the kingdom level. Interestingly, bacteria are not the only organisms that have a 16S gene. But the GreenGenes database mostly covers bacteria. The ones not classified on the kingdom level are thus likely sequences coming from 16S genes from another Domain, such as archaea. However, many of these very long branches are still classified as bacteria, and in many cases get classified as mitochondria or choloroplasts. First, of all it's really hard to distinguish mitochondria and chloroplast based on the V4 region alone. These organelles have diverged a lot from other bacteria, however they are still related to bacteria, which is somewhat consistent with [endosymbiotic theory](https://en.wikipedia.org/wiki/Symbiogenesis). Mitochondria, chloroplasts, and plastids often still contain their own genomes, seperate from their eukaryotic host cells (i.e. plants, animals, fungi, protists), and their internal 16S genes are still close enough to bacteria that the V4 primers will amplify them.\n",
        "\n",
        "**2)** In most cases you would expect that taxonomy should match the tree that you obtain by aligning the V4 sequences. After all taxonomies like NCBI, GreenGenes, or SILVA are often built by aligning the the core proteins, full 16S gene, or other parts of the genome. These molecular taxonomies are often designed to be monophyletic. However, initial taxonomic classification of microbes in the 19th and 20th centuries was based upon phenotypic assays. And as we see in the metabolic modeling course, sometimes unrelated taxa can show similar metabolic functionality. This is why you will find, for instance, genera like `Prevotella_1`, `Prevotella_2`, and `Prevotella_3` in SILVA. The sequences from organisms classified originally as Prevotella fell into polyphyletic groups (i.e. they didn't all share the same common ancestor) and were thus split up further into monopyletic clades. When building trees from amplicon sequencing data, you will often see examples of polyphyletic taxa.  You can further explore these polyphyletic taxa by coloring nodes by their environment, and you will often see that many of those polyphyletic groups are found in distinct environments.\n",
        "\n",
        "**3)** Shearing the tree you will that there is a Prevotella sub-tree with quite a lot of structure to it. Though, of course, no Prevotella ASVs are completely distant from each other there are some branches with quite a bit of phylogenetic distance between them, forming distinct clades. You will find both cases. Some ASVs in the Prevotella genus can be found in individuals from all 3 ethnic groups but many of them appear only in a single ethnic group and in many examples only in a single individual. So if there is already notable differences in the V4 region what does that mean for the rest of the genome? Well it *may* show quite some differences as well meaning that those clades are not the exact same organisms. So they may very well behave differently from the Prevotella strains observed in North-American and European populations, having distinct or maybe even no connection to human health in those specific ethnic groups. Thus one can not generally extrapolate microbiome-health connection from overrepresented groups to other populations and communities.1"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tFuoGkWtptmR"
      },
      "source": [
        "<br><br><br>\n",
        "\n",
        "---\n",
        "\n",
        "# 🦠 Space for your own exploration\n",
        "\n",
        "Feel free to use this space to run some additional analyses. If you rerun all the previous cells in the notebook (`Runtime > Run All`) you should be all set to start working on the previously generated artifacts. Remember that you can get help in the Slack channel and also in the [Qiime 2 documentation](https://docs.qiime2.org)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 20,
      "metadata": {
        "id": "pq919yDYqYoe"
      },
      "outputs": [],
      "source": [
        "# You can add more code cells with the \"+ Code\" button on the top right\n"
      ]
    }
  ],
  "metadata": {
    "colab": {
      "collapsed_sections": [],
      "include_colab_link": true,
      "provenance": []
    },
    "kernelspec": {
      "display_name": "Python 3.10.7 ('micom')",
      "language": "python",
      "name": "python3"
    },
    "language_info": {
      "name": "python",
      "version": "3.10.7"
    },
    "vscode": {
      "interpreter": {
        "hash": "c991a7ed881363492957ff225bb30af9d5174cd8515a21cbef71fcaa303e4050"
      }
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}