LucaCappelletti94/keras_synthetic_genome_sequence

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "from keras_synthetic_genome_sequence import MultivariateGapWindowsSequence\n",
    "from keras_synthetic_genome_sequence.utils import get_gaps_statistics\n",
    "from ucsc_genomes_downloader import Genome\n",
    "from ucsc_genomes_downloader.utils import tessellate_bed\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from scipy.stats import pearsonr, ks_2samp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Downloading chromosomes for genome hg19', layout=Layout(f…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Loading chromosomes for genome hg19', layout=Layout(flex=…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    }
   ],
   "source": [
    "assembly = Genome(\"hg19\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "window_size = 1000\n",
    "batch_size = 100000\n",
    "max_gap_size = 3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Rendering gaps in hg19', layout=Layout(flex='2'), max=25.…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Rendering sequences in hg19', layout=Layout(flex='2'), ma…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    }
   ],
   "source": [
    "number, mean, covariance = get_gaps_statistics(\n",
    "    genome=assembly,\n",
    "    max_gap_size=max_gap_size,\n",
    "    window_size=window_size\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Rendering gaps in hg19', layout=Layout(flex='2'), max=25.…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Tessellating windows', layout=Layout(flex='2'), max=385.0…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    }
   ],
   "source": [
    "ground_truth = tessellate_bed(assembly.filled(), window_size=window_size)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Rendering sequences in hg19', layout=Layout(flex='2'), ma…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Converting nucleotides to numeric classes', layout=Layout…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(FloatProgress(value=0.0, description='Generating synthetic gaps', layout=Layout(flex='2'), max=…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r"
     ]
    }
   ],
   "source": [
    "gap_sequence = MultivariateGapWindowsSequence(\n",
    "    assembly=assembly,\n",
    "    bed=ground_truth,\n",
    "    gaps_mean=mean,\n",
    "    gaps_covariance=covariance,\n",
    "    batch_size=batch_size\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "X, y = gap_sequence[0]\n",
    "synthetic_mean = np.isclose(X, 0.25).all(axis=-1).mean(axis=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAACUIAAAJECAYAAADZ38xZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAewgAAHsIBbtB1PgAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOydd7gkRdWH398GclqQnBZUEEVAJCuwoKKCCigqIMKC4cOEoJ/iZ1xUFAMoBsyyKlFQQBFECSsgUSVJVlhyzrDkPd8fp4apOzs90zPTM3Pv3fM+Tz+3p291dXV3VfU5p06dkpkRBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEwlpkw7AIEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRD0SjhCBUEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEw5glHqCAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIxjzhCBUEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEwZgnHKGCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIBjzhCNUEARBEARBEARBEARBEARBEARBEARBEARBEARjnnCECoIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgCIJgzBOOUEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQjHnCESoIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgjFPOEIFQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQTDmCUeoIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAjGPOEIFQRBEARBEARBEARBEARBEARBEARBEARBEATBmCccoYIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgGPOEI1QQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBGOecIQKgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgmDME45QQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRCMecIRKgiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCMU84QgVBEARBEARBEARBEARBEARBEARBEARBEARBMOYJR6ggCIIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgCMY84QgVBG2QNFOSpW36sMvTC5KmZvcye9jlCUYHklaR9E1Jl0l6RNLcrJ5MG3b5grGBpFnzY72RdGq659slLdzk/9HvVoykfbNn+pphlycIgvmTrB+yivOdL7+nML70rmB8IWltST+UdLWkx/L2L2nqsMsXjA0kzZ7f6o2kiZIuT/f8T0lqkmZa9lxmDaGY4w5Jh6Tn+ZSkNYddniAImtMvfaLNNcdNnytpenYvM4ddnk4Z6+UPQNLCkg6QdK6k+yU9l73TGcMuXzA2kDRjfqw3kv433fNcSRsUpJnv9Kd+IunlWT/12WGXZ1CEI1QwbpG0laQjJP1D0n2SnpH0pKR7kwHm2NTZbiop2kIwXyJpI+BK4FPABsASwDzGySAI5kXSDsAO6ecMM3tymOWpAjlrS9pD0uGS/i5pThXOXJJWkPRFSedLukfS05LulnSRpM9KWrlkVj8Hbkj7349veBCMLRocfWrb2zrM49tN8pjRpyIHQTAfI+ktwOXAh4GXA4sNt0RBMKbYF1g/7R9oZgMb7O8XyblrPUnvk/SjZHN8pgrHAklrJiemS5Id82n5hJu/Sfq4pKVLZnUI8CCwIPCdbssTBEGh7pJvc+VO0jdLOlnS/0haYtjlDoKgv0haCrgYOAzYElgGmDjUQgXBGEHSisAX089jzOzyYZanKiStJGlHSV+V9GdJDzTIDFO7zHehJF/8UdJt8skOD0q6StJ3Jb2qTD5mdg3wq/Tzs5JW6aY8Y41Jwy5AEFSNpLWBXwJbNPn3ZGAhYFlgQ2DXdPwWYOogyhcEowVJAn4DTEmHHgbOBu4B5qZjdwyhaEEw6knON99MP28BZg6vNNWQHCPPwh0iq857P9wg3xg1a/m0bQp8RtL+ZvbLVnmZ2XOSvoY/81cB78H7siAIxi57AX8ok1DSRGD3/hanGtJg6Nbp5zZmNmt4pQmKUDYL38xiQkDwApIWxQ2FC6VDdwHnA/cBtXrz6BCKFgSjHkmLAzPSzwvM7MwhFqcSJO0EHA0sUnG+E4Cv4hPUGm31K6dtK+D/JL3fzE5tlZ+ZPSzpe/jzf5ukLc3svCrLHATBCwh3kl4MH1vYEfiKpP8xs5OGWbAgCPrKN4BXpv3ngTNx+/Cz6dglwyhUEIwRZgCL42OQXxluUapB0p3Ain3Idxquf6zU8K8F8XHddYH9kuz/aTN7pk2WB+M22EWBg4D3VVrgUUg4QgXjCknrAecA+Syp+4BLgbvxjnVpfCbnWtSjoi01wGIGwWhhU+Blaf8+4OVmdv8QyxMEY4nd8W8JwGFm9twwC1MRi9EfJ6hDgAOzQ48B5wJ34gL7VsByuAL0C0kLmtmP2mR7DK4orQp8SdKx4+QdBMH8ylskTTGzh0qkfQN9MC4EQRA04a3UbQtXAxuPhwigQTAgPg68KO1/Y5gFqZClqNgJKnEUsFv2+wFcX7ofnzSyNbBk2j9Z0k7tnKGA7+OOVYviTlZbt04eBEEJLmVe54YJeN+wHvCKdGxZ4ARJ7zCzUwZYviAIBoCkSfikzBqvM7O/Das8QTCWkLQGsHf6eYqZXT/M8lRIP5yg3gKcRN2X5zngQnyljAWBjYG1cafsjwPLS9q9VRReM7tJ0gl4kJg9JR1iZjdWXfbRRDhCBeMGSZOB46gbKu8BPgqcZGbPN0m/NPA2YA/g1UX5mtl0YHrFxQ2C0cCG2f4p4QQV9IKZTRt2GQZMzbHnCcZBNKgGHgL+gRv5LsUF6kO6ySgJ7LkT1G+A/czs4SzNAsBngS+lQ9+XdLGZ/asoXzN7VtJPcWeoFwO74DJAEARji2twp9IFcCW8nRMkwJ5Nzh83zIff0xcIvSsYheT60rHhBBX0gplNHXYZBoWkhXBjPMCtQDunnbHGPdR1pUuBN1K/346Q9FFGOkF9C/hS3t9IWgw4FPggvuzOcZLWMrM7i/I1swclHYfP8t5K0mZmdlE3ZQyC4AVOM7MZRf+UtDlul1gNb6tHSDq9KDpDRCLtDTObyfizxwVjg7VwR2OA/4QTVNAL6bsyY8jFGCSfxFduAvjhMAvSB54CLqeuI9wDnNFNRpJWwiNB1fx4/gns1ui0JOkdwJH4BPNdgYuAw9tk/8OUdhL+PvbtpoxjhQntkwTBmGFHYJ20/xS+/MSJzZygwI0CZjbTzF4PbDSoQgbBKGJKtn/X0EoRBGMMSa/Hw44CnGBm42VJlCuAl5rZ0ma2nZl9zsxOxoX2bjk42z/dzPbMnaAAzOyZpPR9Kx2aSLlZ40dSX8azq4GHIAiGzrHUQ8fv2SohgKQlgJ3Sz8uBq/pUriAIAgh9KQi6ZQ/q0aB+aWZzWyUeQ/wZWN3MVjCzt5rZl83sdODhdic2IzmMfT479CMz+3Sj06WZPW5m/wOckA7VlrJox8+z/dCXgqDPmNmF+CStGisB04ZTmiAI+kjoCEHQBZKWoj757Sbg7OGVpnI2BBY3s83NbD8z+w0eualbDqS+cscdwHbNIjeZ2e9wp6YaX5C0ZKuMzex8oBaJ672SprRKP9YJR6hgPLFdtv8HM7u27Ilm9t8+lCcIRjuTs/3xYpgMgkGQr518/NBKUTFm9pCZ/aeq/CStjYeHr/HZNqd8FXg87b9eUssoL2Z2B/D39HMzSeu0Sh8EwajkfuD0tL+ZpJe2Sf9OYOG0/6u+lSoIgsAJfSkIuiPXl347tFJUjJndbWa3Vpjl1vhyd+CO4TPapP+/bH+vNJhUSIoAVSvvzu3SB0HQO2Z2KT64WyPsFEEw/ggdIQi6493Uo6md0GoJt7GGmV1mZs9VmOU7s/1vm9mDLa59GnBe+rkMI5fuLKKmoy3CSEeqcUc4QgXjiVWy/ZurylTSTEmWtukFaWZkaWakY5Mk7SnpTEl3SHpa0l2STk5LBXVShtUlfVvSNZIel/SQpCslHSxptZRmWlaGWT3edn7tZSR9UtJfJd0m6SlJD6ey/FBSX6JpydlN0unpuT0l6RZJp0raRZJSulnZfU9rkd+SKb+fSLpY0v2SnpH0qKT/SDpG0jslte0XJU3PrjkzHZsk6b2S/iLp9vS+75B0iqS3d3DfL5P0TUkXZWV8JJXxIklHSNopRUTomLzs1JehAvhSdk/z1PeCe54oadd0jzdJejL9fyeakN7pzpJ+JemGdF9PpXp1sqS95Gtsl72XhSV9QtKFkh6Q9ETK95eSNsnSvXBPBflMzdLMLnnt2dk5U0ukX1TShyT9MdXjOZIek3RjKu+2JfKY5x2k4zunfG9N9e7eVA/3qLWTskhaU96fnZvq71OprDeld/QxScsVnFuqLWbpK6sP/W43DddaFI9ACPAYXc5ckLSRpJ+ne58j79cvkfTZdI2y+UyWtG96/vemdniTpOMlvSFL11GdrYjNsv27zOzyVolTZK2/Z4feUeIaJ2X7e3RQtiAIRg+/zvbbRYWq/f854JgymRd9P1uk71guaJJHTe7YOjt8juaVtebRL1p9T+Uyee1/pcNrS9oqO++eZt9XSctJ2jt9ly+T9KCkZ+Vy/3WSjpT0xpLXa6YbLSzpfUlGuTV9q03SBtl5bfWulG6ypDemb/85ku5MMsSTcjn8NEkfly8lVJTHC7pTw/Fm76jlt1PSqpK+IOm8VJan0/O7TK7DrVXmuXVK1TLAAOrAIpI+nJ5TTb+7VdKxKiELZ3lvLOkHkv4ll5+eTWW9QdL5kg5L9WPh9rm1LjuwV/avI5vUi2lt7rltvW+49mS5Tvnb9C4fk+s4N6fntLNUXr6XNEXSF9P7fDjld62k7yk5nKtEn6cubA35cyqZvme7Q8E7qMwuk/JbT9IhcpvC3emdPi7p+tT23qeCGbhl22KWvrL60O9203CtNanrAdd3MkmxIZ9tJR2X7v0puc5/rqSPSprcPocX8llU0qfkemHXdoM+ketL/zSze1slTpM4a7O3JwNvLXGNk9PfBRk5oBIEQf+4O9svtO908a18haRvpe/6/embdqdcFjxQ0jI9l7z5dd+Y+swb5Hb0J+W2zZPksmPpPjnlt4WkX6fv4lPpe/p3uc1xsZRmnm96k3w60vPSOV3bPdP5lcnL/aSqZ5zS9ax7ZXk1lTslbSMfm/lvyvf+JJ/sJ4+eWOaeXyTpf+UyX62Mc1Jd/adcdtpb0spl8mtVduCc7F9b5225sT62uOfXyu3C18lt2Cbpuy2uv7Gk70i6XNJ9chn0bkl/S+2/o4gu8nGwP2XP6hZ1Me6mzuXbUjp/ll7qcfygxTuoxDaf8loita8/pmfyuOp99FmSviTpFQXnlmqLDedUUh/63W6asHu2f3JhqtZlXkXSVyRdIe97n0jt6PuSVu8wr+0l/U7elz0l/y6cJf++TE5pOqqzVZDuY8Xs0OlFaTNOy/ZjTCXHzGKLbVxswJ8AS9vxFeY7M8t3ekGaGVmaGcDK+GCutdh+CUwocf334hEyivJ5BB+Un5Ydm1WQ19QszewS1/4IHuK71X3MBX4BLFDhM5+CC5StrnsKHhpwVnZsWkF+b8eXS2yVX227DJjapnzTs/QzgRVwj9tW+Z6Gh0Zsle8MfBZgmXIe1eWznV4y/xH1vck9rwScW3DeTk2uu156tu2ueR3w8hL3sS5wY5t6+eWU9oXjVbSLdM7s7Jx29eWdeJjcdvf+R2DJDurdkng7aJXn6cDCJe5nQeAHJevfM83qMiXaYj/qAwNoNw3Xe1uW32kl0o+oX4DwZQyeb1HOm4A1S+S9Or6cXat7/iluJC9dZ9vUvVJtJJ336ey8i0qe85PsnFkl0m+Qpb+i1/cbW2yx9X9r+F7sCywAPJB+3wyo4Lyp+PfdgFPTseOyvGYUnJf3YTNLlG9Ev90iXaF8UfKbVNumt3g+0xr+t2L2/XgOWL7kM/9plufhTf6/X8qvTHnPApZpc70Z+XvBZ8L/uyC/DbLzZhY9lyzNqngksTJlvQ94Q0E+0zp8T1Ob5DEB+DLwZJtzn8WXim1at7tsR5XKAAOoA2sDV7fJ90hgcos8JzFSTmi3fbXLZzujg2tMa3HPpep9Q538T4lrXgisXOI+puEDsUX5PI33wVOzY037PErYGjrpI5ukrcTu0OQdVGmXWQr/5sxtk58BdxfkMTtL07QtVl0fGFC7adKf1PI7omRdfaF+4XJBuzL/E3hRibxfhcsWrerVQZ3W2TZ1r1QbSecdkZ13XMlzzsjOmVki/U5Z+lN6fb+xxTa/bYyUzWeUPOem7Jx9WqQr1e+kvvx7tJfVHgL2apPXiD63TdrlgDNLfDtuADYq8VwEfJfW39JrgZc19KtNnzsd6HlUY/esWl4uXf4O6mvVz7gS3SvLb2qWfjauq/yoTb7X0942vCPwYMly3t7ls51aMv8R77PJPS/Q4p6/2+S6U4ATS1zzIWCXEvexBCNliWbbHyg/7jY7SzO1xPVnZumnt0lbyfhBk3dQmW0+5b9vB/XvTU3On5H9v2lb7Ed9YADtpknZa33oo8DEEueMqF+4XNtKb5wD7FAi3wXxVT5a3fMF+Bho6Tpbot9o20bSeZs0nLdQiXN2y9I/SwvbSkov6rbY5ymhW43VrXTkjSAYA+TL+bxF0jrW5cyzHlkM+DPuqDEHd5C5DVgc2AZXIgD2xgW5bxRlJOmdeEdbi1L0HO58chMujGyNh9H+Le2XHOoISd8B9s8OPQBcBNwJLIQblNbFO8x9gJUk7WBmPYUDlXv5nwFsnB2+DTgff55rA1vgzghHlsx2OfzjBnA7cA1uFJ6Dv6918DVchQ+qnydpAzN7oETek3Hv5U3xD8bfcQedxYCtqHvuvhk4XdK2ZvZMYyaS9mdkhKb78eddc6BZGldS1gEmlrrr5lwL/DDtb0L9OV8KXNIkbTMWxAXiV+N18gK8/S2EP8cRSNoKd/KpReN5DvgHXv+fxQWC16bz1wYukLR5UfuVzzA9k3oIeYArccF4YrqvtfD1cMu8w74h6QDgULxugUcQuhCv0xOBl+PvQMBbgL9J2sLM5rTJeiLwO+B1uIJ+AfBf/BluCayW0r0JOAz4UIsyLgb8Bdg8OzwHb3O3p7KtjL/vZfA633UdrLI+DLDd5Lwh2z+/i/O/BHwx7V8OXIXf9wbU288awMmSXm1mzzbLRNKyeDSqNbPD1+Nt+Tnglfg7+wBe74ZBPkPduji/6UyZBq7AlaclgPUkrWBmd7c5JwiCUYSZPSPpt9QH5LcC/tYk6Z7U+5VfN/n/aKIma+2MG07A5cU7mqTtZDnvuySdDbwe/67tChze6hxJCwC7ZIeOapJsJerfyZtSme7DJxIshX9Tan3ytsCZkjYzs6dLFHsZXDdaLeV3HnALrhtt1uK8IhZNeYIb965O+T2OG5XXSPkuBLwIOE3S1mZ2QUM+d1B/Tx/Jjv+Q5jya/5A0ETec5TPt7gIuBu7FdYFNgRfjg1efBZYFPljmJlvRJxmgn3VgSXxSyJq43DoLX6ppaVw/rs1WnY6HZX93QT7fZuTzuwPXX+7H+4ZlcNl67RJlasUl1OvB63BZEnxA67qGtM3aNHRY75PefzT1ZTaewmXa2biOuRYuq09K518oaWMzu6fZxeVRbk6lHoXC8Lp5Db686Gvxga0fMVLnHzh9tDtUaZdZCW9zed16GNf978Lf2+p4u1s8lbtrKq4Pg2o3Ob3qSz/B+4O5eL29DreJbUa9nBvissD2RZnIo/H9lfo3A9xmcAXe370av/cvSrq/i3JWwSD0pfOy/ddJmmhmz3dxrSAISpC+wWtkh84rSlsyvwm4/e9t2eEHcXnqQfx7vg0uBy8FzJQ0xcwKI8uUvO7y+Hfuxdnh/+L98tN4/7lpOv5SPPrtm8zs7xRzOPCx7Pej+Pf1PlwW3QaXu07DJ35WQoV2z37Ky1VR9TOuSvcq4hu4HQB8AkHN6WVDXIYDl3vOTvbymxozkEcNPRFeGOt+krrc9DRuL3wx/n4WKVmuZjxKXUdYGXfIAJdZT2pIe1GLfL5D/Z6vwuWSZ/H7HCHfSloBf3/5EpvX4vbkx3B59rX4s18K+K2k95rZ0c0unGwDpwGvyQ7fgfdTc9J1NsMjTs5scQ99p+rxpAYqsc2ncn6PkW3ueVwvvxHvG5ZN+U5N/+9aT6iyPgyw3eRsS70PvagLefR1uJ4wEbcnXIi3yzVwJ99JuK77W0nrmtnNzTKRJHw8Pf+u3od/Vx/B39XW+Dfj94z0OxgUHa3y0oRJeJ9ydVECMzNJF+BjkhPw53t8j9cdnQzbEyu22KramHdW773AAcCKPeY7kzYen4z02q1FHpoJLN2QbhF8GY9a2seARQvyXJ6RHrkXAWs0pJmEf7Tz6xo9RoTCDYx5GfelycxLXHi+PUv76Qre49ey/J4DPkrD7GlcqbimyX1PK8jzrcBngJe0uO4auKG0ltfPW6SdnqV7Ov39F7BWQ7oJeESUfBbGQU3ym8TI2RWfocBjFx8w2LuiZ53X2xlt0ub3XJs9M4vms+MXzPZXYORM5GNoMmM11fffZ+mupMArHBf4aukeALZrkmZXXIB/OktrvbSLhnNmZ+fM8wxSmtdRn1nwDPA5mrR3XBjOZ8g3nTXb8A5q9f60xueZ6tO3srRzi8qY0ufRNJ7D+5Rm5ZyAt/mTaRK5inKzRCqrDwyv3VySXfMtJdLn9evp9D7+A2zSJO07U12ppd+zRb75t+QJmszywJWfO9I187ZQWB9K1L1SbSSdl39L7ix5Tt4PG+Vmeucz/d/W6zuOLbbY+rs1fC/2Tcc2z479ouC8WiTIh0izoRilEaEK7nVaF89nnnPwpbpq/7+0RH47Z+mvL0izDy5zt4oqsh5u0Kvl9fkWaWdk6Wpy4wmNfTouW0zOfs/MzptekPfq+Iz4TSiI4oIb7XJZ6IaitGXfZZNzvpyddw/wrmbXwJ3Q8hmL76qgDVUuA/S5DtSu/1cadHPcEPyd/B0AezfJ70VZXXoutYOi6HEr4sbo91XwrNvWyQrq/SvSe6yd+x0a7Agp3ZqMjELcNDIpPmnl2izdzTREacANq/vjukpeP2YX5DktSzOr5LMr00dWanegP3aZSfhAaS3tHNx5ch69Ax8QfCtwUkFes7N8itpiZfVhiO3mnqxc65ZIn9ev2nu7BHhZk3r78bxuAVsV5KmG53MPzb+pO+IDKCMiiHdxz3ndK9VG0nn5t+TCkudcl53zeMlz7sjOWa/XdxxbbPPTRgcRoXAHgluy9L9vk77MtzKPtG2448iCDWlWYGSEl2eBTQvyy/vcWS2ue1qW7glg9yZpNsKdo2rpbgWWKshvu4b7+AWwWEOapXGHkvx7UPjcKannUZ3ds2p5uVT5O6ir/XjGq1Oh7sVIXbtm/7yf5lFy3szIsbFzaCLHpPdVS3MiMKXg2gviDtQ/ruBZl2pHTe65Fg3nVmDLZmVsqI/5GMg/aRJ5DdenvkR9/OlxGsYQs7QzsvzmAp9sfFe4I1CtXed6wrSCPGdnaaaWeHYzs/TTC9JUOp5E/2zz+2bpDHciWaUg7bq4o2Kzcaz8vRS1xUrrA8NpN9/Mrvntkufk9eupdD970Hy8ONcXf9kizw82vLev0qDb4U7G52fXbVlnW1xrasO12raRdN6aDeetXeKczzScUyZC3MFZ+sN6fcejdRt6AWKLrcqNuiCXb3NxQ8Gv8RCmmwCTOshzZruOruFjZcAxLfJbCBd2amnfXZAuFyBvoUCRSGm/2XD9WQXp8o53dkGaxfEBpppw1tSwlKVfh/pyEPcDi/Tw/pZu+LB8okXalbNythTIOrj+ZOpLTDxJsQAwveG6t9Ei3C3uAGNZvo2G2HWz/58/wPaS19sZbdI23vOVlFty7RfZOT9rk3YiIwW6edoGI5W6uTRRGrK0724os3XbLpqcMzs7Z2qT/0/AFb9amve0yW8F6gbjZ2giNDd5B+dS0Jfhht/cYefAgnSvb8hz1x7q06x2bbHK+jCMdpOea778TaFzZUH9qvWTK7VIn/f9pxekeUVDnm9vkd/6zBv6e5462+Ye8rpXqo2k89ZruO48y8A0pF8CH4jKz1mrxHV+nqX/4iDqQmyxxdb9RhNHqHT8+nTsERpkDDwaaO2cn2TH50dHqMUYOUjesp9kZNj0L/T47pakHnXxToqd1mc09OVnUG7pqZnZOdMrqGv5kgNv7uVdNqkjeUj3l7VJv212jWvoYYk8BiwDVFgH/kWLcO7Az7K0d9Ag4+KzFGv/73mp4w7ut3Sd7KHen5Wd87k2aRelPiHIaDLIyUjD7hzgpS3yO7ChzLML0k3L0swq+exativ6YHdo8g6qsMu8P0vzLC30zxLPZHa7tlhlfRhGu8H12to159JiKcOC+mW4Hr1Yi/QnZGl/VJBmhyzNc8DGLfLbvuH6Tetsm3vI616pNpLOy5ddfwZYrk36FzeWteQzzpe2KhzQiy222ObdGCmbX4Ivr5ZvR+BR/BqXwz2LJs40DXm3+1Y22kgObZHXgoy0A55dkC7vc2cVpNmm4V4KJ5zhcnHu9N/UJsNI56BTKHbMncy8y9rOKEg7PUszsyBNZXbPDupMWXm5bfk7vG7lz7jD67fVvZjXRvo8sEWLPLdk5ATz7ZukqU2SfYoW8kPF77htO2pxz09Qztb53uycy9rdG77UW6F8hEcImpOlKZQzcXnj8YZyTytIOztLM7XEfc3M0k8vSFP1eFLjO6jCNj8FtwW0lElL1qcZ7dpiH+rDMNpNPvn5/SXPyevXXJo4TWZpc/n/MZqMm+GTXO7M0n2vRX6LM+9S5U3rbIs8Gute2zaSzhMjAwDsX+Kccxuu9cES5+yRpW8qN4yHrbbcVhCMF/bAjSI5wsMjvhf3ur0YeEjSMZK27kMZngE+UfRPM3sKODY7tHFjmhT2dq/s0Awze7jFNQ/ClY4q2AcXjMCF8HNbJTYPN/mr9HMZfCmubtmd+hJ2N9FiqQ8zuwMXSirDPMxlLVTkQvhM6jJ80Vovo/dN/KNdy/e9Df9fItu/r6e+vw8AACAASURBVOQ1h82BZvZkqwRp2Y73pJ+P4BHaCjEPh5kv8fieJsnen+0fb2aFIabN7Hh6DEHdA2/FQ0MDnGUFIWlrmC/l9Z30czIeVaAdB5jZcwX5GSOXjpynn0l8Mts/3syOK3HdruhDfRhGu1mBkSFsi5ZDacXXzOzOFv//ZbZf9N72zvbPNbPfF2VmZldQfhnRqrkKj0JQ4+A26T+PD/DnLF7iOvl7WKMwVRAEo53fpL9LUA8zX2PPbH+0L4vXV8zscXyJ4hrN5CUAJC2JG4NqtJRHSlz7Eeph/1fEl8Qow/7W4/LZXZJ//15fYb4fpx7S/Ztm1rhc2gjM7GzcKQbcmeNVPVx7qDJAD3Xgk0kPLuJT+KAA+LIdOzT8fyzqS23rvaT1cUc5cOePQ1qlN7Mn8AgyNdrpS4eb2Y0tsjyUkbLaIOm33aFnu0wi15cObaV/9kof6sMw2k0ui99rZs90kcdn0reuiE71paPN7NKizMzsNDxi3TCYRX0J08n4QFgrvt7kWOhLQTA4NsajAubbh3B7dm1JtPvw6JavS3JTL+xO3UZyD24zaYr58msfzQ5tI6nbZU//J9v/o5n9oSihmc3GV3eosW9adugFJK2LR48CH2z9RLJbNsvvWVwurIqB2T1r9CAvd82Qn3GNbnSvo6zFMnpJ7jomO/SBJslq8s6cNvLDaOEHZnZDiXS5HPuxEvf2derjg7ul8cWc3fFlw8AnA3yzKCMz+y9wWIkyVk6fxpMaqcI2/0HqMtgt9H+58arrwzDaTS6D3t7F+aea2Z9b/P80PJIY+LfzZU3SvAnvl8Gd/Vp9Vx+jvoTiQEn996nZof+VtHRReknb446jOaEjJMIRKhhXmNkTZvYuPGrMn/GZX81YDNgNmCXpFElTKizG+cmpoRWXZftTm/z/5fj6seCzDk9slVkygJ1ctoBt2D7bL6scnJ3tl3Ueasa0bP+31n6d2I4HcyQtJelNkj4p6WuSvifpB7UNrzs1NiiR5dP4mrKFJAUjF9q3aUhya7a/raR1GN08RH0wpxWvp+7YdmpJoepifHYCNK9L07L93zT5fyPDGiztdzu6ycz+2SZNy35G0oKMfJ7fb5Nfr1RdH4bRbpbP9ue0cwYsoNFZt5Hr8NnuAMtIanQMgpHv7Zgm/2+kp4HvbklCe24Q217SryQtlaeTtICkL9HcELNwk2ON3J/tr9B5SYMgGCX8BjfaQub4lL5X704//2tmfx90wUYhR2X7rQx9u1B34L3AzG5ql7Gk5SS9TdKBkg6R9P0GWXmjLHkZWfnK5MBQOZImS3qtpI9I+oqk7zaU9SMdlrUso0Vf6osM0Ic6cBs+2F9ImvRzSnaolb70Dkmj/Xtftt7ndemEEvovtKhLkhYHXp0daqkvpUkVZepRP+h3O+rZLiNpdUYa0H/QJr9eqbQ+MJx2k+tLrSaLFfEUI43+zWhnTwPIJz2OZn3pUXzJoRofkvQNSfnkGyQtJunH+HItjYS+FASji2WBX0g6qoLxhm2z/ePa2aDM7BJ8QlqNRnmqLPl5vyxMVedIPFIH+ABzowPWtGz/4uRkUUhyjJld4rot6afdsw/ycq9My/b78oz7pHuVsdv/Ktuf1uhoR13emSJp95LXHSZt5V5JK1J/hneY2fntzklO/hemn0viKynk5O36+DRe1Yphjan0YzypkSps8/mkjJ8lZ9S+0Kf6MIx206ue0PK9pXGIK7JDU5skm5btn5pk8VacRL0uDJpv4L4B4Ksj/UXSSxsTSdqZ5v1K6AiJScMuQBD0AzP7K/BXScvgBpAt8PVtN8Q7/py3AedJ2jx5efbKVe2TjOjoG8sDI4XFa0uW61I8pGqvbJ7t7ympcUZ+M1bJ9lft4dr5fV/SLrGZ3SLpPupOY4VIWgWfUbkLdWGqHS8qkeaq5IjWjguz/RGzwM3sdkkX4PV0CeAfko7GP7TnV1Qvq+TykrP687q0VlKKylAbBJ0iadHa85W0MiPf9cUl8rqo5DWrJr/3HSSVUQDzvqBdO6qqn6kZWOdQ7nn2QqX1YUjtZtFsvxsh+BEzu61VAjMzSQ9RF1aXxGcoAC9EDFwvO6VtX4mvG254hMSBYmY/l7QddaP9nsDOkmbhocKnAFtRV4b+iEdUq1HmPebvYtHCVEEQjGqSXHcuLru/QdIKaRD7bdSjhpRxgp4fOAOfbb4s8GJJm5lZM5lnj2z/qCb/fwFJL8cNLW+mHu2oHWVk5XaO2x0jaWF81ue+JctAB+naXXsZYK3s0AGSms62biCfCd6VvtRvGaCPdeDiohnpDVyIz1SGeaNmXYzPtF0df35XS/oVLjdc1KVzej8pW+9z+XjzkvJx/i4b69J61Cc8PgqUccYaDfpSP+wOVehLm2X7N5pZNzOXO6Hq+jCMdtOrvnR9iShSLd9bsv3kfVNhNKgO0/SLL+ODk1uk358G3ifpb/jgxPK4bFSThUJfCoLhcZCZzWg8KGlRfMD1zXgbrkU02UDSlmb2UJfXy+WhspNB/g68Mu1v2OkFk+11uexQYbSeGmZ2n6QbqDsPb4g7EtToyN6fuJRiZ9eyVG737KO83Ct9e8Z91L2Mcu/kYur6zFJ4mfOIpsdTjwh0lKRd07GzzeyukuUdFM9STkbNZUJ1YEN/cba/KnBl9juvI22fu5n9R9IDeCTWQVLp+EETerbNJzbN9s8pWcZu6Ud9GEa76VVPqHrcvcz485OS/k1xZLC+YWbXSjqA+mScVwPXpPGwG/Ax7o2pf3tvxJ2Sa87IneoIC0maMKRo7n0lHKGCcY35cmW/T1vNePxqfFDgA9Q/Zq/Al+zZr4LLlgl7m3tcT27y/1xYLGts62aJphEk7+Y8ZF7jEm5l6GW2S7f33dIRStKr8LXZOy1bmfCBt7ZPAvhM6BrNyrsPLjStCCyC188PAM9Lugpf4/V04EwrWA5tgJQNq79Str8x3QkMU6gvk5E/tznWejnCGi0F2z6S33sZo34j7epqFf1M7oV/2wDqVdX1AYbbbsoM6DVSNix6q3e3RMOxtn2lmT0u6RHqxvNB8x48NO1HccPF4ow03tc4Fh8IyP9XZtnXgTt4BUHQN36ND/ZNxPuOQ6lHhzLCEQrwKC6Sjqe+/MV7aHBmSAPBW6Wfz9IigqmkN+LReMpOFqhRRlaudDmmNLP+bDqfWV2mrGVYseH3h7vIo1t9qW8yQJ/rQM/6kpk9K2kP4E/4c1gaXybhAOAZSf/C5b4/AeeVdLzqJ93oS9syMvJDGRrrUv7cbi/5HAauLw3I7lC1vtQ2ol4FVFofRkG76Yu+lO6r9rOZXTu3K80xswdLXLdne1q3mNkzqQ/+BfCudHgZ4O2NSYHv4hGzavrSXNzpsR2hLwVBH0kD7lfjDqfH4k4GK+PjDYcxcrnOTsi/67eUPGd2tt+NA05+zSfN7N4OrlsbjG287lDGOajY7tlneblX+vKM+6x7PVQm2o+ZPdqgzyzLSEeog3Ebwmvw791b04akm4HzgDOBU0pEf+k3D5Wshys17H+kKGELWukJnehng3aE6sf4QU7PtnlJSzAy2k6/9YR+1Idht5t+jav0a9x94I5QAGb2Q0lz8IiGi+K6z1bUbX01/oUHAMmd8mJMJRFL4wXzFWY218wuNbOP4w5Reaj0DyQP954vU0EeebjFOYWpRlImKlE7mnnJdkovDpaV33cKg/s76h/7e4CD8DCIq+IfkAlmJjMTI5XTMn1kN+VcWNKIWSNmdj2wPvAdIDfUTcSVjf1wh45bJDVbD3uQlPXYrro+DatddEOv996uHVXRz+RK6SDWgq68fxlCu8nr0yJdnF/19wFGf1vAzJ41s/3wMLyH42FqHwaewZXv3wJvNLPdGRmG9Tm8z25HvnTE0O4zCIJKOIF6v7anpOWohxw/v8zSbvMReYSnd0tqlB12py7Lnl7kQC5pWXzmX82gfzNwIG4QWwn/3uWy8kHZ6WVk5aojjvyQuiH+aeAneNSwl+KyzaSsrGtk51Vl+ximvtQXGWAAdaCbcs4zeJKWAFgfX6IlT7sAHrnn08DfgOtSePhhMih9qTESwVjRlwbRjuZHfWmeyBRDaDejTV8aze3gBczscTN7Nz7j/yd4NLdH8e/cTcBMYHMz+wSwWnbqnSUd2EJfCoIBYWZ3MFJW2kPdL02a92dl225LeapP12x33WH1y5V9xwcgL/dKv55xP3WvsuVsLOuI+mVmc/DIigcAjUsCroFPrPo1cJekb1U0/tctwxpTgbEjH/VbT6haR4D+6wn9GFMZRrsJPaELzOxIPBLe5/FJLPfhDl/34I6q78f1ujsYOXGvjKNXriM8NR6jQUFEhArmY1JouU8CR6dDC+GenecOr1Qv0M1HoYrw1o2d+lJmVtZLugqewGcqQnX3/Q7qgvjtwEZm1mpgvVMlsZtyPmlmzzcmMLP7gE9IOhAPr7kVvq7xa6g/l5WAn0p6ZXIsGM3k9Wl/Mzu8x/xyoXKQ7aIZ7RS6J6gLqRuY2RWtEg+JPDxms7Wuq6bq+gAMvN3kzrOLSFp4CEuxNPbTizQ51oyhL4FgZtcA+7dJls+wuCqtad6OfFbTaAt7HQRBB5jZY5JOxp141sOXHajprL8eUDHGxGQhM7tY0o24EXpZ4A2482+N92T7rZbF+wB1meUyYKs2M3MHMaO5KWmpjl3Tz+eB7cysle7Wj7Lm39yHzayXaLi9XBuqkwH6XQe6kdubhnE3s9n4klEfxZeR2grYEjf81YyzawG/l/RJMzusw7IOmvz97WRmp/SY36jQl1Ik7lYM2+5QlmHqS1XUB2Dg7SbXlwaxFFAzhmVP65m0zG275Spzfanskn6hLwXBYDkj25+ER904vot8Hqcuo5Xtp9rKUyWu2SyvXq47rH65yu/4aNeZKn/GA9C9OnGEaFmvzexZPGridyWtS13W2RKPzla73v8CW0raZhQurZ2Tv8+TzawKZ/W8PxmmfNSJnlDZ+EHFNNbBxeivM1Q/6sMw2s3d1Otg6AkdYGb341G8Di5KI2ljfLILuNPYP0pkPV/oCGPCyBsEfeT0ht+NSx0Mi/uz/VVKnrNy+yStMbOHce/+Gi/tNc8O6cd9vy7b/04bJyiA1Utet8Zq7ZPMk+7+wlS8ED3lfDP7mpltjwsGb8Jnatb4WPq4jWbyZ11FXcqXmFhEUpnQrKuWSJOHzCzrINzOE7/qe+8HeRlXbRJJop/Xq/yZDKjd3MPI2Ttl+6kqeYSRdbZtGdLyI1XMHhkE22T7F5Q8J/8OzK6uKEEQDInc4Wl6+vsUHi2qGzr9zo+V/hLqEzrAl/4GIBmx1ks/HwX+2CKPXFb+aonlCTqVlatkW+qhu09rY4iH/pQ1l2eWSrPDB0W/ZIB+14F+6EtPmtlZZvYlM9sWX7LhHXjUyRpfTwM4o5l+6kurKls/rAX90Jda1rlRYHcoS/5+1ihM1Z/r9UNfGkS7mZ3tLytpgaKEfSTvPxZJy/q0Y7T3FTmhLwXB6KdxMLFbmTT/rpeVp/JrtZSnSlxzYUllB6tbXXco4xxUa/cc7TpTP55xv3WvKZLaOk+lZchy2bKdnvBvMzvCzHYzs1WAV+HLz9bYlO6WFhsk/ZAJu+lPBq4nMAbGVNJScfn4QL/1hL4/kwG1m9nZ/jDGVGB436NBkOsI15ScaDRf6AjhCBXM7zRGmni6aarBc3m2v04yYLejKqeYS7L9N1aUZ1ny+96kXWJJqzPSa7UZ+Rq6V5coQ+P6qu14paQyXsGbZfv/6uQCycHjDGA74KrsX2/tJJ8hcHG233NdSuGlc6F9s6K0HabJvfintBswkLQa9UhDRVR6733icup94CK4QNtPBvpM+tFu0rIDeV5rd5tXD2WY21CGtn0lvhTsqF/zOYWL3y479JuSp74s27+8MFUQBGOFM4E7G46d0kO0kPw7X8aJ+pVdXqcZVYTubkUe6WnHTCbdIzt+YpvoeqVl5bS082s6K2Kl9Fuub4uZ3YUv61pju6K0fbh2v2SAfteBzUo65PSiLz1pZr/HjX+1iDQLMHrl8BpVy8dXArVw9osDLy9xTqf6UlX96DDtDmXJI/OsJanfBvtB60uVtxszuxu4N/0UHmlqoJjZbUC+HGyZvnK0TzIDQNJGwDrp5/PAsSVPDX0pCAZLY5SJbpeauSzb36LkObmc1pE8BS/YXu/NDrW9bpqomtvHGq/bkb0/UUW/XKXdc7TrTP14xv3WvUS5d7IpdX3mYXxZwtKY2eVm9n7gp9nht3WSxxDIZdBXVCSD5nWkrfwv6SWUk/ur1hPGwpgKjCzntn2+Vj/qQ0v61G7yyRcDH1NJdDr+vDCwbv+KUyl7ZvsxppIRjlDB/M4GDb9vbZpq8FxN3eFjMvDOVonToMdOFV371Gx/X0kLFaasnlnZ/ruS4tCK97T5P4xUNluGO5T0ajpXtBYE3tUm38n4Ei81zunwGgCY2TPAX7NDy3eTzwA5A3gu7b9E0lsqyHNWtr9HUaKMvdolSF78D6afi9DeWNvyfSfydrS7pOVKnDNQzOxpRtbFj/b5kv2oD23pQ7vJB23W7zGvbpmV7e9elCijTF85GjgIqPX7l5vZxa0SA6SB1VwhuaQobRAEY4O0fPAxDYd7WRYvN5SuX8Iho8x3viy5A9LkCvMFwMz+S90otiiwU7q/3bJkrZbFgw5kZVzfWKGjQlZLJ3L9Iow0BLXihfeU5PZ2/Cnb37+kk09VzMr2q5IB+l0HVgGmtUogaSlgx+xQt/rSQ4yMkDLa9aVcZ3idpJ4cMc3sMeCf2aGW+lKKjFCmHs2m7tj5khITtTrVlwZtdyiFmd0CXJsd6nf0gErrQ1n60G5Gg76URwUeT/rSV7P9PyaHhZYkB4XaYPbjlBvMDoKgNzZs+N22rRZwdra/a7tvpaQNqUeFhS7lqYbzppdIvxf1ccY7gesb/j8r299U0pqtMpO0BTC1xHVbUrHdc7TrTLOy/aqecb90r5z3lkgzPduflSapdkMuZ41qHcHMbmakDLp/BdnmbeHdJfTetmMqidze0jjWOoLk0N0uetJQxg+6IF9p6AOSFuzXhfpUH8pSZbsZDTrCrGz/LSniXCt2pr6U+KhF0m7AK9LPp4EjS56aywxtx2HGKuEIFYwbJH1C0us7SD8J+HJ26B5Giddjmu2bD/jMSMbhIr4EtPp/J/wE964HN1ofUda4L+lFJZyXWnEM9ahcawL7tbjWysCnSuR5U7a/Y1GiJLD/tOj/bfhKm2XaPk1duXiahgEpSVMkle2P89Cl9xWmGgUko1x+rz8uG2Zf0oSCpUbykJzvlrRlizzeDby2VGFHCmLTW+S5CvB/JfL7HfCftL8IcFTJgTUkLVYyylgVHJbt7ypp18KUPVJ1fRhiu8mdqsrWr6rJhdmtJBWuDS5pfWDv/hepNyS9E3h/+mm06P8bWJ96hLYr0yz0IAjGPgfjzum17Ywe8rqW+izFFWkRwUfSDsAOPVyrkTwiRb/Caeff1j3wmbi1797tjBwMbkZZWXlZ4DvdFLBC8rLu0GZ5i0Mpb6jr9D0dikfiANgI18VKkaIf9kI/ZIBB1IHD2gzefQt35gNfTiZ3NqONrtXIWNKXLqFujBWuM7QzxnpiaYGCJb9yfWk/Sa2WT/gkJZZySBNHaoOak2jhNCLpVcAH2uXJcO0OnZDrS59spX/2StX1YYjtZrTpS+9ptTy6pO2BN/S/SL0h6ZPUIyM8iduZypDX2bOTw3kQBP3lgGzf6N4h6RhG6jGFMmdaivT72aFzzKzRIaksP8n2d5b0phbXXRX4fH5uo6OKmV1F3VFbuFzY9Juf5PtvdVXq5lRl9xzVOlOfnnG/dK+cPZJTVlOS3JVP8vl5w/8XLOGgX2PM6AiJb2T7H+9w3LOZznkM9eXcVqWFHCHpxYzsx1qRj6kUOk+l+nN4u8z6NJ7UD36GO5iDLwv53T5fr7L6MMR2czZ1O8qmbfqUfnE69Ui4iwFfKUooX7rzy0X/Hy1Iehkj29ZXzezeovTZeaIe9XEucFYfijcqCEeoYDyxCfBXSf+U9DFJKxYllLQu3unlH4xvJAek0cKh1A2DqwF/ljTCSClpkqQv4A5BlSzrl5YdyQWdvYE/pg51HuRsLukHwC304CFrZg8wstP+lqSPNArukl6BG9eWov19517Le0r6ZKPRVB7q8y/4jJ0nOiz2M/iAyV8ajcxJ+PoUIz+o30r3mbMjcKOkTzW+4yyvhSTtD7wjO3xah2UdBp/FBzPAn9OlknYpcmCRtLKkjwPXAe9ukuQv1Af0BJwsaZ4BzaTYHom/nzLkkSc+IekdjQkkbZauPYWR61/PQzIufoi6cPcG4Nw2Btj1JH0dj0zX77WlATCzM4ETskNHSfpicgxsLN8ESdtIOklSu/W8i6iyPgyr3ZxFXXF8TTI2DRQz+zfw2+zQbyTt0phO0mvw+51I+bZQOZJ+I2las/csd/w7CDiaulx6uJmdVzL7fP3rUwtTBUEwpjCzh83sH9nW9aCdmT3HyG/dzySNWC4qybPvxfvWKpfKzpcx26XsIH+HHEddLnk98Insf8eU0G/yvvMzkuaJICOfXf433GDaqaxcJWcDc9L+i4GZapgsImkJST8F9qV8WfP31DaSTYrElUfl+JKkmSoIUS9poqTXS/o1XSxR0nDtfsgA/a4Dz+Czg//QqKMnI+xh1J2hAb6Q2m3OxyRdLunDRXq+pMUlfQN3TgOXw3txohwUH6NuRF8PuKSVYVvSSyR9Dp993WzZlV8BN6T9RYAz5LOv8zyUZOyv0Z2+dIikeRxcJL0Z19faztQfpt2hQ2ZSj5Y0GbfLfFhNJrnInZHeKumkHq5XZX0YVrvJ+5R+LxVSxGnAhWl/InCqpGmNiSTtiH9Hq/z2d4Tcse/nRbYCSctI+iHw7ezw/5nZjSUvEfpSEAwISUtJ+gnw1uzwMd1O2EqOyLnM+RlJX2m0Q0laHjiJ+oDmc5SbxFl03XMYGe3kBPnktRHInZ/Pwu2kALcB3yvI9gvZ/o7Az9UwEC9paVxv24LqxjmqsnuOBZ2p6mfcL92rxrO4HfAPkuZZAk3ugHcK9WXxzmVem+6KwG2SDpXUdJmrJEO+iZEODWNhTOUo6lHhJgF/knSgCiZPJ/vqbpLOZqRTJOA2FkbKEl+RB5YYYatN7fov+CSVMnrCcdSjh20u6RDNO/a2Ct6GyrbtqseTKidFVD0wO7SvpONb2AReIelwNRnHKkmV9WEo7SY9s7+nn4szhKWpk50hv6f9JH21UbeTO/mehvd9w9QTviJp14Jv1kRJ78H7xpoD4KWMdJprxQbA0mn/IjO7v+cCj1KG4XEXBP1mw7R9T9Js3LB9Py5cTcGNSY1rkJ5EEwFhmJjZXZL2BY6lvmbyDZLOxT3ylwC2xr3tn8VnX9S8+Xty6DKzmfIQqjUBegdge0n/Bv4NPIoLQysDr6K6aFQAM4DXAa/GDVY/AD4t6e+48L0WbuCbAPweX3t463TuPPdtZmdI+ltKI1zg+4ikfwGPAC/FhbCJeKjiw4FvdlDeE/EP4qbAtZLOwyMBLYbPyM/X076QkQpszprput+UdCtwJb4mu/BwuptRVywBjjazC+bJZZSR6vGOuODwIlzQOgG4V9LFeCS2Cfh7XBd/DoWDhGZmkt6HP8tl8Y/1GZIuxyO6TcSdImttfH/KeeQfi8+IXh9YADgx1ZFanuvhdR28ju6Ne/u3uvczJX0I+FHKYzPckH0jcBnwEG7AXwEXPIa1fN778XvZJJXzIOpt7jb8fayMG8drs4q7Gsituj4whHZjZk9IOhmfkbQ43l+d3vqsvvAx/J1NxfvjEyRdj8/EeQ5fc702oHEY7gxWq7OF3whJP8/Oq7F0tr9Sam+NvN/M/lGQ7TvxKCUPSPoH3tcKNw69FsgjRByJt8Wy5JEwju7gvCAI5i++ghvEFsX7nsuTfFiTqbfAJx48D/wPDTNNe+D3wNfxPm8H4EpJF1Cf2Q1wXIv+sy1m9oCkM4C34Pr927J/t1sWD3yQ/xO4jL0g7ljzWeAKfMm4dal/F67AB8jLRqGoFDN7SNK3gS+mQ+8B3pxkiDtwuWIa/p6fBz6MO4W043dAbab7IcngdzUjDV4HJ8NdjYPwb3Bt1ute+Izmy3AD7ON43Vodly9rgxCNEyK6oWoZYCb9rQNH4IMxbwBuljQLd/xfGh+gz+WM35nZL+bJwVkf+CHwA0n/xfXS+3HnlJXwdpwbgg8xs9s6KOdQMLN/y0PZH487Lq2NT/C6DTdk3o/rJ8viz6CpcT3L7ylJe+EDk4vgEywukXQRHiFvITxCzKrplAMpF7nge/gg10q4/n9u0heuS3luBNScmPamRDj+IdsdSmFmz8kjDZ+N2w4WwevhwZLOxwdpJuNt7NV4u3+kh+tVWh8YQrsxs5skXQhsDqwlaR0zu7bdeVViZnMl7YMPtiyN69rnJB3/Slzn3AioOUbvR33wvqUjn6TTGGnngZFLIG1UoC9tb2Z3Njk+CXgf8D5Jd+ERPe7G3/NU/B3l9vuDzKxtRIWMndLfZ3D7VRAE3bO9pBc1Ob4I3l43Z6R94wZGTlLohm/jdpOac9XngQ9JOge3K66Ky1P5skyfMrNel7jZG+9DX4zLsb9N9syL8f5kHdzmVrPZPQHslpwt5sHMTpd0BC6fA+yDTxQ5G48yslK6j0Vw595TqC//1OvE9SrsnjMZ5TpT1c+4j7pXjTtxffkA3NH8StxeDj62ly8RfC+wT8GyeEvh7+YTkh5MedyB63LL4Tb9fALtDZSITDRszOx5Se/CgwG8CpcLDgG+mOT6W/G2OAWvly/H5Ttw/bYZB+M6Wa3tHgockMaz5uCy/BbpfyenvLduntUL5bxF0o+p17sDgd3SGOZTeB/ymlT+s3DZueXy3X0YP+gLZnaEPODGlUCDigAAIABJREFUh9KhdwHvkHQpXs+ewmX2V1FfLabbJeCrrg/DajdH42Om4DLqhS3S9osf43pnLRr854APJjvFI/i72hp/fhcD/6W+1HarMZV9cX05p3EC/WmSGh0Mf2xmPy7IdmP8u/9k0mP+i7+j5fG2mssklwFvMrOWARwy5p8xFTOLLbZxseHh12/CjRZltzm40W1Si3xnZumnF6SZkaWZUaKs07L0s9qk3QtXJoru4RHqhuXasVMK8pqapZldopzvwj+AZZ/nxcCCFbzLpfEZFK2u9QfcyPj37NirCvJbHjcmtcrvalxAmJ4dm1mQ34g0uDB2QZv8zwCWKMhvF/wjWuYZP48bMidX8JxL19syz6XN+asDZ3ZQl+4G3tgiv/XwD3/R+XPxASvy423KuEaJPL+KC9azs+NT2+S7DS4slr33fwMr9foO6KC94w5ZP8UHz9qV70lg8SZ5zMrSTOt3fWBI7SZd+3Vl30Un7yE7p1T9SnX2qjb3/jNccM/78ikt8pzVJr+irfCd44pfu/MfwQd2J3TwHlZK79bwmQs9v9vYYout/1tDP7Nvj3kdl+U1o03aN9Fepn572X47P7fNdb/Spv+b3uL5TCv5HN7dJN8rOniOa9FaBjLgfHyAYEa7Z14mTcF5M4ueS5ZmIm5gb1XWh3CjWtl3OQl3cmiV59SCcz8KPNjm3No2lwI9rYu6X6kM0O86gBvVr22T/68pkNVwJ+kyz9hww+AXK3rObetkr/U+O3994B8d3OfNwAYt8tsWHyQoOv8ZfMBianassJ2kPDfEB9JaPfsPp7QvHC9x75XYHTp9B3Rml1kaH6wrU77bC/KYnaVp2qdUWR8YUrtJ135f2XfRyXvIzilVv3DntFta3PdcfFb45OzYI23ynN0iv1Zb03eOO1GVOf8uYNcO38Om2fnHVfV+Y4ttftro3kZyCrBCm7zL9mWT8Anc7Wx2D9NeXind5+L29LNK3OuNwMYlnuWEdB+tbHnX4XLjwdmxAwrym56lmdnm2lXYPauWl0uXv4P6WvUzrlT3akyDf39/WqJ+rVuQ38qUsznWtnNo0y5LPudpWZ7t2lHb51Ki7v4ID4RQ5h7n4JEji/JbkvY2+T+ldLOyY9Na5LlQOqdVnn/EnW9mZsfa9VerlyhrvjUdT+rmHdCZ3P5x3J7Urnxzge2anD8jS9O0v6iyPjCkdpPVv8dTvjcBquo9pPSl6leqsye2ue8L8LGHo7NjO7fIc0ab/Iq2wncO/LnE+bVxr6Zjzy3yrtln5gBLV/F+R+sWEaGCcYOZ/Qxf6mJd3GNzM1yoWx3vYIXPvr4bnwF2NnCijZzZO+ows1+lGev7AdvjMz2exWcs/BH3GL01zVKs0XT2RRfX/q2kU4BdgTfiHqjL4jNBnsCNldcC5wGnmdkNRXl1eN0H5WHLdwfeS3325z34u5uJzxY2eUjXGkWzTu6Rrzf9/nQv6+KzH+4FrsdnWx5tZnOKwkG2Ke9dkrbGvdnfg8+KeRE+KPIPXJkp8sTHzE6Uh6rfDveQXx/3ZK/NeH0Ed6Q5H/i1mV3TaRmHjZndArxe0uZ4dJit8Lo8BVdCH8AVm3/g4Vdn2bxLYuT5XSnplbjH/buoz8i5E3eO+4l1GPnHzG6WtB7ujPH2hjzPA35kaUaVOljZxszOkS/xsDPuab4ZbuxcAhc07sGV0AuA082s2ezRvmJmT+Ke74cBe+KOPlNxg/8zuNH1SnzWwfFm9lhBVmWv13N9GGa7MbOzJF2Fz056p6SPmy/vMVBSnd0Q79t2w505F8Pf16XAz8zsrwCSapGxDJ9dP0i2wpdrmoYP3C6Hz4i4B3/Pp+D1qtO1xvchW06vkpIGQTBuMbM/p+/x/+Jy7aq4weBWXKb+UZKpp1Z83S+k2cb74LOEl8fl0Co5Be/bl8iOlYkGBYCZ3SAPgf8RXAZaG++n78adbY4BTjCPjFJZobvBfJnEvSSdAHwQH+Sdghvgb8WfxS/N7M6y7zLd1xvxd/QO/Pu+NPPO3mt27g8k/QrXWd6AyyPL4oa1x4Db8QkXs3B9qZIIRVXLAP2uA2Z2nXzZp31w2f0l+Hu7F5eBf2pmZ7U4/1BJv8Of8Rb4O5qK1/m5uB54La7j/zrJmmMKM7sCjySzHT6Y9Brq0Zeexh2QbgAuwifZXGjJglmQ39mS1sF1m51xOXkCXifPwvu8f3fS55nZv1I/+gk8MsUaWZ5/BY7oRuYelt2hwzI+CLw91ePdcbl2FbweP4k/g8txY3XPUXeqqA9DbjdH4xERlwX2kfRla79Ua+WY2T8lvQJ3+tsFj+q1IF6nzsftaRfLl5WqUYk9rQPuxSdPbYvrTavh+pLwPvgaPIr9iebLZHXCB7L90JeCoH88jdug/oP3y0ebWU/LIecke9jHUtSVfXCb3ap4lPIH8e/Babj8V0X00dp17wFel+TkXfHIVCvgziv34tEnTgaOshIRKNJ34GOSjsejZmyJ60aP4A5Gx+Fy/ONl7P0d3kvPds+xoDNV/Yz7oXs15P8s/l5OwJ2oN8EnnM/BbeUn4N/qpwrOv0PSMvg3dEvcAfoluPyxAK6P3YLrR8ebL5U4pkh190PypYz3wO91LTwa0gT83d6ERyI7C/hzK3kh2a9fn8YS98InOtT0sqtwx7ffpnG3smV8StJbcL20lueSKc8r8LG8EzvJM+Vb6XhSvzCzwyUdhTs3vhHXzWuReu7H5e2/4XWw7NLGRdfquT4Ms92Y2SOSjsQnlK2B98UDb5epT9lF0g5437Mp/s4ewN/X0cBvzOzZqr9HHfJhvE5ti9er5XFd7n7cP+DPwLFmdl0nmcqXua9Fc/5N0nXHLWphNwmCYAwh6WB8/VxwL99DhlmeQSBpYdyYPwkXkJdIAnq/rzudeqj/X5nZ9H5fM+geSS986MxsuKN3wbhBvgZzbZB5f+tseYKBIumluGEM4Hoze1mr9GMBSZNwI85quIK39jCU3SAIgiAYzQxLBpA0A/hS+nmQmc0YxHWDzkmDVjenn7eY2dShFSYYV0j6HB5RGWAnMztlmOVphaQ34ANoAGeY2ZtapR8LJCfY2/Ali84zs63anBIEQTBqSBNJtkg/Nzezi4ZZnvHIMJ5xyJ1ji7RU2Nbp5zZmNmt4pQnGC5LWwINTTAZONrOd25wyVCTdQX1Z7BXN7O5hlqcKJB2DOy4+D7x8GBONBsmE9kmCIBjtyF2pd8kOXTqssgyYt8MLke3+NQgnqCAIgsSx+Oxc8PW0R3OUzXdl++Pl+7A77gQFPsAaTlBBEARBMC/jUQYIgmBscDg+WxngwGEWpATjsa/8GO4EBfD5YRYkCIKgEySthkfnAI/WdMUQizMuiWccBMGwMLObgV+mnzumiMOjEkm1qLzgy5+PByeoNajrPr8e705QEI5QQTBe2A8PgQgeyvVvQyzLQJC0FPCV7NCxwypLEATzHync9KfTz9WAvYdYnEIkrY4vA1VjzPeVyemsFgHxcjpY/ikIgiAI5hfGowwQBMHYwcweB2akn5unqEujDkkb4cu41DhuWGWpCklL4nZCgD+Y2bnDLE8QBEFZ0mTv7wIT06FT0lJQQUXEMw6CYBQwA1+CT8AXhluU5kiaDByWHRov9pTP4f3/HOpRvMc14QgVBKMYSVtI+qmk9Qv+v1hadiDvkL8z1iNjSDpC0nRJixX8f1PgPHwdWXDnrxgID4JgoJjZn4A/pZ8z0nKdA0PS8ZLeIWnBgv+/ETgfWCodugpfO3qs8z5g7bT/0eSUFgRBEATzDfOxDBAEwdjix9SjTBySBl8HhqQzJG0naWKT/02QtBvwV3xpDoDTzOzqQZaxT3wGWAZ4GjhgyGUJgiAAQNJHJX1e0soF/18NOBGoLZM0l5FjHkEb4hkHQTAWSJGVvpx+7iZpg0FeX9KXJe0naZmC/68DnAlskg7NAY4YVPn6Rbqv6enn18zstiEWZ2CM5mVcgiCABYAPAB+QdAse+eJe3GNzVXwd50Wz9OcyPoTXlwMfAn4k6TLgP/jHZkngVdQHwAGeA/Yxs0cHXsogCOZ7zOwtQ7z8pngo08ck/QuYDTyFG703ob50HMATwJ7jwWnIzH4C/GTY5QiCIAiCITJfygBBEIwtzOx5YKADGw1sl7YHJf0TuB14FlgO2BxYPkt7D/DBgZewD5jZ/wH/N+xyBEEQNPAiPPrEQZKuBq4BHgYWwVe6eDUjxyu/bmYXDbyUY5t4xkEQjAnM7NvAt4d0+dXwSFSHSroCuB54FFgcWBdYD49WVWN/M5s96EJWjZldy3zoFzTf3XAQjGFWT1sRxwHvS4am8cJCuHFq84L/3wvsZWYxuzkIgvmZxYGt09aMm4B3mdnlgytSEARBEAQDIGSAIAiC9iwNtFqa73LgHWZ2x4DKEwRBMD8zAXhl2prxNHCQmX19cEUad8QzDoIgaM8k3EH01QX/fxTYz8x+NbgiBVUTjlBBMLo5F9gW2B6f2bsi7tm/OPAIPpvtXOA3ZnbpsArZB3bFQ7ROw2crLIvPbp4LPICHVv8zcKSZzRlSGYMgCIbNNsDbgS2BNfHvwzLAM8B9wD+BU4FjzOzZYRUyCIIgCILKCRkgCIKgPesCOwGvwScWvgiYgkccvxe4CDgZOMnMbFiFDIIgmE/4Fh6h6A14tI3l8H55AeBB4EbgbODnZnb7sAo5xolnHARB0J79gVn42PvL8fHnZfEoUA8CV+NL4/3czB4cUhmDilDoeUEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQjHUmDLsAQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEvRKOUEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQjHnCESoIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgjFPOEIFQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQTDmCUeoIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAjGPOEIFQRBEARBEARBEARBEARBEARBEARBEARBEATBmCccoYIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgCIIgGPNMGnYBgv4jaUHglennfcDzQyxOEARBEARB0B0TgWXT/lVm9vQwCxMENULfCIIgCIIgGBeEvhGMSkLfCIIgCIIgGBcMVN8IR6j5g1cClw67EEEQBEEQBEFlbAz8Y9iFCIJE6BtBEARBEATji9A3gtFE6BtBEARBEATji77rG7E0XhAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEQRAEY56ICDV/cF9t55JLLmHFFVccZlmCIAiCIAiCLrjrrrvYZJP/Z+8+o22ryvuPf38X6SBFghSp9iiCEVAQIomIBSSCUYoJAVvAWFAc0SQmUdMUxRIFGwpY0D8EuYLECIqBKKiIEEGIgF6wgCBVilLuff4v1jrezeH0e/bad+/z/Yyxxlx77znXfO7gzXlYz5xz57GPv5qqr9Qx8w1JkqQhZ76hlZj5hiRJ0pDrOt+wEGph+N2Z2ZtuuimPetSjBhmLJEmSVtzS6btInTHfkCRJGi3mG1qZmG9IkiSNlr7nGx6NJ0mSJEmSJEmSJEmSJGnoWQglSZIkSZIkSZIkSZIkaehZCCVJkiRJkiRJkiRJkiRp6FkIJUmSJEmSJEmSJEmSJGnoWQglSZIkSZIkSZIkSZIkaehZCCVJkiRJkiRJkiRJkiRp6C3YQqgkGyfZJ8k7k3wlyc1Jqr1O7NOcByb5apIbkvw2ybVJPpPkGf2YT5IkSdLgmHNIkiRJ6hfzDUmSJGliDxt0AAN0Y1cTJVkDOBXYZ9xPW7XXwUneXlX/1FVMkiRJkvrOnEOSJElSv5hvSJIkSRNYsDtCjfMz4Ow+Pv+TLE8QvgG8CNgZeAXwY5r/Du9M8so+xiBJkiSNpHZV8mqDjmMa5hySJEmS+sV8Q5IkSWot5B2h3glcBFxUVTcm2RpYMt+TJHkWcHD78Uxgv6pa2n6+KMkZwMXAlsDRSf6jqm6f7zgkSZKkEXYycFuSk4FPVdUlgw6oZc4hSZIkqV/MNyRJkqQJLNgdoarqH6vqy1XV7+1j/7ptlwKv6UkQxuK4GXhL+3EDmhUUkiRJkmZnA+A1wPeSXJLktUk2GGRA5hySJEmS+sV8Q5IkSZrYgi2E6kKSdYBntx/PqaqfT9L1i8Cv2/v9+x6YJEmSNFr2BU4HHgACbA98ELg+yeeT7DXI4PrJnEOSJElSv5hvSJIkaRhZCNVfOwOrt/fnTdapqu4Dvj02Jsmq/Q5MkiRJGhXtKugXA5sBRwGX0RRErQ68FPhKkuuSvCPJNgMMtR/MOSRJkiT1i/mGJEmSho6FUP31xJ77/5um79jvDwMe259wJEmSpNFVVbdU1furantgR+A44HaaoqgtgLcBVyf5epKDk6w+xeOGhTmHJEmSpH4x35AkSdLQedigAxhxW/TcT7Zl7JifjRt3xUwnSfKoabpsMtNnSZIkSaOgqr4PfD/Jm4AXAYcBz6FZDPJHwB7Ah5N8Hjihqr43qFhXUN9zDvMNSZIkacEy35AkSdLQsRCqv9btub9rmr5399yvM8t5fjZ9F0mSJGnhaY9oOAU4JcnmwF8ArwMeCawPHA4cnuRy4CPAiVX120HFOwdd5BzmG5IkSdLCZL4hSZKkoWMhVH+t0XN/3zR97+25X7MPsUjSSNv6rWf97v7ad+09wEgkSSujJGsBewJ7ARsDRXNkXtou2wHHAv+Y5LVVddpAAp09cw5J6oD5hiRpgTLfkKSOjOUc5huStOIshOqv3pXkq03Td/We+9/Mcp4tpvl9E+CiWT5TkiRJGnpJdqM5Fu8lwNpjXwO3AZ8DPkNTBPVyYFeanaJOSbJ3Vf1X9xHPWhc5h/mGJEmStDCZb0iSJGnoWAjVX3f23E+3FezaPffTbTH7IFU15dncSab6WZIkSRopPUfgHQo8euxrml2gzgOOB06rqrEVyxcBn0qyK01h1DbA3wHDUAjV95zDfEOSJElasMw3JEmSNHQWDTqAEdf7B/yjpunbu+rBM7ElSZKkWUiyWpIDkvwXcC3wT8BjaAqgbgTeDTyuqv64qk7uKYL6naq6ADiq/bhdN5GvMHMOSZIkSf1iviFJkqSh445Q/XVFz/0Tpuk79vsDwDX9CUeSJEkaPUmOBQ4C1hv7ClgKfBX4BPDlqlo6w8f9sG3Xndcg+8ecQ5IkSVK/mG9IkiRp6LgjVH9dBNzX3j9rsk5JVgOeMTamqu6brK8kSZKkhzgCWJ+mAOo64B+Arapqn6r60iyKoADuBX7aPmcYmHNIkiRJ6hfzDUmSJA0dC6H6qKruBL7eftwzyWRbx+4PPLy9P73vgUmSJEmj5X7gVOC5VbVtVf1zVV0/lwdV1U+rauuq2nZ+Q+wPcw5JkiRJ/WK+IUmSpGFkIdQKSHJokmqvt0/S7b1t+zDg2CSrjHvGRsC724+3A8f3JVhJkiRpdG1WVQdU1TmDDmS+mXNIkiRJ6hfzDUmSJI2ihw06gEFJshvwmJ6vNuq5f0ySQ3v7V9WJc5mnqs5N8gXgQGBf4JwkHwCuB7YD/g7Ysu3+1qq6bS7zSJIkSQtVVd0y6BgmYs4hSZIkjYYkBwJfXJmOfDPfkCRJkia2YAuhgFcCfzHJb89sr14nrsBcL6fZFvYFwB+1V69lwD9V1cdWYA5JkiRpwUoy9j/eb6yqe6fpuwawMTRH4fUxLHMOSZIkaTScDNyW5GTgU1V1yaADwnxDkiRJmpBH43Wgqn5TVXsDLwPOAW4C7gN+RpNA7VZVbx9chJIkSdLwaldCXwtcBqw1gyFrApcDP0mycx9D64w5hyRJktR3GwCvAb6X5JIkr02ywaCD6oL5hiRJkobJgt0RqqoOBQ5dwWecyCxWUVTVyTRJgSRJkqT5c0DbLp7JMQxVdVuS02hWTx8IfLcfQZlzSJIkSSNjX+AwYB9gVWB74IPAe5IsBk6oqrO7DMh8Q5IkSZqYO0JJkiRJGna7AEWzMnmmxl5S7Db/4UiSJEkaJVX15ap6MbAZcBTNbrQBVgdeCnwlyXVJ3pFkmwGGKkmSJC14FkJJkiRJGnZbtO2PZjHmmrbdfJ5jkSRJkjSiquqWqnp/VW0P7AgcB9xOUxS1BfA24OokX09ycJLVBxiuJEmStCBZCCVJkiRp2K3XtktnMWas7yPmORZJkiRJC0BVfb+qXgtsSnPk9ldpdqpdBPwR8BnghiTHJtlxcJFKkiRJC4uFUJIkSZKG3c1tu+0sxoz1vW2eY5EkSZK0gFTVfVV1SlU9H9iKZleoG2l2iVofOBz4TpL/TXJ4kjUGGK4kSZI08iyEkiRJkjTsLm3bA2Yx5sC2vXyeY5EkSZK0ACVZC9gT2AvYmGZ3KGgKogJsBxwLLEny4oEEKUmSJC0AFkJJkiRJGnZfonmxsH+Sl0zXOclLgf1pXkws7nNskiRJkkZYkt2SfBL4JfApYHea/OR24MPA04FXAhe03z8SOCXJ8wYTsSRJkjTaLISSJEmSNOxOBK6lealwcpJjkmwxvlOSLZK8D/gcTRHUz4DjO4xTkiRJ0ghIsnmSv01yFXAecCiwTvvzecCfAZtV1eur6qKq+lRV7QbsBiyhyV3+bgChS5IkSSPvYYMOQJIkSZJWRFXdn2R/4Hyalw9HAkcm+SlwA03R02bAlu2QAHcB+1XVvQMIWZIkSdKQSbIasB9wGPBsmoXmaX/+JXAScHxV/XiyZ1TVBUmOAr5Ic1SeJEmSpHlmIZQkSZKkoVdVlyZ5BvBZYIf26614cPHTmIuBP6+q/+swREmSJElDKsmxwEHAemNfAUuBrwKfAL5cVUtn+Lgftu268xqkJEmSJMBCKEmSJEkjoqquAP4gyV7A3sBTgY3an28Gvg+cWVVfH1CIkiRJkobTET331wKfAj5VVdfP4Vn3Aj8Fls1DXJIkSZLGsRBKkiRJ0kipqrOBswcdhyRJkqSRcT+wmObou3NW5EFV9VNg6/kISpIkSdJDWQglSZIkSZIkSZI0uc2q6pZBByFJkiRpeosGHYAkSZIkSZIkSdLKyiIoSZIkaXi4I5QkSZKkkZJkXWAbYF1glen6V9X5fQ9KkiRJ0lBLsmV7e2NV3TtN3zWAjeF3R+FJkiRJ6oiFUJIkSZJGQpJXAa8BnjKLYYV5kSRJkqQpJNkNOB+4E9gamLIQClgTuBxYK8muVfXd/kYoSZIkaYxH40mSJEkaaklWSbIY+ChNEVRmeUmSJEnSVA5o28VVddt0nds+p9G8gzmwn4FJkiRJejBXPkuSJEkadocD+7b3NwInABcDtwLLBhWUJEmSpJGxC81usufMYszZwF8Au/UlIkmSJEkTshBKkiRJ0rA7pG2vAHafyQptSZIkSZqFLdr2R7MYc03bbj7PsUiSJEmagkfjSZIkSRp2T6RZnf1PFkFJkiRJ6oP12nbpLMaM9X3EPMciSZIkaQoWQkmSJEkaFbNZnS1JkiRJM3Vz2247izFjfV2sIUmSJHWok0KoJG79KkmSJKlfrm7bDQcahSRJkqRRdWnbHjCLMQe27eXzHIskSZKkKXS1I9S1Sc5K8uIkq3Y0pyRJkqSF4QtAgH0GHYgkSZKkkfQlmpxj/yQvma5zkpcC+9Mc4b24z7FJkiRJ6tFVIdQqwPOAU4BfJHlfku06mluSJEnSaPt34AfAEUl2H3QwkiRJkkbOicC1NMVQJyc5JskW4zsl2SLJ+4DP0RRB/Qw4vsM4JUmSpAWvq0Ko9wO/okkSNgLeAFya5KIkf5nk4R3FIUmSJGnEVNW9wF7AxcA5SY5OskOSNQYcmiRJkqQRUFX30+zwdBfNwu8jaU7CWJLkgiTfSrKEpljqDW2fu4H92nxFkiRJUkc6KYSqqqOAzYEXAWcAS2mKop4GHAfckOQzSf64i3gkSZIkjY4kS4EbgF2A1YCjaIqi7k6ydJrrgUHGLkmSJGk4VNWlwDOAS2nebwTYCnh6+/1WPd9fDOxcVZcMJlpJkiRp4XpYVxNV1VKaIqgzkvwecAhwKPAkYE3gYODgJNcCJwAnVdXPuopPkiRJ0tDKNJ8lSZIkaYVV1RXAHyTZC9gbeCrNKRgANwPfB86sqq8PKERJkiRpweusEKpXVf0KOAY4JslOwMuBA4D1gW2AdwBvT/I14JPA4nbrWUmSJEka7x2DDkCSJEnSwlFVZwNnDzoOSZIkSQ81kEKoXlV1EXBRkiOBF9PsEvXHNMf2Pae9bkvyWeDj7YoLSZIkSQKgqiyEkiRJkiRJkiRJLBp0AGOq6l7gAuBCmi1ki+XnaW8IvA64LMniJNsOLFBJkiRJkiRJkiRJkiRJK52B7wiVZE3gT4HDgD9kefETwGXAZ4HtaHaLWhN4IfCHSZ5ZVVd2H7EkSZIkSZIkSVqokqwLbAOsC6wyXf+qOr/vQUmSJEkCBlgIlWRXmuKnlwLrjH0N3AV8ATi+qr7b0/+1wGuBvwPWA/6ZpjhKkiRJkn6n3UF2F2ATYC3gI1V182CjkiRJkjTskrwKeA3wlFkMK1aCRemSJEnSQtHpH99JNgMOAQ4FHjv2ddt+Bzge+EJV3T1+bFX9GvjXJLcBxwK79j1gSZIkSUMjyVOBDwC7jfvpNJrjt8f6/RXwj8AdwO9X1f2dBSlJkiRp6CRZhSaveOHYVwMMR5IkSdIUOimESvISmt2fngMsYnmScCvwGZrdn344w8ed17Ybz2uQkiRJkoZWkr2B/wBW48EvJWqC7icB7wIeAewDnN73ACVJkiQNs8OBfdv7G4ETgItp3nEsG1RQkiRJkh6qqx2h/h/NC4i07bk0uz99sarum+Wz7p3n2CRJkiQNsSSbAJ8HVgd+CLwZ+CZw50T9q+quJIuBg4HnYyGUJEmSpKkd0rZXALtX1W2DDEaSJEnS5Lo8Gu8GmpXXx1fVkhV4zrXANvMSkSQp+o1pAAAgAElEQVRJkqRR8EZgHeA6mpcStwMkU55W8d/Ay4Cn9Ts4SZIkSUPviTSLvP/JIihJkiRp5dZVIdSfAGdV1QpvEVtVS2lecEiSJEkSwHNpXkocM1YENQM/atut+xKRJEmSpFH0o+m7SJIkSRqkRV1MUlVnzkcRlCRJkiRNYGzH2O/OYszYsXnrzHMskiRJkkbP1W274UCjkCRJkjStTgqhJEmSJKmPVm3b+2cxZv22vXueY5EkSZI0er4ABNhn0IFIkiRJmlonhVBJHpfkJ0muSrLZDPpvnuSaJD9OslUXMUqSJEkaWr9s222m7PVgu7Ttz+c5FkmSJEmj59+BHwBHJNl90MFIkiRJmlxXO0IdAGwNXFNV10/Xuap+AVzVjjmwr5FJkiRJGnbfatv9ZtI5yVrA4UAB5/crKEmSJEmjoaruBfYCLgbOSXJ0kh2SrDHg0CRJkiSN01Uh1HNpXjKcOYsxX6LZavYFfYlIkiRJ0qg4iSZ3OCjJXlN1TLIOcAqwZfvVJ/scmyRJkqQhl2QpcAPNzrKrAUfRFEXdnWTpNNcDg4xdkiRJWmi6KoQae8nwg1mMuXzcWEmSJEl6iKr6GrCYJr85I8l7kuzc02XDJE9P8vfAj4Dn0yzU+HRVXdJ9xJIkSZKGTHqu8Z9nckmSJEnqyMM6mmfjtr1rFmPG+m4yz7FIkiRJGj1/BnwZ2AN4U3tV+9t5Pf3GXkJ8neZ4PEmSJEmazjsGHYAkSZKkmemqEOoOYCOaoqb/neGYsQKoe/oSkSRJkqSRUVX3JNkTeCNNEdSmk3S9FXgvcHRVLesqPkmSJEnDq6oshJIkSZKGRFeFUFfTFEI9D/jqDMc8v21/3JeIJEmSJI2UtrDpmCQfBHYGdqTZnXYV4BbgEuCbVXXv4KKUJEmSJEmSJEn90lUh1FeBXYFXJ/l4VV05VeckTwJeRXOUxX91EJ8kSZKkEVFVDwAXtJckSZIkSZIkSVogFnU0z0eAu4E1gHOTvHCyjkn2Bb4GrAn8Bji2kwglSZIkSZIkSZKmkWTbJC9LclSSv0+y0aBjkiRJktToZEeoqro5yeHAZ2iOplicZAnwP8ANNDs/bQbsDmwDpP3uiKq6sYsYJUmSJEmSJEmSJpPkqcAHgN3G/XQacHNPv78C/hG4A/j9qrq/syAlSZKkBa6ro/Goqs8lWQU4DlgL2Jam6KlX2vZumiKoz3YVnyRJkqThlOSQFRlfVZ+er1gkSZIkjaYkewP/AazG8ncZ0CzqHu8k4F3AI4B9gNP7HqAkSZIkoMNCKGheMCQ5B3g98ALgySxPGJYBlwFnAh92JyhJkiRJM3QiE798mIkCLISSJEmSNKkkmwCfB1YHfgi8GfgmcOdE/avqriSLgYOB52MhlCRJktSZTguhAKrqBuBvgL9J8jBgw/anW6vqga7jkSRJkjQSMn0XSZIkSZqTNwLrANcBu1fV7QDJlGnIfwMvA57W7+AkSZIkLdd5IVSvtvDppkHGIEmSJGnojT9yeyJrA4+neRGxP3AB8GqaY7klSZIkaSrPpdlN9pixIqgZ+FHbbt2XiCRJkiRNaKCFUJIkSZK0oqrquhl2vQI4PclBwGeADwF79i0wSZIkSaNibPHFd2cxZuzYvHXmORZJkiRJU1g06AAkSZIkqUtV9Xng08AewF8ONhpJkiRJQ2DVtr1/FmPWb1t3oZUkSZI61GkhVJInJnl/ku8luTXJ/UmWTnM90GWMkiRJkhaEU4AAhw44DkmSJEkrv1+27UyO5R6zS9v+fJ5jkSRJkjSFzgqhkrwJuBR4PfAHNKshVqF5+TDdJUmSJEnzaexFxuMHGoUkSZKkYfCttt1vJp2TrAUcDhRwfr+CkiRJkvRQD+tikiTPA97bfizg28DFwK3Asi5imEqSLWkKtPYGtgTuBa6hWSV+XFXdMw9zbAUcATwHeDSwNs0Z4f8HfAX4aFX9akXnkSRJkjQjW7btqlP2mgfmG5IkSdLQOwl4GXBQks9U1dmTdUyyDvAFmr/9C/hkv4Mz55AkSZKW66QQCjiybW8D9q2qb03VuUtJ9gY+B6zX8/VawE7t9cokL6iqn6zAHAcDn2if22sDmu1xdwHekOSlVXXuXOeRJEmSNL0kqwJ/3X68ps9zmW9IkiRJQ66qvpZkMfAi4IwkHwJO7emyYZKnA3vR7AS1CU0R1Ker6pJ+xmbOIUmSJD1YV4VQO9L80f/OlawIanuaFRFrAXcB/wZ8A1gTOBB4Fc1RGWcl2amq7prDHLsAn6Y5BnAZzcqRLwHX06zM+AvghcAjaBKoJ1fVtSv2L5MkSZIWjnb183QW0fxP+h2B1wJPpslRvtDHuMw3JEmSpNHxZ8CXgT2AN7VXtb+d19Mvbft1mqKovjHnkCRJkh6qq0KosVUC3+xovpn6AE1sDwB7VdWFPb+dm+Rq4GjgCTRJzTvnMMff0iQIAK+rquN6frsIOC3JMe3z127b189hHkmSJGmhWjKHMQEuBN4/z7H0Mt+QJEmSRkRV3ZNkT+CNNH9XbzpJ11uB9wJHV9WyPodlziFJkiSNs6ijeX7Rtqt1NN+0kuxEs3ID4JPjEoQxxwBXtvdHtkdozNYz2/aWcQlCr97kY9c5zCFJkiQtZJnldRvwLmDPqrq3LwGZb0iSJEkjp6qWVdUxwFbAbsCRwL8C76Y5fvs5wOZV9a5+F0GZc0iSJEkT62pHqDOBN9D8wTzRH+OD8KKe+xMm6lBVy5J8mmY72Q1okopzZjnPWPHXpKvUq+qOJDcDGwGrz/L5kiRJ0kJ32Az6LAPupPm7/PKqWtrfkMw3JEmSpFFVVQ8AF7TXoJhzSJIkSRPoqhDqvcAhwFFJPltVv+xo3qns3rZ3AxdP0a/3bO/dmH2ScBXwVGCbyTokeThNgjDWX5IkSdIMVdVJg45hAuYbkiRJkvrJnEOSJEmaQCdH41XV9cCf0JwjfUGSF3Qx7zSe2LbXtKs3JvN/E4yZjY+17SOSHD5Jn7+foL8kSZKk4WW+IUmSJKmfzDkkSZKkCXSyI1SSc9vbW4HHAWcmuR24GrhnmuFVVc+e53jWYPnqhJ9PM/ltSe4G1ga2mMN0x9OszHgZcGySpwFnADcAWwJ/BuzX9n13VZ092wmSPGqaLpvM9pmSJEmS5sZ8Q5IkSRotSQ5ZkfFV9en5igVGK+cw35AkSdJ86+povD2A6vkcmvOod55iTLX9aoo+c7Vuz/1dM+g/liSsM9uJqmop8GdJzgDeCryyvXp9A3jXXF5KtH42x3GSJEmS5p/5hiRJkjRaTmTu7yoKmNdCKEYr5zDfkCRJ0rzqqhDqfPpT0DRXa/Tc3zeD/ve27ZpzmSzJE4CDge0m6bILcEiSy6rqhrnMIUmSJC1UPTvQzqcV2ZnWfEOSJEkaPRl0AD3MOSRJkqRJdFIIVVV7dDHPLPy25361GfRfvW1/M9uJkuxOs03s+sB1wNuAc2iOCXwksC/wzzTbyj4ryV5VdeUsp5luO9tNgItm+UxJkiRpWOzB5DvKjr2smOz78eZjZ1rzDUmSJGm0bDODPmsDj6f523t/4ALg1TS7Mc23Uco5zDckSZI0r7raEWplc2fP/Uy2gl27bWeyxezvJFkd+DxNgvBL4BlV9cueLj8HjktyHvA94FE0W+TuNJt5qmrKM8CTlWmhiiRJkjTvxnag3Qx4bPtdAdcCN9IUNm0MbM3yIqergH6tVDbfkCRJkkZIVV03w65XAKcnOQj4DPAhYM8+hDQyOYf5hiRJkubbokEHMAhV9Vvg5vbjo6bqm2QDlicJsz2r+nnA5u39h8YlCL3x/BD4bPtxxyTbz3IeSZIkacFqd6B9N/B7wK+BtwCPrKpHV9WuVbVLVT2aZrXyW9o+GwPvrqo/muxagXjMNyRJkqQFrKo+T1MQtAfwl314vjmHJEmSNImBFUIl2TzJ05LsnmRO51KvoLGtWR+TZKqdsZ4wwZiZemLP/fen6XvxJHNKkiRJmkKSxwGn0uz29Myqek9V3Ty+X1XdXFXvAZ7Z9j2lHdsP5huSJEnSwnYKTd5xaJ+eb84hSZIkTaDTQqgk6yZ5R5KfAT8Fvgv8N+PO105yYJJTknyij+F8s23XBp42Rb9n9dx/a5ZzPNBzP90xhKtOMk6SJEnS1I6i+bv+6HYl8pSq6grgaJojJN7cp5jMNyRJkqSFbWz3pMf36fnmHJIkSdIEOiuESvIYmhUDb6PZSjXtNZELgf2AlyfZrU8hLe65P2yiDkkWAYe0H28HvjHLOZb03O8+Td/eZGTJpL0kSZIkjfccoIBzZzFm7G/7Pec/HMB8Q5IkSVrotmzbVafsNXfmHJIkSdIEOimESrI6cBbwaOAemtXX+0zWv6quY/kf5Pv2I6aq+i7wP+3HVyTZZYJuR7F869cPVtX9vT8m2SNJtdeJE4z/Os2/F+CIJNtNFEuS59MUfgH8Arh05v8SSZIkacHbdA5jqm03mc9Afvdw8w1JkiRpwUqyKvDX7cdr+jGHOYckSZI0sem2Mp0vhwOPBe4Gdq+qSwGSyTaEAuArNKuzJ/rjfb68gWYr2DWBs5P8K00B1prAgcCr235XAcfM9uFVdXuSdwHvBNYFLkjyIeAc4DbgkcCfAK9ieVHaW6tq2Zz/RZIkSdLCczuwMc0K5O/McMwebXtHPwJqmW9IkiRJIyDJltP3YhGwAbAj8FrgyTQLML7Qx9DMOSRJkqRxuiqE2p/mD/4PjhVBzcAP2vax/QkJquqSJAcAnwUeDvzrBN2uAvauqjvnOM0/AxvSJCTrAH/TXuPdD/xtVX12jvNIkiRJC9U3gRcDb02yuKqumqpzkscBb6HJUb7Zr6DMNyRJkqSRMZej3gJcCLx/nmP5HXMOSZIk6aE6ORoP+P22PXsWY25p2/XnOZYHqaozgafQJCNX0WzzejvwPZqXI0+tqjlvXVuNNwI7AR8FLgfuBJbSrD6/GHgf8OSqeu8K/FMkSZKkhep9wDJgPeDbSY5MsuH4Tkk2SPIG4AKaPGMZc1gVPRvmG5IkSdJIyCyv24B3AXtW1b39DMycQ5IkSXqwrnaEWrdtZ3PsxBpte/+UveZBVV0HvKm9ZjPuv2mSmpn0vZgmIZAkSZI0j6rqwiRvAd5DUwx1DPDeJEuAm2h2fnoksA3LX0wAvKWqvt1BfOYbkiRJ0nA7bAZ9ltEUCC0BLq+qpf0NaTlzDkmSJGm5rgqhbgE2oXn5MFPbte2N8x+OJEmSpFFSVcckuRb4d2BTmv+Z/2hg27ZL7//cvwF4XVV9sdMgJUmSJA2lqjpp0DFIkiRJmpmuCqEuBZ4HPJuZH4/3cpqV29/pV1CSJEmSRkdVnZbkDOBFwJ40iys2oCmCuhW4DPgasLiq+r7zrCRJkiRJkiRJ6lZXhVBfBJ4P/GWS49ptWieV5B+Bp9MUQv2/DuKTJEmSNALaAqdT20uSJEmSJEmSJC0gizqa50TgSmBd4LwkeyfpPZqikixKsnuSM4F/oCmCuqiqzugoRkmSJEmSJEmSJEmSJElDqpMdoapqaZJ9gW8BWwJnAPf0dDkTeCSwVvs5wPXAS7qIT5IkSdJoaRdebEiTY1xfVUsHHJIkSZKkIZXk3D48tqrq2X14riRJkrSgdXU0HlX14yQ7AJ8A9gbWbn8KsO247mcDh1XVDV3FJ0mSJGm4JVkFOAQ4DNgJWI1mp9mnAFf09NsH+EPgjqr6lwGEKkmSJGm47EGTW6Rte42dfjHZ9+NN9hxJkiRJ86CzQiiAqvol8MIkTwL+BNgR2BhYBbgFuAT4UlV9r8u4JEmSJA23JBsDi4GnM/kLhzFLaHaprSRnVdWl/Y5PkiRJ0lA7n6ZwaTPgse13BVwL3EiTg2wMbM3yIqerABd7S5IkSR3rtBBqTFX9EPjhIOaWJEmSNFqSLKIpbNoZWAacSvOi4sMT9a+qHya5EHgGsB9gIZQkSZKkSVXVHkmeB5wM/Br4F+CEqrq5t1+SjWh2qP1bmsKoI6vqv7qOV5IkSVrIFg06AEmSJElaQYfQFEHdD+xdVQdW1XHTjDmTZqX2bv0OTpIkSdJwS/I4mgUXAZ5ZVe8ZXwQFUFU3V9V7gGe2fU9px0qSJEnqiIVQkiRJkobdQTRHT3ysqr46wzGXtO3j+xOSJEmSpBFyFLA2cHR74sWUquoK4GhgHeDNfY5NkiRJUo9OjsZLsuWKjK+qn85XLJIkSZJGzg5te8YsxtzUto+Y51gkSZIkjZ7n0Cy+OHcWY77RtnvOfziSJEmSJtNJIRSwZAXGFt3FKUmSJGn4rN+2N03Z68FWbdtl8xyLJEmSpNGz6RzGVNtuMp+BSJIkSZpaV0fjZQUvSZIkSZrMbW07m92dxo7E+9U8xyJJkiRp9Nzets+axZg92vaO+Q1FkiRJ0lS62mnpsBn0WZvmZcSLgc2AC4BP9DMoSZIkSSPhCpoXErux/PiJ6RxMs0L74n4FJUmSJGlkfJPm3cVbkyyuqqum6pzkccBbaHKOb3YQnyRJkqRWJ4VQVXXSTPsmeTPw78CrgQuq6q/7FpgkSZKkUXAGzWrr1yQ5tqpunapzksOA59K8lDi9/+FJkiRJGnLvA/YD1gO+neSdwKfH5x5JNgAOAf6e5gjvpcAxHccqSZIkLWhdHY03Y1V1f1UdAZwPHJXkuYOOSZIkSdJK7WPA9cDGwDlJnjRRpyRbJPkQzc6zBVwNnNxZlJIkSZKGUlVdSLPDU2iKoY4BbkpydZJvJflmkqtpjt5+H7BhO/QtVfXtgQQtSZIkLVArXSFUj4/QJBWvG3QgkiRJklZeVfUbmtXZ9wA7AD9IckVPl48muRK4FngNTR50F/CnVbWs43AlSZIkDaGqOgZ4CfBLmncXi4BHA88AdmnvF7W/3UCTb7xvMNFKkiRJC1cnR+PN0dVtu+NAo5AkSZK00quqi5LsCnwW2A54Qs/Pz6R5GTHmSuCAqrq8wxAlSZIkDbmqOi3JGcCLgD1pco8NaPKNW4HLgK8Bi6vq/oEFKkmSJC1gK3Mh1HrjWkmSJEmaVFVdBmyfZG/gT2gWVWwMrALcAlwCnAGc5k5QkiRJkuaiLXA6tb0kSZIkrWRW5kKov2jbGwYahSRJkqShUlVnAWcNOg5JkiRJkiRJktStla4QKsljgaNoCqEK+M/BRiRJkiRpZZbkkPb2R1X1nYEGI0mSJGlBSBJgQ2At4PqqWjrgkCRJkiTRUSFUkp/MoNsiYH1g3Z7vbgL+pS9BSZIkSRoVJ9IsojgIsBBKkiRJUl8kWQU4BDgM2AlYjSYXeQpwRU+/fYA/BO6oKt9xSJIkSR3qakeorecw5tvAYVXl0XiSJEmSpnIH8HDg6kEHIkmSJGk0JdkYWAw8Hcg03ZcAZwCV5KyqurTf8UmSJElqdFUIddIM+iwD7qRJEM4zMZAkSZI0Q0uA7YENBh2IJEmSpNGTZBFNYdPONO8yTgXOBz48Uf+q+mGSC4FnAPsBvu+QJEmSOtJJIVRVHdbFPJIkSZIWpNOBHYAXAucOOBZJkiRJo+cQmiKo+4F9q+qrAEkmLIRqnQnsAuzW//AkSZIkjVk06AAkSZIkaQV9ELgOOCLJHw86GEmSJEkj5yCggI+NFUHNwCVt+/j+hCRJkiRpIhZCSZIkSRpqVfVr4DnA/wFfTfLxJHsk2TBJBhyeJEmSpOG3Q9ueMYsxN7XtI+Y5FkmSJElT6ORoPEmSJEnqlyRLez8Cr2ivsd+nGl5VZV4kSZIkaSrrt+1NU/Z6sFXbdtk8xyJJkiRpCp38D/8k/9CP51bVO/vxXEmSJElDZXylk7tASZIkSZpPtwG/x+x2dxo7Eu9X8x+OJEmSpMl0tfL57TTnZ883C6EkSZIkvWPQAUiSJEkaaVcAzwJ2A74xwzEH07wXubhfQUmSJEl6qC6PgBhblV1Mv0J7Jn0kSZIkiaqyEEqSJElSP50B7AG8JsmxVXXrVJ2THAY8l+Zdx+n9D0+SJEnSmEVdTFJVi4BtgYtoCpy+ArwE2ApYo722Av60/S1t322qatFkVxexS5IkSVo5JHlfe2086FgkSZIkLSgfA64HNgbOSfKkiTol2SLJh4BP0BRBXQ2c3FmUkiRJkrrZESrJesDZwDbAIVX12Qm6/ay9vpjkZcBJwNeS7FhVd3QRpyRJkqSV2pE0LxOOB24a+zLJEmAZ8NyqumZAsUmSJEkaUVX1myT7AecCOwA/SPKjni4fTfJ7wOPazwHuBP60qpZ1G60kSZK0sHW1q9IbgccAn5ikCOpBqupzNCsmHg0c1efYJEmSJA23rYCtgdUGHIckSZKkEVVVFwG7ApfTFDo9oefnZwKPb78PcCXwzKq6vOs4JUmSpIWuq0KoF9Os3D51FmNOadv95z8cSZIkSUPonrbdaKBRSJIkSVqQquqyqtoeeCHNTrWXAjfQ7Fh7BfA54ABgO4ugJEmSpMHo5Gg8mtXZALM54m6s71bzG4okSZKkIXUNsB1wSJL/qaoadECSJEmSFp6qOgs4a9BxSJIkSXqorgqh7m/b7YDvz3DMduPGSpIkSVrYTgeeAhwGPD/JT3hwvnBCkrtn+cyqqmfPV4CSJEmSRk+SQ9rbH1XVdwYajCRJkqQpdVUI9b/As4C3JDm1qu6ZqnOStYC30Byn94MO4pMkSZK08ns3sBewC7Bpe40JsNMsnlXtGHeVkiRJkjSdE2lyh4MAC6EkSZKklVhXhVDH0xRCPR747ySvrqpLJ+qYZHvg48ATaBKLj3cUoyRJkqSVWFX9NsmzgJcAewKbA6vT5BoFXAzMdkcoSZIkSZrOHcDDgasHHYgkSZKkqXVSCFVVn0uyH7A/8DTg4iSXARcBN9G8tHgkzQru7XqGfrGqTu4iRkmSJEkrv6p6APh8ewGQZFl7e2hVXTGQwCRJkiSNsiXA9sAGgw5EkiRJ0tS62hEK4ADgA8ARwCLgKTy46GnM2PEUHwbe1Fl0kiRJkiRJkiRJD3U6sAPwQuDcAcciSZIkaQqLupqoqpZW1euApwIfYfkWsum5rml/e2pVvb5d7S1JkiRJU9kG2Ba4atCBSJIkSRpJHwSuA45I8seDDkaSJEnS5LrcEQqAqroM+CuAJKsD69MUQd1WVfd2HY8kSZKk4VZV1w06BkmSJEmjq6p+neQ5wH8AX01yAnAy8AOadxs10AAlSZIk/U7nhVC92sKnGwcZgyRJkiRJkiRJ0mSSLO39CLyivcZ+n2p4VdVA38VIkiRJC4l/fEuSJEmSJEmSJE1ufKXTlJVPkiRJkgan80KoJI8BDgF2ATYB1gSeV1XX9PR5MrAlcHdVndd1jJIkSZIkSZIkSa13DDoASZIkSTPTWSFUkkXAu4EjgUUsXzFRwGrjum8BfBl4IMk2VfWLruKUJEmSJEmSJEkaU1UWQkmSJElDYlGHc30MeBOwCnA98B+TdayqrwA/afv+aSfRSZIkSZIkSZKkBSvJ+9pr40HHIkmSJGluOimESrIH8Ir2478CW1fVS6cZdirNrlF/1MfQJEmSJEmSJEmSoDnR4g3ARr1fJlmS5MdJHjOYsCRJkiTNVFdH4x3etv9ZVW+b4Zjvtu2T+hCPJEmSJEmSJEnSTGwFFLDaoAORJEmSNLWujsbbhSZJ+OQsxvy8bTeZ/3AkSZIkSZIkSZIe5J623WjKXpIkSZJWWl0VQo2dp71kFmMeaNtV5zkWSZIkSZIkSZKk8a5p20OSZKCRSJIkSZqTro7G+w3NlrFrzWLMlm172/yHI0mSJGlUJFkVeGz78cdVde+439cA/gV4Kc3K7iXAcVX14U4DlSRJkrSyOx14CnAY8PwkPwHu7/n9hCR3z/KZVVXPnq8AJUmSJE2tq0KoJcAOwFOBC2c4Zp+2vaIvEUmSJEkaFfsBnwduBraY4PfTgb2AsRXdTwA+mOSxVfWGbkKUJEmSNATeTZM77AJs2l5jAuw0i2dVO6bmLTpJkiRJ0+qqEOpsmiKoVyf5aFUtm6pzkqcBf06TIPxXB/FJkiRJGl7PpXnB8MWquq/3hyR7t78X8HPgImBnYHPgtUm+UFUzXawhSZIkaYRV1W+TPAt4CbAnTd6wOvAsmpziYmC2O0JJkiRJ6lBXhVAfBl4PbAd8IsnhVXX/RB2TvBj4KM1RencAH+8oRkmSJEnD6Q9oXkqcP8Fvh7XtVcDOVXVnkvWAC2h2hnolM9+1VpIkSdKIq6oHaHac/fzYd0nGFncfWlWeYiFJkiStxBZ1MUlV/YKmECrAocBPkhzX0+UVST6S5GrgFOARNC8yXl1Vd3QRoyRJkqShtXHb/qT3yySLaFZxF/DhqroToM0xPkyTn+zaYZySJEmSJEmSJKmPOimEAqiqT9Kstv4NzXayf8nys7GPBF4NPJrmZcS9wMur6tQuYkuyZZL3Jrkyyd1Jbk3y3SRvTrLWPM+1c5Lj2rl+neSuJD9OclaSNyX5vfmcT5IkSVoANmrb3477fgfg4e39WeN+u7xtt+hXUGPMNyRJkqShtw2wLc1Osysdcw5JkiRpua6OxgOgqj6V5Gyawqd9gceM6/IL4AzgPVV1bRcxJdkb+BywXs/XawE7tdcrk7ygqn4y0fhZzLM6zarzV9AUe/Xatr1eQLOKffGKzCVJkiQtMPfR5DYbjfv+D9v25xPkF3e27Sp9jMt8Q5IkSRoBVXXdoGOYjDmHJEmS9GCdFkIBVNXPgTcDb07ycJpjLFYBbqmqm7uMJcn2NEfxrQXcBfwb8A1gTeBA4FXA44GzkuxUVXfNcZ7VgNOB57df/Q/waeBK4AFgK2B74CVz/sdIkiRJC9e1wO8DTwe+3vP9C2l2oT1/gjEbtu2v+hWU+YYkSZKkfjLnkCRJkh6qk0KoJFu2t3dV1a1j31fVr4FfdxHDJD5AkyA8AOxVVRf2/HZukquBo4EnAG8C3jnHed7G8gThzVV1zLjfv0OTrFWPHFYAACAASURBVPxdklXnOIckSZK0UH0DeBLwuiSnV9WVSfYF9mh//88Jxjy5bW/oY1zmG5IkSZL6yZxDkiRJGmdRR/NcCyyhWYGwUkiyE8tfjHxyXIIw5hiaFQ0AR87lD/gk2wJvbT+eOEGC8CBVdf9s55AkSZIWuA/RHI+3MXB5kptpViuH5vjt0yYYsxfNblHf60dA5huSJEmS+smcQ5IkSZpYV4VQv2nbizqabyZe1HN/wkQdqmoZzfauABuwPKmYjVcDq9K8ZJnragtJkiRJk6iqq4E/B+6hKX7asG1vBw6qqvt6+yfZBHhO+/H/s3fncbdXY+PHP9fppEmlNIpKJzQpGoSfUgo9QoNS5ic04hHx4PHjUX7hMRbhkbFkKCpRUnGopFElMhQNCnVOdUrTaTrX74+1trPbZ9/DPvfe3/vs+/68X6/9Wt+99lrf73X+u6/zvdZa5wwoLPMNSZIkSYNkziFJkiR10cjReJRV2LOApRp63nhsV9t7gV+PMu7ctuvn0fuLktaZ2Jdl5vUAETEDeAIlebglM+8fabIkSZKksWXm9yLiXGBXYC3KkXc/bD+au83mwLfr9ewBhWS+IUmSJGmQzDkkSZKkLpraEers2j6voeeNx8a1/XNmPjzKuD92mTMuEbE6sEH9emFErBQRRwG3ATcB1wH/jIhzI2LXXu4tSZIk6dEyc05mfj0zP5qZ3xihCIrMPDsz96uffw4oHPMNSZIkSYNkziFJkiR10dSOUEcD/w68KyK+k5l/a+i5XUXEssBq9evNo43NzHkRcS+wAvCkHh+1Sdv1/cDllJ2x2s0Etge2j4jPZOY7e3wGEfHEMYas1es9JUmSJC0e8w1JkiRJgzSVcg7zDUmSJPVbI4VQmXltRLwaOAG4KCLeA3w/Mx9s4vldrNh2fc84xreShMf2+JxV264PBZYBfgW8D7gUWBb4N+CTwNrAOyLimsz83x6fc1OP4yVJkqQpLSLWBDZj4d/kdwC/y8xbG3i8+YYkSZKkQZpKOYf5hiRJkvqqkUKoiJhdL+cCTwa+CXw1Iq4F5gGPjDI9M3OnPoe0bNv1eIqxHqjtcj0+Z4W262Uo53TvlJnza9/9wLcj4lLgijr+iIg4zjO1JUmSpN5ExAzgQOAQHr1yuX3M74EvAF/KzAUDCsV8Q5IkSdIgmXNIkiRJI2jqaLwdgGz7HpQ/mjcbZU7WcTnKmMU1v+36MeMYv0xte/3DfX7H9/e3JQj/UnfM+iLwLmB1YGfgRz08Z6ztbNeirM6QJEmSpqS6A9QZwDNbXSMM3QQ4BnhjRLwsM28ZQDjmG5IkSdIUEhFLA0+pX/+SmQ90/L4scCTwSsqRddcDX8jMYwYU0lTKOcw3JEmS1FdNFUKdx2AKmhbX3W3X49kKtrXqYTxbzI70nAeBn48y9ixKkgCwDT28mMjMUc8AjxjpHZAkSZI0/OpLh9nARpQCqLnAScAlwK21bw3K39mvrNdbAT+NiK06X2L0gfmGJEmSNLXsAXwHuI3uhTunAi9i4YKMjYCjI+Ipmfn2AcQzZXIO8w1JkiT1WyOFUJm5QxPPGa/MnB8Rt1FWZjxxtLERsQoLk4Rez6puH39rZo62RW372DV6fI4kSZI0nb0D2Jiy+OKrwKGZeW+Xcd+MiPcCnwH2r3PeAXysn8GYb0iSJElTzospRU6ndP7dHRG71t8TuJmye9GzgHWAt0bEdzPzwn4GY84hSZIkjWxGP28WES+vnxXGHj3p/lDbDSNitIKwjbrMGa9rgYfq9VJjjG3//eEenyNJkiRNZ/tQXjqck5n7j1AEBUBm3peZBwJnU15k7DugmMw3JEmSpKljS0rOcV6X3/ar7TXAppn5CmAzFv59/+YBxWTOIUmSJHXR10Io4AfAKcB67Z0R8bX6WbvPz5uIX9Z2BcqxGCN5ftv1Bb08IDMfAlorPdYco0BsVtv133p5jiRJkjTNbVjbL/QwpzV21qijFp/5hiRJkjR1tHY4uq69MyJmADtTiqSOycy7ATLzLuAYyuKL5w4oJnMOSZIkqYt+F0LBwjOw2/078AZglQE8b3H9oO16v24DahLz+vr1TkY//3okJ9d2KWC3Ucbt2XZ9/mI8R5IkSZquHqhtL8c8tMaOdrTDRJhvSJIkSVPHarWd39H/DGClen1Gx2+/q+2TBhSTOYckSZLURb8LoVovIB7b5/v2XWZewsI/xt8UEc/pMuwwYON6fXRd/fAvEbFDRGT9fGOER30NmFOvPxIRa3YOiIgdgNfVr7+jx1UZkiRJ0jT3x9r28oKhNfaPo45aTOYbkiRJ0pTSWkCxWkf/9rW9OTNv6Pjt7tqOdaTcYjHnkCRJkrrrdyFUa7vT7fp830F5O3A/MBM4OyLeFxHPjogdI+JLwMfruGuATy3OAzLzHuA/KFvjrgdcGhGHRMTWEfG8iDgS+DElGXoYOCgzc2L/LEmSJGla+QZlZ9qDephzEOVv9OMHEVBlviFJkiRNDTfUdtuO/pdR/hY/r8ucVWs7d0AxgTmHJEmStIiZfb7fz4D9KasCnkX547p9hcEhETGn68xRZOYRfYqv875XRMQ+wAmU7Ws/0mXYNcCurbO9F/M5J0bEasCnKSvPP99l2D3AazPTlRKSJElSb75COYbhxRHxBeCdmdl5ZAUAEbEM5QXALsBZmfmlQQVlviFJkiRNGT8HNgXeFhGnZuYfIuLlwA719x93mbNZbf8xqKDMOSRJkqRF9bsQ6v9RXkA8Htir47cADl7M+w6kEAogM38UEZtTVk7sCjyRss3tn4HvAcdk5n19eM7nI+IXwFuAFwLrAI8A1wE/AY7KzIElRJIkSdIUth3lP+RXBQ4Edo+Ik4BLKUc4JLAmsA2wN7BW/e1TEbF91zsCmdltVXdPzDckSZKkKeFzwAHAGsDvImIesArlvcfNwMld5ryIkotcNsjAzDkkSZKkR+trIVRm3hQRWwIfAHai/CH8GMof+1E/S5zMvBF4Z/30Mu8X9PBvysyrgUN6Ck6SJEnSWH5ByTla1gTeNsacrYGzRvk96VO+ZL4hSZIkDbfMvDYiXgd8DViBhcfe3Qm8KjMfbB8fEWtRioUAzmkgPnMOSZIkqer3jlBk5k2UlRH/EhELKC8Snp6Zv+/3MyVJkiRNe0vkogtJkiRJU0Nmfi8izqXsurQW5ci7H2bmHV2Gbw58u17PbihESZIkSQygEEqSJEmSGrbjZAcgSZIkaerLzDnA18cx7mzg7MFHJEmSJKlTU4VQ+9X25oaeJ0mSJGmayMxzJzsGSZIkSZIkSZI0+RophMrM45p4jiRJkiRJkiRJ0qBFxJrAZsCqtesO4HeZeevkRSVJkiTJo/EkSZIkSZIkSZLGEBEzgAOBQ4BNRhjze+ALwJcyc0GD4UmSJEkCZkx2AJIkSZIkSZIkSUuyugPUJcAxlCKoGOGzSR1zcUSsNTnRSpIkSdOXO0JJkiRJmjIiYhbwcmALYDVgOcrLiJFkZu7URGySJEmShlNELAvMBjai5BdzgZMohVG31r41gG2AV9brrYCfRsRWmfnAZMQtSZIkTUcWQkmSJEkaehGxPPB54HUsWvgUQHbpo0u/JEmSJHV6B7AxJX/4KnBoZt7bZdw3I+K9wGeA/eucdwAfaypQSZIkabrzaDxJkiRJQy0iAjgVeD0lx7kd+A0Li53OB64GHm7r+xNwLnBeo8FKkiRJGkb7UIqgzsnM/UcoggIgM+/LzAOBsyn5x74NxShJkiQJC6EkSZIkDb+9gRfW68OBtShFUQBk5vMzc3NgVcpq7Hvr9Qcyc8eGY5UkSZI0fDas7Rd6mNMaO6vPsUiSJEkahYVQkiRJkobdq2t7YWYenpkL6HLkXWbem5lHAzsBKwKnRMQTGoxTkiRJ0nB6oLY39TCnNfbBPsciSZIkaRQWQkmSJEkadltTCp++PJ7BmXkp8EVgNeA/BhiXJEmSpKnhj7V9Ug9zWmP/OOooSZIkSX1lIZQkSZKkYbdaba9r63uodRERy3WZc0ZtXzqooCRJkiRNGd8AAjiohzkHURZsHD+IgCRJkiR1t8QUQkXEyyLimxFxZkR8ISKeOdkxSZIkSRoKD9f27ra+9uu1usy5q7a9rOiWJEmSND19BTgLeHF9f7HsSAMjYpmIOAbYBTg7M7/UVJCSJEmSYGYTD4mIHYETgfnA5pl5Z8fvHwb+q2PamyNiv8z8VhMxSpIkSRpafwdmAau39d0C3A8sC2wJXN8xZ8PaNpITSZIkSRpq2wGfBlYFDgR2j4iTgEuBOZSdn9YEtgH2pizGuBT4VERsP9JNM/O8AcctSZIkTTtN/af/SyjHVXy/SxHU5pQiqKhd84BVamzHRsQvM/PGhuKUJEmSNHx+QymEejpwDkBmZkRcDDwfOAQ4uTU4ImYC76hfr202VEmSJElD6BeUYqeWNYG3jTFna8ouUiNJXJghSZIk9V1TR+M9j/JH/TldfjuYUgQ1D9gqMx8PPAu4g7J6u5cztyVJkiRNP7MpOcUuHf1fq/07RMS5EfGWiHg3cDEl50jgpEYjlSRJkjSsYgAfSZIkSX3W1GqDtWr7xy6/vZTyAuLzmXkFQGZeVs/Q/iCwM/C+RqKUJEmSNIxOBY4BdoyIDTLzOoDMPCEiXk0pkHpe/bS7knK8hSRJkiSNZsfJDkCSJEnS+DRVCLVGbe9q74yIWcA6lEKoUzrmnF/bDQcbmiRJkqRhlpm3RMTSQGTmIx0/7wH8X+BNLFygcSfwLeD9mTm/uUglSZIkDaPMPHeyY5AkSZI0Pk0VQrW2eF25o3+72t6VmVd2/HZ7bZcfWFSSJEmSpoTMXDBC/wPAB4APRMSqlBxobmZmk/FJkiRJkiRJkqTBm9HQc26p7cYd/S+u7QVd5qxQ23kDiUiSJEnStJKZd2TmHIugJEmSJEmSJEmampoqhLqIsivUwRGxPEBEbADsRjkW75wuc55a21u6/CZJkiRJkiRJkiRJkiRJ/9LU0XhfAfYFNgd+FxGXA9sDywL3Ad/uMmf72v6+kQglSZIkSZIkSZJGERGzgJcDWwCrActRFoKPJDNzpyZikyRJktRQIVRmzo6Io4BDgfWB9ViYGLw7M29rHx8RyzL6blGSJEmSBEBEzF6MaQnMB+4CrqXsYntWZi7oZ2ySJEmSpoZ62sXngdexaOFTUHKMzj669EuSJEkaoKZ2hCIz31lfUOwNrAX8Azg+M7u9tHg58E/KSwkLoSRJkiSNZgfKy4VeXj506781Ig7LzO/0PUJJkiRJQysiAjgV2JmSS9wG3Aw8g5JTnA+sAjwNWLr2/Qm4ZTLilSRJkqazxgqhADLzdOD0cYw7CThp8BFJkiRJmgLOo7xoWBt4au1L4Dpgbv2+OrABC4ulrgFuBVaqc5ajLNg4ISKelJkfbyx6SZIkSUu6vYEXUnKJw4EPA5sAVwFk5vMBImIF4M3191WB/TPzl5MRsCRJkjRdzZjsACRJkiRpIjJzB+AjlGKnO4C3A6tl5lMy87n18xRgNcpx3fPq2I9m5jOBlYF9KCu6AzgyIjZp/l8iSZIkaQn16tpemJmH1yO1FznyLjPvzcyjgZ2AFYFTIuIJDcYpSZIkTXuTWggVETMjYvX6aXR3KkmSJElTQ0TMAr5PeRHxnMz8XGbO6xyXmfMy87PAc+rYkyLiqZn5cGZ+D9geuJOSJx3S3L9AkiRJ0hJua0oO8eXxDM7MS4EvUhZj/McA45IkSZLUofFCqIjYNCI+FxF/AOZTzsi+BZgfEX+ov23WdFySJEmShta7KKutP5aZ1441uI75OPDYOrfVfwPwJcquUDsOJFJJkiRJw2i12l7X1vdQ6yIilusy54zavnRQQUmSJElaVGOFUBExIyKOAq6krK5+Wn1+1M+M2ncIcEVEfCYiPLpPkiRJ0lheRFmdfX4Pc86t7c4d/bNru85Eg5IkSZI0ZTxc27vb+tqv1+oy567aPmkgEUmSJEnqqsnj6L4LvIJS9ARwNXAJcGvtWwPYBtgMWIqyXewTgH0ajFGSJEnS8HnCBOZ2vrCYU9tlJnBPSZIkSVPL34FZwOptfbcA9wPLAlsC13fM2bC2Tb6HkSRJkqa9RnZciohXA3vVr78Bts3Mp2fmmzLzvzLzffV6c+BZwBWU4qi9ImLfJmKUJEmSNLRaK62362HO9h1zW1ao7e0TikiSJEnSVPKb2j691ZGZCVxcvx7SPjgiZgLvqF/HPL5bkiRJUv80dfTc/rW9BnheZl460sDMvIzyUuJPlGKoAwcfniRJkqQhdgEld3hvRMwaa3BEbAi8l3Kc3q86ft60trf2NUJJkiRJw2w2JefYpaP/a7V/h4g4NyLeEhHvphRIPYuSc5zUaKSSJEnSNNdUIdTmlD/4/ycz7x1rcB3zP/XrFoMMTJIkSdLQO4qSb6wCXFRfPqzcOSgiHhcRbwUurGMT+HTHsJfSvUBKkiRJ0vR1KiVP2DEiNmh1ZuYJwE8oxVDPAz4LfAx4Rh1yJYvmHJIkSZIGqKlCqMfU9qoe5rTGLt3nWCRJkiRNIZl5PvBBysuHVSkvH26LiGsi4pf1cw0wFzgaeHyd+oHMvKB1n7qb1K71Pmc2+W+QJEmStOTKzFso7yqWzczrOn7eAziSsqts1M9dwOeBHTNzfpOxSpIkSdPdzIaecyOwMbDIquxRrNQ2V5IkSZJGlJlHRsR1wGeANYClgA2B1lF50TZ8DnBoZn634x5/obkcSZIkSdIQycwFI/Q/AHwA+EBErErJKeZmZjYZnyRJkqSiqf/kP5mSCLwC+Pk45+xF2Wr21EEFJUmSJGnqyMzvRMTJlBXZOwGbUY7AA5gHXA38DDi1vqyQJEmSpL7JzDsmOwZJkiRpumuqEOrTwGuAAyPivMw8abTBEbEXcCBwPfDJBuKTJEmSNAVk5oPAifUjSZIkSZIkSZKmkRlNPCQz7wJ2Bi4HvhMRP4iI3SNinYhYOiJm1uvdI+JUykuLy4Gd6lxJkiRJkiRJkiRJkiRJGlEjO0JFxCPtX4GX1c+IU4CtgesiYqQxmZlN7WglSZIkSZIkSZKmoYiYvRjTEpgP3AVcC1wEnJWZC/oZmyRJkqRHa6qQqLOaacTqph7HSJIkSZIkSZIkDdIOlMKmqG271ruM8fTfGhGHZeZ3+h6hJEmSJKC5QqjDG3qOJEmSJEmSJElSP51HKWhaG3hq7UvgOmBu/b46sAELi6WuAW4FVqpzlgPWAk6IiCdl5scbi16SJEmaRhophMpMC6EkSZIkSZIkSdLQycwdIuKFwHeBOyiLv0/IzHnt4yJiFeB1wAcphVGHZuZPImImsAfwKeCJwJERcXpm/r7Jf4ckSZI0HcyY7AAkSZIkSZIkSZKWVBExC/g+Zaen52Tm5zqLoAAyc15mfhZ4Th17UkQ8NTMfzszvAdsDd1LezRzS3L9AkiRJmj4shJIkSZIkSZIkSRrZu4AVgY9l5rVjDa5jPg48ts5t9d8AfIlyfN6OA4lUkiRJmuYshJIkSZIkSZIkSRrZiyg7PJ3fw5xza7tzR//s2q4z0aAkSZIkLWpm0w+MiACeAWwBrAYsR1n9MKLMPKKB0CRJkiRJkiRJkjo9YQJz1+r4Pqe2y0zgnpIkSZJG0GghVES8AfhvYL0ep1oIJUmSJEmSJEmSJsNdwOrAdsDF45yzfdvcdivU9vY+xCVJkiSpQ2NH40XEkcDXgPUpO0CN9qHLd0mSJEkal4hYJyK2iojtImK5yY5HkiRJ0lC7gPKu4r0RMWuswRGxIfBeynF6v+r4edPa3trXCCVJkiQBDRVCRcS2wPvq13MoR+NtWb8nsBTlmLxdgNMoCcUvgbUzs7FiLUmSJEnDKyJWjIjDI+Im4K/AJcAvgCd3jNs3Ik6KiC9PQpiSJEmShs9RlHcZqwAXRcRbImLlzkER8biIeCtwYR2bwKc7hr2U7gVSkiRJkvqgqaPxDq7tjcCumflwRLRWPZCZCdwBnA2cHREHA58HfhIR22bmgw3FKUmSJGkI1RXXZwIb8OhdZbPL8AuBbwIzIuK4zPxlAyFKkiRJGlKZeX5EfBD4MLAq8FngqIi4HphTh61BWYQxg4U5yQcy84LWfepuUrvW389sKHxJkiRpWmlqt6XnUl5AfDYzHx5rcGZ+ETgZ2Bw4ZMCxSZIkSRpiEbEMcAYwC7gP+DhllXVXmXkj8PP69eUDD1CSJEnS0MvMI4HXAHMphUxLARsCz6mfDWtfUIqjXp2ZH+24x18yc2ZmzsjMHzcZvyRJkjRdNLUj1Nq1vbqtb0HrIiKWzsyHOuZ8E3gFsA9l21lJkiRJ6uYg4CnAvcB2mXklQESMNudMYGfKCwtJkiRJGlNmficiTgb2AHYCNqMcgQcwj/IO5GfAqZn5wOREKUmSJE1vTRVCLV3bOW1997Rdrw78vWPOTbXdcFBBSZIkSZoS9qTsQHt0qwhqHK6q7VMGE5IkSZKkqSgzHwROrB9JkiRJS5imjsabW9uV2vpuBR6p1xt3mdPaRWrFQQUlSZIkaUrYpLZn9zDn9to+rs+xSJIkSZIkSZKkSdJUIVTrSLyNWh111USrf58uc15T286doiRJkiSpXWvxxF09zFm2tp1HdEuSJEmSJEmSpCHVVCHU+UAAO3b0n1j73xgRR0TEphGxTUQcA7yKcrzFmYMOLiLWjYhPRsQfIuLeiLgjIi6JiHdFxPIDeubyEXFdRGT93DCI50iSJEnTQGt3pzV7mPP02t7a51gWYb4hSZIkaZDMOSRJkqSFmiqE+kFtXxoR7cfjHQ3cUON4P3AVcBFwcP19HvDRQQYWEbvW5x5G2bFqeWAVYBvgE8DlEbHBAB59BPDkAdxXkiRJmm6urO1OPcx5I2XhxcX9D2ch8w1JkiRJg2TOIUmSJD1aI4VQmXk1ZTeoPYCZbf331f4LKDtDtX9+B+yUmTcPKq6I2AI4CVgZuIdSjPVcyguUL9dhTwPOiIjH9vG5zwQOBeYDd/frvpIkSdI0dQolhzgwItYba3BE/Dewbf164qCCMt+QJEmSNEjmHJIkSdKimtoRisw8NzN/lpl3dPTfmJnbARsDewH7Altl5uaZ+ZsBh3UUZXXEw8CLMvMjmXlhZs7OzAOA/6zjNgLe2Y8HRsRSlARkKeAjwB2jz5AkSZI0hm8AfwBWBM6NiF0jItp+z4iYERHbRcSPgA9SdoO6NDN/OMC4zDckSZIkDZI5hyRJktShsUKosWTmnzLzlMw8KTOvGPTzImIbYIf69auZeWGXYZ+ivFABODQilu7Do98ObAX8CfifPtxPkiRJmtYy8xHg5cBcYF3gh8A/24b8CLgL+AXwEsruUf8A9h5UTOYbkiRJkgbJnEOSJEnqbokphJoEu7ddf73bgMxcABxfv67CwqRisdRjOo6oXw/OzAcncj9JkiRJRWb+BXgGcAal0GmF+lMAG9TvrWO4zwa2ycybBhiS+YYkSZKkQTLnkCRJkrqYOdkBTKLtansv8OtRxp3bdv084JwJPPMLlBcw38zMn0/gPpIkSZI6ZOYtwMsiYlNgN2BrYA3KkQ23A1cAp2XmZQ2EY74hSZIkaZDMOSRJkqQuGimEiogPLsa0BOZTjrC4Fvh1Zv5z9Ck92bi2f87Mh0cZ98cuc3oWEftSjuGYB7xrce8jSZIkaXSZeTVw9SSHYb4hSZIkaZDMOSRJkqQumtoR6kOUwqaJeCgiTgPen5l/nsiNImJZYLX69ebRxmbmvIi4l7LK4UmL+bxVgKPq1/dm5pzFuc8o93/iGEPW6ufzJEmSJI3MfEOSJEnSIE2lnMN8Q5IkSf3W5NF4Udtsux5JtzGPAfYCdo2I3TLzZxOIZcW263vGMb6VJDx2MZ/3CWBN4ELgy4t5j9HcNIB7SpIkSUMhIhYAC4DNM/P345wzi7Lz7ILM7HdeZL4hSZIkDaGIWHcQ983Mv/b5llMp5zDfkCRJUl81UgiVmTNqAvFd4NnAqcDxwGXAXErh0xrA1sAbgN2Bi4FXA3cAm9XrA4Dlge9HxIaZeftihrRs2/WD4xj/QG2X6/VBEbE98EbgYeCgzJzozliSJEmSFjXWYot+zxuN+YYkSZI0nK4fwD2T/r+LMeeQJEmSRtBIIVRErAicDWwA7J2ZJ3cZ9rf6OS0iXkEpmjob2DozLwAuiIgfAqcDKwFvAY5YzJDmt10/Zhzjl6nt/b08JCKWAY6lvFw5OjOv6mV+D8baznYt4NIBPVuSJEkaRu071vab+YYkSZI0nAaxUGIQplLOYb4hSZKkvmrqaLxDgacCnxuhCOpRMvPkiPgi8FbgMOCDtf+siPgW8Hrg31j8Qqi7267HsxXsCrUdzxaz7d4PPI2yteuHepw7bpk56hngEcOSu0mSJEmNeXxt7x3Avc03JEmSpOG032QHME5TJucw35AkSVK/NVUItTdlpfWpPcw5hVIItSe1EKo6jVIIteHiBpOZ8yPiNmA14ImjjY2IVViYJPR6VvV7avtT4KUj/MHeuvcKEbFvvZ6TmbN7fJYkSZI03Y1rd6eIWAF4W/36l74HYb4hSZIkDaXMPG6yYxgPcw5JkiRpZE0VQj25tv/sYU5r7Hod/TfWdqUJRQR/ALYDNoyImZn58AjjNuqY04vWlrT7MfZKktWA79TrcwGTBEmSJKmLiLhuhJ/OjoiHxpi+DLAGMINSOPWjfsbWxnxDkiRJ0iCZc0iSJEldzGjoOa2XEU/vYU5rbOeLjFbMd04oIvhlbVcAthpl3PPbri+Y4DMlSZIkTdz6HR+AANbp8lvnZ21gqTr+YuDjA4rRfEOSJEnSIJlzSJIkSV00tSPUbykrE94VESdm5vzRBkfEcsC7KSu0f9vx86zazp1gTD8A3lev96O8BOmMYwblGD4ohVc/7+UBmTnm4dURcQNl16sbM3P9Xu4vSZIkTVOdx1W8gZI7tR5JngAAIABJREFU/JDRF0wkMB/4B/ArYHZmjus4vcVgviFJkiRpkMw5JEmSpC6aKoT6GqUQahNgdkTsn5lXdxsYEZsBX65jE/hqx5Cda/9VEwkoMy+JiPNrXG+KiOMy88KOYYcBG9frozPzUbtTRcQOLEwcjsvMf59ITJIkSZLGlpmPOpIhIt5QL9+fmb+fhJAWYb4hSZIkTU0REcAzgC0ox8EtR9lxdkSZeUS/4zDnkCRJkrprpBAqM4+LiL2AXYFtgasi4krg18CcOmwNYGtK8tByemYe3/oSESsDe1OSijP7ENrbKVvBLgecHREfofzRvxywL3BAHXcN8Kk+PE+SJElS/x1e2zmjjmqe+YYkSZI0hdRFGP9N2QGpF30vhKrMOSRJkqQOTe0IBbAncAzwZkoh0zPqp1NQdnw6Fnhbx28zgd3q9WUTDSgzr4iIfYATgJWAj3QZdg2wa2bePdHnSZIkSeq/zDx87FHNM9+QJEmSpo6IOBJ4L2Ps/lTlOMdNiDmHJEmStKgZTT0oMx/KzAMpuz59CfgLJRFo/1xHKYDaJjMP6tymNTNvz8xz6+fePsX1I2Bz4DOUhOA+ylnZlwHvAZ6ZmX/ux7MkSZIkTS/mG5IkSdLwi4htgffVr+dQFnlvWb8nsBTlmLxdgNMo7zt+CaydmQN9D2POIUmSJD1akztCAWWFAnAwQEQsAzyu/nRnZj7QdDw1phuBd9ZPL/N+wQRXdWTm+hOZL0mSJGlREbE+5UXEcozxN3tmnjfIWMw3JEmSpKF3cG1vpOyu9HBEbNr6MTMTuAM4m3JE3cHA54GfRMS2mfngIIMz55AkSZIWarwQql0tfLp1MmOQJEmSNDVExNOA/wJeTjkWYjySSc6LJEmSJC3xnkvJHT6bmQ+PNTgzvxgRLwD2BA4BjhpwfJIkSZKqxo7GkyRJkqRBiYjdgcuB1wIrs+gx3KN9JEmSJGk0a9f26ra+Ba2LiFi6y5xvUvKNfQYYlyRJkqQOrnyWJEmSNNQi4knACZRj8P4GfAK4DziWsmp7Z2AVYGvg9cATgF8CHwIeaT5iSZIkSUOmVeg0p63vnrbr1YG/d8y5qbYbDiooSZIkSYuyEEqSJEnSsPsPYHngbmDbzPx7RGza+jEzf14vT4mIDwNfpazKflNmvqbxaCVJkiQNm7mUBRXtR3DfSllYMQPYmEULoVq7SK048OgkSZIk/YtH40mSJEkadjtTdn76QmZ2vnx4lMy8n3J83hXAvhHxigbikyRJkjTcWkfibdTqyMwH2/q7HX/XWnQxao4iSZIkqb8shJIkSZI07Nav7a/a+rJ1ERGP2gk3MxcAnwUCeOOgg5MkSZI09M6n5A87dvSfWPvfGBFHRMSmEbFNRBwDvIqSl5zZbKiSJEnS9GYhlCRJkqRht0Jtb2rru6/teuUuc1ort7cYSESSJEmSppIf1PalEdF+PN7RwA2Udy3vB64CLgIOrr/PAz7aUIySJEmSsBBKkiRJ0vC7q7bLtvXd3nY9q8uc1suL1QYSkSRJkqQpIzOvpuwGtQcws63/vtp/AWVnqPbP74CdMvPmxgOWJEmSprGZYw+RJEmSpCXan4DnABtQVl+TmXdHxI3AusCLgEs65uxc2zubClKSJEnS8MrMc0fovxHYLiKeBmxKee9ybWZe0WR8kiRJkgp3hJIkSZI07C6s7bM7+k+nrMR+d0S8oNUZEXsBhwJJWbktSZIkSROSmX/KzFMy8ySLoCRJkqTJYyGUJEmSpGH3Y0rB054RsVRb/yeA+4DHAudExNyI+CdwIrAcsKCOkSRJkiRJkiRJU0DjhVARsWNEHB8R10bE3RHxcERs0jFmu4g4JCJe23R8kiRJkobOL4DDga8D67Q6M/OvwN7AXZRCqcdTiqICeADYPzMvajpYSZIkScMlIpaKiO3rZ+VxjH9ca3wT8UmSJElaaGZTD4qI5YHjgD1bXbXNLsMfAY4BMiIuzsxrGwhRkiRJ0hDKzKQUQnX77cyI2JBSELUpJQe6FjgpM//WXJSSJEmShtjOwJnA7cC64xj/IHAysGpE7JyZPx9kcJIkSZIWanJHqBMpRVABXAp8cqSBmfkr4Lf16ysGH5okSZKkqSoz78jML2Xmf2TmIZn5GYugJEmSJPXglbX9XmbeP9bgzLyP8k4kgH0GGZgkSZKkR2ukECoi9gB2rV8PyMxnZ+Z/jjHtFEqS8PyBBidJkiRJkiRJkjSyrSinW8zuYU5rF6hn9T8cSZIkSSNpakeoN9T2hMz8yjjn/Lq2Gw8gHkmSJEmSJEmSpPF4Ym2v72HODbVdp7+hSJIkSRrNzIaesw1ltcSJPcz5R21X7384kiRJkoZNRKzbus7Mv3brXxzt95IkSZKkLpavbfQwpzV2pT7HIkmSJGkUTRVCPb62f1uMuU3tWiVJkiRpydZafZ08OpfpZVV2p857SZIkSVKnuZRdoTYCLhvnnI1qe/tAIpIkSZLUVVNFRnfXtpeVD7Nqa5IgSZIkCcqK6tZnpP7F+UiSJEnSaC6l5A6v72HOv1MWXlw+iIAkSZIkddfUyudrgW2BZwHnj3POK2r7m4FEJEmSJGnY7NdjvyRJkiT1w/eBPYGdIuKwzPzUaIMj4jDgBZRCqO81EJ8kSZKkqqlCqB8DzwYOiYjPZ+b80QZHxC6UQqgETm8gPkmSJElLuMw8rpd+SZIkSeqTE4H3AFsAH4+I5wJHARdl5kMAETETeA7wDmA3yvuN3wEnTErEkiRJ0jTV1NF4xwB3AusDp0TE47sNiohl60qJU2pstwBfbyhGSZIkSZIkSZKkR8nMBPYA/kE5Im934BfA3RHx94j4G3BP7dutjvkHsFudK0mSJKkhjRRCZeadwGuBBcCLgb9GxI/bhnwgIs4E5gAfB5YFHgJeM9buUZIkSZKmt4i4rn7eOtmxSJIkSZqaMvMG4JnAabUrgMcAawFr1+uov50CbFnnSJIkSWpQU0fjkZk/joiXAN8E1gB2oWwNC/DK2raShNuAV2XmL5qKT5IkSdLQeiKwFHDlZAciSZIkaerKzDnAHhHxNOAllMKo1erPtwGXA2dk5rWTFKIkSZI07TVWCAWQmedExAbAfpTtYbcGHld/vg+4Avgh8L+ZeXeTsUmSJEkaWrcA6wD3T3YgkiRJkqa+zPwT8KfJjkOSJEnSohothALIzPuAz9cPETETWCozH2g6FkmSJElTwsXAnsCmwK8nORZJkiRJkiRJkjRJZkx2AJn5sEVQkiRJkibgi5Rjtt8REUtPdjCSJEmSJEmSJGlyNL4jlCRJkiT1U2bOjoiPAu8DTo+IN2fmTZMdlyRJkqThEhHrtq4z86/d+hdH+70kSZIkDZaFUJIkSZKGWkR8EHgAuAp4IXBdRFxQv88DHhltfmYeMfAgJUmSJA2D62ubPPr9yfVdxo5X570kSZIkDVBf//iOiNn9vF+VmbnTAO4rSZIkaWr4EOXlArVdCtiufsbDQihJkiRJUI7c7qVfkiRJ0hKm36sQdqC8eOhHUtC6T441UJIkSdK015mD+KJCkiRJUq/267FfkiRJ0hKm34VQ52HhkiRJkqQGZeaMyY5BkiRJ0vDLzON66ZckSZK05OlrIVRm7tDP+0mSJEmSJEmSJEmSJEnSeLhyWpIkSZIkSZIkSZIkSdLQsxBKkiRJkiRJkiRJkiRJ0tDr69F4kiRJkiRJkiRJU0lEfHAxpiUwH7gLuBb4dWb+s6+BSZIkSVpEI4VQEbHuYkz7V5KQmQ/2OSRJkiRJU0xEPAZ4DbA7sAWwGrDcGNMyM10gIkmSJGk0H6K8s5iIhyLiNOD9mfnniYckSZIkqZum/sP/+olMjoibgYuAb2Tmmf0JSZIkSdJUERFPBX4APA2ISQ5HkiRJ0tTTyjOSsXOObmMeA+wF7BoRu2Xmz/ocnyRJkiRgRkPPiQl+nkRJEE6PiLMiYpWG4pYkSZK0hIuIFYAzgY0oLxx+AHy5/pzAh4FjKIsrWn2/Ag4Hjmg0WEmSJElDJzNnAOtTcooATgX2oLy7WBZYpl7vQclHArgYmAWsAmwHfBFYACwPfD8iHt/oP0KSJEmaJpraEWq/2h4EbEs58u4s4DJgbv1tdWBr4MWUpOES4FhgJWAzYDfK0RY7UxKJ5zcUuyRJkqQl20HAk4FHgBdn5uyI2BTYHyAz/7s1MCKeAZwAPBv4bmYeMwnxSpIkSRoiEbEicDawAbB3Zp7cZdjf6ue0iHgF8N06Z+vMvAC4ICJ+CJxOee/xFlyYIUmSJPVdIztCZeZxwHOAZwGnAetn5h6ZeWRmHls/R2bmHpRVFT+qY/9PZh6dmfsD6wLHU1ZSPC8i9mkidkmSJElLvJdRdnk6KTNnjzYwM68EdgTmAJ+OiK0aiE+SJEnScDsUeCrwxRGKoB6ljvkiZUeow9r6zwK+RXnP8W+DCVWSJEma3hophIqI3YEDgAuBPTNzzkhj6297ULaY3S8iXln75wNvBC6vQ/cdaNCSJEmShsUmtT21248REe3fM3Mu8GnKDrlvHWxokiRJkqaAvSmLL7rmHCM4pbZ7dvSfVtsNJxqUJEmSpEU1UggFHEJJEo7OzBxrcB1zFGVVxAFt/QuAL9f+rQcTqiRJkqQh87ja3tjW90Db9WO7zLmgth65LUmSJGksT67tP3uY0xq7Xkd/K29ZaUIRSZIkSeqqqUKozWv75x7mtMY+vaP/qto+fkIRSZIkSZoq7qtt+6KLO9uu1+0ypzV2rYFEJEmSJGkqeai2ne8rRtMa+1BHf+u9zJ1IkiRJ6rumCqFaKxtW72FOa+yKHf3317YzeZAkSZI0PV1f2ye0OjLzNuCO+vX/dJmzVW0fHGBckiRJkqaG31JOqnhXRCw71uCIWA54N2UBxm87fp5V27l9jVCSJEkS0Fwh1E21fU0Pc15X27929K9RW5MESZIkSQCX1bbz+OyfUV5WvDsi/rWjbESsD7yH8lLiygbikyRJkjTcvlbbTYDZEbHpSAMjYjNgdh0L8NWOITtTcpGrkCRJktR3TRVC/ZDyAuK1EXHoWIMj4h2Uoqmsc9ttW9sbkSRJkiQ4h5JvvLyj/7O13QC4JiK+FxFnAL8Bnlh/O7aZECVJkiQNq8w8DjiDkndsC1wVEb+OiGMj4v/Vz7ERcTkl33hWnXp6Zh7fuk9ErAzsXe9zZrP/CkmSJGl6mNnQcz4GvAF4PPCpiHgNcDzwa2BOHbMGZQX364Ata9/cOrfdvpQCqXMGHLMkSZKk4XA6cB6wVETMysy/AGTmBRFxBPBBYBVgzzo+avv1zPx249FKkiRJGkZ7AscAb6bkFM+on05BeYdxLPC2jt9mArvV68uQJEmS1HeNFEJl5u0R8ULgJ8CalEKnLUeZEsAtwC6Zece/OiM2AC6pn5MHF7EkSZKkYZGZ9wE7jPDbhyLifMrLik0pOdC1wPGZaU4hSZIkaVwy8yHgwIj4X+AAYCdgw45hf6Ec0X1sZl7e5R63A+cOOlZJkiRpOmtqRygy8zcRsTFwOPB6YOURht5F2S3qQ5k5r+Me1wH7DTRQSZIkSVNKZv6M8jJCkiRJkiYkM68ADgaIiGWAx9Wf7szMByYtMEmSJElAg4VQAJl5J/D2iPhPYCtgM8oRFQDzgKuBy0wWJEmSJEmSJEnSkqy+y7h1suOQJEmStFCjhVAtNTn4Vf1IkiRJ0mKLiLcAJ2bmbZMdiyRJkiRJkiRJmjyTUgglSZIkSX30OeAzEXEO8C3gB5l53yTHJEmSJGkKioilgN2BnSmnXqxaf7oD+B3wU0pO8sjkRChJkiRNbxZCSZIkSZoKZgK71M99EXEapSjqbF9ASJIkSeqHiNgFOBZYp727tgk8FzgAuDkiDsjMsxoOUZIkSZr2Gi+EiogtgO2ADYAVgaXGmJKZ+aaBByZJkiRpWG0LvBrYB1gLWAF4Vf3cFhEnAt/OzIsmL0RJkiRJwywiXgd8nVL41Cp+ugG4pX5fE1ivXj8JOCMi3pCZ32o+WkmSJGn6aqwQKiKeBnwNeHYv0yirKCyEkiRJktRVZl4KXBoRhwEvAF4D7AGsBKwOvAV4S0TcAJwAfCcz/zhJ4UqSJEkaMhGxHmUnqBnAvcBHga9k5pyOcasDbwbeBzwW+HJEnJ+Zf204ZEmSJGnamtHEQyJiHeA8ShFUa7XEvcDNwF9H+dxY20HHt25EfDIi/hAR90bEHRFxSUS8KyKWn+C9l42I3SLicxFxcb33Q7W9MCI+FBFr9+vfIkmSJE1XmbkgM3+amftRVmO/EjgNeIiSgzwZ+L/A1RFxWUQc2sTf4uYbkiRJ0tB7O7AMcA+wXWZ+pLMICiAz52bmRymnYtxT57x90MGZc0iSJEkLNbUj1PspK7ET+Arwycy8pqFnjyoidgW+Bazc1r08sE39vDkiXpKZ1y3GvTcHfkk5ArDTKpTCsGcD74yIN2fmSb0+Q5IkSdKiMvMB4PvA9yPiccBelOPztqcsCNkSeCbwceAxg4rDfEOSJEmaEl5Eeb/xicy8cqzBmfmbiPgk8CHgxcBhgwrMnEOSJEl6tEZ2hAJ2oSQJx2fmAUtQEdQWwEmUBOEeSsHWc4GdgC/XYU+jnOX92MV4xEosTBAuoGyH+0LKS5cXA18CHqljvh0R/7Z4/xJJkiRJI8nMOzPzK5n5AmA94D3AnZRdopYa1HPNNyRJkqQpY93a/rSHOed0zO07cw5JkiRpUU3tCPWE2h7f0PPG6yjKyoiHgRdl5oVtv82OiGspK8Q3At4JHNHj/RdQkpDDM/P3XX4/OyLOBE6lvID5XEQ8JTOzx+dIkiRJGkNEbAa8BngVj14tPSjmG5IkSdLU0FpA8UgPc1pjB7kg3ZxDkiRJ6tDUjlDzantnQ88bU0RsA+xQv361I0Fo+RTwh3p9aEQs3cszMvNXmbnPCAnC/2fvvsMmqeq8/7+/w0hGchBRRzGAYkBBQUUGMYLZR8F1ZUFdXNf1p64J3VVh3fXBrI9ZVx1QjGsAxYAJUAygomsAEWUUkCCigEOG7++POu1d9HS8p/vu6u7367rO1VVdp06drtN90x/mdFWrzvHAZ8vqzsB9hjmGJEmSpO4i4vYR8fKI+F/gp8DLqH6RHcDVwCfGdFzzhiRJkjQ7LiyPDxxin1bdP4y4L4CZQ5IkSepmqSZC/bA83nWJjjeIJ9SWP9ypQmbezMJVrLZkIVSM2rdqyzuP6RiSJEnSXIiIrSLiORFxKvBb4HXAblSTn24Gvgo8A9g+M58+pm6YNyRJkqTZ8S2qPHFEROzYr3JE7AQcASTwzTH1ycwhSZIkdbBUE6H+H1VIOHyJjjeIfcrjGuBHPeqdUlt+8Jj6skFt+eYxHUOSJEmaWRGxUUQcHBEnABcB7wYeRJV5AjgdeAGwY2Y+OjOPy8w1Y+ySeUOSJEmaHe+g+i69LfCDiHhKRKzXXiki1ouIpwLfA7Yr+7xzTH0yc0iSJEkdLF+Kg2Tm1yLiDcDLIuI9wP+XmTcsxbF72LU8npuZN/aod3aHfUZt3y7HkyRJktRHRBxL9WvoTVpPlcdfAx8DjsvMc5e4W+YNSZIkaUZk5s8j4lXAfwE7Ut1i+y8RcSZwCdWVn3YAdge2YCGTvCozfz6mbpk5JEmSpA6WZCJURBwC/BL4LtVVoR4bEf9D9YX46n77Z+ax/eoM2Z8NgW3K6gV9jv3niFhD9Y8qtxtlP0pf7g0cWFZ/0ete2z3a2KlPlR2G7pgkSZI0Pf6+tnwp8Engo5l5xiQ6Y96QJEmSZk9m/t+IuAJ4A7Ax1a3m9mur1poAdTXw0sx8zzj6MkuZw7whSZKkUVuSiVDAKqpfRLTcBnj+gPsmC/ewHpXNast/HaB+KyRsOspORMQGwH8DrUvovnKRTZ0/mh5JkiRJU2kN8DngOOBrmTnpWzGYNyRJkqQZlJnvjohPAYcBDwN2A7Yqmy8Hfg58HfhwZl42xq7MUuYwb0iSJGmklmoiFCz8EqIJNqwtXz9A/evK40Yj7sc7gT3K8jGZecKI25ckSZLmwXaZec2kO1Fj3pAkSZJmVJng9MZSJsXMIUmSJHWxVBOh7rhExxnUtbXl9Qeov0F5HNk/rkTEK4Bnl9UfAc9bh+b6Xc52B2AitwWRJEmSxq1hk6DAvCFJkiRpvGYpc5g3JEmSNFJLMhEqM3+3FMcZwlW15UEuBbtJeRzkErN9RcRzgNeV1V8Bj87MNYttLzN73gM8okkX45IkSZJmnnlDkiRJ0jjNTOYwb0iSJGnUlvLWeI2RmddGxGXANsBOvepGxJYshIR1vld1RDwNeHdZ/R3wsMz847q2K0mSJM26iDiktZyZx3Z6fjHqbY2CeUOSJEnSOJk5JEmSpO7mciJUcRawD3DniFiemTd2qbdL2z6LFhGPA44FlgEXAfv3+7WDJEmSpL9ZBWQpx3Z4fjHa2xoV84YkzagVR5x4i/XVRx84oZ5IkkYtIn47hmYzM3ceQ7tmDkmSJKmDeZ4I9R2qkLAJcD/gB13q7VtbPm2xB4uI/YFPUZ3zPwEPz8zfLLY9SZIkaU51uy9C0+6XYN6QJEmSps+KMbS52B9t9GPmkCRJkjoY6USoiPhmWczM3L/D84txi7ZG6PPAK8ryYXQICRGxDGjdZuMvwLcWc6CIeCBwPLABcCXwyMz8xWLakiRJkubYHYd8fpLMG5IkSdL0OWbSHRiCmUOSJEnqYNRXhFpZHtt/4bCyPDfMr7Rb9cfya4nMPD0ivk31i4lnRcQxmfm9tmovBnYty2/PzBvqGyNiJQvB4ZjMPLT9OBFxH+BEql9lrAEOyMwfjeyFSJIkSXMiM383zPOTZN6QJEmSpk9mHjbpPgzKzCFJkiR1NuqJUKfSeeJSt+cn7QVUl4LdCDgpIl5H9aV/I+Bg4PBS7xzgzcM2HhE7A18FtihP/TtwRUTs1mO3SzPz0mGPJUmSJKlxzBuSJEmSxsnMIUmSJLUZ6USozFw5zPOTlplnRsRBwEeBWwOv61DtHODAzLxqEYfYB9iutv7WAfY5CjhyEceSJEmS5lJEfIjqhxf/npkXDbjPtsDrqW7F/axx9Mu8IUmSJGmczBySJEnS2pZNugOTlplfAO5F9QX+HOBqqntl/xB4ObB7Zp47uR5KkiRJ6uPQUrYcYp9b1/YbG/OGJEmSNJuisnNE7FnKzhGx5P/mYuaQJEmSbmnUt8abSpn5O+BfSxlmv5OB6LF9FbBqHbomSZIkacqZNyRJkqTZERGPBP4FWAls3Lb56oj4FvDOzDxpqfpk5pAkSZIWNOaKUBGxdUQM8wtuSZIkSVqsDcvjdRPthSRJkqSpEBHrR8THgS8BBwCbUE0iqpdNgAOBL0fExyJi/Un1V5IkSZpXE70iVERsD7wWeBLlNhYRcSVwPPDqzPz9BLsnSZIkaXY9qDxeMtFeSJIkSZoWHwOeSDXh6Ubga8APgIvLc9sD9wceDtwKOIjq32CeOonOSpIkSfNq5BOhImIH4Mdl9bWZ+Z4u9e4EnArchlteenVz4BnAYyNi/8z8yaj7KEmSJGl6RcSru2z654i4tM/uGwA7A48DEjhtlH2TJEmSNHsi4kCqH3Qn8C3gmeV2dJ3q3h74EPBQ4MkRcUBmfmnJOitJkiTNuXFcEWpfYAfgeuBTPep9Atixtn4+8Afg7sBmVFeI+nhE3DMzbxxDPyVJkiRNpyOp/gGiLoDnDtFGANcCbxxRnyRJkiTNrkPL40+BR2XmDd0qZubvI+LRVFeLujdwGNXt9CRJkiQtgWVjaHNlefxWZv6pU4WIeAywB9U/XlxOFRzukJl7U02i+nCpelfgyWPooyRJkqTpFrWSpcQA5TpgNXAcsHdm/nSpOy5JkiRp6uxFlTne3GsSVEup8yaqDLLXmPsmSZIkqWYcV4S6N1Ug+FqPOk+vLb84M09qrWTmNRHxbKqJUrsBjwc+OYZ+SpIkSZpCmXmLH3RExM1UGWS3zPzlZHolSWq6FUec+Lfl1UcfOFT9pT62JKlxti2Pw+SNs8vjNiPuiyRJkqQexnFFqO3LY69fVq8sj1cAH2vfmJlJdQ/toJpYJUmSJEnd/L6U6yfdEUmSJEkzaU153HqIfbYqj1ePuC+SJEmSehjHRKjtyuNlnTZGxJ2oJksl8O0el5E9szzuONruSZIkSZolmbkiM++YmedOui+SJEmSZtKvyuNBQ+xzcNu+kiRJkpbAOCZCtW63t36X7Q+oLf+oRzt/KY+brHOPJEmSJEmSJEmSFucEqjtYHBYRh/arXOocRvWD8M+PtWeSJEmSbmEcE6FaV4K6a5fte9eWf9ijnc3K47Xr3CNJkiRJkiRJkqTFeQdwEdVkqA9GxIkR8eSI2CkiblXKTuW5LwEfLHX/ALxzgv2WJEmS5s7y/lWG9lOq29k9GTiuviEiAnhsWb0ZOK1HO3coj5eMuoOSJEmSpk9EfHMMzWZm7j+GdiVJkiTNiMxcExGPAb4ObAk8qpRuAvgz8JjMvHoJuihJkiSpGMdEqOOBA4DHR8QhmXlsbdtLqSY4JfCNzLyiRzutK0d5/2xJkiRJACupskSMoK1WOzmCtiRJkiTNuMw8MyLuCbwdeAKwXpeqNwGfA16UmRcuVf8kSZIkVcYxEeo44JXA7YEPR8TzgHOBXYF71+q9pVsD5cpRT6D6R4nvj6GPkiRJkqbPqThxSZIkSdKEZOYfgKdExA7AfsBuwFZl8+XAz4GTM/OiCXVRkiRJmnsjnwiVmVdHxEHAScCtgT1KgYVfbn8oM0/q0cwBwG2p/pHj66PuoyRJkqTpk5krJ90HSZIkScrMi4GPT7ofkiRJkta2bByNZubpwP2ATwPXUE2ACuB3wEuAw/s08aryeHFmekUoSZIkSZIkSZIkSZIkST2N49Z4AGTmb4CDImIZsC1wfWb+ecDd9y+PN46lc5IkSZIkSZIkSZIkSZJM3gQGAAAgAElEQVRmytgmQrVk5s3AJUPus2ZM3ZEkSZIkSZIkSVpLRDyktZyZp3Z6fjHqbUmSJEkar7FPhJIkSZKkcYqI26/L/pn5+1H1RZIkSdJUOxnIUpZ3eH4x2tuSJEmSNEZ++ZYkSZI07c5bh339RwlJkiRJdTHk85IkSZIaxP/hL0mSJGna+Q8SkiRJkkZhvyGflyRJktQwToSSJEmSNO0OG6DOJsDdgCcDOwLfBT4wzk5JksZjxREn/m159dEHTrAnlVZ/mtAXSdK6ycxThnlekqRRqGec+jKYMyRpMZwIJUmSJGmqZeYxg9aNiJcA/w84HPhuZr5sbB2TJEmSJEmSJElLatmkOyBJkiRJSyUzb8jM5wKnAi+OiEdOuk+SJEmSJEmSJGk0nAglSZIkaR69Bwjg+ZPuiCRJkqTpFxFbRMS2ERGT7oskSZI0z5wIJUmSJGke/bo87jHRXkiSJElqrIhYLyJ2i4g9ImK7Dts3iIgjI+J84E/AxcBfI+IzEXGPJe+wJEmSJCdCSZIkSZpLm7c9SpIkSRIAUXkNcBnwU+AHwEUR8b2IeECpcyvgy8CrgB2prjgbwEbAE4DTI2L/SfRfkiRJmmfLJ90BSZIkSZqAfyiPF020F5IkSZKa6APAYWW5fqu7BwBfjoi9gH8EVpbnL6e66uxy4O5Uk6E2Ao6LiLtl5hVL0WlJkiRJXhFKkiRJ0hyJiLtExHupJkIl8KUJd0mSJElSg0TEvsAzy+q1wP8AbyqP11BdVfaFpc4NwOHAtpm5d2buCWwDvLnsvy1w6JJ1XpIkSZJXhJIkSZI03SLitwNUWwZsAWxWe+5S4L/G0ilJkiRJ06p1JahLgYdm5lmtDRGxC/BNqqtBLQPemJn/Xd85M68BXhoRuwGPBA4E3r4UHZckTc6KI0782/Lqow8cqv5SH1uSZp0ToSRJkiRNuxWL2Of7wGGZ6a3xJEmSJNXtTXX12LfWJ0EBZObZEfFW4PWlzkd6tHMM1USoe4yro5IkSZLW5kQoSZIkSdPumAHq3AxcBZwHnJKZPxlvlyRJkiRNqR3L4/e6bK8/f26Pdn5dHrda5x5JkiRJGpgToSRJkiRNtcw8rH8tSZIkSRrIJlRXe7q8y/a/tBYy87oe7VxbHtcfUb8kSZIkDWDZpDsgSZIkSZIkSZLUMDnk85IkSZIawIlQkiRJkiRJkiRJkiRJkqaeE6EkSZIkSZIkSZIkSZIkTb3lk+6AJEmSJA0qIh4y6jYz89RRtylJkiRp6v1zRFza4fntWgsR8eoe+2/XY5skSZKkMXEilCRJkqRpcjKQI2wvMRdJkiRJWttze2xrZZLXLEVHJElLa8URJwKw+ugDJ9yTZvVFkqaF/8NfkiRJ0jSKSXdAkjQdWv9wMK76nfZbzD9SrOv+kqSRMm9IkiRJU8qJUJIkSZKm0TXA8cDXgJsn3BdJkiRJs2O/SXdAkiRJ0uI5EUqSJEnSNLkK2AzYCDgIWAl8DPhIZv50gv2SJEmSNAMy85RJ90GSJEnS4i2bdAckSZIkaQjbA08DvgTcBOwAvAj4cUT8NCJeEhE7TrKDkiRJkiRJkiRpMpwIJUmSJGlqZOa1mfnJzHwMcFuqSVBnAgHcE3g98LuI+FpEPCMiNplgdyVJkiRJkiRJ0hJyIpQkSZKkqZSZf8zMt2fmHsA9qCZBXQCsB+wPrAIuiYiPRMQjIyIm11tJkiRJkiRJkjRuyyfdAUmSJElaV5l5FvAK4BURsRI4BHgScGvg6cDfUSZFZebLJ9ZRSZIkSZIkzYwVR5w41vrd9l199IFLtq8kTRuvCCVJkiRppmTmyZn5TGAHqglQXwZuKuvPn2TfJEmSJEmSJEnS+DgRSpIkSdKsSuDm8pgT7oskSZIkSZIkSRozb40nSZIkaaZExL7AM4D/A2zWehq4CPjIpPolSZIkSZIkSZLGy4lQkiRJkqZeROxKNfnp6cBOraeBq4HPAccC38jMmyfTQ0mSJEmSJEmSNG5OhJIkSZI0lSJiO+BpVBOgdm89TXU7vG9RTX76bGaumUwPJUmSJEmSJEnSUnIilCRJkqSpEREbAk+gmvz0cGA9qslPAD+nuvXdcZn5h8n0UJIkSZIkSZIkTYoToSRJkiRNk0uBTcpyABcDHwc+kpk/mVivJEmSJEmSJEnSxDkRSpIkSdI02RRI4FrgBOAk4CbgXhFxr8U0mJnHjq57kiRJkiRJs2vFESd23bb66AMH3rdb3XVpf1Dd+tF+7F59GeYYozBo31rPrz76wL516trP7SBjtRSG7Uf99UuaX06EkiRJkjSNNgSeWsq6SMCJUJIkSZIkSZIkzQAnQkmSJEmaNjHpDkhSUy31r6dH2W63YwzyC+32Xzsvtk+j/tX2uhyjKb/AliRJkiRJmiZOhJIkSZI0TfabdAckaZYtZvJNEyfsDDLZaJyTnta17WH7Py23qmhKnyVJkiRJ0uxyIpQkSZKkqZGZp0y6D5IkSZIkSZIkqZmWTboDTRARt4+IN0XEWRGxJiIuj4jTI+IlEbHxCI/zqIj4bERcEBHXlcfPRsSjRnUMSZIkSc1i3pAkSZI0TmYOSZIkacHcXxEqIg4EjgM2rz29MbBnKc+OiAMy87frcIwA3gsc3rbptsATgSdGxPuBf8rMXOxxJEmSJDWLeUOSJEnSOJk5JEmSpFua6ytCRcS9gU9RBYS/Av8GPBDYH/hAqXY34MSI2HQdDvWfLASEM4GnAfcvj2eW5w8HXrsOx5AkSZLUIOYNSZIkSeNk5pAkSZLWNu9XhHob1S8jbgQekZnfq237ZkT8GngDsAvwr8B/DHuAiLgz8LKy+kPgIZl5TVk/IyJOAE4B9gBeHhEfzszfLOrVSJIkSWoS84YkSZKkcTJzqLFWHHHi35ZXH33gBHvSLO3nZZLnqX7sfn3oVLdbvU6vq9P+g7Y5yD71frfqDHp+F9OPfsbR5rj7MMwYd9M6x4t5/Z3eJ+1j1mts2/cdpJ/thvlMjEuT/kY0Xf39ILWb2ytCRcSewMqy+sG2gNDyZuCssvzCiLjVIg71IhYmnD2/FhAAyMyrgeeX1eXACxdxDEmSJEkNYt6QJEmSNE5mDkmSJKmzuZ0IBTyhtvzhThUy82bg2LK6JQuhYiDlvtmPL6tnZ+b3uxzn+8CvWv0q+0mSJEmaXuYNSZIkSeNk5pAkSZI6mOeJUPuUxzXAj3rUO6W2/OAhj3FH4LYd2ul1nJ2AFUMeR5IkSVKzmDckSZIkjZOZQ5IkSepgnidC7Voez83MG3vUO7vDPsMeo72dUR9HkiRJUrOYNyRJkiSNk5lDkiRJ6mB5/yqzJyI2BLYpqxf0qpuZf46INcAmwO2GPFS9fs/jAOd32a+viNipT5XWLza46KKLhmlakqbGjVde9rflCy7o9ydXkqZP2/e49SbVD/Vn3pA0SfXvxe0G+Z48yPfqdT3GILr1o9uxL7jggp79Wtf6k9bpvDYlAw3Tj6b0WdLazBvTZZYyh3ljdvnf/c7az8skz1On78OLyQCd2uj1OsehU2bodX6nKQs0waDnrlVvMee30/uk23tnXd9jw7zPJ/m5nPTfiKarf9bVfEudNyIzx32MxomIbYFLy+onM/PgPvUvAbYDfp6Z9xziOC8F3lBWH52ZX+lR99HAl8rqSzLzzUMcZ/4GUZIkab7tmZk/nHQn1Jl5Q5IkSVPOvNFws5Q5zBuSJElzZ+x5Y15vjbdhbfn6AepfVx43GuNxrqstD3scSZIkSc1h3pAkSZI0TmYOSZIkqYu5vDUecG1tef0B6m9QHq8Z43E2qC0Pe5x+l5ldH9iF6hcifwRuGrL9xdoBOKMs7wlcvETH1Wg5jrPBcZx+juFscBxnw6TGcT1g27L8syU6phbHvLE0/Js6GxzH2eA4zgbHcTY4jtPPvKFBzFLmMG9onBzH2eA4zgbHcTY4jrNhEuO4pHljXidCXVVb3nSA+puUx7+O8Tib1JaHOk5mDnLjy98O0+YoRER99eIB+6mGcRxng+M4/RzD2eA4zoYJj+PvlvBYWjzzxhLwb+pscBxng+M4GxzH2eA4Tj/zhgY0M5nDvKFxchxng+M4GxzH2eA4zoYJjuOS5Y25vDVeZl4LXFZWd+pVNyK2ZOEL/PlDHqr+hul5HG75q4dhjyNJkiSpIcwbkiRJksbJzCFJkiR1N5cToYqzyuOdI6LXlbF26bDPoH7ZpZ1RH0eSJElSs5g3JEmSJI2TmUOSJEnqYJ4nQn2nPG4C3K9HvX1ry6cNeYzzgD90aKeTh5THC4HVQx5HkiRJUrOYNyRJkiSNk5lDkiRJ6mCeJ0J9vrZ8WKcKEbEMOKSs/gX41jAHyMwEji+ru0TEXl2OsxcLv5Y4vuwnSZIkaXqZNyRJkiSNk5lDkiRJ6mBuJ0Jl5unAt8vqsyJi7w7VXgzsWpbfnpk31DdGxMqIyFJWdTnU24Aby/I7ImKjtjY2At5RVm8s9SVJkiRNMfOGJEmSpHEyc0iSJEmdze1EqOIFwDXAcuCkiHhFROwVEftFxPuAN5R65wBvXswBMvMc4E1ldQ/gtIg4KCL2iIiDqC5Fu0fZ/sbM/PViX4wkSZKkRjFvSJIkSRonM4ckSZLUZvmkOzBJmXlm+aL+UeDWwOs6VDsHODAzr1qHQ/0bsB3wTGB34BMd6nwQ+Pd1OIYkSZKkBjFvSJIkSRonM4ckSZK0tvBWzRARd6D65cSBwE7A9cC5wKeBd2bm1V32W8nCPbWPycxD+xznAOBwYE9gG+Ay4AzgfZn55XV+IZIkSZIax7whSZIkaZzMHJIkSdICJ0JJkiRJkiRJkiRJkiRJmnrLJt0BSZIkSZIkSZIkSZIkSVpXToSSJEmSJEmSJEmSJEmSNPWcCCVJkiRJkiRJkiRJkiRp6jkRSpIkSZIkSZIkSZIkSdLUcyKUJEmSJEmSJEmSJEmSpKnnRChJkiRJkiRJkiRJkiRJU8+JUJIkSZIkSZIkSZIkSZKmnhOhJEmSJEmSJEmSJEmSJE09J0JJkiRJkiRJkiRJkiRJmnpOhNJYRMTtI+JNEXFWRKyJiMsj4vSIeElEbDzp/s2iiLhvRLwyIr4cEedHxHUR8deIOCciVkXEPkO296iI+GxEXFDauqCsP2qINpZHxHMi4tSI+GNEXBMR50bEeyPi7sO/yvkVEW+IiKyVlQPs4xg2QERsExEvi4jTIuLiMhZ/iIgfRMQbI2LvAdrYOyI+EhGrI+LaiLgoIr4SEQcP2ZeDI+KrZf9rS3sfiYi9Fv8KZ19ErB8Rzyrn/KLa39dfRcSHBj1/fiZHLyK2i4jHRMR/lP/+XVb7O7lqEe01ZowiYuuIOCoifhoRV0TElWX5qIjYetjXJs2aMG8suTBvzLwwc0ylMG9MvTBvNFaYN6S5FeaNiQgzx0wL88ZUCvPG1AvzRmOFeWO8MtNiGWkBDgT+AmSXcjZwp0n3c5YKcEqP810vxwLr92krgPf1aed9QPRpZ2vg+z3auBZ45qTP3TQU4N7ADW3nb6Vj2PwCPAW4rM9YfL5PG68Gbuqx/wnAhn3a2BD4Qo82bgJeNenz1cQC3A743wH+vr6l22fKz+RYx6fXOV01RDuNGiNgT+APPdq5ENhj0uffYplUwbwxiXNu3pjxgpljKgvmjakvmDcaXfqc01VDtNOoMcK8YbH0LJg3JnXezRwzXDBvTGXBvDH1BfNGo0ufc7pqiHYaNUY0JG9MfIAts1WovsysKW/iq4BXAnsDDwXeX3uDnwVsOun+zkoBzq394Xgb8OTyR2Yv4EXABbVz/7E+bf1Xre6PgYNLWweX9da2/+zRxnrcMrh8BngUcH/g+cAl5fkbgUdO+vw1uVBdue/0cr4uqZ3TlY5hswtwCAtf8C8BjgQeBtwXOKCcx5OAT/do49m1MTgXeGYZy8cD36xt+2ifvhxXq/vNsv+epb1za9uePenz1qQCLOeWIeGnwD+Uv60PB44C/lrb/tIu7fiZHN8Y1b88/x74am191RDtNGaMgNsCF5e6NwCvB/Yp5fUs/E+ji4HbTnoMLJalLpg3JnXezRszXDBzTGXBvDH1BfNG40vtfJg3LJY5KZg3JnnuzRwzWjBvTGXBvDH1BfNG40vtfJg3xnF+Jz3AltkqwLdqb+y9O2x/ae0D9OpJ93dWCvBF4KnAel22bwP8qnbu9+lS7861P0BnABu1bd+4PN8a4527tHNo7Vjv6nKcK8r2c4Dlkz6HTS3AC1kI16+rndeVjmFzC7Ar1czoBE4FNu9Rt+MvmIAtgD+XNn4HbNO2fT2qX0u0xukhXdrZt1bnhPa/E+Xvw+/K9suBLSZ9/ppSqP6HS+vcfbfT31jgfsD1tfO3vG27n8nxjtFRwGOA7cv6itp5WjVgG40aI2BVrZ2ndNj+lNr2D016DCyWpS6YNyZ13s0bM1wwc0xdwbwxEwXzRuML5g3zhmXuCuaNSZ57M8eMFswbU1cwb8xEwbzR+IJ5Y6x5Y+IDbJmdQjWzsPXGfW+XOsuAX9b+oN5q0v2el1L+kLbG5+1d6ryrVmevLnX2qtV5R5c6v6iN8cZd6hxRa+fJkz4/TSxUl6y8qpyjlVQz7vuFBMewAQX4ejkvf6TtC/4QbdT/x8rBXersRDX7OoEvdKlzYtl+I7BTlzoH14714kmfv6YUqsvBts7LY3vU+2yt3m5t2/xMLu2Yraidg1UD7tOYMQK2r32mv9Kjz1+pfa63n/R5t1iWqmDeaHTBvDGVBTPHVBbMGzNRMG9MXcG8YbHMdMG80fiCmWPqCuaNqSyYN2aiYN6YuoJ5Y6RlGdLoPKG2/OFOFTLzZqp7OANsSfXFR0vj5Nryzu0bIyKoLicJcHZmfr9TI+X5X5XVJ5T96u3cBbh7Wf1kZl7dpT+rastP6tnz+fVuYFPgmMw8uV9lx7AZImIXYP+y+s7MvGyRTbX+pl5J9UV0LZl5AVUoAXh4RGza1pdNa335WqnfyWfLccCxrFu/tvzbHvV+U1veoLXgZ7L5GjhGj6P6NRR0+S7V1s56ZR9pXpg3mu3k2rJ5Y3qYOaaMeWOmmDdmXAPHyLwh9WbeaL6Ta8tmjulg3pgy5o2ZYt6YcQ0co0blDSdCaZT2KY9rgB/1qHdKbfnB4+uO2tT/g3dzh+13pLpvJ9xyjDppbd+JanZq3T4d6q0lMy+munQe+D5YS0Q8leoXLpdTzZwfhGPYDE+pLX+6tRARW0bEXSJi634NRMT6VPfdBfheZl7fo3prjDag+uVa3f1Z+OLaayyvB1pfkO4fEbfq18c5cU5t+U496rX+x0sCv64972ey+Zo2RgO1g9+lNL/MG81m3pgyZo6pZd6YHeaN2de0MTJvSL2ZN5rPzDFFzBtTy7wxO8wbs69pY9SovOFEKI3SruXx3My8sUe9szvso/Hbt7Z8doftu/bZTpft7WO4mHZuFxGb9Kk7NyJiC+DtZfXlmfnHAXd1DJthr/J4BXBWRDw9In5KFfjOAS6LiN9GxGvaf+FQcxdgeVle6rFcXo4v+DgLvyR5eUSs114hInYHDiyrn8jMK2ub/Uw2X9PGqNXOFSVUdJSZF7Hw3vS7lOaJeaPZzBtTxMwx1cwbs8O8MfuaNkbmDak380bzmTmmhHljqpk3Zod5Y/Y1bYwalTecCKWRiIgNgW3KardLEwKQmX+m+lUFVPcH1phFxDKq+3a2fKpDtfpY9BxD4Pwu+y22naCagarKG4AdgO8CHxxiP8ewGVqXj1wNvAP4KHCvtjp3pLof+vciYscObUxyLDu1M5dKQD8UuAZ4EHBGRBwSEXtFxMMi4jVUM9fXB34C/GtbE34mm69pY9Rqp18b9Xb8vGoumDeazbwxlcwc08u8MSPMG3OhaWNk3pC6MG80n5lj6pg3ppd5Y0aYN+ZC08aoUXnDiVAalc1qy38doH4rKHSbLazRehELl6H8XGb+sEOdYcZwTW25fQxH1c5ciogHA88GbgT+KTNziN0dw2bYqjzuAjwP+AvwT8B2wIZUl3f9cqmzG/DpEuTrHMuGyMzPAXtQBfb7AMcA3wO+RhX2rqYKCA/uMMPdcWy+po1Rqx2/S0lrM280m3ljipg5pp55Y4aYN2Ze08bIvCF1Z95oPjPHlDBvTD3zxgwxb8y8po1Ro/KGE6E0KhvWlnvd67XluvK40Rj6opqI2Bc4uqxeCjy3S9VhxvC62nL7GI6qnblT7pv8fqpZtG/NzJ8N2YRj2AytS0FuANwEPDoz35eZf8zM60pIfwwLYeGBwJPa2nAsG6LcT/zvgMdSfTbbbQ88DVjZYZvj2HxNG6NWO36XktZm3mgo88Z0MXPMBPPGDDFvzLymjZF5Q+rOvNFgZo7pYd6YCeaNGWLemHlNG6NG5Q0nQmlUrq0trz9A/Q3K4zVj6IuKiLgH8Dmqe+JeBzw1My/pUn2YMdygttw+hqNqZx69kupeqL8HjlrE/o5hM9TP36cz8/vtFTLzZuCltaee1qMNx3JCyv2Nvw78G7A11SWdd6U6V5sDjwC+Q/UrmC9ExAvamnAcm69pY9Rqx+9S0trMGw1k3phKZo7pZ96YEeaNudC0MTJvSN2ZNxrKzDF1zBvTz7wxI8wbc6FpY9SovOFEKI3KVbXlQS5h1ppRPMil0bQIEXFH4CRgS6pZ20/LzFN67DLMGG5SW24fw1G1M1ciYhfgFWX1+Zm5plf9LhzDZqifvy93q5SZvwAuLKt79mjDsZyco4CHlOVnZebLM/PszLw+M6/MzK8B+wHfovo1xVsion6/dMex+Zo2Rq12/C4lrc280TDmjelj5pgZ5o3ZYd6YfU0bI/OG1J15o4HMHNPFvDEzzBuzw7wx+5o2Ro3KG8vH1bDmS2ZeGxGXAdsAO/WqGxFbsvDmPn/cfZtHEbEj1SzfHYEEnlnuA9vLBbXlnmMI3K623D6G7e1cNkA72bbfPHoR1QzZ3wIbR8TBHersVlt+aETsUJa/UEKFY9gM5wOtsel3Ts4Hbkt1f+26cY3lDxfZztyJiAAOK6vnZOYxnepl5o0R8SqqX04sK/u8qGz2M9l8TRujC6guR9yvL/V25v7zqvlg3mgW88bUMnPMBvPGDDBvzI2mjZF5Q+rCvNE8Zo6pZN6YDeaNGWDemBtNG6NG5Q0nQmmUzgL2Ae4cEcsz88Yu9XZp20cjFBHbAF8D7lSeen5mHjvArr+sLe/Stdba29vHsL2dnwzQzvmL/HXALGldAvBOwMcHqP+q2vIdgTU4hk3xCxZ+AbFen7qt7e1/L8+h+pXTeox2LAdp50bg3D5158H2wFZl+cw+dX9UW66fZz+Tzde0MfolcD9g84jYITMv7tRARNwGuHWXvkizzLzRAOaNqWbmmA3mjdlg3pgPTRsj84bUm3mjIcwcU8u8MRvMG7PBvDEfmjZGjcob3hpPo/Sd8rgJ1Zu8m31ry6eNrzvzJyI2B74K3L08dURmvmvA3c8D/lCW9+1VkYVLKV4IrG7b9p3actd2ykz/u5ZV3wej4Rg2w6m15Z371G2F+QvrT2bm9cDpZXXviOh1P93WGF3H2r+IOAO4vq3eWkr7e7X2Kcefd/Xw1m/i+K267OdnsvmaNkYDtYPfpTS/zBsTZt4QjmMTmDdmg3ljPjRtjMwbUm/mjQYwc8w9x3DyzBuzwbwxH5o2Ro3KG06E0ih9vrZ8WKcKEbEMOKSs/oXqvqMagYjYGDgRuG956r8y8/WD7p+ZCRxfVneJiL061SvPt2Z7Hl/2q7dzDguzN59a+tXJobXlfpe0nXmZeWhmRq9CdT/flv1q21aXNhzDZjgBuKEsP6lbpYjYF9i6rH67Q5XW39Rbd2snInYCHlZWv5GZ9fv4Uta/UVYfVup38iQWZl87lpXLgSvL8t4R0Sss1L+0ndda8DPZfA0coxOAm8tyx+9Sbe3cXPaR5oV5Y4LMG9PPzDEzzBuzwbwxBxo4RuYNqTfzxoSZOaabeWNmmDdmg3ljDjRwjJqVNzLTYhlZoZopnFT/kdy7w/aXlu0JHDnp/s5Kobrv8ldr5/Zti2znrmXskmqm9UZt2zcqz7fG+C5d2nlmrS/v7LB9Z+CKsv1cYPmkz+E0FODI2nld6Rg2twDvrp2/gzts34zqcqStOnt2qLMV1f9QSarZ2Vu3bV+P6gtCq439uvTlobU6xwPrtW3fBvhd2f5nYMtJn7+mFOBjtXP3mi51tqS6XHCr3iPatvuZXNoxW1E7T6sG3KdRYwQcW2vn/3TY/pRhX6PFMksF88akzrt5Y04KZo6pKJg3ZqJg3pi6gnnDYpn5gnljkufezDEHBfPGVBTMGzNRMG9MXcG8MdrzOekBtcxWAXYHri5v3quAV1BdjnA/4H21N/avgM0m3d9ZKcBnauf2G8A9gd16lLv2aOv/1tr6MXAQsEd5/HFt2+t6tLEe1eXvWnX/B3gkcH/gX4BLyvM3AY+e9PmblsIAIcExbEYBtmXhy/cNwDvK38H7Uc10Pqt2bt/do53n1OqdSzWDeg/gccA3a9s+1qc/H6/V/WbZf4/S3rm1bc+Z9LlrUqGaIb+mdn5OAJ5c/lu3N/Ci2jgn8PUu7fiZHN8YPbh8plrlJbVz9J22bYf2aKcxYwTcDri09vfj6PI6H1yWW6HmUmCnSY+BxbLUBfPGpM67eWNOCmaOqSiYN2aiYN5ofMG8Yd6wzF3BvDHJc2/mmIOCeWMqCuaNmSiYNxpfMG+MNW9MfIAts1eAx7IwI7BT+RVw50n3c5ZKj3Pdrazu0dYy4IN99v9vYFmfPm1DdQ/gbm1cB/zjpM/dNBUGDwmOYQMKsCvw6z7j8EHgVn3aOYrq8pDd2jgR2LBPGxuVet3auAl/xdbt3D0M+OMAf1e/QZdfm/iZHOv4rBpgbP5WerTTqDECHgBc1KOdi4AHTPr8WyyTKpg3JnHOB/5bi3ljqgtmjqkpmGsKy6oAABSfSURBVDdmomDeaHTBvDHxMbBYJlEwb0zqvA/89xYzx9QWzBtTUzBvzETBvNHognljvOd30gNsmc0C3AF4C1UoWEN1OcIzgJcBG0+6f7NWhvkjSZ+QUGvzAKp7+F5Y/rBdWNYHnokLLAeeS3V/4MuAa4DfAO8H7jHp8zZthQFDgmPYnAJsQjWD+/vAn8o4nA98gi6Xeu3SzgOB44DflzYuAU4CnjZkf/6u7HdJaef3pd21LvVtucV527r89+tbVLPUr6f6deBvgU8CjwdigHb8TI5+bFYN89+/aRojqtDxWuBnVL9CvQr43/Lc1oO2Y7HMasG8sdTn27wxJwUzx1QVzBszUTBvNLZg3rBY5rZg3pjEOTdzzEHBvDFVBfPGTBTMG40tmDfGWqJ0RJIkSZIkSZIkSZIkSZKm1rJJd0CSJEmSJEmSJEmSJEmS1pUToSRJkiRJkiRJkiRJkiRNPSdCSZIkSZIkSZIkSZIkSZp6ToSSJEmSJEmSJEmSJEmSNPWcCCVJkiRJkiRJkiRJkiRp6jkRSpIkSZIkSZIkSZIkSdLUcyKUJEmSJEmSJEmSJEmSpKnnRChJkiRJkiRJkiRJkiRJU8+JUJIkSZIkSZIkSZIkSZKmnhOhJEmSJEmSJEmSJEmSJE09J0JJkiRJkiRJkiRJkiRJmnpOhJIkSZIkSZIkSZIkSZI09ZwIJUmSJEmSJEmSJEmSJGnqORFKkiRJkiRJkiRJkiRJ0tRzIpQkSZIkSZIkSZIkSZKkqedEKEmSJEmSJEmSJEmSJElTz4lQkjQhEXFyRGREnDzpvoxCRKworycj4tB1aOfQWjsrRtZBDSwiVtbGYOU6trWqtLN6NL2TJEnSoMwcXdsxc0yYmUOSJGn6mTe6tmPemDDzhiQnQklqvLYvLBkRnxhgn9YXk1yKPkqSJEmaXmYOSZIkSeNi3pAkaWk5EUrSNHpqRNxz0p2YF7VwduSk+6JmiIgjDeGSJGnGmTmWkJlD7cwckiRpxpk3lpB5Q+3MG9LsWz7pDkjSIgRwFPCkSXdECzJzNdXYaMpl5smMaCwz81Dg0FG0JUmStITMHA1k5pgdZg5JkjTnzBsNZN6YHeYNSV4RStK0uaw8PjEi7jvRnkiSJEmaRWYOSZIkSeNi3pAkacycCCVp2rwduK4sHzXJjkiSJEmaSWYOSZIkSeNi3pAkacycCCVp2lwAvL8sPyYiHrCYRiJidbn/76o+9VaVeqv71Ns2Il4dEadFxKURcV1EnF/WXx0Rd1tMP0vbW0bEv0fE9yListL2HyLi+IgY6NK5EXHPiHh/RPw6Iq6OiKsi4hcR8daIWNFln9Vt90d+Te1e2q2yqlZ/Re35Q/u8nqMj4uyIuKacr69HxFMGOiEL7SyPiGdFxJfK+biunJ9TI+KFEbFhn/3vFxEfjIhzImJNRFxbxuxHEfGuiHhcRAx96dSIOLR2HlZExAYR8ZKI+HFEXBERV0bEDyLieRGx3gDtbRsR/xkRZ0bEX0o/V0fERyLiwQPs/9CI+HhEnFfO99Vl/+9HxJsi4qEd9llZew0r218b8Jrac+3viay/p4b4DA39Hi37rfW+i4iHR8QXIuLi8r44LyLeExE79TtfkiRJmDnMHAvtmDnMHGYOSZI0auYN80arHfOGecO8IY1LZlosFkujC7ASyFIOBW4DXF3Wv9pln1WtfbpsX122r+pz7FY7q3vUeTrw11ofO5W19gdOLttO7tH2AcCf+7T9RWDTHm28Aripx/7XAof0OEe9yqpa/RX1cerSl7sDf+jR3gfLGLfWV3RpZ2fgF336dg5wly77v6jPOWmVrue1x/mu93934Ic92v82sFmPth4BXNGnj+8ElnXZ/y0DvMbL+nzmVnZ5bb3Kito+q+j/GVrUe7TT+w44ukc7lwK7LsXfLYvFYrFYLNNVMHOYOdZux8xh5uj4vsPMYbFYLBaLZciCecO8sXY75g3zRsf3HeYNi2UkZTmSNGUy86KIeA/wr8AjIuLBmfmdSfQlIg4Bjimr1wIfAL4MXAxsCtwLeCxwl0W0/XDgBGA9qi/s7wF+AFwJ3BY4CPh74MDShyd3aOOfgdeV1T8CrwdOK20+DHgpsAmwKiIuy8wv1XZ/BLA+8LOy/h7g3W2H+PMQr2dz4KtUIQ/gk6XflwJ3pRrPZwL37NPObcpr2B64iurXM18HLgE2L/1+AdU5/0pE3Dczr6jtfy/gTVRXRTyP6kv2T4DLqcbsLsB+wBMHfW09vA+4X4fX+iJgT+DBwHHA4zq8zvsAX6AagxuAdwHHA2uowscRwB2B55XnXt62/2PKcQD+l2r8zqIKHZsDuwAPB/Ye4vV8nir0/DPw3PJcp/G6cNAG1/E92u4fgQcCp1Cd+3OALYBDStkW+BDDvWZJkjSHzBxmDswcZo7OzBySJGmdmTfMG5g3zBudmTekUZn0TCyLxWLpV2j7tUR5blsWfqHwzQ77rGrt06XN1azjryWAHam+nCXVF9TderSzU4fnTqbLryWovhRdXLZ/Fdi4S7v/WDs3+7dt27bWvwuB23XYf/faebwAuFWHOq32j+xzrla0j1Pb9jfXtr+iw/ZbldfacdZ9rd4XyrbfA3fq0pf663pt27b/KM//Fdi+x+vZnC6/QuhzHg5tew2dXuty4Cu1Ogd2qHN62XYj8IgO27dk4RcjNwH3aNt+bOu9S+9f02zV5zO3ssP2I1vbBzgfvT5D6/webXvfJVVojA7tfKBWZ/dhx9VisVgsFstsF8wcZo5b1jNzpJmjy/vOzGGxWCwWi2XognnDvHHLeuaNNG90ed+ZNyyWEZVlSNIUysw/Us1yB9gvIvabQDeeD2xclp+TmT/vVjEzLxiy7cOofg1wLfCMzLy6S7sfoPoy2dqnvY1W/16cmed32P9M4P+W1dsCTxiynwOJiA1q/ftfqhnx7X25AXgW1S8DurWzG/CYsvovmfnbTvXK63pXWX1m2+YdyuM5mXlJt2Nl5hWZeXO37QPq9lpvBJ7Nwmv95/r2iLg/1a8pAP47M0/q0MafgcPL6rL2Nlh4nT/OzL9262BmXt7nNYzTqN+jFwHPz8zssO1NteV9FtFXSZI0Z8wcf2vXzNGBmQMwc7Qzc0iSpIGZN/7WrnmjA/MGYN5oZ96QhuBEKEnT7I1Ulw0FeO0Ejn9geTyP6nKeo/T48nhKZl7ap+6p5bH9UpgPK49/AT7TY///7rDPqN2PanY/wDHdvnyXMLXWF+Ka1nm5GjixzzFb52XHiLhd7fmLyuPdy5fxcRr0ta6MiPVqm+vj8MFujWfmaVSXgm3fBxZe50MiYufBu7ykRv0e/Z/MvK7Thsz8FdWvLgDuNHAPJUnSvDNzVMwcnZk5KmYOzBySJGlRzBsV80Zn5o2KeQPzhjQsJ0JJmlqZ+SfgbWX1QRHxyKU6dkTcCtitrH67y+zsdbFHeXxkRGSvAryk1N2hrY1W/84sv0ToqPxiYHXbPqNWv8fyGX3qnt5jW+u8bAzc2Oe8fLG2X/3cfJzqVwobAKdFxBci4p8i4h4REQO+nkEN+lo35pZfXFvjcD1wZp82flAe7xIR69eeP7Y8bg38PCI+ERGHRcSd+7S3lEb9Hj27z/Fa93vfbKDeSZKkuWfmMHOYOQAzR52ZQ5IkjYx5w7xh3gDMG3XmDWlEnAgladq9hWqmNVT3RV4qWwGtL5QX9ao4rBJAtljErhu3rW9VHrteGrXm4rZ9Rm3L2nK/X3/06u92izz+385NZp4NPI3qC+NyqsvQvgf4OXBpRHwkIkZ1WdFhXutWHZYvL5eY7aU1dkHtPGfmN4B/Aa4BNgQOAj4E/DoiLoiI90bEvfu0PW6jfo92vLxyTeuXK+v1rCVJknRLZo4FZo7uzBxmDjBzSJKk4Zk3Fpg3ujNvmDfAvCENbPmkOyBJ6yIz/xIRb6EKCPePiMdk5hf77Tfqboy4vfoXmE+x7pfEHaR/o/6VQK/2+/WnV19a5+Y84HFDHP+8+kpmfiYivk71xfmRVPdT3hbYBvh74O8j4hjgmet4D+11ea2D7N+zjcx8V0R8Gvg74OHAg4DNqe5D/Rzg8Ih4XWb++wDHGacmvEclSZI6MnMMpAnf58wcnZk5Kk14j0qSJK3FvDGQJnyXM290Zt6oNOE9KqnGiVCSZsHbgBdQXR7zPyKi332VYWHWdL8r423S5fnLSxvLgB0H6eSgMvPaiLiaaob/Fpn580U2dTlwG9a+nGwn29f2GYd6u9sD5/So2+sXEX+qtXH2AL8k6CozrwDeXwoRcXeq4PF8qjH9B6pLtr59scdguNd6eYflrSNieZ/X2Rq7ZOGyqH9T7r/+NuBtEbEMuA/wJOB5VL/K+beIOCMzR30P+EE06T0qSZLUi5mjsyZ9nzNzdGbmaM57VJIkqRvzRmdN+i5n3ujMvNGc96ikGm+NJ2nqZeZVwBvL6u7AEwbY7aryuGXPWnC3Lse8geoyowD7jOG+y617Jj8oItovBzuoVv92L5ei7SgitgPu0LbPqP2strxnn7q9trfOy8ZUM/9HJjN/mZlHA3sBa8rTT13HZgd9rVcDv6093xqH9ane073cvzz+OjOv71UxM2/OzB+XX0fsX9s07Osc1S+EmvQelSRJ6srM0VWTvs+ZOXpvN3NM/j0qSZLUkXmjqyZ9lzNv9N5u3pj8e1RSjROhJM2Kd7Jwr+Kj6H+JydZlRO/b7Qt+ROwG3LNHG18oj3cEHj9gPwd1QnnchGpW+2J8vTxuATy5R71nsXC+vt5h+7XlcYNF9gPgRyzM5H9Gj3N+W+ARPdqpz+h/2Tr0p6vMPJ+FXzhss47NDfpaT87Mm2qb6+PwrG6NR8TewN077NNXZv6YhTEZ9nW23hNExLq8L0b1HpUkSVoKZo61mTkWwcwxEDOHJEmaN+aNtZk3FsG8MRDzhjTjnAglaSZk5hrg9WX1nsABfXY5pTzuCDytfWNEbAZ8qE8b72RhVv37SqjoKCJ26tNWu/cCl5Xl10bEo3tVjogHRcRD2p7+MNUsfIA3R8TtOux3b+CVZfVC4PMdmr+oPO48SMc7yczrSn+gumzpSzv0ZTnwAapfCHRr5wzgpLJ6QEQc1eu4EbEiIp7W9twTImKLHvvcDtilrJ7Xrd6ABn2t76lvz8zTgTPK6rMj4uEd2tgceF9Zvbm9jYg4KCI26taxiNiDhV8LDfs6L6otL/p9wejeo5IkSWNn5jBzdGLmMHNIkiSNgnnDvNGJecO8IWlxlk+6A5I0Qu8BXkJ1P95+s78/ChwJ3Br4YETcGfgq1eUw9wD+Fbgt1WVKO162MzMvjojnAsdS3Qf59Ij4APBl4GJgU2A3qnsy340hvkxl5pXly+2XqX6l8MWI+AzwGeA3pdptgPsBTwTuRXXf51NrbfwxIl4KvIsqDP0wIo4GvgusBzyM6gvspuV1H14uh9vuu1S/CHlcRDwHOI2F2fJXlvszD+I/qC5PuhPw+oi4D9W5uxS4K9U535Pqy3Gvy60eBvywvP5XR8QjqQLdz0q/ti7n41HAQ6m+VH68tv8LgeOius/6N4GzgCuovjDvQXUeW1+ub/HFexF+2OG13qW81tblXr+QmV/ssO/hwA+ogsSJEfEOql/o/JXqPXkEcKdS900d7rP+euC9EXE81fviHKpQuzXw4PI6AW6iCizD+G5t+a0R8V9UwaF1OdnVg9zbfITvUUmSpKVi5jBzmDkWmDkkSZJGy7xh3jBvLDBvSFq8zLRYLJZGF2Al1ReEBA7tU/dfanWz+jPXte5TgBvb65dyTdm+qqyv7tHOP1DN+O7UTnbbHzi5bDu5R9sPZeHLV79ySJc2Xkn1RbDbftd227fsf59Sp9O+q2r1VvQbJ+AefV7Ph4BDa+srurRzB+D0Ac/Lh7qc917lRuAVi3y/1vu/O/DjHsf5DrBZj7YeQRVgevX1ncCyDvuuHuB1XtNp7LnlZ25ll759ske7K2r1VtH/M7To9+gg77sO52RVr3oWi8Vi+f/bu3vVqoIoDKDffZIgWIidIIiCii9hoRYilnkBSW0ExcLGJpWCIFjYioV2FgoSULGw0aRUUAwoRnFbzAkx3MTERA1zXQt2lcy+5/4Us+Gbc5T6/ypmDjPHeJ+pmDlWysxh5lBKKaXUDirmDfPGeJ+pmDdWyrxh3lDqj5dH4wGTZi7J4lb+saruJDmS5G6St0mWh7U3khwc/r6VPjfSTkJcTHtO9Ieh10LaJnAmyYnfehervR8MvaeT3EvbYC+nbZoW026hOpNkX1Xd3KDHbNpmdS7tpMXntNT8yyTXfrV2WD+f5HDaiYOFJF+2816GXi/SBoXLSV4Nvd4leZjkVFWd22KfN0kOpZ0UuZ1229NPSb6mfZePklxNcjzjz58+meR02uZ1Pu1ky7e0UwjPk1xPcqCqLm3zbf7sfdpv7MLwWkvD6zxJO61wvKqWNlpcVfeT7E0yO6z/mPaZLSS5leRoVU1X1fd1lh9Lcj5tM/8s7XP5NvR4muRKkv2/+u43cSbtGeaP0waZ9a5hS3b6GwUA+MfMHOM9zBxrmTnMHAAA22XeGO9h3ljLvGHeADYxqqrdvgYAmBij0ehsVp8VvqeqXu/e1QAAAJPGzAEAAPwt5g1gErgjFAAAAAAAAAAA0D1BKAAAAAAAAAAAoHuCUAAAAAAAAAAAQPcEoQAAAAAAAAAAgO4JQgEAAAAAAAAAAN0bVdVuXwMAAAAAAAAAAMCOuCMUAAAAAAAAAADQPUEoAAAAAAAAAACge4JQAAAAAAAAAABA9wShAAAAAAAAAACA7glCAQAAAAAAAAAA3ROEAgAAAAAAAAAAuicIBQAAAAAAAAAAdE8QCgAAAAAAAAAA6J4gFAAAAAAAAAAA0D1BKAAAAAAAAAAAoHuCUAAAAAAAAAAAQPcEoQAAAAAAAAAAgO4JQgEAAAAAAAAAAN0ThAIAAAAAAAAAALonCAUAAAAAAAAAAHRPEAoAAAAAAAAAAOieIBQAAAAAAAAAANA9QSgAAAAAAAAAAKB7glAAAAAAAAAAAED3fgBdCgzyRf0QkwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 2400x600 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Gaps number 51\n",
      "wilcoxon WilcoxonResult(statistic=2652.0, pvalue=4.108463998517357e-13)\n",
      "ks_2samp Ks_2sampResult(statistic=0.114, pvalue=4.441088228134397e-06)\n"
     ]
    }
   ],
   "source": [
    "fig, axes = plt.subplots(nrows=1, ncols=3, figsize=(12, 3), dpi=200)\n",
    "x = range(window_size)\n",
    "\n",
    "single_gap = np.zeros((window_size,))\n",
    "single_gap[window_size//2] = 1\n",
    "\n",
    "ax1, ax2, ax3 = axes.ravel()\n",
    "\n",
    "ax1.bar(x, single_gap, width=5)\n",
    "ax1.set_title(f\"Single gaps frequencies ({assembly.assembly})\")\n",
    "ax1.set_ylabel(\"Single gaps frequency\")\n",
    "ax1.set_xlabel(\"Nucleotides position\")\n",
    "\n",
    "ax2.bar(x, synthetic_mean, width=5)\n",
    "ax2.set_title(f\"Multivariate gaps frequencies ({assembly.assembly})\")\n",
    "ax2.set_ylabel(\"Multivariate gaps frequency\")\n",
    "ax2.set_xlabel(\"Nucleotides position\")\n",
    "\n",
    "ax3.bar(x, mean, width=5)\n",
    "ax3.set_title(f\"Biological gaps frequencies ({assembly.assembly})\")\n",
    "ax3.set_ylabel(\"Biological gaps frequency\")\n",
    "ax3.set_xlabel(\"Nucleotides position\")\n",
    "\n",
    "fig.tight_layout()\n",
    "fig.savefig(\"hg19_distributions.jpg\")\n",
    "plt.show()\n",
    "\n",
    "print(\"Gaps number\", number)\n",
    "for test in (wilcoxon, ks_2samp):\n",
    "    print(test.__name__, test(mean, synthetic_mean))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}