{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Q1 vs Q2 contact analysis\n", "\n", "Here we calculate a Q1 vs Q2 plot, where Q1 refers to fraction of native contacts along a trajectory with reference to the first frame, and Q2 represents the fraction of native contacts with reference to the last.\n", "\n", "**Last executed:** Feb 06, 2020 with MDAnalysis 0.20.2-dev0\n", "\n", "**Last updated:** January 2020\n", "\n", "**Minimum version of MDAnalysis:** 0.17.0\n", "\n", "**Packages required:**\n", " \n", "* MDAnalysis (Michaud-Agrawal *et al.*, 2011, Gowers *et al.*, 2016)\n", "* MDAnalysisTests\n", " \n", "**Optional packages for molecular and data visualisation:**\n", "\n", "* [matplotlib](https://matplotlib.org)\n", "* [pandas](https://pandas.pydata.org)\n", "\n", "
\n", " \n", "**Note**\n", "\n", "The `contacts.q1q2` function uses the `radius_cut_q` method to calculate the fraction of native contacts for a conformation by determining that atoms *i* and *j* are in contact if they are within a given radius (#best_native_2013)\n", "\n", "
" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import MDAnalysis as mda\n", "from MDAnalysis.tests.datafiles import PSF, DCD\n", "from MDAnalysis.analysis import contacts\n", "\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Loading files\n", "\n", "The test files we will be working with here feature adenylate kinase (AdK), a phosophotransferase enzyme. (Beckstein *et al.*, 2009) The trajectory ``DCD`` samples a transition from a closed to an open conformation." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "u = mda.Universe(PSF, DCD)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculating Q1 vs Q2" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We choose to calculate contacts for all the alpha-carbons in the protein, and define the contact radius cutoff at 8 Angstrom. `contacts.q1q2` returns a `contacts.Contacts` object, which we can run directly." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "q1q2 = contacts.q1q2(u, 'name CA', radius=8).run()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The data is in `q1q2.timeseries`. The first column of the data is always the frame number." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
FrameQ1Q2
00.01.0000000.946494
11.00.9809260.949262
22.00.9736600.952952
33.00.9727520.951107
44.00.9700270.948339
\n", "
" ], "text/plain": [ " Frame Q1 Q2\n", "0 0.0 1.000000 0.946494\n", "1 1.0 0.980926 0.949262\n", "2 2.0 0.973660 0.952952\n", "3 3.0 0.972752 0.951107\n", "4 4.0 0.970027 0.948339" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "q1q2_df = pd.DataFrame(q1q2.timeseries, \n", " columns=['Frame', \n", " 'Q1', \n", " 'Q2'])\n", "q1q2_df.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plotting" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can plot the fraction of native contacts over time." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "scrolled": false }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Fraction of native contacts')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydeXhU1fn4P+9k3wghC1tCEkIIuxAQwYXFlUVFcKUWrbZabeluW/1qbau1ta3aTbS1rb8WrQuitirIIqKgiOw7BEICJIFASMi+Z87vjzOTTJJZbiAhC+fzPPeZmXPvufdMGO57312UUhgMBoPB0BJbZy/AYDAYDF0TIyAMBoPB4BYjIAwGg8HgFiMgDAaDweAWIyAMBoPB4Bb/zl5AexETE6OSkpI6exkGg8HQrdi6detppVSsu309RkAkJSWxZcuWzl6GwWAwdCtE5KinfcbEZDAYDAa3GAFhMBgMBrcYAWEwGAwGt/QYH4TBYDB0FHV1deTm5lJdXd3ZSzlrgoODiY+PJyAgwPIcIyAMBoPBB7m5uURERJCUlISIdPZy2oxSisLCQnJzc0lOTrY8r8NMTCLysoicEpE9HvaLiPxZRDJFZJeIpLvsu1tEDjm2uztqjQaDwWCF6upqoqOju6VwABARoqOj26wBdaQP4l/ADC/7ZwKpju1+4EUAEekD/By4BJgI/FxEojpwnQaDweCT7iocnJzN+jtMQCil1gFFXg6ZAyxWmo1AbxHpD1wHrFZKFSmlzgCr8S5oADhZWk1BWU17LN1gMBguDLLXe93dmVFMA4Ecl8+5jjFP460QkftFZIuIbDlVVkPOmcoOW6zBYDB0Nrm5ucyZM4fU1FQGDx7MwoULqampobCwkOnTpxMeHs7ChQutn/DQSq+7u3WYq1LqJaXUBKXUBIDiytrOXpLBYDB0CEop5s2bx0033cShQ4c4dOgQVVVV/OQnPyE4OJgnn3ySZ555pm0nLcr2urszBUQekODyOd4x5mncJ0UVde22OIPBYOhKfPzxxwQHB3PPPfcA4Ofnxx/+8AcWL16MUorLL7+c4ODgtp20KMvr7s4Mc30PWCgib6Ad0iVKqRMishL4tYtj+lrgESsnNBqEwWDoaH75/l72HS9t13OOGNCLn98w0usxe/fuZfz48c3GevXqRVJSEpmZmYwdO7ZtF7XbfWoQHSYgROR1YBoQIyK56MikAACl1F+B5cAsIBOoBO5x7CsSkSeBzY5TPaGU8ubsbuSMERAGg8FgjfJ8qK/yekiHCQil1Hwf+xXwbQ/7XgZebsv1/G3CmUpjYjIYDB2Lryf9jmLEiBEsXbq02VhpaSn5+fmkpaW1/YSFh30e0q2d1K742cSYmAwGQ4/lqquuorKyksWLFwPQ0NDAj370IxYuXEhISEjbT+jD/wA9SED424QzxkltMBh6KCLCu+++y9KlS0lNTSU6Ohqbzcajjz4K6J44P/zhD/nXv/5FfHw8+/bt837Coiywea/L1GNqMfnZbMYHYTAYejQJCQm89957AGzYsIH58+ezbds20tPTOXLkSNtOVpQFUUlAocdDepCAECMgDAbDBcOll17K0aMem8H5pigb+gwGtno8pGeZmCrr0L5vg8FgMHhEKa1B9Bns9bAeIyD8bEJtvZ2quobOXorBYDB0bcpPQV3FhSUgABPqajAYDL5wRjBdKALC3ykgKowfwmAwGLzSKCC8Nw/qMQLCqUEUGw3CYDAYvFOUBeIHvQd5PazHCAh/m/4qRSaSyWAw9FA8lftevXo148ePZ/To0YwfP56PP/7Y+4mKsrRw8POeB9FjBISfn1ODMALCYDD0PLyV+46JieH9999n9+7d/Pvf/2bBggXeT2Yhggl6koBo9EEYE5PBYOh5eCv3nZqayoABAwAYOXIkVVVV1NR46LCplEsOhHd6TKKcABFB/iZZzmAwdCwfPgz5u9v3nP1Gw8ynvR5itdz322+/TXp6OkFBQe5PVFkENSUXloAAiAoLNCYmg8FwwbJ3715++tOfsmrVKs8HWQxxhZ4mIEIDTB6EwWDoWHw86XcUvsp95+bmMnfuXBYvXkxKSornE7VBQHSoD0JEZohIhohkisjDbvYnisgaEdklIp+ISLzLvt+KyB7HdruV6/UONRqEwWDomXgr911TU8Ps2bN5+umnueyyy7yfqCgLEIhK9HnNDhMQIuIHLAJmAiOA+SIyosVhzwCLlVJjgCeA3zjmzgbSgbHodqQPiUgvX9eMCg0wYa4Gg6FH4q3c9/PPP09mZiZPPPEEY8eOZezYsZw6dcr9iYqyIDIB/D34KFzoSBPTRCBTKZUF4Og9PQdwLVI+Avih4/1a4L8u4+uUUvVAvYjsAmYAS7xdsHdoIMUmislgMPRQPJX7fuyxx3jsscesnaQoy2cGtZOONDENBHJcPuc6xlzZCcxzvJ8LRIhItGN8hoiEikgMMB1IaHkBEblfRLaIyJaCggKiQgMpq6mnrsHe7l/GYDAYuhLOct/p6enWJykFhYe6hICwwkPAVBHZDkwF8oAGpdQqYDmwAXgd+AJoVaZVKfWSUmqCUmpCbGwsUWE6K9CU2zAYDAY31FdDdQmkXGnp8I4UEHk0f+qPd4w1opQ6rpSap5QaBzzqGCt2vD6llBqrlLoGneZw0NcFo0IDAZNNbTAY2p/u3mtGKQW15RAaA0NnWprTkQJiM5AqIskiEgjcAbzneoCIxIiIcw2PAC87xv0cpiZEZAwwBvAS2KtxCggT6mowGNqT4OBgCgsLu62QUEpRWHCK4MK9cNEd4B9oaV6HOamVUvUishBYCfgBLyul9orIE8AWpdR7wDTgNyKigHXAtx3TA4D1IgJQCnzV4bD2Su9QbWIqMiW/DQZDOxIfH09ubi4FBQWdvZSzJrjqBPFbn4ZvrLQ8p0MT5ZRSy9G+BNexx13eLwWWuplXjY5kahNRYcbEZDAY2p+AgACSk605drskSsGiuyEuFeKGWZ7W2U7qdiXKoUEYE5PBYDC4kLsZTmfAOB9VXlvQowRESIAfgf42o0EYDAaDK9sWQ0AYjJrn+1gXepSAEBFHPSYjIAwGgwGAmnLY+y6MnAtBEW2a2qMEBOhIJmNiMhgMBgebXtLhreltMy9BTxUQJorJYDAY4OgG+PhXMOx6SLikzdN9CggR+Z6I9BLNP0Vkm4hce1aLPQ9EhRkTk8FgMFB2Et76GkQlwU0vgE4baBNWNIh7lVKlwLVAFLAA6JyC6BbQJb+NiclgMFzANNTD0nuguhRufwWCI8/qNFYEhFPszAJeUUrtdRnrckSFBlBcVddtMx4NBsMFTEM9fP4n/fR/Lnz8JBz9HG78M/QdedansSIgtorIKrSAWCkiEUCXLZcaFRpIg11RWu0z8dpgMBi6FrvegNWP69ezpbIINr4AF30Fxtx2Tsuxkkn9dXTjniylVKWjRtI953TVDsS1YF9kSEAnr8ZgMBgsUl8Ln/xWvz914OzPs2sJNNTC5G+d85KsaBCrlVLbXKqsFgJ/OOcrdxDOkt8m1NVgMHQrti+GkmO62uqpfb6Pd4dSsP0V6D8W+o0+5yV5FBAiEiwifYAYEYkSkT6OLYnWjX+6DL2dFV1NqKvBYOgu1FXBumdg0GQYczsUZIC9VQsc35zYASf3nFXOgzu8mZi+CXwfGABspckxXQo83y5X7wCaSn4bAWEwGLoJW16GshNw8z/gzBGor9Kv0SltO8+2V8A/GEbd0i7L8igglFJ/Av4kIt9RSv2lXa52HjAF+wwGQ7eiphzWPweDp0HS5RAQosdP7W+bgKirgt1LYcQcCOndLkuz4oOwi0jj1RzmpnP3fnQQvYIDsIkp+W0wGLoJW/4Jladh+mP6c6yjHPep/Z7nKAXvfw/eXAAluXps33tQU9Lmiq3esCIg7nM6qPW61BngPisnF5EZIpIhIpki8rCb/YkiskZEdonIJyIS77LvdyKyV0T2i8ifRaylAdpsQmSIyaY2GAzdhD1vQ/zFkHCx/hwYprOfvTmqN/wFtv4LMpbDoktg0991xdaoZK2FtBNWBISf681ZRPwAn/3qHMctAmaim//MF5GWTYCeARYrpcYATwC/ccy9FLgM3Wp0FHAxMNXCWgGIjQgiv6Ta6uEGg8HQOZTkwomdulaSK3EjPGsQRz6Dj36hTUnf2aqFy/KH4OhnMO6rZ1VSwxNWBMQK4E0RuUpErgJed4z5YiKQqZTKUkrVAm8Ac1ocMwL42PF+rct+BQSjBVEQugWp5dTC4f17sfd4qdXDDQaDoXM44Gi4OWx28/G44VB4SOdGuFKWD2/dA30Gw43Pa01jwbtw04sweDqk39Wuy7MiIH6Kvnk/6NjWAD+xMG8gkOPyOZfW4bE7AWcHi7lAhIhEK6W+cFzzhGNbqZTyYpBrzqgBkZwoqeZ0eY3VKQaDwXD+yVgG0akQk9p8PG4E2OuhMLNprKFeC4fackd9pV56XATGfgXu+i+Ex7Xr8nwKCKWUXSn1olLqFsf2N6XUWQTouuUhYKqIbEebkPKABhEZAgwH4tFC5UoRuaLlZBG5X0S2iMgW12biowbqwlRGizAYDF2WqmJtLho2q/W+uOH61dUPsf89OLYBZj/XtL+DsVLuO1VElorIPhHJcm4Wzp0HJLh8jneMNaKUOq6UmqeUGgc86hgrRmsTG5VS5UqpcuBDYHLLCyilXlJKTVBKTYiNjW0cHzFAS9Y9eSUWlmkwGAydQOZHWkto6X8ArVXY/Jv7IbYthsgEnUh3nrBiYvp/wItAPTAdWAy8amHeZiBVRJJFJBC4A3jP9QARiRER5xoeAV52vD+G1iz8RSQArV1YNjFFhgSQGB1qBITBYOi6HPgAwuJg4ITW+/wDIXpIk4AoPgZZn8DYO8F2/vq8WblSiFJqDSBKqaNKqV8As33MQSlVDywEVqJv7kuUUntF5AkRudFx2DQgQ0QOAn2BpxzjS4HDwG60n2KnUup9619Lm5n2HDcCwmAwdEHqa+DQR5A2w/MNP254k4lp+3/067g7z8/6HFip5lrjeMo/JCIL0WaicCsnV0otB5a3GHvc5f1StDBoOa8BXerjrBk1IJJlu05QXFnbWJ/JYDAYugRH1kNtGaR5edaOGwF7/ws1ZbDjPzrTuveg87VCwJoG8T0gFPguMB74KtC+sVQdwKiB2g9hHNUGg6HLcWA5BITBYC/pXXHDAQVf/g1KctqtAF9bsCIgkhzO4lyl1D1KqZuB8yvGzoJRA3Qkk/FDGAyGLsWZI9r/MOTKprpL7ohz5BV/9gcIiXLvzO5grAiIRyyOdSmiwgIZ2DuE3UZAGAyGroC9Ab5YBC9MhtoKuORB78dHJenKrLXlOnLJP+i8LNMVjz4IEZmJbjM6UET+7LKrFzqiqcszaqDJqDYYDF2AklxYcjfkbYHU6+D65yAy3vscmx/EpulSHOO+en7W2QJvTurjwBbgRnQ/CCdlwA86clHtxeiBkazce5Ky6joigk37UYPB0Emsf0438rn5nzDqZuv1kobOgIj+7dId7mzw1g9iJ7BTRF5TSnXL5gojXTKqJw2O7uTVGAyGCxKlIONDGHI1jG5jI5/p/9cxa7KIFR/ERBFZLSIHHVnU2RYzqTsd46g2GAydzvHtUHa8U5zM54qVPIh/ok1KW4H2qsF0XoiNCKJvryDjhzAYDJ1HxnIQPxh6XWevpM1Y0SBKlFIfKqVOKaUKnVuHr6ydGD0w0kQyGQyG5lQWwe9TdSJaR3NgOQyaDKF9Ov5a7YwVAbFWRH4vIpNFJN25dfjK2omRAyI5XFDO5iNFnb0Ug8HQVTi0CipOwc7XO/Y6Rdlwaq/7iq3dACsmpkscr64VpRRwZfsvp/2ZO24gS7fmcutfv+DOSwbx05nD6GUimgyGC5sDy/Tr4bVQUw5BlqoHtZ0MR6WhtO4pIKz0g5juZusWwgEgKSaMVT+Ywr2XJfP6pmNc89ynbD16ptVxdrvihU8yOXK6ohNWaTAYzht11ZC5RmcqN9TA4Y99zzlbDiyHuJHQJ7njrtGBWOkHESkizzkb84jIsyISeT4W116EBfnz+A0jePdbl+Fvs/Hw27tosKtmx3y4J5/frchg6dbcTlqlwWA4Z0pPwGd/1M14PJH9KdRVwFU/1yUsMpZ7PvZcqCzSDX66qXkJrPkgXkYnx93m2ErRPSK6HRcl9ObR2cM5dKqc/25v6l1U32Dn2dUZAGQXGg3CYOi2bP4HfPRzWHQJ7PfQIeDAMgiMgJTpOqv54ArdzrO9ObgSlL3bmpfAmoBIUUr9XCmV5dh+CQzu6IV1FDNG9mPUwF784aOD1NbbAXhnex5ZBRX0CvbnqBEQBkP35cRO3XUtLBbe/Cq8uQAqTjftt9u1QBhyla5tNGwWVJ2BY1+0/1oOfAARA2DAuPY/93nCioCoEpHLnR9E5DKgquOW1LHYbMJD16aRe6aKNzcfo6a+gT99dIiL4iO5adxAjpyuRCnl+0QGg6Hrkb8Lkq6A+9fClT/TT/Gvz4f6Wr0/byuUn4Rhjj4MKVeBX1D7m5lyNmlBNPwG62U1uiBWBMSDwCIROSIiR4DngQesnFxEZohIhohkisjDbvYnisgaEdklIp+ISLxjfLqI7HDZqkXkpjZ8L69MHRrLxOQ+/PnjTF7+7Ah5xVX8+LphJMeEUV5TT2FFbXtdymAwnC/K8vXNv/8Y8AuAKQ/B3L9C7iZY7ehTlrFMJ62lXqM/B4XrRjwHlumSGFaoKfe+v+K0LswXGQ/Tu3zha69YiWLaoZS6CBgDjFFKjXPUafKKiPgBi4CZwAhgvoiMaHHYM8BipdQY4AngN45rrlVKjVVKjUWH01YCq9rwvXytjR9fl0ZBWQ2/XXGAyYOjuWxINEkxYQAmkslg6I6c2KVf+1/UNDZqHkz6Fnz5IuxeqgVB0uXaOe1k2CwoPtrU3tMbB1fBbxPhyGfu99sbYOm9UFUEt73S/DrdECtRTL8Wkd5KqVKlVKmIRInIryyceyKQ6fBb1AJvAHNaHDMCcMaYrXWzH+AW4EOlVKWFa1rm4qQ+TE+LBeCh69IQEZKitYDINgLCYOh+5DueW/uOaj5+zROQMAn+9204fbDJvORk6ExAmnIjPGG3w5pfgr0e1jzhXuNY+5SOkpr9rNZkujlWTEwzlVKNMWNKqTPoPhG+GAjkuHzOdYy5shOY53g/F4gQkZZlV+8A3KY7isj9zvDbgoICC0tqzm/mjeHFO9MZn6ilfHxUCH424Whhu8oig8FwPjixE/oMhuBezcf9AuDW/wdBEfpz2szm+yP6QvwE336Iff/VJbsHT4ecLyHzo+b7Mz6E9c9C+l2d1r+hvbEiIPxEpLGVkYiEAO3V2ughYKqIbAemAnm4FAQUkf7AaGClu8lKqZeUUhOUUhNiY2PbfPF+kcHMHN2/8XOAn434qBAT6mowdEdO7GpuXnKl1wD4yhKY8TT0dtMxOXkK5O/WSXTuaKiHtb+G2GEw/w19jo9/1aRFFGXDO9+EfmNg5u/b5/t0AawIiP8Aa0Tk6yLydWA18G8L8/KABJfP8Y6xRpRSx5VS85RS44BHHWOuGS63Ae+ez34USdFhJtTVYOjKlOTBiv+DOpdgyqoz2o/Qz4tZZ2A6TPLQ5rPfGG068uSH2P0WFB6C6Y9CQDBMewRO7NBmqboqWLIABLj9Fb2/h2DFSf1b4FfAcMf2pFLqdxbOvRlIFZFkEQlEm4recz1ARGJExLmGR9BJea7Mx4N5qaNIjgkzoa4GQ1fmo5/DxkWw552msfzd+tWTBuEL57z8Xa33NdTBJ7/RQmT4DXps9G0Qnap9Dsse0tef93fdR7oHYUWDQCm1Qin1kGNza+5xM6ceWIg2D+0Hliil9orIEyJyo+OwaUCGiBwE+gJPOeeLSBJaA/nU4ndpFxKjQ02oq8HQVTm5T0cjAWx/pWn8hMNBfbYCIioJgiKbzuPK9le1dnLlz5pyGvz8dQjrqX2w41WY8uNu2e/BF1aquZ41SqnlwPIWY4+7vF8KLPUw9witndodjmuoa0x4kHMtrNp3kqlDYwkO8DvfSzIYDE4++bV2Nk+4Fz7/I5w+BDGp2v/QayCExZzdeUV01JE7AbH1X9B/bFPuhJMRcyHhJR3KOq175zt4wpIGcSHhLtR1y9EzfPOVrc3qNxkMhvPM8e26vtLkb+vcBvFr0iJO7PTuf7BCvzFwcm/zukxVZ/S502a1zoi22eCe5TD/dbD1zAdHSwJCREJEJK2jF9MVcBfqunJPPgCHTvnIoDQYDB3H2l/rp/VJD+rQ1KEzYMfrUF2iHchna15y0v8iqK/WuRJOjnwOKB3l5A6bX7cupeELK4lyNwA7gBWOz2NF5D3vs7ovLUNdlVKs3KcFxOECIyAMhk7h2Je6C9xl34NgR7eBcV/VXeHWP6erpp5rYppzvqujOnsdBITCwPHndu5uihUN4hforOhi0KU3gO7Z/cIirqGuB/LLyCmqIjjARqbRIAyGzuHTpyEsDibe3zSWei2E94WNL+jP56pBRKeCf0hzP0T2Ot1P2j/w3M7dTbEiIOqUUiUtxnp0DKhrqOvKvfmIwK3jE8grrqKqtsH3CQwGQ/tRlK27vk28DwLDmsb9/OGi+dBQCyF9tJP6XPDzh74jm2o6lZ+Cgv2ezUsXAFYExF4R+Qo6ozpVRP4CbOjgdXUqrqGuq/aeZPygKC4Z3AelTJ0mg+G8s+M/IDYYe2frfeMW6Nf+F7WPL6D/RdrEZLdr7QGMgPDBd4CRQA3wGlACfL8jF9XZOENdPzt0mn0nSrl2ZF9SYnVTc+OHMBjOI/YG2PGa7tsQ6UZDiBkCVzwEF3+jfa7XfwzUlELxES0ggiLP3XTVjbGSBzFMKfUojlIYFwLOUNeX1mUBcO2IfvSLDEYE44cwGM4nhz+G0jy47teej7nqZ+13PacwOLFTC4iky3tsCKsVrGgQz4rIfhF5UkRG+T68++MMdd13opS0vhEkxYQRHOBHQlSo0SAMhvPJtsUQGn3++jrHjQCbv66xdCYbkq84P9ftolipxTQdmA4UAH8Tkd0i8liHr6wTCfCzkRAVAsB1I/s2jqfEhnG4wPggDIY2UXi4qeVnW6g4rUtoj7nj/EUR+QdB7HDY+67+fAH7H8B6LaZ8pdSf0a1GdwCP+5jS7Ul0mJmuHdmvcSwlNpysgnLsdmtBXNV1DXyw67gp/Ge4cCkvgBcmNYWitoWdb4C9DtIXtP+6vNHfUdk1NEYLiwsYK4lyw0XkFyKyG3BGMMV3+Mo6mUmDoxk9MJKRA5qaj6TEhVNTbyevuMrLzCY+2HWCha9t5/PMwo5apsHQtTmyXoehtmyu4wuldBmNgRMg7jzfpJ1+iOQrdDmNCxgr3/5ldJLcdUqpaUqpF5VSpzp4XZ3Og9NSeP87lyMuoXND4nQkU2YLP4Qnv8ShU2UArNyb30GrNBi6OM5Q0ZxNzfs3AJQeh98maUd0S45vg4ID5197AF2YDy548xJY80FMVkr9USl1/HwsqCvTGOrqEsm0NuMUVz37KRuzWmsJh09pf8XqfSctm6UMhh5F9qe6NEZDjRYSrux/XxfD2/lG63n7P9DF+Ibf2HpfR5MwEW552X3exQWGRwEhIkscr7tFZJfLtltE3HTV6Pn0CQskKjSgmaP6PxuPAbArt7jV8VkF5YQE+JFfWs2uvJbJ6AZDD6c4B4qyYNK39c3eqU04ObBMvx5coZvyuJKxHJIug9A+52etrojAqJu1w/oCx5sG8T3H6/XADS6b87NPRGSGiGSISKaIPOxmf6KIrHEInk9EJN5l3yARWeUIsd3naCDU6aTEhjealE6WVrM2Q1vbMvKbm5lq6+0cLarklvHx+NmEVcbMZLjQOLJevw6brdt9ugqIqjNw5DMdVlpdAkddijMUHtbmpbTZ53e9hlZ4FBBKqROOt99SSh113YBv+TqxiPgBi4CZwAhgvoiMaHHYM8BipdQY4AngNy77FgO/V0oNRxcL7BJ+j5TY8EYT09KtuTTYFYNjwsg4WdrsuKOFFTTYFeMTo7gkuY/xQxguPLLX6RyGuBHanp+3FWq0X45Dq0E1wIynwT9YawxOnJrFsPOU+2DwiBUn9TVuxmZamDcRyFRKZSmlaoE3gDktjhkBOD1Ua537HYLEXym1GkApVa6UqqQLMCQunMKKWooqanlj8zEmD47mymFxHDpZToOLn8GpZaTEhnPdyH4cLqgwWdiGCwelHJnIjkig5ClaIBz9Qu8/sAzC++n9g6fBgeV6Dmhh0Xc09B7UWas3OPDmg3jQEdqa1sIHkQ1Y8UEMBHJcPufSuoXoTmCe4/1cIEJEooGhQLGIvCMi20Xk9w6NpOUa7xeRLSKypaCgwMKSzp2UOJ0f8erGo+QUVXHHxATS+kVQU29vLBEONPopBseGcc0InWy3ap/RIgwXCEVZukSGMxIo4RLwC9RO6/oaHfaaNkMLj7RZUHIMTu7RyXE5X2qzlKHT8aZBvIb2NbxHcx/EeKXUV9vp+g8BU0VkOzAVyAMa0DWirnDsvxgYDHyt5WSl1EtKqQlKqQmxsbHttCTvOCOZ/vbpYSJDArhuZD/S+kUAkJFf1njc4VPl9I8MJizInwG9QxgTH8mqvSfPyxoNhk4n+1P9mjxVvwaEaCGRvQ6y10NteZOPIW0mIFqLOLhCN/8x5qUugTcfRIlS6ohSar7D71CF7gMRLiJWdL88IMHlc7xjzPUax5VS85RS43AUA1RKFaO1jR0O81Q98F8gvS1frKOIjwol0M9GRW0Dc8cNJDjAj9S4CER0cyEnmQXljXkTANeO6MuOnGJOllZ3xrINhvNL9jqIGADRKU1jyVMgf7dOgAsIa9IuwuMg/mI48IE2PUUmnHt/aUO7YKnlqIgcArKBT4EjwIcWzr0ZSBWRZBEJBO5AayOu544REecaHkEn5Tnn9hYRp1pwJbDPwjU7HD+bkOwoB37HRC3/QgL9SIoO4+BJLSCUUhw+Vd6obQBc5yjZsWqf0SIMPRy7XWsJyVOa92hIngIo2PdfSL0aAoKb9g2brfswZH6kNYoe3Oe5O2HFSf0rYBJwUCmVDFwFbPQ1yQ4UqgcAACAASURBVPHkvxBYCewHliil9orIEyLizH6ZBmSIyEGgL/CUY24D2ry0xuEHEeDvbfliHckVqTFMT4tlWL+mMhxD+4Y3mphOltZQUdtASmxT96shceEM6hPKuoPnx1diMHQaBfuh8nTrTOQB6VpzgNYhrE6fQ0Ot8T90Iaz0g6hTShWKiE1EbEqptSLyRysnV0otB5a3GHvc5f1SYKmHuauBLqlnPnZ9y2hdSOvXi9X7TlJd19AsgsmJiDAhKYp1BwtQSjUr4WEw9BhqK+DzP+n3LUtl+wdC4mQ4vBZSWwRHxqTqntAVpyDxsvOzVoNPrAiIYhEJB9YB/xGRU4Cped2CYf0isCvdUMgZzurqgwAYnxjFO9vyOFZU2Vgt1mDoMWR9Au99F4qPwuSF7sNUp/0fjLjJfYb0rN/rpDm/gA5fqsEaVgTEHKAa+AFwJxCJTmozuOCMZDqQX8bhgnIigvyJjWieqp8+KAqAbcfOGAFh6FmsfBS+eB76pMDXlulObO6IH683d6RM77j1Gc4KnwJCKeWqLfy7A9fSrUnsE0qgv42M/FIOF5QzOC68lRlpaN8IwoP82Xr0DHPH9fiK6YYLhcoiLRxG3QJzntchrYYegZUopnkickhESkSkVETKRKTU17wLDX8/G6lx4WScLOfwqYpmDmonfjZhbEJvth1tXdjPYOi25G3Vr+PvNsKhh2Eliul3wI1KqUilVC+lVIRSqpfPWRcgaf0i2JlTTH5pdTMHtSvpiVEcyC+lvKb+PK/OYOggcjeD2HSUkqFHYUVAnFRK7e/wlXQFygvA3nDW09P6RlBSpcsWt3RQO0kf1Bu7gl05Rosw9BByNkHcSAhy/5s3dF+sCIgtIvKmiMx3mJvmicg839O6GdUl8KeLYNPZp1s4HdWARw1inMNRvfXombO+jsHQZbDbtYkpfkJnr8TQAVgREL2ASuBamveE6FnkboG6Cji08qxP4Uyc87cJidGhbo+JDAkgNS6cbceMgDD0AE5nQE2p7sJm6HFYiWK653wspNPJ3aJfj34B9bU6qccKpzPh1bkwZxF9k64gMiSA6PBAAvw8y97xiVF8uCcfu11hs+lIpxV78rkoIZL+kcbJZ+hG5G7Wr/EXd+46DB2CFQ3iwiB3k3a01VdB3hbr8/YsheJjsPRepOwEVw/vy7ShcV6npA+KoqSqjqzTOoL4rS05PPDqVr7/xg6UMr2rDd2InE0Q3Buih3T2SgwdQM8XEKcPwQc/gGovkbl2u9Yght+ghUTL3rneOLBM/+eoq4Ild/PsvOE8fkPrUhyupCc6EuaOnmFPXgmP/ncP0WGBfJldxGeZp61f29AxKAVrfw0HV3X2SjqXklx49wGo8hJQkbtFaw+mdEyPxFvDoO85Xrt3YZRPfwtbXob/faupY1VLCjOhuhiGXAP9L7IuIIqP6QqU6XfBjX/RWsjqx31OGxwTRmRIAGsOnOSbr2wlJiyQD757OQN7h/D7lRlGi+hsPv+T/t2s+Knn38yFwMYXYefrsOM19/urS3TvaON/6LF480HcA/wJ+AtdpBdDm6k6A/ve0zVh9r+vsz0v/U7r43I36deEiVB4CL54AWorIdC9o7mRDEfV87TZEDNEq9tfvgj9RsE4zz2VbDYhfVBvVu49SaCfjbcemEz/yBC+d3UqP1m6i5V7TzJjVD+3c+0n96NKj+OXepWVv4ChrWSvhzW/1L+Zoiw4+rnnshG+yN8D+9/zvD/p8tYVT7sK9bVaOIDu3zDpwdZaQt5WQBn/Qw/Gm4lpv6MPRMuWo7tFxErL0c5n91JoqIHbFsPwG2H1z+HI562Py90MwZG6mmTyFLDXQY7PiubavBQzVAsHgGuegMTL4X/fhrfu0XkVHpiQpIuVPTFnJBcl9AZg3riBDI4N49lVGc36W7uS+cZPqPvPfG3SMrQvpSdg6T3aZPiNNRDUC7a9cvbnW/Gw1kQ8be8+0HU1lIMfQmUhDLseTu2DvG2tj8nZDAgM9FBbydDt8dZRbj667WcmzVuOXu947fpsWwz9RsOAcTBnEUQl6RtAWYve0DmbYeAE3R930GSw+fs2M1Wd0U+XaS6tEf0DYcG7MP0x3R1r0cWwa4nb6XdNTmTxvRO5Y2JTxUt/Pxs/uiaNQ6fK+d+OPLfz+pTsJ5gaag9+bOUvYPBGSR5krmna3vqa1hxve0V3ORt1M+z7nzaluHLmiO6r7I3KIji6Aa74EfyipPU25wXds/nEjo76dufGtld0R7gb/wL+IbB9cetjcjdB3HAINoUVeipendRKqXyl1EXACSDCsR13tCD1iYjMEJEMEckUkYfd7E8UkTUOzeQTEYl32dcgIjscmxc93QMndmr/wLi79OfgXnD7q9pZvcJlKTVl+gnJqSYHhun3vgTEodVgr9dPWK74B8LUH8MDn2nt4p374Pj2VtMjggOYMrR1H+2Zo/oxckAv/rTmUGtfREUhMXatlVTubvufxOBCfS38vxnw6rymLWcj3PhniBumj0lfoKPadru0LMnZBH8ZD6/d5j3r/tAqUA2tG+M4GTpDB0QcWO5+f2dSkgeH18DYr+iy3CPnwu63da8HJ87ADpMg16OxUqxvKnAIWAS8ABwUEZ+GUxHxc8yZCYwA5otIy/CeZ4DFSqkx6BLiv3HZV6WUGuvYbqStbHsF/IJgzK1NY31HwORvw953dW9caLKjJrjYUZOn6Ju665Oj3d78/AeWQXhfz+p1bBrMfwPET/s/LGKzCXPHDeRoYSWFFbXN9lXlaDX/uOpDSPbqcyoLcsGzfbEOMrj+j3DvKr0t3Aqjb2k6ZkC6LiGx3WFmKi+AJXfrh4isT+CTpz2f/8AyCO+ntVd3hEVrbTWjCwqIna+Bsjf50dIXQG2Z1qacOAM74o2DuidjJcz1OeBapdRUpdQU4DrgDxbmTQQylVJZSqla4A10bwlXRgBOW8laN/vPjroq2L1Eh62GRDXfd+lCCIrUYYzQlOgz0OVJKHmK/g9y5HNoqIP1z8Jv4rVvoeqMNi9kfqSfAm1e/oShfSDx0jY/JTrrOB12NB5yUpqlq2a+UD+HoJrCpuQ+Q9uoq4J1z+gb9PivwaBL9BbTIpZfRN8cj2/XGunb90JVEdz9AYz9Kqz7HRx0k3lfV61NVmkzvf8+0mbByT3aZNVVsNth+6uQdAX0SdZjgyZrv4zTH1NRqB35YCKYejhWBESAUirD+UEpdRCw0vJpIJDj8jnXMebKTsBZ12kuECEi0Y7PwSKyRUQ2ishN7i4gIvc7jtlSUODiEN7/gX76T1/QelJIlI5kylgOuVu1/yEmDUJ6Nx0TfzH4B8O2f8Pfp8OaJ3Rk0o7XYdEl8NEvoba8tXnJHcNm6x69RVm+j3XQKCAKmjfuazi+kxx7LO/bL6Mef+3nMLSdLS9D2Qm48jHf8ftjbge/QHjtdm12nP0s9B8Ds5+BvqPhnfvhTAuLa/Y6XbbF1+9jmMN/1ZXMTEc/0wIr/a6mMRGtTRzboEOAF12sBeOVj2lN2dBjsVqs7x8iMs2x/R1or0fXh4CpIrIdmArkAU67SaJSagLwFeCPIpLScrJS6iWl1ASl1ITYWBd7/rZ/Q+9ESPJgCZv0AIRGw8dPag0ioUWYnn8QDJoEB1dA+Sntu/j6KrjvY21W2rhIN1+3EqKY1vabwIDIEIIDbI29rZ2EFu5hr0oiddBAdvmP7prmia5OTTmsfw4GT7MWvhraRwv5shP6puk0uwSEwO2LdRTSkgVaa3By4AMIjGjdk7klfQZD7HDv/47HNmrhVHueuvxueklr2MNbxKFc9BVtLl39uA72+OY6mPLj87MmQ6dhRUA8COwDvuvY9jnGfJEHJLh8jneMNaKUOq6UmqeUGgc86hgrdrzmOV6zgE8AD8bcFhzdAEfWw8Xf8KzeB0XA5T+ArLXaZOAujvuKh+DyH8K3v2z6zzJgrBYSM56GGb+BgGDf64lKhL6j2nQzt9mEwTHhjb2tAagupXfVMbIDhjAmvjfLasdpO3DBQcvnNQCb/gaVp3WkmVWmPwaXfQ9m/r75eJ/BMO9v2vz04U/0mN2uHyyGXKUfNHwxbLb+zVYWtd5XegLeXKDPd8xC2PW5svXf2l926cLWjX8i+urf/axn4OurtT/P0OOxUqyvBu2HeK6N594MpIpIMlow3IHWBhoRkRigSCllBx4BXnaMRwGVSqkaxzGXoRsX+VosfPwr/ZR/8Te8Hzvh67DheSjPd+9oS77C/ROgX4BOGmoLabNg/TPadhsW7ft4ICUunO2uFV9P7gGguNcwUuLCeL52HD8LBjKWQezQtq3nQqXqjDaRDJ3RWmv0RswQnePijrSZOpR1/bPaHh+TBuUn9Y3fCsMcv42DK2Hs/Kbxhjodkl1bDoj2Nw3xkhxZUwYb/6p/z076jdY+Fisc3w7LfwwpV+rv445L7rd2LkOPocNqMSml6oGFwEpgP7BEKbVXRJ4QEWdU0jQgQ0QOAn2Bpxzjw9GmrZ1o5/XTSql9Pi+a9YnOTbjiId9Z0IGhcN1TOrGto+2ow2Zrp/fBFZanDIkNJ6+4iqpah8XthM5NrIsbTUpsOCeIpqzPKB0tY/CN3Q7vPqhNNVe2QXuwwvRHtblx2Y/gs+e0KSb1Gmtz+4+DiP5a0Lvy0S/g2Bdww58hbkRTtr87Dq2GRZNg7VM6Qm/vu7o8xrIf+c7XAK29LLkLwmJh3j/A5mdt7YYej08N4lxQSi0HlrcYe9zl/VJgqZt5G4DRbb7gx7+CXvG6N64VRt/SPKyxo+h/kV5XxnIYd6elKSlxYSgF2acrGDGgF3V52ylWkUT3T2xsRnQoagrph1/UiX8R7ktzGBx8/kedHTzjt/rJuj2x+cHNL8PfrtD/xslTWkfPeZxr01rIzjcdSZUCxUd1WZiJ9+sw7SPrdYip3d7cbFpdCssfgl1vas3l66uaoop2L4W3v67NkH1Her5+Qx28+01tzrp3hWUN13Bh0HOquVaX6jLdU39izfZ7PhHRN4HMNTpT1wJOIeB0VDfk7WSPPYmk6DBiwgPpFezPOv9L9blfvg6yPu2w5Xd7sj7VAQkj58El3+yYa4THwq3/1hFPo25u29yR83TU0zv3wTvf0GtNmATXOhTqhIk656Aws/m8z56D3W/BlJ/AA+ubh5zGDdevp7x0Cz6+Q0fpHVqlfWom6c3QAp8ahIgMBX4MJLoer5S6sgPX1XbKjkPUEJ392RUZNgs2/12bwYbNar2/vgZsAY1PiMkxYYigHdV11QSeOchedT3TokMREYbEhbOxzAZ3vQfvfxcW3wjjFsC1v2oesnuhU3oclt6r62zd+JeOLUs96BJ46JCu69UWkq+AH+xtHgkVlQR+jv9uziCK3E3N/U37P9DaypWPtj5ndKouGXPKjWW2rkon+W34C4TF6NIiI9qei2ro+VjRIN4CtgGPoQWFc+ta1FXBtEe0E7krkni5Dh88+GHrfQ11uh/2xkWNQ8EBfiREhWoN4tQ+bKqBPfZkkmLCAK1hHC6o0DeXBzfoKJsdr+mQS0MTqx/Xv43bX4Eg933C25WQ3mcnhCLjtTPcufm5PLtFp2qh40zqBN3npPCQ51Ie/oE6uc2dBrHiYW1yGztfR+kZ4WDwgBUBUa+UelEptUkptdW5dfjK2kpk/PnxJ5wt/oH6Zu6uxlPeNh1nn/lRs+GU2DAtBE7sBCA/NJXwIH3jSIkLp6CshpKqOh2SeM0T2ryWvV7nbhg0x76EtBndO6HLZtOZ/jkuAsIZnJA20/O8uOHuNYjDH+vQ7TmLrPtKDBckVgTE+yLyLRHpLyJ9nFuHr6ythMV2/eiL5Ck6S9Vd5i3ozG6X+kopseFkFZSjTuyiQsIIiE5utg9onkw37HpANfWpuNCpLIKSY9BvTGev5NxJmKhv9jVl+vOBZTr4oXeC5zlxI/TvzTXJruykrkGVMKlDl2voGVgREHejTUobgK2OzRQBOhucmdcttYhsh4O5tgwKGquakBIXTk29ndrc7ewniaSYJhOJ23pNfUfqRjcm9FWT72hb0v+izl1HexB/MaB0ccnyU9rc5Mm85MTpqC440DTmNFOZJj8GC/gUEEqpZDfb4POxuB5H7DCt6bgKiLoqXULaWZLDJd59SFw4fjTgX7CP7XWJjf4HgISoEAL8pHm9JhF908j6RJeUsMrGv8LRL87yS3VhTvQgAeGsGpyz2aEhKvfBDq7EObKdXf0QuZt0MERP+JsYOhwr5b4DROS7IrLUsS0UkS7qCe7iiGgtIntdUyexnC9117vxX4OQPs0ckSmx4YyXg/jZa9hmTyUpuklA+PvZSIoOa1WviWGz9fkOW2woVFUMKx/RSVY9jRM7ITJB11Pq7oT01g8YuZu1hhg5SJdw8UZUkm724yogcjZr4WClTIzhgseKielFYDy6F8QLjvcvduSiejTJU3Q5hNOH9OfsdTrzNvFSrfa7OCL7hAVyQ9B2avFnnX0MidHNs8N1JFMLATFosnY8Wq39dHSDzvL2VA/IldLjcOQza+c9n5Tkamd0S/J39Qz/g5P4CbqpkTNU2le0lM1PO+edjuqGOl1Sw5iXDBaxIiAuVkrdrZT62LHdA5hf2NnS6Idw+B2y12nzQVCErg90OkPXDAJQiqv9trKhYSQVhDQzMYE2QR0trKS23qWZkZ8/pF6ny3o01Ptej9PcpRp0wpQnlIK374NXb7ZWvuF88uFPYfGc5kmINeVaCPckU0r8RF3GvqHGeq2nuBFNGsTJPbpDXlvqUBkuaKwIiAbXUtsiMpimktyGthKVrM0e2esc2d/bmoSG88kuzxFFXHCA/g0nWG0fT0x4UGOIq5OUuDAa7IpjRS1KQQ+bpYXMMQt+hex1+voR/b07t7M/1b0C6qubbPtdgdpKnaFeX6Wr8zo5uRdQundDT8H5+wjuDYMutTYnbrgOoa4satJOTRc4g0WsCIgfA2sdPaM/RXeA81Du0eATpx/iyHpdWFA1NAmIgeN1n2Lnf2RHQ6DVDeNJjmldfDCtr24W/9Sy/ZwoqWrakXKVbrfqy8xUXgCn9ureCM5SIK7ZvE6cVXJDY/Rn14StzibrEy0coHnPDUfuSI/SIGKH6T4mw2Y3T6TzhtNRXXBA/7uF99M5QwaDBaxEMa0BUtG9IL4DpCml1nqfZfBK8hT9hL/xBX0jd9bQCYpwVO50CojllPQZwymiSIwOa3WaEQN68djs4XyRVcg1z61j8RdHsNuVzhgePE1rBE5nuDuOrHesZ6q+6dRVNJm+XDm0Sq/pqp9p7cdbZdHzTcYynaE+Yo42qznzSPJ3aoEW0b9z19ee2GzwjY903SSrNNZk2qf/3RIu7thyI4YehUcBISJXOl7nAbOBIY5ttmPMcLYkOfpMZK/TwsG1OUv8BF37vyQPjm/DPlSHMibHtBYQAN+4YjCrvj+VsQm9efx/e3loqePJedgsXRX0V32btg9/2nxy9jrd+az/WL2mwIjWZia7XWsPUUkw9k5t5ugqvbDtDZCxQpfWHjFHNwLKcQivEzu19tDTboZ9Bret1lOvAVqAZq/TSXPGQW1oA940iKmO1xvcbBaaMRs8EjlQ18kB/fTuSvxEqCmBDX8GoPe4OTw1dxS3TvBsFhgUHcorX5/I/ImDeH/ncUqr62D0rbo21aQHdAXTQZN1O0nXDnTZ6yDpMm2u8A+C1Kt1jL3dxel94H0dDeSscxV/MZTk6PLQnU3OJi0Uhs2CIdfo+P6MZdqJfupAz/I/nC0iWotwmt+M/8HQBjwKCKXUzx1vn1BK3eO6AU9aObmIzBCRDBHJFJGH3exPFJE1IrLL4eOIb7G/l4jkisjzbflS3QKn36FlX2unuWnzPyEqGYkbzp2XJBIX4T1uXUS4OX0gdQ2KtQdOQWAYTHtY12i65gm4+R8QEAqfOMwTJblQdLj59dNmQ8UpXTYd9BPnR7+EmKFa4LiuryuYmTKWaaEw5BoI7qVrXR1YpqN27HU+/Q//7/NsnluV0dSUqacSN1z/PWz+um1uNyenqJIfLtlBRY2FKD3DOWHFSf22m7FWTX5aIiJ+wCJgJjACmC8iLRvZPgMsVkqNAZ4AWhpXnwTcVLfrAYz/Goz9KgxMbz7eJ0VHqdjrtF+gDSaScYOiiAkPYtW+k613hsXoVql734H8PbqoHzQXEKnX6JvI/vfhixfghcm6febM3zXVueo3Wvc86GxHtVL6qTh5ihYOoP9eRVm6RwJ4zYEorqzlNx8e4M8fZ3LtHz9l/aGC87DoTsLpqO43unWv6W7I29tyeWdbHqvd/c4N7Yo3H8QwEbkZiBSReS7b1wAraZgTgUylVJZSqhZ4A5jT4pgR6Kgo0K1FG/eLyHh0G1IvwfndmP4XwU2LWpcnt9ma7MRpPkoptMDPJlwzIo5PDpyipt7NU/Hkhdp+vfbX2rwUGg1xLt3GQnpD4mXavLXyEUi6HL61EVKmNx3jH6R9FjmdLCBOH9QakGu5Ceffa/M/IaiXDin2wLvb86itt/PkTaMIsNlY8M9N/GjJTuoa7B7nuCOroJx3t+e63bcnr4QPdh1v0/k6BKejuoeYlzZkFgKwal++jyMN54o3DSIN7WvoTXP/Qzpwn4VzDwRyXD7nOsZc2Qk4Hd5zgQgRiRYRG/As8JC3C4jI/SKyRUS2FBT0oCfAETdqAZJwSZunXjuyHxW1DY3/iZoR0hsu/Y42zex/XzumbS1+AhPu0QX/5v0DvrLEfbXQhIlwYgfU17Z5fe1GY7lrFwHRawAMGKfDXvuNbv3dHCileGNTDhfFR7JgUiLLv3cF908ZzNvbclm51/pNRynFD5fs5IdLdlJd11ogP/9xJgtf286KPZ18I+s/RgtLq8l1XZjK2nq255zB3yZ8klHg9u9uaD+8+SD+5/A3XN/CB/FdR8/o9uAhYKqIbEc7xfPQSXjfApYrpdw/mjWt8SWl1ASl1ITY2Nh2WlIXIP0u+OY667HuLlyaEk14kL/np6tLHtCaQ21Zo3np4MkyzlQ4bvYj58L3d+teyJ7MW/ETdMLcyd1tXh9K6bIYmWuaNl8lPtxxYBkMSNdCwRVnhVMv/oftOcVknCzjjomDAN2c6aczhjGwdwhvbMrxOK8lH+0/xY6cYpTSdvGWHCnUCYwPvbVTdwZ03Xe6wu2cDiE4Er63AwZP9X1sB7Art7jd/Dybsouoa1B87dIkKmsb+DzzdLuc1+AeKz6IB0SksYeliESJyMsW5uUBro+f8Y6xRpRSx5VS85RS44BHHWPFwGRgoYgcQfsp7hKRpy1c84InyN+PaWmxrN53kga7mxyIoAi44iGdkJcyHaUUt//tC57+8EDrYz3hNFWcTbjrx0/Cy9fCq/Oatve+07ZzVBZpR7q7ZjkjbtTfzYv29eamHEID/bjhoibh4mcTbpuQwGeZpy3duO12xTMrMwgL1L6Z7NMVrfYfKaxg9uj+BPrbeODVrZTX1FNd18DvVhzg6uc+5cH/dL2+W+1NQVkNc1/YwL+/ONIu59twuJBAPxvfuSqViCB/Vu01foiOxIqAGOO4aQOglDoDjLMwbzOQKiLJIhII3AG853qAiMQ4zEkAjwAvO65xp1JqkFIqCa1lLFZKtYqCMrjn2pH9OF1ey/ZjZ9wfMOlB+O4O6DOYvOIqzlTWsflIG57iIwdCxICmnAOrZHwI65+Fi+bDvav0Nvo23UmvtsL3fCdOB3mim3ITsWnwvV06L8IN5TX1vL/rODeMGdCqdMmtE+KxCby52bcW8f6u42ScLOORWdq+f7SwuVA5VVZDdZ2dSSnRPD9/HFkF5TzwylZm/HEdL3xymAG9Q9h3vJSy6joLX7j7su3YGRrsij15Je1yvs8zTzNuUG8iQwKYNiyOj/Z7eBAytAtWBIRNRBr7Ejq6yfm0fSil6oGFwEpgP7BEKbVXRJ4QEWcT3GlAhogcRDuke2DN6fPPtLRYAvzEfTQTaNNRVCIAGfm6Q1nW6QqKKtrgU0i4uCnUNX8P/ONqXTDPE0XZ8M43tenn+j/CoEv0Nu5Oba463Ibk/NzNugLuAA/PKb0TPJrH3ttxnMraBm6f2Nq3MqB3CFOHxvLW1hzqvTir6xrsPLf6IMP79+IrEwfROzSA7MLmAs6pUSRHh3HpkBh+OmMYn2Wexq7g1a9fwlNzR2FXsDOnfW6cXZVtR/VDivN3di6cqahl34lSLhuiS75cN7IvhRW1bD3q4UHIcM5YERDPAl+IyJMi8it0Z7nfWTm5Umq5UmqoUipFKfWUY+xxpdR7jvdLlVKpjmO+oZRqVSZUKfUvpdRC61/J0Cs4gEtTYli5Nx/lrdQGkHGy6T+uR43DHfETdevKDx+Gl6bqAoNZn7Rupwq6KdKSBfqmfdvi5r0IEi/TNnKr5clBay59R+pcjzby5uZjpPWNYFxCb7f775g4iJOlNXyS4TnoYcmWHI4WVvLj64ZiswlJ0WEcbSEgnP4HZ4n2+6cM5o37J7Hy+1O4PDWGsQm9EdFP2D0Z5/fLPl3hPrKuDXyRVYhScNmQaACmDo0l0M/GqjYEFhjahpVaTIuBm4GTQD4wTyn1SkcvzHBuXDuyL0cLK/nrp1len4Yz8suICQ/C3yY+n8Re+/IYX/3Hl1roOENxv3xRJ9F9fbXjhG5u9Kt+Bvm7Yd7fdckOV/wCIPXa5nWUvGFv0BVwE9oesrk7t4SduSXcfnEC4kHDuHJYHDHhQbzhMDOVVdfx+P/2MOrnKxn+sxUM/9kKfvbfPYxPjGJ6WhwASdGhHDnd3MR0pLCCQD8bA3rrvAMRYdLgaEIcPouI4ACGxkX06Kff2no7O3NLGBAZTL1dkVXQXIieqajl2j98ys6cYg9naM7nmacJC/RjTLwW7hHBAVw2JJqV+3w/CBnODisaBEqpvcAStA+hXEQGdeiqDOfM3HEDuXp4HL9dcYC5HacEyAAAIABJREFUL2zwaAPOyC9jTHwkIwb08vk0u3RrDp9lnuZ4SbWuPHvpd+HOt2HuX3VkU+yw1rWcqoph+yuQfjcMvdb9idNmQWWh7q7ni4IDOgLrLGoKPbc6g8iQAG7xUrYkwM/GrRPiWZtxijc2HeOa59bxysajXDOiLwsmJ7JgciL3XTGY3948plHIJMWEcbykqlnI5ZHTFQyKDsXP5jnRMT0xiu3HzugCiz2QfSdKqa23c9vF2pzX0sy0MauQgyfL+fSgtRD1DYcLuWRwNAF+Tbeta0f2I6eoigPtYMIytMZKy9EbReQQkA18ChwBPuzgdRnOkdBAf/5+1wQWfSWdEyXVzFn0Oa9vOtbsmLoGO4cLyhnaN4L0QVHszCnxqG2UVdexM1cLmW1Hz+gQ3Guf1PWbnKTNat2Zbvdb2scw4V7Pix1ytS6Z4a0fhROnY7yNAmLzkSLWZhTw4LQUegV775h7+4QEGuyKh9/ZTe/QAN791mX84fax/N+s4fzfrOE8Mms4Q+LCG49PjglrFep6tLCSpBYdAFuSPqg3pdX1rbsC9hCc2tEt4+MJ8JNm5kxoMj9Z8U/kFVeRfbqCS1Oim41fPbwvIvCt/2zj1r9u4Na/buB7b2yn3JThaBesaBBPApOAg0qpZOAqYGOHrsrQLogIs8f0Z80PpzJqQC/+sT6r2f7s0xXUNSiG9YsgPTGKqroGj09im7KLGqNFPJpFhs1u3Zlu+ys6ac1bDaDgXjpG31d5ctChtaHRuqqpRZRS/H5FBrERQdw9Ocnn8UkxYfzwmqE8MnMY73/ncsZ68Fc4cZZiP+KIZHKGuCa5KdHuyvhEHfvRU/0Q246dYWDvEOKjQkmJDW8lCJy/owP5pT7P5cx3cDqoncRGBPHtaUPoHxlMgJ8Nf5uND3ad4Mdv7TRmp3bAioCoU0oVoqOZbI5eEBM6eF2GdiQyNIDrxwzgcEEF+SVNDYGcwiCtX0TjzcrTzf/zzEKC/G2MTejt2Zk9IF03pHFqAid26bLb4xb4XmTaLDiTrU1I3sjdpLWHNtSo+vRgAZuOFPHdK4c0+gB88d2rUvnm1JRm5gxPJDsFhCNy6WRZNdV1dhI9lGhvnBcTRlRoQI/1Q2w7eoZ0x+8qrV9EMwFRU9/AnrxSAv1sHCms9JkRvSHzNNFhgaT1jWi176Hr0njtvkm8dt8kXr9/Ej+dkcaHe/J5aV2WmzMZ2oIVAVEsIuHoonn/EZE/AW0IWjd0BS51RH5sONyUeXowvwx/m5ASG86AyGD69gry+DS74fBpLk7qw+SUaPYeL3X/H9pma96ZbvsruiGSsxKsN5wlM7yZmSqLdA2meOvPJ3a74vcrM4iPCuH2izvGdRYZGtAs1NXpsE72oUGICOmDoth2zJqTtjtxvLiKEyXVpA/S2ldavwjyiqsa8z725JVS22DnulH9aLArr2a2qtoG1hw4xdShsdi8+HSc3HfFYGaP7s9vVxxgg8m0PiesCIg5QCXwA2AFcBhdk8nQjRjerxd9wgL53KVG04H8MpJjwgj0tyEijE+Mcvs0W1BWw4H8Mi4dEs34QVHU2xW7cj3E7zs70x1aBbuWwPDrIbSP7wX26q81EG8CIm+bfm1D0bllu0+w93gpP7h6KIH+lmIyzgrXUNeWIa7eSE+MIvNUOcWVnVjXqgNwPmg4NVPnk/9Bhx/CqYXOn+jege3Kh3tOUFZdz60T3NQFc4OI8LtbxpASG87C17d3aEmTzFPlHC+u8n1gF2JXbnFTaR0feP0f4yjZ/YFSyq6UqldK/Vsp9WeHycnQjbDZhMmDo9lw+HSjbTbjZClp/ZpU9vRBUeSeqeJUWfO+1E6t47KUGMY5ngg9mkWSp0BgOKx4GKqLrZmXnAy/Ho5v0z0o3PXGzt2ky2i0LJHuhuq6Bn6/8gA/eHMHaX0juGlcyzqR7YtrqGvLEFdvOP+e2y2GenYXth0tJjjAxvD+uhS783fmNGtuPXqGhD4hTEzqQ6CfzauAeGNTDknRoUwabOFBw0FYkD9/XTCeuno7s/60nlc3Hm33aDG7XbHgn1/yf++eRU2yTqKsuo5b/voFP1yyw9LxXgWEUqoBsItIG3ocGroq/7+9O4+PsroXP/75zmTfV0JISAJJIETWEEAWAWUR1ypuWL0u1Z9V3Leqrfb2emvtgmtrtdW63KtXVKQuxQoKqAjIkrALgYQlC4FsZCOBkOT8/phnJpPMkExIJhOS83698mLmeZ4kh8kzz5nnnO/5fqekRFJcdYIDZcc5frKRgop60uw7COuk6aHWF6t1ueWE+HkxMi6UyCBfhkQFnn5i1cvXEpVUXWTJCtu2Yl57Jt1pKWv6/fPw2lQ4tL71/sJNltoGvo7j0PbW55Vz0UtreGV1HpePHcT7d5zbbrhpd7APdXUlxNVqTHwYZpPYVhz3FVn5xxgdH2abw4kL8yfI14ucIzUopcjOP0ZGQjheZhPJA4IcIpys8kpr2XiwgusmJJx27crpJEcH8fm90xgVH8qTn+zk2r+tZ9PBCnKO1JBzpIbckpoudRrbCisprjpB1sGzJ1R5dU4pDY3NrM4pZbML6XVcueeuBXaIyD9E5GXrV5dbqvW4qcmWCJC1uWW2W/1hdpN+5wwKwcdscrj4r80r49yhkbYL3riEMLIPHTt9lIg1rfTYG0+bctspn0C44q9w41JoaoC35sGyh+FEtaUMamFWh+GtK3cf5adv/EBTs+Ld2ybx/LVjiQj0cb0NZygpsiXU9WBZxyGuVoG+XqQNDO5TkUwnTjXx4+EqMhJsGXoQEYbFWCKZiirrOVp90rY/rc0Etr0PNhXgZRKuGn9md4BJUYG8d/sk/nT1aHJLa7nmtfVc+OJ3XPjid8x+/juuem3dGacBsaayqTnZyL6SsyNUecWuI0QF+RAd7Msfl+d0GOnlSj7ppcaXdpZLjAwgLsyftbnltk92aQNDbPt9vcyMig9t9Wk2v7yOwmP1/L/zWsJKxyeGszS7iPyKOluIZysjLoMZj8GkO86soSmz4K71sOq3sOE1yPkSJt9tqdXdTgdxsOw4D3ywlfTYED66czIBPp1Pl36mkoyIpf1lxzlUcZzzUqM6+I4W4xPD+TirkKZm5fY7nZ6ws6iKU03KNv9gNXxgCP/eWWwbnhxvF+H0zy1FVNWdIjSgZY1KQ2MzH2cVMmvEgA5L7rZHRLgmczAXpA1g44EKrJfEkuoTvLwql0teXsOdM5K554IU/Lxdi3IDy8U2MTKAQ+V1ZOcfazVc2xudbGzim5xSLhsTS3psCE99uovv9rU/id9eRbkEAGPeweGrm9uu9QARYUpyJOv3l7O7uJoAHzPx4a3HyccnhrOtsJLPth1GKcVa6/xDSssCJesnv9POQ3j7w/m/BP9w5/td4RsEF/3eksLDN8hS4Q5g8EROnGrikpfX8NiS7bbJ3bqGRu58NwuzSXjtxvE92jlAS8TSxgMVLoW42stICLcUecrrXRE3uw5Xcd4fV3V6Etb6ydo6v2KVNjCYyrpTLN91BH9vs2140zqB3XaY6evdRyk/3mCr29FVkUG+XDQqlouNr1umDuHrh2Zw+ZhB/GV1Lte//oPToaJfLNnG3f+X3WpbbkkteaXHuXVKEhGBPmccqnyk6gR3vZvF/L+upaGxc9UMO2tdXjm1JxuZmz6Q6yYkEB/uz5+Wtx9W3t79/yfWByLirC61dhaamhJFVf0plu0oJjUm2CFs8JYpSaTHhnDf+1v42dubWLa9mAHBviRHt6wcHhYTTJCvV88MiwyeYCmeNPOXMPJqiEhmZ1EVuw5X88HmAmY//y2fbzvME0t3kHO0hpcXjGNwhGvDO93JGur6TU4J0HGIq73Z6TEMjQrkoQ+3cbTayeS8hyzbXkxBRT0r95S4/D2rc0p4fc1+rsqIJyrIt9U+63Dmil1HGTM4FC/jLtb6ybttB7F4UwGDQv2Ynuq+YmARgT48f91Ynp0/ii35lXzepkTstoJKPtxcyLLtxa1CZq0FueaeM5CMhLBOvxeam5Ulhcvz3/L17qNk51fywab8jr+xC1bsOkqgj5nJyZH4eJl4cPYwdha1v0ixvQ7C/srh+rJVrVezpiooq20gzcmio0Fh/ixdOJWnLk1nw4EKvs8tY2pKVKsJQrNJGDs4jCxjMru+oYlFy3P488p97mm0ly/MfAyu/geYWuZI/udnE4kN9efe97fw6dbDPDxnGNOHea6yYFJkIHmlroe4WgUZETfHTzay8L1st3+SdNXaPEuwoqtrCfLL63hg8VbSBobw2ytGOuy33jE0NqtW8xOxoX4E+3mRY7eiuqCijjX7Srkmc3CPDLtdlzmYtIHBvPDV3lZ1yRetyCEi0IeBIX78aUXLmP2KXUcZHR/KoDB/MhLD2V963OXQ0fqGJn76xg889clORg8O5euHZjAhKZyXV+V2W+W9RctzePeHlszKTc2Kr348ysy0AbZhtCvGxbVKGeNMex2EOs1j7Sw2IMSPYTGWk+J0Y6Zmk3DbtCGseHA6N0xK4GdThzgck5EYTs6Rar7cWczcF7/lL6tz+XM3nuDtyTp0jMTIAKYPi+afC6fw60vTuWtmMgtnprj9d7fHOjHtaoirvWExwfzx6tFkHTrGM8t+dEfzOqWq/hQ7Cisxm4T1+8s7jNKpb2ji5+9moZTitRsznK5YDw/0YUCw5a7Cfn5CRBwmqv+8ah/eZhMLnNTtcAeTSXhk7nAOltexJMtS6XhdXhlr9pWxcGYy981KZUt+JSt3l3Ck6gRbCyq58JyBQMuQ65aCju8ilFI8sXQ7Gw5U8Pv5o3j3tkkkRgby6IVplNac7JbKe8dPNvLat3k89elOWyLErQXHKKs9ydz0GNtxZpPwyNxh7f6s9jqIMSJSLSI1wGjjcbWI1IhIx8lTABGZJyI5IpIrIg4V4UQkUURWish2EflGROLttmeLyFYR2SUid7ry+zTXTDGimdI6mFSLDw/gmStHMSreMco5IyGMZgV3vpuNt8nEfbNSaWhqZvOhM6gv3QmWEMlK25vSy2ziZ9OG8Ni8NJdW2bqTdaLa1RDXti4dPYjbpw3hnfWHeGDxFp76ZCdPfbKT51fktJuy3R027C+nWcG1mfFU1p3ix+L23/JPfrKTPUeqeWnBOOeBCwbrh5JxCa3np4bFBLPHCIHNLallSVYh/3FuIrGhnetou2LWiAFkJITx0tf7jHU0OcSG+nHjuYlckxlPUmQAi1bktAwvGRdba6hy23mIjzYXsDqn9fDcO+sO8snWwzw4exgLJraE7k4cEsHM4dG8+k0e1V2sMrjxYAWNzYpQf2/uX2xZKLh811G8zcL5aQNaHTtvZGy7P+u0HYRSyqyUClFKBSulvIzH1uchp/s+K2OR3SvARUA6cL2IpLc5bBGWcqKjgaeBZ43txcBkpdRYYBLwuIi0qU6vnakrx8UxJj7U6YXfVZlJEWQmhnPvBSl8cf95/Hz6ULxM0mqltjsUHquntOakbc1Gb2JNzudqiKszj1+UxiWjYvluXxnLdhTz2bbDvLwql693uz4P0B3W5ZXj522y3ZWtbWeYKetQBR9nF7JwZrLDBaitS0bFcunoWIfQ47SBwdScaORI9Qle+Gov/t5mFs5M7vp/pBNEhEcvTONI9QnufDeLLfmV3DcrFT9vM95mEw/OGcaeIzU8t2IvQ6MCbcMz/j5m0mNDWq0f2lFYxaNLtnPrW5u4+71sSmpOsOlgBb9dtpvZIwZwz/mOd7uPzB1OVf0p3uhiDql1uWX4mE0svuNcmpoVd76bxZc7jzA5OarDTMZtuTPUYyKQq5TaDyAii7Gk7bC/f04HHjIer8aYGFdK2Q/m+eJi3QrNNWMGh/HpPdO69DOCfL1YclfrmtDjEsLcHoljnX/ISGg/w6onWO8gOsri2h4vs4lXbmhZKd7Y1My0P6zmg035zBs5sMttdNXaXEvurcERAaQMCOL73DJ+PsPxgq2U4o9f5hAV5MPdTi56bS2YmOA0Kmm4EW69NLuIZTuKue+CFCLbTHL3hMnJkZyXGsU3OaUkRQZw9fiW2iGXjR7Eq9/ksedIDQsmti46NT4xnA82WUrVeplNLFqRQ1iANzdPTuLVb/JYs68Ub7OJ+HB/nrt2rNO73ZFxoVwyKpY3vj/AjecmMiDkzEJ71+aWk5EYRtrAEF68biy3vbMZgJ/P6PxUsjsvvHGAffX3QmObvW3AfOPxlUCwiEQCiMhgEdlu/Iw/KKUOo/VqU5Kj2FFURVVd126R25N16BiBPmanWT09LTk6kGA/L4fhk67wMgoYfbu3tMdy/pRUn2BfSS3TjNTa01Ki2HSwwmnJ0O9zy9hwoIJ7zk/pUmix9e/5wld7CQvw5vbpnouLefTC4fiYTfxiXlqrbL4mk/D4RWl4m4XLRrce0BiXEGZLl79hfznf7i1l4cxkHpwzjH8/cB4jYkM42djM3/4jk1D/03+Kf3DOMBqbFBe/vMYWat4ZFda63cYw8qwRMTw0Zxj+3mbm2M0/uMrTn8wfAWaIyBZgBlAENAEopQqMoacU4GYRcfjficgdIrJZRDaXlrpWlUpzn6kpUShlqR18Okoph1xPnZGdf4wxg8NsIZK9SbCfN1lPzuHiUd37Sf/azMEoLLWwe8I6I3rJWnthSnIkJ041s6VN1lmlLJly48L8uX5S19YqhAZ4MzDEUpr0rhkdF3Vyp9HxYWz9zzlcPMpxfH7m8AHs+M2FjIxrPTxrX9tj0YocYkJ8ucmoPZIcHcTiO85l85OzO1xMlzIgiE/unkpcmD/3vb+F297ZTFEnPhisN/52U+zqZtw3K5Xsp+ac0WJDd77LigD7EIR4Y5uNUuqwUmq+Umoc8CtjW2XbY4CdwHltf4FS6u9KqUylVGZ0tOfCGzWLsYPD8Pc2n3aYKb+8jpve3Mik3608ozTMdQ2N7C6ucVih25tYM+N2p8ERAUxLieKjzYW2ok3utDa3jLAAb9KNRHuThkZiEsdw1+W7jrK9sIoHZqfi6+X6CuTTGRkX2urC6knt3Q05W20dF+bPgGBfXl+zn00Hj3HvBamtjhMRl1dppw8KsYWa/7C/nDnPf8ub3x9w6W+/Nq+MIF8vxrSZX3S1Dkpb7uwgNgGpIjJERHyABVhqWtuISJSIWNvwBPCmsT1eRPyNx+HANCDHjW3VuoGPl4mJQyIcJjQbm5r527d5zH3xW7bkVxIR4MMfXMgD09a2giqa2sTQ9xcLJiRQVFnPmn3uvVNWSrE2t4zJQyNt4+Sh/t6Mig+zrYsAS1z9cytySI4O5MpuypT77PxRLF049YwvZp5kTZdfUFFPQkQA17qYmvx07EPNJw6J4Ol//cj8V9exu4NosnW5ZUwaEtFtd9hu6yCUUo3APcByYDfwoVJql4g8LSKXG4fNBHJEZC8QAzxjbB8BbBCRbVjqYC9SSp09OXX7sakpka0q1zU2NXPr25t49t97OC81mq8ems5j89LYVlDJV0ZKBldZJ6jbpnDoD+akxxAR6MPijZZhpvLakzz0wVbmvfid0+JN76w7yLV/W++wvSMHy+s4XHWi1RAFwNTkSLYVVFJ7spGq+lM8/vF29pXU8vDc4d12MYoO9iWuk+tHehPrne2Dc1K7rfZIfHgAb90ygZcWjKWwoo7L/vw9f/xyj9O/eVFlPQfL6xz+dl3h1oQ1SqkvgC/abPu13eMlwBIn3/cVMNqdbdPcwzpuvS6vjPkZ8fzhyz2s2VfGb68YyQ2TLHHf8zPieO3bPJ5bsZdZI2JcXjOQfegYydGBhAW4Pztrb+PjZeKqjDjeWnuQt9ce4KWV+zhmBANkHzrmcFFYml3ItsIqak6cIrgT4/m22s/Jka22T02J4q/f5PHcihyWbS+mrPYkP58xlHnn9FxkVW93TeZggv28uHxM99YeERF+MjaO6anRPPPFbv76TR5f7Cjmd/NH2dY0gX3d7sjT/ahO630zfdpZzb5y3b+2H+b1NQe4aXIiN56baBub9zKbeGjuMHKO1vD5NteC06w1BHrz/IO7XTchgcZmxW8+/5EhUYF8fNcUzCaxJVS0qqo/xY4iS8W/Q+WO1dQ+2VLEY0u2O+R9Ona8gaXZhcSG+jGkTbLB8Ynh+HiZeGvtQaKDffn07mk8cdEIjy9O7E1C/b25bkKC21KDhAf6sOiaMbx3+yQU8NPXN/DoR9tsCSvX5ZYRFeS8bveZ6tmUl1qfZ61ct2rPUf69s5jxieE8eUnb9ZFw8chY0mPzeP6rvVw8KrbDW/IDZcc5VneqX84/WKUMCOLRC4cTFuDN9RMSMJmEMfGhrM0t59ELW46zroIGS3W7thE3/7cxn40HKvhiRzGPXZTGTycm8Pn2w/zX5z9SXX+K31x+jsNEu5+3mV9dPIKmZsVNkxN7ZRRZfzE1JYrlD0znpZX7+Pt3+1mdU8KvLzuHdXnlTEmO6tYgCd1BaN1uSkoky3YUExXky19vyHB68TeZhEcvHM6tb29i8ab8DiNXso0Qy/58BwE4LEabmhLFK6tzqT5xyhYaui6vHF8vEycbmzlYdtzhZ+SV1DJzeDQNjc08+clO/rxqH0erTzJmcBi/nz/KVia0rZunJHX7/0c7M37eZh6bl8ZlowfxxNLt3Pf+FqB7h5dADzFpbjB7RAwj40J49cYMYtpZDTpzeDSTh0by7Bd72HOk/eiML3ceIbpN2nHNsjixWcGG/S05sNbmljFxSAQxIb4cbDPEdOx4A+XHG5iaHGWrtBbm78N/XpbO0rumnLZz0Hon+5DYcQlhXJDW+cVw7dEdhNbtYkL8+Ne95zEhqf0i8yLCSwvGEuTnxZ3/m0VVvfMV2EerT7A6p4Srx8frMe82MhLD8PM22SYoraugp6ZEkRQZ6HAHsb/MUhozeUCgrdLa8genc+vUIX2iml1/ZA2J/efCqUQHd296Et1BaB41IMSPV2/IoPBYPQ9/uNVpWumPNhfQ1Ky4roux5X2Rr5eZCUkRtsWJtlXQyUYHUd66g8grsTzXd2KaK3QHoXlcZlIET14ygq93l/CX1bmt9jU3Kz7YXMDkoZG2ZHhaa1OSo9h7tJaSmhMtq6AHhZAUFUhZbQM1dumjc0tr8fEyER/e81X3tLOP7iC0XuHmKUlcOS6OF77e2yqH/rq8cgoq6nuscMzZyDoxuT6vnHV55UweGonZJAyJsnQC9qGueSW1DI0K1MNJmkt0B6H1CiLC764cRdrAEB5YvJV846K2eFM+YQHetupdmqNzBoUS4ufFexvyKaqsty2asxbusR9myiut1cNLmst0B6H1Gv4+Zl67MQOlLEVODlfWs2LXUa4cF+dyorP+yGwSJidHsvGAJZLJugraWhfbOlF9srGJ/Io6kqP1UJ3mGt1BaL1KYmQgLy0Yx+4j1Vz16joamppZMKFrqaT7A2uKE/tV0AE+Xq1CXQ+V19GsILmDQvWaZqU7CK3XOT9tAPfPSqW46gTjEsI6zKGvtdQZb7uS1j7UNbfECHHVQ0yai/RKaq1Xuu+CVARhxnBd58MVydGB3HtBikNZ0qTIQFbusWTNzTM6iKF6iElzke4gtF7JZBLun53q6WacNUSEh+cOd9huH+qaV1pLXJh/l0qDav2LHmLStD7MPtQ1r/S4nn/QOkV3EJrWh1lDXfeXHTdCXPXwkuY6t3YQIjJPRHJEJFdEHneyP1FEVorIdhH5RkTije1jRWS9iOwy9l3nznZqWl9lDXX9YX85dQ1NeoJa6xS3dRAiYgZeAS4C0oHrRaRtYYBFwP8opUYDTwPPGtvrgJuUUucA84AXRaT/1ZnUtC6yhrqu2m1ZnZ6ih5i0TnDnHcREIFcptV8p1QAsBn7S5ph0YJXxeLV1v1Jqr1Jqn/H4MFAC6HAWTTsDSZGBHDGqx+k7CK0z3NlBxAEFds8LjW32tgHzjcdXAsEi0qrihYhMBHyAvLa/QETuEJHNIrK5tLS02xquaX1JkjEPEeLnRVRQ/6vnrZ05T09SPwLMEJEtwAygCGiy7hSRWOB/gVuVUs1tv1kp9XelVKZSKjM6Wt9gaJoz1iy4yQOCurUcpdb3uTMgugiwT8EZb2yzMYaP5gOISBBwlVKq0ngeAiwDfqWU+sGN7dS0Ps0a6pqih5e0TnLnHcQmIFVEhoiID7AA+Mz+ABGJEhFrG54A3jS2+wD/xDKBvcSNbdS0Ps8a6qrXQGid5bYOQinVCNwDLAd2Ax8qpXaJyNMicrlx2EwgR0T2AjHAM8b2a4HpwC0istX4GuuutmpaXzYsJpiFM5O5fMwgTzdFO8uIUo4lHs9GmZmZavPmzZ5uhqZp2llFRLKUUpnO9nl6klrTNE3rpXQHoWmapjmlOwhN0zTNKd1BaJqmaU7pDkLTNE1zSncQmqZpmlO6g9A0TdOc0h2Epmma5lSfWSgnIjVAjqfb0ctEAWWebkQvpF8XR/o1cdRfXpNEpZTTbKd9qXp5zulWA/ZXIrJZvyaO9OviSL8mjvRrooeYNE3TtNPQHYSmaZrmVF/qIP7u6Qb0Qvo1cU6/Lo70a+Ko378mfWaSWtM0TetefekOQtM0TetGuoPQNE3TnOoTHYSIzBORHBHJFZHHPd0eTxCRwSKyWkR+FJFdInK/sT1CRL4SkX3Gv+GebmtPExGziGwRkX8Zz4eIyAbjfPnAKHHbb4hImIgsEZE9IrJbRCb39/NERB403jc7ReR9EfHr7+cJ9IEOQkTMwCvARUA6cL2IpHu2VR7RCDyslEoHzgXuNl6Hx4GVSqlUYKXxvL+5H0vZW6s/AC8opVKAY8BtHmmV57wEfKmUSgPGYHlt+u15IiJxwH1AplJqJGAGFqDPk7O/gwAmArlKqf1KqQZgMfATD7epxymlipVS2cbjGixv+jgsr8U7xmHvAFd4poWeISLxwCXAG8ZzAS4AlhiH9KvXRERCsdTnonAkAAADu0lEQVR7/weAUqpBKVVJPz9PsCwa9hcRLyAAKKYfnydWfaGDiAMK7J4XGtv6LRFJAsYBG4AYpVSxsesIEOOhZnnKi8AvgGbjeSRQqZRqNJ73t/NlCFAKvGUMu70hIoH04/NEKVUELALysXQMVUAW/fs8AfpGB6HZEZEg4GPgAaVUtf0+ZYlp7jdxzSJyKVCilMrydFt6ES8gA3hVKTUOOE6b4aR+eJ6EY7mDGgIMAgKBeR5tVC/RFzqIImCw3fN4Y1u/IyLeWDqH95RSS43NR0Uk1tgfC5R4qn0eMBW4XEQOYhl6vADL+HuYMZQA/e98KQQKlVIbjOdLsHQY/fk8mQ0cUEqVKqVOAUuxnDv9+TwB+kYHsQlINSIOfLBMLn3m4Tb1OGNs/R/AbqXU83a7PgNuNh7fDHza023zFKXUE0qpeKVUEpbzYpVS6gZgNXC1cVh/e02OAAUiMtzYNAv4kX58nmAZWjpXRAKM95H1Nem354lVn1hJLSIXYxlrNgNvKqWe8XCTepyITAPWADtoGW//JZZ5iA+BBOAQcK1SqsIjjfQgEZkJPKKUulREhmK5o4gAtgA3KqVOerJ9PUlExmKZtPcB9gO3Yvmw2G/PExH5L+A6LNGAW4Dbscw59NvzBPpIB6FpmqZ1v74wxKRpmqa5ge4gNE3TNKd0B6FpmqY5pTsITdM0zSndQWiapmlOeXV8iKZpACLShCWM2OoKpdRBDzVH09xOh7lqmotEpFYpFdTOfi+73D2adtbTQ0ya1gUicouIfCYiq4CVIhIkIitFJFtEdojIT4zjkoz6C2+LyF4ReU9EZovIWqMGw0TjuEAReVNENhrJ9PpdZmKt99B3EJrmojZDTAeUUleKyC3Ab4HRSqkKa7popVS1iEQBPwCpQCKQiyXL7i4sKWK2YakxcDlwq1LqChH5HfCjUupdEQkDNgLjlFLHe+5/qmkWeg5C01xXr5Qa62T7V3ZpKQT4nYhMx5LyJI6W1NkHlFI7AERkF5YCPUpEdgBJxjFzsSQYfMR47ocl/YV9wSNN6xG6g9C0rrP/dH8DEA2MV0qdMjLJ+hn77PP4NNs9b6blvSjAVUqpHPc1V9Nco+cgNK17hWKpQXFKRM7HMrTUGcuBe42soojIuO5uoKa5SncQmta93gMyjWGjm4A9nfz+/wa8ge3GMNR/d3P7NM1lepJa0zRNc0rfQWiapmlO6Q5C0zRNc0p3EJqmaZpTuoPQNE3TnNIdhKZpmuaU7iA0TdM0p3QHoWmapjn1/wHA1ERKdTzv/AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "q1q2_df.plot(x='Frame')\n", "plt.ylabel('Fraction of native contacts')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, we can create a Q1 vs Q2 plot to characterise the transition of AdK from its opened to closed position. " ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Q2')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEGCAYAAACgt3iRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3zU9f3A8dcnew+yd0IIkLBCCEMQZDgQRQRRwT1aWqu1Vv05WqvWarUWF2odtQ60SJWh4GDJlB3CTIBMMi57z0tufH5/3BETCBAgl7uQz/PxyIPkO+7eR+De9/m8P0NIKVEURVGUi2Fn7QAURVGU3k8lE0VRFOWiqWSiKIqiXDSVTBRFUZSLppKJoiiKctEcrB1Ad/H395fR0dHWDkNRFKVX2bdvX4WUMuBiH+eSSSbR0dGkpKRYOwxFUZReRQiR1x2Po7q5FEVRlIumkomiKIpy0VQyURRFUS7aJVMzURRFsRSdTkdhYSFardbaoVwwFxcXwsPDcXR0tMjjq2SiKIpyDoWFhXh6ehIdHY0QwtrhnDcpJZWVlRQWFhITE2OR57BYN5cQ4mMhRJkQ4sgZzgshxCIhRJYQ4pAQIqndubuFEJnmr7stFaOiKEpXaLVa/Pz8emUiARBC4OfnZ9GWlSVrJp8C089y/logzvy1AHgPQAjRD3gOGAuMAZ4TQvhaME5FUZRz6q2J5CRLx2+xZCKl3ApUneWSWcBiabIL8BFChADXAOullFVSympgPWdPSpe8lBNVvL7uuLXDUBRFOSNrjuYKAwra/VxoPnam46cRQiwQQqQIIVLKy8stFqi17curZtHGLHLKG6wdiqIoVlRYWMisWbOIi4ujf//+PPTQQ7S0tLB+/XpGjRrFsGHDGDVqFBs3buzx2Hr10GAp5YdSymQpZXJAwEWvBmCzZo8Mw07A8tRCa4eiKIqVSCmZM2cON954I5mZmWRmZtLc3MwTTzyBv78/q1ev5vDhw3z22WfceeedPR6fNZOJBoho93O4+diZjvdZgV4uXDEwgBWpGgxGtTOmovRFGzduxMXFhXvvvRcAe3t73njjDRYvXkxcXByhoaEADBkyhObmZlpaWno0PmsODV4FPCSEWIqp2F4rpSwWQqwF/t6u6H418LS1grQVc0dF8OCSVHZkVzAx7tJthSmKrfvr6jTSi+q69TETQr14buaQs16TlpbGqFGjOhzz8vIiOjqarKwsEhMTAVi+fDlJSUk4Ozt3a4znYrFkIoT4EpgM+AshCjGN0HIEkFK+D/wAzACygCbgXvO5KiHE34C95od6QUp5tkJ+nzAtPhBvV0eW7StUyURRlE6lpaXx5JNPsm7duh5/boslEynl/HOcl8CDZzj3MfCxJeLqrVwc7blhRChfpRRQp9Xh5WKZWayKopzduVoQlpKQkMCyZcs6HKurq6OkpIRBgwZRWFjI7NmzWbx4MbGxsT0eX68uwPc1c0eF06I38sOhYmuHoihKD5s2bRpNTU0sXrwYAIPBwGOPPdY2ouu6667jlVdeYcKECVaJTyWTXmR4uDcDAj1Ytk+N6lKUvkYIwcqVK1m2bBlxcXH4+flhZ2fHn//8Z9555x2ysrJ44YUXSExMJDExkbKysh6NT63N1YsIIZg7KpxXfjxGbkUjMf7u1g5JUZQeFBERwapVqwDYsWMH8+fPJzU1lWeeeYZnnnnGqrGplkkv0zbnRLVOFKVPGz9+PHl5eSQlJZ374h6gkkkvE+TlwqSBASxPLVRzThRFsRkqmfRCc0eFU1yrZWd2pbVDUZQ+wzQAtfeydPwqmfRCV8YH4eXiwLJ9Bee+WFGUi+bi4kJlZWWvTSgn9zNxcXGx2HOoAnwv5OJozw2JoSzbV6jmnChKDwgPD6ewsJDevKDsyZ0WLUUlk15q7qgIvtiVzw+Hipk3JtLa4SjKJc3R0dFiOxReKlQ3Vy81Qs05URTFhqhk0kudnHOSkldNbkWjtcNRFKWPU8mkF1NzThRFsRUqmfRiQV4uTIxTc04URbE+lUx6OTXnRFEUW6CSSS93VUIQnmrOiaIoVqaSSS93cp+TNWkl1Gt11g5HUZQ+SiWTS8DcUeFodUZ+OKz2OVEUxTpUMrkEJEb4EBvgruacKIpiNSqZXAJMc04i2HuimhNqzomiKFagksklom3OSapqnSiK0vMsmkyEENOFEMeFEFlCiKc6OR8lhPhJCHFICLFZCBHe7tw/hBBHzF+3WjLOS0Gwt3nOyb5CjGrOiaIoPcxiyUQIYQ+8C1wLJADzhRAJp1y2EFgspRwOvAC8bL73OiAJSATGAo8LIbwsFeulYu6ocIpqtezMUXNOFEXpWZZsmYwBsqSUOVLKVmApMOuUaxKAjebvN7U7nwBslVLqpZSNwCFgugVjvST8MudEdXUpitKzLJlMwoD2M+kKzcfaOwjMMX8/G/AUQviZj08XQrgJIfyBKUCEBWO9JDg72DGuvx8/HimmsUVv7XAURelDrF2Afxy4QgixH7gC0AAGKeU64AdgB/AlsBMwnHqzEGKBECJFCJHSmzet6Q6ZpfXc+Z89rE8vJdrPHSGsHZGiKH2JJTfH0tCxNRFuPtZGSlmEuWUihPAAbpJS1pjPvQS8ZD63BMg49QmklB8CHwIkJyf3yapzbZOONzZk8PmuPNyd7HluZgJ3jIvC0d7anxMURelLLJlM9gJxQogYTElkHnBb+wvMXVhVUkoj8DTwsfm4PeAjpawUQgwHhgPrLBhrr2MwSpbuzWfh2uPUNOuYPyaSx64aiJ+Hs7VDUxSlD7JYMpFS6oUQDwFrAXvgYyllmhDiBSBFSrkKmAy8LISQwFbgQfPtjsA2YeqrqQPukFKqIoDZntwqnl+VRnpxHWOi+/HcDQkMCfW2dliKovRhQspLo3coOTlZpqSkWDsMiyqqaeblH4+x+mARod4uPD0jnuuHhyBUgURRlAskhNgnpUy+2MexZDeX0k20OgMfbs3hX5uzkBIenhbHA1fE4upkb+3QFEVRAJVMbJqUkjVHSnjph6MUVjczY1gwT18bT0Q/N2uHpiiK0oFKJjbqWEkdf12Vzs6cSgYHe7Lk12MZH+tv7bAURVE6pZKJjalpauX19Rl8sSsPL1dH/jZrCPPHROKghvoqimLDVDKxEQajZMmefF5bd5y6Zh13jIvi0asG4uPmZO3QFEVRzkklExuxYHEKPx0rY2xMP/46awiDg9W6loqi9B6q78RGDAs3zROp1+pxsFO/FkVRehf1rmUjHrlyIB/dlUxxbTMz3/6Zr/YWcKnMAVIU5dKnkokNuTIhiB//MInECB+eWH6Ih5ceoE6rs3ZYiqIo56SSiY0J9nbhi1+N5f+uGcQPh4u5btE29udXWzssRVGUs1LJxEbsyK4gq6weAHs7wYNTBvDVby7DaISb39/Jp9tzrRyhoijKmalkYgMMRsm9n+xlxqKfWbzzRLtaiSSinyt6o2Tp3oKzPYSiKIpVqaHBNqCwuokWvZEQbxee/TaNV9ccJ9zXlWMl9fRzd+L/rhnEnZdFWTtMRVGUM1LJxAbklDcCMCcpjHc3ZdPQoudYST1XxgexaH4ibk7q16Qoim1T3VxWZjRK3tuSDcC7m7IJ9Xbh1uQIwnxc+elYKe9vzkZvMFo5SkVRlLNTH3mtxGCUfH+4mHc3ZnG81FR4f/Wm4dw4MgwnBzsaW/Q8vyqNRRuz2J5dyVvzEgn3VasFK4pim1TLpIfpDEa+Tingqte38PCX+zGYi+0jIny4ZXQETg6mX4m7swP/vHkEb81L5HhJPTPe2sYPh4utGbqiKMoZqWTSQ1r0Bv67O48pCzfzf8sO4eJoz3u3J7HukUkEeDoTF+jR6X2zEsP44eGJxAR48Lv/pvL0isM0txp6OHpFUZSzU91cFiKlZHduFYODPVmRquHDrTmU1GlJjPDhrzcMYergQIQQ1Gt1lNe30D/A/YyPFennxrLfXsZr6zJ4f0s2e09U8d7tScQFefbgK1IURTkzlUwsoLqxld/9N5WdOZVtx8bG9GPhzSOYMMCvw57tJ0dy9ffvvGVykqO9HU9dO5gJA/y479O9vPVTJu/clmSZF6AoinKeVDLpZmvTSvjN5/s6HPvV5TE8c31Cp9fnVDQAEHuWlkl7sQEe6AySxAifiwtUURSlG1m0ZiKEmC6EOC6EyBJCPNXJ+SghxE9CiENCiM1CiPB2514VQqQJIY4KIRaJ9h/nbdRfvjnSlkhCvF1Y9tvLGBrmxcr9GioaWjq9J6e8EXs7QaRf10ZqbcssB2DSwIDuCVpRFKUbWCyZCCHsgXeBa4EEYL4Q4tSP5wuBxVLK4cALwMvme8cDE4DhwFBgNHCFpWLtLj5ujgAsf+Aydj49jeTofrx+SyL1LXqeXnG40yXlc8obifB1xdnBvkvPsTWzgmAvlzMW7BVFUazBki2TMUCWlDJHStkKLAVmnXJNArDR/P2mducl4AI4Ac6AI1BqwVi7xWNXD+LEK9cxKqpf27GBQZ48cc0g1qeXsmxf4Wn3ZJc30D+ga4nBYJT8nFnBxDh/ekFDTVGUPsSSySQMaL86YaH5WHsHgTnm72cDnkIIPynlTkzJpdj8tVZKefTUJxBCLBBCpAghUsrLy7v9BXSX+ybEmLbjXZ1OTVNr23GjUXKispH+/l2rlxzW1FLbrGOi6uJSFMXGWHueyePAFUKI/Zi6sTSAQQgxAIgHwjEloKlCiImn3iyl/FBKmSylTA4IsN03WDs7wV+uT6ChRc+PR0rajhfVNqPVGbvcMtmaUY4QcPkAf0uFqiiKckEsmUw0QES7n8PNx9pIKYuklHOklCOBP5uP1WBqpeySUjZIKRuAH4HLLBirxQ0J9aK/vzurDxa1Hcs+OSy4iyO5tmWWMyzMm37uThaJUVEU5UJZMpnsBeKEEDFCCCdgHrCq/QVCCH8hxMkYngY+Nn+fj6nF4iCEcMTUajmtm6s3EUJw/YhQduZUUlanBSCn3DQsuCvJpF6rIzW/holxqlWiKIrtsVgykVLqgYeAtZgSwVdSyjQhxAtCiBvMl00GjgshMoAg4CXz8WVANnAYU13loJRytaVi7Skzh4cgJW1rbOWUN+Lp7ECAh/M5792RXYnBKJkYZ7vdeYqi9F0WnbQopfwB+OGUY8+2+34ZpsRx6n0G4DeWjM0a4oI8GRzsyepDxdwzIYacigb6B7h3aWTWtsxy3J3sSYr07YFIFUVRzo+1C/B9zswRoezLq6awuomc8sYuF9+3ZVZwWaxf26rCiqIotkS9M/WwmcNDAfg6pZDiWm2XhgXnVTaSV9mkurgURbFZam2uHhbp58aICB8+2Z4LQOxZZrJLKUnNr+aDLTmAWkJFURTbpZKJFcwcHsKL35sGp3U2kiuvspEVqRq+OaAhr7IJF0c77psQQ3QX1+9SFEXpaSqZWMF15mQiBET7mZJJbZOO1YeKWLlfw768aoSAy/r78dCUAUwfGoyni6OVo1YURTkzlUysIMTblTEx/Sit0+LiaFrg8eYPdpBR2kBcoAdPTh/MrMRQQn1crRypoihK16hkYiWv3zKCeq2+7Wd/D2f0Rsm6P05SizgqitLrqNFcVhLu60Z8iFfbz5MGBpBT3kh5fef7niiKotgylUxsxMllUrZlVlg5EkVRlPOnkomNiA/2wt/Dia2ZtruUvqIoypmoZGIj7OwEE+MC+DmzAqPx9B0Zz8RolPxz7bG2RSMVRVGsQSUTGzIxzp/KxlbSi+u6fM8hTS3vbspm47EyC0amKIpydiqZ2JDLzXWT8+nq2pahusUURbE+lUxsSKCnC/EhXmzL6HoRXhXsFUWxBSqZ2JhJcf6k5FXR2KI/57WmDbOqeyAqRVGUs1PJxMZMjAtAZ5Dszq0857U7syvRn0exXkpJVlkDDV1IVIqiKOdDJRMbkxzti4ujHVu70NW1NbMcB7tzz5YvqmnmX5uzuPqNrVz5+hYW/ZTZHaEqiqK0Ucup2BgXR3vGxvh1qQi/LbOCpEhf9pyoOu1cvVbHj0dKWJmqYVduJVLCqChfXB3tqdfqLBG6oih9mEomNmhkpA9bMsoprG4i3LfzZedPbpg1Nym8QzKpaGjhb9+lszatBK3OSJSfG3+YFsfskWFE+bkz5qUNPfUyFEXpQ1QysSHpRXW8uymLH44U4+Zkj95w5nrIVvMorkkDA3htfQZgqok89tVBduVUcnNyOLNHhpMU6aMWjlQUxeIsmkyEENOBtwB74CMp5SunnI8CPgYCgCrgDilloRBiCvBGu0sHA/OklN9YMl5rOVBQwzsbM9lwtAwPZwd+NzmW+ybE4OfhfMZ7tmWUE+bjSnS7bX+X7MlnS0Y5f5s1hDsvi+6ByBVFUUwslkyEEPbAu8BVQCGwVwixSkqZ3u6yhcBiKeVnQoipwMvAnVLKTUCi+XH6AVnAOkvFai17cqt4e2Nmh7kiH941ivGx/me9T2cwsiO7kpkjQjnZ6MirbGLZvkImxvlzx7goS4atKIpyGku2TMYAWVLKHAAhxFJgFtA+mSQAj5q/3wR01vKYC/wopWyyYKw9RkrJ9qxKFm3MZE9uFf4eTjx97WCKapr5bGcet/17N2Nj+nHf5TFcGR+EfSejtQ4U1NDQomdS3C9J57+783B3duDVucNVt5aiKD3OkskkDCho93MhMPaUaw4CczB1hc0GPIUQflLK9pMs5gGvd/YEQogFwAKAyMjIbgrbMqSUbDpexqKfsjhQUEOwlwvPzUxg3uhIXJ1Muy2O7e/HU8sPsTu3it25VUT0c+We8THckhzeYdvebRnl2AkYP+CXZGKU8LdZQwnxVrszKorS86w9z+Rx4AohxH7gCkADGE6eFEKEAMOAtZ3dLKX8UEqZLKVMDggI6Il4L9iO7Eru+zSFAwU1ODvYcUNiKGE+rh2G6c4YFsIPf5hIcpQvAAVVzfztu3TG/f0nnl+VxomKRoxGyY9HSkiM8MHb1REHO4GdgBnDgpmVGGqtl6coSh9nyZaJBoho93O4+VgbKWURppYJQggP4CYpZU27S24BVkope/3EiNHR/Xj1puHsL6jmYEEt//k5lw+35gAQ7OXC8HBvRkT4MCLchw/vSubT7bm8vSkLKSHM15X/7s7js50nkOYBXm/NSwTAzcmBZQ+MJz7YS3VvKYpiNZZMJnuBOCFEDKYkMg+4rf0FQgh/oEpKaQSexjSyq7355uO9npODHbeMjuCW0ab82txqIL24loMFtRwqrOFgYS3r0kvbro/xdyfQ05nSuhYyShu4Z3w0rk72vLc5G4B3N2XRojNyQ2IoSZG+Pf56cisaSS+qY8awYJXEFEWxXDKRUuqFEA9h6qKyBz6WUqYJIV4AUqSUq4DJwMtCCAlsBR48eb8QIhpTy2aLpWK0Jlcne0ZF9WNUVL+2Y7VNOg5rajlYWMPBghoOFda2nft0x4kO95+oaOKJ5Yf4x5pj3D42kjvGRRHo5WLRmKsaW/nuUBErUjUcKKjBycGO8bHT8HV3sujzKopi+4SUXV8o0JYlJyfLlJQUa4fR7UrrtBwoqOGP/ztAU6vhjNc52guuHx7KfRNiGBbufcbrxry0gWnxgbw8Z3iXnr9Fb2Dj0TJW7New6VhZ28KSlw/w59mZCQwM8jy/F6Qoik0RQuyTUiZf7OOoGfA2LsjLhWuGBPPMdQn8aeXhtuNjYvoxLMybgwU1HCmqRaszsnK/hpX7NUyM8+ffdyXj4mh/Qc8ppSQlr5oVqRq+P1REnfaXVYYj+rnyzHUJXJ0QpLq3FEVp06VkIoRwPLUILoTwl1KqnZl6gFZn4O2NmSRF+rBgUiy//WIf4/r78ehVAwHQG4xkljVwqLCGAwW11Gl1nc5POZfcikZzQiqkoKq5wzlXR3semjqA+y+PueAkpSjKpeusycS8rMnngIsQIhVYIKU8YT69DkiybHgKwJd78imu1fKX6xN4btUR4gI9+N3k2LbzDvZ2xId4ER/ixa2jz++xq0/WQfZr2J9fgxAwPtaP2AAPUvOqqdPquTExlKeujSfY27I1GUVReq9ztUxeBa4xF87nAuuFEHdKKXcBqo+jBzS16nl3UzaX9fdjzZESKhta+c/doy+qddCiN7LmSDErUjVsOl6GziAZHOzJ09cOJjbAg/e3ZLM9q5JhYd48f0NCh0ECiqIonTlXMnGSUqYBSCmXCSGOAiuEEE8Cl0bl3sYt3plHRUMLM0eE8Mn2Ezx61UCGhp25wH4mUkr25VVTVt/CilQNK1I1BHg6c8/4aGaPDCfIy5mF647zyppj9HNz4h83DePmURHYXUB3maIofc+5kolOCBEspSwBMLdQpgHfAbFnv1W5WPVaHe9vySY+xIuV+zWMiPDp0L3VFaV1WpbszuebAxryKn9Z3uyz+8YwIdYPiSlhvbkhg+ZWA/dPiOHhK+Pward8i6IoyrmcK5k8BQQJIeqAAeZjxzEtffKQJQNT4OOfT1DTpMPBTtDcauD1W0bgYN/1FXDWHCnhyeWHqNPqmBDrz8NT43j22yPckBjKFQMD2JpRzgvfpZNV1sCkgQE8e30CAwI9LPiKFEW5VJ0rmWzFVDe5C8jFVCcJAt6WUr4khEiUUh6wcIx9Uk1TKx9tMy23UtHQyvMzE4gN6NobvVZn4MXv0/liVz7Dwrx5a14i/c33/mPNMYpqtPzqsxQ2HC0lys+Nj+5KZlp8oBrqqyjKBTtXMlkIuAJRUsp6ACGEF7BQCPEeMB2IsWyIfdO/t+VQ32Ka3zE+1o+7urjZVUZpPb9fsp/jpfUsmNSfx68ehJNDx9bMloxy3J3seXL6YO67PBpnBzXUV1GUi3OuZDIDiJPtpslLKeuEEA8AFcC1lgyur6poaOGT7ScA8HR24J83jzhnIVxKyZd7CnjhuzQ8nB349N7RTB4UeNp1VyUE0ao38n/XDLL48iuKovQd50omRtnJeitSSoMQotw8RFjpZu9vzm5bOuX5G4YQ5nP2PUpqm3Q8teIQPx4pYWKcP6/dMoJAz84TxUuzh3V7vIqiKOdKJulCiLuklIvbHxRC3AEctVxYfVdJrZbPd+UBcM2QIOYkhZ31+pQTVfxh6QFK67Q8fe1gfj2xvxrOqyhKjztXMnkQ07yS+4B95mPJmOoosy0ZWF/17qYsWvRG/D2c+PvsYWcsihuMknc3ZfHmhgzCfd1Y/sB4RkT49HC0iqIoJmdNJlJKDTBWCDEVGGI+/IOU8ieLR9YHFVY3sXRvPgB/nz0MPw/nTq8rrm3mkaUH2J1bxazEUF68cWiHbX27S0OLHnsh2rYVVhRFOZMuLfQopdwIbLRwLH3e2z9loTNI5o4K5+ohwZ1esy6thCeWH6JVb2ThzSO4KSms24f0ltVp+fe2HL7Ylc91w0NYePOIbn18RVEuPWoJehuRW9HIstRCwnxceXZmwmnntToDL/9wlM925jEk1Iu3549smzvSXYpqmvlgSzZf7i2gVW9kTHQ/Hpoy4Nw3KorS56lkYiPe2pCBwSj5583DT1vKJK+ykd98vo9jJfXcf3kMT0wf1C1zQwqqmtiZXUlipA8f/5zL8tRCdAaJj5sjL84aytxR4aqYryhKl6hkYiMqGlp5eOoAxsf6dzhe26Tjnk/2Ut3Uyif3jGbK4NPnjpwPvcHIpuPl/Hd3HpuPl592fk5SGH+eEX/Geo2iKEpnVDKxEZ/fP+a02ofeYOTBJakUVjfx5a/HkRx94UvBF9c287+9BSzdU0BJnfa08zH+7rx041DGD/Dv5G5FUZSzU8nERnRWRH/x+6P8nFXBqzcNv6BEYjBKtmaWs2R3Pj8dLcXYyaYBTvZ2PDA5lgcmx6odFBVFuWAWTSZCiOnAW4A98JGU8pVTzkcBHwMBQBVwh5Sy0HwuEvgIiMC0d8qMdrs8XvKW7snn0x0nuP/yGG4ZHXFe95bVa/k6pZAv9+RTWN18xuvGxvTjpdnD1ErBiqJcNIslEyGEPfAucBVQCOwVQqySUqa3u2whsFhK+Zl5LsvLwJ3mc4uBl6SU64UQHoDRUrHamj25Vfzl2yNMGhjA09cO7tI9RqNkZ04l/92dx7q0UvSdNEMc7QU6g8TXzZE/zYhn7qhwtVKwoijdwpItkzFAlpQyB0AIsRSYBbRPJgnAo+bvNwHfmK9NAByklOsBpJQNFozTphRUNfHbL/YR0c+Nt+ePPOf+JVWNrSzbV8CS3fmcaLf5VXuTBwWQVlRHeX0Lc0eF86cZ8fRzd7JE+Iqi9FGWTCZhQEG7nwuBsadccxCYg6krbDbgKYTwAwYCNUKIFZiWuN8APCWlNFgwXqtrbNHz68Up6A1GProrGW/Xzme1SynZk1vFkj35/Hi4hFZD5422P0yL41hJHWvTSukf4M6ieSO5LNbPki9BUZQ+ytoF+MeBd4QQ92DaiEsDGDDFNREYCeQD/wPuAf7T/mYhxAJgAUBkZGRPxWwRRqPkj/87QEZpPZ/eO6bTCYm1TTqWpxayZE8+WWUNeLo4EBvowdHiug7XvTxnGC06A6+tz6BFZ+SPVw7kt5P7q31LFEWxGEsmEw2m4vlJ4eZjbaSURZhaJpjrIjdJKWuEEIXAgXZdZN8A4zglmUgpPwQ+BEhOTu5krFLv8caGDNall/Ls9QlMGhjQdlxKyf6CGv67K5/vDhXRojcyIsKH28dG8t/d+R0SyWs3j2BwiCfPfHOE/fk1jI/148Ubh3b7THlFUZRTWTKZ7AXihBAxmJLIPOC29hcIIfyBKimlEXga08iuk/f6CCECpJTlwFQgxYKxWtXqg0W8vTGLW5MjuHdCNAD1Wh3fHCjiv7vyOFZSj7uTPTeNCicxwocnlh3iYEFN2/1v3prIVQlBLNqYyRPLD+Ht6sjrt4xg9sjuX7dLURSlMxZLJlJKvRDiIWAtpqHBH0sp04QQLwApUspVwGTgZSGExNTN9aD5XoMQ4nHgJ2F6N9wH/NtSsVrT4cJaHv/6IKOjffnbjUM5oqljyZ48vj1QRFOrgYQQL16aPZTBwV7c9N4OluzOb7t30fyRzBwewubj5Vzz5lYKq5u5JTmcp6+Nx1cV2BVF6UGik40Ue6Xk5GSZktK7Gi9ldVpueGc7rQYj930hkD8AACAASURBVF8ew9q0Eg4V1uLiaMcNI0K5bWwUHs72XPn61g73vXtbEjOGBVNW38ILq9P5/nAxsQHu/H32MMb2VwV2RVG6TgixT0qZfLGPY+0CfJ+l1RlY8Pm+tqVN/rn2OAODPHh+ZgKzk8KpbGhh6mtbOtzzr9uTmDEsBINR8sWuPF5dc5wWg5HHrhrIgitUgV1RFOtRycRK3tuczYGCGpwc7LhuWAi3jY0kOcoXIQRVja0dEsl7tydx7bAQimubeW9zNstTC8kqa2DCAD9evHEYMf7uVnwliqIoKplYzdTBgfh5ODFzeOhp9Q1vV0d+dXkMIyN9uWJQAGuOlHD7R7vYkV2JlJAU6cNb8xK5YUSoKrArimITVDKxkhERPmfcs93eTjBlcCBfpxTw+NcHadYZiOznxsNT45g9Moxo1RJRFMXGqGRig74/VMyDS1LxdnVkTlIYc5LCSIr0Va0QRVFslkomNqapVc+L36eTEOLFygfHq6K6oii9gkomNuZfm7IprtWyaP5IlUgURek1zr4krdKj8iob+XBrDjcmhjL6InZVVBRF6WkqmdiQv313FEd7wdMz4q0diqIoynlRycRGbD5exoajpfx+WhxBXi7WDkdRzqm2Scen23O5/aNdpBfVnfsG5ZKmaiY2oFVv5IXV6cT4u7ct9KgotkhKSUpeNV/uzuf7w8W06I0IYRrOrvRtKpnYgE+255JT0cgn945WRfc+rrFFj4ujvc29OVc3trI8tZClewvIKmvAw9mBkZE+7MqpYv6YSAYFe1o7RMXKVDKxstI6LYt+yuTK+ECmDAq0djhKD6lubCWrvIHM0gayyhrILKsnu6yBolotc0aG8fqtidYOESklu3KqWLr3lx09R0b68OpNw7lueAi/XpyCt6sjj189yNqhKjZAJRMre+XHY+gMkmeuS7B2KEo3k1JSWtdCVlkDWWX1ZJY1mL9voLKxte06V0d7XJ3sqWpsxdPFgdvHWXfX0IqGFpbvK+R/ewvIqWjE08WB+WMimDcmkvgQLwB+PFzMjuxKXpg1hH5quwMFlUysKuVEFSv3a3hwSmyXl0iRUvLB1hyGhHoxMS7g3DecIjW/GgEkRvioGfXdxGiUFFY3k1Ve366l0UB2WQP1Lfq267xcHIgL8uTK+CDigjyIDfQgLtCDklot93+WQqCnM5/dN6btDbunX8OO7Eq+3JvPurQSdAZJcpQvD04ZwIxhIbg6/dL92txq4MXvjzI42JPbxnRMfI0tejYfL2fCAD983FSS6UtUMrESg1Hy3Ko0QrxdeHDKgHNebzRKnvn2SIfNsU68ct15P+/Tyw9zvLSeGH93bkwMY/bIMCL93M77cfoincFIXmVjh4SRVdZAdnkDLXpj23UBns4MCPBgdlIYAwI92r4CPJxPS+A/HS3lwSWphHi7svi+MUT069nfRVm9lq9TTK2Q/KomvF0duXNcNPPHRBAX1Hkd5L0t2WhqmvnfgnE42HccEPrUisOsPliEk70dUwYHMHtkOFMGB6haYB+gkomVLN2bT1pRHW/PH4mb05l/DXqDkUe/Osiqg0Vtx4aEerF0wbgLet5JA/05XlqPu7M9b/6UwRsbMhgd7cvskeFcNywEbzfHC3rcS0lzq4HsclOSaF/TyKtsQm/8ZTO5MB9X4oI8GB/rx4BAD+KCPBgQ4Nnlv8OvUwp4asVhEkK8+OTe0fh7OFvqJXVgMEq2ZZazdE8BG46WojdKxsb047GrB3LNkGBcHM/8xl9Q1cT7W7KZOSL0tI3YVh8sYvXBIu4ZH42DneDbg0WsTSvF29WR64eHqDXmLnFqp0UrqGlqZcrCzQwM8mTpgnGd/udq0Rv4zef72Hy8vO3Y2Jh+fHrvmA5dDufrUGENN7yznVdvGs7lcf58c0DDilQNWWUNZ7xn2mDTwAB7O8GjVw9kcHDPd8NcrOMl9XxzQEO0nxu3jj69JtHYoufJ5Yc4WFhDYXUzp/63CPNxJTnal+TofiSG+xAb6H7WDwHn8sGWbF7+8RiXD/Dn/TtH4eFs+c91JbVavk4pYOneAjQ1zfRzd2LuqHBuHR1BbIBHlx7jN5+nsDWjgo2PX0GIt2vb8dI6Lde8uZUoP3eW//YyHOzt0BuM/JxVwcr9GtamlaDVGYnyc2P2SFOLOMpPrX5tC9ROi73Ya+syqG3W8fwNQ05LJE2teu76zx5S8qrbjl2VEMS7tyXh5HDxc0yHhXkT5efG6kNF3DI6gt9NHsADV8RyRFPHbf/e1aGP/6SfjpUB4O/hTFOr4aJj6Cll9VpWHShiRaqG9GLTpLoHJsd2eq3eKKlsaMVoBDdHexpPeZ2ammY0B5r59kARjvYCHzcnfFwd8XVzwtvNEV+39t874evmiLerE77ujvi4OuHj5oiLoz1Go+TlH4/y7225XD88hNdvSeyW3+uZGIySLRllLNldwMZjpRglTBjgx9MzBnNVQtB5dT9tyyxnbVop/3fNoA6JRErJk8sPodUZeP2WEW1dXw72dkweFMjkQYE0tOhZc6SElfsLeeunTN7ckMmoKF9mjwzj+uEhqr5yCVAtkx6WXlTH9W9v485xUfx11tC247XNOm5+fwcZpb+0EOYkhfHPuSO6fc7BP9ce473N2ez585VtXStanYF/bc5m0U+Znd4zYYAfL904zOb3UmluNbAuvYQVqRq2ZZZzslfqsv5+/N/0QSRF+nbpcVr1RmqaW6lt0lHdpKO66eT3rVQ36ahtbqW6UUdNcys17Y63tqudnMrRXqAz/PL/bfqQYPp5/JKUfNwc8TEnIp+TP7s6nlaX6ApNTTNf7S3gq5QCimu1+Hs4cXNyBLcmR1zQ77BVb+Tat7aiN0rW/XFShyT05Z58nl5xmOdnJnDPhJhzPlZxbTPf7C9i5f5CMkobVH3FyrqrZaKSSQ+SUnLrB7vILKtn8+NT8HZzpLKhhZlv/0xRrbbtunvGR/Ps9QnYWWji2rGSOqa/uY1HrxrIv7fmdNoauWd8NJMG+lNW18I3BzTsyqkCoH1DysfVkQ2PXoFfD/X1n4nRKNmVU8nyVA1rjhR3aFUMC/PmiemDuHyAf4/01Te3GqhuMiWYGnOCqWlupbhGyzubstquGxXlS83J65p1GIyd/z+0E/DeHaO4ZkjwOZ9bbzCy8VgZX+7JZ0tGORKYGBfA/NERTIsPuqgW0Efbcnjx+6P85+5kpsUHtR3//lAxj399kKQoHz6/b+x5/ZuVUpJWVMfK/Rq+PVBERUMLPm6m+srskeEkRaoRhz2hV3RzCSGmA28B9sBHUspXTjkfBXwMBABVwB1SykLzOQNw2HxpvpTyBkvG2hNWHSxiz4kq/j57GE06PROe20hDuzfy308dwKNXDbT4f6CTXRSvr8847dxLs4cyf3RkhzeFeWMi0dQ0s2JfIR/9nEttsw57O8GsxDCrdk9klNazIlXDtwc0FLdLxgD9/d15/JpBXDs0uEffkFyd7HF1ciXU55duoOrGVu77bC92Al6aPYz5pwynlVJS36Jva/mU1bXw9x+PklPeSFygJ8PCvM/6nAVVTfzP3Aopq28h0NOZB6cM4JbkiG4ZHVZWr+XNDZlMGRTQlkiaWvW8sDqdpXsLGBnpwxu3JJ73hx8hBEPDvBka5s3T1w5mW1YFK1M1LNtXyBe78on2c+NGVV/pNSzWMhFC2AMZwFVAIbAXmC+lTG93zdfAd1LKz4QQU4F7pZR3ms81SCm7VhXE9lsmjS16pr22hcZWPfXaji2Bp68dzG+u6Lwvv7vUNulYl17C94eLOxT1Afw9nLhpVDhf7MyjsdXA1QlB/H5qHMPCTW9ieoORr1IKeWNDBuX1LUwfEsz/TR/U5aJtdzpZB1m5X0NaUR32doJ+7k6U17cAEOLtwiNXxnFTUvgFdQ91t6KaZu76eA/5VU0smjeS6UPP3sLIrWjk4S/3c1hTy53jovjzdfGdjq7SGYxsSC9lyZ58fs6qQACTBwUyb3QEUwcHthXAa5t1lNRpOVhQy6BgT0ZFda2br73HvjrIqoMa1v3xCmL83UkvquP3X6aSU9HI7ybH8siVA3Hsxr/req3OXF/RsDOnEikhOcqX2UlhXD8sVI047GY2380lhLgMeF5KeY3556cBpJQvt7smDZgupSwQpo+PtVJKL/O5SyaZ5JQ3MP3NbbQaOvanvzR7KLePjbLY89Y0tbIuvZQfDhfzc2ZFh2Gt7b1/RxLTh4ZQ09TKJ9tP8Mn2XOq0eq4YGEB8iBfr00vILm9kVJQvf5oxmFFRPbvXysk6yMr9GrZlVmAwSoaFeZMY4UNGaT27c6vwdXPkwSkDuGNc1FmHtvakzNJ67vp4Dw1aPf++O5lxpwylPdWK1EL+8s0RHOzteHXu8LauLaPR1HKpaWrlQEEN72/J4Whxx1V6EyN8qGhoobC6+azPkfvyjPNqqe3Lq+am93bw2ytieXL6IBbvzOOlH47i4+rIG7cmMmGAf5cf60IU1TTzzQENK1M1ZJaZ6itTBwcyOymMKYMCLTp4oa/oDclkLqZE8Svzz3cCY6WUD7W7ZgmwW0r5lhBiDrAc8JdSVgoh9MABQA+8IqX8ppPnWAAsAIiMjByVl5dnkddyoaSUPPttGp/v6hjXW/MSmZUYZpHnrG5sNbdAStiRdXoCiejnyh+vHMisxDBmvfszRzR1bcnkpHqtjse+Osi69NK2Y/eMj+a5mQk91mV0sg6yYr+GNUdKaGjRE+rtwo0jw7gs1o/vDxXz9b5CnB3s+NXE/vx6YgyeLrbziXVfXjX3f7YXR3s7Prt3DAmhXkgpadYZTAX9xlZqm03dWprqZl7+8VjbvfEhXrg62lHTrKOmSUdVu6VXLkRsgDuJEb7cNjbyvFomRqNk1rvbKavXsvyB8Ty/Kp0NR0uZOjiQf84d3qO1spP1lRWpGlYd1FDR0KrqK92kV9RMuuBx4B0hxD3AVkADnKyeRkkpNUKI/sBGIcRhKWV2+5ullB8CH4KpZdJzYZ+b3mBkwJ9/7HDso7uSuTIh6Ax3/KJeq+PzXXncfVk07l2Yf1DV2MratBJ+MK+XZDBKIvq5MirKl6PFddRp9QR7ufD7aaZ+9JNdEjOHh3JE0/ETbm5FI/9ce4x16aUdRh99uuMEhwpreOPWRIv2X2eW1rNiv4Zv92soqtXi4ezAjGHBzB4ZzsAgD97fks39n6WAhLsui+LBKQN6bLJfV+3OqeTuT/a0zav44/8OtBXlT22dnirYywU7Yaq9hPi48v2h4rNeHx/iRWyAOzH+7kT7uRMT4E6Mnzu+3bBe1lcpBRzW1DJ/TCRz39tJVWMrz16fwL0Tonv8jbt9feVPM0z1lRWpGr5O+aW+Mi0+iBBvF4K9XQjyciHYy4VAL2c1OqyHWDKZaICIdj+Hm4+1kVIWAXMAhBAewE1SyhrzOY35zxwhxGZgJNAhmdiyRRt/Gbnzxf1juTyua90BUkoe//oga9NKGRLqzRUDO19/q6KhpS2B7MqpwmCURPm5sWBSf/r7u7M2rYQNR8vwc3fimeviO+3+mZUYxsfbc6lp0vHvrTm89MPRDufbD2MFOFZSz4nKpm5PJuX1Law6aBoqekRjqoNMivPnqRnxXBUfhFFK/vNzLgsWp9DYqmdOUjiPXBlHuK9tLgOTVd6Aj6sT3v0c2+agjHTzaZuD4uPqiI+bI8v2FbLhaBkOdoLP7hvTaZdRdeMuqhpbifZzJ9rfnf7+pj+j/d06XZ6lu9Q26Xh17XHAtFpDjJ87H909nqHnGAzQExzs7ZgyyLTKdr1Wx49HSliZquHzXXmdDs32dXM0JRdvF4I8XQjydiHIy5lgL1PSCfJywc/dyWKjJ/sKS3ZzOWAqwE/DlET2ArdJKdPaXeMPVEkpjUKIlwCDlPJZIYQv0CSlbDFfsxOY1b54fypbqpkc0dRy47vbmTEshEXzR57Xvf/anMWra0z/iT+9dzST2y1LX17fwpq0En44VMzu3EqMEmL83ZkxLJgZw0JwcbTnjfUZfHeoGE8XB34zqT+3jo6kqrGV/Kom8quaKDB/5Vc1UVDdhFZ3+n++/v7uJEX5EtnPjch+bkT0cyOin2u3vnlpdQbWpZeyIrWwQx1k9sgwZo4IJcDTmRa9gSW783lnYxaVja1cMySIx68edMY1o3qL8voWHvv6IFszyrlmSBD/uGm4zU3ae35VGp/uOAHALcnhPDdzSJdaydYkpWwbcFBSq6WsrsX0fZ2WMvOfpXUtVDS0nLbCgaO9INDT1JJpn2SCvZ3bElCwl4vN/x1cCJvv5pJS6oUQDwFrMQ0N/lhKmSaEeAFIkVKuAiYDLwshJKZurgfNt8cDHwghjJi2Fn7lbImkO32xK48Ptl54A6hFZ6TMPLJo0/EyJr66se1cgIczSxdcdsai4c+ZFSxce5y4QA8yzcublNVrWXPE1ALZk1uFUUL/AHcenDKAa4eGMDjYk9T8an7331TyKpvaHivQ05lPd+SxcF3H4b/uTvZE+rkT4uNKbkVj2/EBgR4smjeShFDLLZViNEp25VayMlXDj+3qIAsm9WfOyLC2JGEwSpbtK+SN9Rloapq5rL8fT0wfxMguTji0ZVsyynnsqwPUa/W8eONQbh8baXN9/cW1zXy+Kw9PZwdenD3UYvW97iaEeWUCN6ezLvmjMxgpr2+htE5r/jIlndJaLaX1WjJK69mWWdFh2P5JHs4OplZN+1aOp3Nb11qQlwsBns7dOrqtt1CTFk+xLbOcz3acoLC6GU1N82nDeDtjbycI9nIhzNeVPbmmyX3+Hk5Migtgd24VmhrTCJuZI0J589bETme0a2qamfn2z/h7OPHbK2J59KuDgGmSoJSmyWtxgZ4MDfPG1cmOgqpm9uVVd/oPPtzX1dSi8HUj0s/UsjjZynC0F3yxK5+PtuVQ2djKmJh+PHHNIJKjLTdCK6vs5HyQIjQ1zXg4O3Dt0GBmJ4UxLsavrXtBSsm69FIWrj1OZllDj084tKRWvZGF647z4dYcBgZ58Pb8JJvdnbBeq+P9LdncmhzZp1eUbmjRt0s4Wkpqf0lAptaO6edTB7kIYVp66GRXWqC5fnOyhhNsbuV4uzraxL9rmx/N1dMs1c1Vp9WhqW42fdWYv6qbKTT/WdHQ0ul9Id4uHSbSLZjUn/smxBDs7XLatVqdgakLN1NUq8Xfw/mMj3k2r940nLH9+xHq49rpp6LaJh2f7jjBx9tNkw4nDQzgoSkDGBNjmSRS2dDCt+b5IIc1tW11kNlJ4VwVH3TaYpU7sit4dc1xDhTU0D/Ancev7vkJh5ZyoqKRh5fu51BhLbePjeQv1yfYzPBl5eIYjZKqplZKaju2ck52q5XUaimrb+l0RJ6zg525NePcNmAgyOuXLrWTxy39b0Ulk1NYq2ai1RkoMieZO/+z56zXOtoLDj539WmrzT61/BBL9xac1/MKAbNGhPLIlQPPutZSZUML//k5l8U782ho0XNlfBC/nzqAERE+5/V85+P7Q8U8teIQ9Vo9Q8O8mD0ynBvMdZBTHS6s5dW1x9iWWWFzEw67wzf7Nfx55WHs7QSvzh3eYQi20ne06A1tLZmTtZtSc7Jp39rprIbp7epoSi6ddKmdTDp+Hs4XvIafzddM+goXR3v6B3gQ4+/OQ1MGdFh/aftTUwnydKakTsujXx0kr7IR11M+ZSzdk8/SvQVcGR/EfROiiejnRrC3Cz9nVdBwli62+BBPBgSeuZukrE7LB1tzWLI7H63ewIxhITw4eYBFayLNrQZe+C6NL/cUkBjhw8tzhp1x18Ds8gZeX5fB94eL8XVzPOOIs96qoUXPs98eYUWqhtHRvrw5byRh7ZZYUfoWZwd780CWM3cbSimp0+rbdaudUtOp03K8pI7y+hZOnX9sbycI9HQ2d6md0r3m/Usrx8PZwWKtfZVMukFuRSN/XZ3G5uPlBHu5oDMYqWxs5cllh3jvjiRCvF3JKK3nyvigDr/IgwU1PPttGhPj/PngzlEdPllMaTeK63y06A289P1Rlu4twGCUzBoRyu+mDGBAoGWXPjlaXMfvv9xPdnkDD0yO5dGrOl9io7i2mbc2ZPL1vkJcHOx4eFqcTUw41JkXSZw86OJXrT1QUMMjS/eTX9XEH6bF8fupAy6ZlpZiOUIIvF0d8XZ1ZOBZRizqze8vJbXaU7rUWiir15JT3sjO7ErqOvkw6uZk3260mjNBnXS7XyiVTC5CU6uedzdl8e+tuTg52PHMdfHcMjqCL3fn89q6DH7OquBYST1O9nbUNOmY2G6uSWVDCw98sY8AT2cWzRvZbcvMf3ugiMU787h5VDgPTR1g8QXypJR8viuPF78/irerI5/f1/mcmqrGVt7bnMVnO/NscsLhR9ty+ceaY6z/46QLHnpcr9WxcO1xFu/KI9jLhSW/HnfOJVQU5Xw52Nu1dXONOMt1Ta36tu60zgYRpORVU1Z3/vXZM8bVbY/Uh0gp+fFICS9+l05RrZY5I8P441UDWZ9eytSFW6hoaGFinD+PXDmQUVG+vP1TJkLA5eZJaXqDkYeX7qeisZUVD4zvltnKJx0rrsfF0Y5/3DTc4pOwqhtbeWL5IdanlzJ5UAALbx5xWnJobNHzn59z+XBrDk02OuGwtE7L2xszGRLqdUEtOCkla46U8PzqNMrqW7hrXBSPXTMILxta3kXpe9ycHIjxdyDmLDVVKSV2f++e51PJ5DxlldXz/Kp0fs6qID7Ei4U3jyC7vIG57++gtK6Fcf378a/bkzqMktqWWcHQUO+2tYwWrstge1Ylr84d3u0zijNK6xkY5GnxRLIrp5JHlh6gsrGFZ66L574JMR2eszdNOHz5h6M0tRr4zRWx592fXFDVxHOr0th4rIz4EC8+uDOZRAsOblCU7tSd9ROVTLqooUXPop8y+fjnXNyc7PnL9Qk4Odjxf8sOoalpJjnKlzduTWR8bMcunnqtjtT8ahZM6g/AmiPFvL8lm9vGRnJLckRnT3VRjpfWn3EJlu6gNxhZ9FMmb2/KItrPnZV3T+iQEA1Gycr9ml4z4XDviSq+OVBEuK8rM86xPHx7OoORT7bn8sZ6086Uz1wXzz3jo1VtROmzVDI5Byklqw4W8dL3Rymrb+GmpHDigjz4dEcuBVXNjDCPWpoY1/nEup3ZleiNkolxAWSVNfDYVwdJjPDhuZkJ3R5rVWMr5fUtDLLQp//C6iYeWXqAlLxqbkoK54VZvyyx0dmEw1duGmbTEw4NRslz35pW97n/8pguJ4L9+dX8aeURjhbXcWV8IM/fMMSmuu0UxRpUMjmLYyV1PPttGntyqxgS6sWMYSFsPl7G8tRChoZ58dd7hjBlUOBZ3yy3ZVbg5mTP4GBP5r6/AxdHe967I8kiK5lmlNYDMPAiZ1ZrdQbu/2wv27MquWNcJC/eOIwfDxfz5PJDGCW8eWsiN47suMTGx9tP8LfvTCve3Dshmr9cZ7lth7vLl3vySS+uw9vVsUutxDqtjn+uOc4Xu/MI8nTh/TtGcc2QIJtNlorSk1Qy6URts4431mfw+a48PJwdmDDAj+IaLZ/uOMHgYE8+uHMUVyd07U1ka2Y54/r78aeVh8mtaOSLX41t2za3u51MJhfaMtHqDPx6cQrbMivajiVG+PKnlYdZsjufEeHeLJo/stMRYlfGB5KmqWVNWgmfbD/B+vRSZpu3XO1vhR0Zz6W6sZWF60wLat45LuqsC/hJKfnhcAl/XZ1GeUMLd18WzWNXD7T6cGZFsSUqmbRjNEpW7Nfwyo9HqWhoxd/DGXs72J5VSVygB+/elsS1Q4O7/Ik7r7KRvMompIT8qib+NGPwaTWVrmhq1aOpbj5n8fp4ST1eLqaF6LqiudXAkaJaBgd78uCS/WzN+GU737mjwrlvQgyP/G8/GaUN/OaK/jx21aAzLlIZ5efO67cm8mKrnrVpJaxI1fDupize3phFYoQPc5LCuH54KP26ceTaxXht/XFqmnQ4Odhx9/joM15XUNXEX749wubj5QwN8+Kju5MZHq4K7IpyKpVMzI5oann22yOk5tcApuVKKhpa6O/vzlvz4rl+eOh5zwXZav6En1/VxHXDQvj1xP4XFNvyfYU8vzqdNX+YeNaEklFaz6Bgz7O2mAwndzBM1bA8tfC083NGhvHq3OF8ubeA2f/ajqeLI4vvG8OkLhb13ZwcmD0ynNkjwymt0/LtAQ0rUjU8+20aL6xOZ/KgQOYkhTF1cKDVZrunFdWyZHc+ADclhXe6zIvOYOQ/P+fy5oYM7ITgL9cncPdlUarArihn0OeTSU2Tqbvjv7vzO+xxEOHrxh+mxTErMfSC30C2mT/pDwj04B9zh19w3/rgEC8MRsk/1hzjo7tHd3qNlJLjJfXMHBHa6fnjJfWs2F/It/uLKKnTnnZ+VmIor84dTnOrgQeXpLI2rZRJAwN47eYRnb7ZdkWQlwsLJsWyYFIsR4vrWLlfwzf7NWw4WoqXiwPXDQ9lTlIYyVG+PVZ3kFLy11XpGKXpA8OvJ8acds2+vGr+vPIwx0rquSohiL/eMIRQtRSKopxVn00mBqPkq5QCXl1zjOomXdvxMB9XHp42gDlJ4Re1J4HeYGRndiUezg58cOcoPC5iU53ECB88nB3YcLSMXTmVnc6qLq1roU6r77CseVm9llUHiliRqiG9uO60ewCuHx7CP+eOwNXJnt05lTzyvwNUNLTw5xnx3H95TLcV0eNDvIgP8eLJ6YPZkW3acvWb/Rq+3JNPRD9XZieGMTsp/KwTrLrDqoNF7Dlh2ibgqvigDvWc2mYdr645xpI9+QR7ufDBnaO4ZkjXhwsrSl/WJ5PJgYIanv32CIcKa9uOBXu58NBU0x7pZ6oLnI+Te1b/amIMsRdZgHa0t2N8rB/r0kv5+w9H+eZ3E057kz9uLr5H9HNr61rallmOUdJpy2LqmmArIAAAEy5JREFU4EDenJeIl4sjeoORN9Zn8PbGTCL7ubH8gfEWqwvY2wkmxgUwMS6AF2801VdW7tfwzqYsFm3MYmSkD3NGmuor3bkyAJhm4/+93dbEv7nC1O0opeS7Q8X8dXU6VY0t3Ds+hkevHnhRHwAUpa/pU0vQVza08Oqa4/wv5Zfl3gM8nXlwcizzxkTa9Iq1n+/K+//27jw8qvpc4Pj3zUYIEAgQKCSQBAEBBQmgbCqL4ga1ilbFWlDautXb1T7V61JLrxcseilivbe0tbR9rG1dEBcEZZNWUFlkCSELxAAJCQQCZoHs7/3jHOIkJCEyc5IJvJ/nmSdnzvll5p0zk3lzzu/83h9PvJkCwPMzkrnR53RWTY3ynT9vYm36lx3ocV3aMzIhhh05x8n2mYERYOsTU2o7wg8eP8mP/r6NT7MLmZ4cx5ybLm6VL1Hf/pW0/GLCQ8XpX0mOY/LgHgG5lPqZFWn87zpnFs1RCTG89sA49h89wePLUlifUcDQuM7MnT40KOY5N6al2Hwm9TSVTKprlJc/2cezK9NrK2l26xDBAxMvaDNlz/cdLWXC/HWAM5Pi6p9OYN/RE+4Mhrm1E3HdNiqe4X1ieD81n3U+yQXgmiE92V94os7kW0VlVYSHCL+66WKmj4hvsdfTlNSDRSz9LIc3tx2koLic6Mgwpl3Sm+nJcYw8y/6Vz4+Ucs2CDwkRobyqhhe/NYLso6UsXJVJWIjw8LUXMnNsYsAKbhrTVlgyqaexZLI5u5Anl+2q7TPoEhXOfVdewKxxCadNUhXsJsxfW2eed6B2BsO16QV0ahdGckJMnUt86xuVEFOnXyUiLISZYxM976s4G1XVNXy09yhLt+awctchTlZW07drFDclxzE9Oa7JScHqu+dPn7Ip+xgd2oVyrLSShG5RZB4u4dqLevLUjRd5NvbHmGBnk2OdweHiMua9l8YbW3MBiI4M494r+zFrXGKbG2xWVV3DuzvzTkskP7xqAHeNSSD3+EnWphdQXF7VYCJJ7BblXq4b16bm9A4LDWHCwFgmDIyltLyKFSlO/8qiNZk8vzqTEX27cPOIeKYN7dVk/8qatEOsTS9g0oWxtacCS8ur+P3MUUwZ0rOlXo4x5zRPj0xE5DpgIRAK/EFV59XbngC8BMQChcBdqprjsz0aSAXeVNWHmnquU0cmldU1/HlDNr9ZlUlJeRWd2oUx+/IkZl+eROf2bSuJnLJsWy4//Pu22vt3jenLy5/sJ7lPF9qFhbIx62iDv/f1S3pzz/hEkvt0OadKfuR/8WX/Svohp39lkjt+ZdKguv0r5VXVXLtgPaEhwph+3Xjl0/3MHp/Ej6cMbHLUuzHni6A/zSUioUAGMAXIATYBM1Q11afNq8A7qvpnEZkM3KOq3/bZvhA30TQnmSz6x0p+8VYKGYdKiIoI5Z7xiXzvin50iQqOUddnq6S8ikde38E7O/IAp3M99/jJRttPGBjL4pkjPan/FUxUldS8IpZuzWXZdqd/pXP7cKYN68X0EXGM6BvDi+v2Mn9lOn+ZfRlD4zpTXFbVpo7OjPFaW0gmY4GnVPVa9/6jAKo616fNLuA6VT0gzr/OX6hqtLttJPAzYAUw6kzJpGe/Idr+tvlEhocwa2wi917Zr3b+kHOBqvLgy1t5LyW/0TYd24Xx+gPj6vSJnC98+1dW7MqnrLKGhG5RHC5yJipbPNPvvxVjzkltoc8kDjjgcz8HGF2vzXZgOs6psJuBTiLSDTgGPAfcBVzd2BOIyL3AvQCRX7uA749P4v6J/ejRKXDzGre2mhqntPsLazNJyW144OEpmx67mvYR5/bRSGN8+1dKavtXcjhZUc0T0wJf7t8YU1drnzR+GHhBRO4G1gO5QDXwILBcVXOaOtevqouBxeCc5nrSgzlCWkt1jfLOjoP8du0eMg6VANC5fTh9u0axM/cLJg/qwZq0w7XtrxrU47xNJPV1bBfGrSPjuXVkcFzqbMz5wMtkkgv4ThIR766rpaoHcY5MEJGOwC2qetw9RXaFiDwIdAQiRKREVR/xMF6/LVyVyZq0Q81uP2FgLD+55sLT1heVVXLLixvIPFxCeKhw7UU9uTk5nkmDYjlWWsmYuatZk3aYSxNjSO4bw+L1WX7PYWKMMf7wMplsAgaISBJOErkDuNO3gYh0x+lcrwEexbmyC1X9lk+bu3H6TII6kWzZV8iCVRkMjetMt45Nd/gfKDzB3oJShvRueKR1RGgIIxNimDUukWnDetVeQPDKp/t59I2dAPTo1I5XvjeGGoWKqhpuGh7X4GMZY0xL8CyZqGqViDwErMS5NPglVd0lInOAzar6FjARmCsiinOa6/texeMlVeXpd3fTo1M7/nHfmCYHQ56sqGbqon/Ru3Mkj1w/qME2keGhzLtlWO396hpl/Lw1tdV+fzplIP9x1YDa7U/deFGAXokxxpwdT/tMVHU5sLzeuid9ll8DXjvDYywBlngQXsCsSMln6/7jzJs+9Iyj6p9ZkUZWQSkvf3d0s8a9LPnoc556u/ZqatY9PPErjfw2xpiW0Nod8G1eZXUNz6xIY0CPjmfs8P1ozxGWbMjm7nGJjO/f9IyLJeVVXPyLlXXWfT73hnNq8KEx5txh08b56W+f7Cf76AkevWFQk5NoFZVV8rNXt9MvtgM/v67h01vgnDJbti23TiKZcVkfsudNtURijAladmTih6KyShauzmRsv25MurBHk21/+VYqh4rLef2BcY1ewptVUMLsJZvqlIzf8Mhkm+XPGBP0LJn44f/W7aWwtIL/vGFwk0cNK1LyeX1rDj+Y3J/hfU6fdOpkRTUvrtvDojV7atctuP0Sbk62cRLGmLbBkslZyvviJH/89+d8Y3hvhsY3PpnSkZJyHlu6k4vjonlo8oDTtq9KPcTPX9/B0dIKABK6RfHyd0cTH2P1o4wxbYclk7P03PsZqMLDDQw6PEVVefSNnRSXV/HKbcPrTAd8oPAEv3w7lVW7vxzk+IPJ/fnBVQOa7HsxxphgZMnkLKQeLOL1rTl89/Ik+nRt/AjitS05fJB6iMenDmZgT2eEenlVNb9fn8WiNXsor6oBoHvHCBbNGMHYC7q1SPzGGBNolkzOwrwVaURHhvPQpNNPW52Sc+wEc95OZXRSV2aPTwLg35lHeHJZCllHSmvbXT24J/NvHdbk5E7GGBPsLJl8Rf/KLGB9RgGPTx1M56iGBx3W1Cg/e3UHNao8+81LOFxczq/eTeVddz4ScKbLfXzqYL49JsEu+TXGtHmWTL6Ciqoa/nt5GvEx7fn22IRG2y3ZkM3GrKP8100Xs3JXPgs+yKDMPaUF0L9HRxbNSGZwr+iWCNsYYzxnyeQMVJVtB46z9LNc3t5+kGMnKnl+RnKjsxjuOVzMMyvS6N6xHX/duI/0Q8X079GRY6UVHC2tYMZlfXhi2pAzll0xxpi2xL7RGnGg8ARvfpbL0s9yyTpSSkRYCFOG9OSbI+OZ2MgAxcrqGn7yz+2UV9VQXlJORKgwbVgvPkwvAIEX7kxm2rDeLfxKjDHGe5ZM6tmUXcj8Fel8ml0IwOikrtw3oR/XD+1FdGTThRl/9+FeduR8QXio8K3RCRwpKeedHXmM6NuFhXckN3nllzHGtGWWTOrZsOdobSIBSMsv5o2tuezOK2ZIr2gG9erEwJ6diAw//TRX3hdlXDkwlltHxvObVRl8fqSUhyb150dX29gRY8y5TVS1tWMIiFGjRunmzZsD8lhFZZWk5xeTlldEal4xaflFpOcXc6KiGoAQgaTuHRjcK9q9dWJwr2i+Fh3Jkg3ZzF2eRkyHcBbcPpxxFzRdHdgYY1qTiGxR1VH+Po4dmTQgOjKcSxO7cmli19p1NTXKvsITpOUVsTuviN35xWw7cJx3fC73jYoI5URFNVcN6sH8b15CVxs7Yow5T1gyaaaQECGpeweSunfg+qG9atf7HsWk5RdzcVxn7ri0j40dMcacVyyZ+KmhoxhjjDnfWK+wMcYYv1kyMcYY4zdLJsYYY/zmaTIRketEJF1E9ojIIw1sTxCR1SKyQ0TWiUi8z/qtIrJNRHaJyP1exmmMMcY/niUTEQkFfgtcDwwBZojIkHrNngX+oqrDgDnAXHd9HjBWVYcDo4FHRMTqkBhjTJDy8sjkMmCPqmapagXwd+Ab9doMAda4y2tPbVfVClUtd9e38zhOY4wxfvLySzoOOOBzP8dd52s7MN1dvhnoJCLdAESkj4jscB/jGVU9WP8JROReEdksIpsLCgoC/gKMMcY0T2v/x/8wMEFEPgMmALlANYCqHnBPf/UHZolIz/q/rKqLVXWUqo6KjY1tybiNMcb48HLQYi7Qx+d+vLuulnu0MR1ARDoCt6jq8fptRCQFuAJ4rbEn27JlyxER2ddEPN2BI1/pFbQ8i9F/wR4fWIyBYjEGxoWBeBAvk8kmYICIJOEkkTuAO30biEh3oFBVa4BHgZfc9fHAUVU9KSIxwOXAgqaeTFWbPDQRkc2BKGbmJYvRf8EeH1iMgWIxBoaIBKRCrmenuVS1CngIWAnsBv6pqrtEZI6I3Og2mwiki0gG0BN42l0/GPhERLYDHwLPqupOr2I1xhjjH09rc6nqcmB5vXVP+iy/RgOnrlT1A2CYl7EZY4wJnNbugG9Ji1s7gGawGP0X7PGBxRgoFmNgBCTGc2ZyLGOMMa3nfDoyMcYY4xFLJsYYY/zW5pOJH8Ukh4vIRreQ5A4RuT0IY2yxgpdnG6PP9mgRyRGRF4IxRhGpdvfjNhF5K0hj7Csi74vIbhFJFZHEYIpRRCb57MNtIlImIjcFU4zutl+7fy+7ReR5EW+mPfUzxmdEJMW9efLdIyIvichhccbpNbRd3P2zx41xhM+2WSKS6d5mNesJVbXN3oBQYC/QD4jAKc8ypF6bV4FZ7vJk4K/u8kBggLvcG6e4ZJcgizECaOcudwSygd7BFKPP9oXA34AXgu29du+XBPPn0b2/Dpji835HBVuMPm26AoXBFiMwDvjIfYxQYCMwMchinAp8gHM1bQecMXnRHsR4JTACSGlk+w3Ae4AAY4BPfN7bLPdnjLscc6bna+tHJv4Uk8xQ1Ux3+SBwGPCiJktbKHh51jECiMhInHFC73sUn98xtpCzjlGcitph6lwWj6qWqOqJYIqxnluB94IwRgUicf8RA8KBQ0EW4xBgvapWqWopsAO4LtABqup6nITfmG/gVG1XVf0Y6CIivYBrgQ9UtVBVj+EkvjPG19aTiV/FJE8RkctwPnx7gy1GaUbBy9aMUURCgOdw6qx5yd/3OlKcoqAfe3Vqxs8YBwLHReQNEflMROaLM41DMMXo6w7gFQ/iAz9iVNWNOF/cee5tparuDqYY3fXXiUiUOFVAJlG39FRLaew1NOe1naatJ5PmaLSYJICbif8K3KNOWZfW4FfBy1aO8UFguarmtFJcvpp6rxPUKWtxJ/AbEbkgyGIMw6k/9zBwKc7pk7uDLEag9m9mKE51i9bSYIwi0h+ngkY8zhfgZBG5IphiVNX3cQZzb8BJyBvx2b9tlacj4FuAX8UkRSQaeBd4zD3MC7oYfdtIMwpetnSMIjIWuEJEHsQ5zx8hIiWqelqHZGvF6G7LdX9micg6IJnAH4n6sx9zgG2qmuVuexPnPPYfgyVGnya3AUtVtTLAsfkdo4h8D/hYVUvcbe8BY4F/BUuM7rancctHicjfgIwAx9ccjb2GXJxSV77r153x0QLd6dOSN5xkmAUk8WUn2EX12nQHQtzlp4E57nIEsBr4URDHGA+0d5djcD5wQ4Mpxnpt7sa7Dnh/9mMMX17I0B3IpF5naRDEGOq2j3Xv/wn4fjDF6LP9Y2CSF+9zAPbj7cAq9zHC3b/xrwdZjKFAN3d5GJCC01/mxb5MpPEO+KnU7YD/1F3fFfjc/buJcZe7nvG5vPpAtNQN54qEDJz/Mh9z180BbnSXb3W/PDKAP/h8qdwFVALbfG7DgyzGKTidc9vdn/cG236s9xh341Ey8XM/jgN2uvtxJ/CdYIux3vu9E1gCRARhjIk4/7mGeLUP/XyvQ4Hf4RSXTQX+JwhjjHRjS8VJzF5977yC029UidPv8R3gfuB+d7vgTK2+1/3MjfL53dnAHvd2T3Oez8qpGGOM8dv50AFvjDHGY5ZMjDHG+M2SiTHGGL9ZMjHGGOM3SybGGGP8ZsnEmAASkXgRWeZWW80SkRdEpJ1bematiJSIh5WVjWktlkyMCRC31PkbwJuqOgAYALQHfg2UAU/gfQ0zY1qFJRNjAmcyUKaqfwJQ1Wrgx8BMnCmy/42TVIw551gyMSZwLgK2+K5Q1SKceWj6t0ZAxrQUSybGGGP8ZsnEmMBJBUb6rnArU38NSG+ViIxpIZZMjAmc1UCUiMwEcCe3eg6n+OXJVo3MGI9ZoUdjAkhE+uBUYh2MMw30P1T1PndbNhCNU7L8OHCNqqa2UqjGBJQlE2M8IiLjcMqA36yqW1s7HmO8ZMnEGGOM36zPxBhjjN8smRhjjPGbJRNjjDF+s2RijDHGb5ZMjDHG+M2SiTHGGL/9P2bHDGfHog15AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "q1q2_df.plot(x='Q1', y='Q2')\n", "plt.ylabel('Q2')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## References\n", "\n", "[1] Oliver Beckstein, Elizabeth J. Denning, Juan R. Perilla, and Thomas B. Woolf.\n", "Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of OpenClosed Transitions.\n", "Journal of Molecular Biology, 394(1):160–176, November 2009.\n", "00107.\n", "URL: https://linkinghub.elsevier.com/retrieve/pii/S0022283609011164, doi:10.1016/j.jmb.2009.09.009.\n", "\n", "[2] Richard J. Gowers, Max Linke, Jonathan Barnoud, Tyler J. E. Reddy, Manuel N. Melo, Sean L. Seyler, Jan Domański, David L. Dotson, Sébastien Buchoux, Ian M. Kenney, and Oliver Beckstein.\n", "MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics Simulations.\n", "Proceedings of the 15th Python in Science Conference, pages 98–105, 2016.\n", "00152.\n", "URL: https://conference.scipy.org/proceedings/scipy2016/oliver_beckstein.html, doi:10.25080/Majora-629e541a-00e.\n", "\n", "[3] Naveen Michaud-Agrawal, Elizabeth J. Denning, Thomas B. Woolf, and Oliver Beckstein.\n", "MDAnalysis: A toolkit for the analysis of molecular dynamics simulations.\n", "Journal of Computational Chemistry, 32(10):2319–2327, July 2011.\n", "00778.\n", "URL: http://doi.wiley.com/10.1002/jcc.21787, doi:10.1002/jcc.21787." ] } ], "metadata": { "kernelspec": { "display_name": "Python (mda0170)", "language": "python", "name": "mda0170" }, "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.6.7" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": false, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": true } }, "nbformat": 4, "nbformat_minor": 2 }