transcript_novelty_barplot.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "1ec70275-884c-4a33-b16e-3522257702d9",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message in system(\"timedatectl\", intern = TRUE):\n",
      "\"running command 'timedatectl' had status 1\"\n",
      "-- \u001b[1mAttaching packages\u001b[22m --------------------------------------- tidyverse 1.3.1 --\n",
      "\n",
      "\u001b[32mv\u001b[39m \u001b[34mggplot2\u001b[39m 3.3.6     \u001b[32mv\u001b[39m \u001b[34mpurrr  \u001b[39m 0.3.4\n",
      "\u001b[32mv\u001b[39m \u001b[34mtibble \u001b[39m 3.1.7     \u001b[32mv\u001b[39m \u001b[34mdplyr  \u001b[39m 1.0.9\n",
      "\u001b[32mv\u001b[39m \u001b[34mtidyr  \u001b[39m 1.2.0     \u001b[32mv\u001b[39m \u001b[34mstringr\u001b[39m 1.4.0\n",
      "\u001b[32mv\u001b[39m \u001b[34mreadr  \u001b[39m 2.1.2     \u001b[32mv\u001b[39m \u001b[34mforcats\u001b[39m 0.5.1\n",
      "\n",
      "-- \u001b[1mConflicts\u001b[22m ------------------------------------------ tidyverse_conflicts() --\n",
      "\u001b[31mx\u001b[39m \u001b[34mdplyr\u001b[39m::\u001b[32mfilter()\u001b[39m masks \u001b[34mstats\u001b[39m::filter()\n",
      "\u001b[31mx\u001b[39m \u001b[34mdplyr\u001b[39m::\u001b[32mlag()\u001b[39m    masks \u001b[34mstats\u001b[39m::lag()\n",
      "\n"
     ]
    }
   ],
   "source": [
    "library(tidyverse)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "a203ba6a-b1a0-4de9-af41-b00a7207dff8",
   "metadata": {},
   "outputs": [],
   "source": [
    "read_annot = read.delim(\"data/all_talon_abundance_filtered.tsv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "eb2c2489-185c-4185-8c4d-3c871f5b1147",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAIAAAByhViMAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nOzdeUAUdePH8e8Awhq35C2KCngkamamZo9oiFpq5oMJnpBnmokPWqIl0EGl\nxINaKWiCilee9Wg+hheaZpplaGIeaCGapAhyKsvu74/9PRuB4Ii7DA3v11/DzOzsZ9ZaPnzn\nkvR6vQAAAMDfn4XSAQAAAGAaFDsAAACVoNgBAACoBMUOAABAJSh2AAAAKkGxAwAAUAmKHQAA\ngEpQ7AAAAFSCYlfjFBcX37p1q7CwUOkgCrt161Zubq7SKRSWm5t769atWn4X8aKiolu3bt29\ne1fpIEoqKSm5detWQUGB0kEUlp2dnZOTo3QKoEazUjoAytLr9SUlJbX8d7kQoqSkRJIkpVMo\nTKfTlZSUKJ1CYYb/I5ROoTDDh6DT6ZQOojCdTsfXAlA5RuwAAABUgmIHAACgEhQ7AAAAlaDY\nAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAA\nqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATF\nDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAA\nQCUodgAAACpBsQMAAFAJih0AAIBKWJl164cOHfryyy8zMjLu3Lnj4uLyzDPP+Pv716lTx7D0\n+++/X7NmzZUrVxwdHX18fAICAiRJMtMiAAAA1bMMDw8339YzMjKaNWv23HPP+fr6NmjQ4PPP\nP79582a3bt2EEL/88sv8+fN79uw5bdo0V1fX1atXFxcXd+zY0RyL/l5KSkru3LlTp04dYwOu\nnQoKCiwsLDQajdJBlHTnzh2dTvfII4/U5j9RtFptcXGxjY2NpaWl0lkUo9PpioqKrKysrK2t\nlc6ipMLCQkmS6tatq3QQoOYy74hdz549jdNt2rT59ddfU1JSDD9u3bq1adOmkydPFkK0aNHi\n2rVrX3zxxfDhw21sbEy+6GF2QYqf9TAv/3vRB0UpHQEAAFRdNZ1jp9Pp0tLSTp482alTJ8Oc\n1NTULl26GFfo0qVLUVFRWlqaORbdk14ek30EfwcyP5PqUQMjVT8+BAM+BD3/Mej1+hr/IVT7\ndzZwD+YdsRNCFBcXDx8+3PAfva+v76RJk4QQer0+Ozvb2dnZuJphOisry+SLKgp2+/bt4uJi\nE+6pCty8eVPpCH+h1WprWiRFVPKfce2Rm5urdATlFRUVFRUVKZ1CeTX2a8He3v4hjxEBD8/s\nxc7KymrRokXFxcXnz59PTEx0cHAYO3asud9UDisrs+/7306NOquvuLhYkqRa/s+k1Wr1en2N\n+nepfjqdrqSkxMrKqjafaKjX67VarYWFRW0+0VAIodVqRQ3+9raw4EYTUJ7Z//eQJKlFixZC\nCHd3dwsLi08//XTYsGF2dnZOTk63bt0yrmaYrlevniRJpl1UUTBbW1sT7qY6ODo6Kh3hTzdu\n3LC0tKxRkapfTk5OcXGxg4NDbe40hYWF+fn5jzzySG2+bkCr1WZnZ1tbW9vZ2SmdRUlZWVmS\nJNXyrwWgctX654Vh+MHwJ1e7du1++OEH46IffvhBo9G0atXKHIsAAABqA/MWu7i4uAMHDqSm\npv7888/btm1bvXp1165dnZychBDDhg3LyMiIjY399ddf9+/fv23btiFDhhjOTjD5IgAAgNpA\nMuuFPKtXr/7uu+8yMzMtLCwaNGjQu3fvwYMHG8vW8ePHExMT09PTDfcTHjlypPF4k8kXVRm3\nO1HKjRs3rKysDH8G1FqGQ7EuLi4cinVwcOBQrEaj4VCsJEmlr5MDUIZ5i50KUOyUQrETFDsh\nBMVOCEGx+x+KHXBfXMIDAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0A\nAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBK\nUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwA\nAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABU\ngmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIH\nAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACg\nEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7\nAAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAA\nlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDY\nAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAA\nqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATF\nDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAA\nQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUo\ndgAAACpBsQMAAFAJK6UDKEav1ysdocapgZ9JDYxU/Wr5h2DYfb1eX5s/B+O+1+YPwajGfgiS\nJCkdAajFxS4/P1+r1SqdombJyclROsJfaLXamhapmpWUlAghbt++rXQQJel0OiFEQUFBYWGh\n0lkUY6gyd+7cqeXfWoZ+X2O/FmxtbevUqaN0CtR2tbfY2dnZKR2hxnFyclI6wp9u3LhhZWVV\noyJVv5ycnOLiYkdHx9o8ElBYWJifn29ra2ttba10FsVotdrs7GwbG5ta/sWVlZUlSVIt/1oA\nKsc5dgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUo\ndgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAA\nACpBsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpB\nsQMAAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMA\nAFAJih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJ\nih0AAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0A\nAIBKUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBK\nUOwAAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwA\nAABUgmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABU\ngmIHAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUgmIH\nAACgEhQ7AAAAlaDYAQAAqATFDgAAQCUodgAAACpBsQMAAFAJih0AAIBKUOwAAABUQlaxy83N\nzcjIMP6YkZHxxhtvTJw48eDBg2YLBgAAgAdjJWelqVOnnj179vjx40KIgoKCHj16pKenCyHi\n4+MPHTrUo0cP82YEAACADLJG7A4fPjx48GDD9MaNG9PT09evX3/p0iV3d/eFCxeaMx4AAADk\nklXsrl+/7urqaphOSkpq3769v7+/m5vb+PHjDcN4AAAAUJysYidJUklJiWH68OHDvXv3Nkw3\naNAgMzPTXNEAAADwIGQVuxYtWhw4cEAIcfz48d9++61Pnz6G+RkZGfXq1TNfOAAAAMgn6+KJ\n0aNHz5079+rVq6mpqS4uLgMGDDDMP3HihIeHhznjAQAAQC5ZI3azZ88OCQm5fPly48aNP//8\nc3t7eyFEVlbWjh07jIdlAQAAoCxZI3ZWVlZRUVFRUVGlZ9arV+/OnTvmSQUAAIAHxpMnAAAA\nVKLCEbuEhASZmwgMDDRJFAAAADyMCotdUFCQzE1Q7AAAAGqCCotdUlJSdeYAAADAQ6qw2Pn4\n+FRnDgAAADykB7h4QqvVnjhxYteuXdnZ2eYLBAAAgKqRW+zWr1/frFmzrl27Pvfcc2fPnhVC\nXL16tUGDBomJieaMBwAAALlkFbvdu3ePGjWqWbNmCxcuNM5s0qRJx44dt2zZYrZsAAAAeACy\nil1kZGTnzp2PHj366quvlp7fo0ePn376yTzBAAAA8GBkFbsTJ06MHj3ayqrslRbNmze/du2a\nGVIBAADggcl6pFhJSYmNjU35+ZmZmXXq1KnkhXv27ElOTr58+fKdO3eaNGny/PPP9+vXz7j0\n+++/X7NmzZUrVxwdHX18fAICAiRJMtMiAAAA1ZNV7Dw9Pb/55ptp06aVnqnX67/88ssOHTpU\n8sJ9+/Y99thjL7zwwiOPPHLkyJElS5ZotdqBAwcKIX755Zd333134MCB//rXvy5evPjpp5/q\ndLrRo0ebYxEAAEBtIKvYjRs3bvbs2b6+vgEBAYY5eXl5s2bNOnbsWGxsbCUvjIyMNE63b9/+\n0qVLhw8fNhS7rVu3Nm3adPLkyUKIFi1aXLt27Ysvvhg+fLiNjY3JF1XxswEAAPhbkXWO3Wuv\nvda/f/+XX365RYsWQoixY8e6uLjExsYOHjx4woQJ8t/s7t27jo6OhunU1NQuXboYF3Xp0qWo\nqCgtLc0ciwAAAGoDWSN2VlZW//nPf2JjY1evXl1UVHTt2rUOHTqMHTv21VdftbCQeye8PXv2\nXLhwYdKkSUIIvV6fnZ3t7OxsXGqYzsrKMvmiivLk5OQUFxfLDF9L3LhxQ+kIf6HVamtaJEXc\nvHlT6QjKu337ttIRlFdUVFRUVKR0CuXV2K8Fe3t7jhFBcbKKnRDC0tJy6tSpU6dOrdrbHDp0\naNmyZTNnzvTw8KjaFkzO0tJSr9crnaJmKX/hs4K0Wq2oYZGqX0lJiV6vr+Ufgk6n0+l0lpaW\ntflaKL1eX1JSIkmSpaWl0lmUVFJSIoSosR+C/JEOwHyq4xfGrl27Pvvss1mzZnXv3t0wR5Ik\nJyenW7duGdcxTNerV8/kiypKZWdnZ6odVA0nJyelI/zpxo0bVlZWNSpS9TOMKzs6OtbmTlNY\nWJifn29ra2ttba10FsVotdrs7GwbG5ta/sWVlZVl+LZXOghQc5n9z4sNGzbEx8e/9dZbxlZn\n0K5dux9++MH44w8//KDRaFq1amWORQAAALVBhcVOI1slW1++fPnGjRuDgoLs7e3T0tLS0tLS\n09MNi4YNG5aRkREbG/vrr7/u379/27ZtQ4YMMZydYPJFAAAAtYFU0Xlmfn5+pX88c+ZMampq\n06ZN27RpI0nS2bNnMzIy2rVr1759+82bN1e09VGjRuXm5pae06hRo7i4OMP08ePHExMT09PT\nDfcTHjlypPF4k8kXVZkUP+sht/A3og+KUjrCnzgUK/53KNbFxYVDsQ4ODhyK1Wg0HIqVJKn0\ndXIAyqiw2JV25MiR/v37L168eNy4cYaTQ3U63cqVK4ODg7/++uuePXuaP6diKHZKodgJip0Q\ngmInhKDY/Q/FDrgvWefYzZkzZ9y4cUFBQcZLfiwsLCZMmDB27NjQ0FBzxgMAAIBcsordiRMn\nOnXqVH7+448//v3335s6EgAAAKpCVrGztrb+8ccfy88/ceIEVycAAADUELKK3aBBg2JjY+Pi\n4gz3jBVCaLXaZcuWLV++fPDgweaMBwAAALlk3aB44cKF33333eTJk+fNm+fh4aHX68+fP3/z\n5s02bdosWLDA3BEBAAAgh6wRu0aNGp04cSI8PLxJkyYpKSmnTp1q2rRpRETE999/37BhQ3NH\nBAAAgBxyHylmb28fFhYWFhZm1jQAAACosgd7Vuzt27cvX74shHBzc3NwcDBLIgAAAFSJ3GfF\nnj17tn///s7Ozp06derUqZOzs/OAAQN++eUXs4YDAACAfLJG7C5cuNCzZ89bt2716NHDy8tL\nCHH69Ondu3f36NHj2LFj7u7uZg4JAACA+5NV7ObPn19QULB7925fX1/jzK+//nrIkCFhYWFr\n1641WzwAAADIJetQ7J49e6ZOnVq61QkhfH19X3nllT179pgnGAAAUIk9e/ZIkpSQkKB0ELNT\nfE9lFbvs7GwPD4/y8z08PLKzs00dCQAA/L9z586Fh4enpKQoHaS6qWPHq38vZBW7Jk2aHDly\npPz8I0eONGnSxNSRAADA/zt37lxERMTfvd/07du3sLBwzJgx8l+ijh2v/r2QVeyGDRuWmJj4\nwQcfFBUVGeYUFRVFRkauXbt22LBh5owHAADur6CgQOkI92YIZmFhodFoLC0tlY6jfrKK3fz5\n8zt06BAaGlq/fv3HH3+8c+fO9evXnzdvnpeX11tvvWXuiAAA1E7h4eGGZ7KPGTNGkiRJkry9\nvYUQmzdvliTp888/j4iI8PDwsLa2fvvtt4UQOTk5b7755lNPPfXoo4/a2Ni0atVq1qxZeXl5\nxg0aXrhly5YPP/zQ09PTxsamefPm7733nl6vN66j1WoXLFjg5eVlb29vb2/v4eERGBiYm5tb\neoWYmJgnnnjC1tbW3t6+Y8eOxucXVBSszJlnhtU2bNgwb948Nzc3GxsbDw+PmJiY++54JeTs\nmhAiOzs7JCSkZcuWNjY2DRs2HDVq1IULFwyLkpKSJEl67733ymw5MDDQ0tIyPT3duPvR0dGd\nO3euW7euvb29t7f3119/Lf+fb//+/ZIkvfPOO2VWHjNmjJWV1ZUrVyrfzfuSdVWsk5PT0aNH\no6Kitm7dev78eUmSWrVq9c9//jMkJMTW1vYhEwAAgHsKDAy0sbGZO3fu3Llz+/XrJ4RwcnIy\nLn3jjTeaNm0aGRnZqFGjOnXqCCHS09Pj4uL8/PwCAgKsra0PHjwYHR197Nix5ORkSZKML3z9\n9dc9PT0XL17s5OS0YsWKN99808XFZcqUKYaloaGhUVFRI0eOfO211ywsLH799dcdO3bcvn3b\n3t5eCKHVagcNGrR79+7evXvPnz/fwcHh7NmzmzZtioiIqCTYPc2aNeuJJ57YvHmznZ1dQkLC\nzJkzr1+//v777993xytR+a7l5+f/4x//OHXq1KhRo3r27Hn+/PmlS5fu2rXr22+/bdOmzbPP\nPuvq6rpq1ap58+YZN5ifn79ly5a+ffu6uroKIUpKSoYMGbJ79+7hw4dPmDChqKgoMTFxwIAB\na9euDQgIkPPP17lz5zZt2nz22Wfz5s2zsPj/8bXs7OwtW7YMHDiwWbNmcnazEnKfPGFra8sj\nxQAAqE5ubm6G28e2a9eu/JCVtbX1gQMHrKz+/FXu4eGRkZFh7FJTp07t2LHjvHnz9u7d6+Pj\nY1ytXr16X331laHqdevW7eDBg4sXLza2ny1btvTp06f0vcwMo24GH3/88e7du6dPn75o0SJj\nWdTpdJUHuydra+stW7YYVvvggw9+/fXXBQsWjB8/3t3dvfIdr0Tlu/bRRx+dOnXqvffemzt3\nrmHOwIED+/fvP2PGjP/+978WFhZjxoyJjIw8cuRIz549DSts2rQpLy8vMDDQ8GNsbOyuXbtW\nrlwZFBRkmDNjxozu3bvPnDlz+PDhZXa5or2YNGlSSEjInj17jPcbSUxMLCwsnDRpkvw9rYjc\nJ08AAIAaJSgoqEyTsLGxMba64uLioqKiF198UQhx9OjR0qsZjgwapi0sLLp27Xrx4kVjOXNy\nckpNTT1+/Pg93zQxMbFu3bqRkZGlhwCNI08VBbunwMDA0qtNnDhRp9Nt3779vi+sROW7tmXL\nFjs7u3/961/G9X19fXv06JGUlHT79m1DJCHEqlWrjCskJCQ4ODgYryhYvXp1gwYNAgICiv6n\npKQkICDg+vXrP/30k8yQgYGBGo1m+fLlxjnLly9v1qzZc889V+UdN3qAZ8X+8ccfFy9evHnz\nZpnD1YMGDXr4HAAA4IG0bNmy/MyEhIS4uLiffvqp9OUUWVlZpdcxHFU0cnBwuHv3bm5urqOj\noxAiKirqpZde6tatW/PmzXv16uXj4zNixIhHHnnEsPK5c+fc3d3t7OweNFh5rVu3Lv1jq1at\nhBAXL16U89qKVL5raWlprVu31mg0pdfx8vL69ttvL1++3LFjRw8Pj549e27cuHHRokUajeby\n5csHDx4cP3583bp1DSunpqbevn3b+GNpmZmZMkPWq1fPz8/v888//+OPP+rXr3/s2LGUlJT5\n8+eb5OISWcXu1q1b06ZN27hxY5mxVoMyPQ8AAFQDGxubMnOio6NDQkIGDx68YsWKJk2a2NjY\n3Lx5c9CgQWV+fZcebDMy/jbv27fvpUuX/vvf/+7fvz85OXndunVhYWHffvtt06ZNDavd8+WV\nB7unO3fulP/xvhuvXOW7Jid8YGDgpEmTtm/f7u/vv2rVKr1ebzwOK4TQ6XQeHh6rV68u/8K2\nbdvKzzllypTExMTVq1eHhIQsX77cwsJi/Pjx8l9eCVnF7pVXXtm4ceOLL77o7e1dr149k7wx\nAAC4rwcqOp999lnLli2/+OIL46sOHTpUhTe1t7cfPnz48OHDhRAbNmwICAhYvHjxhx9+KIRo\n06bNmTNn8vLyKh+0k+P06dPlfzSM24mHbnj31Lp16wsXLhQVFZUetDt9+rSFhYWbm5vhxxEj\nRsyYMWPVqlUjRoxYvXq1h4fH008/bVzZ09Pz9OnTHTp0kLn7Fe3F008/3aFDhxUrVkyaNGnD\nhg39+/dv3rx51XesFFnn2O3cuXP06NFbt2597bXXRpdjkhwAAKA8w7WoZY6lVsTCwkKv15eU\nlBh+LCkpiYyMfNB3LPNe3bt3Lz1z9OjRhYWFZW52VrVjd/Hx8b///rthuri4+KOPPpIk6YUX\nXjDMeaAdl2nYsGF5eXml76uyZ8+eI0eO+Pj4ODg4GOY4ODi8+OKLSUlJGzduTEtLGzduXOkt\njB079u7du7NmzSqzy1evXr3nO1ayF5MnTz579uyrr76al5dnkssmDGSN2FlaWnbt2tVUbwkA\nAGTq1KmTRqNZsmSJtbW1k5NTgwYN+vbtW9HKfn5+4eHhAwcOfOmll3Jzczds2FCFytWkSZNB\ngwY98cQTTZs2zczMXLFihaWlpfGhEdOmTduxY0dMTMzJkycHDhzo4OBw/vz53bt3lxl+k6N1\n69ZPPfXUlClT7Ozs1q1bd/To0dmzZxsfYfpAOy7TrFmzNm/eHBoa+vPPPxtvd+Ls7Lxo0aLS\nqwUGBq5bt27y5MkWFhZjx44tvWjatGl79uyJjY398ccfX3jhhfr166enp3/77bc//fTTPc+x\nq2QvxowZ88Ybb6xevbpx48YmvFxBVrHz9vb+4YcfTPWWAABAJkdHx3Xr1kVERAQHB9+5c6d3\n796V9Jt58+ZZWVnFx8e/+uqrDRs29PPze+2112ReymAUEhJy4MCB6OjonJycBg0aPPnkk/Hx\n8T169DAsrVOnzq5du2JiYtasWRMWFlanTp2WLVsaDto+qLlz5168eHHZsmVXrlxxdXX96KOP\nZs6cWbUdl8nW1vbQoUNvv/321q1bN27c6OTk9OKLL7799tvu7u6lVzPc0C49Pd3Hx6fM1RhW\nVlZffPHF8uXLExIS3n//fa1W26hRo86dO0dHR9/zHSvZC0dHxxEjRsTHx7/88styLiKWSZLT\n5S9cuNCrV6/w8PBJkyaVuaRZ9aT4WUpHqD76oCilI/zpxo0bVlZWMu9IqVY5OTnFxcUuLi7m\nONfk76KwsDA/P9/BwcHa2lrpLIrRarXZ2dkajebhz2r6W8vKypIkydnZWekgeFibN28ePnz4\ntm3bhg4dqnQWJU2ZMmX58uUXL140nuH38GQ1RHd3908++WT48OGzZ89u0aJFmV558uRJU6UB\nAACoDW7dupWYmNi/f38Ttjohs9h9/vnnAQEBer2+bt26Wq1Wq9WaMAEAAEDl9Hp9mdujlGZt\nbf03OqJ48uTJU6dOrVy5sqCgwPgMDFORVezCwsJcXV137tz52GOPmfbtAQAA7uvnn382PJ7r\nntavX+/v71+deR5GYmLiRx991KxZs08++aRXr16m3bisYnfp0qXw8HBaHQAAeHh+fn4Permu\nu7t7RU85E+UeYlHDRUVFRUWZ66R2WcWuefPmd+/eNVMCAACAymk0Gu68JoesA9LTp09fvXp1\nXl6eudMAAACgymSN2Lm6ujZs2NDLy2vKlCmtW7cuc1VsLb9WGQAAoIaQVexefPFFw8ScOXPK\nL63ag0QAAABgWrKK3aZNm8ydAwAAAA9JVrHz8/Mzdw4AAAA8pL/N3fwAAABQOZM9dBYAAJic\naR9ZXqOeCQ5zYMQOAABAJSh2AAAAKkGxAwAAUAmKHQAAgErIKnbe3t4nT54sP3/fvn3e3t4m\nTgQAAIAqkVXskpOTs7Ozy8/PzMxMTk42dSQAAABUxUMdis3OztZoNKaKAgAAgIdR2X3sUlJS\nUlJSDNNJSUlXrlwpvTQrK2vJkiXt2rUzYzoAAADIVlmx27p1a0REhGE6MjKy/Ap169bdsGGD\nWXIBAAAlBAcHb9682Tiak5mZ+fzzz+fm5u7atatly5bKZsN9VVbsRo4c2bVrVzgZ/BgAACAA\nSURBVCHE4MGDIyMjvby8jIskSbK3t+/cubODg4PZMwIAACWcO3duwIABTZs23b17d7169ZSO\ng/urrNh5enp6enoKIcLCwgICAtzc3KopFAAAUNqRI0eGDBny7LPPrl692sbGxjDTz8/P0tKy\nXbt2K1euvH37dq9eveLi4po0aSKE0Ol077///ooVKzIyMpo1azZlypTZs2dLkrR3715fX9/M\nzEwXFxchRMOGDSVJ+v3334UQN2/erF+//t69e/v06VPJliGfrIsnwsPDaXUAANQeW7ZsefbZ\nZwMDAzds2GBsdQbbt2+3srI6d+7c5cuXb926NX36dMP8BQsWREZGvvXWWz///POcOXPCwsIW\nLVokhOjVq5e1tfX+/fuFEKdOncrLy8vNzT19+rQQYt++fRqNpmfPnpVvGfJVNmJXxvnz57/4\n4ou0tDS9Xt+6deuhQ4e6u7ubLxkAAFDE9evXX3rppZCQkAULFpRf6unp+eabbwohNBrNtGnT\npkyZIoTQ6/UffPBBSEjIyy+/LITw8PBIS0t79913g4ODbWxsnn766b179/r5+e3du/cf//iH\nXq/fu3dvhw4d9u7d26tXL2NxvOeW8UBkjdjp9frXX3+9TZs2s2fPXrp06bJly2bPnt2mTZvQ\n0FBz5wMAANXM2dm5b9++K1eu/P7778svbdu2rXG6QYMGubm5BQUFV65cycnJ+cc//mFc1Lt3\n75s3b167dk0I4ePjs3fvXiHE3r17n3322Weffbb0j5Vv2Qz7p2ayit2///3vhQsX+vv77927\n99KlS2fPnt28eXO3bt0++OCDmJgYc0cEAADVydraeseOHf/4xz/69u1b/kkElpaWZebodDq9\nXi+EkCSpzCLDnGefffb8+fNpaWnJyck+Pj4+Pj7JyclpaWkXLlzw8fGpfMsm2aPaQ1ax+/TT\nT2fMmLFu3bq+ffu6ubm1adPmn//85zfffNO7d+9PPvnE3BEBAEA1s7Gx2bRp07BhwwYOHLhz\n5877ru/q6uro6Fi6BSYnJ7u4uDRq1EgI8cQTTzg7O7///vs2NjadOnXq3LmztbX1+++/X69e\nvccff9yMu1H7yCp2v/3227hx48rMtLS0HD169G+//WaGVAAAQGGWlpbx8fETJ0588cUX73vb\nWkmSQkNDP/roo/j4+PPnz8fFxS1atMhwwpwQwsLCwtvbOyEhoW/fvpIkSZLUp0+fhISEPn36\nWFg81EOwUIasiycaNWqUl5dXfn5eXp6rq6upIwEAgBpBkqRFixY5OTmNGjUqNzd34sSJlaw8\ne/ZsrVYbERGRkZHh6uoaERExY8YM41IfH59t27YZD7z6+Phs2rSp9Al2MAnJcFC8cm+88cbJ\nkyd37txpZfVnEczMzOzZs+eUKVNmzZplzoQKk+LVvHdl6IOilI7wpxs3blhZWTk5OSkdREk5\nOTnFxcUuLi7lT1upPQoLC/Pz8x0cHKytrZXOohitVmt4NrednZ3SWZSUlZUlSZKzs7PSQaqV\naX8N1ajveZiDrBG7xx9/fO3atW3atBk3blzr1q3v3Llz6tSp+Ph4d3f3Vq1abd++3bjm0KFD\nzRYVAAAAlZFV7AICAgwTYWFhpeefOHHin//8Z+k5csb/AAAAYA6yit2mTZvMnQMAAAAPSVax\n8/PzM3cOAAAAPCSuMQYAAFCJCkfsEhIShBBjxoyxtLQ0TFckMDDQxKEAAADw4CosdkFBQUII\nf39/S0tLw3RFKHYAAAA1QYXFLikpSQhhuHGUYRoAAFQz7jyHB1JhsSv9UN7S0wAAAKiZ7n9V\nbEFBwdtvvz1s2LBu3bpVQyAAAGAkhfzHhFvTfzTYhFtDDXT/q2Lr1q0bHR1dXFxcDWkAAABQ\nZfcvdpIkNW/e/Nq1a9WQBgAAAFUm6z52Y8aMiYmJ0Wq15k4DAACAKpP15Il27dolJCQ89thj\nQUFBLVu2tLGxKb106NCh5skGAACAByCr2I0YMcIwERoaWn6pXq83ZSIAAABUiaxit2nTJnPn\nAAAAwEOSVez8/PzMnQMAAAAPSdbFEwAAoJYIDg5u1qyZYbqgoGD+/Pmenp5169Z1cXF58skn\n33nnHeNqkiSVucfthQsXLCwsJEk6evRodeeGEEJmsZs/f36HDh3KnEun0+nat28fERFhnmAA\nAEBhEydO/OSTT0JDQ7/77ruvvvpqypQpV69eNS5t2LBhSkrKjz/+aJwTFxfXvHlzJZLi/8kq\ndtu2bevfv78kSX95pYVFv379tm7dap5gAABASXq9ftu2bcHBwUFBQR07dnzqqafGjx+/dOlS\n4wpOTk5+fn6xsbGGH+/evZuQkDBhwgSF8kIImcXu0qVLHh4e5ee3bdv28uXLJk4EAABqAEmS\nmjRpcvDgwaysrIrWmTx58rp16/Ly8oQQW7du1Wg0AwcOrMaMKEtWsdPpdLdv3y4///bt2zxq\nDAAAtVq1atWlS5caNmzYpUuXV155Zdu2bTqdrvQKzzzzjKur67p164QQsbGxEyZMsLS0VCgs\nhJBZ7Nq2bbtr164yM/V6/a5duzw9Pc2QCgAAKO/pp58+f/78sWPHJk6cmJeXN3LkyN69e9+5\nc6f0OpMnT46Lizt37tw333wzfvx4paLCQFaxGz169IEDB2bOnGkYaxVC5OXlzZgxIzk5ecyY\nMeaMBwAAlCRJ0uOPP/7KK6+sWbPm66+//uabbzZu3Fh6hbFjx545c2bq1KnPPfdc06ZNlcoJ\nA1n3sZs+ffpXX30VExMTGxvr4eGh1+svXLhQWFjo6+v72muvmTsiAACoCVq3bi2EuH79eumZ\nTk5OL7300qpVq3bs2KFQLvxJ1ohdnTp1du3aFR0d3b59+7S0tMuXLz/22GMxMTE7d+6sU6eO\nuSMCAABFdO7cOSoqat++fT///POuXbtGjBhhbW09aNCgMqstW7bsjz/+eO655xQJidJkjdgJ\nIerUqTNz5syZM2eaNQ0AAKg5XnrppS+//HLBggU5OTn169d/6qmnDh482K5duzKraTQajUaj\nSEKUIZW57TDKkOJnKR2h+uiDopSO8KcbN25YWVk5OTkpHURJOTk5xcXFLi4uZe4iWasUFhbm\n5+c7ODhYW1srnUUxWq02Oztbo9HY2dkpnUVJWVlZkiQ5OzsrHaRaSSH/MeHW9B8NNuHWUAPJ\nOhSbm5ubkZFh/DEjI+ONN96YOHHiwYMHzRYMAAAAD0bWodipU6eePXv2+PHjQoiCgoIePXqk\np6cLIeLj4w8dOtSjRw/zZgQAAIAMskbsDh8+PHjw/w/ebty4MT09ff369ZcuXXJ3d1+4cKE5\n4wEAAEAuWcXu+vXrrq6uhumkpKT27dv7+/u7ubmNHz/eMIwHAAAAxckqdpIklZSUGKYPHz7c\nu3dvw3SDBg0yMzPNFQ0AAAAPQlaxa9GixYEDB4QQx48f/+233/r06WOYn5GRUa9ePfOFAwAA\ngHyyLp4YPXr03Llzr169mpqa6uLiMmDAAMP8EydOeHh4mDMeAAAA5JI1Yjd79uyQkJDLly83\nbtz4888/t7e3F0JkZWXt2LHDeFgWAAAAypI1YmdlZRUVFRUV9Ze719arV+/OnTvmSQUAAITg\nlsJ4QLJG7AAAAFDzyX1WrIFOp8vNzS3zFLJa/tAnAADM58w4Uz5RsP0qniOqcrKKnU6ni42N\nXbx4cVpa2t27d8ss5WmzAAAANYGsYvfuu++GhYV5eHgMGzbM0dHR3JkAAABQBbKK3fLly8eP\nHx8XF2dhwTl5AAAANZTcR4pNnDiRVgcAAFCTyepqzZs3z8nJMXcUAAAAPAxZxe7ll19evHgx\nF0kAAADUZLLOsfP09Pzss8+6d+8+atQoV1dXSfrLpddDhw41TzYAAAA8AFnFbvjw4UKItLS0\nY8eOlV/KSB4AAKoRHBy8aNEiPz+/TZs2GWe6ubn5+/t/8MEHhhU2b9585coVw6KCgoJ///vf\nmzZtOn/+vLW1tbu7+/PPPx8cHMxtbhUhq9iV/qcFAADqptFoNm/efOjQoWeeeabyNbOzs3v3\n7v3bb7+FhoZ269bNycnpzJkzK1assLOzmzVrVvWkRWmyip2fn5+5cwAAgBqiRYsWXl5eM2fO\nPH78eJnzr8qYPXv2+fPnU1JS3N3dDXM6d+48cuTIP/74o1qSoizuYAIAAMpasGDB6dOn16xZ\nU8k6xcXF69evHzdunLHVGdWvX9+c6VAhuc+K1ev1e/bs+e6777KysnQ6XelFMTExZggGAAAU\n07Jly+Dg4Llz5/r5+T3yyCP3XCc9PT0/P9/Ly6uas6ESsopdbm7uwIEDDx8+fM+lFDsAANRn\n7ty5CQkJCxYsCA8Pv+cKXD1ZA8k6FBsWFvbtt99GRkaeOXNGCLFjx47k5GRfX98nn3zy8uXL\n5g0IAACU4ODg8M477yxcuNB4AWwZrq6utra2p06dquZgqISsEbtt27a99NJLoaGhRUVFQggX\nF5fu3bs//fTT3bt3//jjjxcuXGjmkGbB3xnl1cDPpAZGqn61/EMw7L5er6/Nn4Nx32vzh2BU\nYz+Eyi8y+JsaP378J598Mnfu3Hsutba29vf3T0hICAkJKXOa3R9//MFpdoqQVewyMjIMFzwb\nHhdbXFwshLC0tPT391+yZMnftNjl5eVptVqlU9Qs2dnZSkf4C61WW9MiVbOSkhIhRC1/oJ/h\nt3h+fn5BQYHSWRRj+BDu3Llj+PqttQz9vsZ+Ldja2lpbWyudwsQsLCyio6N9fHwq2rUFCxYc\nPXr0ySefDA0NfeqppxwdHQ23O3nuuee43YkiZBU7W1tbw7eJtbW1RqO5evWqYb6Dg8Pvv/9u\nxnTmZG9vr3SEGsfZ2VnpCH+6ceOGlZVVLb+/ZU5OTnFxsZOTkypHAmQqLCzMz8+3s7NT369M\n+Qx/5NjY2NjZ2SmdRUlZWVmSJNWob6raoG/fvoMHD/7yyy/vubRevXrfffdddHT02rVrw8PD\nbWxsPDw8AgICJk6cWM05YSDJGdN+4oknnnrqqU8//VQI0b1798aNG2/durWkpGTAgAGXL1++\ncOGC+XMqRoqvRX9w6IOilI7wJ4qd+F+xc3Fxodg5ODhQ7DQaDcWuFha7M+NM+b9/+1U19EA2\nTEXWxRO+vr5btmwxDNpNmDBh+/bt7u7uHh4ee/fuDQoKMnNCAAAAyCKr2M2ZM2fv3r2G29dN\nmDAhKirK8IdjeHj4nDlzzJwQAAAAssg6x87R0dHR0dH4Y0hISEhIiNkiAQAAoCruP2JXUFAw\nZ86cY8eOVUMaAAAAVNn9i13dunWjo6Nr+TX2AAAANd/9i50kSc2bN7927Vo1pAEAAECVybp4\nYsyYMTExMdzOFwAAoCaTdfFEu3btEhISHnvssaCgoJYtW9rY2JReOnToUPNkAwCgtuPOc3gg\nsordiBEjDBOhoaHll9bYx/YBAADUKrKK3aZNm8ydAwAAAA9JVrHz8/Mzdw4AAFDetQRTPlKs\ncSAH2VRO1sUT3t7eJ0+eLD9/37593t7eJk4EAACAKpFV7JKTk7Ozs8vPz8zMTE5ONnUkAAAA\nVIWsYleR7OxsjUZjqigAAAB4GJWdY5eSkpKSkmKYTkpKunLlSumlWVlZS5YsadeunRnTAQAA\nQLbKit3WrVsjIiIM05GRkeVXqFu37oYNG8ySCwAAAA+osmI3cuTIrl27CiEGDx4cGRnp5eVl\nXCRJkr29fefOnR0cHMyeEQAAADJUVuw8PT09PT2FEGFhYQEBAW5ubtUUCgAAAA9O1sUT4eHh\ntDoAAGqD4OBgSZKGDx9eeqabm9ucOXNkriCEKCgoeO+99zp37mxra+vs7Pzkk0+Gh4ff8w4b\nMC1ZNygGAAC1h0aj2bx586FDh5555pkqrJCdnd27d+/ffvstNDS0W7duTk5OZ86cWbFihZ2d\n3axZs8ycvbaj2AEAgL9o0aKFl5fXzJkzjx8/Lkn3ePRF5SvMnj37/PnzKSkp7u7uhjmdO3ce\nOXLkH3/8Yfbotd5D3ccOAACo0oIFC06fPr1mzZoHXaG4uHj9+vXjxo0ztjqj+vXrmz4o/opi\nBwAAymrZsmVwcPDcuXMLCgoeaIX09PT8/PzSd9JAdaLYAQCAe5g7d65Wq12wYMEDraDX680f\nDRWSVey8vb1PnjxZfv6+ffu8vb1NnAgAANQADg4O77zzzsKFC8s8eqryFVxdXW1tbU+dOlVd\nMfEXsopdcnLyPS9RzszMTE5ONnUkAABQI4wfP97Dw2Pu3LnyV7C2tvb3909ISLhw4UKZlbl4\noho81FWx2dnZGo3GVFEAAECNYmFhER0d7ePjY21tLX+FBQsWHD169MknnwwNDX3qqaccHR0N\ntzt57rnnuN2JuVVW7FJSUlJSUgzTSUlJZUZis7KylixZ0q5dOzOmAwAAiurbt+/gwYO//PJL\n+SvUq1fvu+++i46OXrt2bXh4uI2NjYeHR0BAwMSJE6slcq0mVXKSY3h4eERERCUvrlu37oYN\nG4YMGWKGYDWFFF+L/rbQB0UpHeFPN27csLKycnJyUjqIknJycoqLi11cXO55H6laorCwMD8/\n38HBoaIBg9pAq9UajpDY2dkpnUVJWVlZkiQ5OzsrHaRaXUsw5f/+jQO5skHlKhuxGzlyZNeu\nXYUQgwcPjoyMLH3psiRJ9vb2nTt3dnBwMHtGAAAAyFBZsfP09PT09BRChIWFBQQE8LhYAACA\nmkzWxRPh4eFmjgEAAICHJet2J/Pnz+/QoUOZs/F0Ol379u0rPwkPAAAA1UZWsdu2bVv//v3L\nnL5tYWHRr1+/rVu3micYAAAAHoysYnfp0iUPD4/y89u2bXv58mUTJwIAAECVyCp2Op3u9u3b\n5effvn27uLjY1JEAAABQFbIunmjbtu2uXbtef/310jP1ev2uXbsMl80CAABz4M5zeCCyRuxG\njx594MCBmTNn5uXlGebk5eXNmDEjOTl5zJgx5owHAAAAuWSN2E2fPv2rr76KiYmJjY318PDQ\n6/UXLlwoLCz09fV97bXXzB0RAAAAcsgasatTp86uXbuio6Pbt2+flpZ2+fLlxx57LCYmZufO\nnXXq1DF3RAAAai3JpJTeG5idrBE7IUSdOnVmzpw5c+ZMs6YBAABAlckasQMAAEDNJ3fEzkCn\n0+Xm5pZ5BIWTk5NJIwEAAKAqZBU7nU4XGxu7ePHitLS0u3fvlllapucBAABAEbKK3bvvvhsW\nFubh4TFs2DBHR0dzZwIAAEAVyCp2y5cvHz9+fFxcnIUF5+QBAADUULKK2vXr1ydOnEirAwAA\n8gUGBg4aNEjpFLWLrBG75s2b5+TkmDsKAACoIXJzc5s2barVatPT011cXOS8ZPTo0Xl5edu3\nbzfO6dWrl/GZVagesgbhXn755cWLF3ORBAAAtURiYmKnTp169uwZHx9f5Y1MmDAhODjYhKlw\nX7KKnaenZ2pqavfu3RcvXrxt27btf2XuiAAAoJotW7YsKCgoKCgoNja29MiOn5/fiBEjwsPD\nmzdv7uTkNGjQoKtXrwohJkyYsHbt2i+++MLwiItly5aJvx6K/eabb3r16uXg4GBnZ9ehQ4et\nW7cat7lq1aqOHTtqNBo3N7cZM2bk5uZW/l5V21otIetQ7PDhw4UQaWlpx44dK7+UkTwAANTk\n8OHDFy9efOmllywtLadNm7Znz55+/foZl27fvt3Ly+vcuXNFRUXPP//89OnTt2zZsmLFiqKi\nojKHYo2Ki4sHDRo0ceLE1atXW1hYnDlzxsbGxrDo448/joiIiImJ6dGjx7Vr14KDg4OCgjZv\n3lzJe1V5a7WBrGK3adMmc+cAAAA1xNKlS4cPH25nZyeE8Pf3X7ZsWeli5+np+eabbwohNBrN\ntGnTpkyZct8NZmdn5+TkDBw4sFWrVkIINzc3w/ySkpLw8PAPP/xw1KhRQohWrVrFxcV16dLl\n+vXrDRs2rOi9qry12kBWsfPz8zN3DgAAUBPcuHFj8+bNSUlJhh9ffvnlp59++urVq02aNDHM\nadu2rXHlBg0a5ObmFhQUPPLII5Vss379+oGBgQMGDPD29u7du/eQIUO8vLyEEJcuXbp58+aE\nCRMmTJhQev0LFy4Yqtg936vKW6sNuIMJAAD4U3x8/J07d/r06WNlZWVlZdWzZ0+tVrtixQrj\nCpaWlmVeotPp5Gz2xIkTvr6+33zzzeOPP75gwQIhRElJiRBi165d+r96+umnK3+vqm2tNqhw\nxC4hIUEIMWbMGEtLS8N0RQIDA00cCgAAKEGv18fGxs6cOfPll182zly/fv3y5cvnzZtXvmaV\nZm1trdVqK1nBy8vLy8tr1qxZc+fOjY2Nff3111u3bu3s7Lx9+/YBAwY8aFTTbk01Kix2QUFB\nQgh/f39LS0vDdEUodgAAqMPXX3998eLFyZMnt2nTxjhzypQp77///o4dO1544YVKXtu6des9\ne/akpqbWr1/fzs5Oo9EYF/3yyy8JCQlDhgxp3rx5ZmbmgQMHOnToIISwsrKKiIgIDg5+9NFH\nR4wYYWNjk5qaumXLltWrV1fyRqbdmspUWOwMB9etra2N0wAAQN2WLl3auXPn0q1OCOHq6tqz\nZ89ly5ZVXuwmT578zTffdO/e/fbt20uXLi19UYWtrW1qauqaNWv++OMPFxcXX1/fqKgow6Lp\n06fXr18/Ojr6o48+qlOnTqtWrYYOHVp5SNNuTWUkblZSOSl+ltIRqo8+KErpCH+6ceOGlZWV\nk5OT0kGUlJOTU1xc7OLiIkmS0lkUU1hYmJ+f7+DgYPg7s3bSarXZ2dkajcZwlWKtlZWVJUmS\ns7Oz0kGqlWn/9+eXvurJunjC29v75MmT5efv27fP29vbxIkAAABQJbKKXXJycnZ2dvn5mZmZ\nycnJpo4EAACAqnio250YDg2YKgoAAAAeRmU3KE5JSUlJSTFMJyUlXblypfTSrKysJUuWtGvX\nzozpAAAAIFtlxW7r1q0RERGG6cjIyPIr1K1bd8OGDWbJBQAAgAdUWbEbOXJk165dhRCDBw+O\njIw0PK/DQJIke3v7zp07Ozg4mD0jAAAAZKis2Hl6enp6egohwsLCAgICjA/ZBQAAQA1UWbEz\nCg8PN3MMAABwD9x5Dg9E1lWxubm5GRkZxh8zMjLeeOONiRMnHjx40GzBAAAA8GBkjdhNnTr1\n7Nmzx48fF0IUFBT06NEjPT1dCBEfH3/o0KEePXqYNyMAALXVibmmvEjxiUh/E24NNZCsEbvD\nhw8PHjzYML1x48b09PT169dfunTJ3d194cKF5owHAAAAuWQVu+vXr7u6uhqmk5KS2rdv7+/v\n7+bmNn78eMMwHgAAABQnq9hJklRSUmKYPnz4cO/evQ3TDRo0yMzMNFc0AAAAPAhZxa5FixYH\nDhwQQhw/fvy3337r06ePYX5GRka9evXMFw4AAADyybp4YvTo0XPnzr169WpqaqqLi8uAAQMM\n80+cOOHh4WHOeAAAAJBL1ojd7NmzQ0JCLl++3Lhx488//9ze3l4IkZWVtWPHDuNhWQAAAChL\n1oidlZVVVFRUVFRU6Zn16tW7c+eOeVIBAADggckasQMAAKjJAgMDBw0apHQK5VU4YpeQkCCE\nGDNmjKWlpWG6IoGBgSYOBQAAFBIcHLx58+YrV65Uvtro0aPz8vK2b99ePanuq1evXnl5eUqn\nUF6FxS4oKEgI4e/vb2lpaZiuCMUOAADIVFxcXKdOHZNvdsKECSbf5t9RhcUuKSlJCGFtbW2c\nBgAAtYqfn5+lpWW7du1Wrlx5+/btXr16xcXFNWnSZMKECWvXrhVCSJIkhFi6dOmUKVOEEKtW\nrfroo4/OnTvXqFGjF1544d133zVccOnn5ydJUqtWrRITE69fv56Xl2dpafnGG2+sXr26qKho\n4MCB3t7eM2bM0Gq1hvetZDv3zCOECAwMvHHjxo4dOwxbWLNmTXR0dGpqqq2tbdeuXRMTE+vX\nr6/ER1jdKix2Pj4+95wGAAC1x/bt2728vM6dO1dUVPT8889Pnz59y5YtK1asKCoqKnMo9uOP\nP46IiIiJienRo8e1a9eCg4ODgoI2b95sWPrFF1+EhIT88ssvQggbG5uwsLD4+Pi4uLiuXbvu\n3Llz/vz5MrdzzzxlMi9evHjWrFkRERFDhw7VarUHDhwwVkbVk3VVLAAAqJ08PT3ffPNNIYRG\no5k2bZphZK68kpKS8PDwDz/8cNSoUUKIVq1axcXFdenS5fr16w0bNhRCuLq6vvfeexYWFkII\nnU4XExPz1ltvDR8+XAjx6quvHj9+3DAEeN/t3DePYQuvvfZaaGioYY6Xl5dZPpoaiatiAQBA\nhdq2bWucbtCgQW5ubkFBQfnVLl26dPPmzQkTJkj/06VLFyHEhQsXDCt06NDB0OqEEOnp6bm5\nuT179jS+3Dh93+3cN09aWtqtW7d8fX0fetf/liocsdNoNDI3UVRUZKIwAACgZrG0tCwzR6fT\nlV/N8Ez5Xbt2GR9PVUbdunXLzDGcn2eg1+tlbue+eQybKr3xWqXCYlfmZjBnzpxJTU1t2rRp\nmzZtJEk6e/ZsRkZGu3bt2rdvb/6QAACgZrG2ti594lrr1q2dnZ23b99eUSErzdXV1d7e/vDh\nw8aBum+//bYK27knwxZ2797dr1+/qm3hb63CYmc8S1EIceTIkf79+69cuXLcuHHGo+MrV64M\nDg5esWJFdcQEAAA1SevWrffs2ZOamlq/fn07OzuNRhMREREcHPzoo4+OGDHCxsYmNTV1y5Yt\nq1evLv9aCwuL4ODgyMjIli1bPvHEEzt37tyxY4dhjM3Kykr+du7J0tJy/vz5r7/+uouLy9Ch\nQ/V6/YEDB1588cXGjRubcv9rKlnn2M2ZM2fcuHFBQUHGo+MWFhYTJkwYO3as8cxEAABQe0ye\nPPmxxx7r3r17/fr1DQ8ymD59+tq1a7/++utu3bp17do1LCysVatWFb38SyQnRQAAIABJREFU\nrbfeCgoKmjx5speXV3Jy8pw5c4zngD3Qdu4pODh46dKl69at69Sp0zPPPPOf//zHyqq2XC0q\nGY9qV8LW1jYmJmbixIll5i9fvjw4ODg/P9882WoEKX6W0hGqjz4o6v4rVZcbN25YWVk5OTkp\nHURJOTk5xcXFLi4utfZkESFEYWFhfn6+g4OD4baatZNWq83OztZoNHZ2dkpnUVJWVpYkSc7O\nzkoHqVYn5m4w4daeiPQ34dZM6NVXX/32229PnDihdJC/PVkF1tra+scffyw//8SJEzY2NqaO\nBAAAVO7ixYuHDh3q3bu3hYXFV199tXz58iVLligdSg1kHYodNGhQbGxsXFyc8TRJrVa7bNmy\n5cuXDx482JzxAACACkmStGzZsk6dOrVv3z42NnbZsmWTJk1SOpQayBqxW7hw4XfffTd58uR5\n8+Z5eHjo9frz58/fvHmzTZs2CxYsMHdEAACgMq1atTp69KjSKVRI1ohdo0aNTpw4ER4e3qRJ\nk5SUlFOnTjVt2jQiIuL777833AYaAAAAipN7kYi9vX1YWFhYWJhZ0wAAAKDKeKQYAACASlDs\nAAAAVKK23K8PAIC/oxp75znUTIzYAQAAqATFDgAAQCUqLHZdu3bdt2+fYToxMfH69evVFQkA\nAABVUWGxO3HiRFZWlmF6zJgxqamp1RUJAAAAVVFhsWvUqNGFCxeqMwoAAAAeRoVXxfbr12/+\n/Pn79+93dnYWQrz99tvLli2755obNmwwVzoAAADIVmGxi46OliQpKSnp999/F0Ls37+/ojUp\ndgAAADVBhYdiH3300VWrVl29elWn0wkh9u/fr69ANaYFAABAhWTd7mTatGlNmzY1dxQAAAA8\nDFlPnvj4448NE7dv3758+bIQws3NzcHBwXyxAAAA8KDk3qD47Nmz/fv3d3Z27tSpU6dOnZyd\nnQcMGPDLL7+YNRwAAADkkzVid+HChZ49e966datHjx5eXl5CiNOnT+/evbtHjx7Hjh1zd3c3\nc0gAAADcn6xiN3/+/IKCgt27d/v6+hpnfv3110OGDAkLC1u7dq3Z4gEAAEAuWYdi9+zZM3Xq\n1NKtTgjh6+v7yiuv7NmzxzzBAAAA8GBkFbvs7GwPD4/y8z08PLKzs00dCQAAAFUh61BskyZN\njhw58sorr5SZf+TIkSZNmlTywnPnzm3ZsuXixYuZmZn9+vWbPn166aXff//9mjVrrly54ujo\n6OPjExAQIEmSmRYBAAConqwRu2HDhiUmJn7wwQdFRUWGOUVFRZGRkWvXrh02bFglLywqKmrc\nuPGYMWMaN25cZtEvv/zy7rvvtm/fPjo6evTo0Vu3bjWeq2fyRQAAALWB3IsnkpKSQkND33vv\nPXd3d71ef/Hixby8PC8vr7feequSF3bs2LFjx47i/9q787ioyv7/49cZlhlkVxCXFAFxxSiX\ncE1v81a5U1KT3AntLm27TdPbb1mJZT6s1DSy9NbClUhx777dQ+uB+FBvE0XRVNQHLqSJICjr\nzPz+OI/mNzfiiASc4eL1/GvOOdec87kuGXxzrnPOCLFp06YymzZt2tS0adOJEycKIfz9/a9f\nv75169bIyEi9Xl/lmyo1MgAAALVMhc7YeXl5HTp0KCYmJjAw8Ny5cxcuXAgMDJw9e3ZKSoqX\nl1flDpyent6xY0fLYseOHQsLCzMyMqpjEwAAQF1QoTN2QghXV9dZs2bNmjWrSo5qNptzcnK8\nvb0ta9TX2dnZVb7pQTXcuXOnpKSkSrojjVu3bmldwv8oLS21t5JqmPpdzDZ+jOuOvLw8rUvQ\nXmFhYVFRkdZVaEn9RNjtrwV3d3dnZ2etq0BdV9FgJx+dTufg4KB1FfbFrgaktLRU2FlJNc9o\nNJrN5jo+CCaTyWw263S6unwvlNlsNhqNiqLU8R8Go9Eo7PjXQl3+EYX90CbYKYri5eV1+/Zt\nyxr1df369at804NqcHNzq8ouScHGxLry9vaarERD5gWDtS7h/8vNzS0pKfH09KzL/2EUFBTc\nvXvX1dW1Lp8LKS0tzcnJ0ev1dfwXV3Z2tvrbXutCAPtV0e+KrXJt27Y9duyYZfHYsWMGgyEw\nMLA6NgEAANQF1RvsiouLMzIyMjIyiouL8/PzMzIyLl68qG4aNmzY1atXly1bdvny5aSkpM2b\nN0dERKh3sFb5JgAAgLpAUa9FrSYZGRlvvfWW9RqdTrdlyxb19ZEjR9auXZuZmak+T3j06NGW\n+aYq31RpSty0P7mHWsQ8fv6DNjEVqwl1KrZBgwZMxXp4eDAVazAYmIpVFMX6PjkAZVRvsJMA\nwU5FsNMEwU4Q7IQQBLs/EOyAh9LsGjsAAABUrQoFu1u3brm5ue3fv7+aiwEAAEDl2Qp26hOD\nhBBms/nu3bvqc8Xu3Lnj4+OTnJxcE9UBAACgwmw9x87b27tbt249e/Zs3769ZaXJZLp16xbf\n2QAAAGBvbAW7mTNnHjhwYMGCBbm5ueri0aNHw8LCaqo2AAAAPAJbU7EzZsz4z3/+k52drV5d\nZzQa586d27dvXyHE3LlzlyxZcurUKW6qBQAAsBMPv3lCp9OpU7Hz5s27devWzp07xR8PqAsJ\nCfHz86v2GgEAAFABtqZiGzVq1KtXr169elmusXNyclKnYlesWNGlS5fk5OQDBw7URJkAAAB4\nGFvBLioqKjk5+Z///GdRUZEQYs6cOWfPnu3atau61dXVtX///v3796+JMgEAAPAwtoLdp59+\nKoQoKipKSkoKDw/PysqaMmWKej/s4sWLr1y50rt372bNmtVQpQAAALDp4dfY6fX6zp07CyG+\n/PLL33//PT4+Xghx+PDhqKio5s2bBwYGVnuNAAAAqABbZ+zu5+HhER4eLoRYt25d+/bt9+/f\nz9dRAAAA2IkKBTsXF5fJkyc3b97ceqWvr29kZGRkZGT1FAYAAIBHU6Fg5+rqumjRIvW1Xq+f\nOHFikyZNqrMqAAAAPLJHm4oVQri4uCxdurQ6SgEAAMCf8fCbJwAAAFArEOwAAAAkQbADAACQ\nBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4A\nAEASBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIE\nOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAA\nSRDsAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwA\nAAAkQbADAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRB\nsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAA\nkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAACQBMEO\nAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEAS\nBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAA\nAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAASRDs\nAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASjloXoJk7d+6UlJRo\nXYV9uXXrltYlaM+uBsFsNgshsrOztS5Ee3l5eVqXoL3CwsKioiKtq9CS+omwqw+pNXd3d2dn\nZ62rQF1Xd4Odh4eH1iXYnQYNGmhdgvbsahByc3NLSkrq16+vKIrWtWimoKDg7t27dfy/zNLS\n0pycHIPB4ObmpnUtWsrOzlYUxdvbW+tCAPvFVCwAAIAkCHYAAACSINgBAABIgmAHAAAgCYId\nAACAJAh2AAAAkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEiCYAcAACAJgh0AAIAk\nCHYAAACSINgBAABIgmAHAAAgCYIdAACAJAh2AAAAkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAA\nAJIg2AEAAEiCYAcAACAJgh0AAIAkCHYAAACSINgBAABIwlHrAoBa4/SLiibHvVHjR2y3ylzj\nxwQAVAHO2AEAAEiCYAcAACAJgh0AAIAkCHYAAACSINgBAABIgmAHAAAgCYIdAACAJAh2AAAA\nkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEiCYAcAACAJgh0AAIAkCHYAAACSINgB\nAABIgmAHAAAgCYIdAACAJAh2AAAAkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEiC\nYAcAACAJgh0AAIAkCHYAAACSINgBAABIgmAHAAAgCYIdAACAJAh2AAAAkiDYAQAASIJgBwAA\nIAmCHQAAgCQIdgAAAJIg2AEAAEiCYAcAACAJgh0AAIAkCHYAAACSINgBAABIgmAHAAAgCYId\nAACAJAh2AAAAkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEjCUesCANQm11cqmhz3\nlhYHbRxt1uKwAFB5nLEDAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAA\nAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAACQBMEOAABAEgQ7AAAASRDs\nAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAk\nQbADAACQBMEOAABAEgQ7AAAASRDsAAAAJEGwAwAAkATBDgAAQBIEOwAAAEkQ7AAAACRBsAMA\nAJCEo9YFVKWjR4+uWbPmypUrnp6e/fr1GzVqlKIoWhcFQDZ16heL2WzWugQAj0CeM3Znz56d\nM2dOu3btFi5cOHbs2E2bNq1bt07rogAAAGqOPGfsNm3a1LRp04kTJwoh/P39r1+/vnXr1sjI\nSL1er3VpACCh/76boHUJNaTT3JFalwBUlDxn7NLT0zt27GhZ7NixY2FhYUZGhoYlAQAA1CRJ\nztiZzeacnBxvb2/LGvV1dnb2g96Sl5dXWlpaE8XVHrdv39a6BO0xCIJB+APjIBgEIUSFB8HV\n1dXZ2bm6iwFskyTYVZ/sITNr+Ihms9lkMul0Oru6QDv7vZ41fESj0SiEcHBwqOHj2tB40QP/\nTqgmdjgIhudqehDs8BNh4y/GaqIOgqIoOp0dTbMETh9Qw0c0mUxCCLsaBMDeSBLsFEXx8vKy\n/qNKfV2/fv0HvcXd3b0mKnt0xcXFd+7cMRgM9erV07oWLf3++++Ojo5eXl5aF6Kl3NzckpIS\nLy8v+8k0Na+goODu3btubm51+VxIaWlpTk6OXq93c3PTuhYtZWdnK4piPTkDoAx5/u5p27bt\nsWPHLIvHjh0zGAyBgYEalgQAAFCT5Al2w4YNu3r16rJlyy5fvpyUlLR58+aIiAhuiQUAAHWH\nJFOxQojWrVvPnDlz7dq1u3bt8vT0HDp06OjRo7UuCgAAoObIE+yEEF26dOnSpYvWVQAAAGhD\nnqlYAACAOo5gBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEiCYAcAACAJgh0AAIAkCHYAAACSINgB\nAABIgmAHAAAgCYIdAACAJAh2AAAAkiDYAQAASIJgBwAAIAmCHQAAgCQIdgAAAJIg2AEAAEiC\nYAcAACAJgh0AAIAkCHYAAACSINgBAABIgmAHAAAgCYIdAACAJAh2AAAAkiDYAQAASIJgBwAA\nIAmCHQAAgCQIdgAAAJIg2AEAAEhCMZvNWteA/2EymYqLix0dHR0dHbWuRUuFhYU6nc7Z2Vnr\nQrRUXFxsMpkMBoPWhWiptLS0tLTUycnJwcFB61o0Yzabi4qKHBwcnJyctK5FS0VFRUIIvV6v\ndSGA/SLYAQAASIKpWAAAAEkQ7AAAACRBsAMAAJAEwQ4AAEASBDsAAABJEOwAAAAkQbADAHux\naNGiDz/8UOsq6gpGG1Kq04/ArQ7Lly8/ePBgXFycupibmzt79uyCgoKYmBg/Pz9ta7Mr1gNV\nVFSUmJj4888///7773q93s/PLywsbMSIEWqz7du3BwcHL1iwwPLe69evT5o0yWw2f/bZZ61b\nt9asD1VB7WCPHj1mzJhhWfn3v/+9V69eL774orjvJ6qoqGjr1q3JycnXrl1zdHRs3Lhxly5d\nIiIiXF1dtelAZVWk47YbCDsbjYKCgujoaKPRGBcX5+7uXpG3LFy4sKCgYObMmZY17dq1Kyws\nrLYaq0uZn9IHub+/2qqlow3YRrCrRlevXo2JiWnQoMHs2bMr+Iu+bvryyy+PHTs2fvz4li1b\nFhcXX758+fz585atXl5ely5dysjICAwMVNfs2rXL19f3xo0bGtVbxZydnZOTk0+dOtW+fXvb\nLe/evfvOO+/cvHkzMjIyODjY1dU1MzNz9+7dBoNh6NChNVNtFXpox203sLfR2L9/f0BAgJOT\n0969eytdQP/+/au2KjkYjcbq+N4RRhtSIthVlzNnznz00UehoaFTpkyxfAvQvHnzdDpds2bN\n9uzZc+/evXbt2r3xxhv169cXQpjN5g0bNuzevfvWrVs+Pj7h4eFDhw5VFCU1NfWDDz5Yu3at\nGg3HjRunKMrq1auFEHl5eWPHjp0zZ06HDh1s7NnOmc3mlJSUyMjIfv36qWtatWr117/+1dLA\n1dX1iSee2Llz52uvvSaEKC0t3bdv36BBg9atW6dNxVXN19e3RYsW33zzzYIFCxRFsdEyLi7u\n2rVrsbGxjRs3VtcEBgb27t07Nze3RiqtYg/tuO0G9jYaO3bsiIiIcHR0/O6774YMGWIp+EGf\nzdjY2P379wshIiIihBCvvvpqeHj4okWL7ty588EHHwghTp8+vWrVqkuXLpnN5oYNG44ZM6Zb\nt27qPn/88cfNmzdfu3bN29s7LCxs7NixLi4uNo5Vub1VziP113ZfFEXx8/Pbv39/Tk7O+vXr\nHRwcVq5c+eOPP5aUlHTq1CkkJGT58uVbtmyp9JhYj7YQIikpaevWrZmZmQaDoWXLllOnTvX0\n9Kz0OABa4Rq7anHw4MH33nvvmWeemT59epnvdjx06JBOp1u6dOk333yTn5+/bNkydf3GjRs3\nbNgwcuTIJUuWDB8+PD4+fvv27UKIdu3aOTk5nThxQghx+fLlwsLCgoKCy5cvCyFOnDjh5OTU\npk0b23u2c4qi1K9fPy0tLS8v70FtBgwYcODAAXXSJCUlxcnJqVOnTjVYY7WLjo6+fPlyUlKS\njTZGo/Gnn3565plnLDnGovb+9/PQjj+ogb2NRnp6elZWVs+ePbt3756bm5uammq9tdzP5ptv\nvtmnT5+wsLBt27Zt27ZNTTkWRqPxww8/bNOmzeLFi7/88svo6Oh69eqpm/79739/++23zz//\n/JIlS95+++309PTFixfbPlal91Y5Fe9vBfuSkJDg5OSUkJCwd+/e1157LTY2NiQkJD4+3tKy\nEmNSxvbt27/44ouePXsuWrRo7ty5Xbp0MZlMf3IcAE0Q7KpeTk7OJ598MmjQoAkTJtx/mqFp\n06YjRoxwdnZ2dXV99tlnjx8/LoQwm82JiYlDhgzp169fkyZNBgwYMHjw4O+//14I4eTk1LZt\nWzXYpaamtm/fvn379pZFNfbZ2HOtMGXKlN9++y0qKuqtt976+uuvU1JSynyFcfv27X18fA4c\nOCCE2LlzZ//+/XU6qX50/fz8IiIi1qxZo37Heblu3rxZWFjo7+9fk4VVt4d2/EEN7G00duzY\n0aNHD4PB4Ozs/PTTT+/YscN6ayU+m/n5+ffu3evUqVOjRo0aNmzYuXPn0NBQIYTJZIqPj3/x\nxRf79OnTqFGjtm3bvvHGGwcPHszJybFxrErvrXIq2N+HHt3X13fcuHEGg8FgMAghtm3bFhkZ\n2aNHDz8/v2effbZLly4V3M9D6zGZTN99993gwYOHDx/erFkzf3//QYMGeXt7/5lBALTCVGzV\nc3Nza9GixZ49e3r27NmyZcsyWx977DHLay8vr4KCgqKiory8vHv37oWEhFg2hYSEbNy48fbt\n297e3qGhoXv37hVCpKamhoaGms3m1NTUwYMHp6amWl8jUu6e9Xp9dfWz6rRt23bZsmUXL148\ne/Zsenr6ggULWrZs+dFHH1mf7Bw4cOCuXbtCQkJOnz49derUWjr5aENkZOS+ffs2bdo0atQo\nrWupUQ/tuP2PzJ07d5KTky33V/br12/GjBnZ2dmWayEq8dn09PR85plnZs2a1aFDh5CQkLCw\nMDXF/vbbb3l5ebGxsbGxsdbtr1+/7uXl9aBjVXpvlVPB/j706M2bN7f8bXzz5s2CgoK2bdta\nmrVp00ad263cmFjXk5WVlZ+f/+STT1a6y4D9INhVPUdHxw8++GD+/PkzZ858//33reOaEOL+\nU01ms7nMCaoyQkNDV69enZWVlZaWNmbMGCHE+vXrs7Kyrl+/rv7ZbWPPf6onNUhRlMDAwMDA\nwPDw8AEDBrzzzjs///xz3759LQ369u27atWqpUuXdurUqUGDBvIFu3r16o0ZM2bFihXW1xda\n8/HxMRgM6iy8TB7a8XIb2NVo7N27t6SkxPpmT5PJtHv37pEjR6qLlftsTp48eciQIceOHUtN\nTY2Pj4+Kiho2bJg6PxgTE9OxY8dy3/WgY1Vub5VTwf4+9Oj3Z0HrORDLPis9JmUWbV/hCtQW\nUs1n2Q8nJ6cZM2Z069YtJibm6NGjD23v4+NTr169tLQ0y5q0tDR3d3d1LqBly5Zubm6JiYlO\nTk4BAQEBAQGOjo6JiYnu7u5BQUHV2A2NqFdNlZkMcnV17dmzZ2pq6sCBAzWqq9r179+/SZMm\na9asKXero6Njr1699u3bd/369TKbanvMtd3xchvYz2iYzeadO3c+99xzi61ERkbu3r37oRdp\nOTo6Go1GGw38/f2HDh0aExMzbNiwnTt3CiEaNWrk5uZ26NChSpRatXurhDL9faSj+/r6uri4\nnD592rLm7NmzldhPuRo3buzm5nbs2LFK7wGwH5yxqy46nW7y5Mmurq4ff/zx1KlTe/XqZaOx\noiiRkZEJCQl+fn7t2rU7efLk9u3bo6KiLFs7dOiwb9++bt26qX9TqothYWFy/Ik5efLkPn36\nBAUFeXl53bx5c8OGDY6OjpYLaCxee+21CRMmSPzgGEVRXnrppffff9/RsfwPZnR09NmzZ6dO\nnRoZGdmqVat69eqpD/jo3LlzbXzcicVDO15uAzsZjV9++SUrK2vgwIFNmza1rAwPD09MTDxy\n5EhYWJiN9zZq1Oj48eOZmZmenp7q9XmWTVevXt23b99TTz3l6+ubm5ublpbWvHlzIYSDg8Po\n0aOXL1/u4eHRq1cvR0fHK1euHDx4cMqUKTYOVLV7q7T7+1vxoyuKEhERsWHDBj8/v5YtWx49\nevTIkSPqL8A/3wudTjdy5Mi4uDgPDw/1n+zkyZPdunXjMjvURgS7aqQoyssvv+zq6jp//vx7\n9+4NGDDARuNhw4YZjcaEhAT1cSejR48ePHiwZWtoaGhKSopl4vWJJ55ITk62noet1Xr27Hn4\n8OFNmzbdvXvX09OzVatW8+bNa9asWZlmzs7O1v/zSenxxx/v0qXL4cOHy93q7u4+f/78rVu3\n7t+/Pz4+3snJqUmTJk8//bTtH61awXbHy21gJ6OxY8eOwMBA61QnhPDx8WnTps2OHTtsB7uB\nAweePn16+vTp9+7dszz+Q6XX6zMzM5OSknJzc93d3Z988skJEyaomwYNGuTp6blly5YtW7Y4\nODg0atSoa9eutous2r1V2v39faSjjxw5srCw8KuvvlIfd/L8888nJCRUVS8iIiJcXFy2bdu2\nbt26evXqBQcH9+jR40/1FtCIUosuwwIAQLVs2bIzZ858/vnnWhcC2BfO2AEAaoGsrKxTp06F\nhIQoivLf//53165dEydO1LoowO5wxg4AUAtkZWXNnz8/MzPTZDI1adJk8ODBlq+rAWBBsAMA\nAJAEjzsBAACQBMEOAABAEgQ7AAAASRDsgLpl7969iqKsXLlS60KqXd3pKQBYEOyAivr1119j\nYmJOnDihdSE1TY6Oy9ELALCNYAdU1K+//jp79uzangz69u1bUFAwbty4ir9Fjo7L0QsAsI1g\nB1Sxe/fuaV1C+dTCdDqdwWBwcHDQuhwAQNUj2AEVEhMTo35777hx4xRFURSlT58+QojExERF\nUdavXz979uzg4GBnZ+cPP/xQCJGbm/vee++FhYX5+Pjo9frAwMBp06bl5+dbdqi+cePGjZ98\n8kmrVq30en3z5s0//vhj60dLlpaWfvrppx06dHB3d3d3dw8ODo6Ojs7Ly7NusGjRok6dOrm6\nurq7uz/++OOzZs2y3v/9hZW58kxtlpCQMHPmzBYtWuj1+uDg4EWLFj204zZUpGtCiJycnLff\nfjsgIECv1/v5+Y0ZM+b8+fPqpj179iiK8vHHH5fZc3R0tIODQ2ZmpqX7CxcufOKJJ1xcXNzd\n3fv06bN79+6K//MlJSUpivLRRx+VaTxu3Dj1i+RtdxMA7JEZQAVcvHhx7ty5Qoh33303KSkp\nKSnpl19+MZvNGzZsEEK0aNGiR48e69ev/+mnn1JSUsxm88mTJ319fV999dXPP/98yZIlI0aM\nUBSlV69eJpNJ3aH6xsDAwIEDB+7YsSMlJeWll14SQnz99deWg06bNk0IMXr06H/9618rVqx4\n//33n3zyyStXrqhbS0pK1C+8792797x587766qt//OMfbdu2td7SLhkAAAAHr0lEQVT//YXt\n2bNHCBEXF2fdrGnTphEREUeOHElPT58xY4YQ4v/+7/9sd9yGinQtPz+/Q4cOQogxY8YsWbLk\nrbfe0uv13t7eZ86cMZvNRqOxWbNmwcHB1rvNz893c3Pr16+fulhaWhoeHq7T6UaMGBEbG/vZ\nZ5+FhoYqihIfH682sO7pg3rRunVrf39/o9FoOcrt27ddXFwGDRpUgR8KALA7BDugorZv3y6E\nWLNmjfVKNcS0atWqpKTEen1hYWFxcbH1GvX80549e6zf2LlzZ0vUMxqNwcHBlmRmNpsDAgL+\n8pe/PKge9evP33zzTcse1J3YLqzcYBcQEGDdbOTIkTqd7ty5czY6bkNFujZ79mwhhHoaT7Vr\n1y4hxIABA9TFd999VwiRnJxsaRAXFyeEWLt2rbq4ZMkSIcS3335raVBcXNyxY0c/Pz+1L2V6\nWm4vFixYIITYtWuXZU1sbKwQYtu2bRXsLADYFaZigSowfvx4R0dH6zV6vd7JyUl9XVJSUlhY\nOHToUCHEoUOHrJupM4Pqa51O17lz5wsXLphMJnWNl5dXenr6kSNHyj3o2rVrXVxc5s6da9mD\nuhPbhZUrOjrautnLL79sMpm2bNny0DfaYLtrGzdudHNzmzp1qqV9//79u3XrtmfPnjt37qgl\nCSFWrVplabBy5UoPD49hw4api6tXr27YsOGoUaMK/2A0GkeNGvXbb7+lpqZWsMjo6GiDwbB8\n+XLLmuXLlz/22GN/+9vfKt1xANAQwQ6oAgEBAfevXLlyZffu3V1dXZ2dnV1cXNq1ayeEyM7O\ntm7TrFkz60UPD4/i4mLLVXTz588vKSl56qmn/P39x4wZExcXZ31nxq+//tqyZUs3N7dHLex+\nQUFB1ouBgYFCiAsXLlTkvQ9iu2sZGRlBQUEGg8G6TYcOHUwm06VLl4QQwcHB3bt3//777wsL\nC4UQly5d+umnn1544QUXFxe1cXp6+o0bN1z+1/Tp04UQN27cqGCR9evXHz58+LZt227evCmE\nOHz48IkTJyZMmMDNJQBqKYIdUAX0en2ZNQsXLhw/fryPj8+KFSv279+fkpLyww8/CCEsp6xU\n1ifbLMx/3GTQt2/fixcvrl+//tlnnz1+/PiECRPatGlz9epVS7Ny3267sHIVFRXdv/jQndtm\nu2sVKT46Ojo3N1c9cbhq1Sqz2ayexlOZTKbg4OCU8nTr1q3idU6aNKm4uHj16tVCiOXLl+t0\nOvWKQACojR4+RwNA9UhB55tvvgkICNi6davlXT///HMlDuru7h4ZGRkZGSmESEhIGDVq1Bdf\nfPHJJ58IIVq3bn369Gn1loJK7NlaWlra/YvqeTvxpxNeuYKCgs6fP19YWGh90i4tLU2n07Vo\n0UJdHDFixOTJk1etWjVixIjVq1cHBwf36NHD0rhVq1ZpaWkhISEV7P6DetGjR4+QkJAVK1a8\n8sorCQkJAwYMaN68eeU7BgCa4owdUFHu7u7ivrnUB9HpdGaz2Wg0qotGo1G9K/ORlDlW165d\nrVeOHTu2oKDg/ffft25j/t9HilRQXFxcVlaW+rqkpGTBggWKojz33HPqmkfqeAUNGzYsPz/f\n+rkqe/fuPXjwYL9+/Tw8PNQ1Hh4eQ4cO3bNnz/fff5+RkfHiiy9a7yEqKqq4uHjatGllunzt\n2rVyj2ijFxMnTjxz5swbb7yRn5//yiuv/MmuAYCGOGMHVFRoaKjBYIiNjXV2dvby8mrYsGHf\nvn0f1Hj48OExMTHh4eEvvPBCXl5eQkJCJSJXkyZNBg0a1KlTp6ZNm964cWPFihUODg6WL414\n/fXXf/jhh0WLFh0/fjw8PNzDw+PcuXO7du0qc/qtIoKCgsLCwiZNmuTm5hYfH3/o0KHp06cH\nBwdXouMVNG3atMTExHfeeefUqVPdu3c/d+7c119/7e3tvXjxYutm0dHR8fHxEydO1Ol0UVFR\n1ptef/31vXv3Llu27Jdffnnuued8fX0zMzNTUlJSU1PLvcbORi/GjRs3Y8aM1atXN27ceNCg\nQX+yawCgJc3uxwVqoU2bNoWGhqoXrvXu3dv8x6M9Nm/eXKZlSUnJnDlzgoKCnJ2dmzVrNmXK\nlIsXLwohJk+erDYo940TJ04UQty+fVtdfPfdd7t37+7j4+Pk5NS0adMhQ4YcPHjQun1xcbH6\nBGODwaA+oDgmJsbG/s0PeNzJxo0bP/3008DAQGdn56CgoAULFlg/QqXcjttQka6Zzebbt29P\nmTLF39/fycnJ19d31KhRlgesWKgPtBNCWB5fV2br0qVLu3bt6ubmZjAYWrRoMWTIEMsDTcr0\n1HYvxo8fL4SYOXOm7a4BgJ1TzJWauAEgh8TExMjIyM2bNw8ZMkTrWrQ0adKk5cuXX7hwwXKF\nHwDURlxjB6Cuu3379tq1awcMGECqA1DbcY0dgEdjNpvLPB7FmrOzc5mHJNuz48ePnzx58ttv\nv7137576XRcAUKvVmt+/AOzEqVOnXB5s/fr1Whf4CNauXRsVFXX+/PklS5b07NlT63IA4M/i\nGjsAj6awsNDGjbdBQUHe3t41WQ8AwIJgBwAAIAmmYgEAACRBsAMAAJAEwQ4AAEASBDsAAABJ\nEOwAAAAkQbADAACQxP8DJYlbDfTDYkkAAAAASUVORK5CYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "read_annot$transcript_novelty <- factor(read_annot$transcript_novelty, \n",
    "                                        levels = c(\"Known\", \"ISM\", \"NIC\",\n",
    "                                                   \"NNC\",\"Antisense\",\"Intergenic\",\"Genomic\"))\n",
    "\n",
    "ggplot(data = read_annot) + \n",
    "    geom_bar(mapping = aes(x = transcript_novelty, fill = transcript_novelty)) + \n",
    "    theme(axis.text.x = element_text(angle = 45, hjust = 1)) + \n",
    "    scale_fill_manual(values = c(\"#009E73\",\"#0072B2\",\"#D55E00\",\"#E69F00\",\"#000000\",\"#CC79A7\",\"#F0E442\")) +\n",
    "    ylab(\"# of distinct transcript models\") + theme_minimal()\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "4.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}