{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# [Prody](http://prody.csb.pitt.edu) is an open-source Python package for protein structural dynamics analysis with an expansive and well-written API." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "scrolled": true }, "outputs": [], "source": [ "import prody\n", "prot = prody.parsePDB('3UB5')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The default visualization in ProDy uses matplotlib." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAGJCAYAAACpeBXYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAC8z0lEQVR4nOz9d5xjV33/j7/uVe9leu87M9v77syCDbbBxgZsbGpMMISYQGxsyjc4lCS/D7EhTgg4gCkm9FCT0EI3Ni7Yu8benT6z03tT14y6dO/5/TF776rOSJqrtnufj8c+7NVKulfSved13p0ihBCIiIiIiIgIBF3oExARERERubIQhUVERERERFBEYRERERERERRRWEREREREBEUUFhERERERQRGFRUREREREUERhERERERERFFFYREREREQERRQWERERERFBEYVFRERERERQRGEREREREREUUVhERERERARFFBYREREREUERhUVERERERFBEYRERERERERRRWEREREREBEUUFhERERERQRGFRUREREREUERhERERERERFFFYREREREQERRQWERERERFBEYVFRERERERQRGEREREREREUUVhERERERARFFBYREREREUERhUVERERERFBEYRERERERERRRWEREREREBEUUFhERERERQRGFRUREREREUERhERERERERFFFYREREREQERRQWERERERFBEYVFRERERERQRGEREREREREUUVhERERERARFWugTELm6IISAZVkEAgFIpVJIJBJIJBJQFFXoUxMREREIihBCCn0SIlcHhBCEw2EwDAO/3w8AoCgKNE1DKpWKQiMicoUgCotIXmBZlhcViqIQCoV48WBZFoQQEEJEoRERuQIQhUUkpxBCwDAMIpEIWJYFTdO85RIvFtylyLIsWJYFIFo0IiKliCgsIjkj2vUFbIkERVG89cI9tt3ruT8sy/LPpWkaMpmMFxlRaEREigtRWERyAsMwCIfDvJUSvfCzLMu7wjIRhHihAS5bNDKZjLdo4o8nIiKSX0RhEREUQggikQgikQgAJBWPbIUl2bFEoRERKT7EdGMRweBcXNGLfC4X9Oj3l0gkMUITCAT450gkkpgYjSg0IiK5RbRYRHYNZzGEw2E+s2u7hVsoiyWd84q3aGiaTkgGEIVGRERYRGER2RWpAvTbkS9hiSdeZAghotCIiOQA0RUmkjXRtSmZLMaFWrQ5IaPprU5GnNAwDIMXX3wR7e3tMBgMMVlnUqk07wIoIlLqiMIikjHJalMyXXiLYaGOFhquroaiKP6zcX+Pj9GIQiMisj2isIhkRLzr60pyG3EiE2/RRCKRGOGRSCQxdTTc80VERLYQhUUkbbarTbkSSeU6E4VGRGR7RGER2ZHo2hQu4H2li0oydhIaAEnbz4hCI3K1IQqLyLawLItIJHJFur52SyqhCYfDMVlvotCIXG2IwiKSlExrU0odIT5bMqHhvsPo3mjRQsNlnYmIXEmIwiKSQDptWUR2hou/cEQLTTKLJjrrTESklBGFRSSG+LkpV5PbJte1wukIDU3TCckAotCIlBqisIgAEKY2RSQz0hUacRaNSKkhCosIH3AeGRmBXC5HW1ubuHgVgGihiR56FgqFEAwGRaERKRlEYbnK4RYubjxw9EAtkcIR3bUZSC40Pp8P4XAYlZWVotCIFBWisFylcK4vLuuLc31xDRpzfezNzU0olUrI5fKcH+9KIJnQuFwuuN1uGAwGhEIhAMnraEShEck3orBchaRqy8LNo88l4XAYw8PDsFqtYFkWOp0OJpMJJpMJRqMxJuYgkpro1GaZTBYzIiAYDMa4zqIbaoqxM5F8IArLVUa0lRKfRpxri8XtdqO/vx8ajQa9vb1gWRYulwtOpxPj4+MIBoPQ6/W80HCdhkWSE70JSHfomThdUyQfiMJylZBObUquhIUQgvn5eUxOTqKtrQ3Nzc18e5jq6mpUV1cDAPx+P5xOJ5xOJ1ZWVhCJRGA0Gnmh0el04iIYBbc5SIY4XVOkkIjCchUQPzI4lRWQC1cY5/pyu904fvw4TCZTymOoVCqoVCrU1taCEAKfzweHwwGn04n5+XkAiBEajUZzVS+C2wlLPKmEhmVZUWhEBEcUliuY6LqIdGpTKIoSVFhcLhcGBgag1WrR29ubUaCeoihoNBpoNBo0NDTwAX+n0wm73Y7p6WlIpVJeZEwmE1QqVdbnWoqLZybCEs92QhMMBhEIBMTpmiJZIwrLFUo2c1NomhbEFUYIwdzcHKamptDe3o7m5uZdL0YURUGv10Ov16OpqQksy2JjYwMOhwOrq6sYHx+HQqGIERqFQrHrz3K1EO8ajZ6uyTBMymQAsd2PSDJEYbkC2c3I4N1aLKFQCENDQ9jc3MSJEydgNBp39X6poGkaRqORf3+GYfhEgMXFRYyOjkKtVsNkMsFsNsNoNEImk237nrnOiBOa3VgsO7HdGGdxuqbITojCcgWx27YsuxUWp9OJgYEB6PX6jF1fu0UikaCsrAxlZWUAtmI7nNDMzMzA6/XGpDYbDAZIpaV9+edSWOLJZOhZvOtM5OqjtO8sER4hRgZnmxVGCMHs7Cymp6fR0dGBpqamtI6dy4VRJpOhoqICFRUVAIBgMMhnnCVLbS41a4WjUNaBOF1TZDtEYbkCEGpkcDZZYaFQCIODg/B6vTh58iQMBkNWx841CoVi29TmUCiE6elpbGxswGw2l0RqczGJoThdUyQaUVhKmPjalN1m7GRqsTgcDgwMDMBoNKK3t3fHGEYxEZ/a/Pzzz0Ov12NzcxMLCwsAij+1OZ+usEzZabrm9PQ0GhsbodFoRKG5AhGFpUSJr00RImiarsVCCMHMzAxmZmawZ88eNDY2Fu0Clw7cAlhRUQGz2Zw0tVkikSSkNhf6M3M93kqBeKFZW1tDXV2dOF3zCkUUlhIjlyOD0wneB4NBDA4Owu/3F7XrK1PiW9skS212Op1YX1/HxMRE0aQ2l+rCSwjhRYT7uzhd88pBFJYSIj5AL3Rq504Wi91ux+DgIEwmEw4fPlxSrq/dEJ3a3NLSsm1qM/cnH99NMcVYMoFzi0VbW+J0zSsLUVhKhGxrUzIhVYyFEILp6WnMzs6is7MTDQ0NV/UNvV1q8+zsLIaHh6HVamE2m3Oa2lzMMZbtiHbfpkKcrlnaiMJS5ORzZHAyV1i06+vUqVPQ6/U5OXYpE5/aHAqF4HQ64XA4ctq1uVSFhbvGMvkOxOmapYUoLEWMELUpmRBvsdjtdgwMDKCsrAxHjhwRbNd9pVdny+VyVFVVoaqqCkDyrs0Gg4HvCqDVaksmCC8E6Vgs25HOdE1RaAqLKCxFCsMwWF5ehkqlyltNBRdjIYRgamoKc3Nz6OrqQn19/VVxQ+YqZpGsazMnNAsLCyCExMRn0k1tvposlu3YTmiCwaA4XbMAiMJSZETXpiwsLKCqqipv7ifOYnnxxRcRDAZx+vRp6HS6vBz7aiG6a3N9fT0IIfB4PHA4HBmnNpeqsLAsm1OrNVpo4mfRxAuNOF0zN4jCUkSwLItIJBLj+spn5s/GxgZCoRCUSiWOHj1a8r20SgGKoqDT6aDT6VKmNsvlcl5kzGYzn9pc6sKSL9IdeiZO1xQOceUoAlLVpkgkkpyOCuZgWZZ3fdE0jYMHD+b8mCLJSZba7Ha74XQ6sbS0hLGxMT61ORgM5uX6EJpCF3amKzTi0LPsEYWlwGw3Mlio+SjbEQgEMDAwgHA4jMOHD2NgYCCnx+MQb9D0kEgkMJvNMJvNaGtri0lt9nq92NjYgMVi4S0ao9FY9JZmvi2Wndhu6JkoNNlR3FfgFU50bUp0uwuOXAuL1WrF4OAgKisrsXfv3rzugEu1uK/QRKc2+/1+GAwGqFQqOJ1OTE5OIhAIxKQ26/X6mHqQYqDQFstO7DRd02q1wuv1or6+XhSaFIjCUgDSrU3JlbCwLIvJyUksLCxg7969qKurA3C5jqVUffe7oVQ/r0wmSzu12WQyQafTFXxRLzaLZSfiEw38fj9sNhtqamoSpmuKQrOFKCx5JpPalGzno2yH3+/HwMAAIpEIenp6oNVqY47HnePVekOUEsmsvnRSm41GI98VoBBdm7nNVKnCsizvFgNSj3G+mqdrisKSR7gCrnQr6Gma5gVICCwWC4aGhlBVVYXu7u4EF0l0i3OR4menDUCq1GauK8DMzAxoms571+Zid4XtRLwwZjtd80oWGlFY8gC3m+GyvtI1kWma5luK7waWZTExMYHFxUXs27cPtbW1SZ8XbbGIFD+ZWpbRqc2NjY1gWRabm5twOBxJU5tNJhOUSqXg511qrrB4GIbZNm6VidBcqdM1RWHJMbtpyyJEjMXn82FgYAAsy6K3txcajSblc7nzKsUU1quR3bosaZqGwWCAwWBISG1eXl6OSW3mMs7kcvmuz/tKcYWlS7pCcyVN1xSFJYfsdmTwboVlfX0dQ0NDqKmpQVdX1443Qz5dYcW4Y73aLbXo1GYAiEQicLlccDgcmJ2dhdfrhVar3XVqc6nH8HYrjKmEhuvcDJT+GGdRWHKAUCODsxUWlmUxPj6O5eVl7Nu3DzU1NWm9TnSFlRa5XqClUinKy8tRXl4O4HLX5ujUZp1OF9O1OZ2dfKlbLAzDCFortJ3QcK7w97znPXjrW9+KO+64Q7Dj5hJRWAQmfmTwbm6gbITF5/Ohv78fANDb2wu1Wp32a7kLXHSFlQb53vkn69rMFWuOjo6mndp8JQTvc1kblExopqamBE3kyTWisAhEdFsWoeamZCosa2trGB4eRm1tLbq6urK6eZPNZMkFpewKKSYK+T1yqc01NTXbpjZzQqPVavmNSyn//tywvXxBURR8Pt+28dFiQxQWAcjV3JR0hYVhGIyPj2NlZQX79+9HdXV11sfMl7CI7J5i+p22S23mJmtyqc0sy5Z0Ie5OWWFCQwjh41ulgigsuySXI4PTERav14uBgQFQFJWx6yvbYwoBV0xWyi6RQlPMC/N2qc2rq6sIBAJ4/vnnc57anAsKESPyer2ixXI1kI+RwTst8qurqxgZGUFdXR06OzsFudjzYbE4nU709/cjFArBYDDwmUj5GmiWjGJdoLejmIUlnujUZq5QuKqqik9tvnjxIpRKJT8aQKjU5lyQ6xhLMrxeb0nNRhKFJQvyNTI4lbAwDIOLFy9ibW0NBw8eRGVlpWDHzGXwnhCC+fl5TE5Ooq2tDUajEW63Gw6HAwsLCwDALywmk2nX1tfVQKkISzSEkJSpzU6nE3Nzc/B4PIKkNueCfFvaoVAI4XBYdIVdyey2NiUTkgmLx+PBwMAAaJpGb28vVCqV4MfMhcUSiUQwPDwMp9OJ48ePQ6fT8TdLXV0dCCEJVeAKhYJffEwmE2QymeDnVcoUU4wlE5IF73OV2pwL8m2xeL1eABCF5UpEqNqUTIgXlpWVFYyMjKCxsREdHR052TXlwhW2ubmJ/v5+KJVKnDlzBnK5nP8eo4+r1+uh1+vR3NwMhmH44ry5uTkMDw9Dp9PxIlPIhaVYKCVXWDScxbId8anNgUCAF5qxsTGEw2F+PADnRs2XFZFvi8Xj8QCAGGO50mAYBoFAgM8vz1fzOE5YGIbB2NgY1tfXcejQIUFdX/EI7QrjxLC5uRnt7e1pf28SiQRlZWUoKysDsLWDdTgccDgc/MISHZ/hUlmvJkpVWLJJN1YqlaipqeFTm7nxAA6HA0tLS2BZNmlqc67OP98Wi0ajKalEF1FYtoGrTQmFQvjDH/6Aa6+9VnDX03ZQFAWGYXD27FlIpVKcOXMm55kzQrnCWJbFxYsXsbq6KogYyuVyVFdXo7q6OqZmgrNouFRWTmiy+Z5K0bVUqsKy25YoarUaarWad6MmS22OFhq1Wi3Id8Ul7eRzkfd6vYKdf74QhSUFyQL0+a5It1gsIISgsrIS7e3tebmYhXCF+f1+9Pf3gxCCnp4ewYPw8TUT8ams4+PjUCqVvMgYjcYrMj5TikIICF95nyq12el0wmKxYGpqCjKZTJDUZu47L4TFUkqIwpKEZLUp+RSWSCSC0dFRWK1WAEBHR0fediu7dYXZbDYMDAyguro6rcaXQhDfpTe6eeLMzAx8Pl9CfKaU3AqpuJpcYZkQfT1w8brors3Rqc3cn3RTm6M3mvnC4/EUZCDbbhCFJYrtalMkEkleevVwgW65XI6TJ0/iueeey6tPN1tXGCEE09PTmJ2djRl3nIp8Nk8MBoN8fIYb18stKGazueR2gxylKiz57hW2XWrz/Pw8RkZGoNFo+I3HdqnN3KYrnxaLz+crqYwwQBQWnp1qU4Se5pjs+NwMjObmZrS1tfEXcT6FJRuLJRQKYXBwED6fD6dPny66Qi6FQhET+PV6vXA4HLw/XiKRgGVZOJ1O6PV6KBSKQp9y2pSisBS6V1iy1GZOaHZKbWYYJu+THz0eT8nVdInCgvRqU7jFJxdEIhGMjIzAbrfjyJEj/AXPkc/YTqYxFrfbjb6+Puj1evT09BR9LIOiKGi1Wmi1Wt4f73a7MTw8DLvdjuXlZahUqpj4TLEU5sVTqjGWYmubL5fLUVlZySeYxKc2cx0iuNhMIdq5iBZLCRFdm7LTyOBcWSwbGxt8jUdvb29MULEQEx3TdYURQrC4uIjx8XG0t7ejubk54zG5xQCXTaZUKtHc3Ayj0cjHZ6ampuD3+6HX63k3iV6vL5pFUXSF5YZUqc1OpxOLi4tgGAYDAwN5SW0GRGEpKViWRSQSSbsti9AWS/TC3NLSgra2toTj5ztpgDvmTsfjkgvsdjuOHTvG+66zOVaxIZPJUFFRgYqKCgBIWi8RHZ8pdBpoMX6HO1FoV1gmxKc2c3VUJpOJd6VSFBWTCCD0NVFqLfOBq1BYouemcDu+dC4CIS2W6PYmR48e5YsAUx23mCwWj8eD/v5+yGQy9PT0lExH2mzhZo7U1tby9RIOhwN2ux3T09N8GivnOstn40TRYsk/LMtCKpWisbExIbXZarViamoKUqk0ZvOx23uE65tWSlxVwhLfliWTIJxQFgvn+lKpVOjt7d0xUFwIiyWVsHCDxBoaGnLWUqaYia6XaGpqikljXVxcxOjoKJ9dxMVncpl0UcoxllIURCAxkSZZavPGxgafgcjVVGWT2szh9Xpz2m0jF1w1whJdmxI99jNddmuxEEKwsLCAiYkJtLa2orW1NW1LqdCuMJZlMTExgaWlJRw4cIDv33S1E53G2tbWhnA4zLvNxsfHEQwG+bYzXHxGqAWVE5VSXKCLLXifCTtV3UskEl5AgNSpzdFCs1NyiJhuXIQINTdlNxZLOBzG8PAwXC5XxjGJQrvCAoEABgYGEIlE0NPTU3K+3kzYrQUgk8lisov8fj9fPxM9FoBzkahUqqyFoZSFpdRdYZlYodulNk9PT8Pn8/HNNFM1VxWD90WGkHNTsi2QdLvd6O/vh0aj4Tv7ZkIhXWF2ux0DAwMoLy/Hvn37rvpuwpmiUqlQV1eXMBbAarVicnISCoUiZv5MNvGZUhSWUnaF7bZPWCapzTRNo7q6WvCWLsvLy3jggQfwm9/8Bj6fD+3t7fjmN7+J48ePA9haN//pn/4JX/va1+ByuXDmzBl8+ctfRkdHR9rHuGKFhWseKdTclEwX+PihVi0tLVmdQyGEhWEYzMzMYHp6Gl1dXaivry/ZhaBYSDUWINpFotVqedfaTmMBSjW+AlxdFstObJfa/Mgjj+AXv/gFamtr8eSTT+LQoUM4dOjQro7vdDpx5swZvPKVr8RvfvMbVFRUYHJyknfdAcC//uu/4vOf/zy+/e1vo6WlBf/wD/+AG2+8EaOjo2knIlxxwsK5vrisL6HmpkgkEgSDwbSeGw6HMTQ0hI2NDRw/fjzmR8uUXE50TAYhhO9RdvLkSRgMhrwd+2oi1ViA6Hkj0fGZ+LHNpewKu5otlu2IT23+6le/irvvvhvvf//7MTExgVe+8pWQSCT467/+a/zrv/5rVsd4+OGH0dDQgG9+85v8Yy0tLfz/E0LwyCOP4BOf+ARuvfVWAMB3vvMdVFVV4Wc/+xne+ta3pnWcK0pYcjkyON3gvcvlQn9/P3Q6HXp7e3edfppPi2VjYwOrq6uQSqWCnLtI+sSPBYiOz0SPBeBcZ1yHg1JboAkhosWSJhKJBCdPnoTb7cZjjz2Gnp4eXLhwARsbG1m/5y9+8QvceOONeNOb3oSnn34adXV1+Nu//VvcfffdAIDZ2Vmsra3hhhtu4F9jMBhw6tQpnD179uoTllyPDN4peE8IwdzcHKamprKqRE9FvoRlaWkJY2NjMBgMUKlUoqgUkOida/xYgLW1NX5sM7DVTbqsrKzoW+lwcJZWqQoLwzB5/665GItUKsXJkyd39V4zMzP48pe/jA996EP42Mc+hhdffBH33Xcf5HI57rrrLqytrQFAQuZnVVUV/2/pUPLCkq+RwdtZLKFQCENDQ9jc3MSJEydgNBoFPW4uhYVhGIyOjsJiseDIkSNwuVzw+/05O16xUsw7/2RjAWw2G0ZHRzE3N4fR0VF+LAAXnynWhZu7lov5+96OfKdKc01ThWrsyrIsjh8/jk996lMAgCNHjmB4eBhf+cpXcNdddwlyDKDEhYWrTeEu1lz+4KksFqfTiYGBAej1epw5c0bw3YxQEx2T4fP50NfXB4lEwk+ndLvdJR0YvhrgKruBrThYOBzm4zMjIyOIRCIwGo18fKaYxjZfCRZLPrMjQ6EQGIYRLN24pqYGe/fujXmsu7sb//u//wsAqK6uBgCsr6+jpqaGf876+joOHz6c9nFKUlii27LkyvUVT7zFQgjB7Owspqen0dHRgaamppxZSrmwWCwWCwYHB1FXV4fOzk7+Ri/EpEyR7KEoKmEsgM/n4+Mz3JheTmSEaDGyG64EiyWfwuLxeABAMGE5c+YMxsfHYx6bmJhAU1MTgK1AfnV1NZ544gleSDY2NvDCCy/gfe97X9rHKTlhyWWAfjuiLRZu/ojX68155pTQCz3LspiamsL8/Dz2798fsysBhBlNLJJ7Uv1G0WObGxoawLJsQosRlUoVUz+Tz7EApS4s+Z537/F4+JibEHzwgx9Eb28vPvWpT+HNb34z/vznP+Oxxx7DY489BmDrd/nABz6ABx98EB0dHXy6cW1tLW677ba0j1NSwpJsZHC+4CwWh8OBgYEBGI1G9Pb25jyQJ6SwBINBDAwMIBgMoqenJ+kuKN/pzSLZkW66MU3TMBqNMBqNaG1tRSQS4eskpqen4ff7Y+IzuR4LIGQJQCHIt8XCtXMR6vs6ceIEfvrTn+KjH/0oPvnJT6KlpQWPPPII7rzzTv45H/nIR+D1evGe97wHLpcLL3vZy/Db3/42I0u3JISFq02ZnZ1FKBRK2mI+19A0jVAohPPnz2PPnj1obGzMyzkIJSxOpxP9/f0wm804evRoyl1qLmM6xU4pfe5sOxtLpdKYsQCBQICPzwwNDYFl2Zj4jNCz1ku5hgUojMUidBv+1772tXjta1+b8t8pisInP/lJfPKTn8z6GEUvLNGur1AoBK/Xm/cLMxgMYnx8nO+Xlc+iQSGaX3IdANIRRNEVVjoIcR8olUrU1tbGjAVwOp38WACpVBoTn9nt2OZSrmEB8p8VJnQ7l3xR1MISX5silUpzOnc+GXa7HYODg9DpdHzaZz6haRrhcDir10YiEQwNDcHtdqedBi26wkqDXMxiiR4LED22mRtyNjY2xo8FMJlMWY1tvhIslnwH74W2GvNBUQpLqtqUfAoLIQTT09OYnZ1FZ2cnysrK8Oyzz+Z9uFK2rrDNzU309fXxc1/SLXi8ml1hpUQ+fqPoav/osQBOpxOTk5MIBAL82Gaz2cxvvrajlFvmA/k//1KcHgkUobDE16ZED+PKtsNwpnBB7kAggFOnTkGv1/N9wkpBWFZWVjAyMoLm5ma0t7dnPIteFJbipxDTI1ONBeAGnRFCYrLNksUGrgRXWCEsllKjaIQlnZHB+RAWrlV8WVlZTJCbu5jyHbzLxDXFsizGxsawtraGw4cP8wHaTBDrWEqHQrtHko0FiB7RK5PJYuIzcrm8pF1h3BqVT2EpxVksQJEIS3xtSqqRwbkUFkIIpqamMDc3h+7ubtTV1cWcAycm+e4VlO5C7/f70d/fD0IIenp6ss57z5fF4na7sby8zHfwLXSvq1Jb7IrNqoweCxA9tpkbcjY6OgqtVgulUgmWZfMeqxCC6Nq5fFGK0yOBIhCWTGpTciUs3JTEUCiE06dPJ+3Lw51bvnfz6QiL1WrF4OAgqqur0dXVtasbNh+fcWlpCaOjo6ioqMDc3BxGRkZifPW5rqW4EiiEKywTosc2A1tFxU6nE8vLywgEAnjmmWf4TQUXnynmzwNcLu7Md7qx6ArLgGxGBkskEj6gLxRWqxVDQ0MoLy/HsWPHts1y2W3qbzZsJyzRCQZ79+5FXV3dro+XS4uFZVlcvHgRq6urOHLkCL+YcLNIHA4HX0vBLTjcCN98UGxWwHYUu7DEI5fLUVVVBZZlQQhBd3c3H5+JHtscHZ8pNgphsXi93pKciVQQYcm2LQuXFSbETRXd2qS7uxv19fU7vmY3c++zJZWwcG1lfD5fSisr2+PlYoENBoPo7+/na4GUSiVCoRAAJPS68ng8sNvtWF9fx8TEBJRKJcrKymA2m7NKcb1SKSVh4SCEQCKRJB0L4HQ6+d9coVDwGwuTyVRwVymAtDfAQuLz+QTZMOabgtyh3A+TKpaSCs7Fs9sAGuf6CofDKVubJKNYLBZumJjBYEBPT4+gN10uLBaXy4W+vj6YzWbs379/W4GOrqVobm5GJBKBy+WC3W7nU1w5F0pZWVlRde7NJ6VkXUWTLHgfPRYgemwz10RzeHiYbztjMpl2HNucy3PP93G9Xm9RWm87UbCtXzY74+jMrGx/YC4eUVlZib1792b0PoW2WAghWFxcxPj4uKDDxFIdTwgWFxdx8eLFrDtAS6VSlJeXo7y8HMDlFFe73Y75+Xm+cy/3Z7eV4aVCqbnCONJJN44f2xwMBuF0OuFwODA6OopIJBITn8nX5iLfGaGAmBWWF7gfNRKJZDzhkGVZTE5OYmFhAfv27UNtbW3Gx89XHU003EIfiUQwMjICh8OBY8eO8UFRoRHKYmFZFqOjo1hfX8fRo0f5RWK3RKe4RnfuXV5extjYGLRaLb/gGI3GKzYJoFSFJZt0Y4VCETO2mRsL4HQ6Y8Y2c797rsYCFMpiEYUlx1AUldXi7vf7MTAwwPv3s/2hClHjwbnfzp07B5lMht7e3pzuyoXICgsEAujr6wMhBL29vUmD70KNbY7u3MsNvIre2UYvOEI38ys0pfhZdlu5vt1YgNXVVYyPj0OpVMZsLoRyFRfCYhHTjTMk25si07YuFosFQ0NDqKqqQnd39652HIWwWBwOB8LhMOrr69HR0ZHzC3u3wXun04m+vj5UVFTs6GoUemGUyWSoqqpCVVUVP9KVc5tNT0/zBXtlZWVFExDOllKNsQhdeR+9uQDAx+QcDkfCWAAuPpPt8fNde8Ndw2K6cR5Id3FnWRbj4+NYWlpKOtAqG/JpsUSfP03T6OzszMtxOVdYpq6W6PhPZ2cnGhoaCrqjpigKWq0WWq0WjY2NCQHhkZERfsEpKysruYX6anKFZUJ8TC4QCPDxmZWVFd6K5SzZTBo8FqLPmcfjESzjM59ckcLi8/kwMDAAlmXR29srmOLny2IJBALo7+8HwzA4evQoXnrppZwfk4O7yTJZuBiGwejoKKxWa07jP7shWUDYbrfz8ZlwOIzFxUUwDJPX2plsKVVhyXevMKVSGZPKzlmxTqcTMzMzkEqlMe7S7dzMYowlfUrOFbbT4r6+vo6hoSHU1tais7NT0AshHxYL16uMcyVxvdPytZBwN326O3i/34++vj5QFIXe3t6CzlPPBIVCETOH5MUXX4RSqcTa2homJiagUqli6iiKsf1IKQoLy7IFq0OKt2KjxwJwyR9qtTomPhN9rvmOsXCJCqKw5IFU1fdcVffKygr279+P6urqnBw7l73KZmdnMT09ja6uLtTX14OiKP54+dotRVssO8GJIJe6XaoZWNxIhoqKClRXV/Pjex0OByYmJhAMBvmpikKltzKEYDQSgINlYKYl2CtVQpKh67EUKaYmlNFjAQAgHA7z7tL4sQAmkynvMRa/3w+WZUVhyQfJgvc+nw/9/f0AgN7e3pwVFOWqQDIcDmNwcBAejwcnT56MaeHALdb5Epbo46UieiplV1cXGhoacn5e+YBbrOPH93LprQ6HA3NzczF9sLiuvZnwfMiLr3ltsJPL11IZJcHdmnL0ytNz24quMOGRyWQxv7vf7+c3GEtLS4hEIlAoFFhcXMxLlqHX6wUAUVgyQShX2NraGoaHh1FXV4fOzs6cXrQSiSTraY6p2NjYQF9fH7RaLXp7exMyldJZ6IVkJ4uFYRgMDw/D4XDg+PHj/G7vSoKwBGyABeNjwPpZSKQSVOurUVdTB0IR3n2yuLjId+3lkgCis44YQjAc9mMoEgAAHJAq4WYZfMZnTTimnTD4F886/l5blZa4lKqwFJPFshMqlQoqlYp3l46MjCAUCsFms/FZhtHxmUw3GDvh8XggkUhKxr0cTclZLJywMAyD8fFxrKys4MCBA6iqqsr5sYW2WLhxr62trWhtbU16w3GPFYOw+Hw+9PX1QSKR8P2+ShXCEoTtYbA+FoyXgdlmBu2gYR+wg/WzQApPE62kQWtoVGgqUK2tBikj8Eg9cHgcGBkZAcMwMBqNWKgw4gdKFp6oN/pxGuf1n14bTsnUGbnFSolSnSDJ1dBxNVPRYwG4DQY3tpmLz+zWw+Dz+aBWq0vy+ypJYQkEAjh37hxoms6p6yvZsYVY4KOzqHaqSuf6qeVTWJIdz2azYWBgADU1Nejq6irJi50jtB7CxksbIMHLi74GW1YCi+2/ZzbAgg2wiNgvx/kklATNB5rRfaYbXq8XT7pt+E9FZEucMhQIG2EwGgnggGz7rLRStViK2RW2E9Hu6GRjAbj4zPj4OILBIN92xmQyQa/XZ/x7eTyeknSDASXoCvP5fLBarWhqasKePXvyepEKYbFE7/rTzaLKd8V/dFsXQgjm5uYwNTWVdhfoYoYNsdh4cQMkJGDwmwDeIS8UNQqoNBr8NGwHSPaLvoPd+RorVWEpJVdYPNtlhcnlcn5sMyEkJj4TPRaAc52pVKodv4dSLY4ESshiYRgGY2NjsNlsMBgM6Orqyvs57NZi4VKhM40HFUJYuP5kw8PDcLlcCUkFpUpgLiCsqHAQgPEwGJeHYoLy2WCmd3ahlKqwXCkWy3ZQFMWPBYge2+xwOGCxWDA5OQmFQhEzfyZZfMbn82VUwFlMFFRY0m146PF4MDAwAIlEgra2NrhcrtyfXBKytViiG2Bm0wUg38JC0zTfX00mk6Gnpycn/cnyNQaZg7AE/ml/zt6fDbKw+3aX3EERgsoNDxijfMdFrBQXnCvVYtmO6LHN8WMB5ufnMTIyEtM8lRsLUKrTI4ESsFhWVlYwMjKCxsZGdHR0YHV1Ne/9ujiysViCwSA/9jjbBpj5FhZCCAYGBlBfX5/zTLt8EloNbQXmk0BAACUg18shUUtAq2n+v/aAD+f9K3iGXYSUUKBA4ZS1Cqes1aBxeZHcfGkTEmMEOJz9ORKKwtkpCyo2J6A2q/kBZ/E711KuYynV60molP/4LhDcBFWn04mxsTE4HA589rOfRXNzM7/mCP2d/cu//As++tGP4v7778cjjzwCYKvjx4c//GH88Ic/RDAYxI033ogvfelLWSVGFa2wcAFui8WCw4cP87nlhWgEyZGpxeJ0OtHf3w+z2YyjR49mXXGcL2EhhGBmZgaRSARtbW3o6OjI+THzSSprRXdSh9H1UdTU1sBYYwRLCJaCGxjxLGPUa8FqyLP1RBoIXcryeqZ6BcaQAt3u2PY1e1wSmAIUnMrsF35pyIQ9G5UIm8KwOq1865Ho2hnRFZZ/clV5L5fLY8YCrK2t4TWveQ1+/etfY2hoCJWVlbj++utxww034M4779x1stKLL76Ir371qzh48GDM4x/84Afxq1/9Cv/93/8Ng8GAe++9F7fffjuee+65jI9RlK4wj8eD/v5+yGQynDlzJibAXUhhSffY0QWEe/bsQWNj464WgVyNC44mEolgaGgIbrcbCoVCsPkpxULEHUHYluimojU0FHUKwArMRjbw3LoDo14LNpjgju85oXclCAsNCm+bUuBL+wJAlj+5IUQBBJAvyXHgNQcACeB2u2G327GwsIDR0VEoFApIJBI4nc5ddezNN6XsCstHkTJFUaipqcFHPvIRRCIRdHV14d5778Xjjz+OH/zgB/jLv/zLXb2/x+PBnXfeia997Wt48MEH+cfdbje+/vWv4/vf/z6uu+46AMA3v/lNdHd349y5czh9+nRGxyk6i4Wr7WhqakJ7e3vCDZOqpUs+SMcVFr1Anzhxgm/nvRtybbF4PB709fVBqVSit7cX586dK1lXSypSWSuqVhUCbARPKl1YC1iAQPrvqQ8nL4g7ZpPhb0eA7x4IYZNk8LsRwByksMe1tXiRCAEJEki0kpjWI6FQCBMTE3C73RgeHgbLsglzZ4oV0WJJH6/XC71ej56eHvT09OAf//Efd/2e99xzD2655RbccMMNMcJy/vx5hMNh3HDDDfxjXV1daGxsxNmzZ0tXWCKRCEZHR2Gz2WJcX/EUsytsc3MTfX19UKlU6O3tFawSN5fCYrFYMDg4iIaGBn7eSyEGmgkJYQkYLwNmg0FkMwJmk0FwOYkFIgGUTUo85Z7DmjSU0TFkDI1T1tT96I7ZZHilpAIXleHEynvvpcr76I37JR1/65TictyG2irIjEcul/P9yvbu3QuPxwOHwwGr1cpnHHGxGZPJVLCmj8kQLZb08Xg8gnoOfvjDH+LChQt48cUXE/5tbW0Ncrk8YSNcVVWFtbW1jI9VcFcYsLUg9/f3Qy6X71jbUWhXWKoFl0syaG5uRnt7u6A3Ty4WekIIpqenMTs7m5Cplu9sLSEI28Pwz/oRcUXAeBjsUOcIYEtUaDmNpeBGxsd701w71Mz2t4+EpnBIrsYh+ZYFQQhBaCWEv51W4gftwZg4jDlI4a1TChyzXW7pIyuXgZImv464GAtFUdDpdNDpdGhqagLDMHz9BDfoimukWFZWBp1OV9CFvVSD94SQvFssPp8PTU1NgrzX4uIi7r//fjz++ON56ZhRUGHhhkNdvHgRzc3NaGtr2/GHk0qlYFm2IMFLLtYRfXOwLIuxsTGsra1ta2nt9rhCCkt008vTp08nDBIqNWEJrgaxcTZzcVC1blW3V8u1GMzgdTU+DRp8qYcvERBYlX64hyMoJyrI9DJI9BIEF4MIrYVwDDIcsUkxYWTglhMYQlvuLzouKCOvjLV4CUsQcUZAIgRSqxQmmwmOZQcUDQqo96hB0VttR6IHXfn9fr6B5uLiIiiK4t1mZWVlOR1znYxSdYVx90M+LRav1yuYW/P8+fOwWCw4evQo/xjDMHjmmWfwxS9+Eb/73e/47gHRVsv6+npWneILKizj4+NYXFzEkSNH+BthJ7gflmGYvJv43LE5YfH7/ejv7wchBD09PTntqiyUsHg8Hly4cAFqtRo9PT1J3XX5bCEjBL5xX8avkVXKINVvXT+nNynMvziM9cYauKp2dj3sc6UeZPZC+Rqeqlnm/97o0eL2yXYo2NgFiQaFLtf216+sfMt6IYQguBiEd9QL1seCgMCj8CBAB6AMKVE2WgaKpqDek3j9qVQq1NXVoa6uDizLYnNzE3a7HSsrK7h48aLg/a12olRdYZyXJN8xFqFaulx//fUYGhqKeexd73oXurq68MADD6ChoQEymQxPPPEE7rjjDgBb6/PCwgJ6enoyPl5BhaW+vh6NjY0Z7Zq4Cz8SieRdWPiutZfcDYODg6iurkZ3d3dOLzihFvq1tTUMDQ2hqakJHR0dKW/wfGShcQhheZJIZudKySjoDm9ZHLLvfAe1992HD7EsWJrCf338b/Dcbddt+/omT6K14pIF8dWu4YTHF7QenC+zoNea+Wjszb5NqNpU8E/7wWxsLWwryhUMG4cRkF7OMlBGlDi0egj79+zf9v1omobBYIDBYEBrayvC4TCcTifsdjsuXryIcDgcM3dG6KpvbmBdKVos3P1XqsKi0+mwf3/s9aHRaFBWVsY//u53vxsf+tCHYDabodfr8f73vx89PT0ZB+6BAguLTqfLOMOLpumYAVj5hLuoZmZmsLS0hL1796Kuri4vx92NsBBCMDk5ifn5eRw4cGBH07bULBZlkxLeIW/qJ1CARCeBVCeF1CSFslEJWkmDnVyA4r77QHGLBktw50OPYaTn0LaWy9f3jEIblqEyoEKlX42ghEFfWWIrfI411Tbntg3MBgNPn4f/+4pyBS+VJY6pDkgCeEHyArQOLZrNzWm/v0wmi+lvFT13ZmZmBjKZLKZ2Jn6kQ6Zwm5VStVgoisp7jCWflfef+9znQNM07rjjjpgCyWwonnSRDEg27CsfcLNY1tfXk8YmcsVuhCUcDmNgYAA+ny/tyv98WixCoGpTIbQeQtiSvJ2Kqk0F7UEt2DCLsDUM75gXIUsI8pf6YYj7XiUsi8rFtR1dYh5ZGB5ZGDO6nWM71f7ExYFSUFC3q+EdSU90CAiGjZcsovh1+dLfzy2cQ6OpETR1abw0S0DCBLRi58WQoihoNBpoNBo0NDTEtIXn2o7odDo+20yv12e8yBZi1y8U+c4II4TA6/XmdI156qmnYv6uVCrx6KOP4tFHH931e5eksBQiM8zlcqG/vx8UReHgwYN5ExUge2Hh0p81Gg16enrS3nGWWvCeoinoT+rhetoFZjPxuvBP+RG2hhHZiMTMWQlXNYNQNKioWhOGpmFpEG6stToixUlbbEsMZbMSmv0a0HI6RljCFIuQhIEqIk0I5NsV9hj3VzK8IS/+OPVH7K3aC71FD/+oHyRCIDFIoOnWQF4jT9taiG8LHwwGeWtmaGiIr53hhEal2r7NP3DZYilFYcl3RhggrCss3xRFunGm5FNYuMy18fFxtLe3Y3Z2tiDZaJkKy+rqKoaHh7NKfy7FOhZaTsPQa4Dzj86k3Ysj7kSXK2Ouge2dD6H82x8HxbIgNI1f/t29aQXw00HO0HjfxQOQEhqUlIKyWQllixJS3eXbTmKQgHEzeKnMgmerlhGSsJCwFMxBJcqCSv6/IOllvc055zDnnIMyosR+6X7URmrBuBlsnNuArEIG7UEtpIbMb3uFQoGamhrU1NSAEAKPxwO73Y719XVMTExAqVTyQpSqdoa7pkrRFVaINGlRWPJMvoQlEolgZGQEDocDx44dg9lsxuLiYt4XXZqm045FcZ2UFxcXcejQIVRWVmZ8vHxYLMFgEGNjY5BKpSgvLxckI0mikcDQY4DrWVdadSwAsHntm+E78HLI1ucRrmqCrFkJYGZX5wEANy8244DrskCRCIF/yo/IZgTqDjVkFTKwXhbMBoNpnRtP1C7yz2VoAqvKD6vqcrcAdTiIJg/SJiAJ4KWyl3Dcfhy1gVoAQNgahvMJJ5QtSmj2atJykSUjunamubkZkUgELpcLdrsdU1NTCAQC/JArs9nM186Ueowln64wlmXzHmMRkpIVlly3deH6lXFFm1zmmtDjidMhXQsiFAphYGAAgUAAp0+fznq3k2thcbvd6Ovrg16vB8uyfEZSdH1FtqnbsjIZdMd02HxxM+3XMOYaMOatrK09GwTGoAIuxeVKfRlF45C2BicMdbCGvJhftuNFWepq5PddPJCy3Ut4PQz3unvLapACYbB4vHZhx3P0SeUIUzSkhE2vBRkFgADDxmHUrNWAinpVYDaA4FIQ6m41VK0qUPTuFnpuc5CsdmZ+fh40TcNsNvPdAkpRWPJtsfh8PhBC8upyFxLRFZYELi2Xa9UffUEJNZ44E9IRlo2NDX6x7unp2VUqdi6zwjgXXVtbGz8EiaIo+Hw+2O122Gw2TE1NQalUxrQlyWS3qGxQgvEw8I0l1rcEaQYyloZUKYGsUgY2wMYE/WlQOGmrwu/rFqCKSHDUUYkbTnZBJ1OAsAR1Syq0DyrRE6nEF04MgsStNd0uU1JRYUHgkYUhZ2goWSnvmnuhcg1ueRrtZCgK62oj6ryOtL8HUEBAGoBdYUd5MLZOjIQJvINeBGYC0BzUQFEtXKFkfO3MxsYGHA4H1tbWQAjBn//854TZI8VOvi0Wr3cr9ia6wvJIrrLCWJbF+Pg4lpeXcfDgwaRzCAqROLBTlhbXTqa1tRWtra273hHmIiuMS3leWFjgXXThcJhP4+QykhobG2NcKxMTEwiFQjAajfwMi3TGuqq71ACFLXEhACNh8ZuWBYyo7Pxz3l59CAcVVbD/0h7z2v3OLRfWPqcZciJB8PebCHo2wEyHQAJb34scwHsXu/CHkyuYbLwc/1jQbIKAxFgIVoUfP2ucgUO5FXw/ZanCy9frsCkL4YWK9PswbcpVsDI6VAbSt8YAIEAHIHGsQrY+h3BVM2+dAVtTLzee34CyVQntIa3g1gRN0zAajTAajaioqEBfXx+amprgcDgwNjbG185wmwi1Wl2UFk2+s8J8Ph9kMlneOyMIRUkKSy4W90AggP7+fjAMg97e3pSumGJyhUULoZDtZIS2WCKRSEwLmZ12YdGulej6Crvdjunpacjlcl5kUlkzFEVB06WBqlUFm9+Lz1gSZ0r819oAlKwEr9O2oMWj58VARmgccVz+Ltn1CNi5RNer2i/B659uwC+uXeTFxSuL4FzFGsI0i01ZCEtqb4xbDQBeqFzHusoHGhQidKKAm6RKKGkZVkOJAmJX6mAKeCAFSbsrv/nCOTT+60dBERaEomF750PYvPbNMc8JzAQgr5RDUZu7hYxzJ1VVVaGqqopPqY3+bbnaGe633W3tjFDkOyvM4/EUrcimg+gKA2C32zEwMICKigrs3bt3251JsbjC4idTChnkE9Ji8fl8uHDhAhQKBU6fPp1xx+dk9RVck8XJyUkEAoEYa4a7Gd2RAAY219C3uYrFoDvl+wdoBv/dMoVKvwqnrNXocptiUn0JIWCXk8fzKFAgIHjli9WYqt/g3WLPVK/s+LnmdCmsDgL0GBrx0sZy8n+nKKxrTKjzOkCww8gXAqjCcnR/5v/jU6opwqL82x+H78DLYywXAIg4IzkVlviqe4qioNVqodVq0djYGDOyd3Z2FsPDwwkNNAuVqpzvGEspjyUGrnKLhZuYODMzg+7ubtTX1+/4mmKwWLjgt9Fo3NVkylQIFby32+3o7+9HbW2tYCOO45sscrEZu92OmZkZODU0ZgwsVqkMBqsAsKj8+L/GWTwTWsYrVuvRtbE1+wQeAmwzxp4CBb1PjnqLBovV2VXYx70hfm2fgHHdjj0Lq7Ak6V+2KVdhmZjRuulFSJJiINmln+/ImBKSuOuVYlnI1ucThEVqzO1ysNPiHD+yNxgMwm63w+FwYHl5GYSQmLkz6dTOCEUhYiylGl8BikBYslnEJBIJgsGdJ/xtR3SH35MnT8JgMKR97EJaLNwgtPb2djQ3N+fEVN6tKyy69iddwc4WtVoNtVqNhoYG/Mkxh+fsF5M+z7huR2WKhToatzyEnzfNILjUhEPO8rT7kGn8wt1KZ372JN7+0FdBsyRl/7JT9j1YVm1iXrsObdgPQ8gPaVShp5JRYr9rP6ooKqEIlFAUWEXsokxraMhrhJkflIpMG1AqFArU1taitrYWhBBsbm7ySQATExNQqVQxDTRz2TuwEFlhoissz+zWYnG73ejv74dWq0Vvb29GftxCWSwMw2B0dBSrq6sZdYPO9nhc+5pM4cYIrK+v4/jx4/zUw2QIedOwhOBJ12zSf0tnoY6nr8yKQ85yUFIK6UiLVyVM+rtx3c6fK5C8f5k+JIeMpTFmcgFQwCdTwKIyQB0JocutxRFHLcqCZVBUK6DYr4D3U5+F5uMf2ioCBUARgrp/fmNMrEW7T7vrtOOd2E0DSoqioNfrodfr+doZziU6MTGBYDDI186UlZXxqc1CwbJsXpveejwe0WLJN7vJCuN2/NlmUBUiK4xhGPh8PlAUldP2/BzZusJCoRD6+vrAMAx6enp2dFUImXnmigSSzqlPZ6FOhl1xqThRSwEypHSHEQBeVRhLlem5wa6VN6JhVIXKgBq/rZ/HqDE2fbhyYZU/V474/mUHneV4smYp9o0pCj6ZAmNGFq9cM6Ps1WWQaC+NeXjvu7Dsb0TdP98BiitSjIq1UO0NkNfl1loBhG2ZL5VKUVFRwSesRDfQjK6d4YRmt9NcGYYRbCJsOni9XjHGshuydYVlurhzO36r1YqjR49mPfKTpmmEQpmNsd0NLpcLo6OjoCgKp0+fzoufNxtX2ObmJs6fPw+j0YgDBw7kvTbBEU4+kyWdhToZEZogSDNQsBLQddKkWWHAVvBcU63ErYutsCn9+FPVatLnSUDhncwBVI4pwPq3vtvXLbbghLUSL1ZYcNHgAEsBlsYasDQVc87R/ctMQQWW1R4EJcmvf7+UwfyJAGq0l79/SkKBDvp4UeEfvxRrUb6hOy8ul1y2zOdcovX19WBZlm+gyW0ktVptjNssmwaaYowlfQouLNmQqbB4vV709/dDIpHsOPo4nWPnK8bCTdesq6vD2tpa3i7sTLPCuIJSoepossGeQlhOozchxsCm2WjyP/eM4KblJrTBALoZW9lh0ZaLDKDrpKAMEnRumLCuSj1w7LrlepQ7ZGDjes1UBzR47WIzTlmr8FT1MmargP/6+N/gzoceg4RlwdA0vvfx9/AiaAjJMbtDR+XRlQGcmZkH29YGUleHwGIgacNNQtOgDnZAVpaflN58DfmiaRomkwkmkwltbW0IhUK8NTM6OopIJBKTBJBOLCPf6caixVIAMmnpsr6+jqGhIdTV1QmSmZSPGAvLshgdHcX6+jqf9bW6mnwnnAvStSIJIZiensbs7GzKgtJ8YQ/7kz5uDirhvvGvYPjdN7bqOGgatrsehELfAGD7yZMeWRj/0zyFWxabsR9loA004CEgEbI1i14b257EI00dl3qhYh2VATXqfBrMaTcxqXfBJQ8mFYnnbrsOIz2HULm4BktDdYxllTJN+RLR8SRC03B98GFsHrgdG1VlmPjQB3Dd5x6BhGXBXvoelK9o3/b9hKRQ8+7lcjmqq6tRXV0dUztjs9n4uqjoBprJYq6FyAoThWUXZLODScdi4ZoxLiwsYP/+/aipyXyCX6pj59JiCQQC6OvrAyEEvb29UKlU2NzczGsmWjqusEgkgqGhIWxsbOR1Nk0qkrnCrvnJH9H1qa9uCQoouG66G+5X3wXGXINXrG7gR62Tab3372sX0OTRQReRAzrqUhsuAkm5BKzt8vfU6NVhyGxP+h4b8hC+3zoOkubl7qoqy7jLcnw8iWJZGD/3AJ7/zz34bUcYoa7T+MP1j/KCtV9zELfp8rcEFMP0yJ1qZ5LNneHuh3xbLEIVPBeCggtLNuwUvI8vHhTSV5nL4L3T6URfX19CoWa+29jv5Arz+/24cOECpFIpenp68hrUTIUtTliM63a87dNfuRysBoHh91+H+9V3AQCavXqYggo4FckD/tGpyWEJi6erl/HapRYAAKtnwfpZULZYldjnMmNd5cNL5Zak55iuqGRLsngSxbK4GL6AkGQfgFjB+hOWUbdhxgl97qegAsU57z6+diYQCPBus8XFrY7TJpMJwWAwr0k7Xq8XLS0teTue0JSksGy3uDscDgwMDMBsNuekeDAXi3x03UdnZycaGhpibkBuoRdiPnw6bOcKczgc6O/vR3V1Nbq6una1ixPyszjiXGGpFtnowsBqvzpBWFKlJo+YHDhir0CdXwt6gwaNxM8doQiUjDDXW+96DYZMNmzK00/7DqqVYCkKNEke+E/GL6xjOKCpglKS+6WgUK6wTFAqlTG1M1wDTbvdjosXL2JhYYHPNBNi1EMqSrllPlAEwpKtK4wQEnOhEkIwNzeHqamppIuzUAhtsURnq3EzX+KJ/oz5Cn4mE08umaCzsxONjY05P4908TNh+NjYBThZdhWhaYSrmvi/6+K6EO+Umvxf7eO4fa4NrZsGsNRWt2KvNAyPNAyPLIQnauNSgLOkwaPFMXslFKwEf4xPK04BL4iE8K1e4gP/yfCzEYx6LTiqrxXk3LcjX9evUFAUBYPBAIPBgJWVFXR1dYFhGDgcDoyPjyMYDMJoNPLxGSFrZ8Q6lgLA7RK4TI1wOIzh4WG43W6cOHECRqMxZ8cWMnjv9/vR19cHiqK2zVbjhCVfO754i4WbmbK6uppS/AoFSwj+z5ZYbe+qKkvIrvrhx/4GoW4Wp6wRKFkpDjjL8OeKdf416aQm/6R5Oqef51WrDThsqwANCsdslZjWubGg3T5gb1y34+0PfgVcP0sKAENRePhbD2J+387B+YWgG0eRe2EpBYslFVwdi06nQ2VlJQgh8Pv9fMuZubm5mHHOZrN5Vy7iXM+7zzUlLSyRSIQPdqtUKvT29ubc3y9U8J5rfFlZWYm9e/due8NFC0s+iBaWUCiE/v5+Pl4lZHGmELu7P28s4c8pGjYmz65aQ1+ZFSdslXDJYuuRdqohyRVmqQp/1/QyyGgJrCPWmBb3b0UH5rSb8EsiqFlxoXphDUtNVfjlwSAslyZMXveDXyO+SbKEECj8sW4+fUiOjSSzXxYCrlx9tBiKIXifLfF1LBRFxbQTiq6dWVhYwOjoKHQ6XczcmUw+O9fSpVQpuLBks7hQFAWJRILV1VVMT09nNdc9W3ZrsRBCMD8/j8nJSXR1daGhoSGtYwL5ExbOFba5uYkLFy5Ar9fnJF4lBCOe5IFyjmTZVUEJk1DIyAXsf/PO23DTN38GCSFpuZKEwBHx49cD53BYUo+qP/4S5d/6ON/i3vrOB9Fy7Vuge/rH/OO1FI3mdz6Ip29/JYYwiFf91y8T3pMBEFRd6lRMgGP2StyAZvyyaQ4j3tjvbDm4gQhhIaViFz42yCLijEBWLttKr94l+S4yFIp4t3sytqudGRkZAcMwCQ00U61XXEq0aLHkGZZl+c7EQs4hSYfdWCwMw2B4eBgOh2PHPlrRcONc82mxhMNhnDt3Lq+inQ1GWfbFrhwxtR/YciWxFIVnL1k8+eAljQ+HXljnxQPg2q58AuePVOCmb8c+XvntT+DUgWdw1KaICdZzSAD8/Ts/jp9+5B5UHL0b9UEtzDcY0ODTYwSxwsIQgmWvGzV+DSLOCCLOCMLOMFjf1vEM1xggL9+9J4Bl2aKZr5IJ3H2XiSjG1854PB44HA5YrVZMTk5CoVDEzJ2J37SVeuV9ydmlPp8P586dAyEE+/bty3uuN7ebz7QNDXfefr8fvb29aYtK/HFzDSEEa2trCAQCOHDgADo6OnIqKrvtF/ZKUyvUdPaLVULtx6XHaULwip/8AZ++5X0487Mnd3WO6eCnGPxe8kxMZTwA0CyLjdlnQcX99lyGm6y8HYRKfhvTLMHt//olNC1tAgzg+IMDhgvJf8ux8ytwP+uGd9iL4HKQFxVga06LEJRa8J6D81DspoGmTqdDU1MTjhw5gmuuuQZ79uwBTdOYnp7Gs88+i/Pnz2N2dhazs7OIRCLw+XyCCcunP/1pnDhxgo8P3XbbbRgfH495TiAQwD333MM38Lzjjjuwvr6e4h13pqSExWq14uzZszAYDFCr1QVxzXC7lkwWeZvNhrNnz8JkMuHkyZNZjRvNh7AwDIOBgQFYrVYoFApUV+c2tiBEE0qzTIUPNPZk/fpkAftoaALc+dBjMK4nL3wUEi7GEw1L05g+0gU2blHjMtwYcw1s73wIJMWixwkQAIABKgLJG4OuqlI30RRKWEo1eM/dd0KdOzdTaM+ePTh9+jROnz6N6upqeDwe/M3f/A2am5thMBjwq1/9CktLu880fPrpp3HPPffg3LlzePzxxxEOh/HqV78aXu/l3/yDH/wg/u///g///d//jaeffhorKyu4/fbbsz5mwX/ldHYw3Lz0/v5+dHd3Y9++fTmbe78T3MWVzrEJIZidnUVfXx86Ozt3DNLvdNxcV/y/8MILCAQC2L9/f0ntLM0yNT7d9ioYpJkLdrLFPB4uMyzXcJlsbPR3T1jUTy7g2fs/zIvHVluah/h6nM1r34yFzzyDtb/9fMLAYk6AbAo/nqhZxH8bnsGeF4cThHJVnVpYws7sRijEU6rBey77NFf3hEqlQl1dHQ4cOIBf/vKXeOyxx+D3+/E///M/aG5uxv79+/Gzn/0s6/f/7W9/i3e+853Yt28fDh06hG9961tYWFjA+fPnAWyNEfn617+Oz372s7juuutw7NgxfPOb38Tzzz+Pc+fOZXXMoo+xhEIhDAwMwO/3x7QOyaRfmJCka7FEIhEMDw/D5XJlNEgsFbkUFq7in8tQ29jYELSl/XZw1di7XXBktAT/0PJK2EI+/Mv8M2m/Lj4tOdm4X4aiLgfCc8xWTOfyd89ZTBcfeQILn7kVsvV53lKJOUdzDXwnb4HN70X5tz++NXuFptH/vo/i/467saxZwZmfPYn/X4q5NF5p6nuJ9bJggyxoxe5+o2KsvE+HfCYdyOVyHDlyBC6XC88++yx8Ph+efPJJtLa2CnYMt3trVDdXNnD+/HmEw2HccMMN/HO6urrQ2NiIs2fP4vTp0xkfo6iFxeVyob+/HwaDAb29vTGur0LMRQEuB9K3OzY3510mk6Gnpycr11c8uRIWrq34nj170NjYyC/yuXa7cZk2hBBEIhH+uLsVmXK5Gv/WfiPOuhfxU+toWkO6uLTk637wa7zqu/93qRcY+P9KCMHfv/PjaQ0IiybdqZXRVC6sJqYOsyxGw+dhaH4lWsynINlm0r3zFW/EcM9BWDyjGNijg73GDMC7Y/FnWSAqCYICpHoppKatPzKTDJRcmKywUrZY8oXP54NcLodMJoPJZMIdd9wh2HuzLIsPfOADOHPmDPbv3w9gqzu5XC5PqP+rqqrC2lp2lnrBhSXZDoYQgoWFBUxMTKCjowNNTU0JzyuUKwzYPjPMarVicHBQ0DnvwO7HBcfDsizGx8exsrKSMJFSqJn3yeBa03CiIpPJwLIsWJaN+T1pmub/ZAK1vAzJ9DTOtLWhq/ka/NI2gUFPejfHq773S943HC0uQPoDwjiymVoJpK6lealbA1fVNFQRCQ46ytFrrYGc3dpFExDMaTcxanRgQu9ESMICaIp53+2KPyNl1biZtEF7SLslJgYpKInwlkWpBu/znSbt8Xig0Why8l3dc889GB4exp/+9CfB3zuaots+RCIRDAwMYGZmBsePH085171QFkuqY3Mt5Lk4UHd3t6C7HCGtiHA4jPPnz8Nms+H06dMJY45zldrMCQrDMCCEQCKRQCqVQi6X83+kUinfGy0SiSAUCiEcDoNhmB3PSfad70Czbx/Ur30tNHv3ovI/vorblltx9/i+Hc+tdWA8sbdY3HNSxVqM6/aYuEUq6yCdBADONcdwsby4Whq/lMELlev4RcMsCAgCdAT/1TqO36vPITT1NNQ2a9L3TZUYsKf9KP6u/WXoOlILVZsKMrMsJ6ICiBZLunDCIjT33nsvfvnLX+KPf/wj6uvr+cerq6sRCoXgcrlinr++vp51Ak/BLRbg8g7Z4/Ggr68PCoUCvb2927qQChVjARKLJLkW8m63G6dOnYJer8/JMYVY7D0eDy5cuACNRoOenp6kmXWZDvpKB0IILyjcMeKPCcTGsKL/7GTNUMvLUNx3H5+WSxEC5Sf/CbJxL8yv+RtU+dUpB3Gd+emTePtDX9nxM7A0hWpSDrw4zLu3klkm1rrKrKZWcmw3j4VjWu/GlM6Nee0mWh7/vx2to4RYEk1j7qN/D2lQhpeePweTycR3+d3NILztKNXgfb4tFi7VWCiLhRCC97///fjpT3+Kp556KqFr8rFjxyCTyfDEE0/wbrfx8XEsLCygpye7jMuiEBYAWF1dxfDwMBobG9HR0bHjBSiRSPI6Ijj+2Nwiz4mhUqnMaUsZIYTFarViYGCA/45TXbic0AvluuBEhduxpvOe0cLBuc04YeKsGe5cLREfxl74Ld4QX+sBoPzH/wbrmdegdVOfVFje8aIevQ89BipOR+OD+ATA6jW9+Iu/vR/UpUFZP7v3L3DbF7+XYJk8/K0HE9xZLE1DQevQ9eIw1tKIuaQzj+WpmiWUrzrx19vETjgq/SqoT9yFhe+8A1UmJ9jWVlTU1aH80obObrdjfX0dExMTUKlUvMhkM8Y3FaUavC/E9Egh27ncc889+P73v4+f//zn0Ol0fNzEYDBApVLBYDDg3e9+Nz70oQ/xM2je//73o6enJ6vAPVAEwkIIwejoKJaXl3Ho0CFUVlam9bpCusI4i8VisWBwcBANDQ1pieFuj5mtsER3ft63bx9qa7dvOMjd/EIISzaiEk8ya4ZhGFiCHjzunMagdx2GSiVupZAQ+KYJwR8lz8LMXp/0vd3usYSiRAA495qX4+TvnovJFKt96rmomAuLN3zhewlz5CUsi/b+i/jJ++/EG77wfb4B5uirrsc97/kI36rlxx99H56849qMv4toHIogyq0zKa0jYqrBXpcZ3S4zKoJb9SsEwJq1CuajZkhwuXhPp9OhubkZkUiEbxM/NjYWM8a3rKwMKlXyOph0EC2W9BC6s/GXv/xlAMArXvGKmMe/+c1v4p3vfCcA4HOf+xxomsYdd9yBYDCIG2+8EV/60peyPmbBhYWiKH63n4lKF1pYVldXYbPZBJ1OudMxsxEWhmEwMjICu92edtpzdJv+3RDtyspEVAghYEAgAZX0Nc5IAL+3T+JF9xKf9eWqKsNP7ns77viP/4qxNBiaxmJLBYYqExtVGtftcAVXQECBissfO/XbP+E3f/UGvOYbP4npGhwNRQgIRcWICwHw5s9+ByxN4SfvvxPze9sQVCnw9+/8+OWhY4TFm/7ly7jwsv277kOWLNjP0jRezpxB2XgrqBQZZI7fOaA/rQeJELABdutPcOu/0oAUFaEKVCmrIGmWwCV1wWKxYHJykrdmzGZzxvNIRIslPYQeS5zOfaxUKvHoo4/i0UcfFeSYBRcWAGhra8t40SxUVlg4HIbX64XX683rSN5s4h5c52cA6OnpSdt3Hm2xZAMXpI+uWE5nQWEJwW9s4zjrWoSXCYECIKMkkNMSyCgJZDQNGSXBanATbJJE4sff8XoAwO1f2HJPbddEMqY/GEWBJbGZLDQheM03f5pgAUXD0DQef+9duPEr396qG0FsFtkbvvB9fOyXjybNyKJTxFwyTVFONh7gl393Lw6o23Z87ca5jW3/nQkyYAYY1L6iFk1NTYhEInA6nfzQq3A4zMdmzGbzjhvDUg3e5/u8S33IF1AkwpINhbBYuMA3RVFobm7Oa/fRTC0Wt9uNCxcuoKysDPv27ctoZ7mbbso7Bem343nXPP5gvzzvhAAIEQahDH7nx9/xerx445ltA98J/cEuWR4g8Ys/ARvnXosfovXcbdfh7E09eN8PZ1D9/X+JeT3nkkq3HX+2KcrxwX65vh4HpnZ8WdoEl4KQmWWQSqWoqKhARUUF34E3urGiUqmMic3EX3Ol6gpjGKakXWGFoCiEJdspkvnMCltbW8PQ0BCamppieuzki0yEZWVlBSMjI2hvb0+Zrr0d2Vos0ZZKNvGUnVrg74SMorFHU471oBIT2+z2k1oQhCSM9SU0jcW3fgANP/jc5ccp4Hdvfx2efNvNvGit1Rnx0qt6cMsP6Jh4DSceyayKn77/L1C5sNW631VVtmMB405EB/vLA8Jm9FGyxN+RoihotVpotVo0Njby1gw3XTEUCsFoNPJCo1arS9YVlm+LRWhXWCEoCmHJhnxZLFyfsvn5eRw4cADV1dUYGhrKWwt7jnSEhRCCiYkJLC4u7mqcAHfzZ/IZhQjSq+jsLscyixOvsEdw4PDL8JXQMmyKQNLncW6moFqZ1ILoe+VJHHviHF8cOXrD9dAcfwOoHzwCrs0KRYAbvvcrPPm2m2Pee6BVjlPvfCimncrFv7wHrQNbXWRHeg5tucUW19A0Os2764RKUY6GJsIt3pSMgqp154B9vDXj8/lgt9ths9kwNTUFpVIJhmHgdruhUqlKai4LwzCCdM9IF5/Pl/MGsLmmZIUlHzGWcDiMgYEB+Hw+9PT08OapkOOJ02UnYeEKS7nYz25Maa5tTboWS7ZB+nja1WXo21zd+YmXoEHhzt9dwJlPPMwv5m0ffw/WkriP4t1M526+Bqd+/WyMBXH7F77Hx0goAF1/eAIvHL4GzXFZY8kW/EXNJtavux2+Ay+HbH0eitkh7P32w9h3SZBYisJ/feJvMNJzCB983ydjLJO/fOgxvPTQtwWbXmlR+WFR+lAZuBzziJ5KGd9nLCUUYOg1ZNwjjKIoaDQaaDQaNDY28nPih4aGMDc3h8nJSd6a4WIzxWzJiBZL5hSFsGTrCsvl4r65uYm+vj6+kDB6QJFQ44kzgaZphMPJu8x6vV5cuHABKpUq4VyzJR1h4WpKoudV7GaB6NZWAklGQGgkMrSqzNiIBDEfNUb3bqoGRy6JCrDVIj6Z+yiZm+nUr5/Fw996EAp/EJaG6pQtT4aNDvRQl11chKbhPHQINZQWEasfs+UqEAogFHChzIKXMVup3DUPvz0m04wmBG9/6DF87aH7Eiv8WRY1bgrWux5Cxbc/AfpSncwfPnQ/KpStKNsAGIpgVe1FUJLeNf/NjjG8ca4dbZuGmOmThKJhe+dD2Lz2zdu/AQXoT+shK9v9tSSRSFBWtvV7HDt2jBcau92O6elpyOVy3mVmMpmKzprJd4yl1Id8AUUiLNnACUsu+g9xxZotLS1oa2tLeH+JRIJgMJji1bkhlcVis9kwMDCAuro6dHZ2CvZd7NTWJT5Iz1k5u8EkU6FarsVayBPzOHPpN14OxmYx/enFP+Aou7M1kUo0FP4gJo5fbveSzGKYPrwH3/3Ee/D2hx6D/cwZTL3//QhVVqIOQB2A4zYfnjeGsSQL4rmqVYwaHej50xSakmSt0SwLUFRCLIelKPzvMRY3BO5A4MA1fAfjDnMNOha2nkNA4JGG8UTLEsaVzrS+z6erlrFn3pdkKuXH4Tvw8qSWC2fZyK/fB0WNcEP0uOtEIpFAqVRCo9GgoaEBDMPwsZnJyUkEAoGE2EyhrZl8Wyxi8L6ARBfLCbWbYFkWk5OTWFxc3LZYM1/THLc7JiEE8/PzmJycxN69e1FXVyf48VJZLLsN0m9Ht7YSa45YYQmwEQxuJvboWm6oAkvTWwv2Jdgk7qN0srKSBdi5VOXnbrsOtpe9DEfldQlVIVK5Ei/3q/AsNrAkC8KpCGKkOojXIbERH0vTsNdWAElEZ0XjxQ9qJvBWdEBtroFXEsaq2oUVtRereh9WFV4EqcwsdKvKD8ZqTygA5YZ/xQtLjGXzbzSCn/88wu94R0bHTAV37SbbpJWXl/P96rjYjMPhwMzMDGQyWYw1U4jhfvm2WIScHlkoikJYsnWFAcL96Nzcl0AgsGOMohCpztHCwrIsRkdHYbVaceLEiYR210KQyhUmRJB+O7o1FfijYyat57qqyvDEJz+EG/7pc6AYBkQiwbP//HcJwe7tRCOaVD26KAJ0KWsAFkD856VpgBAcC2qxLA2i9+eXYjmIbQvDUhT+6+PvgcIXSNodoHJxDRPHy/CFvYOQsTTC9GUxMK7b0ZRh+30O+6ub0fIwHTPamBv+FY3EsRpr2bAsFPffj8j114MIsGlJN/1crVZDrVbz1ozL5YLdbsfU1BQCgQAMBgMvNLnqAByPGGPJnKIQlmzgfuhIJLLr/lwbGxvo6+uDXq9P2Zgx/tiFsliCwSD6+vrAsmxGRY/ZHi8aoYL029GiNkNBSxBk0xPu4Tfegpe9+T2gpqdB2trgUQYA68WE58WLBgDsiWomyZGsR1cFI4OGSBJL7y9BURQ0RIIWmz+2PgYASwHf//u/xtDLj/FpxTtZT9Gikm1tC8fTVQwOfP7zUNx/Py++jvd8KsFaka3PJVo2DAN6ZgaMAMKSymLZDi42w8VnfD4fH5uZnZ2FTCbjW82YzeacWTP5tFi4+qB81sjlgpIVFoqiBMkM42o+Wltb0dramtaFXwiLhaIohEIhnD17FiaTCfv378/pxR5tsQgdpN8OKUVjj7ocQ54kUfwUkLo6flfdGdhIKizAZdHIdLFWkfR2q5V2T5L6GGC9uQ6uqjIYJAr8zdHXY+BTqzj4sQe3tZ6A1O3345MTdCEZ9myYcL48sQ5o0m9H+B3vQOT660HPzIBtbYXEWA36KRfY4GUhCVc1g1CxdThEIgEr0ORCroZlN9cNZ83U19fz1gznMhsZGcmZNSNaLJlTFMKS7QWwmwWeG3S1vLyccc1HIdKNNzc34Xa7sWfPHrS0tOTcBRDd4VjoIP1O3FDWjmHPelrTH0ncs2qVerzM2IQ/ueaTPn+nxfo6byMqGvT4kWOYf42fSs86tZRpU1ojKlqKd9ceg0mqgvHu+/GrMycwOfpSyu4AwPbDubjX1Pg0uHG5EYaQPKmwAABDWEjq6njLQwJAd1IH97Puy88x18AWXYcjkSD4H/8hiBsMEL7qPtqa6ejogN/v52Mzc3NzMf9uMpl2lSkpxlgypyiEJVuyFZZgMIiBgQGEQiH09PRkvDvIZ7oxIQRTU1NYWlqCWq0WdPb1dnDiycVThJhLny6NKiPeXX8cT9inwRKCeqUBx/S1+NHaENbjMsaS5Re8XtGJljU17GE/qv1qfLf9Iu9e2mmxNtll6NKbUavQYSW4CQCwSsLwUgzUhE7a1JGAwEexmC1XpYzlPNDUi2rFlnuDZVlce+DlcNdUYMK9kPJ7SNVg8iRO48xkPbRhGdTMzgvmRa8N+7SxiSih9cSRE5vXvhm+Ay+HucEBan+HYKIC5L4BpUqlQn19Perr68GyLB+bmZ2dxcjICPR6PS80mc46yafFEolEEAgERGEpJNkIi9vtRl9fH4xGI44ePZqVXzZfFkskEsHg4CA8Hg/27NmD5eXEDr25gqKonAbpd2Kftgr7tFUAtsQ1MB8A42GAhHBarEiwERaOpxyojahRi60CwWvX6vCH2kUA22eIHbFXoM1twB9WJrFStcn/u8Fix6pjHW11+0EoEiMunMV0XuEBoVInAHx96kW8u+4Yqo0GftbM7dX7IYlQeNqb3Lrikg7e/tBjoC8VgNrvegh1qj3YiGxAI9OAMJc/xz6nGSMmR8L7BNnY+ic2yMI/7U96TPmJVuDQwbSsxUzIZ58wmqZhNpthNpsBAH6/n4/NzM/PQyKRxMRmtrNmuAzIfFksHs/WxkmMsRSQTPuFLS0tYWxsLOseWtHHzbXF4vP5cOHCBSgUCpw+fRputzuvVpJUKsXExAQqKytRUVEBg8FQkAaCwbUgNgc3EXFHgI7E5S7kCMNP+aGsV8I37cNm32bCc47ZK1Hj02BBu4kwVYM/3/sATn7xX/lCxNH3fgxvc7wC+siWas1pL9fLRMdjLNdeg4mPfByM5nJ7Dx/F4rzCgyXZ5bqmZAkANokfj82/iPdOHoW5Ww8sLCHwuyGcUdejrWYvFjWbeLxuMeHcOaF637MaUBWtYMw1oFU05mXzeMUrXoHATADewa3edWcsNUmFRb8kBfZe/rtvwgeJNbESX1Yug+ZAbnz7hewTplKpUFdXh7q6OrAsC7fbDbvdjrm5OYyOjkKv1/NCo9PpYs6Tcwfn69r3+baG0YkWiwDkOsbCsiwuXryI1dVVHDlyJGHGe6bk2mKx2+3o7+9HTU0Nurq6+B1uroUlOqZy8OBBfpc3MDAAQghfb1BeXi5Idf92hJ1hbA5uImS57LKhkvTACjEM3OfccMOd8G/R1Po1qPVrQKtoKN96Pzxv/0vQszNwucqhUVQCUfsTfViOZXgT4jGVTz+Dsueex2d++XUwRiP8FAurJIx0W3M5FUE871zAq/7+GZR/6+PQX6qEV7/zIVRc+2Z0bpjww5YJ2JSxvc5cVWWYP7IHDb5LrrQwC8i2rmtFs2LLmnMzMIWU6F2vwfNVl9vidLlMqIhczhxkAyyk3/g2qr+RWImvPagFRV/Z8+5pmobJZILJZEJ7ezsCgQB/nS8sLICmad6SKSsr49emfFksXq8XSqWy6LoPZEpRCAuQXguReNLJCgsGg+jv70ckEkFPT48gIz9zabEsLCxgfHwcXV1daGho4B/PxRz6aOJnqMhkMlRXV6O6uhqEELjdbthsNszNzWFkZARGoxHl5eWoqKgQtDqa8TLYHN5EYCGxkaQ2IoMVsS4cq9IPFgR0qlzgKEzXmCCvlF86Vz2wrxnGCAvPoAe+6ctji09bqzGhdyWPx0QikM7PYKZiH7LhIjOCt6WohNeYa3DjchO+1zae8Dqr0s8LC0JAma4MDMOAoiio9qngeX7LhfJySy3qfVqsqbwwB5XYs2FECCEEV4OQV8kR+NMUKr+RePzIDddDahSu0j6eYm2Zr1QqUVtbi9ra2hhrZmFhAWNjY7zl4PF4YDAYcm51eTyevNXn5JKiEZZs2Mlicblc6Ovrg9lsFjQ9l7NYhGwnw7IsxsbGsL6+juPHj8NkMiUcM1diFl30mCxIT1EUjEYjjEYj2tvb4ff7YbPZYLPZMD09DYVCgYqKCpSXl8NkMmW1gBCWwDPqgXfcu1WImIQanwazustuKq5b8WbECIO25fJ7gcArjcCi9MGhCELO0ujYMALPOGG+zgx52eVADS2loT+qh6JWAfdLbrB+FmUBFU5ZqzHaaBOsMSSHYWVp20r4ikjyTsIWZayg1tprATsgqZKAKqMgq5MhvLwVS2nx6NHi0cc8f+PsBmgVDXnfxaTH18rWAHRdfmx5GfT0NNi2NkGC+KXQMj/amgG2NqWrq6vY3NzE4OAgKIqKic3stn4uGZywlDolLyypYiyLi4u4ePEiOjo60NTUJOhFzQmUUMISCoXQ19fHW1XJ5ornSliyqaRXqVRoaGjgq6O59ujDw8NgGAZlZWW80KRz8xFC4P6zG4HF5O3uOWr8l63NmFoUisbg334U51/3aqwrfbCq/PBJY6+LZ8MreNNcO/AkoGpWQdWqgsws4z+volqBhd4Q/rA6hXXae6kYcueK/e0mPt641Ijf1cdmfSVLHuAq4WkFjYYbK2FeGocjHGeZqRKD7ZvnN1HxmgrQShqSwxI41h0gkdRWLetnEa5MUq9CS0Ad6OAD9rLvfAeK++7jO0YL0dqlWC2W7VAoFCgvL8f8/Dxe9rKXYWNjg7dm4mMzer1ekLWAmx5Z7CK8E0UjLNm4wpJZLFy7E4vFgqNHj/JVu0IS3U5mtzfL5uYmLly4AL1ej2PHjqXMUsuFsAjRnkUikaCyshKVlZUghGBzcxNWqxWLi4v8zceJTKo0T8+wZ0dRAYDWugoA04m1KITFgS9/Gj+4pR4uXfLf2yML43+bpnHXVBcwB/jn/JBoJFA2KqFsUGJG6sb3LQMgcUZtqiwvYPuqeHVYikPOclw0ODGvu5xQkCzTy3bXQ2DMNai6uQKUlEKtQp8oLAo/COIy0sIEkc0I5Ao5ZGoZ9Ef1cP95+1hTQr0KTcP38Gd5q4RaXuZFBRCutUspWCzJ4DLCaJrmrfa2tjYEg0E+NrO0tMRbM5zQZGvNXAnFkUARCUs2SCSSmFby3Ix3QkjKnb8QcGLCMMyugtjr6+sYHBxM2UU5/pgsywpiJXFBeqHbs1AUBb1eD71ez998NpsNVqsVMzMzkMvlfFyGa4/um/XBezG9iZz0xTD03fK0CgeTvt65imHrOHpwBoy5BoyXgXfMC++YF/+3ZxwGV6z1EW2NRHdBBnYutPTJIvjM/gtgk3ytz912HRyHevHmFyUIVzVdbq9ySdRqFToMx3Ue0NisiEw+A2XZnsvPpwBadXljo2pSQaKRwDvuRXA1mKzXJYDL9Sqy9XlQhzqgvnkP3zCTnp6O6SsGCNPapViC95mSavOoUChQU1ODmpoasCyLzc1NXmTGxsag0+n4uplMrJkrobMxUOLCEh28dzqd6OvrQ0VFBfbu3ZvTrApuIc7WgiCEYHp6GrOzszh48CCqqqrSOib32t2IQHyQPpc1KgqFgk/z5Nqj22w2jI2NIRQKoUZVA9OKKWnRYSqqveq0Z8hHE21dxM8kWVf60PybXyQMAjv962dStn1JR9ySiQrHSp0eAc+BmMdC1hAUlQrUKGLjI/Hnbn3ng/Bc+xZo92oh1cTewvJyOeTlcjB+Bv5ZP3yzPrC+xOuUMdeAMddA26Pl3ck0TYO0tIDQcU0rBWjtUoquMCC97uk0TcNgMMBgMKC1tRWhUAh2u50XGgAxsZntplFyrrBSp6SFhYuxcJlUnZ2daGhoyIvJnW3KcSQSwfDwMFwuF06fPp12IRR3U+5m5xffniWfN3p0e/TOzk5srG7Ad86XNIV4O+p8WkzUpNetmCPeuojOxAIAt/13ePunvsJ3HaZZgp5fPs3LXbIeXdmIWzTacKKl6/yTE9W3V6NWqYsZo5xw7t/6OBxHVNDqbwaQfHcrUUmg3auFpluD0FoIvhlfghWj7lZDVabirwlCCMJVVfB+9rPQfPjDfNNKIVq7lKorLBt3t1wu560ZQggfm+GsGa1WG2PNRL+/6AoTmGwuOoqisLGxgY2NjaSZVLkkm5Rjv9+PCxcuQCqVore3NyM/bLSwZEOu291nAhtkEewLgmKSn8Oafg0mygSFO3Fnd9BRhucrV1PGPlQRCSoDalQEVHjpUu+sZNYFxbIw/P7bMPzu62giid9p/JlJWBY3uyX4ebUMXhJOuxV/KnThJL89A4RsIci/8hg+/fA/XkpMiB0KBmy12e+674Ng6Q9h4TOfRsV77k35e1IUBUWNAooaBRg/g8BSAIyPgbxSDkW1AhRFxcw2YhgGkbvuguu66/imlWxtLehIhK+nyoZStVh22yeMoqgEa4aLzQwNDYEQArPZzMcf8zE98tFHH8W//du/YW1tDYcOHcIXvvAFnDx5UtBjFI2wZIrf78fMzAwYhsHLX/7ynLWPT0WmFgvnqquqqkJ3d3fGN9luhIUQEuPuKKSokAiB7Vk7zrrWYWEDqKSVOCErh+TSOam71GiobYBjyZFUWJSsFHdNdeHZqhVsarWg29txMKBC5ZwalX4VtBEZKFAgIBgy2RGUMMkzsSgKht/9J6hUw8wQKy6ElqDFtg/3BKvAUAQSQuHJl1XiY79MHtjfCX0yYQHg/t8RNP7LP/HnRROScC4cNEvQ8P99DF8+3IBXHLgGrWrztseUqCTQdKTeDUcLB9vSAtLcDMIwoKLGJXDP45qRpnsdl6rFInRsSC6Xx9SHcbGZp59+Gvfffz/q6+tRVVWFZ599FqdPnxa8EPlHP/oRPvShD+ErX/kKTp06hUceeQQ33ngjxsfHUw42zIbS20JgqzL97Nmz0Gg0UCgUeRcVIDOLZWlpCS+99BLa29uxb9++rC5U7kbORFi4eEp0zU2hb+4f/WkUvdO/wp3uZ/HBzRdxp/tZXOP4LX4XXIayUQnd/q2gZ8ehDsjMyW8qU0iJ1y+24s6ZTty62Ioea83WbPeInI/X2BUBfj48Z10wXJyKpvHsW29NKSoMTcN2zetBop5vu+tBRHTVoEBBeqkZZZNHB1dVGSaO78t4AJcuiSsMSDEXBQCT4neTsCyCkxN4dOEcLmwI10uOpmlIJBLI5XIolUooFApIpVK+UJdhGEQiEYTDYd4S3o5SDd7nsk8Yl+zS0tKCu+66C+Pj4zh48CD8fj/e+MY3oqKiAu9617sEPeZnP/tZ3H333XjXu96FvXv34itf+QrUajW+8Y1vCHqcorFY0lnwosfxdnV1QavVYmBgIA9nl0g6FgvXmn9lZUWQ1OdMhCWfQfp0+fn0NP569OmEZKV11o97Nl6AwWzEbZSRfzwYDILOcu+zrI7NNONcZ7Xz61hp2kqWeNmPfp7QOfhrn7ofwbYjeK33DHy3fYSfP8+Ya/h58Fx/rQavDhRBypYuXS4T2jYN+FXDXMK/JXWFIflcFJam8aXHHobctoK7P/pIjGuMi+uwIPi1dQJHdLU5+Z1jrJlL2YmcoETfB9HPixaSUnaF5eu8q6qqIJfL8da3vhUf/ehH0dfXh6mpKcHePxQK4fz58/joRz/KP0bTNG644QacPXtWsOMAJWSxMAyDwcFBzM7O4vjx42hoaCjIwC2OnSwW7ke02+3o6ekRpJ4m3VqW6Juee12hRYVhWfzdn/6UNAOWe+wjz/8JDMuC8TOYe2YOtDf7y3NJ7Ul4zFVVhtGTe/kmkdFWDEPT+K+Pvwfjr3g5TkeObj1mrkGg+zQYcw10T/8YjR++BrUPvx2NH74Guqd/DAUrQac7Ma5nYBR4q6MTb2a6Ia1ObplI2eSfjaszibaW7Hc9hJukb0Bl91vxP3//tzHnHB3XsYd98DKJ7fCFZidrJhKJJFgzpewKy2ffLm56JE3TOHbsGN7ylrcI9t42mw0MwyRkoVZVVWFtbU2w4wBFZLFsh8/nQ19fHyQSCXp7e/l0vUIKy3YWi8fjwYULF6DVanH69GnBRqamIyy5qE8RgudWV7HsTV2vQgAse7347W/HccxrghK7c2+uaBKFJeGc4hIAGmWd+KvJOqiZ2N8rYR58VFbZ9dIGhBQsZlRuqGkZThsb8OryDijorfcIrlsBZ+Kxf9Y0jTPrtThhq4IkLoISXWfCW0sADjnLEen6AJ56/PWYmH0Rc7WJXZQ3IkFopanTWXNBvDUT/Ye7RyKRCL8ZKyXLJZ8WCyBmhQlOqgXQZrNhYGAgptMvB3ehCtmzK11SWSwWiwWDg4NoampCe3u7oOe1k7AUU+ZXPGs+385PArBs9+CYcnfZfWGKhUMR3PmJ2LJi5Pp6vHalEXW+5Nk4SeMel/p7ac01eNNUO6QtchiOGiCjY3e31GIkaUZwhCZ4umYZY0YHblpuQo0/djHh6kzikRIaJ/afhnZfO/qXzyf8+2YeLJbtSOYys9vtsFgs2LNnD59EwgX+8zlALhtYlhVsY5gOuZweWV5eDolEgvX12OLb9fV1VFdn3/8uGUX7ixJCMDs7i76+PnR2dmLv3r0JF2B0a5V8E2+xEEIwMzODgYEB7Nu3Dx0dHYIv7NsJC7c7LJSoEELAeJmYwVPRVKfZVbqSVkLiWIVy7CwkjtWdX5AESZq1MXKGxvUrDbhrqjulqACX4x7RcP29OCKzIYRnY8WMjbBoWNZAwUhgXLdjz4vDMK7bY/7fovLju20X8YeaRQTpna9jiUYCqVYKXQqrZDOyc2ucfEHTNBwOB4aHh9HZ2Yna2lrIZLKY+zYSiSAUCiESieRt3lAm5NtiyWUTSrlcjmPHjuGJJ57gH2NZFk888QR6enoEPVbRWCzRRBcRnjx5EgaDIenzuAs0EonkdVfBHZu7ERiGwfDwMJxOJ06dOgW9Xr/Dq7MjmbBEz1DhnpNPUSEMgeMpB8KOS611JIBuvw7qdjXYEAsSImBDLI5LzKhVa7Dq8yaNs1AAqmkVrjv7R1R9K3FWSCbQoHDQUYZBsz3lc/Y6zXjlWj20kZ3TOZP11+L6e0WzObgJZaMStGwr1hCyhGAKKXDPNwbR8dgnL9WlbH1amsRW9J8vt2BC78QtS81o8qa+fig5Bf+iHwpp8t94M5KepZYP1tfXMTw8jP379/N+/WQJANw1XIzWTD5jLIQQPsaSKz70oQ/hrrvuwvHjx3Hy5Ek88sgj8Hq9gmefFY2wcIshNzlRJpOhp6dn2/YH3MVXSIslEAjgwoULoGl6x/MV4pjRjToLWUlPCMHGhQ34Z+K67jLA5sAmNgcSJzl+Qrof9+AFUIhtY8Utkf+EBlR96w1JYxnJ3ELbcf1KA4whBYaNdjgVQT5zqzygxC1DCjTPriFcpUz7fZPFPeIhEYKwLbzVtn/Ui7AjDIljFS1f/eeouhTwnz6+on9THsaPm6fwjukuVAWSW3gRZwTuc25EKBbYn+Q8C+wK41hZWcHFixdx8OBBVFQkznnhrtX44szoFPno5+6mOHM35NtiyXVLl7e85S2wWq34x3/8R6ytreHw4cP47W9/m1ZbqUwoGmEBAKvVisHBQdTW1qKzszOtHzSdYV+5QCKRwOfz4ezZsygvL8+6PiUToi2WQgXpCSEIrgThet6V8WtvVNThUf0p/HNgCKuhyzGXalqFf9AexOtmFradVZIJciJBj7UGPdYahCkWTkUAFKHQ8vtfoOJbn8jKIkoV94jG+afYSH2y+Ew0Cf3FaIKXytdxy1JLytcAgIzQkDM0QpLY994oAotlcXERk5OTOHz4MD93fie2S2febXHmbsh3skE+Ku/vvfde3HvvvTk9RtEIi8fjQX9/P/bt24fa2tq0X1eozDCv1wu73Y7Ozk7B572kIrrDcSHiKYQlcL/oTjrdMV1uVNThBnkt/qRahm0ziDqi5yvvw1V04qwQiSQmlpENMkKjMqCGxLHKiwqwO4soXZLVpUSTrL/YlM6d1lRMTUSGkCRWSArtCpudncXc3ByOHj0Ko9GY1Xuksma4a58TmXy4zHbb0iUTQqEQQqFQTl1h+aJogvc6nQ7XXHNNRqICbD/sKxcQQnDx4kU4HA6YTCY0NzfnbWHn3G+FCtIHlgO7EhUOCUXh2kA97pC14bS8gm/nklDDIZEg9IUvoOyvDkGi2f3NvV12V65I+EwUBXLp8xKaxi//7t6ElOGAlMFSGunSLJW7UdWZQgjB1NQU5ufncfz48axFJRk0TUMmk0GhUEAul0Mul+ctASCfMRbvpXR8sW2+wGQTn8inxRIOhzEwMAC/34+mpib4/YlT/XIFl1K9sbGBiooKyGSyvGd+RZy5F/CYWSH7OqC/fS9oioLhlAGOJx27eu9k1kN8dlc8ygYlDKcM8E35tuJGWazlnle/BZFXXA8yOskfi4vV6Bo0ACYTXjOpd6HRm3rnyoBgQ5YYTzHLcjODaDsIIRgfH4fFYsGJEydyGiPYyZoROgEgnzEWTljEOhaBEWqKZC7gih41Gg1Onz6NlZUVeDw77yqFgHN9VVZWYm5uDs888wzMZjM/mTFXA83iSdW7S2j4WAYLyGb9ULeqIS+TQ9OtgXcsvaFgqd43neyuaAKLAQQWA5BXyzMWlYgkAv1ePXQdOlCSKmz0NYCZ8vHnAgANXgIFI+H7mnFM6Vy4brU+6awaSk4h3EInbSVTJk8vrVsoCCEYHR2F0+nEiRMn8nYtcqRTnBn9vExFIt8Wi1qtLngmnBAUlbBkQz6C91arFQMDA2hoaMCePXv4VuP5yLuPjqfU1NSgtrYWPp8PVqsV6+vrGB8fh1arRUVFBSoqKqDT6XI3uKtOAXWbGr7p9IodhWCzfxPyKjmkGikUVYpdCQuQXnZXMkJrmWVbrWpW4ZQ7YfKbUL6yNTVTf0QPZb0SgeUAQpYQIu4IJKDQuqnHmDE26O9ShGBTBFARVEHZoIS8Ug55pRwSjQQURW01nFxJPG65LH+7XZZlMTw8DI/Hg+PHjxekGWw0yRIAOJHJ1prJp8XC1bAUU2FztpS8sOTSYiGEYG5uDlNTUwlJBdkO+sr0+MniKWq1Gk1NTWhqakI4HObH/87Pz0MqlfIiYzabBb0pKIqC/qge6k41wrYwws4wwq4wIs4ISCQ3/n7CELj+5ELZq8vgXxDG9ZhOdtduYCUsDr7qIAKhQMwGQKPRbP02LRUoO1QGNsgiZAlhn6MSY0n6vkzpXagOaGE8bYx5vH9jFT9YHUx67HxZLFzvvmAwiOPHj2c94z1XCJXOnM+ssCtleiRQZMJSTK4whmEwMjICu92etEgz1xZLtEm/XZBeJpPFzN52Op2wWq0YGxtDOBxGWVkZLzRCzXaQaqSQaqRQNW25Pbiq+4gzAu+E93KxpEBENiJY/5/1nZ9YJNAMjYg9Ak21BhqNBs3NzfwGwGaz8XVP3ETNQ3vr8dOZCTBxvrZzFWtQBGV4NWuGlJaAEIKnnbP4uWUs6XEloFCRB4slEomgv78fLMvi2LFjgs8MyQWp0plTdQGPHlORL1eYx+OBWq0WLZZiIBdZYYFAAH19fQCAnp6epCZ+riyW3VTS0zTNjzzt7OzE5uYmrFYrFhYWMDo6CqPRyIuMOs0WK+lAURQkGgkYHyO4qAgFJaMg1Unzdn60InaXG78BcLvdsFqtmJuYg2RDgtpaFRZVsS7GkITF4+p5XJi14rUVXZjxO/Cscy7lMY8b6qCU5PaWDofDfEPYo0eP5r3jhRDsZM1wIhOd1pwP8lHDki9K76qIQ2hhcbvduHDhAsrKyrBv376Uu5VcWCzxF/ZuBnNxQ4T0ej3a2toQCGy5ZaxWKyYnJ6FWq3mRMRgMu7p52CAL51knwtYiFRU5BVWTCr7J3MSG4ue0UHIKEl3qXS5FUVAH1aiwVkC/rgcI0O2wYrFuIenz7WEfvr1yYdtzOKCtwu1VSUrxBSQUCuHChQtQKBQ4ePBgXtvJ55JU1szKygokki1LMRQK5bw4M5d9wvJNUQlLNoubVCpFMChMUdjKygpGRkbQ3t6+Y32K0C646HhKLi5cpVKJhoYGNDQ0IBKJwG63w2q1oq+vj3fLVFRUoKysLOMFwzvuLVpRAQDDCQP8c7lJDdc9/WO+pX50Jb/lpxYoG5TQHdRBor60Mw6y8M364J/xg/HGXjv7nGacq1jDhjzzlizXmJpxa+VWWnau4FoXabVa7N+//4rIXEoG97mWl5exsLCAI0eOQC6X56U4U4yxFBFCLPCEEExOTmJhYQGHDx9O2tsonnSHbqV7/HwWPUqlUlRVVaGqqgosy8LlcsFqtWJiYgLBYJCPy5SXl6dVWySEe8lwyoDgalCQAsxo1O1qKGuVCDvCCC4LW5W+3ZwWxlyDwGIAbIiFtlsL34wPgaUAkOKSkRMJ/nK6C4/XLmDC4Err+BSAWyv34lrz9u1fdgvXv89kMmHv3r1XRAxgOxYWFjA9PY2jR4/GxFbTmZy5G5ERXWFFxG6FJRKJYGBgAF6vF6dPn077hxXKYin0DBWapmE2m2E2m7Fnzx54vV5YrVYsLy9jbGwMer2ed5mlSoUMJSnUywSJYxX43QVoTncjAOHqIKQGKXQHt4oMlfXKtFKVFXWKtAVou0p+LusstB6CYz29wk5tRIY3LLRhUb2JJ2uWsKZO7bqjCXByU41qKgQHHDAajTlzz1y4cAFVVVV8qv2VTCpRAWJjM9GdzYUqzhSFpYjYzQLP7cQUCgV6enoyym7hOg1nm47IBem5i7IYBnNRFAWtVgutVouWlhYEg0E+lXlmZgYKhYIXGaPRCIqiMDc3h7nAHLrRndUx411JoSza5KfC9DITGD8D70Uv/PPpucIi7gj0x/TYOL+x43OzqeRPhwafDu+Y7sJ40waeMi/BHTdjRU3L8K7ao9D7WNhsNgwNDYFlWZSVlfGZZkKk/25sbODChQuor69HW1tbwa/PXLOdqMTD3fNCFmeKwpIjsrlwsxUWu92O/v5+1NXVYc+ePRmLQ3RGSaavFTJIn0sUCgXq6upQV1cHhmH4uAy3kMlkMoTDYew37AfjyPw32MmVtKtzr1Ngc3ATgcXMXGuMh0FwLbhVyLi0/WuzqeSPh1bSUNQqtmJAUcYPBQonzPXoaW3DU45ZPOecxyYTRJvKjDfXHEClXAtogcrKShBC+AzAxcVFjI6OpmVpbofL5UJfXx9aWlrQ3Nyc0WtLkfn5eczMzKQlKskQojjT6/WivLx89x+mCCgqYcmGTLPCCCFYWFjAxMQEuru7UV9fn9VxuYuCYZiMUi4LOUNlN0gkElRWVqKyshLhcBgXLlyA3++HOWIGM7uDqNBIGltIx5WULbuJp2Ty2mwr+eVVcqjb1FDUKEDRW1lrrhdcYH3s5X9vVoOiKdxY3oFXlbWDIWzC6GMgMQMwGAzCarXCZrNhZmYGcrmcj5mZTKYdkzMcDgf6+/vR0dGBhoaGtL+LUmW3ohJPtsWZosVSRGTS0oVlWYyOjsJiseD48eMwmbKfrR69O0mXaFEpVitlJ7gaH7lcjjNnzsDzkgdBV+qFWH9Cj4gjkrQNTK5cSekir5IjtL77wViZVPLLzDIYThkg1cbeevJyOSpeU4GQNQRKSkFmjm0ySlMUaCq9bD2FQoH6+nrU19eDYZiEotnoPnPxyRmcRdrV1ZVxp/FSRGhRSUa6s2bW1tZyOigwnxSVsOTSFRYMBtHf3w+GYdDb27vrvkacKZuuqBU6SC8Em5ub6OvrQ1lZGbq7u7c+hyT155jTz6F8sBzaYPJdmBCupGxQNiih6dJAZpSB8TGwPW4DCeWnBb3hRKKocFA0BUWVsAuLRCLh4y6EEHg8HthsNj45Q6fT8anmXq8Xo6OjMaOEr2Tm5uYwOzubU1GJJ5U109fXh/Pnz+PUqVN5OY9cU1TCkg2csHBWQDK4IKTRaMSBAwcEK+xKp0gyOkhfyqLCBYmbmprQ0tLCfwZNh2YrRTjJuty80bzj+2brSkqGrFwG/RE97I8nn3Uvq5BB3axGcC2IzaFNSLVSEIbkTVQUdQpI9YW75SiKgk6ng06nQ0tLC0KhEJ+cMTc3xycAcBumK6UAMhmcqBw7dgx6vb5g50HTNIaHh3H77bfjk5/8JB544IGCnYuQXBHCwi3eyRbstbU1DA0NobW1Fa2trYIu6jtZLKn6EJUaS0tLGB8fx969e1FTE7vwy0wy6A7rsNmXOON+O5xyJ9SUGvKQfMuVVF4DiqZA0Vs7d9Dg4w07QgHafVpouraC1KoWFfyziVlgYWsYbqub/3sI+ZkPT6toaPZooG7Pb0v7nZDL5aitreULZjs6OhAMBjE+Po5gMAiz2cxbM4XuXCwkxSIqADA6OorXve51uP/++/HAAw+U5PqQjKISlmxdYcBWPUp0iiU30W5ubg6HDh1CZWWlYOcZfexUFkupBumjIYRgenoai4uLOHLkSMr55eo2NSKuCL+Yx7c4iYEGjKeMMJQZYLPZsGpZhcPhgFqz1WKmsrISer0eFEXB8axjx3b1Eq0EhlMGyM2Xf3vdQR2CK0GwwRw1CU2RjBCP1CSFZo8GynrlllgWIdwo4WPHjvHuoD179iSMZuA7M1dU8L9PKVJMojI+Po7Xvva1uPvuu/GP//iPJfudJqOohCUbokeUckQiEQwNDWFjYwOnT5/O2QzpVBbLlRBPYVkWIyMjcLlcOHHixLbZKhRFQX9Ej8hGBMqffC9pixMO8yvMkJdtiQDXYiYcDvOpzFzn30pTJUzr2ydXqFpU0B3WgZbGijYtp6HZq8nYikqX8u4AAr8fgl9Vn9J1Z3q5CfIqedH+9tzGa2VlBcePH4+5RyiKgkaTXmfmsrKykmlEGS2ihRaVqakpvPa1r8Xb3/52PPjgg0V7nWQLRTLtU59DuGZvmfL444+jp6cHWq2WL3qUy+U4fPhwTudEnD17Fi0tLaiuruYfuxJEJRQKYWBgACzL4vDhw2lnqrDTC9Ac2peQ5bXwmWfAmGugO6yDpmP7Xkhcixn3iBuqtRRV+BRgPG2Esv6ye4YQgpAlhMBSAGFbGJGN3IxRTtUbLJrKWytBy4vXQuVGCVutVhw9ejSj/lTRnZltNht8Ph9MJhOfZSZk12whmZ2dxfz8PI4ePVpwUZmbm8NNN92EW2+9Ff/xH/9Rkt6MnSiqrUa2izAXwHc4HOjr60NNTQ26urpy/oPFWyzpzlApZnw+H/r6+qDRaDJOdJAuzW5blyJR7fxeNE3D7PNB+vwggsbG5BYBAVxnXWmfl1CkW9BJIgQorrlXPFzKvcvlwvHjxzMeJUzTNEwmE0wmE+8y4xIAJiYm+K7Z5eXlMBgMRbFocqJy7NixnHkv0mVxcRE333wzbr755itWVIAiE5ZskUgkWFlZwdLSErq6uvJW1BWdkXYlBOldLhf6+/tRU1OTVV8o0t4OQtOg2OR1KelkREm+/W3I770Xaja5RbBt/CbHpFvQGXaF+Y7GxQTLshgaGoLX68WJEycEqZlQq9VobGxEY2NjTNfsgYEBAOAbmpaVlRVkIFgxicrq6ipuueUWXHfddXj00UevWFEBilBYMp0iybIsIpEIlpeXcezYsZQB5lzACUupB+kBYH19nR8Z0NjYmNV7kLo6hL74Rcjf/35QDBNTl6LuUO8oLNTyMuT33ssLU7xFkMwNtZWqnB+hSbegsxgD9QzDYGBgAKFQKGejhKO7ZhNC4Ha7YbPZMDs7i+HhYRiNRj7LLB/t4YtJVNbW1nDzzTejp6cHX/va167oVG6gyGIswJZ/P91TCoVC6O/vh8vlQmdnJ5qa8lOxzTEwMACNRoPGxsaStVK4FjfT09PYv3+/INlz1PIyqOlphMobEVJWQmqUQmaS7fj90E8/DeXNNyc8vvLA9xCuakLjh6+JXdQpCiAABZIy3iE0uqd/nFDQGX1MmVkG83XmoroWuFHChBAcPny4IJaD3+/nEwAcDgeUSiUvMrnozDwzM4OFhYWiEBWr1Yqbb74Z+/fvx/e+972SSXbYDSUrLJubm7hw4QL0ej3C4TBqa2uz7vuVDVwAdGlpiU+TzWZIViFhWRYTExNYX1/H4cOH81Z9nApqeRnKrq4EV9rCZ56BbH0OtQ+/fdvXRycK5JItd1xsQSetoqFqUUHbpd22G0G+4fq6SaVSHD58uCiuT66hKRebEbozczGJit1uxy233IK2tjb8+Mc/LoioF4Kik850XGHr6+sYHBxES0sL2tra0NfXl5P586ng4iktLS0oLy+PGZLF7cIqKiqK+iLiUrL9fj9OnjyZcRA3F2znSgOQ4IaKR6gGljvB9QajJBSU9QqomlSQVxZfanEwGMSFCxegUqlw8ODBonHTRjc0FbozczGJisvlwq233orGxkb86Ec/Kur1QGiKzmIJh8PbFh3OzMxgZmYGBw4c4NN8BwYGoNPp0NramtNz4yr8ORGLbiRJCIHX64XFYoHFYoHH4+HTMCsqKopi4eYIBoPo6+uDVCrFoUOHiu6Cp5aXET43BqetLEYkttxQnwDFMiAUDRACKqqXTL4sFg7tAS20XcXZjTYQCOD8+fPQ6/XYt29f0YjKTgQCAd5lZrfbM+rMzBXzFoOobGxs4NZbb4XRaMTPf/7zK6pzQTqUjLAwDIOhoSG4XK6EXPTh4WEoFAp0dHTk7LziK+l36k7s9/thtVphtVrhdDqh1Wr5XVo28zGEwuPxoK+vjx8zW6wLzkbfBnxTiR2RKw4G4f35AIKGBqiHnt023pFrpCYpym8ovvkZPp8P58+f55uFFpsllS7RnZltNtu2nZk5UTl+/HjBW897PB7cfvvtkMvl+NWvflVUm8p8UXTCEolEEtxafr8ffX19kEgkOHLkSIIPdmxsDADQ3Z3dFMOdiJ6jkE2QPrrZn81mg0Kh4EXGYDDk7ca32+0YHBxEY2Oj4H3ThMb2exsi7tgiR1pBo+J1FVj/n3X+sWTxjnxBSSlUvaG4ugB7PB6cP38eNTU16OjoKOrfOBOiOzNbrVZsbGzwnZmDwSDW19d37BCRD3w+H+644w4QQvDrX/+64OdTKIpeWJxOJ/r6+lBZWZlyhz0xMYFQKIT9+/cLfj5CV9JzgUuLxQKbzQaKovjgv9lszpkFsbKygrGxMXR3dxf9nA02xMLyc0vSf6PVdPrNKXOMvEIO8yvyl96+E1wX74aGhqLfOOwWbrM2NzcHr9cLmUyGyspKVFRUwGw2FyRJwe/34y1veQt8Ph9++9vfFrzCv5AUXfA+mqWlJYyNjaGzsxMNDQ0pbxSpVAq/P72Z5pmQi0r66MAl177EYrHwQ5jKy8v5DDMhYh9cXGphYWHbRpLFxHbtWFKJSr4LJ4mCQLm/ePzmV9soYZlMBr/fj3A4jFOnTvH9zArVmTkYDOLtb3873G43Hn/88ataVIAitFgYhkEoFML4+DhWVlZw+PBhlJWVbfua+fl52O12HD16VJBz2C5Inyu47BiLxQKr1Qqv1wuz2czvwrKpkubadzidThw5cqRkzPKwK5xypkoy0unftRskWgn0h/Wg5BR8mz64Nl1Y965j07MJg8EQk8FUCOx2OwYGBrBnz568ptwXCq7rNlcUHX1dE0L4zsw2mw0ulwtarZYXmVx0Zg6FQvjLv/xLLC0t4Q9/+MOO69XVQNEJi9/vx/nz5xEMBnH06NG0mtotLS1hdXUVJ06c2PXxMw3S5wqfz8eLjNvthl6v50UmnQUsHA5jYGAAkUgER44cKamRp4QQOP/k3LFlPrBlqSQUTgqQHUbJKFTcXAFKSgFU8j52XAaT1Wrli/44kTEajXm5biwWC4aHh9Hd3Z0wK+dKZDtRSUZ0Z2abzSZ4Z+ZwOIy/+qu/wsTEBJ588klUVFTs6v2uFIrOFba0tASJRILTp0+n/aNLJBJEIrvvZhsdT+FGDxcKtVqN5uZmNDc3IxgM8hlmU1NTUKvVvMgk24FxyQ4qlQpHjhwpiqK4TJGoJWm5t9Lt35Up6jb1jh2KlUplzGz5+D5Z3C45V63lV1dXMTo6igMHDuRk3lCxkamoAFsus5qaGtTU1MR0Zp6ensbQ0FCMyyzT7K1IJIL3vve9GBsbE0UljqKzWCKRCCKRSEa7PYvFgomJCbzsZS/L+ril0u6ea/THBf+5mE1FRQVMJhM2NzfR39+PqqoqdHZ2Fu3nSMn8Ery/7Af9Uh/M//2v27u3KMBY64HpxiNx1foSLHzm6eyFhQYqbqmARJmdIBNC4HK5+M1AIBCISZMVwue/tLSEiYkJHDp06KpwvWQjKjsR3ZnZ6XRm1JmZYRjcc889OHfuHJ566qmiT4jJN0UnLCzLIhwOZ/Qau92O4eFhXHvttVkds1REJR6WZeF0OvmiTO4zVFdXo6urq+R6EtH/+U0oPnDflpgAiP4Vkrm3lI1K+Kp9CH/hK+j66ucEq2dRtalgOCpcexuv18uLjNvthk6n411mWq024+ttfn4eMzMzOHz4MEym7YehXQnEDyXLRSwrujOzzWYDkLozM8uyuP/++/HUU0/hj3/8Y9ZNW69krghhcbvdOH/+PK677rqMXleIIH2umJ+fx9TUFMrKyuD1ehEIBPgbo6KiIqcDz4SAGZ+H9tj+bVu2rDzwPQS6Tyf9N8HqWSig/KZySLW5EeX4miauspyzOLfbJXMZftyo6EL3dssH+RCVZMfkOjNziTRarRa///3v8frXvx7f/e538etf/xpPPfUUWlpacn4+pUhpbWlTwLWvz4T4GSqlKiqEEExMTGB1dRXHjx/nFxuPxwOr1cqnbBuNRr5eJl+VwNTyMqipqa05LXV1KZ8XsoXg/98+6LYRlWTt6aPh+nelA62iwfqTH0t3UJczUQEAuVyO2tpa1NbW8sPprFYrhoeHwbJsTFwmepdMCMHk5CT/O5dKht9uKISoAFtrgdFohNFoRHt7O/x+Py5evIinnnoK//7v/w6JRIK3ve1tmJ2dRX19fdG1RCoGik5YslncowdupfP6YgrS7wauzY3X68XJkydjMui0Wi20Wi1aWloQCARgtVphsVgwOTkJrVbLi0w2rph04AZ2ce6p0Be/COauuwBECY5Wi/CMA25XOVDWlDjrBFvuMEJLYHv3g7uuT5FXymG+dquOh7AEbIAF42PA+BiQCIHMLIPMmL9FQiKR8NYKIQQbGxuwWq38/JLokb/z8/Ow2Wx5XWALSaFEJRkqlQqHDx/GNddcg+XlZTzwwAPo7+/HnXfeiWAwiMXFxaviN8mEonOFZTP3PhQK4cknn8QNN9ywY1yhVOMp8QSDQfT394Om6YxmbHDplxaLJabJX2VlpWApsknb30skCIyNgf7DHy4LDi4Jx6XgPICY3l+ON34EwZYDCNc0gZbTkMzP7KoAklbQKL+pvKjn0XNwveYsFgucTicoikJ9fT1qampyUotRTHCisrq6imPHjhV80SaE4NOf/jQee+wxPPnkk3yHD5ZlcfHiRezdu7eg51eMXBHCwrIsfv/73+OVr3zltvUaV4qocI0kjUbjrjrXcq4Yrl4GQEx7mWzTlFMN7Ap+97uQ33VXjOBwcMF5AAmxknQLICk5BbXUAaVvEe7NSkR01QnP0XRpoDtQ2M636RI9SrihoQFOpxN2ux0SiYR3mRWqfUmu4Fx+a2trRSMq//7v/47/+I//wBNPPIHDhw8X9HxKhSvCFcbFR1LFWbggvdDtWQqBw+HAwMAAGhoa0NbWtqvPEe+K4drLXLx4EeFwGGVlZaisrER5eXlGfmTS3g5C0wkWCyEkqagAl2tPAt2nYywSiWOVFxUgcVwxJIC6VQ1lnRKqX34f8vveD4ploU+RHRZYDpSEsHCjhMPhME6cOAGZTIaGhgY+E9Bqtcb8TpzLrNiTNLajGEXl85//PB555BH87ne/E0UlA4pOWLKBoqiUAfz4IH0piwpXENfV1YW6bYLh2UBRFEwmE0wmE/bs2QOPxwOLxYL5+XmMjIzAZDLx9TI71WEkDOySSBD6whdATp9OEBz+NSmC8zsVQBqOG6BqVIFaXuZFhXtO+bc/cVmALsF404/FFYpIJIK+vj4AwLFjx2LcuzRNo6ysDGVlZejs7OSTNLghWcXQYiYbokXl+PHjaXXcyPX5fOUrX8HDDz+M3/zmN4J09biaKEphSWeKZDzJhCW+PUupBukJIZidncX8/HxavdN2C0VR0Ol00Ol0aGtrg9/vh8Viwfr6OsbHx6HT6WJmyySDuesuBG64AdT0NEhbG58VFiM44ILzsZMiowlXNScG9TkRogCZacuSoqamEgSLYpnECnwWICECSlGcwhIKhdDX1weZTIZDhw5t6+aK/p1aW1v5FjMWiwVTU1NQqVT8ZiCf4xkyhctsXF9fLxpR+cY3voH/9//+H371q1+hp6enoOdTihRdjAVIf+59NM888wz27t2L8vKtwUtXSjyFZVmMjY3BbrfjyJEjBZ+MFwqF+GI/u90OpVLJi0y6QWVqeXlLcDQaBMftcG9WbBuQ35ocmTjQS7NXA90+Hf+eCQkDFI319z6CYMdRMOYavkWM9k1HIO1u3vV3ITTcKGG1Wo0DBw7saiOUrOCPs2TKysqKJi5TjKLy3e9+F3/3d3+HX/ziF3jlK19Z0PMpVa4YYXnuuefQ3t6OqqoqEEL43mGlWp8CbGVwDQ4OIhwO4/Dhw0U33jR68bJarXzMprKycsdiPwAgDIHl/ywg4Z1/6/gCSHW7GrrDupjfVvLtb1+2hi49ThECQtHw9N4G7fM/20oAiEt/Lga45qtGo1HwyZ7RPbLiW8xk2zlbCDhRsVgsOHbsWFGIyo9+9CPcd999+MlPfoJXv/rVBT2fUqYohWW7ufepeOGFF/h0zCshSB/dSPLAgQNF354lOqjMtZeJni0Tf/4hWwiuc66UhYo7UfXGKt5lGnaEwXgZ0DIaEtcqZEPnoXz3O2OtF8S1iLmU/rxd4Wa+8Hq9uHDhAsrLy9HV1ZXzazZVi5l8js0uNlEBgP/93//Fe9/7Xvz4xz/GLbfcUujTKWmKe7XKAJqmY6ZPlrKobGxs8FMzOzs7SyI2FB9U5or9pqenMTw8zO+QKysrIZPI4PyTMy1LJRkS9ZYbh7AE7hfdCCwEov5VCeVFGrXx8Za496AYZssdV2Bh2dzcxIULF/I6Slij0UCj0aC5uTnGtTk7OwuFQhHT+j8X114xisovfvELvPe978X3vvc9UVQE4IqwWAgh6OvrA8MwaGpqgtlsLllRsVqtGBoaQmtrK5qamkr2c0TD7ZAtFgs2NjZQKa9E5fru2ryXvaoM/hk/fNO+hH9LOqMFxWexuN1u9PX1obGxES0tLQX/raNbzFit1m1bzGQLIQTj4+OwWq04fvx43toLbcdvfvMbvOMd78C3vvUtvOlNbyr06VwRFKWwxM+93w4unXhjYwOLi4t8oR8XUM7lHHmhWVxcxOTkJPbt24eqqqpCn05OCAaDsC5ZIb2wO2OZVtBgg6k3H/EBf0/PbdCe/dmlv0sQ+uIXChpjcTqd6O/v5zcQxUZ0ixmuESPXYiab2SXcexabqDzxxBN429vehsceewx/8Rd/UejTuWIoaWFJlvkVPUee8/VzLphiyoaJhsvh50YxG43GQp9SznE840BoPbMOC5kSH/Dn/k4d7IDhtn05PfZ22Gw2DA4OltQoYW52icVigcvlgkaj4VOZdTrdjtZWMYrKM888gze+8Y344he/iLvuuqvgFuOVRMkKCxeg3y5Iz7W/5kQmFArxAeXy8vKiCIgzDIPh4WFsbm6mPYr5SmBjeQO+5xPdWPlAd1gHTUdhigctFguGhoawd+/ekh0lzPWb41LOozs4JPMQFKOoPPfcc7jjjjvwmc98BnfffbcoKgJTlMLCMEzKUcPZzlAhhMDj8WB9fR0WiwV+vx9msxlVVVWoqKgoSOvrUCiE/v5+UBSFQ4cOlXQ7jkzgkhNa/a2QbQr3vWsPaEFLaTD+rY7FjJ8B62PB+BngktdM3XYpTZnO/0KyurqKsbEx7N+//4oZJRydDWi1WhNazMhkMly8eJHvzFwMovLCCy/gtttuw0MPPYR77rlHFJUcUFLCImQlvdfr5avJPR4P37KksrIyL3n9Xq8XfX190Ov12LdvX1G66HIB1+uspaUFdcY62B+3Jz6J2ur/FVgJpJ2OrOnWQLc/efEoIQQkTEBJqYIICnB1jBLmNm9cU1OPxwOZTAaWZXHo0CGYzeZCnyIuXLiA173udfiHf/gHfPCDHxRFJUeUjLBwQXqGYQRPJeZallgsFrjdbhgMBl5kcrHDcjqdGBgYQF1dHdrb26+ai3ttbQ0jIyPo7u7mZ4S7zroQWArs8MrtiRgiMPQYcjZbZrfMzc1hdnYWR44cuSriZ8DW/To8PAy73Q6tVguXy8XPlC9Ui5mBgQHccsst+MhHPoIHHnigKK+VK4WiFJb48cT5bM8SDAZ5kXE6nWn1xcoErpFkZ2dnyQRuhWBhYQFTU1M4cOAAKioq+McjGxHYfmfL+n2JmsBSbYHNZYNCoeB/q2LojRU9Svjo0aPQ6/UFPZ98QQjBxYsXYbfbcezYMahUqoQWMxRFxaQy59piHxkZwWte8xrcd999+Id/+IeCXxtXOkUvLOkE6XMFVzzGDcXiMmGymbxICOF3rgcPHuR7ml3pEEIwPT2NpaWllBlvrj+7EJjP3GqhFTTM15kh1UrBMAzsdjssFgu/cEXPlsl3yjlXBMi1gL8aRgkDW597bGwMDoeDF5V4uMxNLi4TDAZz2mLm4sWLeM1rXoO//uu/xoMPPiiKSh4oWmHh+oVlGqTPFZFIhBcZm21rd1xVVZVW80Vu0pzNZiuKRpL5gvvcXAPNVItrxBOB7be2rSrGNJDoJJBXyKHp1CSdTx+dcs4FlPOZDcgtrtyO/WrJ9IsWlePHj6fV244QEtNiZmNjA3q9Pqb1/27u+8nJSbzmNa/BnXfeiYcffjhvG4xPf/rT+MlPfoKLFy9CpVKht7cXDz/8MDo7O/nnBAIBfPjDH8YPf/hDBINB3HjjjfjSl750RdSwFaWwMAyDQCDAB+kLLSrxcLvj9fV12Gw2SCQSVFZWoqqqKmG8byQSweDgIILBII4cOVJ0jSRzBcMwGBoags/nw9GjR3f83N5xLzYHN2MfpACpXgqpUQqZSQaZSQapQQpalv7iQAjB5uYmLzJerxdms5mvwRB6d8yyLEZGRrCxsYFjx45dNb93NqKSjGAwGJPKvJsWM7Ozs7jpppvwhje8AY888kherdabbroJb33rW3HixAlEIhF87GMfw/DwMEZHR3mX+vve9z786le/wre+9S0YDAbce++9oGkazz33XN7OM1cUpbDce++9WFhYwK233oqbb765qGd8syzLj/e1WCy8C6aqqgoqlQoDAwNQKBQ4ePBgUdTN5INwOIz+/n4QQnDkyJG0U7lD9hBC1hBoGb0lJkYZKImwv7vP5+NFxu12Q6/X8yKz2xgaJ6Z+vx/Hjh27atLHhRKVeLZrMbOT5bmwsICbbroJN910E770pS8VvPuG1WpFZWUlnn76aVxzzTVwu92oqKjA97//fbzxjW8EsOWy6+7uxtmzZ3H69OmCnu9uKUphGRkZwY9+9CP85Cc/wdTUFK677jq8/vWvx2tf+1qYTKaiFhnOBbO2toZwOAyVSoWOjg6Ul5dfFSnFgUAAfX19UCqVOHjwYFF/5mAwGDNbRq1W8yKT6WaGYRj09/eDYZiMxLTUyZWoJDtOdOt/n88XE5eJPu7KygpuuukmXHvttXjssceK4hqcmppCR0cHhoaGsH//fjz55JO4/vrr4XQ6Y+KOTU1N+MAHPoAPfvCDhTtZAShKYeHgskv+53/+Bz/96U8xPDyMa665Brfddhte97rXoby8vChFxmazYWBgADU1NZBKpUVb9S80XPt3s9mM7u7ugu8SM4HLWuJiaJx7s6KiYsfZMpyFRlEUDh8+fEX+tskghGB0dBROpzOnopIMn8/Hi4zL5cLy8jIGBwdx3XXX4eMf/zhOnTqFb/3/2zvzsKbO9P3fgChuoKwqiKBSsKisBa1osSKIIIkV21pq1Vb9uqF20VqsM844bcfK1aLWqm3HpY5bJRFcALGytAzUUXZQKKhohSZE2QQDgeT9/eHvnCGCK0nOCbyf6+ofnpwmjxje+7zLc98HDvBCVFQqFcLCwlBXV4eMjAwAwJEjR7Bo0SK0tLSo3evj44OpU6di69atXJSqMXgtLO1hThcxIpOTk4OXX34ZQqEQYWFhGDJkCC9E5vbt2ygtLVWz7Ois69/CwoIduLrD0y3j1NsdenOYbnJmeVOlUrFPxg/PPBUKBXJyctjlTj4MZLqAEZW6ujrO95JaW1uRmpqK7du3IysrC4aGhliwYAGEQiH8/f05CzJjWL58ORITE5GRkcG2GFBh4SGEENy8eRNisRhisRi//fYbfH19IRAIIBAIYGdnp/OBjRCC8vJyVFZWws3NDYMHD37kvUx3cnV1Ndv1z1jLcP1L8DwwpoqjR4+Gvb091+VolPZLMNXV1WhubmYtS8zMzFBQUID+/ft3OUpYn+CTqDDcvXsXISEhcHR0xNKlS5GQkIBTp07B2toa2dnZnNW1atUqxMfH45dffoGjoyN7nS6F8RxCCKqqqiAWiyESifCf//wHHh4eEAqFEAgEcHBw0LrIKJVK9iSQh4fHM20Cy+VydibT0NAAMzMzVmT44Kv0JJiGT1dXVwwZMoTrcrROY2MjZDIZJBIJGhsb0bt3b4wYMYI9rNHd4aOo1NXVITQ0FLa2thCJROyhCUIIJBIJJ2afhBBERkbi5MmTSEtLg5OTk9rrzOb90aNHMWfOHABAaWkpXFxc6OY93yCEQCqVIi4uDiKRCOnp6XB1dYVAIIBQKNRKQp9CoUB+fj4IIXB3d+/SSaDm5mb2yVgbXf+a5ubNm7h27Vq39r/qjKamJmRnZ8Pc3BympqaQyWSora3FgAED2KZMvtrLdAVCCIqLi1FfX88bUWloaEBYWBjMzc0RFxfHi5oAYMWKFThy5Aji4+PVelfMzMzYB5Dly5cjISEBBw4cgKmpKSIjIwEAmZmZnNSsSbqVsLSHEIKamhpWZC5cuIAXXngBYWFhmD17NsaMGdPlX/z79+8jNzcXAwYMwNixYzW6vq7Jrn9N0z4/xsPDA2ZmZpzVomuYKOFhw4ap7SUxVvLtG2gZkXm4t0kfaS8q3t7evFiybWxsxOzZs2FiYoIzZ87wasb4qH/v/fv3Y+HChQD+1yB59OhRtQbJ7jDz77bC0h5mnfzUqVMQi8U4d+4c7O3tIRAIMHv27OdaH6+rq0NeXh6GDRum9azyhwctExMTVmR03eOjUqnYk0Cenp68nElpi/r6euTk5MDBwUFtvfxhmP4Lpl8GgJq9jL5t8DOiwjR98kFU7t+/zy4hnT17tsdY5ugLPUJYHqahoQFnz56FWCxGYmIirK2t2eUyLy+vJ4qMVCpFcXExnJycMHz4cB1V/QClUqkmMr169WJFRttPxkqlEgUFBWhuboanpycvBhhdUVNTg7y8vGc+oEAIUUs0ZfJKmGPnfD8RyEdRkcvleP311yGXy5GUlNRjzD31iR4pLO1pampCUlISRCIRzp49i0GDBiEsLAwCgQC+vr5qT5fMabTr1693cOnlApVKxfZeyGQyta7/J/VePCtMKJmhoSHc3Nx4PyBqEubUm7OzM2xtbZ/7fTrLK2FygB5u8uMDfBSVlpYWzJs3DzU1NUhOTu4xMQT6Ro8XlvbI5XKcP38eIpEIp0+fhomJCcLCwiAUCuHt7Y0VK1Zg/PjxWLx4Me+ektp3/UulUrb3wtrausu25HK5HLm5uejfv7/G95L4jlQqRVFRkVZOvT2cA8SnwxpMnsq9e/d4IyoKhQLz58/H7du3ceHCBV4Eh1E6hwrLI1AoFLhw4QJEIhFOnjwJQ0ND9OnTB3/7298QHh7O6yd2Zk+JEZmuuPs2NjYiJycHlpaWGjnwoE9UVVWhpKREJ7PT9oc1ampqON9HKy4u5pWotLa24t1338Xvv/+O1NTUHhM7oa9QYXkCVVVVCAkJgYGBATw8PJCQkACFQoGQkBAIhUJMnTqVF794j6K9u++zdv3X1dUhNzcX9vb2GDlyZI8SlT/++ANlZWVwd3fX+ZNx+1AsmUwGIyMjdvap6SXOh+GjqLS1tWHp0qUoKChAampqt7CV7+5QYXkMTFa3t7c39u7di969e0OpVCIjI4Odydy7dw/BwcEQCoUICAjg1ZHHzni46/9RFvIymQyFhYWcHFDgGj5FCTP2MsxsRqlUsrNPCwsLjfqS8VFUlEolVq5cid9++w1paWlspDWF31BheQLXrl175NO6SqXCb7/9xoqMTCZDUFAQBAIBgoKCeH8EkrGQf7jrX6VS4dq1axg7dmyPejpsn3bp5eXFu0A2QggaGhpYkZHL5WoPBl1pzlWpVCgqKkJjYyO8vb15YfmvUqmwZs0apKWlITU1tdvZBXVnqLBoCJVKhZycHMTGxkIsFuP27duYPn06mykzcOBAXi8lNTc3o7q6Grdu3YJcLke/fv0wbNgw2NjY9IgERCZKWCqVwsvLi/PN86eBSV5s/2DAiMyz/JsxotLU1MSbHBmVSoV169YhMTERqampj+0bovAPKixaQKVSobCwkBWZa9euYdq0aWymDB87sQkhKC0thVQqxbhx41gPs5qaGrbr38bGpstRsXykvf27p6enXgopky3DbP4z/2ZWVlaPfajhq6hERUVBLBYjLS0No0eP5rokyjNChUXLMEFIjN1/cXExXnnlFQiFQoSGhvIiU4YZXBoaGjoMrHzq+tcG7ZeAniZCWR9obW1Vy5YxNjZmRaZ9vC9fRWXz5s04fPgwUlNT4eLiwnVJlOeACosOYaz1GZHJzc3FpEmT2EwZGxsbnQ/UbW1tyM/PR2trKzw9PR87uLTv+pfJZOyApa9+WIyTQEtLyxP/7voK00TLnDAjhLC5MhKJBPfv3+eNqBBC8Pnnn+P7779HamoqXF1duS6J8pxQYeEIpotfJBJBLBbj4sWLmDBhApspY2trq/WBmgmpMjY2hpub2zOdMOqs658RGW0fidUEjKD2pChhpr9JKpWisrKSPWHGxDRw+TMghCA6Oho7duxASkoK3NzcOKuF0nWosPAAQggqKyvZTJnMzEx4enqyIqONTBm5XI7s7GyYmppi7NixXRICpuufyZVhnor5arrY2tqK3NxcGBkZPbOg6jvM/l9TUxPGjBnDJmUy9jLMv5sulwQJIdixYwe2bduG5ORkeHt76+yzKdqBCgvPYDJlTp48yWbKjBs3jjXJ1ETsL2P9bmNjA2dnZ42KVvunYsZ0kXkq1nTfxfPQU6OEgf+JSmfLX3K5nN38r6urw4ABA9TsZbQ1eyaEYM+ePdiyZQuSkpL0PuCK8gAqLDyGEIK7d+8iPj4esbGxSElJgbOzM+tf9jwWKzU1NcjPz4eDg4PW0zW70vWvDZqbm5GTk8Pm5/B9uU6TMKIil8ufuJ+kUCjYvbS7d+/CxMSEncmYmZlp7DtDCMG+ffuwceNGJCQkwM/PTyPv21V27dqFbdu2QSKRwM3NDTt37oSPjw/XZekVVFj0BMZ+ncmUSU5OxogRI9hMmacZKBlDRRcXly659D4vj+r6t7a21vrmMbP0N3jwYLz44ot6d9CgKzyLqDyMUqlUO2HGOGgzy5zPK86EEBw6dAjr1q3D6dOn4e/v/1zvo2mOHz+Od955B3v27IGvry9iYmJw4sQJlJaWwtramuvy9AYqLHpKQ0MDzpw5A7FYjKSkJNjY2LDLZZ6enh1+4RnvKz7Y/QMdu/4HDRrEioym1/eZKGFra2uNL/3xHZVKxWboeHl5dWmW2N5BWyaTPbe5KSEEx44dw5o1a3Dy5ElMnz79uWvSNL6+vnjppZfwzTffAHjwdx4+fDgiIyOxYcMGjqvTH6iwdAOampqQmJjIZsoMHjxYze4/KioKgwYNwsqVKzn3vuoMpuufWd83NTVlRaarzYr37t1DdnY27OzsMGrUKCoqGqL9MqdMJkNTU9MjfeceRiQSYdmyZfjpp58QEhKisZq6ikKhQL9+/RAbGwuhUMheX7BgAerq6hAfH89dcXpGzzkO043p378/wsPDER4eDrlcjuTkZIhEIoSHh6Nv375oaWnBli1beOtdZmJiAnt7e9jb20OhULAiU15e3qWuf8ad+UlRwt0RbYoK8CDT3dTUFKamphg9ejQ7A/3zzz9RUlLCPhxYWVmp2eOcOnUKy5Ytw+HDh3klKsCDQDelUtnBH8/GxgYlJSUcVaWf0BlLN6W5uRlvvfUWLl26hClTpuDcuXMwMjJCaGgoZs+ejcmTJ/O+d6MrXf9MlHBPdGfWtqg8CcZeRiaT4e7du0hISIBSqcSoUaPw97//HQcPHsTcuXN1WtPTUFVVBVtbW2RmZmLixIns9fXr1yM9PR0XL17ksDr9gs5Yuilvv/02qqqqkJubC0tLS7S2tiI9PR2xsbFYvHgxWltbERoaCqFQCH9/f15YpD+MsbExhg4diqFDh6p1/WdnZz+265+x/HdxcelxNusqlQr5+flQKBSciAoA9OnTB3Z2drCzs2OzZfbv34+DBw9i4MCBSE9Ph7m5OaZMmcKrhxtLS0sYGRlBKpWqXZdKpRpPD+3u0BlLN6WwsBAjR47s1KVXqVTi119/Ze3+GxsbMXPmTAiFQkybNo33mTKP6/pXKBS4cuVKj7P8B9RFxdPTkzeDdnp6OubOnYuYmBjY2dkhLi4O8fHxCA0Nxffff891eWr4+vrCx8cHO3fuBPDgZ2pvb49Vq1bRzftngApLD0epVKplyty5cwczZsxgM2X4bh/PBGFVV1dDIpGgra0N5ubmsLe3h4WFRY/pVeGrqGRkZGDOnDn46quvsHjxYnZmqVKp2NOAfOL48eNYsGAB9u7dCx8fH8TExOCnn35CSUlJj3tQ6QpUWCgsKpUK2dnZrElmZWUlmykTHBwMU1NTrkt8JLdu3UJZWRmcnJwgl8t52fWvLfgqKhcvXoRQKMRnn32GlStX6s2JvG+++YZtkHR3d8eOHTvg6+vLdVl6BRUWSqcwG8BMpsz169cxbdo0CAQChISE8MrN+MaNG6ioqICnpyfMzMwAqB+HlUqlaG5u5rTrX1vwVVSys7MRFhaGv/zlL1i7di1vvisU3UCFhfJEmCAsZiZz5coV+Pv7s5kyFhYWnAwcTAxBVVUVPD09HxslzGXXv7ZgbP/5Jir5+fkICQnBxx9/jPXr11NR6YFQYaE8E4QQlJWVsSKTl5cHPz8/CIVCzJo1S2eZMkziZXV19TNHCeuy619bKJVK5Ofno62tjVe2/8XFxQgODsbq1auxadMmKio9FCoslOeGEIKKigo2U+a///0vJk6cCIFAgLCwMK1lyrSPEvby8urSKTZtdv1rC76KSklJCYKDg7FkyRJs2bKFikoPpscLS0tLC3x9fZGfn4/c3Fy4u7uzrxUUFGDlypW4dOkSrKysEBkZifXr13NXLI8hhOD27dsQi8UQi8XIzMyEl5cXmykzYsQIjQw02owSbt/1z+TG29jYaN06/lloLyqenp68OZBQVlaG4OBgREREYOvWrT3mNB6lc3q8sKxZswZlZWVITExUE5aGhga88MILCAgIwCeffILCwkK8++67iImJwdKlS7ktmucQQiCRSNhMmV9++QXjx49nReZ5M2V0GSXMdP1LpVLWOp6xlhk4cCAnIsNXUblx4wZmzJiB1157DV9//TUVFUrPFpbExER88MEHEIlEcHV1VROW3bt3Y+PGjZBIJOwAtmHDBsTFxVHfoGeAEII7d+6wmTKpqalwdnZmReZpM2Xa2tqQl5cHQgjc3d11uvzTvutfJpOxXf82NjYazSd5Ug18FJVbt24hKCgIM2fOxK5du6ioUAD0YGGRSqXw8vJCXFwcLC0t4ejoqCYs77zzDhoaGhAXF8f+P6mpqXj11VdRU1ODwYMHc1O4HkMIQW1tLZspc/78eTg4OLCZMq6urp0OTO2jhN3d3TlNfVQqlaipqem063/w4MFaGViVSiXy8vKgVCp5JSpVVVUICgrC1KlTsXfv3h6Vxkl5PPz4huoYQggWLlyIZcuWwdvbGxUVFR3ukUgkHRxxmc5biURCheU5MDAwgLm5ORYuXIiFCxeivr6ezZR59dVXMXToUDZTxsPDA4aGhqisrER6ejpefPFFjBs3jvPBy8jICFZWVrCyslLr+i8qKgIhhA3B0lTXPyMqKpWKV6IikUgQEhKCSZMmUVGhdIAf31INsWHDBmzduvWx91y9ehXJycm4d+8ePvnkEx1VRukMMzMzREREICIiAo2NjWymTEhICMzNzREYGIiEhARMnDgRb775Ju+WWQwNDWFhYQELCwu4uLigvr4eUqkUJSUlaG1tZUWGMTd8VtqLioeHB29Epbq6GqGhofDw8MC+ffuoqFA60K2Wwhib7scxcuRIvP766zh9+rTa2rhSqYSRkREiIiJw8OBBuhTGIffv38fBgwexbt06WFhYQKlUYtasWRAKhZg4cSJvBthHwXT9S6VSVFdXo7m5WS1p8Wn2h/gqKnfv3kVISAicnJxw7Ngx3hx1pvCLbiUsT8utW7fQ0NDA/plZK46NjYWvry/s7OzYzXupVMr+8kRFRUEsFtPNey1z5coVBAQE4M0338Rnn32GlJQUiEQinDp1CkZGRpg1axZmz54NPz8/3g9shBA0NTWx1jLtkxYf1fXPV1Gpra3FrFmzYGtrC5FIpLeOBRTt0yOF5WEqKio6bN7X19fD2dkZgYGB+Pjjj1FUVIR3330XX3/9NT1urGWOHTuGq1evYvPmzWqzytbWVqSlpSE2Nhbx8fFoa2tDaGgoBAIBbzNlHuZJXf98FZWGhgaEhYXB3NwccXFxeuNQoAt+/PFHvP/++6iqqlL7DgqFQgwcOBCHDh3isDpuoMKCzoUFUG+QtLS0RGRkJD7++GPuCqWwtLW1ISMjA7GxsYiLi0NjYyNCQkIgEAj0IlMG6Nj1P3DgQLS2tqJXr17w9vbmjag0NjZCKBSib9++OHPmjF78bHWJXC7H0KFD8f3337PJmNXV1bC1tUVycjKmTp3KcYW6hwoLRe9RKpXIyspiM2VqamrYTJnAwEDeZ8oADwannJwcKBQKKJVKDBgwgJ3JDBgwgLO6mpqaMGfOHBgYGODs2bOc1sJnVqxYgYqKCiQkJAAAvvrqK+zatQvl5eW8cGzQNVRYKN0KlUqFy5cvsyaZVVVVCAwMZDNlHueAzBVKpRK5ubkAAA8PD6hUKshkMlRXV+Pu3bvo27cvKzK67PqXy+V4/fXXIZfLkZSUxOs8Hq7Jzc3FSy+9hJs3b8LW1hbjx4/H3LlzsWnTJq5L4wQqLJRuC5NVwojMjRs31DJldNU1/zgeFpWHj+4ymfFSqRR37tzRWdd/S0sL5s2bh5qaGiQnJ3Oe9FhRUYEtW7YgJSUFEokEw4YNw9tvv42NGzeqHSLg0t/Py8sL4eHhCAwMhI+PDyoqKjB8+HCdfDbfoMJC6REQQlBcXMyKzNWrVzF16lQIhUKEhIRwkinzJFHp7H6m67+6uhqGhoasyAwaNEhjfT4KhQLz589HZWUlfv75Z5ibm2vkfbtCUlISjh8/jnnz5mH06NEoKirCkiVLMH/+fERHRwPg3t9v9+7diImJwfTp01FWVoZz585p/TP5ChUWSo+DEILff/+dtfvPz8/H5MmTWbt/a2trrYtMW1sbcnNzYWBg8FSi8jDtu/6rq6s11vXf2tqKRYsWoby8HCkpKbC0tHyu99EF27Ztw+7du3H9+nUA3Pv71dfXY9iwYWhra8OPP/6IN954Q+ufyVf41cpMeSwVFRV477334OjoiL59+2LUqFH461//CoVCoXZfQUEBJk+eDBMTEwwfPhxffvklRxXzEwMDAzg7OyMqKgqXLl1CaWkpgoODcezYMTg5OSE4OBjffvstKisroY3nLkZUDA0Nn0tUgP91/Y8ZMwZTpkyBm5sbevXqhZKSEqSnp6OwsBBSqRRKpfKZ6lq6dClKSkpw/vx5XosK8GAgbz+bysrKwpQpU9SWxoKCglBaWora2lqt12NmZoY5c+ZgwIABEAqFWv88PkOFRY8oKSmBSqXC3r17UVxcjK+//hp79uxBVFQUe09DQwMCAwMxYsQIZGdnY9u2bdi8eTO+++47DivnLwYGBhg5ciTWrVuHzMxMXL9+Ha+99hpOnTqFF198EQEBAdi+fTtu3rypEZFpLyqaMtQ0MDDA4MGD4ezsDD8/Pzanpry8HGlpacjPz8eff/6J1tbWR76HUqnEypUrkZubi59//pn1xeMr5eXl2LlzJ/7v//6PvSaRSDrU3d7fTxdUVlYiIiJCL3qqtAldCtNz+LYc0F0ghODPP//EyZMnIRaL2UwZoVAIgUCAUaNGPfNymTZE5XEwXf+MtQzT9W9jYwMrKyv2+6FSqbB69Wqkp6cjNTUV9vb2Wq2rPU/r7+fi4sL+ubKyEq+88gr8/f3xww8/sNcDAwPh6OiIvXv3steuXLkCV1dXXLlyBWPGjNH8X+D/U1tbi7S0NISHh+PKlStwdnbW2mfpA1RY9JxPP/0USUlJuHz5MgBq968NmEwZRmRSUlLg4uLCioyLi8sTRUbXotIZTNe/VCpFfX09YmJiMHnyZFRWViIlJQWpqakdHL21zdP6+zEiWFVVBX9/f0yYMAEHDhxQ20vi8rvv4OCA2tpabNq0CR999JHWPkdf4EdrL+W5YJYDmFMxALX71wYGBgawsrLC0qVLsWTJEtTW1iI+Ph5isRhffvklRo4cydr9d5YpwwdRAYB+/frBwcEBDg4OqK+vx6RJk3DixAn8/vvvcHd3h0gkwpw5c3QqLkwEwdNQWVmJqVOnwsvLC/v37+/wc544cSI2btyI1tZW1kPu/PnzcHZ21vr3vrPojZ4M3WPhARs2bICBgcFj/3t4GauyshIzZszA3LlzsWTJEo4q73kwmTKLFi3C6dOnIZVKERUVhdLSUrz66qvw8PDApk2bkJOTA5VKhZqaGrz//vts8iVfLOYHDhyIlpYW1NfXIyMjA8uXL8fPP/8MZ2dnxMbGcl1eByorK+Hv7w97e3tER0dDJpNBIpGo7Z289dZb6N27N9577z0UFxfj+PHj2L59Oz744AMOK++Z0KUwHtBdlgN6Oo2NjUhISIBIJEJCQgLMzc3Rq1cvmJmZITExkTdd/4QQfP755/jhhx+QkpICV1dX9rXa2loYGxvzzrrlwIEDWLRoUaevtR/CqL8fP6DCome0Xw7497//3eEJmNr98wOJRAJ/f380NTWhubkZJiYmapkyXM1cCCGIjo7Gjh07kJKSAjc3N07qoHRv6FKYHkGXA/SDhoYGzJ49G46OjigrK8Mff/yB3bt3Qy6XY968eXBycsKaNWuQlpb22CPAmoYQgh07dmD79u04d+4cFRWK1qAzFj2CLgfoB01NTfjiiy/w6aefdsgtaW1tRWpqKpspo1QqERoaCqFQCH9/f62FZxFCsHv3bvzjH/9AUlISJkyYoJXPoVAAKiwUCme0tbXh119/ZTNl7t+/j5CQEISFhSEgIEBjYVqEEPzrX//Cpk2bcPbsWfj5+WnkfSmUR0GFhULhAUqlEpmZmWymTF1dHYKCgiAUChEYGIh+/fo91/sSQnDo0CGsW7cOp0+fhr+/v2YLp1A6gQoLhcIzVCoVLl26xDoxSyQSTJ8+HUKhEDNmzHjq02WEEBw9ehRr165FXFwcAgICtFw5hfIAKiwUCo9RqVTIy8tjRaaiogIBAQEQCASYOXPmYzNZYmNjsXz5cpw4cQIzZ87UceWUngwVFgpFTyCEoKioiBWZ0tJStUwZc3NzVmTi4+OxePFiHDlyBAKBgOPKKT0NKiwUih5CCEFpaSmbKVNQUIApU6ZAIBCgf//+WL16NQ4ePIjw8HCuS6X0QKiwUHTKrl27sG3bNkgkEri5uWHnzp3w8fHhuiy9hhCC69evQyQS4fjx48jJycGePXvULOUpFF1ChYWiM44fP4533nkHe/bsga+vL2JiYnDixAmUlpbC2tqa6/K6BYQQFBQUYPz48TqPWqZQGKiwUHSGr68vXnrpJXzzzTcAHmxMDx8+HJGRkdiwYQPH1VEoFE1BLV0oOkGhUCA7O1vtyKuhoSECAgKQlZXFYWUUCkXTUGGh6IQ7d+5AqVR2Gh2rq9hYCoWiG6iwUCgUCkWjUGGh6ARLS0sYGRlBKpWqXZdKpRgyZAhHVVEoFG1AhYWiE3r37g0vLy9cuHCBvaZSqXDhwgVMnDiRw8ooFIqmoZn3FJ3xwQcfYMGCBfD29oaPjw9iYmLQ1NT0yCgACoWin9AZC0VnvPHGG4iOjsZf/vIXuLu7Iy8vD0lJSR029Cn6SUtLC9zd3WFgYIC8vDy11woKCjB58mSYmJhg+PDh+PLLL7kpkqITqLBQdMqqVatw8+ZNtLS04OLFi/D19eW6JIqGWL9+PYYNG9bhekNDAwIDAzFixAhkZ2dj27Zt2Lx5M7777jsOqqToAios3QCZTIYhQ4bg888/Z69lZmaid+/eansaFIq2SExMRHJyMqKjozu8dvjwYSgUCuzbtw+urq548803sXr1anz11VccVErRBVRYugFWVlbYt28fNm/ejMuXL+PevXuYP38+Vq1ahWnTpnFdHqWbI5VKsWTJEhw6dKjTQLKsrCxMmTJFLXY5KCgIpaWlqK2t1WWpFB1BhaWbMHPmTCxZsgQRERFYtmwZ+vfvjy+++ILrsijdHEIIFi5ciGXLlsHb27vTeyQSSaeNscxrlO4HFZZuRHR0NNra2nDixAkcPnwYffr04bokip6yYcMGGBgYPPa/kpIS7Ny5E/fu3cMnn3zCdckUHkGPG3cjrl27hqqqKqhUKlRUVGDcuHFcl0TRUz788EMsXLjwsfeMHDkSKSkpyMrK6vAQ4+3tjYiICBw8eBBDhgzptDEWAG2O7aZQd+NugkKhgI+PD9zd3eHs7IyYmBgUFhZSO/rn5IsvvoBYLEZJSQn69u2Ll19+GVu3boWzszN7T3NzMz788EMcO3YMLS0tCAoKwrffftujjk/funULDQ0N7J+rqqoQFBSE2NhY+Pr6ws7ODrt378bGjRshlUphbGwMAIiKimJ/vpRuCKF0Cz766CPi4OBA6uvriVKpJH5+fiQkJITrsvSWoKAgsn//flJUVETy8vLIzJkzib29PWlsbGTvWbZsGRk+fDi5cOECuXz5MpkwYQJ5+eWXOayae27cuEEAkNzcXPZaXV0dsbGxIfPnzydFRUXk2LFjpF+/fmTv3r3cFUrRKlRYugGpqamkV69e5Ndff2Wv3bhxg5iampJvv/2Ww8q6D9XV1QQASU9PJ4Q8GCyNjY3JiRMn2HuuXr1KAJCsrCyuyuSczoSFEELy8/OJn58f6dOnD7G1tSX//Oc/uSmQohPoUhiF8hSUl5fDyckJhYWFGDt2LFJSUjBt2jTU1tZi0KBB7H0jRozA2rVr8f7773NXLIXCMfRUGIXyBFQqFdauXYtJkyZh7NixAB4ck+3du7eaqAA0X4ZCAeipMArliaxcuRJFRUXIyMjguhQKRS+gMxYK5TGsWrUKZ86cQWpqKuzs7NjrQ4YMgUKhQF1dndr9NF+GQqHCQqF0CiEEq1atwsmTJ5GSkgJHR0e11728vGBsbKzmxVZaWopbt27RfBlKj4du3lMonbBixQocOXIE8fHxar0rZmZm6Nu3LwBg+fLlSEhIwIEDB2BqaorIyEgADwxAKZSeDBUWCqUTDAwMOr2+f/9+tiOdaZA8evSoWoMkXQqj9HSosFAoFApFo9A9FgqFQqFoFCosFAqFQtEoVFgoFAqFolGosFAoFApFo1BhoVAoFIpGocJCoVAoFI1ChYVCoVAoGoUKC4VCoVA0ChUWCoVCoWgUKiwUCoVC0ShUWCgUCoWiUf4fRHdmy9WmkBcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "prody.showProtein(prot)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "With the latest version of ProDy, `py3Dmol` will be used if it has been imported." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "prody.checkUpdates()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import py3Dmol\n", "prody.showProtein(prot)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Configuration\n", "\n", "A number of keyword arguments are supported." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The background color is configurable." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,backgroundColor=\"black\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The height and width of the viewer can be changed." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": true }, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot, height=400, width=200,backgroundcolor='gray')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The viewer can be zoomed and centered on a given [AtomSelectionSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomSelectionSpec)." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot, zoomto={'resn':'ATP'})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Styles\n", "\n", "Custom styles can be applied by providing [AtomSelectionSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomSelectionSpec) and [AtomStyleSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomStyleSpec) objects.\n", "\n", "A single style may be applied." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,style={'cartoon':{'colorscheme':'chain'},'stick':{'colorscheme':'chain'}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For more complex styling, generate and manipulate the py3Dmol view object directly." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "view = prody.view3D(prot)\n", "view.setStyle({'cartoon':{'colorscheme':'chain'}})\n", "view.setStyle({'hetflag':True},{'stick':{'colorscheme':'magentaCarbon'}})\n", "view.setStyle({'bonds': 0},{'sphere':{'radius':0.5}})\n", "view.addLabel('Chain A',{'inFront':True,'showBackground':False,'fontColor':'blue'},{'chain':'A'})\n", "view.addLabel('Chain P',{'inFront':True,'showBackground':False,'fontColor':'green'},{'chain':'P'})\n", "\n", "view.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# GNM/ANM\n", "\n", "Squared fluctuations (or any other alpha-carbon length numerical vector) can be visualized in a RWB colorscheme. Red values are larger." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "gnm, sel = prody.calcGNM(prot,cutoff=7.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Individual modes may be visualized." ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,data=prody.calcSqFlucts(gnm[0]))" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "Styles can still be applied." ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,data=gnm[1].getArray(),style={'cartoon':{},'stick':{}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Motions\n", "\n", "It is also possible to provide an array of vectors (one per a calpha) that can be used to animate motions. The animate kwarg takes the 3Dmol.js animate options." ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "anm, sel = prody.calcANM(prot)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "scrolled": false }, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,mode=anm[0].getArray(),anim=True,\n", " data=prody.calcSqFlucts(anm[0]),\n", " style={'cartoon':{}})" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'2.4.0'" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.__version__" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10.6" } }, "nbformat": 4, "nbformat_minor": 1 }