Comparison_plot (Sandoval-Velasco, 2019).ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "toc": true
   },
   "source": [
    "<h1>Table of Contents<span class=\"tocSkip\"></span></h1>\n",
    "<div class=\"toc\"><ul class=\"toc-item\"><li><span><a href=\"#Figure-1:-Comparing-BGI-obtained-values-for-3-samples-against-statistics-for-a-set-of-Illumina-samples-(n=320)\" data-toc-modified-id=\"Figure-1:-Comparing-BGI-obtained-values-for-3-samples-against-statistics-for-a-set-of-Illumina-samples-(n=320)-1\"><span class=\"toc-item-num\">1&nbsp;&nbsp;</span>Figure 1: Comparing BGI obtained values for 3 samples against statistics for a set of Illumina samples (n=320)</a></span><ul class=\"toc-item\"><li><span><a href=\"#Formatting-the-data-from-the-4D-genome-unit:\" data-toc-modified-id=\"Formatting-the-data-from-the-4D-genome-unit:-1.1\"><span class=\"toc-item-num\">1.1&nbsp;&nbsp;</span>Formatting the data from the 4D genome unit:</a></span></li><li><span><a href=\"#Read-the-data-from-the-4D-genome-unit.\" data-toc-modified-id=\"Read-the-data-from-the-4D-genome-unit.-1.2\"><span class=\"toc-item-num\">1.2&nbsp;&nbsp;</span>Read the data from the 4D genome unit.</a></span><ul class=\"toc-item\"><li><span><a href=\"#Placing-the-sample-values-on-the-graph\" data-toc-modified-id=\"Placing-the-sample-values-on-the-graph-1.2.1\"><span class=\"toc-item-num\">1.2.1&nbsp;&nbsp;</span>Placing the sample values on the graph</a></span></li></ul></li><li><span><a href=\"#Saving-tabular-data-for-Figure-1\" data-toc-modified-id=\"Saving-tabular-data-for-Figure-1-1.3\"><span class=\"toc-item-num\">1.3&nbsp;&nbsp;</span>Saving tabular data for Figure 1</a></span></li></ul></li></ul></div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Figure 1: Comparing BGI obtained values for 3 samples against statistics for a set of Illumina samples (n=320) "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Formatting the data from the 4D genome unit:\n",
    "\n",
    "This raw data is private and cannot be published yet. We used it with permission from the 4D Genome unit. \n",
    "\n",
    "For any inquiries email to: `martirenom [at] cnag [dot] crg [dot] eu`"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Removing those that were not classified as INHIC (insitu HiC), and also those that had 0 reads mapped (most of them were tests & prototypes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "NOTE: **TO DO IN BASH!**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cd /home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder\n",
    "\n",
    "# check column names\n",
    "head -n1 hic_juanro_dataset.tsv | sed 's/\\t/\\n/g' | cat -n\n",
    "\n",
    "# Removing those that are not insitu hic, have 0 valid pairs or are test samples\n",
    "cat \n",
    "<(head -n1 hic_juanro_dataset.tsv | awk -F'\\t' -v OFS=$'\\t' '{print $1,$45,$46,$47,$48,$67,$69,$70,$71,$72,$73,$74,$75,$76,$77,$78,$79,$98,$99}' ) \\\n",
    "<(awk -F'\\t' -v OFS=$'\\t' '$3 == \"INHIC\" && $99 != 0 && $1 != \"TE_S_T\"' hic_juanro_dataset.tsv \\| \n",
    "  awk -F'\\t' -v OFS=$'\\t' '{print $1,$45,$46,$47,$48,$67,$69,$70,$71,$72,$73,$74,$75,$76,$77,$78,$79,$98,$99}') > hic_juanro_dataset_f1.tsv\n",
    "\n",
    "wc -l hic_juanro_dataset_f1.tsv"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Now run the analysis in R:**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "\n",
      "Attaching package: ‘dplyr’\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:stats’:\n",
      "\n",
      "    filter, lag\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:base’:\n",
      "\n",
      "    intersect, setdiff, setequal, union\n",
      "\n",
      "\n",
      "\n",
      "Attaching package: ‘reshape2’\n",
      "\n",
      "\n",
      "The following object is masked from ‘package:tidyr’:\n",
      "\n",
      "    smiths\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "library(ggplot2)\n",
    "library(dplyr)\n",
    "library(tidyr)\n",
    "library(reshape2)\n",
    "library(scales)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Read the data from the 4D genome unit. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "a=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/hic_juanro_dataset_f1.tsv\",\n",
    "             header=TRUE,\n",
    "             stringsAsFactors = FALSE)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>320</li><li>19</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 320\n",
       "\\item 19\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 320\n",
       "2. 19\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 320  19"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#summary(a)\n",
    "dim(a)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 10</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>EXCLUDED_SELF_CIRCLE</th><th scope=col>EXCLUDED_DANGLING_END</th><th scope=col>EXCLUDED_ERROR</th><th scope=col>EXCLUDED_EXTRA_DANGLING_END</th><th scope=col>EXCLUDED_TOO_CLOSE_FROM_RES</th><th scope=col>EXCLUDED_TOO_SHORT</th><th scope=col>EXCLUDED_TOO_LARGE</th><th scope=col>EXCLUDED_OVER_REPRESENTED</th><th scope=col>EXCLUDED_DUPLICATED</th><th scope=col>EXCLUDED_RANDOM_BREAKS</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>180727;0.0008181000356266231</td><td>13964364;0.06321272795931505 </td><td>86353;0.00039089561812272536</td><td>34461646;0.15599812878182048</td><td>71165405;0.32214566924633875</td><td>13074145;0.05918295822034138</td><td>27994;0.0001267209238095674 </td><td>9420193;0.04264255052598486 </td><td>8891413;0.04024891295750507 </td><td>123792;0.0005603713867340847</td></tr>\n",
       "\t<tr><th scope=row>2</th><td>176738;0.0008162783164235938</td><td>17665621;0.08159005628929424 </td><td>33312;0.00015385408501116204</td><td>38540540;0.17800250713064636</td><td>55525632;0.2564494868523805 </td><td>10077260;0.04654261577568392</td><td>34522;0.00015944256492421159</td><td>11475967;0.05300265377051183</td><td>15688726;0.07245961166308923</td><td>158380;0.0007314904534122191</td></tr>\n",
       "\t<tr><th scope=row>3</th><td>110767;0.0007128736641671585</td><td>12852246;0.08271441583502041 </td><td>718378;0.004623333277213204 </td><td>39532601;0.25442370136347714</td><td>72535121;0.4668211424709363 </td><td>11660154;0.07504235653880081</td><td>24329;0.0001565764476380402 </td><td>9008993;0.05798002880249787 </td><td>18836436;0.12122743372277099</td><td>145135;0.0009340590541307478</td></tr>\n",
       "\t<tr><th scope=row>4</th><td>294303;0.0013325714223721011</td><td>12517567;0.056678158434769865</td><td>231953;0.0010502575207642327</td><td>46079136;0.20864122962116421</td><td>74787736;0.338630594107125  </td><td>11548820;0.05229177920075354</td><td>49365;0.00022351925826579673</td><td>13857740;0.06274631350228425</td><td>36123427;0.1635628807669129 </td><td>265574;0.0012024896889431926</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>176931;0.0008646272499276631</td><td>27806915;0.13588696410138576 </td><td>1399862;0.006840852260701846</td><td>43149998;0.21086561487316616</td><td>84818738;0.41449261112679503</td><td>12416541;0.06067721144651913</td><td>44783;0.00021884577679157714</td><td>10745709;0.05251218170469244</td><td>54506386;0.2663620656113156 </td><td>1193757;0.005833657369211146</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>119490;0.0007923299385199313</td><td>29990011;0.1988616919561642  </td><td>544565;0.0036109729096167574</td><td>39737860;0.2634990055294471 </td><td>61974834;0.4109508445309477 </td><td>9770407;0.06478689411352169 </td><td>35375;0.00023456918215032695</td><td>8266866;0.054817017570780065</td><td>32920930;0.21829641344814596</td><td>845282;0.00560500657035739  </td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 10\n",
       "\\begin{tabular}{r|llllllllll}\n",
       "  & EXCLUDED\\_SELF\\_CIRCLE & EXCLUDED\\_DANGLING\\_END & EXCLUDED\\_ERROR & EXCLUDED\\_EXTRA\\_DANGLING\\_END & EXCLUDED\\_TOO\\_CLOSE\\_FROM\\_RES & EXCLUDED\\_TOO\\_SHORT & EXCLUDED\\_TOO\\_LARGE & EXCLUDED\\_OVER\\_REPRESENTED & EXCLUDED\\_DUPLICATED & EXCLUDED\\_RANDOM\\_BREAKS\\\\\n",
       "  & <chr> & <chr> & <chr> & <chr> & <chr> & <chr> & <chr> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & 180727;0.0008181000356266231 & 13964364;0.06321272795931505  & 86353;0.00039089561812272536 & 34461646;0.15599812878182048 & 71165405;0.32214566924633875 & 13074145;0.05918295822034138 & 27994;0.0001267209238095674  & 9420193;0.04264255052598486  & 8891413;0.04024891295750507  & 123792;0.0005603713867340847\\\\\n",
       "\t2 & 176738;0.0008162783164235938 & 17665621;0.08159005628929424  & 33312;0.00015385408501116204 & 38540540;0.17800250713064636 & 55525632;0.2564494868523805  & 10077260;0.04654261577568392 & 34522;0.00015944256492421159 & 11475967;0.05300265377051183 & 15688726;0.07245961166308923 & 158380;0.0007314904534122191\\\\\n",
       "\t3 & 110767;0.0007128736641671585 & 12852246;0.08271441583502041  & 718378;0.004623333277213204  & 39532601;0.25442370136347714 & 72535121;0.4668211424709363  & 11660154;0.07504235653880081 & 24329;0.0001565764476380402  & 9008993;0.05798002880249787  & 18836436;0.12122743372277099 & 145135;0.0009340590541307478\\\\\n",
       "\t4 & 294303;0.0013325714223721011 & 12517567;0.056678158434769865 & 231953;0.0010502575207642327 & 46079136;0.20864122962116421 & 74787736;0.338630594107125   & 11548820;0.05229177920075354 & 49365;0.00022351925826579673 & 13857740;0.06274631350228425 & 36123427;0.1635628807669129  & 265574;0.0012024896889431926\\\\\n",
       "\t5 & 176931;0.0008646272499276631 & 27806915;0.13588696410138576  & 1399862;0.006840852260701846 & 43149998;0.21086561487316616 & 84818738;0.41449261112679503 & 12416541;0.06067721144651913 & 44783;0.00021884577679157714 & 10745709;0.05251218170469244 & 54506386;0.2663620656113156  & 1193757;0.005833657369211146\\\\\n",
       "\t6 & 119490;0.0007923299385199313 & 29990011;0.1988616919561642   & 544565;0.0036109729096167574 & 39737860;0.2634990055294471  & 61974834;0.4109508445309477  & 9770407;0.06478689411352169  & 35375;0.00023456918215032695 & 8266866;0.054817017570780065 & 32920930;0.21829641344814596 & 845282;0.00560500657035739  \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 10\n",
       "\n",
       "| <!--/--> | EXCLUDED_SELF_CIRCLE &lt;chr&gt; | EXCLUDED_DANGLING_END &lt;chr&gt; | EXCLUDED_ERROR &lt;chr&gt; | EXCLUDED_EXTRA_DANGLING_END &lt;chr&gt; | EXCLUDED_TOO_CLOSE_FROM_RES &lt;chr&gt; | EXCLUDED_TOO_SHORT &lt;chr&gt; | EXCLUDED_TOO_LARGE &lt;chr&gt; | EXCLUDED_OVER_REPRESENTED &lt;chr&gt; | EXCLUDED_DUPLICATED &lt;chr&gt; | EXCLUDED_RANDOM_BREAKS &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|---|---|\n",
       "| 1 | 180727;0.0008181000356266231 | 13964364;0.06321272795931505  | 86353;0.00039089561812272536 | 34461646;0.15599812878182048 | 71165405;0.32214566924633875 | 13074145;0.05918295822034138 | 27994;0.0001267209238095674  | 9420193;0.04264255052598486  | 8891413;0.04024891295750507  | 123792;0.0005603713867340847 |\n",
       "| 2 | 176738;0.0008162783164235938 | 17665621;0.08159005628929424  | 33312;0.00015385408501116204 | 38540540;0.17800250713064636 | 55525632;0.2564494868523805  | 10077260;0.04654261577568392 | 34522;0.00015944256492421159 | 11475967;0.05300265377051183 | 15688726;0.07245961166308923 | 158380;0.0007314904534122191 |\n",
       "| 3 | 110767;0.0007128736641671585 | 12852246;0.08271441583502041  | 718378;0.004623333277213204  | 39532601;0.25442370136347714 | 72535121;0.4668211424709363  | 11660154;0.07504235653880081 | 24329;0.0001565764476380402  | 9008993;0.05798002880249787  | 18836436;0.12122743372277099 | 145135;0.0009340590541307478 |\n",
       "| 4 | 294303;0.0013325714223721011 | 12517567;0.056678158434769865 | 231953;0.0010502575207642327 | 46079136;0.20864122962116421 | 74787736;0.338630594107125   | 11548820;0.05229177920075354 | 49365;0.00022351925826579673 | 13857740;0.06274631350228425 | 36123427;0.1635628807669129  | 265574;0.0012024896889431926 |\n",
       "| 5 | 176931;0.0008646272499276631 | 27806915;0.13588696410138576  | 1399862;0.006840852260701846 | 43149998;0.21086561487316616 | 84818738;0.41449261112679503 | 12416541;0.06067721144651913 | 44783;0.00021884577679157714 | 10745709;0.05251218170469244 | 54506386;0.2663620656113156  | 1193757;0.005833657369211146 |\n",
       "| 6 | 119490;0.0007923299385199313 | 29990011;0.1988616919561642   | 544565;0.0036109729096167574 | 39737860;0.2634990055294471  | 61974834;0.4109508445309477  | 9770407;0.06478689411352169  | 35375;0.00023456918215032695 | 8266866;0.054817017570780065 | 32920930;0.21829641344814596 | 845282;0.00560500657035739   |\n",
       "\n"
      ],
      "text/plain": [
       "  EXCLUDED_SELF_CIRCLE         EXCLUDED_DANGLING_END        \n",
       "1 180727;0.0008181000356266231 13964364;0.06321272795931505 \n",
       "2 176738;0.0008162783164235938 17665621;0.08159005628929424 \n",
       "3 110767;0.0007128736641671585 12852246;0.08271441583502041 \n",
       "4 294303;0.0013325714223721011 12517567;0.056678158434769865\n",
       "5 176931;0.0008646272499276631 27806915;0.13588696410138576 \n",
       "6 119490;0.0007923299385199313 29990011;0.1988616919561642  \n",
       "  EXCLUDED_ERROR               EXCLUDED_EXTRA_DANGLING_END \n",
       "1 86353;0.00039089561812272536 34461646;0.15599812878182048\n",
       "2 33312;0.00015385408501116204 38540540;0.17800250713064636\n",
       "3 718378;0.004623333277213204  39532601;0.25442370136347714\n",
       "4 231953;0.0010502575207642327 46079136;0.20864122962116421\n",
       "5 1399862;0.006840852260701846 43149998;0.21086561487316616\n",
       "6 544565;0.0036109729096167574 39737860;0.2634990055294471 \n",
       "  EXCLUDED_TOO_CLOSE_FROM_RES  EXCLUDED_TOO_SHORT          \n",
       "1 71165405;0.32214566924633875 13074145;0.05918295822034138\n",
       "2 55525632;0.2564494868523805  10077260;0.04654261577568392\n",
       "3 72535121;0.4668211424709363  11660154;0.07504235653880081\n",
       "4 74787736;0.338630594107125   11548820;0.05229177920075354\n",
       "5 84818738;0.41449261112679503 12416541;0.06067721144651913\n",
       "6 61974834;0.4109508445309477  9770407;0.06478689411352169 \n",
       "  EXCLUDED_TOO_LARGE           EXCLUDED_OVER_REPRESENTED   \n",
       "1 27994;0.0001267209238095674  9420193;0.04264255052598486 \n",
       "2 34522;0.00015944256492421159 11475967;0.05300265377051183\n",
       "3 24329;0.0001565764476380402  9008993;0.05798002880249787 \n",
       "4 49365;0.00022351925826579673 13857740;0.06274631350228425\n",
       "5 44783;0.00021884577679157714 10745709;0.05251218170469244\n",
       "6 35375;0.00023456918215032695 8266866;0.054817017570780065\n",
       "  EXCLUDED_DUPLICATED          EXCLUDED_RANDOM_BREAKS      \n",
       "1 8891413;0.04024891295750507  123792;0.0005603713867340847\n",
       "2 15688726;0.07245961166308923 158380;0.0007314904534122191\n",
       "3 18836436;0.12122743372277099 145135;0.0009340590541307478\n",
       "4 36123427;0.1635628807669129  265574;0.0012024896889431926\n",
       "5 54506386;0.2663620656113156  1193757;0.005833657369211146\n",
       "6 32920930;0.21829641344814596 845282;0.00560500657035739  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message in lapply(a[, 2:5], as.double):\n",
      "“NAs introduced by coercion”\n",
      "Warning message in lapply(a[, 2:5], as.double):\n",
      "“NAs introduced by coercion”\n",
      "Warning message in lapply(a[, 2:5], as.double):\n",
      "“NAs introduced by coercion”\n",
      "Warning message in lapply(a[, 2:5], as.double):\n",
      "“NAs introduced by coercion”\n"
     ]
    }
   ],
   "source": [
    "# select only those columns to compute the percentages.\n",
    "chunk=a[,8:17]\n",
    "head(chunk)\n",
    "# compute valid pairs\n",
    "valids=a[,19]/a[,18]\n",
    "# change character for double\n",
    "percs=as.data.frame(lapply(a[,2:5],as.double))\n",
    "# compute percentage over 1\n",
    "percs=percs/100"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Split the columns and create a new column besides each filter class with the percentage reported there\n",
    "for (i in seq(1,dim(chunk)[2],1)){\n",
    "    ca=data.frame(chunk[,i],stringsAsFactors = FALSE)\n",
    "    cn=colnames(chunk)[i]\n",
    "    colnames(ca)=cn\n",
    "    caca=ca %>% separate(colnames(chunk)[i], c(cn, paste0(cn,\"_perc\")),sep=\";\")\n",
    "    px=as.data.frame(caca[,2],stringsAsFactors = FALSE)\n",
    "    colnames(px)=paste0(cn,\"_perc\")\n",
    "    percs=cbind(percs,lapply(px,as.numeric))\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# add valid pairs\n",
    "percs$valid_pairs=valids"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# column names formatting\n",
    "colnames(percs)=c(\"dangling ends r1\",\"ligated r1\",\"dangling ends r2\",\"ligated r2\",\"self circle\",\"dangling ends\",\"error\",\"extra dangling\", \"too close from res\", \"too short\", \"too large\", \"overrepresented\",\"duplicated\", \"random breaks\",\"valid pairs\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>320</li><li>15</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 320\n",
       "\\item 15\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 320\n",
       "2. 15\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 320  15"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# How many we have\n",
    "dim(percs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>316</li><li>15</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 316\n",
       "\\item 15\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 316\n",
       "2. 15\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 316  15"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 15</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>dangling ends r1</th><th scope=col>ligated r1</th><th scope=col>dangling ends r2</th><th scope=col>ligated r2</th><th scope=col>self circle</th><th scope=col>dangling ends</th><th scope=col>error</th><th scope=col>extra dangling</th><th scope=col>too close from res</th><th scope=col>too short</th><th scope=col>too large</th><th scope=col>overrepresented</th><th scope=col>duplicated</th><th scope=col>random breaks</th><th scope=col>valid pairs</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>0.07331</td><td>0.23132</td><td>0.08017</td><td>0.23060</td><td>0.0008181000</td><td>0.06321273</td><td>0.0003908956</td><td>0.1559981</td><td>0.3221457</td><td>0.05918296</td><td>0.0001267209</td><td>0.04264255</td><td>0.04024891</td><td>0.0005603714</td><td>0.8864393</td></tr>\n",
       "\t<tr><th scope=row>2</th><td>0.03291</td><td>0.25579</td><td>0.03056</td><td>0.26364</td><td>0.0008162783</td><td>0.08159006</td><td>0.0001538541</td><td>0.1780025</td><td>0.2564495</td><td>0.04654262</td><td>0.0001594426</td><td>0.05300265</td><td>0.07245961</td><td>0.0007314905</td><td>0.8284620</td></tr>\n",
       "\t<tr><th scope=row>3</th><td>0.06644</td><td>0.27054</td><td>0.06821</td><td>0.26698</td><td>0.0007128737</td><td>0.08271442</td><td>0.0046233333</td><td>0.2544237</td><td>0.4668211</td><td>0.07504236</td><td>0.0001565764</td><td>0.05798003</td><td>0.12122743</td><td>0.0009340591</td><td>0.7867297</td></tr>\n",
       "\t<tr><th scope=row>4</th><td>0.02121</td><td>0.25370</td><td>0.02163</td><td>0.25563</td><td>0.0013325714</td><td>0.05667816</td><td>0.0010502575</td><td>0.2086412</td><td>0.3386306</td><td>0.05229178</td><td>0.0002235193</td><td>0.06274631</td><td>0.16356288</td><td>0.0012024897</td><td>0.8000160</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>0.07891</td><td>0.21146</td><td>0.08121</td><td>0.20882</td><td>0.0008646272</td><td>0.13588696</td><td>0.0068408523</td><td>0.2108656</td><td>0.4144926</td><td>0.06067721</td><td>0.0002188458</td><td>0.05251218</td><td>0.26636207</td><td>0.0058336574</td><td>0.6612235</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>0.03578</td><td>0.27713</td><td>0.03734</td><td>0.26850</td><td>0.0007923299</td><td>0.19886169</td><td>0.0036109729</td><td>0.2634990</td><td>0.4109508</td><td>0.06478689</td><td>0.0002345692</td><td>0.05481702</td><td>0.21829641</td><td>0.0056050066</td><td>0.6483567</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 15\n",
       "\\begin{tabular}{r|lllllllllllllll}\n",
       "  & dangling ends r1 & ligated r1 & dangling ends r2 & ligated r2 & self circle & dangling ends & error & extra dangling & too close from res & too short & too large & overrepresented & duplicated & random breaks & valid pairs\\\\\n",
       "  & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl> & <dbl>\\\\\n",
       "\\hline\n",
       "\t1 & 0.07331 & 0.23132 & 0.08017 & 0.23060 & 0.0008181000 & 0.06321273 & 0.0003908956 & 0.1559981 & 0.3221457 & 0.05918296 & 0.0001267209 & 0.04264255 & 0.04024891 & 0.0005603714 & 0.8864393\\\\\n",
       "\t2 & 0.03291 & 0.25579 & 0.03056 & 0.26364 & 0.0008162783 & 0.08159006 & 0.0001538541 & 0.1780025 & 0.2564495 & 0.04654262 & 0.0001594426 & 0.05300265 & 0.07245961 & 0.0007314905 & 0.8284620\\\\\n",
       "\t3 & 0.06644 & 0.27054 & 0.06821 & 0.26698 & 0.0007128737 & 0.08271442 & 0.0046233333 & 0.2544237 & 0.4668211 & 0.07504236 & 0.0001565764 & 0.05798003 & 0.12122743 & 0.0009340591 & 0.7867297\\\\\n",
       "\t4 & 0.02121 & 0.25370 & 0.02163 & 0.25563 & 0.0013325714 & 0.05667816 & 0.0010502575 & 0.2086412 & 0.3386306 & 0.05229178 & 0.0002235193 & 0.06274631 & 0.16356288 & 0.0012024897 & 0.8000160\\\\\n",
       "\t5 & 0.07891 & 0.21146 & 0.08121 & 0.20882 & 0.0008646272 & 0.13588696 & 0.0068408523 & 0.2108656 & 0.4144926 & 0.06067721 & 0.0002188458 & 0.05251218 & 0.26636207 & 0.0058336574 & 0.6612235\\\\\n",
       "\t6 & 0.03578 & 0.27713 & 0.03734 & 0.26850 & 0.0007923299 & 0.19886169 & 0.0036109729 & 0.2634990 & 0.4109508 & 0.06478689 & 0.0002345692 & 0.05481702 & 0.21829641 & 0.0056050066 & 0.6483567\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 15\n",
       "\n",
       "| <!--/--> | dangling ends r1 &lt;dbl&gt; | ligated r1 &lt;dbl&gt; | dangling ends r2 &lt;dbl&gt; | ligated r2 &lt;dbl&gt; | self circle &lt;dbl&gt; | dangling ends &lt;dbl&gt; | error &lt;dbl&gt; | extra dangling &lt;dbl&gt; | too close from res &lt;dbl&gt; | too short &lt;dbl&gt; | too large &lt;dbl&gt; | overrepresented &lt;dbl&gt; | duplicated &lt;dbl&gt; | random breaks &lt;dbl&gt; | valid pairs &lt;dbl&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|\n",
       "| 1 | 0.07331 | 0.23132 | 0.08017 | 0.23060 | 0.0008181000 | 0.06321273 | 0.0003908956 | 0.1559981 | 0.3221457 | 0.05918296 | 0.0001267209 | 0.04264255 | 0.04024891 | 0.0005603714 | 0.8864393 |\n",
       "| 2 | 0.03291 | 0.25579 | 0.03056 | 0.26364 | 0.0008162783 | 0.08159006 | 0.0001538541 | 0.1780025 | 0.2564495 | 0.04654262 | 0.0001594426 | 0.05300265 | 0.07245961 | 0.0007314905 | 0.8284620 |\n",
       "| 3 | 0.06644 | 0.27054 | 0.06821 | 0.26698 | 0.0007128737 | 0.08271442 | 0.0046233333 | 0.2544237 | 0.4668211 | 0.07504236 | 0.0001565764 | 0.05798003 | 0.12122743 | 0.0009340591 | 0.7867297 |\n",
       "| 4 | 0.02121 | 0.25370 | 0.02163 | 0.25563 | 0.0013325714 | 0.05667816 | 0.0010502575 | 0.2086412 | 0.3386306 | 0.05229178 | 0.0002235193 | 0.06274631 | 0.16356288 | 0.0012024897 | 0.8000160 |\n",
       "| 5 | 0.07891 | 0.21146 | 0.08121 | 0.20882 | 0.0008646272 | 0.13588696 | 0.0068408523 | 0.2108656 | 0.4144926 | 0.06067721 | 0.0002188458 | 0.05251218 | 0.26636207 | 0.0058336574 | 0.6612235 |\n",
       "| 6 | 0.03578 | 0.27713 | 0.03734 | 0.26850 | 0.0007923299 | 0.19886169 | 0.0036109729 | 0.2634990 | 0.4109508 | 0.06478689 | 0.0002345692 | 0.05481702 | 0.21829641 | 0.0056050066 | 0.6483567 |\n",
       "\n"
      ],
      "text/plain": [
       "  dangling ends r1 ligated r1 dangling ends r2 ligated r2 self circle \n",
       "1 0.07331          0.23132    0.08017          0.23060    0.0008181000\n",
       "2 0.03291          0.25579    0.03056          0.26364    0.0008162783\n",
       "3 0.06644          0.27054    0.06821          0.26698    0.0007128737\n",
       "4 0.02121          0.25370    0.02163          0.25563    0.0013325714\n",
       "5 0.07891          0.21146    0.08121          0.20882    0.0008646272\n",
       "6 0.03578          0.27713    0.03734          0.26850    0.0007923299\n",
       "  dangling ends error        extra dangling too close from res too short \n",
       "1 0.06321273    0.0003908956 0.1559981      0.3221457          0.05918296\n",
       "2 0.08159006    0.0001538541 0.1780025      0.2564495          0.04654262\n",
       "3 0.08271442    0.0046233333 0.2544237      0.4668211          0.07504236\n",
       "4 0.05667816    0.0010502575 0.2086412      0.3386306          0.05229178\n",
       "5 0.13588696    0.0068408523 0.2108656      0.4144926          0.06067721\n",
       "6 0.19886169    0.0036109729 0.2634990      0.4109508          0.06478689\n",
       "  too large    overrepresented duplicated random breaks valid pairs\n",
       "1 0.0001267209 0.04264255      0.04024891 0.0005603714  0.8864393  \n",
       "2 0.0001594426 0.05300265      0.07245961 0.0007314905  0.8284620  \n",
       "3 0.0001565764 0.05798003      0.12122743 0.0009340591  0.7867297  \n",
       "4 0.0002235193 0.06274631      0.16356288 0.0012024897  0.8000160  \n",
       "5 0.0002188458 0.05251218      0.26636207 0.0058336574  0.6612235  \n",
       "6 0.0002345692 0.05481702      0.21829641 0.0056050066  0.6483567  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# use only those rows with full info (remove rows with NA's)\n",
    "percs=percs[complete.cases(percs), ]\n",
    "# After removing NA's (4 cases)\n",
    "dim(percs)\n",
    "head(percs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Placing the sample values on the graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# In bash; for each sample (workdir):\n",
    "cd ${wd};\n",
    "# Get the values for the errors in sample\n",
    "tadbit describe . | grep -A 4 \"INTERSECTION_OUTPUTs\" | tail -n1 | cut -d$'|' -f4 | sed 's/,//g' | awk '{print $1}' > tots;\n",
    "# Get the total number of interactions\n",
    "tadbit describe . | grep -A 14 \"FILTER_OUTPUTs\" | tail -n11 | cut -d$'|' -f5  | sed 's/,//g' | awk '{print $1}' > nums\n",
    "# Get the dangling ends and ligated sites\n",
    "tadbit describe . | grep -A6 'MAPPED_INPUTs' | tail -n4 | cut -d$'|' -f9,10 | grep -Eo '[0-9\\.]*' | sed -e 1,3b -e '$!d' > dangs_lig"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# disable scientific notation\n",
    "options(scipen=999)\n",
    "\n",
    "# enable again\n",
    "#options(scipen=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# In R, read both files for all\n",
    "vals10=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/10_Zebra_finch_Oz13/nums\")\n",
    "total10=as.numeric(read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/10_Zebra_finch_Oz13/tots\"))\n",
    "dangs_lig10=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/10_Zebra_finch_Oz13/dangs_lig\")\n",
    "\n",
    "vals18=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/18_Zebra_finch_Mz13/nums\")\n",
    "total18=as.numeric(read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/18_Zebra_finch_Mz13/tots\"))\n",
    "dangs_lig18=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/18_Zebra_finch_Mz13/dangs_lig\")\n",
    "\n",
    "vals19=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/19_Zebra_finch_Mz17/nums\")\n",
    "total19=as.numeric(read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/19_Zebra_finch_Mz17/tots\"))\n",
    "dangs_lig19=read.table(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/19_Zebra_finch_Mz17/dangs_lig\")\n",
    "\n",
    "vals10=vals10/total10\n",
    "vals10=bind_rows(dangs_lig10/100,vals10)\n",
    "rownames(vals10)=colnames(percs)\n",
    "\n",
    "vals18=vals18/total18\n",
    "vals18=bind_rows(dangs_lig18/100,vals18)\n",
    "rownames(vals18)=colnames(percs)\n",
    "\n",
    "vals19=vals19/total19\n",
    "vals19=bind_rows(dangs_lig19/100,vals19)\n",
    "rownames(vals19)=colnames(percs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "No id variables; using all as measure variables\n",
      "\n",
      "No id variables; using all as measure variables\n",
      "\n",
      "No id variables; using all as measure variables\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# Merge all samples\n",
    "colnames(vals10)=c(\"value\")\n",
    "colnames(vals18)=c(\"value\")\n",
    "colnames(vals19)=c(\"value\")\n",
    "\n",
    "a10=melt(vals10)\n",
    "a18=melt(vals18)\n",
    "a19=melt(vals19)\n",
    "\n",
    "a10$variable=as.factor(colnames(percs))\n",
    "a10$Sample=\"Oz13\"\n",
    "a18$variable=as.factor(colnames(percs))\n",
    "a18$Sample=\"Mz13\"\n",
    "a19$variable=as.factor(colnames(percs))\n",
    "a19$Sample=\"Mz17\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "all=bind_rows(a10,a18,a19)\n",
    "\n",
    "all$Sample=as.factor(all$Sample)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 45 × 3</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>variable</th><th scope=col>value</th><th scope=col>Sample</th></tr>\n",
       "\t<tr><th scope=col>&lt;fct&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;fct&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td>dangling ends r1  </td><td>0.030000000000</td><td>Oz13</td></tr>\n",
       "\t<tr><td>ligated r1        </td><td>0.316010000000</td><td>Oz13</td></tr>\n",
       "\t<tr><td>dangling ends r2  </td><td>0.030000000000</td><td>Oz13</td></tr>\n",
       "\t<tr><td>ligated r2        </td><td>0.298200000000</td><td>Oz13</td></tr>\n",
       "\t<tr><td>self circle       </td><td>0.000623176006</td><td>Oz13</td></tr>\n",
       "\t<tr><td>dangling ends     </td><td>0.332374768033</td><td>Oz13</td></tr>\n",
       "\t<tr><td>error             </td><td>0.000489399795</td><td>Oz13</td></tr>\n",
       "\t<tr><td>extra dangling    </td><td>0.288776641185</td><td>Oz13</td></tr>\n",
       "\t<tr><td>too close from res</td><td>0.292769182340</td><td>Oz13</td></tr>\n",
       "\t<tr><td>too short         </td><td>0.050580325471</td><td>Oz13</td></tr>\n",
       "\t<tr><td>too large         </td><td>0.000001149468</td><td>Oz13</td></tr>\n",
       "\t<tr><td>overrepresented   </td><td>0.184060958601</td><td>Oz13</td></tr>\n",
       "\t<tr><td>duplicated        </td><td>0.746469175613</td><td>Oz13</td></tr>\n",
       "\t<tr><td>random breaks     </td><td>0.053636871044</td><td>Oz13</td></tr>\n",
       "\t<tr><td>valid pairs       </td><td>0.221349489967</td><td>Oz13</td></tr>\n",
       "\t<tr><td>dangling ends r1  </td><td>0.050000000000</td><td>Mz13</td></tr>\n",
       "\t<tr><td>ligated r1        </td><td>0.148030000000</td><td>Mz13</td></tr>\n",
       "\t<tr><td>dangling ends r2  </td><td>0.050000000000</td><td>Mz13</td></tr>\n",
       "\t<tr><td>ligated r2        </td><td>0.144440000000</td><td>Mz13</td></tr>\n",
       "\t<tr><td>self circle       </td><td>0.000245294992</td><td>Mz13</td></tr>\n",
       "\t<tr><td>dangling ends     </td><td>0.424299467028</td><td>Mz13</td></tr>\n",
       "\t<tr><td>error             </td><td>0.000266190492</td><td>Mz13</td></tr>\n",
       "\t<tr><td>extra dangling    </td><td>0.337880224272</td><td>Mz13</td></tr>\n",
       "\t<tr><td>too close from res</td><td>0.292011877729</td><td>Mz13</td></tr>\n",
       "\t<tr><td>too short         </td><td>0.043126493631</td><td>Mz13</td></tr>\n",
       "\t<tr><td>too large         </td><td>0.000000454250</td><td>Mz13</td></tr>\n",
       "\t<tr><td>overrepresented   </td><td>0.202854870309</td><td>Mz13</td></tr>\n",
       "\t<tr><td>duplicated        </td><td>0.538828153142</td><td>Mz13</td></tr>\n",
       "\t<tr><td>random breaks     </td><td>0.063488703484</td><td>Mz13</td></tr>\n",
       "\t<tr><td>valid pairs       </td><td>0.383513723573</td><td>Mz13</td></tr>\n",
       "\t<tr><td>dangling ends r1  </td><td>0.020000000000</td><td>Mz17</td></tr>\n",
       "\t<tr><td>ligated r1        </td><td>0.078500000000</td><td>Mz17</td></tr>\n",
       "\t<tr><td>dangling ends r2  </td><td>0.020000000000</td><td>Mz17</td></tr>\n",
       "\t<tr><td>ligated r2        </td><td>0.073100000000</td><td>Mz17</td></tr>\n",
       "\t<tr><td>self circle       </td><td>0.000283001150</td><td>Mz17</td></tr>\n",
       "\t<tr><td>dangling ends     </td><td>0.502507776100</td><td>Mz17</td></tr>\n",
       "\t<tr><td>error             </td><td>0.000377549261</td><td>Mz17</td></tr>\n",
       "\t<tr><td>extra dangling    </td><td>0.332689720112</td><td>Mz17</td></tr>\n",
       "\t<tr><td>too close from res</td><td>0.142757357387</td><td>Mz17</td></tr>\n",
       "\t<tr><td>too short         </td><td>0.028925290269</td><td>Mz17</td></tr>\n",
       "\t<tr><td>too large         </td><td>0.000000000000</td><td>Mz17</td></tr>\n",
       "\t<tr><td>overrepresented   </td><td>0.041307876951</td><td>Mz17</td></tr>\n",
       "\t<tr><td>duplicated        </td><td>0.166845900728</td><td>Mz17</td></tr>\n",
       "\t<tr><td>random breaks     </td><td>0.113966492664</td><td>Mz17</td></tr>\n",
       "\t<tr><td>valid pairs       </td><td>0.717090182176</td><td>Mz17</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 45 × 3\n",
       "\\begin{tabular}{lll}\n",
       " variable & value & Sample\\\\\n",
       " <fct> & <dbl> & <fct>\\\\\n",
       "\\hline\n",
       "\t dangling ends r1   & 0.030000000000 & Oz13\\\\\n",
       "\t ligated r1         & 0.316010000000 & Oz13\\\\\n",
       "\t dangling ends r2   & 0.030000000000 & Oz13\\\\\n",
       "\t ligated r2         & 0.298200000000 & Oz13\\\\\n",
       "\t self circle        & 0.000623176006 & Oz13\\\\\n",
       "\t dangling ends      & 0.332374768033 & Oz13\\\\\n",
       "\t error              & 0.000489399795 & Oz13\\\\\n",
       "\t extra dangling     & 0.288776641185 & Oz13\\\\\n",
       "\t too close from res & 0.292769182340 & Oz13\\\\\n",
       "\t too short          & 0.050580325471 & Oz13\\\\\n",
       "\t too large          & 0.000001149468 & Oz13\\\\\n",
       "\t overrepresented    & 0.184060958601 & Oz13\\\\\n",
       "\t duplicated         & 0.746469175613 & Oz13\\\\\n",
       "\t random breaks      & 0.053636871044 & Oz13\\\\\n",
       "\t valid pairs        & 0.221349489967 & Oz13\\\\\n",
       "\t dangling ends r1   & 0.050000000000 & Mz13\\\\\n",
       "\t ligated r1         & 0.148030000000 & Mz13\\\\\n",
       "\t dangling ends r2   & 0.050000000000 & Mz13\\\\\n",
       "\t ligated r2         & 0.144440000000 & Mz13\\\\\n",
       "\t self circle        & 0.000245294992 & Mz13\\\\\n",
       "\t dangling ends      & 0.424299467028 & Mz13\\\\\n",
       "\t error              & 0.000266190492 & Mz13\\\\\n",
       "\t extra dangling     & 0.337880224272 & Mz13\\\\\n",
       "\t too close from res & 0.292011877729 & Mz13\\\\\n",
       "\t too short          & 0.043126493631 & Mz13\\\\\n",
       "\t too large          & 0.000000454250 & Mz13\\\\\n",
       "\t overrepresented    & 0.202854870309 & Mz13\\\\\n",
       "\t duplicated         & 0.538828153142 & Mz13\\\\\n",
       "\t random breaks      & 0.063488703484 & Mz13\\\\\n",
       "\t valid pairs        & 0.383513723573 & Mz13\\\\\n",
       "\t dangling ends r1   & 0.020000000000 & Mz17\\\\\n",
       "\t ligated r1         & 0.078500000000 & Mz17\\\\\n",
       "\t dangling ends r2   & 0.020000000000 & Mz17\\\\\n",
       "\t ligated r2         & 0.073100000000 & Mz17\\\\\n",
       "\t self circle        & 0.000283001150 & Mz17\\\\\n",
       "\t dangling ends      & 0.502507776100 & Mz17\\\\\n",
       "\t error              & 0.000377549261 & Mz17\\\\\n",
       "\t extra dangling     & 0.332689720112 & Mz17\\\\\n",
       "\t too close from res & 0.142757357387 & Mz17\\\\\n",
       "\t too short          & 0.028925290269 & Mz17\\\\\n",
       "\t too large          & 0.000000000000 & Mz17\\\\\n",
       "\t overrepresented    & 0.041307876951 & Mz17\\\\\n",
       "\t duplicated         & 0.166845900728 & Mz17\\\\\n",
       "\t random breaks      & 0.113966492664 & Mz17\\\\\n",
       "\t valid pairs        & 0.717090182176 & Mz17\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 45 × 3\n",
       "\n",
       "| variable &lt;fct&gt; | value &lt;dbl&gt; | Sample &lt;fct&gt; |\n",
       "|---|---|---|\n",
       "| dangling ends r1   | 0.030000000000 | Oz13 |\n",
       "| ligated r1         | 0.316010000000 | Oz13 |\n",
       "| dangling ends r2   | 0.030000000000 | Oz13 |\n",
       "| ligated r2         | 0.298200000000 | Oz13 |\n",
       "| self circle        | 0.000623176006 | Oz13 |\n",
       "| dangling ends      | 0.332374768033 | Oz13 |\n",
       "| error              | 0.000489399795 | Oz13 |\n",
       "| extra dangling     | 0.288776641185 | Oz13 |\n",
       "| too close from res | 0.292769182340 | Oz13 |\n",
       "| too short          | 0.050580325471 | Oz13 |\n",
       "| too large          | 0.000001149468 | Oz13 |\n",
       "| overrepresented    | 0.184060958601 | Oz13 |\n",
       "| duplicated         | 0.746469175613 | Oz13 |\n",
       "| random breaks      | 0.053636871044 | Oz13 |\n",
       "| valid pairs        | 0.221349489967 | Oz13 |\n",
       "| dangling ends r1   | 0.050000000000 | Mz13 |\n",
       "| ligated r1         | 0.148030000000 | Mz13 |\n",
       "| dangling ends r2   | 0.050000000000 | Mz13 |\n",
       "| ligated r2         | 0.144440000000 | Mz13 |\n",
       "| self circle        | 0.000245294992 | Mz13 |\n",
       "| dangling ends      | 0.424299467028 | Mz13 |\n",
       "| error              | 0.000266190492 | Mz13 |\n",
       "| extra dangling     | 0.337880224272 | Mz13 |\n",
       "| too close from res | 0.292011877729 | Mz13 |\n",
       "| too short          | 0.043126493631 | Mz13 |\n",
       "| too large          | 0.000000454250 | Mz13 |\n",
       "| overrepresented    | 0.202854870309 | Mz13 |\n",
       "| duplicated         | 0.538828153142 | Mz13 |\n",
       "| random breaks      | 0.063488703484 | Mz13 |\n",
       "| valid pairs        | 0.383513723573 | Mz13 |\n",
       "| dangling ends r1   | 0.020000000000 | Mz17 |\n",
       "| ligated r1         | 0.078500000000 | Mz17 |\n",
       "| dangling ends r2   | 0.020000000000 | Mz17 |\n",
       "| ligated r2         | 0.073100000000 | Mz17 |\n",
       "| self circle        | 0.000283001150 | Mz17 |\n",
       "| dangling ends      | 0.502507776100 | Mz17 |\n",
       "| error              | 0.000377549261 | Mz17 |\n",
       "| extra dangling     | 0.332689720112 | Mz17 |\n",
       "| too close from res | 0.142757357387 | Mz17 |\n",
       "| too short          | 0.028925290269 | Mz17 |\n",
       "| too large          | 0.000000000000 | Mz17 |\n",
       "| overrepresented    | 0.041307876951 | Mz17 |\n",
       "| duplicated         | 0.166845900728 | Mz17 |\n",
       "| random breaks      | 0.113966492664 | Mz17 |\n",
       "| valid pairs        | 0.717090182176 | Mz17 |\n",
       "\n"
      ],
      "text/plain": [
       "   variable           value          Sample\n",
       "1  dangling ends r1   0.030000000000 Oz13  \n",
       "2  ligated r1         0.316010000000 Oz13  \n",
       "3  dangling ends r2   0.030000000000 Oz13  \n",
       "4  ligated r2         0.298200000000 Oz13  \n",
       "5  self circle        0.000623176006 Oz13  \n",
       "6  dangling ends      0.332374768033 Oz13  \n",
       "7  error              0.000489399795 Oz13  \n",
       "8  extra dangling     0.288776641185 Oz13  \n",
       "9  too close from res 0.292769182340 Oz13  \n",
       "10 too short          0.050580325471 Oz13  \n",
       "11 too large          0.000001149468 Oz13  \n",
       "12 overrepresented    0.184060958601 Oz13  \n",
       "13 duplicated         0.746469175613 Oz13  \n",
       "14 random breaks      0.053636871044 Oz13  \n",
       "15 valid pairs        0.221349489967 Oz13  \n",
       "16 dangling ends r1   0.050000000000 Mz13  \n",
       "17 ligated r1         0.148030000000 Mz13  \n",
       "18 dangling ends r2   0.050000000000 Mz13  \n",
       "19 ligated r2         0.144440000000 Mz13  \n",
       "20 self circle        0.000245294992 Mz13  \n",
       "21 dangling ends      0.424299467028 Mz13  \n",
       "22 error              0.000266190492 Mz13  \n",
       "23 extra dangling     0.337880224272 Mz13  \n",
       "24 too close from res 0.292011877729 Mz13  \n",
       "25 too short          0.043126493631 Mz13  \n",
       "26 too large          0.000000454250 Mz13  \n",
       "27 overrepresented    0.202854870309 Mz13  \n",
       "28 duplicated         0.538828153142 Mz13  \n",
       "29 random breaks      0.063488703484 Mz13  \n",
       "30 valid pairs        0.383513723573 Mz13  \n",
       "31 dangling ends r1   0.020000000000 Mz17  \n",
       "32 ligated r1         0.078500000000 Mz17  \n",
       "33 dangling ends r2   0.020000000000 Mz17  \n",
       "34 ligated r2         0.073100000000 Mz17  \n",
       "35 self circle        0.000283001150 Mz17  \n",
       "36 dangling ends      0.502507776100 Mz17  \n",
       "37 error              0.000377549261 Mz17  \n",
       "38 extra dangling     0.332689720112 Mz17  \n",
       "39 too close from res 0.142757357387 Mz17  \n",
       "40 too short          0.028925290269 Mz17  \n",
       "41 too large          0.000000000000 Mz17  \n",
       "42 overrepresented    0.041307876951 Mz17  \n",
       "43 duplicated         0.166845900728 Mz17  \n",
       "44 random breaks      0.113966492664 Mz17  \n",
       "45 valid pairs        0.717090182176 Mz17  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "all"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Saving tabular data for Figure 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Background data\n",
    "write.table(percs,\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/Figure1_background_data_316_samples.tsv\",\n",
    "            sep=\"\\t\",row.names=FALSE, quote=FALSE)\n",
    "# Observed 3 points data.\n",
    "write.table(all,\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/Figure1_points_3_samples_data.tsv\",\n",
    "            sep=\"\\t\",row.names=FALSE, quote=FALSE)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "adding..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "No id variables; using all as measure variables\n",
      "\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAIAAAByhViMAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nOzdZ3hUZcLG8efMTCaTSU+GhITeg9I7hBYjUgKIyLpZpSNiX6SsbRFlldVL\nUGm6UixYkBb1BWNEBMRQZGkBlBKQpQTSCSGZTKad98NgTALCzJBwJpP/7wPXnGfOnHMjJrnz\nnCbJsiwAAABQ86mUDgAAAICqQbEDAADwEhQ7AAAAL0GxAwAA8BIUOwAAAC9BsQMAAPASFDsA\nAAAvQbEDAADwEjW42D366KNhYWGnT59WOggAAIBHqMHFrri4+NKlSzabTekgAAAAHqEGFzsA\nAACUR7EDAADwEhQ7AAAAL0GxAwAA8BIUOwAAAC9BsQMAAPASFDsAAAAvQbEDAADwEhQ7AAAA\nL0GxAwAA8BKa27Mbm802duzYxYsXh4aGOkbsdvuqVau2bNlis9l69+49YcIEtVp9g3EAAADc\n2O0odmaz+Ysvvrhy5Ur5wTVr1nzzzTdPPvmkRqNZsmSJSqWaOHHiDcYBAABwY9Ve7DZs2PDh\nhx9ardbyg1arNTk5ecyYMb169RJClJaWLlmy5MEHH9RoNNcd1+l01Z0TAACgpqv2Yte3b992\n7dqdPXv2zTffLBs8f/58QUFBp06dHIudOnUyGo2nTp3y9/e/7vidd95Z3TkBAABqumovdsHB\nwcHBwRaLpfxgfn6+ECI8PNyx6O/vr9PpCgoKSktLrzte9sHs7OxDhw45XhcWFmo0t+kcQQAA\nAM+nTDEqLCz08fEpX8v0ev3ly5ctFst1x8sWf/nll+eee65s0c/P7/YEBgAA8HzKFLuAgACL\nxWKz2cqueDUajQEBAXq9/rrjZR9s1arVCy+84Hj9wQcflM3eAQAAQJli57jpSX5+fp06dYQQ\nJSUlJpMpNDRUr9dfd7zsg9HR0SNHjnS8/vLLLysd4QUAAKjNlLlBcePGjYODgw8ePOhYPHjw\noJ+fX4sWLf5sXJGQAAAANYsyM3ZqtXrIkCGffvpp3bp1VSrVihUr7rnnHsc9Tf5sHAAAADem\n2FWliYmJFovlnXfesdvtvXv3Hj9+/I3HAQAAcGOSLMtKZ3DTmDFjPv300/T09ObNmyudBQAA\nQHnKnGMHAACAKkexAwAA8BIUOwAAAC9BsQMAAPASFDsAAAAvQbEDAADX99tvvykdAa6h2AEA\ngOs7duyY0hHgGoodAACAl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AAICXoNgBAAB4CYodAACA\nl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AAICXoNgBAAB4CYodAACAl6DYAQAAeAmKHQAAgJeg\n2AEAAHgJih0AAICXoNgBAAB4CYodAACAl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AVGA2m3/+\n+efCwkKlgwCAyyh2AFDBihUrnnjiialTpyodBABcRrEDgArsdnvZnwBQs2iUDgAAnmXy5Mld\nu3aNiYlROggAuIwZOwCoQKvVRkdHBwQEKB0EAFxGsQOAClauXDlixAjOsQNQE1HsAKCCvLw8\nIURubq7SQQDAZZxjBwAVTJkypW3btu3atVM6CAC4jGIHABXo9fq7775b6RQA4A4OxQIAAHgJ\nih0AAICXoNgBAAB4CYodAACAl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AAICXoNgBAAB4CYod\nAACAl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AAICXoNgBAAB4CYodAFRw5syZ1157befOnUoH\nAQCXUewAoIKVK1d++eWXc+fOVToIALiMYgcAFfTq1SsoKOiee+5ROggAuEyjdAAA8Czx8fHx\n8fFKpwAAdzBjBwAA4CUodgAAAF6CYgcAAOAlKHYAAABegmIHAADgJSh2AAAAXoJiBwAA4CUo\ndgAAAF6CYgcAAOAlKHYAAABegmIHAADgJSh2AAAAXoJiBwAA4CUodgAAAF6CYgcAAOAlKHYA\nAABegmIHAADgJSh2AAAAXoJiBwAA4CUodgAAAF6CYgcAAOAlKHYAAABegmIHABXk5uZ+9NFH\nx48fVzoIALiMYgcAFaxYsWLx4sXPP/+80kEAwGUUOwCoQJZlIYRWq1U6CAC4jGIHABWYzWYh\nRGFhodJBAMBlGqUDAIBnGTt2rNlsHjRokNJBAMBlFDsAqKBx48avvvqq0ikAwB0cigUAAPAS\nFDsAAAAvQbEDAADwEhQ7AAAAL0GxAwAA8BIUOwCoYMuWLXFxcQsXLlQ6CAC4jGIHABWkpqZe\nuXLlu+++UzoIALiM+9gBQAVjx45Vq9VxcXFKBwEAl1HsAKCCxo0bv/jii0qnAAB3cCgWAADA\nS1DsAAAAvATFDgAAwEtQ7AAAALwExQ4AAMBLUOwAAAC8BMUOAADAS1DsAAAAvATFDgAAwEtQ\n7AAAALwEjxQDAKC2eOn1N5u1inF+fbPZ/PGXG5xf/+L5c8899bjruVBlKHYAANQWDZs069an\nf/Vtf8cP31ffxuEMDsUCAAB4CYodAACAl6DYAQAAeAmKHQAAgJeg2AEAAHgJih0AVGA0Gr//\n/vvs7GylgwCAy7jdCQBU8P7773/22WctW7b8/PPPlc4CAK5hxg4AKggPDxdChIWFKR0EAFxG\nsQOA65AkSekIAOAyih0AVJCXl1f2JwDULJxjBwAVTJkypU2bNu3bt1c6CAC4jGIHABXo9foB\nAwYonQIA3MGhWAAAAC9BsQMAAPASFDsAAAAvQbEDgAp48gSAmouLJwCgguXLl69cuTImJubT\nTz9VOgsAuIYZOwCoQKfTCSH0er3SQQDAZczYAUAFDz/8cO/evZs1a6Z0EABwGcUOACpQqVR3\n3HGH0ikAwB0cigUAAPASFDsAAAAvQbEDAADwEhQ7AAAAL0GxAwAA8BIUOwAAAC9BsQMAAPAS\nFDsAAAAvQbEDAADwEhQ7AAAAL0GxAwAA8BIUOwAAAC+hUWSvO3bseOONNyoN3nXXXVOnTk1K\nSvroo4/KBlUq1VdffXVbwwEAANRMyhS7O+64Y9asWWWLVqt10aJF7dq1E0JkZ2e3b99++PDh\njrckSVIkIQAAQI2jTLELDQ3t2rVr2eKXX37ZrFmzuLg4IUR2dnZMTEz5dwEAAOAMZYpdeTk5\nOevWrZs/f75jci4rK6tHjx5KhwIA4DpkWd5zek/a+bQLBReC/IKa1WnWv1X/YL9gpXM569jh\ntEt5udW3/cuX8qtv43CG8sXu888/79OnT926dYUQsixnZ2fv3bt39erVJpOpdevWEydOjI6O\nLlv5ypUr58+fd7w2mUxqtVqZ0ACA2scu25duX/ru1ncjgyIDdAFmm3n1f1dvP7F96oCpDUIb\nKJ3OKWGGOk1atKy+7f92/Fj1bRzOULjYXbhwITU19b333nMsFhYWlpaWWq3Wp59+2mazrV69\n+sUXX1yyZIler3essHfv3pkzZ5Z9vGwcAIDqtvXY1vd/fL9zo85ajdYxUj+4ftq5tOU/LZ89\nbLZKqgE3moiIir6zY+fq235BPjN2ClO42H311Vddu3Y1GAyOxYCAgI8//jg4OFilUgkhWrRo\nMWHChF27dsXHxztWqFev3siRIx2vt27darFYFIkNAKiFUk+mNghrUNbqhBBCEk0MTZL2Jz3U\n/aGWkdU4EwY4ScliZzabt2/fPm3atLIRtVodGhpathgYGBgREZGb+8fZAC1btnzhhRccr48e\nPWoymW5bWgBAbSbLclZhVqAusNK4SqXy1/pfvHyRYgdPoOS88d69e+12e6dOncpG0tLSZsyY\nceXKFcei0WjMyclp0KBmnLgAAPBikiSpJJVdtl/7ll22a1TKn7MOCGVn7A4ePNiqVSuN5o8M\nMTExWVlZ8+bNGzZsmK+v75o1ayIiIrp166ZgSAAAHJrWaZqelR6kCyo/WGotLS4tbmJoolQq\noDwlZ+wOHz7cqlWr8iO+vr7z58/39fVdsGDBm2++GRoa+uqrr5ZvfgAAKOWeO+7JKMjIL/7j\n+gCz1Xw88/jT8U9Hh0Tf4IPAbaNkZyq7GLa8iIiIsrPoAADwHG3qtVn0t0Vr963dc3pPoC7Q\nbDVfLrk8IXbC2F5jlY4GXMVkGAAAzoqLievQsMORjCMZlzKC/IKaRzRvEdGCp1/Cc1DsAABw\nQag+tE+LPkqnAK6vBtxNEQAAAM6g2AEAAHgJih0AAICXoNgBAAB4CYodAACAl6DYAQAAeAmK\nHQAAcI9dI2er7QVKx8AfuI8dAABwjUo2BpZ+5Wfb72/+UQhxRTu4WBtX7NNf6Vyg2AEAAFdI\nsjnMtCigdJNF3dio7SHZ7VrbbwHFP+TqpwsRqHS62o5DsQAAwAV6S2qg6ZtSzR12KUDIKlnS\n2FThpeo7DMWv+2vzlE5X21HsAACAC3S2w1Z13UoVwi7p7VJQmO8ZpVLBgWIHAABcoJKvyEJ7\n7bgs6bRq4+3Pg/IodgAAwAV24S8J87Xjkig12/W3Pw/Ko9gBAAAXlGraaOQsIdnLD6pko8p+\nOb+koVKp4ECxAwAALijW9r2iHeJr+VUlFwkhC2FT2/N8bb/m6v9RbDEona6243YnAADABbLw\nzdc9ZVXV01n26q27hCwX+8QV6MYVaeOF+EHpdLUdxQ4AALjGLukLfEdLvolqe44saW1SKMcA\nPQTFDgAAuEMWGqsqSukUqIB+DQAVmM3mn3/+ubCwUOkgAOAyih0Az5Kfn79nzx6r1apUgKVL\nlz7xxBNPP/20UgEAwG0UOwAeRJblxx577PHHH1+8eLFSGTQaTdmfAFCz8J0LgAeRJEmr1Qoh\nHH8qYuLEiZ06dWrTpo1SAQDAbRQ7AJ7l/fffP3v2bKtWrZQKoNVqu3XrptTeAeBWUOwAeBa9\nXh8TE6N0CgCokTjHDgAAwEtQ7ACggr17944ePfqLL75QOggAuIxiB8CzWK3Ws2fPKhjgm2++\nOXbs2Mcff6xgBgBwD8UOgGd57LHHRo4cuXTpUqUC9OvXLyoqaujQoUoFAAC3UewAeBBZlvPy\n8oQQRUVFSmXYsmXLxYsXN27cqFQAAHAbV8UC8CCSJC1evPjo0aOxsbFKZcjJyRFClJSUKBUA\nANzGjB0AzxIdHR0fH6/T6ZQKkJCQIISIi4tTKgAAuI0ZOwCoYOjQof369QsMDFQ6CAC4jBk7\nAKhMlmWlIwCAOyh2AFBZamqq0hEAwB0UOwAAAC9BsQMAAPASFDsAAAAvQbEDAADwEhQ7AJ7F\naDQeO3aM61IBwA0UOwCeZcqUKaNHj162bJnSQQCg5qHYAfAgsiw7nhJbUFCgdBYAqHkodgA8\ni81mK/sTAOASih0ADyJJUkREhBAiMjJS6SwAUPPwrFgAnuXdd989e/Zs8+bNlQ4CADUPM3YA\nPItWq6XVAYB7KHYAAABegmIHAADgJSh2AAAAXoJiBwAA4CUodgAAAF6CYgfAs/CsWABwG8UO\ngGfhWbEA4DaKHQAPIsuyyWQSQjieGAsAcAlPngDgQSRJWrBgwdGjR2NjY5XOAgA1D8UOgGeJ\njo6Ojo5WOgUA1EgcigUAoLLCwsJ58+YlJSUpHQRwDTN2AABUlpKS8sUXXwgh+vbtazAYlI5T\nZbIuZGxJ3uD8+kajUa/XO7++xWx2PRSqEsUOAIDKunTpYjAYWrRoER4ernSWqvTiM0+5tH5y\ncvKQIUOqKQyqA8UOAIDKmjZtmpKSonQKwGWcYwcAAOAlKHYAAABegmIHAADgJSh2AABUZrPZ\ntm/ffvLkSaWDAK7h4gkAACrbtGnTrFmz9Hr9N998ExgYqHQcwFnM2AEAUJnjLieBgYG+vr5K\nZwFcwIwdAACVdevWLTk5OSgoSKvVKp0FcAHFDgCA64iIiFA6AuAyDsUCAAB4CYodAACAl6DY\nAQAAeAmKHQAAgJeg2AEAAHgJih0AAICXoNgBAAB4CYodAACVnTlz5oEHHnjhhRdkWVY6C+AC\nih0AAJXt2bPnt99+27RpU15entJZABfw5AkAACobOHDgmTNnmjZtajAYlM4CuIAZOwAAKtNo\nNEFBQQEBAUoHAVzDjB0AAJV9++23S5cu1Wg0Xbp0CQsLUzoO4CyKHQAAlcXExOh0ugYNGoSE\nhCidBXABh2IBAKgsMzPTZDJlZWWZTCalswAuoNgBAFCZWq0u+xOoQTgUCwBAZX379l2+fHn9\n+vX1er3SWQAXUOwAAKhMpVJ16NBB6RSAyzgUCwAA4CUodgAAAF6CYgcAAOAlKHYAAABegmIH\nAADgJSh2AAAAXoJiB3iQw4cPjxw58s0331Q6CACgRqLYAR7kp59+Onv27FdffWU2m5XOohij\n0Xjs2DFZlpUOglqtsLBw3rx5SUlJSgcBXMMNigEPMmLEiPPnz3ft2lWr1SqdRTFTpkw5evTo\nI4888sgjjyidBbXXd99998UXXwgh+vbtazAYlI4DOIsZO8CDREdHz50797777lM6iGJkWXY8\nc724uFjpLKjVIiMjVSpVWFhYaGio0llwcyaT6Z133unWrVt4eHhISEjnzp1nz559+fLl25+k\nf//+PXr0uP37LcOMHQAPIknSggULTpw4ERsbq3QW1Gq5ubl2uz0/P//y5cthYWFKx8GNWCyW\n+Pj4nTt3xsfHP/nkk7Is79u3b86cOZ988smBAweCg4OVDnhbUewAeJbo6Ojo6GilU6C2GzRo\nUG5ubuPGjWl1nm/ZsmU7d+6cN2/e9OnTywa//PLLkSNHzp49+5133lEw2+1HsQMAoDK9Xs9Z\nnjXFrl27hBCV/r3uu+++tm3bpqamKhRKMZxjBwAAajDHRfRpaWmVxlNSUtauXVu2+PXXX/fp\n06dOnToBAQFt27ZduHBh2dX3gwYNGjFixN69ewcMGBASEtKpU6f169ebzeZp06Y1b948MDBw\n8ODBZ86cEUJYrVZJkt57770VK1b07ds3KCioR48eH3zwwZ9lO3v27IMPPti4cePAwMBevXqt\nX7++6v/+FVHsAABADfbXv/5VCBEXFzdmzJgff/yx7HZR0dHRTZo0cbxevnz5iBEjSktLn3nm\nmccff1ylUv39739fuXJl2UaOHj06cuTIPn36zJkzJycnJzExMTY2du/evTNmzHjooYdSUlKe\neOKJspVXrFjx+OOPt2zZ8qmnnioqKpo0adKcOXOuDXbs2LEOHTr8+OOPiYmJ06dPLyoqGjVq\n1IIFC6rxv4UQQq6xRo8eLYRIT09XOggAb/PNN98oHQHwCDXla2HlypWNGzd2FBt/f/+EhIRF\nixbl5+eXrTBgwIDo6OiSkhLHoslkCgwMHDt2rGNx4MCBQoitW7c6Fr/66ishRMeOHS0Wi2Ok\nf//+oaGhsixbLBbHXlJSUhxvFRcX9+zZ08/P7+LFi7Is9+vXr3v37o63hg0b1rBhw7y8PMei\n4yIPnU5XPliVY8YOAIDKbDbb9u3bT548qXQQOGXMmDGnTp3av3//m2++2bt37x9++OGpp55q\n2LDhxo0bHSusX7/++PHjOp3OsZifn2+1WktLS8u2EBUV1b9/f8fr1q1bCyESExM1Gk3ZiNFo\nLFu5Z8+eji4ohNDr9bNmzSopKfn+++/LRyouLt6wYcP48eP1er3JZDKZTFar9eGHHzaZTI6T\nAqsJxQ4AgMo2bdo0bdq0iRMnXrlyRekscIpKperYseOMGTNSUlLOnTv38ssv22y2hx56KC8v\nTwgRGBj422+/vffee48++mj37t0bNGhQUlJS/uPlL39WqVTXHSnTrl278osdO3YUQpw6dar8\nYHp6uhBizpw5fuX87W9/E0Lk5uZW0V/6OrgqFgCAyhw3P9PpdLX5MTA1wuXLl8eNGzdixIjx\n48eXDRoMhtmzZ2s0mn/+8587duwYPnz43Llz//nPfzZo0GDEiBEzZ87s1q1bnz59qiqDWq0W\nQlit1vKDjsVZs2YNGTKk0vplZ/5VB4odAACV9erVa82aNeHh4b6+vkpnwY0EBgampKQUFxeX\nL3YOer1eCBEaGnru3LkXX3xx0qRJy5YtkyTJ8a7NZnN7p4cOHSq/eODAASFEixYtyg86Fu12\ne/kHURw5cmT37t133HGH27u+KQ7FAgBwHU2bNq1tDy2oiVQqVWJi4ubNm9966y273V42npub\n++677wYHB3fp0uXSpUtCiLZt25a1ut27d2dkZMi/3+7EVbt27fruu+8cr0tKSubMmePr6xsf\nH19+neDg4H79+r3//vvHjx93jJhMpvHjx7/yyisBAQHu7dcZzNgBAIAa7I033ti1a9f06dM/\n/PDDTp06GQyGjIyM5OTkK1eufP31135+fjExMQ0bNnzttdfy8/Nbtmy5Z8+eL774Iioqavfu\n3SkpKYMGDXJ1j/Xq1Rs+fPi4ceMiIiK+/vrrI0eOzJ49u379+pVWe+utt/r06RMbG/vXv/61\nXr16a9euPXjw4OrVqyudsVe1mLEDAAA1WGRk5N69e+fMmaPRaL788sv//Oc/aWlp995774ED\nB4YPHy6E0Gq1ycnJHTt2fOedd5577rmcnJy9e/fOmzfPZDLNnz/fjT2OHj166dKl+/fvX7Bg\ngU6nW7Zs2ezZs69drVOnTgcOHOjTp8/XX3/9xhtvBAQEJCcnP/DAA7f6F74hye15SMWNGTPm\n008/TU9Pb968udJZAHiV5OTka893BmohvhYqsVqtPj4+zz777Ouvv650lutjxg4AAMBLUOwA\nAAC8BMUOAIDKbDbbpk2bjh07pnQQwDVcFQsAQGWbNm2aNWuWXq/fuHFjUFCQ0nHgKTQajYdf\nnMCMHQAAlYWHhwshAgMDy54uCtQIzNgBAFBZ586dZ8+e3bRpU7WslWXx+31tAU9HsQMAoLKU\nL7ese2VvhJTf01C/Xh9tQD1VgwE+QU04zAVPR7EDAKCC0kuy/cdmMSV6u3+xoaO6tEC+fNJy\n5Zy91UPa0Bi10umAG6HYAQBQwbnvrZqLkd2GhfhqtWq1WvgKH391cYb9t68tnVuqOTsdnszZ\n/z0TExN//fXXa8e3bNkyZcqUKo0EZVy4cGHSpElz585VOggAKEm2i7xfrPooSe/np1b/MT+n\nj1Sd32K5ctZ+g88CirtJscv93erVq9PT03Mrys7O/vbbbz/55JPbkxXVKjU1NS0tLSkpKScn\nR+ksAKAYm0lc3GE1ll5JTk7esWPHH29IQq2VzIUefasL4CaHYuvUqVP2esSIEdddJy4urioT\nQSF33XXXrl27mjVrVv4fHQBqG5VWjuyhOXE4z/S7sjue2K2y2o/rY+HRblLs3n77bceLZ555\n5oknnmjevHmlFXx8fIYNG1Yt0XB7GQyGsn9uAKi1VBoppJlK2uVbabz0kly3h09gI4odPNpN\nit3UqVMdL7766quHH364Q4cO1R8JAAAlNRjgc/FIcH5Bc98IWafTyUKU5stX/mdr9ZBWraXY\nwaM5e1Xstm3bqjMGAACeIqC+qsuTYaHRAWdSLLlpNrtFju6jiRnjFxXLrSTg6Zz9f7SgoGDa\ntGmbN282Go3Xvpubm1ulqQAAUFJQE1XHabrW432NWXZtsKSPlFQa5upQAzhb7KZPn/7hhx92\n7969Xbt2KhX38AEAeLlffv3ljjvu0IVzR2LUJM4Wuw0bNvz1r39dtWqVxAPzAAC1wLlz51q3\nbs1PPdQszs69FRUV3X333fz/DQAA4LGcLXa9evU6cOBAtUYBAACoWnaLyPqv7bevzCfXWM7/\nYDXlVsEtpuvWrStJUqUnctnt9nr16kmSdOTIkRt/3Gq1GgyGzMzMspGMjIz777/fYDBERkZO\nmDAhLy/P7WzOFrtFixYlJSUtWrTIbDa7vTMAAIDbxpRnT1to2vNKycnVlt82mA+9a0pbaMra\nY62Sja9bt6784u7duy9evHjzSCbTrFmzylc3WZYTExP/97//ffrppytXrtyxY8etPKzV2XPs\nnn/++YYNGz799NMzZsxo2LChr2+FOzfetJwCAADcTrJNHPvEnLXXamirFr+fSmbKk3f/syTu\nff+gJrd0JWhsbOzatWtfeumlspH169fHxsampqbe4FMLFy6cOXNmpTmyU6dOpaamHjhwwHG3\n4FdfffWhhx6yWCw+Pj5uBHP2b2UymUJDQwcOHBgXF9esWbP6FbmxYwAAgOpTkG47+501qPEf\nrU4IoQ2W/CJVF1Itt7jx4cOHHzt27Pjx445FWZaTkpLuv/9+x2JqaqpU0csvvyyESExM3Ldv\n36pVq8pvqqioKCEhISYmxrEYHBwsy7LbB0idnbFLSUlxbwcAXHLx4kWDweDeL2oAgDLGi7I2\n4DqXfWqDpOILt3qmXVhYWHx8/Lp161588UUhxP79+3NycgYOHOh4t3v37jk5OY7XGzdunDJl\nSkJCghAiIiIiIiLCZDKV31SHDh02btwohLDZbGfPnl24cOGAAQP8/f3dC+baPOSlS5eSk5NX\nrlyZlZVVWFhot9vd2yuA60pOTh42bNi4ceNkuQpO7wWA2sxu/5PvpLIQVfEtdtSoUWvXrnW8\nXr9+/dChQ/38/ByLPj4+BoPBYDBcunRp6tSpixYt6tq16003OGDAgKZNm+7evXvx4sVup3Kh\n2L399ttRUVEJCQnjxo07ffr0mjVrGjVqVPZXAnDrHA9xyc3NtVqr5txeAKi1/KNU1mL52g5n\nKZL1davg9m0jRow4cuRIenq6LMvr168fNWpUpRWKiopGjhw5cuTIyZMnO7PBlStX7t69OyEh\noUePHllZWe6lcrbYrVu3btq0ab179y47MNyjR4/IyMgHHniAo7RAVRk1atSsWbMWLVrEoVgA\nuEUhLdX14jRXzlU4umgpkkuy5KheVfA91mAwxMXFrV+//pdffjl37tzgwYPLvyvL8qRJk7Ra\n7ZIlS258G+AzZ84cPHhQCFG/fv3u3bt/8MEHxcXFbpcrZ8+xmz9/frt27VJSUsomEtq0aZOa\nmtqjR4+5c+cOGjTIvd0DKE+v1997773Xf89+RVgzhaQVmrpC8r3+OgCA36k0Il/d10YAACAA\nSURBVGaM9qjNnLnLqg2SJJVkM8mmS/YuL/qFtKyah6OOGjVq6dKlpaWlgwcP9vf3LzuvTgjx\n1ltvbd68ed++fWXHZ//Mzp07p06dmpGRodFohBAmk8lisYSGhroXydlid/jw4ZkzZ2o0mvJH\niHQ63V/+8pc333zTjR0nJSV99NFHZYsqleqrr74SQtjt9lWrVm3ZssVms/Xu3XvChAlqNc/p\nQ+1mKxBXPhOmQ6Jkh5Bl4R8n9HeLwOGuniMLALWNvq6q/d91Ub2sRRfsNpPsV0dlaK8OqF9l\n3zxHjBjx+OOPZ2Zmzps3r/z4tm3bnn322Y8//jggIMBxjo2Pj09wcPB1NxIfH19UVPTwww8/\n+uijFovl9ddfj46O7tevn3uRnC124eHhlS7icMjIyAgMDHRjx9nZ2e3btx8+fLhjsWyWcs2a\nNd98882TTz6p0WiWLFmiUqkmTpzoxvYBLyGbRP58YdwqtE2FX6wQsrDlitw5wn5ZBI9TOhwA\neDqNn4ju62zbcVVkZGTfvn137tzpuOi1zMaNG2022+jRo8tGEhISHJe+XisiIuLbb7997bXX\nhg4dqlKp+vTps3nz5j9rgTfl7F+1Z8+en3zyyT/+8Y/yM4rHjx9fs2aNe6UyOzs7Jiam0kUi\nVqs1OTl5zJgxvXr1EkKUlpYuWbLkwQcf1Ol0buyiZjlw4MCqVavuvffe2NhYpbPAkxRvFsWb\nhK69uHojJkmogoWurch/U+j7CZ/GyqarDkaj8ezZs61ateLh1AA8U/mngW3durXsdePGjR33\nNJg3b16lObzyunTpUunWB3379u3bt2+VZHN2NvKNN94oLi7u2LHjv//9byFEcnLyiy++GBsb\nazKZXn/9dTd2nJWVVadOnUqD58+fLygo6NSpk2OxU6dORqPx1KlTbmy/xvnwww+3bNmyYMEC\npYPAw5h/FZoIISpWHEknVMGi1Duf+PLYY4+NHj162bJlSgcBgJrH2Rm7Ro0a7dq1a+bMmXPm\nzBFC/Otf/xJC9O/ff/78+S1atHB1r7IsZ2dn7927d/Xq1SaTqXXr1hMnToyOjs7PzxdChIeH\nO1bz9/fX6XQFBQWubr8mGjRo0P/+97+yY9PAVfYiIWmvMy75CnvRbU9T7WRZvnLlihCiqMgL\n/3YAUN1cOOocExOzYcOG4uLi9PR0q9XaokULtw8AFxYWlpaWWq3Wp59+2mazrV69+sUXX1yy\nZElhYaGPj4/jqhAHvV5/+fLlssWtW7fOnDmzbNG90/s805AhQ4YMGaJ0CngeVbCQr3N6q5BL\nhNrNa6Y8mSRJixcvPnr0KOckAIAbXCh22dnZK1eu7NKlS//+/YUQCxcuNBqNU6ZMceOK3ICA\ngI8//jg4OFilUgkhWrRoMWHChF27dgUHB1ssFpvNVnYlrNFoDAgIKPtgYGBg69atHa9Pnz5t\ns9lc3TVQw+g6iMJPhaZuha9W+xXh20H4tlMuVjWKjo6Ojo5WOgUA1EjOnmOXmZnZoUOHmTNn\nHj161DGSnp7+/PPPt2vX7syZM67uVa1Wh4aGOlqdECIwMDAiIiI3N9fRER0HZIUQJSUlJpOp\nfHHs0qXLJ79r0qSJ0Wh0dddADaPvJwIThemwsF0SslnIJmHLFKYjwn+Q0EQpHQ4A4FmcLXYz\nZ840Go1btmx59NFHHSOLFi366aefrly58vzzz7u617S0tBkzZjjOpBFCGI3GnJycBg0aNG7c\nODg42HH/ZSHEwYMH/fz83DiHD/AiGhE2VYTNFD6NRMnPQvITvp1E1Hsi8E/uYwwAqMWcPRSb\nmpo6efLkuLi48oO9e/eePHly2UPGnBcTE5OVlTVv3rxhw4b5+vquWbMmIiKiW7duarV6yJAh\nn376ad26dVUq1YoVK+65557acK8T4EYkXxH0FxH0F2F4SUg+QuIrAgBwfc4Wu8uXL/v7+187\nrtfr3Tge6uvrO3/+/OXLly9YsECSpA4dOkybNs1xzURiYqLFYnnnnXfsdnvv3r3Hjx/v6sZR\nc23bti07O/u+++7jSanXp/Keq4UAANXB2WLXpUuXtWvXzpgxo/ylDMXFxevWrevYsaMbO46I\niHjhhReuHZckaezYsWPHjnVjm6jRcnJyZsyYIYQIDAys9ChlAADgDGeL3csvv9yvX7+ePXv+\n/e9/v/POO9Vq9dGjR99+++2jR48uXry4WiOilggNDW3WrFl2dnbLli2VzgJv88bi9+rWq+/8\n+maz+eMvNzi//plTp5596jFfX1/XowFAVXK22PXq1SspKemZZ56ZPHly2WB0dPRnn31W6cQ7\nwD0ajeaLL76w2Wzlb2QIVInwOpHd+vSvvu2r1RruvgTAE7jwE3TYsGGDBg3av39/enq6xWJp\n0aJFp06d9Hp99YVDbSNJEq0OAAC3OftDtEGDBoMGDVq2bFn37t27d+9erZkAAADgBmfvY1en\nTp2ffvpJluVqTQMAAAC3OVvsVq9eLUnSI488wpO5AQAAPJOzh2Kff/75+vXrL1++fMWKFY0b\nNw4LCyv/7t69e6shGwAAAFzg7IxdQUGBLMvx8fF33XVX06ZNQyqq1ogAAACeo27dupIk/frr\nr+UH7XZ7vXr1JEk6cuTIjT9utVoNBkNmZqZjcd26ddI13H5Ag7Mzdps3b3ZvBwAAAIqxZovi\nTcJ6TsgmoY4Qfj2FrlOVbHjdunUvvfRS2eLu3bsvXrx400+ZTKZXXnklLy+vbCQ2NnbDhj9u\nnGmxWCZNmuT2veRcu7XEpUuXdu3alZubO3DgQD8/v4CAAJXK2Tk/AACA28p8TFx6V5QeFKoQ\nIdRCNopL74vwf4jgh4SQbmXDsbGxa9euLV/s1q9fHxsbm5qaeoNPLVy4cObMmWazufxgVFTU\n0KFDyxbnzZvXuXNntx/B5UIte/vtt6OiohISEsaNG3f69Ok1a9Y0atRo7dq17u0YAACgGskm\nUbBcWH4T2tZCEyU0EcKnsfDrKPLfFCU/3+K2hw8ffuzYsePHj1/dlSwnJSXdf//9jsXU1NRK\nh1ZffvllIURiYuK+fftWrVr1Z5s9e/bsv//97/fff1+S3Oydzha7devWTZs2rXfv3mVpevTo\nERkZ+cADD6SkpLi3bwAAgOpSelgYtwhNvQqDkk5o6omSnbe47bCwsPj4+HXr1jkW9+/fn5OT\nM3DgQMdi9+7dc3734YcfarXahIQEIURERESbNm2aN2/+Z5udPXt2YmJi06ZN3Q7mbLGbP39+\nu3btUlJSRowY4Rhp06ZNampq+/bt586d6/buAQAAqoU1U6gCrjOuChDWzFvf/KhRo8qOW65f\nv37o0KF+fn6ORR8fH4PBYDAYLl26NHXq1EWLFnXt2vWmG0xPT1+zZs1zzz13K6mcLXaHDx8e\nOXJkpcc96XS6v/zlL4cOHbqVBABQ3qJFi0aMGHHj81QA4OYkHyFf9yHOdiH53PrmR4wYceTI\nkfT0dFmW169fP2rUqEorFBUVjRw5cuTIkZMnT3Zmg/Pnzx86dGiDBg1uJZWzxS48PNxkMl07\nnpGRERgYeCsJAKCMLMv/93//d/78+e+//17pLABqOJ8mwl4oZGvlcVue8HH/WGcZg8EQFxe3\nfv36X3755dy5c4MHDy7/rizLkyZN0mq1S5YsceaEOZPJ9MUXXzz00EO3mMrZYtezZ89PPvnk\n0qVL5QePHz++Zs2abt263WIIAHCQJOmZZ54ZMGDA6NGjlc4CoIbTthTB44T5WIVuZ8sU2jtF\nwMAq2YPjaOz69esHDx7s7+9f/q233npr8+bN69evLzs+e2PffPONzWYbNGjQLUZy9nYnb7zx\nRocOHTp27Dhu3DghRHJy8oYNG95//32TyfT666/fYgg4mM1mrVardAolGY1Go9FoMBiUDgIl\nDRkyZMiQIUqnAOAFJBEyRQghLn8sVIFCqIW9WOj7iZCJQlO/SnYwYsSIxx9/PDMzc968eeXH\nt23b9uyzz3788ccBAQG5ublCCB8fn+Dg4Bts6vvvv+/evfut1wBnZ+waNWq0a9eutm3bzpkz\nRwjxr3/9a+7cuW3btt2+fXuLFi1uMQSEECtXruzVq9drr72mdBDFmEyme++9d9CgQZy1Wcvl\n5+fv2bPHar3m6AkAuEoVKMKmi+g1wvCyCJ8p6r4rIl4Tvm2ravORkZF9+/bNzc11XPRaZuPG\njTabbfTo0XV+d9NjrNu2bevRo8etR3LhBsUxMTEbNmwoLi5OT0+3Wq0tWrS4cfeES9LS0oQQ\nR48eVTqIYkpKSoxGoxDC8ctN7ZSdnf3BBx906tTpnnvuUTqLYp588skTJ06MHz/+ySefVDoL\nAC8gCd/Wwrd1FW6x7GlgQoitW7eWvW7cuLEsy0KIefPmVZrDK69Lly6O1co7duxYlWS7ebEz\nGo0nTpzIyclp0aJFw4YN/f39O3ToUCX7RnnTp09v06ZN//79lQ6imNDQ0P/85z/Z2dm1+T/C\n+vXr161bt3Hjxv79+9fO4/KyLNvtdiGE408AgEtuVOxkWZ4/f/4rr7xSVFTkGOnVq9eHH37Y\nsmXL25KtdomOjp4wYYLSKRSWl5eXl5dns9lq7aPq+vTp8/333/fo0aN2tjohhCRJ77777smT\nJzt1qpqHOQJArXKjYvf555/PnDkzJCRkypQp4eHhu3bt2rp165AhQ9LS0ipd+gHcupycnBkz\nZgghAgICau25823atElKSlI6hcLCwsK41h4A3HOjYrdgwYLQ0NBDhw7Vr3/14pEXXnjh3//+\n94YNGxITE29LPNQioaGhzZo1y87OZkoYAAD33OiA19GjR++7776yVieEeOKJJ0TVnd8HlKfR\naO69996EhITo6GilswAAUCPdaMauqKiobt265UccP3G5DQGqQ2Zm5ltvvSWEaN68edkjiQEA\ngPNucop6pYdgOPNMDMA9ERERHTt2jIiI6NKli9JZAACokWrptYceaN26df3791+6dKnSQRSj\nUqmWLVuWnJxc/ug/AABw3k3uY5eWlvbRRx/ddHD8+PFVGapW+umnn4qKin744YdHHnlE6SwA\nAKBGukmx27hx48aNG286SLG7dY8++midOnVq7W0+AADArbtRsVu7du1ty4HWrVv/85//VDqF\nwn788cesrKz77rvPx8dH6SwAANQ8Nyp2o0aNum05gNzc3OnTpwshgoKCBg0apHQcAABqHi6e\ngKcICgqKiorSarWNGjVSOgsAADXSTc6xA24brVabmJh4/vx5ih0AAO6h2HkQi8VSm88ty8rK\nevvtt4UQMTEx9957r9JxAACoeTgU60G+//57pSMoqU6dOh06dIiIiOjUqZPSWZRktwpjplx4\n2m6+LCudRRmLFi0aNmzYDz/8oHQQAKh5mLGDp1CpVMuXL1c6hZJkuzj3gyVzh+3iTotKI9nM\nctN7fZqN0vpH1aJfwGRZ3rBhQ35+fmpqanx8vNJxAKCGqUU/MAAPdzbFkvZ2qTHbHtFZY+ig\nrtNBc3Gn7dfl5tJLtWjqTpKkF198ceTIkRMmTFA6CwDvYZftpdbSKtxgaWnpv/71r9jY2KCg\noFatWiUmJh45csSZD1qtVoPBkJmZWTaSkZFx//33GwyGyMjICRMm5OXl3UowZuwAj2C+LF/Y\naQ1urvIJuPpEZpVWBDVR5R6yntuian6/Vtl4t1O/fv369eundAoAXuJwxuENaRsyCjKsNmuI\nPqRH0x5D2w31Ud/SGe25ubl33313RkbGtGnTZs2adfHixTVr1nTu3HnZsmVjx469wQdNJtMr\nr7xSvrrJspyYmGg0Gj/99FNJkp566qkpU6asW7fO7WwUO1zl+L+qcePG99xzj9JZaqOCk7bc\nQ1ZDm8pfkn51VJdP2BWJpKCzZ882bNhQ6RQAarwfT/z4xGdPNAxrGKoPValUuUW5W45tSc9K\nf2bAM7fS7Z599tlLly6lpaVFR0c7RsaPH//cc8898cQTd999d9lgJQsXLpw5c6bZbC4/eOrU\nqdTU1AMHDnTo0EEI8eqrrz700EO3cjElh2JxVUpKytKlS1944YVbnASGe2wlQu1zna9HtY+w\nltz+OApz8ogGANxAgbFg7d61raNaNwhrEKAL0Gv1hgBD+/rtV+9dveXYFrc3e+HChQ8++GDu\n3LnlC5wkSS+//HJwcPCCBQtSU1Olil5++WUhRGJi4r59+1atWlV+a0VFRQkJCTExMY7F4OBg\nWZYrlT+XUOxwVfv27UNCQjp37hwaGqp0ltpIGyTZTNc5l85WKnwCb38cAKjxDmcc/vm3n8P8\nw8oPqlXqeiH19p/d7/ZmHb95Dhw4sNK4n59f3759Dx061L1795zfffjhh1qtNiEhQQgRERHR\npk2b5s2bl/9Uhw4dNm7cqNPpbDbb6dOnFy5cOGDAAH9/f7fjcSgWVzVr1iwlJUWj4X8JZQS3\nUEV20xgv2H3DpbJB2S4XX7Q3GVF7724IAG67ZLzkp/W7dtzPxy+vyP1jU//73//8/PwMBsO1\nbzVq1Oi///2vj4+P49309PSpU6cuWrSoa9euN93sgAEDtm7dGhYWtmfPHrezCWbsUB6tTkEa\nP6nREJ/Lp23FF2RbqbBbZXOhfOm4vV4/n3r9KHYA4DK9Vm+xWa4dN9vM/r7uT4k1bty4pKTk\nuqctnTt3Ljg42PG6qKho5MiRI0eOnDx5sjObXbly5e7duxMSEnr06JGVleV2PIod4Ckiu6r7\nvKOP6KpW+4mc/TZ9XVXMWN82U7Sa6/zCCQC4iZi6MYUlhSaLqdJ4dmF266jWbm+2TZs2Qojv\nvvuu0rjJZNq+fXufPn2EELIsT5o0SavVLlmyRJKk62zld2fOnDl48KAQon79+t27d//ggw+K\ni4tTUlLcjscMDeBBwu5Qh92htluF1Shrg270vQAAcGP1Q+v/Y9A/Fv6wsFXdVoG6QCGE1WY9\nnXe6V/Neg+4c5PZmo6OjJ06c+Pzzz8fFxUVFRTkGZVmePXv2xYsXH3vsMSHEW2+9tXnz5n37\n9vn53eRX8507d06dOjUjI8Nx0MxkMlkslls52Z1iBw9SWlp66tSpO+64Q+kgClNpBK2uah3a\nu+enzZV/va5CNpt15N3cew/wOA91f8hP67ft+LbU9FS1Sl1qLR3dY/TE2Ikh+pBb2ewbb7wR\nHx/fvn376dOnd+jQITMzc/Xq1d9+++3777/fsmXLbdu2Pfvssx9//HFAQEBubq4QwsfHp+wQ\nbSXx8fFFRUUPP/zwo48+arFYXn/99ejo6Fu5lyfFDh7EZrPl5OQonQJeqF2XbrHxA6pv+/t2\nplbfxgG4TaPWPNDlgcFtBmcUZJRaSyMCI+oG1b3xsVFnGAyGPXv2vPHGG19//fWrr74aFhbW\nuXPnffv2OZ51vnHjRpvNNnr06LL1ExISNm7ceN1NRUREfPvtt6+99trQoUNVKlWfPn02b978\nZy3QGRQ7wLOcOHGiXr16t3KtOwCgvEBdYEzdmKrdpq+v70svvfTSSy9d+9a8efPmzZv3Zx/s\n0qWLLFe4uVXfvn379u1bVcG4eALwIElJSQ8++ODYsWPt9lr3tAkAwK2j2AEexGq1CiEsFkul\n3+cAAHAGh2IBDzJq1KimTZs2btxYrVYrnQUAUPNQ7AAPolKpunTponQKAEBNxaFYAAAAL0Gx\nAwAA8BIUOwAAAC9BsQMAAPASFDtcZbfbDx48eOHCBaWDAAAAN1HscNX27dsffvjh0aNHG41G\npbMAAAB3UOxwVV5enhDCZDIpHQQAALiJ+9jhKsczD8xmc3FxsV6vVzoOAEBJxRn2rD22kp/q\n/3qh1C9Siuql0YUzGVQDUOxw1V133bV79+5mzZrVqVNH6SwAACVl7rL+PLtEFyaZiwIvFljN\nV0T2XluLRG1Yax6K4+kodriqTp06b7/9ttIpAAAKKzpv3/1SSWgrtTZIMmWr9BEqfV1hvCCn\nrzZ3nKbTBklKB8SNMK0KAAD+kP1fq1+45ChwERERjkFdHSk3zZZz0KZoNNwcxQ4AAPzBmCNr\nAq4zLacNkEqy7Lc/D1xCsQMAAH9QqYS4Xn+T7bLEKXYej2IHAAD+4F9PVXpZrjQoy8JSKAXU\npzZ4Ov6FAADAH+r20IS3VRsz/5i1k2VRdM5et7fa0J4pO0/HVbEAAOAPvqFSzBjtiVXmrD02\nnwAhZMlcKNeP08SM8VX7ckmsp6PYAQCACoKbqTtO98s7bC3JkSWV0EeqwtupVVSGmoB/JQAA\nUJnGT0R2oyTUPJxjBwAA4CUodgAAAF6CYgcAAOAlKHaABzEajV9//fXx48eVDgIAqJE4LxLw\nIKtXr16yZElYWNg333zj4+OjdBwAQA3DjN1V//nPfyZOnMhMCZTVuHFjIUSDBg00Gn7pAgC4\njB8eQghhMpmWL18uhNi0aVOrVq2UjoPaKy4ubsuWLQEBAZLEXUABAC6j2AkhhE6ne+yxx9LS\n0hISEpTOgtouKChI6QheyFhcvOenbS6sbzTq9Xrn18/MOK/RxLscCwCqGsXuqkmTJikdAUB1\neXriGJfWT05OHjJkSDWFAWqKpKSkpKSkSZMmxcXFKZ0FzuIcOwAAcB1r1qw5duzYunXrlA4C\nF1DsAADAdTz66KM9e/acOHGi0kHgAg7FAgCA6+jfv3///v2VTgHXMGMHAADgJSh2AAAAXoJi\nBwAA4CUodgAAAF6CiycA3A42s1x8QS7Nk3Xhkn89lYoH4QJANaDYAah2mbutZzdZM3dZ1T6y\nzSJF9dI0GuwT0UWtdC4A8DYcigVQvXL22/bMNpVkynU6qsPaaAwd1UUZ9l0vGHMP2ZSOBgDe\nhmIHeDS7bDdbzUqnuAWyOPu9JaChShtydUASQhcmBTZQndtkkWXXN2fLE+bjwnpRyNYqjgoA\nNR+HYgEPlXYubeOhjRkFGTa7LUQf0qtZr4R2CRpVDfuaNebYz222RHSuHFtnUJ1JsbSeoNWF\nO/3rpeWMKFwpLq8SklrINhEwRASNEroeVZwYAGqyGvZDAqglth7f+tSqpxqFNQrxC1GpVDlX\ncn44+sPJ7JN/j/+7Rl2TvmxtJiFJQrq2vKmFpJKsJqc3ZDkv8ucL8zHh10NIPkLYhfmkuDhZ\nRL4r9H2qMjEA1GQ16ScE4PUyMjKWLFnSukPr3Zbdd0bdGaK/evxSr9WH+od+/vPn7Ru0v7v1\n3cqGdIlviCRkYSsVat8K4/ZSOaKr2jdEcnZDRUmi9Behbf77skpoIoWkFlfWCb/uQtJWYWYA\nqLk4xw7wIElJSZs2bXp96ev/Pf3fslbnoFFpokOi95/Zr1Q292iDpBZ/1RZl2MqfTScLUXRe\nDmmh9vF3stjZROkxoalbeVhtEMVbheW3KgoLADUexQ5XXbhwYdKkSXPnzlU6SK1mtVqFEMJX\n6LS6a9/V+ejyivNud6Zb1nSENvxOdcFxW2mBbDWJ0gL58gl7eDt1k+FO38tONgnjtutPy6m0\nwn6lCtMCQI1GscNVqampaWlpSUlJOTk5SmepvWJiYoQQEWERVvt1Lvm02Cz+Wv/bHupW6QxS\nu6d0Lf+mDWqi0viJoCaqlg9q2z/lqwtz+jis5Cf844X92jPyZCGXClXIdT4CALUS59jhqrvv\nvvvIkSPNmjWrU6eO0llqr8GDB3fo0MGoMg5/d7jJYtL5lJu3k0VWYVbrqNbKpXOfb4jUfNSt\nnAanEr7tROknQh1UYdh6UfgPEtqmt5YOALwHxQ5XhYWFzZkzR+kUEFFRUUKI6fdMX7J1SUxk\nTIAuQAhhtVl/y/2tT4s+A9sMVDqgQgJHCvMJYfpZaBoISS/kUmHLFtbzos5rQvAECwC4imIH\neKKxPcfqffQ/nvhxx8kdapW61Fo6vtf4cb3GBemCbv5hr6QOF+HPiytfiNJfRPFWoe8rtL1F\n0ANC21LpZIDX2r59+4YNG0aNGtW9e3els8BZFDslGc3GjIKMYF1wncA6kuT0+UaoBTQqTWK3\nxCFth1y4fMFis9QPrR+qD1U6lNLUYSLkcSHswnZJqAK5xQlQ3d59992TJ08ajUaKXQ1CsVPG\nhYILn/382cnskztP7bTZbQltE4a1H6Z0KHicIL+gIL/aOkX3p1RCHa50BqBWuOeee4qKigYO\nrK1ngNRMXBWrgJwrOW9+92by4WSz1dyrWa/Y5rHp2elPfv7kgewDSqbKyXnmmWcWL16sYAYA\ngOfYsWNHZmZmamqq0kHgAoqdAjYc2rD3zN6WkS21Gq0QQiWpIoMiW9Vt9fPFn41mo1Kptm7d\n+tNPP3300Ufc7gQAIISIiIgQQkRGRiodBC7gUKwCDp8/HB0SXWkwzD9s88nNxzOPd2zYUZFU\nffr0+e6777jdCQDA4bXXXnviiSfq16+vdBC4gGL3B4vF4uPj9K3w3WWz24xmo1Z9nfO+fVQ+\nV0yK3UM/KipqxYoVSu0dAOBpVCoVra7G4VDsVY899lh8fPxtOJNArVIH+AaYLNfcQ18WZrs5\nWB9c3QEAAIC3otgJIYTJZDp06JDRaDx06NBt2F3Hhh0vFFyQRfmnoovsK9ntDO1aRba6DQEA\nAIBX4lCsEELodLrXXnvt6NGjf/vb327D7oZ3GH444/CO9B0NwhoE6ALMVnPOlZyz+WcntphY\n4RFSAAAArmDG7qr+/fs/9thjISG342niQbqgZwc9Oy52XJh/2I6TO3Q+uu5Nu382+bOYsJjb\nsHcAAOCtmLGrRikpKUbjn96+JFyE3+V3V/c7u/uqfX0sPqf2nDp16pTJdM25d3/Oz89v8ODB\nVZG0uphMpry8POfXLykpycvLy8jIcP4jQUFBgYGBrkcDAMALUeyq0fnLxtj4Ac6v37qvCysL\nIXb88L2LiW63ffv2ZZpsfnp/5z+ii2qUdjbzz96VZfuZgt/yirOtdmuo3tAkrLk9O/O+++6r\nirAAANR4FDtUr0ZNm/tX0YzaZeOlL35e/u2hJH/fAElIRrOxV4u4nr5tC4Z9igAAIABJREFU\nqmTjAAB4AYodagabzfrprv/8fGp7u/pdVCqVEELIIj3z11zb6fuK7wv35+GhAABw8cTv9uzZ\n8/nnnxcXFysdRDEXLlyYNGnS3LlzlQ5yfUcvHv7hl41N6jS/2uqEEJKICql/suDk9794+iFp\n1Dg22ZZRkHEi60R+cb4syzf/AAB4BmbshBDCZDLNnDmzuLi4qKjokUceUTqOMlJTU9PS0tLS\n0iZPnuyBTxXLuPS/IL9g6ZpfRQJ9Ak/nnlYkEpxktVrT09Nd+siFCxeOHj3q/PrBwcHR0ZUf\n01eJbBe5B2wF6TbTJdnHXwqop4rsofbxl65dc9vxbSsOr3huz3M+ah+z1Tym55hxPcdd+xhA\nAPBAFDshhPD19a1fv/7JkyebNm2qdBbF3HXXXbt27fLYZ8VabFaVpL52XFJJVrv19ueB8woL\nCzfv3tss5g7nPxJ9Z8fTBX96Rfm1Tuz4eerD42+wgmwTx1aWpq+x+BkktV7YzVJpnq3uHp/W\nk7T6iAq/LWz6ddOMNTMC7AGxzWMlSTJbzZt/3ZxdmP38kOcjAiOcjwQAiqDYCSGEJEmffPJJ\ncXFxbb5xhsFgePvtt5VO8acMAXWKS4uuHS+xlNQNrnv788Aldes1aNKi5U1Xyy/O3Zm+5ULB\nObO1NMzf0K5B1zb1Ozqz/Qtnz9x4hYwfrSfXWw1tVJLP1Sk6/7qanINW7VqpzeO+0u/TdkWl\nRf938P9aR7cuzCmUJEkIodVom0c0P3DuwFcHvnqkby2dzgdQg3CO3VUqlao2tzrPd2e9ju0a\ndMktyio/WGIuvlR6qVezXkqlQhU6nZP+3g9vrP55xeFze09mHd169Nt/rn88ad8nVXKKW/Ze\nq389UdbqhBBCEoGNVKe+MhdfsJeNncg6sePkjhC/yjcqrxtU90jGkVuPAQDVjRk71AyBfsH3\ndRmdtPeT9MxfA/2CJKEymotzr2SObT32zug7lU6HW1VqNa3b+/G5/NPNIq4+fyVEH2YIjPgk\n9b1mdWLaN+x6Kxu3W0RpgawNqHw6nUoj+eglU44cUO/qSHFpsVajvXYLWrW21FpqtVs1Kr5n\nAvBofJNCjXFHdPvw/nV2ntyaffmCXbbXDa7XvkFX6X/ZSudCFThx8Zfd6dva1u9cftBXo4sM\njt5/ZtctFjshyZIkhP0610nINkkqd9wi2C/YZDGJa6YISywlEUERtDoAno/vU7jqt99+e/zx\nx/+fvfuOj6LO/wf+ntle0za9ACmQ0JNAIISOIKJUK3pgwa7n1y+Hd56nHnp33s9TEeU4Pb5W\n0ANRsSAoUo5eAgHpJZAC6T27m+0z8/sjEQJEJJDdzyb7ev7BI5nMfubFBmbf85n5fD7JycmL\nFi3iuDY+Av1BpDFmesY9rbcUorDrEmqbqrUqHV32706r1NVYr/dXzMs5XQxftU/Qai46gMdB\nHqeoj79Q2fWK6jU2dWxBTUHr3SSSSupLpgyYcp0xAAB8AM/YQYu8vLyamprdu3fX1NSwzgIB\nRy5TiJJ4+XZREhW84vrbjx4ht5YJLvOFvjjRRY1nhD5zVKqQC9WeSq66a/BdZ6rO1Nhr7C67\nW3CbHeZjpcfGpI6Zmj71+mMAAHgbeuygxcSJE4uLixMTE/1zuhPo2uJCulkdZo/glssuKuPq\nm2qH97zh+ts39ZcNek5zbr279pAg13GSh1xmKWWmMnH6pVXjkMQhKx5eseCLBWq1evPJzTf2\nuXFC7wm3Zt5qVBuvPwYAgLehsIMWBoNh3rx5rFNAgOoWljQt4571R79Niugl+/lRtmpLRXJk\n2vCUcR1yiNhRclN/WUO+aKsQlUGcsRtv6N72LYu+sX1vS7ltwsQJFoclRBvSIUcHAPANFHYA\nfuTUqVN//vOfhw8f/sQTT7DO4lMcx905dA5x3Nd5n+pUep6TOTy2Qd1zpg/6TVRwXEcdRRXC\nRWbJiNqY6fpycl6Oqg4AOh0UdgB+5O9//3t+fn5+fv6cOXPUajXrOD6lVeruH/Hb8X2nlNYV\n2d32SGN0UmSqUqZinQsAoDNBYQfgR5rHI3McJwgC6yxsxIV0iwvpxjoFAEBnhcIOwI+8+uqr\nb7zxRnZ2tk6nY50FAAA6HxR2AH4kPDz8//2//8c6BQAAdFaYxw4AAACgi0BhBwAAANBFoLAD\nAAAA6CJQ2MEFVVVVDoeDdYqAdvTo0RkzZrz++uusgwAAQKeEwg5a5ObmTpo0acaMGS6Xi3WW\nwLV58+azZ8+uWrUKvwUAALgGGBULLerr64nIarXa7XalUsk6ToCaNm1aaWnp4MGD8SsAAIBr\ngMIOWowfP57juPj4+KCgoI5q02azfbHkXwrFpeusd6Ck2CjvNe57sbGxr7zyCusUAADQWaGw\ngxY8z0+YMKFj21Sr1TdOnaHR6zu22dbq8o96r3Hfq6urW7FixYABA3JyclhnAQCAzgeFnRe5\nXa4mi8Wr7Xuv8Q7B87wpMkpnMHjvEE1nT3uvcd/7/PPPP/jgA51Ot379etyNBQCA9kJh50VH\nDuQd3r/Pe+17PG6adZf32gffGzBggEqlysnJQVUHAADXgFlh53K5Pvzww/379zc2Nqampj70\n0EOxsbFEtGrVqo8++uj8bjzPf/3116xCXqf0Idk548Z7r/0dG9d7r3FgYujQoZs3b/bqU4kA\nANCFMSvs3njjjVOnTj344IPBwcErV6584YUX/vnPf2q12qqqqgEDBkyZMqV5N47jWCUEYAJV\nHQAAXDM2hV1dXd2uXbteeOGFwYMHE9Gzzz47e/bsvXv3jho1qqqqKjU1tXk7AAAAAFw9NhMU\nm83m5OTknj17Nn+rVqtVKlVdXR0RVVZWhoeHM0kFwFxeXt4NN9zw4osvsg4CAACdEpseu+7d\nuy9YsOD8tzt27DCbzWlpaZIkVVVV7du377PPPnM4HGlpaQ888EBMTMz5Pfft2/fWW281f11Y\nWKjVajsqUm5u7unTp6dMmaL35twcAFe2b9++hoaGzZs3u1wujJ8AAID2YjwqVhCE1atXf/TR\nRzfeeGNqampjY6PT6fR4PE899ZQgCJ999tmf/vSnxYsXny/gLBbL8ePHz79cJpN1SAyHwzFv\n3jybzdbY2PjYY491SJsA1+D222+3WCxZWVldrKr7Zvmy7z5f7r32Q8NMdOcM77UPANBZsCzs\niouL33jjjYqKijlz5txyyy1EpNfrP/7446CgIJ7niSglJeX+++/ftWvXuHHjml8yZsyYffta\nJhCZNWtWXl5ehyRRq9VpaWnHjh1LS0vrkAYBro1cLo+MjDR4c+Y/JqbOnNU3I9N77WOEOABA\nM2aF3eHDh+fPn5+env7SSy+FhIQ0b5TJZOe/JiKDwRAREVFTU+ODPP/+97/dbjcGJAJbK1eu\nfPfddw0Gw7p167pYpx0Ac//9738FQbj6/aurqzdu3Hj1kzM4HI7+/fsnJCRcUzqAjsGmsHO7\n3a+99tqECRMefvjh1v9nDh48uGzZsj//+c/NPRY2m626ujo+Pt43qQK8qhMEYceOHTExMcnJ\nyayzBK7myxiXy9VRjxkAwHl2u71H1oir3z92YPvar62qqqmpQWEHbLEp7A4ePNjQ0JCSknL+\nvioRJSQkpKamVlZWvv7665MnT1apVCtXroyIiMjKymISMtBs3Ljxueee0+l03333Xde7FdhZ\n9OjRg4hCQ0MlSWKdBQAAOh82hV1paSkRLVy4sPXGRx555Oabb37jjTfee++9t956i+O4gQMH\nzp07Vy7Hume+0DxCRalUBnjPJVu33Xab3W6/6aab8M8eAACuAZsPj6lTp06dOrXNH0VERDz3\n3HM+zgNENHz48JUrV4aGhqrVatZZApdMJouIiIiMjGQdBAAAOiX0CsAFiYmJrCMAAADAtWOz\n8gQAAAAAdDgUdgAAAABdBAo7AAAAgC4ChR0AwKUwNhwAOikMngAAuNT48eNZRwDwa5IklTaU\nltSXaJXa7mHdjRoj60TQAoVdi8LCwrNnz+bk5GD+MAAAgCsoqin6cOeHq/JWqRVqj+jJTsrO\nTsy+K+suhQxd3eyhiCEicrlcjzzySF1d3VNPPTV79mzWcboOnue3bVinVKmu9gUSCYIgk7dj\nNa1ILRZUBQDwndqm2oUbFh4tOzokaYiMkxFRo6PxrY1vuQTXnOFzWKcDFHZERCSTyVQqFf28\n+gJ0lOzs7Ozs7Kvf32az7d27d9SoUd6LBAAA1+OHIz/sK9rXJ7bP+S16pb53dO8F6xfc2OfG\nuJA4htmAUNg1k8lkn376aWVlZUpKCussAAAA/utExYkIY8QlGzVKjVFlPFFxAoUdcyjsWhiN\nRqOxg5/9rK+tKcw/1bFttlZbXdmBrRUUFDz++OPJycmLFi3iOK4DWwYAgC7D7rK3+SydQq6w\nu+y+zwOXQGHnRbdOvMHpdF79/jt27MjJybn6/VMn3dj+UL8oLy+vpqampqamtrbWZDJ1YMsA\nANBlhOpCT1acDNYGX7Ld5rKF6cOYRILWUNi1eO+99w4fPvz000/36NGjo9ps79KrhYWFaWlp\nHXX09po4cWJxcXFiYiKqOgAA+CVDegz5Iu+LCEOEXHahhKiyVOUk5/SP688wGDRDYUdE5HA4\n3n33XSJKTk7+7W9/yzoOGwaDYd68eaxTAACAXxvVc9TdQ+5enrs8LiTOqDK6RXdtU21ZQ9kf\nJv5Br9KzTgco7IiISK1W33333fv378espAAAAFcgl8mfvuHpPtF9cotyqyxVKrlqaOLQm/vd\nnBKJ0Yd+AYUdEZEkSTt37iwqKiooKEhNTWUdB6BLkcvlxWfybU2Wq3+Jw+FQq9Xt2N9ua38u\nALhGcl5+U7+bbup3E+sg0AYUdkRELpertraWiMrLy1lnAehqjEbjvMceatdL1q5dO2nSJC/l\nAT8hOCVHjSS6SR3OKXQYiQ/QMVDYERGpVKrFixefPn16woQJrLMAAHRxokcq2eAp2+Gp3OMh\noohB8pA0Pmm6UqFHeQdwvVDYtejdu3fv3r1ZpwAA6PrOrHKfWOo29uAiMuXEkeCgwm88jhqp\n/xNqHmsEAlwfnnUA8Bdms/n1119ftWoV6yAA0JVZzorH3nOG9OKVBo44IiKZmoJT+NL/ust3\nelinA+j0UNi12LZt27/+9S+Hw8E6CDPffPPNihUrXnnllZqaGtZZAKDLasgXVEG8THXxVo7U\nJr7hlMgmE0AXgluxRERWq/V3v/udKIonT5586623WMdh49ChQ81fmM1mzFEMAF7iaSJeLl2+\nnVdw7iYUdgDXCz12RERKpbJ5dVS9PnAnVxw6dCgRyeXymJgY1lkAAprNZnv66afnz58vSW0U\nQJ2d0kiCq43toouURgyeALhe6LEjIlIqlatWrTpx4sS4ceNYZ2HGYDAQkUqlcrvd7ZpCDAA6\n1tdff719+3YiuuWWWwYNGsQ6TgcLTZO7LE63TVJoL5RxkkD2KjG0N4ZOAFwvFHYtYmNjY2Nj\nWadgSRTF5j8VCgXrLAABLT4+vvmLqKgotkm8QRPJpf9O9dNCpz6OUxh5nid3k9RUInSfqowY\nLGOdDqDTQ2EHLcxmMxHZ7XaLxYIeOwCGRowY8c9//lOv18fFxbHO4hUJNypUoXzJJrerUarM\n9cSMUHS/RZEwXsHjEwnguuG/EbQYOnRoYmJiUlJSeHg46ywAga75mdcuLHKwLHKwzGOTRA8e\nrQPoSCjsoIVarcb9aADwJbkWJR1AB0NhBy02b968bdu2bdu23Xnnnei0A2BIEISlS5cGBQXN\nmDGDdRbGTlWe2nBsQ2FNodPjjAuJG9x98Kheo3gO8zkA/CIUdi1efvnlgwcPPv/88+np6ayz\nsDFixIh169YlJyejqgNga8eOHYsXLyaigQMHJiYmso7DzI7TOx5e9nBMUEyoLpTn+cKawv/k\n/ufhEQ8/OvpR1HYAvwSFHRGRw+FYvXq1JEmffvppwBZ20dHR77//PusUAEApKSlBQUEGgyGQ\nJ5U0O8wr9q5IjUo16VvmSzeqjSa96d9b/52ekJ6dlM02HoDfwkVPi+Y5PvCEGQAwp1AoNBqN\nVquVyQJ3+o8DZw/sKdhzvqprppQrY4JjdhfsZpUKwP+hx46ISK1WL1iw4MCBA7NmzWKdhaWC\ngoLQ0NDg4GDWQQAC2oYNGyoqKioqKo4cORKw9xBqrDVapfby7TqlrtpS7fs8AJ0FCrsW8+bN\nczgcTqfz6aefZp2FjZ07dz711FOhoaHffvst5rEDYOj89HVhYWFskzCkkqs8oufy7R7Ro1bg\nBAXwi1DYERE1NDQ4nU4iOn78OOsszFgsFiJyOBxYUgyArZycnFdeecVgMCQkJLDOwkzPyJ6N\ntka34FbILloLp9pS3TOyJ6tUAP4PhR0RUXBw8H333bdnz5758+ezzsLM+PHjNRpNTExM86Kx\nAMAKx3ETJkxgnYKxnpE9Hx718Ce7P0mNSm3uohMkobimeEjikIl9J7JOB+C/UNgREXk8nmXL\nlnk8nnfeeefll19mHYcNnudHjhzJOgUAQIs5w+eo5Kq3NrylUWpkvMzmtN2aeet9OfcFa/Ec\nMMAvQmFHRORwOARBIKLGxkbWWQAAgIhIq9Q+PPLhaenTSupLnG5npDGyu6k7ZrADuDIUdkRE\ner1+4cKF+/fvnzNnDussAABwQYQhIsIQwTpFmySlJ18l5MukWg8f4ZL3ItKxjgSAwu5nOTk5\nOTk5rFMAAEAnwEnuYOdHIY73RC5Y5NS8ZOelBpG7z0LjWEeDQIfCDi6wWq0ajSaQ50QFALga\nBtc3wY5lDlmGxCmbt3DkiPEsrZCpiQaxzQYBDg8rtHjzzTcfeOCBY8eOsQ7CzMGDB8ePH3/H\nHXd4PG3MHQUAPuNwOP7617++/fbbrINA2zhyaj27XbKk81UdEUmkdkjxBtleIoFhNgAUdkRE\nDofj008/PXTo0IYNG1hnYaa4uNjtdpeXlzdPaAcArOzcufPrr79eunTpiRMnWGeBNsjFCq1r\nm8gZL9nuoSAjv5M8NUxSATTDrVgiIpVKFRERUV1dbTKZfn3vLio7Ozs9Pb1nz54hISGsswAE\ntIyMjF69egUFBaWkpLDOAm3gJLHN7RJxRETU9k8BfAOFHRGRy+Wy2+2SJFmtVtZZmNm8efOB\nAwcOHDhw//33B3KBC8BccHDwp59+yjoF/CIPH2lTDJORRbx4GKycLBZxiFqG8yewhFuxREQq\nlSonJycuLi4rK4t1FmaCg4M5jtPr9eixAwC4ApHTOuTpCqGg9eN0nORRc+esYjpxiiu8FsDb\n0GNHRORwOH744Qci2rp168CBA1nHYcNms0mS5HK5GhsbQ0NDWccBAPBfZtVtcqnK6PzKw4VJ\nnJqX7DKppkya0iCMZh0NAh0KOyIitVr929/+dv/+/dOmTWOdhZmJEyfW1dUlJiaiqgMAuDKR\n09aq/8emyFF68mVSrcCHO/lepfXxWgmDz4AxFHZERKIofvHFF+Xl5aNGjUpISGAdhw2VSnX/\n/fezTgEA0DlInMImz7bJs1ttq2KWBuBnKOyIiDwej81mIyI/nOlDEqipQnQ1SAoDp4vieeWv\nvwQAAAACEwo7IiKlUmkymaxWa3R0NOssF2k4JZ75ylWy0c0rONEtxYxQdJsojxjsld9aXV3d\nwoULk5KS7r33Xm+0DwAAAN6GUbFERA6H49y5c4IgHD58mHWWCxoLxJOfOBtOCuEZMtMAWcQg\nubVE3PW8vWqvV1aG2LBhw9q1axctWlRdXe2N9q+eJJKjRnI1SmxjAAAAdDrosSMiEgTB7XYT\nUWFhIessFxSvdVuKJV3sheJbHcYRyYp/cJvSZbyc69jDZWdnJyYm9u3bNzw8vGNbvnoeu2Td\nGbr3v46y7W4iih+niB4ujxkupw7+uwJcRnKRaCVZMC53AaBTQ2FHRKRQKBQKhcvliomJYZ2l\nhccu2cpFlenSikYdypVvF5rKJENCBxc78fHxK1eu7Ng220X0UP5SybLBZOgtRg6Wi4JkKRLL\ntjmcj6p6TMG8UOA1jgNk/ZY85WTbSroxpOxFhjtJHsk6FgDAtUBhR0SkVCqXLl2an59/ww03\nsM7SwmOnij2eiEwZXdJbxRHJSHAwiuVNlbme0h9FZZxDYQgmIl7OqcKI1/CH/umMGCzTRaMf\nBbzAvofK55AynmQRpMkhoZ4sn5OriEx/JBmzrmsAgGuGD8sWycnJN910k0LhLz1DSiMXPUzh\ncV7aLScJJHkkZZBXDnrw4EFRZLbKYf1xURXOEX/Ro3UKLac0cnXHhF96FcC1k9xkWUmqJJIn\nEKcmjifeQMre5NxLlq9ZhwMAuBbosfNTvJyCe/KFq93GJL51cWerkLpPUmoivFKRV1RU9OvX\nzxstXw2XRZSpiC4bGcIryW3pxAMp3v5gqUzejv9odru9cOl/rn5/q7lxzl23Y3nfa+E6Q03r\nSTP00u2yWHIcJJIu7S8HAPB7KOz8V/ebFdYSsWqfoI3iZCoSXOSskYJ7yRJnKLiu+HGj0HGi\nq41OZNFFCkMn/gtrdfqcceO91/6ZE8ebZ2GEdpPMxCnbqN54NZGbRBvxurZeBgDgv3Ar1n+p\nQrh+j6uSb1NoI7mag4I6jIu/Ud73MaUhoWv+1kLSZI4aiaSLPmUFB7nMYmiqjFUq6Mp4A4mu\nNrZLTiI58RqfBwIAuF7osWvxzDPPHDx48OWXXx469LL7MuwojVzKXUoiyvg9yVSs03hZ1FB5\nzFj+9Fq1O0hSaDlJIreZGouEfo+pW8/5AtBhFEmkH0fu4kvHSbjLKWgkrnsBoDPCmYuIyOFw\n7Ny5s66uLi8vj3WWtnX5qo6IeDn1vJc3jKlRGriqPE9VnkcTSZnPqDHXCXgLpyTDdHLlk1BB\nJBIRSS5ynSb1ANJPZR0OAOBaoMeOiEitVv/tb387fvz4nXfeyTpLQJNrST+ibshwjb1G5BWc\nKpjjcOkBXqUZSVH/Jss3ZF1NnIw0w0h3AxnvIbm/TGkJANAuKOxajB49evTo0axTABERJyNt\nJAo68BXNMNIMo7BnSKgjeQzxetaBgIhIEISTJ0/27NlT3p5B5QCAj08AACKZiZQ9UdX5j9de\ne2327Nm///3vWQcB6GRQ2AEA+BePx/Phhx+uWrWKdRCWzGYzEVksFtZBADoZdHFfsHbt2kmT\nJrFOAQCBbufOnYsXLyai/v37Jycns47DxrPPPjt69OhBgwaxDgLQyaCwAwDwL2lpaTExMUaj\nMSEhgXUWZoxG44QJE1inAOh8UNgBAPiX8PDwb7/9lnUKAOiU8IwdAAD4nYKCgrfffvv48eOs\ngwB0MuixAwAAv/P2229v3759586dK1as6Kg2T58+/d36TR3V2uVcTsecWfd4r32Aq4HCDgAA\n/M6AAQNyc3PT09M7sM3k5OTxd9/fgQ1eoraqSuVo9ErTkpOaNpDrJIn1JAslRS/S3UCc0ivH\nakUUxbKysri4OG8fCDoQbsUCAIDf+fLLL10u19q1a1kH8QOihWr/TrXzqWkjOY9S0waqeZFq\n/x+JVm8f+fnnn582bdqCBQu8fSDoQOixAwDwO5s3bzYajRkZGayDMGO1WonI6XSyDuIHzP+h\npvWkGnihL0YeQ00/kDyagh/y6pEPHz5MRD/99JNXjwIdCz12AAD+Zd++ffPmzXv44YdLS0tZ\nZ2GmeVbR4cOHsw7CmmQjx0+kSLz485onZSI5DpBk9+rBFQoFEanVaq8eBToWeuzAj6jV6iFD\nhrBOAcBYRUUFEXEc5/F4WGdh5p577lEqldOmTWMdhDVPNdm2kvayApfTkX0LeapI0c17B3/6\n6afXrFkzc+ZM7x0COhwKO/AjPM/j0hCgubCTJKmkpKRbNy9+bPuz1157bfv27Xv27Fm+fDnr\nLExxzR/TEhF38Q9EkiTiFF49+MiRI0eOHOnVQ0CHQ2EHAOBfZs2adejQocjIyJycHNZZmCks\nLCSiqqoq1kFYk0eRfgK5S0gWetF2oY70E0kexSgW+C8UdgDedebE8f27d3qvfYfd1v+5Z73X\nPvieSqV6++23WadgTBAEInK73ayDMCcj3Y1U/SxxfYjXt2wTreQ6Q8GP4EF5uBwKOwDvSkpN\nyxk33nvtnzmBqfmhC0pMTKyoqIiNjWUdxA/obiSxiapfJF5PvIpEJ6kzyPRn0nnxxAKdFwo7\nAAC/U1BQYDKZjEYj6yDMND9uq1KpWAfxBxwZbiXNMHKdIrGe+BBS9iR5NOtU4KdQ2AEA+J26\nujq9Xh/IhV1kZCQRxcTEsA7iN+TRKObgauD2PAAA+J3169cT0datW1kHAehkUNj5EUzFuWTJ\nkpdeeslisbAOAgCMZWVlEVH//v1ZBwHoZHAr1o8E8m0XIqqoqFiyZAkR9e/ff/r06azjAABL\nSUlJRJScnMw6CEAng8IO/EVkZOTIkSPLysqGDh3KOgsAMJabm0tEu3fvZh0EoJNBYQf+guO4\nBQsWsE4BAH6B53kiksvxIQXQPnjGDgD8Tnp6OusIwJjBYCAinU7HOghAJ4PCDgD8TnQ0pnUI\ndM2DqBobG1kHAehkUNgBAIDfSUtL02g0/fr1Yx0EoJNBYQcAAH7nwIEDdrt9//79rIMAdDIo\n7AAAwO+43W4i8ng8rIMAdDIo7DoDTyXZtpHla7L9l9xFrNMAAHhd87RHgwcPZh0EoJPBSHK/\nZ/mamn4gx37idSQ5SDBT6FMU9ABxCtbJAAC85dFHH505c2ZQUBDrIACdDAo7/2bbRtUvkqoP\nqX+e/UGyU8MS4g1kvJtpMgAAAPA7uBXr35rWkbI7yVotNcZpSNmLbNtJamIX65pIdnKdIudR\nEjF/AQD8ihdeeGHcuHGPPfYY6yAAnQx67PzIuXPn4uLiOI5r+V5qIqGa5KGX7scbyL6DXMWk\n6u3jhNdIcpJ5Jdl3kW078TyJLgp+mILvJT6EdTJfqK2u3LR29dX6rtA2AAAgAElEQVTv7/F4\n2jXbfm1VVXqP29ufC3xq9erVde0ZBuDxCIeLy3ie+/VdiYhIFMVEU/CoUaOuJZxfal5M7NCh\nQ6yDBDSXy3X69OnU1NTmhUCgU0Bh5y8WLlz4ySefjB8//u9//3vLJtFNtm2kGUqXn9s5nujX\nPyVWr159tt5y9RkcDsfpT1Zc/f5Oh2PK2JG/vkp343vU8BEpU0mbTUQkuci6ioRqMj1PnKb1\njidOnKitrc3Ozu5KJ5E/PNm+Loe1a9dOmjTJS2GAFZlMlpU9wnvti4JQejDXe+373lNPPfXO\nO+/ceeedrIMEtD/84Q/btm2bMmXKiy++yDoLXC0Udv6ivLyciMrKyi5skhlJP5Fc5y7ttJOc\nJDpJ/utT88tksrGTJndw0FbqqqvM5l+7r+o6RfXvkmYQkbJlC6ckZSrZfiRbDukuVDD19fUP\nPvigw+F4/fXXR48e7a3QANAZDBw4cOzYsSNHjmQdJKAVFxcTUUFBAesg0A4o7PzFM888k5WV\nNWzYsFbbeNIMI9urJDe2+k1J5C6g4AdJFs4g5TVwHiVZyIWqrgVHfCQ5j7cu7DQajUajcTgc\nwcHBPs4IAP7mN7/5TVNT0xdffNF8TxaYaH4sRKPR/Oqe4D+6bGF3+vTp7zZvDwkLu/qXuN2e\nj79qx7NQGrf9jjvuaH+0tplMpltvvfXSrfrJ5C6kxk9JEUOcmiQXeSpJN56CZnfUcb1OshG1\nNTMLpyDxovEfarX6yy+/NJvNcXFxPsoGAP7K6XQSkSAIrIMEtN///vfr1q2bMWMG6yDQDl22\nsBNFsc/AjLjuPbx3iMLcbd5rvAWnotC5pM4i1zFyl5M8jJQ9STuKOK3XD91R+BCSHG1sl+wk\nv7TsNhqNRqOxjZ0BIMA88sgjS5YsaeNyN7AJTqn+hOiokXgFaSK5kJ6yNh7C7jiDBg0aNGiQ\nFw8AXtBlC7suREbaUaTttIPd1OkkNpBoId5wYaPkJk8FqTLYxQLwqYqKis1vL/Be+5IkDR3Y\nz3vt+979998fGRmJgUStNRaI+ctd5Ts9Ci0nipKnSeoxVdnrHqXS6M3iDjobFHbgZfJoCn+V\nqv9I8m4kCyFORmIjuc9S8P2kGco6HICPhIWF3Zw13HvtS4JgLyvyXvvAnKNWPLHUaT0rhg2Q\ntdRxIpX+1yN5pP5PqDEpLZyHwg68T38LySPJupY8FUQiKdMo+CHSTSCv3kIA8CcKhSIi6tdH\nsl8zURBKK895r31grmy7UH9MDEpuVcHxZEzki9Z64sYKoX1k7KKBf0FhBz6hHkzqwUQCkdj2\nWAoAAPhl1nOiKujSi2GOJ4WRs5wVUdjBeSjswJdkRDj7AAC0m+iR2jx98rwkYegwtILCDgDA\n6+Li4kp/2nP1+1dWVhqNxqufP0ySpMTExGuKBp2DNpyvyhPUoRwR2Vy2elu93WWXyeTy+jCV\nScU6HfgRFHYAAF7Xv3//du2/b9++hISEiIgIL+Vhwmq1iqJ49fvb7Xaz2Xz1+/M8r9fr25+r\ncwjPlB3/2Kkx8eX2c8W1xTVNNQpeobYGV2t35JfVPml/NEgTxDoj+AUUdlcieSRnGeeqlHg1\np44jeUCsWd9CkshVQa4KiYiUUZwyijgMdQDwvoPnDuadzTtw/EB8efzApIE5STk6lY51qI7x\n9nsf9ezT9+r39+iCf9yz/+r3P3H40PNzn2p/rs4hJFU28H/Vu19rPGk/pzEqo/juMqvaHlkh\nZRatObVJrZf/bsLvWGcEv4DC7hc5zlL1N9S4U+LVkiRIkpOLvJPCbiFe2fULHHcDVa+S6n4k\nXk1EJDqksBspfAYnxwUhgNdIkvTpnk9f/eHVSEOky+46Vn/sy8Nfju41+n/G/U9McAzrdB0g\nIjqmT3qm99pvqKvzXuP+oNtNiq8qv63YWRUjJtl5jzOkvrF7gVtvSRVSP9j+wW2Zt3UL68Y6\nI7CHwq5tnjqq/Ex0nuV0fSSO54hIdErV33CSQBG3sw7nZZKbKj6VrAdJ14fjZBIRSQLXsEvy\n2KS4RzgO/2QAvGP76e2vrXstPT5do9RUV1fr9XqNWpNXnPfvrf9+8ZYXZTwGHl0vq9V69EDe\n1e8vShLfnlsVjQ0N6UkJ7c/VDgWqn+oHVwvqwtYbFTKFXq0vri1GYQeEwu6X1O+UbPmcpjud\nn2uNV3GaHlT5uRQ8nFN6cToq9qyHJfMuSZfGESc1b+FkkjpBMu8iazYZMrp+hyUAE9vyt8WH\nxGuUrQZMcNQjrMeq/atmZs1MjUplF61j5G7bsm3DOu+1L5PJHrzzSquaTp06tXkJ2qu0cePG\nMWPG8Hw7Jv/19kN+oiRybdWaPMcLIgbHAhEKu1/iKCZF8KUbOSXJdWQvpq5d2NnPkjz40smD\nOZ6TB5PjLBmwDBiAd5Q3lhvUBk7kVeagoGqNjHeRWuR5Xq/Wl9aXdoHCLmvEqJxx473X/o6N\n66+8g0qlUqnaMYBUo9EYjcZ2FXbeFhMcU1RbZFAZWm8URMHqsMaFxrFKBX4FhV3bJA9xbd73\nkHOiW+raSyZIrrb/7pyMxHZc6wJA+8hIHlaYFlucGXImxSU6ZaRs6HekakCeJEpyWVc4V1eU\nluRu23yFHThO0MgaRZI5BaMkcW63R6Fox1+89GzRdSb0f6N6jvos97NQbej5nl2JpDPVZ24f\ndHtyeDLbbOAnusLJwhsUweQoJNlFF0UkEYl2SRHsR1dv3iAPlgR7G9sFuyQPoa5d1AIwlFwx\nzLMhVYxqNCcU2e12pVylOxcfaVVx3XZ3D+vOOl0HuNKQVclO5pVk3032bSRJpLuBdBO+3yq/\nadJkHwbsBIYlDZt347xX170aaYjUKrUe0VNjrRmbOvahkQ+1eYsWAlCXLewkSTqct7fsXPG1\nNhBMRYnktpJMurDNrKQox7GqQtomEhFnqe+AoP7H0I8ve08U7cS3etRHaCJPA+n7dfGiFoAV\nwSEllQ7brt+rlslUpCIiiROdQfVcqXFW8u+7+EPxkovqXifrGlIkkmYYkUSeMqp5PkY/lgiF\n3UU4jpuVPSs9If1gycHS+lKjxphoSsxJ7jpz4sD167KFHRHZbU1Wi+UaX2yyKvrz8gMJktFB\nSg+JPGdTiuEWT9pZwWkmJxGRxuPpwLT+QxlNcU9wpe+QIpzkeomIPFbOVSXFP8mpuvTDhQAM\nmYtEz7GgzN79z1SfqbHUkEicwLlFd8+Y3r1pIHXtB0DsW8n8NWkyiJovHTmShRLXLzHoC3I9\nQUqsqHGpvrF9+8a2Y0ZACChdtrDjOC5rxOi47j2uuQXpJrIeJMtByV1FvIZUcRSSY1JEXGiw\nMHdbRyT1oqKios27917ba1W9w/Tl3RSlBiJy6yxNfc7u+6mGfrpoH4fDfveMadefEwA8duLl\nFB0UHaINqW2qraqt0ml1YcawYHmYjEhwkawLrxrlPEyKqJ+rup/xGreoJ9dhFHYA7dJlC7vr\nx/FkSCdDeie+TI6Oju4z8oYOaqzX5Zsaa2vlcvwTAugASiMnuIgkUivUscGxSrdSr9drNBqX\nWVIG8V25qiMiwUKkvHyzSAoS2rGkGAAQCjv/ceDAgeXLl0+dOjUnJ6ej2lSpVBFRXrx7KpfJ\neHuj99oHCBzG7nzsaHlDvqCNbNVxJVFTmRQ/rqtPTcwbSWpjyL1BryQZlrsBaB8Udi2qzOVV\n5gpB9JiMkTFB8b4fXvTBBx/s2rWrqKioAwu7a2NzNZXVnzU7GoLUITEhCRqllm0egEDAyShx\nmmLrU27BRZpQThI4t5lzFAtRQ+XxE7r6iVo9kMxLSRFD1KqEFW0qeROpBrKLxZ7NZrPZbCaT\niXUQ6Ey6+vniKlgc5i/2frxq3zKNUsMRb3fbxveZckfWfZFBsb6MMW3atNLS0ttuu82XB73c\n1pM/bjmx7kDxbpVc7XDbM3sMG5t287CUMWxTgS+dPn1aEIRevdq4+Q5eFdLTPuHfOxqO5bvN\nZuWBcEVkelCPAQkTFApdJ34g5Kpoh5PhNrJ8Q4oexBuIE0loIHdBrfRwiCwuYIfiO53OqVOn\nNjY2vv/++/369WMdBzqNQC/sRElcvnvJlhPr+sVnKmVKIhJEIa9op8Nte3zcH3Uq7y4O01pq\naurkyZOzsrJ8dsTL7TmzdcEPf+5hSukfP4iISKJqS8Wra5798/S3MroN9UGAN954o6ys7IUX\nXggOvmzdD/CJkpKSu+66i4j+85//9OzZk3WcQOIuofq3NPw2zYBQIll8to2TFlHYc6S/i3Uy\nH5BT6P+SPIEcudS0gYhIf/M/ltR/8f2y8PAf16xZwzoeG1artb6+nogqKytR2MHVC/TC7mT5\n4e8Pruobl3F+gW0ZL+selpxXtCu3YOuYtEk+S7JgwYItW7Zs3Ljxk08+8dlBWxMlcevJdfFh\nPYK0IS2bOArRhglhns3Hv09PGOLt29Pl5eXLly8nopycnBkzrrTgI3iPTNbyHwGTnfqWSA1L\nyLGfVAOav+fkRFIs1fydlImkHsI2nC9wGgr6DQXdQ54q4hQkC/3yhyGiSFVVVayTMRMWFvbu\nu+9WV1ePGYN7JtAOXbmwO7x/369OUJxbvp1zU21VJRGJosTzLR9mnJO25q7T1Vzp8TJZU0eO\nG0hNTd2yZQvDPpIaS+XWkz+mX9YzF6I1bTq+9r4RT4bqwr0aIDo6+uabby4vLx85cqRXDwRX\nEB0dvWLFCkEQUlJSWGcJJK6TZPmatBcXcJyGlHFk2xIQhV0LjuSRzV/Fx8cXFRUZjUa2gdga\nNGgQ6wjQ+XTZwi45OfnBqKhf3c2xr+yoNax7WDQRlZaWxsbENG/XWeQp8bHTx426wms1Gs0V\nftpeDz300B133BEU1JFDwOLi4o5v/ZVVsc+rtFU2NTY2VpRd+gOJXLamozs2mNRtPMA7duzY\n6wzZ2syZM2tra0NDQzuwTWiv5GSsOOlz7lKS6y+dyI2I+CByl7IIxN6rr776+eefT56MlScA\n2qfLFnY8z1/pUs9TRp5SIrFXhNolOpVKJRHJZLLmL4jIKTrjTHE+vljs2KqOiPr379+/f/8r\n7eGpIE8JkYdkkWYh9Ltz34VFhJ1fW7qZzWlLSkq659Z7vL1kTX19/YMPPuhwOF577TXceoDA\nwslIktr6gdhGtRcYkpKSnn32WdYpADqfLlvY/SLRSo0fUf07xGuJ+BGaxmezQtYWV0jyC917\nDrejwlwxNMkXwwWYEW1kXkZ1bxOvIZKR2GQMunNAROj22qK0qLQLixdJVFRbNHPITB8sRHi+\nB1ShUHj7WPBLGhsbJ06cKEnS119/HXUVfd7QMRQJJFpI8hB38TlZqCdNNqNMwN7HH39cUlLy\nv//7v1otpp2CqxVohZ1I9QvJupY0WcQpiUgmeUZE76tq2Plt6cAmt9PsMFsclnO15+aOnzsg\nboAvk3ljguJfJlH9P8m6ijSDiFMTEZFAts339tZyZQP/e2pfdHC0Sq5yepxlDWUje468Y9Ad\n3o9E9fX1DoeDKKAfl2buk08+cbvdRLRo0aK//e1vrOMEDEUiBd1Hlq9IlXZhLjdPHbnLSDeR\naTJg5uzZs4sWLSKikJCQxx9/nHUc6DQCrLBz7Cfz56QZfOHUyclDQobe029nRFTyx9urIgwR\nmd0ysxOzhyUN83E0n05Q7DxK5k9IM7jVPwAZKZKjtNv+OPL2vvHZJypOVFuqIwwRt2fePqnf\nJKPGF7ekMXiiWVJSEsOjz5kzZ8WKFZIkzZ07l2GMwMNR0IMkOcm8kvhg4hQkWUmVSaGPkzKV\ndTZgQ6PRcBwnSVJYWBjrLNCZBFhh5zpFstCLJjcnIiKNJmFKai95w+RJk3w3v8klevfufeDA\ngb59+/riYK5TJAu+/LfvEoNMXOnMrD/6IkNbXnrpJVaH9h9sZwZWq9VbtmyRJOn8vCfgI7Jg\nCnuOdDeTp5hEG8nCSZ1OMnyiB67w8PB33nmnurp6woQJrLNAZxJghZ3kJK7N57fkRG2sVOhL\nx48fdzgcR44c8cXBJGebv3qJZCQ5fBHgFyxZsqS8vHzu3LkGg4FhjEBmtVpvu+02QRBWrFiB\nfgKf40g9kCigF9HyK/HxDJaXPE8UxUcffZSIXC7X1KlTWcWATifAxlvJTCQ2tbFdbCKZdydp\n+1XTp0/v1atX86T/XicPJ9F2+WYZ52D4PlRUVCxZsmT16tXr11/tFC3Q4Y4fP15TU1NfX797\n927WWQAY69OnD8PCbuHChc1fvPrqq6wyQGcUYD12mixSDSCxjvhWE4uIdhIqST2E6FdmM/aq\nMWPG+G6OD3UmqTNJqCJZqxnjJKeKr5PUWaxOY1FRUXjGjrnGxpZpt2tra9kmAQhw50fCSm1P\nhQPQtkDrsYskwxRyHCFXEYkNJJjJU0LOnyj096TqwzqcD/EhZLydnMfIXfjz+1BKjgPFlimk\nzmSY66WXXlqyZInJ1MZMyOAb5zvqdu3axTYJQIA7P44qPJzxDSXoXAKsx46IdBMpJobsm8lV\nRCSQYhBpsgNxpijtaIr7gmybyFVIkovk6RQ6tOSUrdUUdhCIzo+ZQCcBAFs33HDDokWLamtr\nz9+TBbgagVfYEZG6P6nbWI+B7TQTDCjTSJnWekPPnoVY+j3A7d27t/mLU6dOsU0CEOA4jlu+\nfHleXl5iYiLrLNCZBNit2CtiO82EP+jRowcKuwA3e/bs5i9mzJjBNgkAAFwDv+uxE0Vx+fLl\nmzZtEgRh+PDh999/P+bTAi9xuB27CnblV+bX2+qjg6L7xfYbGBPLuY6Rp5xkoaRIIeWlF8oe\nO9WfEBy1okzF6aI5Y5Kss1TCkkANp0RbpSgKkjacD+7Fy1SXRj9Xd07Zwz5hljpIYR89ttZp\n2azSD/15bRLoms7WnS1vKOd5Ps7AR+tEIpFkUaRIoIufynDUSrZyUXCR2sTp43gOfQIA/srv\nCruVK1euWbPmySeflMvlixcv5nn+gQceYB0KuqC6proF6xd8f/j7CGOEQqawu+y99PmPD+qW\nFurkeA1JbhIbKPgxCn6oefU5Iqo/IZz+3FWxR1DqOFEgl0VMmqHsOVOpNPp7cWevlE4sc57b\n6FEaOOLIZZGic+Q9ZyqDEi98Pn9/+Pv3Nzx5Vy/7vUOsFhfnrH6j1vGO0XSrPvL3pOjGMDx4\nSV1T3Xvb3vt418dhGuXk7jUzelRXacPDdCYZ5yLjbAqeQzITEYkuOvOV6+h7TpmSI5kkOqj7\nzcrkOxS6aBR37Xbw4MG9J84olFe7HLYoSm63u/Cr1Vd/iGidChMaBzj/Kuw8Hs/atWtnzZo1\nbNgwInI6nYsXL7777rvVavQZQAf7cMeHW05tyUjI4HmeiBJ1Vfd0b8wr3XvkJ1VGz4zk5D5E\nTmr8iHg9Bd1LRE2l4slPXbYKydS/pQtZ9HDn1nuIoz4Pqfy5367ozNldrzQIxQbe5CBOIiIy\n0slt6vwTNtWNBZzWQ0Ql1pIPT7z64lApWuM5Va2XJEkmlxU12HtZvlAUlRyvvluUrvRRlJGR\n0b17d5/8baBjOD3ON9e/ueXUluzEIbfEHukf7C6zxx+pt8SHygfGD+Sta0k0k+kF4tSnVrhO\nf+ky9ZfLVEREkihV7vW4rNKA36r8/5LG31it1oFZQ3XenIC9MHeb9xqHTsG/CruSkpKGhoaM\njIzmbzMyMmw225kzZ/r0CaS5SMD7yhvLP9j+QVaPrOaqjoh6B52rdxvrGyvkMtfhw4eTk5OJ\nVKTqRfZdZJhGfFDZNk/jadHY40IvBS/njIl0ZpU7fqw8KNl/Hxg4uPasvKRfSD9Oav0pHEO2\n0xStSggZyRHRgV1LJvcO6R1Vc84Wqvn5Q0cl6Y6ba8YnnlQnB9vkV1o9+VjuNhR2ncv2/O3f\nH/k+o1tGD13NoNCCIlu4RFyINqSwpjA2ODbSmHps38J3Pv/c0zRlYsOfTQNk/M+fFRzPGRK4\nmgNC2VZP91uutucJmtnt9rWffqxQKr13iPiwYO81Dp2CfxV2dXV1RHR+ISOdTqdWqxsaGs7v\ncOrUqS+++KL564qKCvTkwa+69957bbZLl9mol9WXa8u3Fm9t/lYpk8ZMqD1dL/EiL3KiRNKG\nDRuaf5QcV7f4lYKi8qBeJdMNzmjXGeslTYUJPaznUoOSvf33uHbmYk9tzRlH3qUTDqscxkMf\nFJRs2k5E+8Tddw0qL69x1zW5Wu/jJGdRSc2+da/vOtnzCocYNiijw2NDx3I6na1XE1mbv5Z3\n8xXlFalxpXU2rqnJ3rydE7ni8uJjeceUvDwhyrHnO7s9ylJZfelKg4KkOLipqdhQ2Xpjt27d\nUN9fWVRUlE5v4OVXex0oSZLb5Vaq2lEIxsbGXlM06Dr8q7Azm80KhUIuv5BKq9WenwqfiEpL\nS1etWnX+W4UC14vwK+bPn+9yuS7ZeKzy2J/W/alfVL/mb1UywWQ6oJCJco3TarZmDcnSa/XN\nP1LTkSefeMxBvc++b3CU83LDpbO7ucvVgsevp3zr329Afr5DGXHptM+eBnlaqsk0NYWI/rmn\nyRTaqNM1BcvVRCRJUvP4aIvTEhoaNm708IHDb77CIbp1w0N4/s7pdJ6/XCGiw02HG9wNQoPQ\nZKyx2V1NPy+16JJctZ5ag8IguDmVQpSTqqGxzk31l7Qmc6ldrvraDRfd9cvKykJhd2V9+/bt\n27cv6xTQxflXYafX691utyAI50fC2mw2vV5/foehQ4d+8803zV//7ne/O3DgAIOU0Kn06NHj\n8o2mWBNtIJVOpVFoiIgjEuV6o8piFxQ9E3vGxcS17Ce5yU7dEoaQPEZIdZVVenTBFz1UJElU\nd1bQRfr1U+ThPQwlvCI0/NJOArNFSE43JQ1KIqLBlsG1VfuCNDYHrzm/gyRJZrc53KDQxA3v\noR3k09DQ0YxG41/+8pfz3y7btezDnR+mRKSYgotMIYfc6pDm7Q22htTY1B6mHvnHfti6X/7W\n//1j7yuOsD6XzvFpK5eih8t7z8FD+gB+x78+kEJCQujnG7JEZLfbHQ5H88ZmGo0m9mcKhUIU\nRTZBoZML04fNHT/3ZPlJt+AmIoko3xIdrGiwuy3xIfE/7yWR+zQF30/yGCKKHCSzVYrCxbek\nbOVS1DB5SKr/PmBHRKYBMlM/mb32ov8sbovkqJMiBrVc2o3pNWZTsVsra9LKnM1bJJLMDnO/\nyDC1cSyps3wdGrxscI/B1eZqm8tWbAvXyZ1qmZuI3ILb4XaE68NJtKf0MD7/129C0mThA2TO\nhov6pEWPZKsUw/r718cHADSTzZ8/n3WGC4KCgn744YfIyMjmRSD27duXl5c3Z86c1jdnz/vq\nq68OHTr01FNPhYaGXv5TgCvrHdPb4XasObzG6rRaHJYDlWaPx3p7qj5MZyDJTWIduc6QdjSF\nPEG8log0EbxST0XfuyWBBI/kaZKsJWJwL77X3SpNuF+PDZSpOXUIbz0nNZwSRY/ksUn2StFa\nIg16TmMa0FKShhvCw4NSFu/YGKl1xGkbRMFBYkPPUGVk+Dil6QlSxF/5ENDpmPSmUF3oF/u/\ncIqqBndwtqnA6XFaHZbMhH4mLZHrFIX9gXQ3yFSc0sCd/twliRyRJLrIZSZzgZgyU9HtRqU/\nDwYHCFicv60IuXz58h9++GHevHk8z7/55pvZ2dlz5sxpc89Zs2Z98skn+fn5ycl+/OA6+Ldj\nZcfOVJ+xOq1h+rB+sX2jlWfJ+RO5y0geSspU0o6+ZHre+uNizSFPU4UoV3OGeD4qW64K6Rwf\nbrYKqTLXYy0RJYG0MVxEhrz1CN9m5+rOHT+3TrLvDlFYQ/URCZEjlIYxzZOZQZe0u2D3phOb\nKhor4rU1gyMaeocHRRhMJI8i7WjS5JzfrbFALN3kbqqQBKekCecjMmRRw+SYoxjAP/ldYSdJ\n0rJly7Zs2SKK4vDhw++7775fWnkChR0AAABAa/41eIKIOI6bPXv2+QUrAQAAAOAqoTMdAAAA\noItAYQcAAADQRaCwAwAAAOgiUNgBAAAAdBEo7AAAAAC6CBR2AAAAAF0ECjsAAACALgKFHQAA\nAEAXgcIOAAAAoItAYQcAAADQRaCwAwAAAOgiUNgBAAAAdBEo7AAAAAC6CBR2AAAAAF0ECjsA\nAACALgKFHQAAAEAXgcIOAAAAoItAYQcAAADQRaCwAwAAAOgiUNgBAAAAdBEo7AAAAAC6CDnr\nANfryJEjjY2NrFMAAADAr0hISAgPD2edoovrxIVdnz59QkJCpk+f3iGtaTQahUJhtVpFUeyQ\nBttLoVBoNBqHw+FyuZgE4DjOYDB4PB6bzcYkABFptVq5XG42m1kFUCqVarXaZrN5PB4mAXie\n1+v1LpfL4XAwCUBEer2eiKxWK6sAarVaqVQ2NTUJgsAkgFwu12q1TqfT6XQyCUBEBoNBFMWm\npiZWAXBKxCmRvHBKfOeddx599NEOaQp+SScu7J599lmFQnHq1KkOaa2wsLCurq5fv35KpbJD\nGmyv2traoqIihlczgiD89NNPRqMxJSWFSQAiOnnypNVqzcjI4DiOSYDKysqSkpLExMSQkBAm\nARwOx9GjR00mU7du3ZgEIKIjR44IgjBgwABWAUpKSiorK1NTU3U6HZMAZrM5Pz8/KioqNjaW\nSQAiOnDggFqtTktLYxWgoKCgvr6e4Smxrq6usLAwPj4+IiKCSYDmU2JQUFBycjKTAOQHp8SK\niorS0tKkpKTg4OAOabB3794d0g5ciQSSJEnSH//4x8zMzLKyMlYB1qxZk5mZ+dlnn7EKYLFY\nMjMzn3zySVYBJEl68MEHMzMz3W43qwAff/xxZmbmxo0bWQUoKCjIzMx86aWXWAWQJGnq1Kk3\n3HADwwBvvPFGZmbmoUOHWAXYs2dPZmbm4sWLWQWQJGnYsGEzZ85kGODZZ5/NzMwsLy9nFYD5\nKdFsNmdmZv72t79lFUCSpDlz5mRmZno8HlYBPvroo8zMzEC2jjMAACAASURBVE2bNrEKANcA\ngycAAAAAuggUdgAAAABdBCdJEusMfuH06dN1dXUDBw5k+IzdmTNnunXrFhkZySSAIAh5eXlB\nQUG9evViEoCIjh07ZrVaBw8ezOqBkvLy8nPnziUnJ4eGhjIJYLfbDx8+bDKZEhMTmQQgooMH\nDwqCkJGRwSrA2bNnKyoqevfu3TyMw/fMZvOJEydiYmLi4uKYBCCivXv3qtXqfv36sQrA/JRY\nU1NTUFDQvXt3hs/Y4ZRYVlZWUlLC8JQI1wCFHQAAAEAXgVuxAAAAAF0ECjsAAACALgKFHQDA\ndXFZret///v9779vq60t27dPYjSjLwAAobADALhOuxcuLNm9O7hbt28ffPC7Rx/dPH8+60QA\nELhQ2AEAXDvR7T79/ffZc+c2FBa6zOa0GTMKN24kDEoDAEY68ZJiwJDTbK45cSIkMVFrMrHO\nAsCSrabG43QqtNqDy5YNeeopURBUQUHEaHKKQCN6PMRxvEzGOgiAH0GPXfvkr11beeiQyGht\ncj9RdfTof265Zc3jj6+YPr1w06bmjS6LhW0qACa0JpM6JGTjc88Z4+Jis7KOffFFtxEjWIcK\nFGd37Fg6duyP8+YdXbmy8exZdJQCEOaxa6+vZs+uPnZMqdPFDBoUO2RIbFZWcLdugXZ1vvrh\nh+11dQNmzy7avLnmxIm7vv5acLmWjhs3/I9/TJ02jXU6AF+ry88v2LCh5+TJ+/7974oDB6Z+\n8IGO0Zy6gaYsL+/Q0qVl+/d77HYi0kdGxg4ZEjd0aOzgweqQENbpANhAYdc+NSdPrrrnntCU\nFHNJSfOpRBcRETtkyLC5c5UGA+t0viCJ4vvDhg156ql+d9/tcTg+nTQpe+5cldG4bu7cGcuW\nmdLSWAcE8CnR4zm4dGnl4cPdRoyIGTRIHxUlY7RUQ8AS3e7jX3216403RFE832k38k9/Sp0+\nnW0wACZwK7Z9dOHhRHTTW2/dt2nTlP/7v4yHHtJHRZ3dvl2h07GO5iMcx6lDQo5+/nnloUNy\ntTp12rQTX31Vsnu3ymgMY7fwDgArh5Yt2//ee2W5uR6HY9+7757dvp11okC0/733Btx778N7\n997x5ZehKSmhSUkRffuyDgXABgq79pGr1aPnz9dFRPAKRVR6+qBHHpn6wQf3rFnD8QHzTnLc\nyOeeczY21p46RUR977qr+vjx4199FTNoUAC9CQBERCSJ4k8ffzxs3jyZShU7ZAgR/fTRR6xD\nBRyPy2Wvq/M4HJIoBnfrFjd0qCooKDQlhXUun8JkinAeRsW2j0Kr7XnLLZdsDLQ7LwkjRtyz\nZo0oikSki4jod/fdP330UfOnGkBAcZrNLqvV43DwcnloYmJkv37Vx4+zDhVwlDpd4vjxh//z\nn9LcXE1oaHleXtqtt7IO5WvNkykmTZjw7YMPms+dS5k0aczLL7MOBWygi6Xd9rz99u4332Sd\ngqU9b7+97913lT/ffc564onx//gHRgKCj9lqa83nzrEdCKkOCgru3n33woVytXrfkiU/ffxx\nwrBhDPMErDEvvZT+wAMeu702Pz/pxhsHP/YY60Q+hckUoTUUdu0jiaLWZDq3axfrIMy08Q5w\nXI+xYzEM0Mdw5+XkN9+smD79P5Mnb/3LX878+KOjvp5BCI6b8Prr0ZmZ1oqKA++/H967d+Yj\njzCIEcBEj+fABx+s/8Mf9FFRty5fPnv9+jEvvxwgQ9nOaz2Z4qDHHgtJSsJkioEMt2LbpzQ3\n9/Cnn9rr63/6+ONekydrQkNZJ/I1vAN+Andewnr1Uup0Lqs1//vvT3zzDRGF9+7d/ze/SZow\nwWcZ9rz9tiQIt7zzjigIotstV6t9dmho1jx4heO42MGDt/71r4k33NBj7FjWoXzt/GSKYT17\nxmZlrfvd73ALJZChx659OJ7XR0dLgpC7aNGyCRO+mDkzf+1a1qF8Cu+AP8CdF9Ht3vqXv/Sd\nOfO+zZvv+vrriD59ojMynGbzxueeO/HVV77J0Lr3mpfJUNX5HgavEJHH6eQVilv+9a+MOXPG\n/+Mfef/3fx67PX3OHNa5gBn02LVPbFZWbFaW22Yrz8sryc0t3bPHbbezDuVTeAf8AZax8jid\ntpoat80miaIuIiI6M7P25Mk7Pv9854IFp777zjcTmKH3mjkMXiGiLS+9NOK550JTUpoHAg99\n+mmVwRBoQ/qgNRR210Kh1SaMGJHQ3NcdSN0k5+Ed8L3Wy2LizotSrz8/EFIbFlaWl5c6dSqv\nUAR362Y+d843GZp7r201NbmLFuUuWhSakjJg1qyUSZN8c3SgVoNX9FFR+5YsOfHVV4njxrEO\n5VuSVLZv3+6FC0c+/7zodlcdOVKSmxueltZt5EjWyYAZrDwB0DkUbdmy+cUXYwYPjs3Kihs6\n1ONwFG7cGMjLWAku14EPPjj9/ff2+vrYwYNHPv/8+j/8oXz//tHz518+J5H3tO697nPnnb0D\nb6INthqKira/+mp5Xh4RJQwfPnr+fJXRyDqU17mamgSns7mTuOi///3xmWdiBg2qOnq0eT2k\n6UuXhvfuzTojMIPC7mq5rNYtL78c1qtX6rRpDYWF0RkZATgfL94Ehn5pWcyQxMSghATcefE4\nHF/Nnp10440ZDzzA7K60JAXUDXGGTn33nTE2dtebb0798ENeJgu0wSunVq/e/NJLoSkpcVlZ\nsUOG5K9Zc3rduuSJE2MyMyMHDAjp0QP/DgMZCrurtfWvfz3z44+jXnwxd/HiwByESHgT/ACW\nxSSiutOnN8+fL7hckf37x2ZlYcX3QGOrrv7kppuIiFcokidOjBsyJDYrqzY/v6mysteUKazT\n+ULN8eMnV68u3bOnobj4/Maekyf3uf12U2oqrrcDHAq7qyK63R+NHj3smWds1dVHVqzoMW5c\n/po1D2zfHlBXRXgT/IHodn96882p06YNfvzxhuLiDc8+S6I49q9/DagFlH6cN+/czp2GmJjG\ns2ebJ/DzcWkrCkLR5s28TBaVnq4OCvLZcaGZJIrVR49u+/vfHQ0N9ro60eMhIo7ne02ZMvL5\n57199NZPuzLXVFVVsnt3aW5uaW6uva6OAu8yDy6HwRNXBYMQCW+Cf7h8Wczqo0cDqqqTRFFl\nNPa+7bbsuXNdFktZXl7Jnj0+XvE995//PLRsGRERx5lSU+OGDOk5eXJwt26+zBDIOJ6P6Ndv\nyP/8jzEmRhMWVr5/f2lurrWiIvPhh31w9LM7drR+2jUoPp7haVBlNLqbmiL7909/4AHR4ynN\nzY3D2icBD4XdVcEgRMKb4B+wLGZpbm7hpk1yjSasV6/47Ozuo0d3Hz3axxmKt2yJy85OnTKl\n/MD/b+/O46Ku9v+BnxmYgWHfQVaVTUxRURZZblhuaRKX1CiXNi2XzKtpmXW5anb79dDSrLya\nloR142LmvmbuC4yAC6DsCALKIKssAzPzmd8fn5wmRCz7es7HmdfzL4Y+PT6n0wzz/pzzPu/3\nhcr09IvJyb1CQhDYUeZ5pz+1d3S0d3Q0tftKrazcBg+uzMi4dvw40ct2ZZISkL5mzZXt2/mf\nLV1cnIOCrFxdrUaNojYAJF4LELZie6K/5F5fVFR65IixHULssunAcBI6mptv5efb9+1r4eRE\n546CpenszN68ueTQoc62Nq/hw6MWLzaqBkpVcnnWV18pcnI4jYYQ4hgQEDx1KuUiI4cWLnQN\nDh780kv8y5aaGpm9Pc6vUMY21VII2a6cWp0cGxsyc2bl2bOWrq5lv/wiMjGZsHGjU1AQtTEg\n8VqAsGLXky5L7sNmzSIiEeXyj0X799t4ejo/9hiTlI4uM+Dg58ekBqYiL2/f7NmqtjZTmWzE\n8uV8y6DO27cpBzScWn0pJaUmJ8cnJsZv7FiJhQXNu+uYSKWhc+aEzpnD5O6stCoUZz7+OHTu\nXLfBg+M2bVK1t7Mqkd1QWtpUUXH9zJmmigp+ncbK1ZXa3bucBqV2XwHK3LChobTU2t29YPfu\n/J07CfW4Knvz5kEvvqif7Uo5JaClpkatVFo6O4ul0uCpU6VWVjaenjSjOr4FTuTixboWOEX7\n9o1YvhwpOmwhsOtJD0vu1GKa3NTU2itXpJaW7sOGeYSHe4SF2fn4UPvYCGEGCCEZn31m6eIy\naPr0a8ePn1292js6WtPZmTJyZPS77/aLj6c2DOZdKY35NGhTeXl1ZqapuXn62rXlJ054hId7\nhIcPefnlyLfeolwi+/z69c1VVXZ9+hTu3VuwezehGE+01dYeX7aMECKWSE6uXMnwNCjzhxzm\nqZZCyHa1cHIys7GpvXpVLBZf2rq15tIlr6gomgNA4rUwYSv2/tguud8qKPhpyhQHf//mykq+\ngJmli4tHeHjkwoXU1qvYzoCW476OjAx/882BL7ygViq/Hzdu+MKFZjY2hxYuTNi6ldrjqZbj\nvn3iifA335R/8cWETZsufP11c2Xl31NS6Nydx/w0KFucRiM2Mbn6009F+/cz3Ir9ado076io\nYbNm6eKJ/gkJdL7R2Z4G1Xdxy5asTZtEIlHo3LmK3Fz6DzmV6elHliwxlcnC3njDa/hwJv3c\njrz7bunPPzv4+emyXaPefpvyGOoKCwkhlenpGevWiSWSCRs3ugYHU7s7p1J9N26csqHBMSBg\n5P/7f4feess9JCT63XepDQC6hRW7P4ThkrulszMh5KnPPpPZ2/PtYqoyMipOn378n/+kMwAe\nwxkQiUTm9vZ527a5DBjgGhzcLz4+f8cOx8BAMxsbx8BAOmMgAuhKyXyJgiF+H5ZfpQv6+9+D\nEhJYdStuLCszt7XVdHYSQqTW1pSPbrA9Daqj5biL334buWiR/IsvPMLDFbm5F5OTKQd2IrHY\nwc9PkZNz/F//IoxSLUcsX27r5VVy6FBdUZHvmDGhs2fTvLv+omn/iRP7PPGEiVRKOfNbLJE8\nvX69LvFa3d4+5NVXaQ4AuoXA7v7YLrmbmpvHLlvGf1zdhgxxGzJk2Ouvazo7aZ48YrzpIBL9\nbenSY//6V11hoWtw8IDExNzUVMWVKz4xMTQngXlXSiGcBmWF34e9duIEIcTSxcUzIsIjPNwj\nNJR+t+L83bsr09Mr09MVublew4d7hIfTrwdblZFRqdFELFhA+TSoDvOHHE6tVuTmSq2tw954\nw8LFRZGbSzm+1w+qnv3hBybptswzQzLWrdNqNBELFjBJvIYemCxbtoz1GITORCptKCsr2LXr\n2rFjZUePlh4+7B0d7U0rlcFEInEMCOjyS8pJ02xngBBi6+MzYPJkez8/E6lUammpbm+/kZ09\nIDHxYfdDbFUojiclKZuapFZW5ra2HuHhDaWl9UVFNy9ccB86NHLRIlMzs4c6AH23q6sby8qa\nysvLjh69vHXrtePHJTKZo3FUsLP28AiePt0rMtK6V6+O5ubrp0+XHjny2yTc9QF5eOz79nXw\n9ZWYm9deuVJ+8mT+jh1WLi4009W1HNdQWlp65MhjkyZRu2kXpmZmJT//XHLokLmdXcft27mp\nqb0ff5xmdtelb7/N3ry5pbraIzy8+vx576ioqLffdg4KopbdpRtAr5CQKz/+KBKJ7Pv0oXNr\nnpbjDi1cOHzBghvZ2cNmz7554UL5yZNBFLMy7n4fSiwsjPw0j3Bgxe4PYbvkrnswonnTLgQ1\nA2Fz5zr3709hF/LuhaJ+8fEjli83s7Gh35XSIyzMIyyM1RYkc2ITE2sPD78xY0JefVXV3n4j\nO7syI4P+JMgcHe19fV0GDBjxwQf1xcX068FWyeU533/f3tBw8dtvAydMYJJbRkSi0atXn/74\n4xtZWRe+/to7Onro669Tu/k9N4JpRXVC2IlmvmgqiPch3AMOTwidluNyU1Pzd+6clJbGeixs\nsJ0BTqOpzcurkssrMzIUly8zrJ3WzTlEI2s5f+Gbb86vX2/l5uYZHs5vxZrb29OchFaFYu+s\nWU0VFYQQB1/fUatW2Xp707m1Tpcyfg7+/oOmTaP5Vmyrq1O3tdl4ehKRiNNoOJWK8kOOsrEx\nZeTI4QsXXkxOnnboUG5qam5aWuKOHcYzAEII0WrTJk1qqqiwcnPzHz+ezwyJXLz4Yd9WV3Xo\ndlXVxeRkRW4uq/ch9ABbsfdRX1x84M03c1NT6woK1EqlhYODqUxGcwBVGRnnv/iiVaEQm5ra\nentL6N6dGP0MiMRiK1dX5/79G69dkzk6RvzjH6q2NmVjo2Ng4MPeCO6im92fvn1pDoA5VVtb\nxalTGpXqVkFByeHDl7ZurTh9Wmpl5eDrS2cAZz7+uKmiImL+fKeAgGsnTypycmhW2+HZeHgE\nxsUNnDLFLTjY3MGh+fp1ym/F3NTUw4sWFezZ01BSounstHJxofwHgflGMPMBEEIyPv9cq1Zb\nurpSzgy5deXKxeTkAYmJRfv31+bleYSHOwYEOPj6tt26Rf9PItwLVuzug3mBCeYP6JgBQsjJ\nDz8s3Lt38PTpw2bN+m9c3IDJk4OnTqW5WiaEYitscSrVfydM6BcfP2zWrFaF4ufFi03MzFoV\niubKSmpvyOQRI8LnzQtKSCCE5KWlnVm1amZGBuWTE5xGcys/XyKT2ffp8+s7kO7CbcWZM0eX\nLiUikaazkz8d7Ny/f/DUqb6jR1MbQ+O1a/xGMCHEOzo6dtkyMxsbandnPgD9TQz6i6b3qjo0\ncMqUgPHjqQ0DeoAcu54IocAE28wqzAAhhGi1xQcPRsyfPyAxkWi1zkFBOT/8EDxtGs0hME+p\nYU7d0dF265aqrU3LcZYuLr2GDq0rKJi8bdvZTz8t3LuXQmCn5Tius7O+uJh/qWxstHByohzV\ndb8XTDGq41Sqkx98MOD557uE178sXapqbaXwf0HXmfSJlSvri4tdBw6keSJVv+3H0//5D5Od\naMI6v40/IRGUkNCl6pBaqaQ5DOgBArueCKfAhMTCwjsmhn5xB8wAIYRTqzUdHR3NzfzSiGNA\ngCI3l+YAiACKrTA/wSO1suo7alTOf/9bJZdbODpWZ2X1e+YZsURi5+PTfP06hQGIxGK/sWPz\n0tIUubnmdnbXz50bOnMmhfvqk3/xhVqpjFm6tLWm5mJKyrGkpPjkZJoDYB5ep69dW5me7jt6\n9O4ZMyh3JhVO2w+RWGzVq1fbrVvyzz+Xf/45/U0MtVKZv2OHyMTEIzTUOzqayZ9l6AECu54I\noQZmd8Oi94COGSCEiCWSviNHZn31VdmxYxZOTtWZmQMSE2kOgBD25xAtnJzyd+5kfjTbzsen\n+MCBmpwcvvfDntdeu5GdHUsrUTjqnXcsnJ1LDh9uVSgGPv88/Vqs5adO6faCZY6OZ1at0nIc\nzVVDtuE1286kMkfH+C1b+LYfxQcOFO7ZQ+60/aAc2DHfxEhfs+bK9u38z5YuLs5BQX5jx/Yd\nNYrmGKAHCOx6wuTzw6nVRCQSSEEgVn9BBDUJhJDHk5KsPTxKDh1S3LzpN3bsMIpBlY5d796s\ndn8EUtrARCodNmvWsFmz+JdqpVLZ2Dhs9mw6mT26NUvdACgTwl4wYRpes+1MKpC2H4QQTqWq\nLy4WSyReUVH0V8s4tbpw376wefMqz561dHUt++WXqvPnQ2bMoDYAuC8EdvfUpboEtc9PxZkz\nx5OS3ENDPcLCPCMibL28GJa00J+E0NmzqTVcF9Qk6L7Rw+bOZTUGHbGJCf14l/nWT31x8fFl\nyzSdna7BwR5hYXyhE1Nzc2oVcISwZsl8L1g/tGUSXls4OZnb2/+ydKljQIBHWNiht97y4f8s\nU+Q2aNBvu5ALFtD/u9RUUbH/jTduV1cTQhx8fUevXm1D989jS02NWqm0dHYWS6XBU6dKraxs\nPD1p1uiG+0K5k3vqvrb4w//8tNXVtdbUVGVmlp84kfe//xXs3l1fUqLp6KBfZ4SwK7EhnEkQ\nQqF/5piX2Dj173/XXL4stbKqkstLjxy5tHWrpbMzze8S5lWHeJ7Dh4tEoltXr7bcvNkvPn7Y\n7NnUVuzu9UEQm5o+NmlSr5AQCn8bRSYmXsOHm9vZhb3xRtZXXzVfvx67fLnU0vJh31ff2VWr\nsr/++vqZM3nbtuXv2nUzO1skFtvTKrhDCDm5cmVzRUXYvHlugwdXnD5dnZXFb81TY2JmduXH\nH83t7Tubm28VFFSfP29ma8uktR3cC8qddI95dQlOpbq6Y8e5Tz7hOE63SEa5zggmgRBSmZ5+\n8oMP2hsahr7+OpNdSP16sJRvrcO2xIaW406uXCm1stI/mt0/IYFet2Jh1Nxhi/kHgfx+A8F9\n2DArNzfKnUk5tTo5NjZk5kzdLqTIxGTCxo0P+xmjs7VV09HBz3napEn94uL4U/mXUlIy1q17\n5cwZmr0NCSF1hYWEkMr09Ix168QSyYSNG12Dg2kOAHqGrdjuCaG6RPbmzYNefDF0zpzG8vIj\nS5YQjqNcZwSTQASwC1mwa1c37RYoYl5ig+HR7PPr1/uPH2/n48M8Xb2+pOR4UpJGpdLfjKY5\nAOYfBPL7tveZGzbQb3vPahfy2tGjx5cvd/D39wwPt3R2risq4ovJtdfV2Xh5mdKNbtUdHZqO\nDlNz8+CpU/s88YSJVGrp4kJzAHBfCOy6x7y6hLqzs72+Xq1UajnOzsfHMyKiNi+P5hIFMeJJ\n0LXNsfH09AgNZfuN7hgYKLW07GxpKTpwIH/XLkK9HizzEhusjmar2toufPNNe0NDwPjxx5cv\nn7xtG8OaO5n/+U99cbGtj0/B7t35O3cS6kvXrsHBfZ980js6WiQWtyoUVXI55Q+CEDq0Wjg5\nmdnY1F69KhaLL23dWnPpEp2GEw5+fo9NnlyVkXH5u+/439QVFFi4uChycsauWUPzKaupomLf\n3LktN24QQuz79h29ejWiOgFCYHcPTKtLEEKklpa6sgIyB4cbWVlBzz5LcwCEGO8kNJWXV2dm\nmpqbp69dW37ihMedpTIGB9BY14MlAiixwWq1TCKTeYaH5+/Ykb9jh1giOfnhh3zpsl93IR/+\ntyn/gOFx56b9J06MXLy48/bt6szMSrmc8tL13RUuzKytaQ6A+QaCWqk8v369patr4IQJV7dv\nv7J9u1gi8X/qKQq3dgoK4tcFWxWKyvT0Krm8Si6vLykhhDReu+Y2ZAiFMfDOf/mlur096u23\nO1taLiYnH1my5NnvvzeqjtWPBAR298SwugQhhFOr7Xx8CCH1xcUmUqnvmDGhs2dTHgNhOgl8\nPk1HU5O5nV1bXV1bXR21SXAPDZ1+9KjYxMTBz6+usLBo376C3bsJizJ+zOvBCqHEhn411Mjo\naCISUYqtRaJxX3yhyMs79e9/65cuo/Y24B8wrp04QQiRWlmZmps79e/vNXx47xEjeo8Y8bDv\nrk8IFS6YbyCc/eQTvq+go79/xdmz/Z99dsirr1Jer7J0ceHL5mk5rr64uEou94yMpDmAxoqK\nwS+++NjkyYQQsalpxrp1qvZ2ms0/4I9AYHcfTKpLEEIub916KSXF1MwsdO5cRW5u79hYKd3n\nY31MJkGXT8PPAOV8mnu1zaG8/YR2C4RtNVSRyGXAgCEvv+zg52fp6nojO7syI4Pa24B/wKjN\ny6uSyyszMhSXL7OqEy6IChdsNxDu6itYfupU9Lvv0hvA74nEYseAAMeAAJo3bSwrs3Jz65rh\nR33VA+4LgV03mJ9DZJ5Nghkgd3ZeOpqaYpcvt+vbl1XbHLbtFjLWrRNLpSEzZrBqtyCEtSJd\nEOkdHf1rWQdabwOxiYlrcLB93751hYUDp0yx9fHJTU3VdHZSfsBglVumj6+ix3ADgXlfQeby\nd++uOHWKEHIrP9/S2bm+uHjMmjX0S2TDfSGw6wbzc4jMs0kwA0Rv54VotXtnzRoweXLwtGk0\nI139qryhc+fy/xdo1oPVVeWdlJbGqt2CINaK7kb3gYdvkPp4UpL8yy/5Bqn9H37pst8y/MLD\n7by9n96wgehVuKCTW6ajXyCayQaCEPoKchrNtePHxSYmbkOGmNvaUr47IWTg88/b9+lTlZFR\nKZc3lJYSQvbNnesRGtovPt6L7o4w9AwFiruhamurOHVKo1LdKigoOXz40tatFadPSy0tHWhV\noTQ1Myv5+eeSQ4fM7ew6bt/OTU3t/fjjNB+RMQNEq/3lvffC580b/PLLhJCbFy+WHT0aPGUK\nvQHcoyqvy8CB1OrBCqEqL/NqqJxafTE5+VJKiqqtzc7Hx0QioXbrXwegUh19//3IRYtuV1Xd\nyMryHz++5NChIa+88rDfALeuXLmYnHzt+PG8tLT8XbuUjY0SS8s+I0b0f/bZ4KlT6RQq12Hy\nVmxVKI4nJTn4+0stLcVisXd0tFarvXX1anNlZd+RIyPmzxebUl0Zyfj883OffML/QS4/dep2\ndbWlq6u5nR21AUgtLZ0CA/s8+WTw1Kl9Royw9fTUqFTVcrmDvz/q2AkKVuy6EsI5RLbZJJgB\nIoCdFy3HmdnY9J84Ub8qr5HU8NPvFGxqZsZ2rUi/dtrJlSvp105j1SC1S4YfwyNEhNFbsdvT\n8QOff55Vr+TyEyc8hw/vFxd348KFyvT0i8nJvUJC+DN2lOky/IKnTdN0dnJqNf0xQA8Q2HXF\n9hxi4d69Nh4e59aseWbLFlbZJJgBIoCdF4ZVeXVY1Rnp0inY0d+fj63pV0MVQq4nwwapfIaf\na3BwyIwZqrY2ymdH9DF5KwrndDzPrndv1+DgvqNG8UmfLTU1soefIaP/lNUtE6mUcv8PuC8E\ndl0xPIfYVlt7fNkyQohYIjm5ciVfNKuuqKi1piYwLu6h3lofZoD3eFKStYdHyaFDips3/caO\nHUZ3yZBVVV59uqbvlKvySq2s3AYPrszIuHb8OCHEytXVIzzcMyLCOHM9xRLJ0+vXlx45EjBh\nQubGjer2dsrnV3gSCwv6Z0d01B0d9YWFpubmXlFR1N6K+lmGQX//u6q9ndXpeEJIQ2lpU0XF\n9TNnmioq+M+Clasrhft2ecqy9fJC1TrhQ6/Ybmg6PVgK/gAAGoxJREFUOy98803xgQPtDQ0e\noaF/e//9n995hz+HGPD00w/vvlqOq83LO/XRR8rGxvb6en59WyQWB8bF/e399x/efe+GGWBC\nv+OFiURCRCL9JYrHnnuuP8Ui1VqOy01N5U9OULupPiF0CiZabdqkSU0VFVZubv7jx/O10yIX\nL6Y3gN9rq6szs7Y2tgWSu7sd2Hp7U7hv9fnzhxct6mxtJYRYurh4RkTwJ8lkDg40eyXzDi9a\nVHHmjJ2PT0NpqZbjCK3PQnVW1uWUlOrsbHV7O2H6lAV/HAK7+1MrlTumT/cdMybk4ScsE0Iq\nMzJs3N1ljo43srOr5PKWmzcjFy1i27YFM0AH/0Xy7A8/XP7uO/2OF7/m9ND9LmHe9J1Tqb4f\nP75ffLx+p+AnVq6k3Fiv8do1PteTEOIdHR27bJmZjQ3NAcCRd96pzsoa+tprfLcDG09Pat0O\nOI1Gv44gp9EQRsvnP02bxlc70mXc9k9IoPZZEMRTFvxh2Ir9Hf0CE7pO26bm5tQWLdRKZWNp\naVN5uUdoqHdUFM3TfzzmM0AIcRs0SNdpwHvBAmrRjK7jO53b3U1QOT3Mm74zb5fc2dJyYsUK\nx8DAJ1aurC8udh040Dgr7HMaza38fIlMZt+nD5NtOIbdDgSSZdhYVmZua6vp7CSESK2tmWTc\nZm/ePOjFF/Wfsiif5YI/DoHd72Ru2NBQWmrt7s6q0zbLIvuEEAHMAGE0CXd3fKdcy4AILKeH\n1ckJHebtkvnqcb6jR++eMYOvHjdixQqaAyACCKpaFYq9s2Y1VVQQQhx8fUetWkVnG1RHCN0O\nfsvwi4xkkmWYv3t3ZXp6ZXq6IjfXa/hwj/Bwp379aFYGZv6UBX8KArvfMC8wwbzIPvMZIOwm\n4T4d36nQ7w2qy+kZ8vLLkW+9RfmLhG/UW5OT4xMT4zd2LMOuG7ZeXiWHDtUVFVFul8ypVMUH\nDkQuXtxYVtbZ3ByUkFC0b9+I5ctpRlfMgypCiPyLL9RKZczSpa01NRdTUo4lJcUnJ9McAPNu\nB91n+NENsh+bNEnV2lqTk9NQVnYjO5t8+SXl523mT1nwpyCw+w3zAhPMi+wznwHCbhJaa2tN\nzMz8x40rPXLE3M6Ofsd3IqTKYd0Xb6O+YmQilYbOmRM6Zw7l+xJ21eP0MQ+qCCHlp06Fz5sX\nlJBACJE5Op5ZtUrLcdTiKk6tVt9ZJ26urGTS7eD8l1+q29uj3n6bz/A7smQJtQw/XqtCsf+N\nN/j43r5v3+glS1pu3PB8+P/tXc5yMXzKgj8Lgd1vmBeYYN6QkfkMEHaT0FRefiMriz8B16pQ\n9B05UiyRaDWa+uJimluQQsjpEULxNuYYVo/TYRtUEUK0HMd1dtYXF/MvlY2NFk5ONAdweevW\ngt27Tc3MQufOrcnJsfPxkVpaVsrl18+ccR00iE5gxzDDj9clvr+8dSud+L7b+szPbNnCqj4z\n/HHGHtjpP5R4hIayTSpiVWRfNwlugwfHbdrENrWL1SR0WS27duyY7gQc5YSetro6dVubjacn\nq8phzIu3qZVK3ekZu969meSWMa8exzyoIoSIxGK/sWPz0tIUubnmdnbXz50bOnMmtbvf/YBx\n/ezZv6ek0Ox2IIQMPybxfatCkZuaOvT11zUqlRDOcsGfYuyBXbcPJR6hofSTivivcyZF9rud\nBCapXQwngdx7tUytVNIZAK9g167z69dbubl53nk3mtvb0wxuzG1t7Xr3Tl+71srNLfOrr/ji\nbdTuTgRwhIgQwmk06s5O3zFjbDw8Iv7xD/rV49gGVTpR77xj4exccvhwq0Ix8PnnaUa3PTxg\nUOt2wDzDj1V8r8v3Pffpp3y+r/+4cSITk/riYiYP/PCnoI4d4R/Frv70U9H+/YqcHFaVii58\n8003X+e0YBJ4XQ4N/LrhQreAXMWZM0eXLiUikaazky9w4Ny/f/DUqb6jR1MbA8PibZxanRwb\nGzJzpu70jMjEZMLGjTSTTYVwaoHcqRNecviwqq3Nd/To8Pnze+js9DAwPpMrgOrQrQpFZXp6\nVUZGpVyubGgghEitrWlm+BFCTn7wQf6uXc79++vi+6GvvUbhvj3V8HvqKfSfEDKjDuz0q0vY\neXuzLfTP6usck6Dv4pYtWZs2iUSi0LlzFbm59Du+cyrVfydM6BcfP2zWrFaF4ufFi03MzFoV\niubKSvp1Z5g06m2uqkp95pkRK1YUHTgQPm9e/s6dNp6eA194geYYjiUlVWdmhsyYwWc1OQUG\n0j+1gEInhHV1aF1SBBGJtBzHL1ZVyuU3s7OHzZkTPGXKQ7pvl/aszON7Qoj+DgblLwV4AEa9\nFdttdQkm+7CcSnXygw8GPP98l6/zX5YuVbW2PtSvc0yCjhAODag7Otpu3VK1tWk5ztLFpdfQ\noXUFBZO3bTv76aeFe/dSDuzEJib0v0KYHyEiAji1IISgivmZXL5V8dP/+Q+TBwzSXVJE8LRp\nFDL87m7POmzWrGGzZj28O/ZAP+E1cuFCIhLRL3sEf5ZRB3bCqS7B8Osck6DD/NAAIURqZaWr\nF2Xh6FidldXvmWfEEomdj0/z9euUB8MEq9MzOkI4tcA8qCKso1stx1k4OeXv3BmxYAGTBwxC\niHd0tFqprJLLC/bsyd+1i9yp1v6wM/ykVlZugwdXZmRcO36csG7PKoSEV/izjDqwI8KoLkFY\nf51jEnjMDw3wRixfbufjU3zgQE1ODt8dcs9rr93Izo5dtoz+YOjj12kiFiygf3qGJ4RTC8yX\nDJlHt1Vyec7337c3NFz89lsmrYoJIY4BAY4BAaFz5igbGg7Mny82NaVTrd196FD3oUN17Vlb\nFIqC3bv5523K+RjMa+bDgzH2wI4wrS6hX2yF7dc5q0m4uwYmy5hGJBq9ejWf03Ph66+9o6OH\nvv46jfv+nolUqr/zolYqlY2Nw2bPDhg/nv5gKNNfpyGE2Hh6Uru1frJp5NtvszoKSgQQVBEB\nRLfMWxU3lpW13LzpNniwqUxmbm/vHBTUcfs25SZazNuzMq+ZDw8GgR3L6hJ31xkZNmeOR2io\nzMGB8tc5q0nouQYm5UnobGk5v369+7BhI5Yvryss9IqMpPxtWl9cfHzZMk1np2twsEdYGP9/\nwdTcfFJaGs1hMMRwnebuZNOhr7/OfxipjYHHPKjiMSx0Qhi1KtY/tZC/e/flrVvFEolbcLBd\n796F+/f70t1/FEJ7ViEkvMIDMOpTsTy2JzEFUmeE4SQIZAYIISdXriw5fPjxpCT5l18yafp+\neNGi62fPWru7N1VUaDmOUN95Ya5KLs/66ivd24DyOk1P9R3ovhUZHoQs3LvXxsPj3Jo1z2zZ\nwiSzjafbkf/tVw+/6tC1Eyd0pxbs+/Rprq6ulsv5QiemMlncV19RXqw68u67pT//7ODnp2vP\nGvX229Tu3nLzpqlM1lpTQ/QSXids3OgaHExtDPBgjD2wY1hdQjh1RlhNgnBmgBDCqVTJsbGR\nixe31dbmpqb2efLJon37Xjl9mlqlCS3HnVy5UmplNXzhws7bt6uzsiozMvonJFB+RhcChm8D\ntVJ5fv36jqamqHfeKdizp7mykuYYmAdVbbW13z31FCFELJH4jR3rGR7uERZWV1TUWlMTGBdH\nbRhajstNTc3fuZPyWnV1VtbllJTq7Gy+QS1/asEjLMzGy8vOx0dqZUVzMIQQTWdn9ubNJYcO\ndba1eQ0fHrV4sdTamtrdLyYn5/7ww9+SkryjogghzZWV9BNe4cEY+1Ysw5OYwqkzwmoShDMD\nRABN36vk8rKjR01lMsfAQK/hw3vHxvaOjaV2d+b0YxqJhYV3TAyTt8HZTz4p3Lt38PTpEpns\n8nffDZg8eVJaGp0xtNXWHl+2jBAilkhOrlzJJKiSOTrGb9ly6qOPlI2NxQcOFO7ZQwgRicWB\ncXE0AztWO/LCObXAM5FKQ+fMCZ0zh/J9eW5DhhTu23dw/vyghISIf/yDZsIr/EUmy4zjqN29\nmEilDWVlBbt2XTt27NqxYyWHD3uEhfk8/nhzZeXtGzce6v6LtYdH8PTpXpGR1r16dTQ3Xz99\nuvTIkctbt147flwikzkGBDy8W3fBahKEMwOEEFNz8ys//VS0b599nz6DX3wxfe1ar4gIb4p9\n329XVzeWlTWVl5cdPfrbJBjHcl1bbe3umTML9uxRNjberq5WK5UyB4eay5dvZGc79etHbxxa\n7S/vvRc+b97gl18mhNy8eLHs6NHgKVPoxPemMpl3ZGTt1asikUiRk8N/FooPHjS3s/P5298o\nDIAQIhKJLF1dbTw9g194Ydjs2W6DBskcHCycnMLnzZNaWtIZAyHkdnV1XWFh682blenpl7du\nLTt2TGJhQe2zoOW4w4sWDUhMjNu0yXfMmBsXLsjs7QckJsocHekMgBBSX1x84M03c1NT6woK\n1EqlhYODqUxG7e48Kze3oIQEqbV13v/+V7h3r3NQkJWbG+UxwIMx9hU7wrS6hEDqjBB2kyCc\nGWDb9F2tVDaUlPiOHh2xYEF7XV3V+fNG1ZBRIAtFnFqt6ejoaG7m07kcAwIUubnU7i4Si10G\nDgyfP9/G3V3m6HgjO7tKLm+5eZNO/yh9nuHh/A+/HZCni8nJCR0hnFrI3LChobTU2t29YPfu\n/J07CfUlw9vV1Teys9UdHSKRKDAu7sqPPx5etGj6zz+jk9gjwdhz7O4+h2gqk+2YPt13zJiQ\nV15h9iam2Jy025OYaqWS8STQbc/aRVtdHeWm76c/+qibKqAjRxrVn9Gi/fud+/e3dHHRxTSR\nixZRzun5ZenSksOHHfz9LZycqjMzByQmRsyfT3MAIIR+zWxPLQgh3TZj3bpLKSn8z1Jra5m9\nvU9MzO/OsoCAGXtgh3OImAHmhND2XgguJie3KhRRixczDGfVSmX211+XHDrUcft279jY6Hfe\nod/Jysgx79dMWJ9aqExPP7JkialMFvbGG17DhzMpznyroIAQIrO3l9nbiyUS+gOAv8Kot2K1\nHGdmY9N/4kT9ByPKFSDZwgwIAaqA8jpv385LS7N2dw+eMqWuqKj44EEHPz/KdUZMzc3D5s4N\nmzuX5k1BRwj9mgnrUwsisdjBz0+Rk3P8X/8ijAruOAUG0rwd/N8y6sDOyM8hEsyAMKAKKC90\n7tzbN26kf/bZ1Z9+aqqoIIT8/c5mEBgJIfRrZo5tiiEYAKPeiu1SDZXJgxFbmAGG9Mvc1xUW\nEuOuAqrIyzv14Yd1RUV8bZGhr70W9OyzFhTPIYIgaLVpkyY1VVRYubn5jx/P92uOXLyY9bCo\nUiuV+Tt2iExMPEJD7Xr3JiIR25xjeOQYdWDHY1gNVSAwA0zol7n3jIiw9fIiIpHRVgHNWLeu\nMj3dd/Roj7CwI0uWqJXKpzdssO/bl9oAOI3mVn6+RCaz79MHX6KstNy82VJTk7lhw42sLEKI\nd3R07LJlZjY2rMdFVfdHqeg2NINHmvEGdt08FRHGhzHpw6MhQ92WufeMiPi1Ua+RUTY1mdva\n8j83lpXtfOWV8HnzghIS6Ny9VaHYO2sWv//r4Os7atUqW29vOrcGfbpuB54REZxKZYQnV3CU\nCv46482xS1+zBk9FmASGhFbmniFOrb66fbuuwoVdnz6jV6+27tWL2gDkX3yhVipjli5tram5\nmJJyLCkpPjmZ2t1Bp0u3A9bDYQBHqeCvE7MeABucWl24b1/YvHnuQ4f6jxvX0dRUdf68sbVM\nwSQIRPbmzYNefPG18+cnb9/u4O/v4OtrbAeTL2/dmr15c7VcrlYqT65cWXb0qPvQodbu7tQG\nUH7qVMiMGUEJCcNmzx6+YIEiL48v/QOUuQ0aNPGHHyIWLCg+ePDHxMSbFy+yHhFtXY5SVZw+\n3XT9OutBwSPGSFfs8FREMAnCIIQy92wxr3Ch5Tius7O+uJh/qWxstHByEomN9KGXIXQ7IISY\nmpk9vWED0TtK5f/UU6wHBY8YIw3sUGCCYBKEQWpp2XfUqJz//rdKLteVuWc9KKqYV7gQicV+\nY8fmpaUpcnPN7eyunzs3dOZMmgMA3pUff9TvdmDt7u4TE2NUUV3GunVajYZv8OAYENDniSeM\n8ygV/EVGdHhCv7oEMdYCE5gEAWJb5p49AVS40HR2Xvjmm5LDh1Vtbb6jR4fPn6/7jAA1tVeu\n3MrPt3Z37xUSQrOhn0BoOS43NTV/585JaWmsxwKPNiMK7LqtLkEIMaoCE5gEEKDGa9dOf/wx\nwwoXqHXCHA4mV6ann/zgg/aGhqGvvx44YQKTTmJgGIwosEN1CYJJAAHjNBomFS4QUgjBsaSk\n6szMkBkz+IPJToGBxnYwuUu5eAd//0HTpqFcPDwAIwrseLrqEhzHkTv/7UZXXQKTAHAHQgoh\nSB4xQle5MC8t7cyqVTMzMozwCAvKxcNfZ4yHJ/jqEqFz5jSWlx9ZsoRwnLFVlyCYBIA7yk+d\n0oUUMkfHM6tWaTnOCEMKhnAwWUdiYeEdE+MdE0MIIUa27AL/V4wusEN1CYJJALgDIYUQ4GBy\n95DuCQ/E6AI7VJcgmASAOxBSCETUO+9YODuXHD7cqlAMfP75Ia++ynpEAI8qo8uxI6guQQjB\nJADcgVonAGBIjDGwAwAo3LvXxsPj3Jo1z2zZgkgOAAwGAjsAMDpttbXfPfUUIUQskfiNHesZ\nHu4RFlZXVNRaUxMYF8d6dAAADw6BHQAYHS3H1eblnfroI2VjY3t9PadWE0JEYnFgXNzf3n+f\n9egAAB4cAjsAMFKVGRk27u4yR8cb2dlVcnnLzZuRixah/woAPNIQ2AEAAAAYCJRrAgAAADAQ\nCOwAAAAADAQCOwAAAAADgcAOAAAAwEAgsAMAAAAwEAjsAAAAAAwEAjsAAAAAA4HADgAAAMBA\nILADAAAAMBAI7AAAAAAMBAI7AAAAAAOBwA4AAADAQCCwAwAAADAQCOwAAAAADAQCOwAAAAAD\ngcAOAAAAwEAgsAMAAAAwEAjsAAAAAAwEAjsAAAAAA4HADgAAAMBAILADAAAAMBAI7AAAAAAM\nBAI7AAAAAAOBwA4AAADAQCCwAwAAADAQCOwAAAAADAQCOwAAAAADgcAOAAAAwEAgsAMAAAAw\nEAjsAAAAAAwEAjsAAAAAA4HADgAAAMBAILADAAAAMBAI7AAAAAAMBAI7AAAAAAOBwA4AAADA\nQCCwAwAAADAQCOwAAAAADAQCOwAAAAADgcAOAAAAwEAgsAMAAAAwEAjsAAAAAAwEAjsAAAAA\nA4HADgAAAMBAILADAAAAMBAI7AAAAAAMBAI7AAAAAAOBwA4AAADAQCCwAwAAADAQCOwAAAAA\nDAQCOwAAAAADgcAOAAAAwEAgsAMAAAAwEAjsAKB7P/74o+j3bGxswsPDf/jhB61Wy3p0AADQ\nDVPWAwAAQXv66acHDRpECNFqtdXV1bt27XrhhReqq6vfeust1kPrydq1axcsWFBbW+vk5MR6\nLAAA9CCwA4CePPvssy+99JLuZWVl5cCBA5ctWzZ79mwLCwt24wIAgG5gKxYA/gRPT8/nnnuu\npaWlsLCQ9VgAAKArBHYA8Oc4OjoSQmpra/mXu3btiomJcXZ2trKyGjhw4Lp163QZeCNHjpw4\ncWJRUdHQoUO9vLzue/3YsWPj4+MzMzNHjRplZ2cXEhKyffv2zs7OhQsX+vn5WVtbP/XUU+Xl\n5bqRVFRUvPDCC71797a2to6MjNy+fTv/+9jY2AULFhBCnJ2dExMTe774XuMEAHgUIbADgD/n\n3LlzhBB/f39CyObNm+Pj4zs6OhYsWDBnzhyxWDx//vyUlBTdxU1NTXFxca2trWPHjv0j11+9\nejUhISEmJmbFihW1tbWJiYlRUVGZmZmLFi2aMmXKwYMH586dy1+Zn58/ePDgEydOJCYmvvXW\nWy0tLRMnTvzss88IIWvXruUv27lzZ1JSUs8XdztOAIBHlRYAoDvbtm0jhGzZsoV/yXFcdXX1\ne++9RwiJjIzUaDRarXbUqFHu7u7t7e38NUql0traevr06fzLJ598khAyY8YMjuP43/R8/Zgx\nYwghx44d41/u3LmTEDJkyBCVSsX/JjY21t7env95woQJ3t7edXV1/EuVSvXkk0+am5vX19dr\ntdo1a9YQQmpra//IxXePEwDgEYUVOwDoycsvv8zXOhGLxe7u7h9++KGvr++WLVvEYjEhZPv2\n7QUFBebm5vzF9fX1arW6o6ND96+LRKI1a9aIRCL+5X2v79WrV2xsLP9zUFAQISQxMdHU1FT3\nm7a2NkJIa2vrnj17XnrpJQsLC6VSqVQq1Wr1jBkzlEolv6Co749c3GWcAACPKJyKBYCe6Mqd\nEEJMTU0DAwPj4uIsLS3531hbW1++fPnMmTOXLl26cOFCVlaWRqPR/9d79+5tZWWle3nf6x0c\nHHQ/87Hj3b8hhBQVFRFCVqxYsWLFii4DvnXrVpff/JGLu4wTAOARhcAOAHrSpdxJF//+97/f\nf/99Ly+v+Pj4xYsXh4WFxcTE6F+gH5b9kev/ILVaTQj55z//OW7cuC7/qE+fPg9wcZdxAgA8\nohDYAcADun79+nvvvffqq69u2rRJt4nZZQXur1zfA/7oBsdxERERul/m5uamp6f379//r1wM\nAPBIQ44dADyghoYGQsjAgQN1UVp6enpVVZX2Hg3H/uz1PbC1tX388cc3btxYUFDA/0apVL70\n0kvLly/X31HlOO6PXwwAYACwYgcAD6hfv37e3t4ffvhhfX19QECAXC5PTU3t1atXenr6wYMH\n764b8mev79mnn34aExMTFRX13HPPeXh4bNu27eLFi//73//4PDyJREIIWbVq1bhx40aMGNHz\nxQAABgN/1ADgAUml0v379w8ZMmTt2rVLliypra3NzMxcvXq1Uqn85JNP/vr1PQsJCblw4UJM\nTMyuXbs+/vhjKyur/fv3T548mf+nEyZMiI2N/fLLL1NTU+97MQCAwRA9wCYIAAAAAAgQVuwA\nAAAADAQCOwAAAAADgcAOAAAAwEAgsAMAAAAwEAjsAAAAAAwEAjsAAAAAA4HADgAAAMBAILAD\nAAAAMBAI7AAAAAAMBAI7AAAAAAOBwA4AAADAQPx/7wWE6MXcrwkAAAAASUVORK5CYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      },
      "text/plain": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "a=ggplot(data = melt(percs),\n",
    "           aes(x=variable, y=value*100)) + \n",
    "geom_boxplot(aes(fill=variable),\n",
    "             fill=alpha('lightblue',0.5),\n",
    "             outlier.size = 0.1,\n",
    "             lwd=0.1) + \n",
    "theme_classic() +\n",
    "\n",
    "geom_jitter(data=all,aes(color=Sample),alpha=0.6,size=2,width = 0.3) +\n",
    "xlab(\"Parameter\") +\n",
    "scale_color_manual(values = c(\"purple\",\"gold\",\"darkgreen\")) +\n",
    "ylab(\"Percent\") +\n",
    "theme(legend.position=\"right\",\n",
    "      axis.text.x = element_text(angle = 75, \n",
    "                                 face = \"bold\", \n",
    "                                 color = \"#993333\", \n",
    "                                 size = 7,\n",
    "                                 margin = margin(t = 30))) # top margin\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<strong>png:</strong> 2"
      ],
      "text/latex": [
       "\\textbf{png:} 2"
      ],
      "text/markdown": [
       "**png:** 2"
      ],
      "text/plain": [
       "png \n",
       "  2 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pdf(\"/home/jrodriguez/scratch/ancient_assembly_hic/paper_test_folder/response_referees//illumina_4DGenome_comparison.pdf\")\n",
    "a\n",
    "dev.off()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**END**"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.6.3"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": true,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}