{ "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": "iVBORw0KGgoAAAANSUhEUgAAAcUAAAE1CAYAAACWU/udAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl4JGd19n1Xb1JLrdE+m0YzGi0zkmZfPAsYM7Yx2A7YmBjbkAQHAgkJJIa83xWcD5sAJnjI64vEwbzgNwHb8BHAGALEYLOMbcDLePYZj2XPqDftrZZ636uru74/pKemqruqu6p60fb8uHQNbnUt3equu855zrkPw/M8KBQKhUKhAIaFPgEKhUKhUBYLVBQpFAqFQpmHiiKFQqFQKPNQUaRQKBQKZR4qihQKhUKhzENFkUKhUCiUeagoUigUCoUyDxVFCoVCoVDmoaJIoVAoFMo8Jo3Pp/Y3FAqFQlmKMGqeRCNFCoVCoVDmoaJIoVAoFMo8VBQpFAqFQpmHiiKFQqFQKPNQUaRQKBQKZR4qihQKhUKhzENFkUKhUCiUeagoUigUCoUyDxVFCoVCoVDmoaJIoVAoFMo8VBQpFAqFQpmHiiKFQqFQKPNQUaRQKBQKZR4qihQKhUKhzENFkUKhUCiUeagoUigUCoUyDxVFCoVCoVDmMS30CVAo1YTneWQyGTAMA4PBAIZRNYybQqGsEKgoUlYERAw5jkMqlQLP8wAAg8EAg8EAk8kEo9Eo/DcVSwplZcKQi4NKND2ZQlloxGLI8zwYhkE6nYb4c8/zPHK/BwaDAUajUfihYkmhLHlUfXmpKFKWJXJiSASNZVnhMaVtyb/i70cymUQ6nUZLSwsVSwpl6aHqS0rTp5RlBc/z4DgOmUxGED6DQVs9GRG4XKFLJpMIBAJYtWoV0um05Hc0sqRQlgdUFCnLAiKGHMcBgKIYJhIJjIyMwGg0wmazob6+HjU1NarEizwnd78kmuQ4joolhbLEoaJIWdLIiaGc4MTjcTidTkQiEWzYsAHZbBaBQADj4+NIpVIwGo2or6+X/FgsFsm+GIbJW3skj4v/FZ8bMCeWLMtKfk/FkkJZnNA1RcqShIih0+nExo0bFcUwFovB6XQiHo9j8+bNaG9vB8dxyGazkudzHId4PI5oNIpYLIZYLAaWZWEymQSRzGaziMfj6O/vL/ncxT8Mw4BlWQQCAaxbt46KJYVSGeiaImX5kc1mhQIaAJicnERXV1fe86LRKBwOB5LJJHp6etDa2lpQWEwmE1atWoVVq1ZJHuc4DrFYDNFoFKFQCOFwGCdOnIDJZBLSr+THbDareg1yAk4i17Vr1yKdTstGliaTSYgwqVhSKJWBiiJlSZArhkqRYSQSgcPhQDqdRnd3N1paWkoSDpPJhMbGRjQ2NsJms8Hj8WDr1q1Ip9NCROn1ehGNRsFxHMxmc55Ymkzqv2Zy66AkopSrms2NKo1Go+J7Q6FQikNFkbKoyWazQjUpoCyGoVAIDocD2WwWPT09aG5uruh5mc1mNDU1oampSfI4y7KCWHo8HkSjUWQyGVgsFtTX10sE02g0qjqW0mvmeV64WSgklkQwqVhSKMWhokhZlKgVQ47jcPr0aTAMg56eHjQ2NlbsnJQKbcRYLBZYLBaJKJMoj4jlxMQE4vE4MpkMampqhGiSZVlkMpmKiKXb7UZXVxcVSwqlCFQUKYsGkiZMp9PIZrMAlC/8fr8fDocDLMti165dFRVDghpRVNqupqYGNTU1aGlpER7neR6pVAqxWAyBQACJRAJnz55FNptFbW2tEFHabDbU1dWp7reUe8/8fj+6u7sFscyFiiWFMgcVRcqCo1YMeZ6Hz+eD0+lETU0N+vv7cfHiRTQ0NCzEaZcMwzCora1FbW0trFYrkskktm/fDp7nkUwmhcjS7/cjHo8jm83CarVKxNJqtZYkloA0ssyFiCQp8qFiSVnuUFGkLBjkYsxxHIaGhtDf368ohrOzs3A6nbBarRgcHITNZgOgP3rTQ7WOxTAMrFYrrFYr2trahMd5nkcikRDEcnZ2FvF4HAAkYllfX192sWRZVvJ8UtRDxZKy3KCiSKk6YjEk/YLBYFDWKcbr9cLlcsFms2HHjh2oq6uTPEePUOm9cFdTgJWOX1dXh7q6OrS3twuPZ7NZiVh6vV4kEgkAQF1dHerr64U+TKvVqvr1qxXLeDyOUCiE9evXU7GkLHmoKFKqhpwYKkWG09PTcLlcaGxsxK5du2C1WmX3aTAYhJTrSsVgMAgRohhiNhCLxZDJZOBwOJBIJGAwGASxJD+1tbW6xTKTySAej4NhmIKRpdx4LiqWlMUGFUVKxVGaWCEnhlNTU3C73WhubsaePXtQW1tbcN/LMX1aLgwGA2w2G2w2G0ZHR7Fjxw4AV0QsFoshFAphcnISyWRSIq7kR40vbKG/Kfl9NptFKpWSPM4wjKx7j1YDdwqlnFBRpFQMOTGUu+Bls1mwLItXXnkFra2t2LdvH2pqalQdQ2+kWGh0lBJLTRSVMBqNaGhoyCtQymQyQgpWiy9ssfdS6QaIHJMYMoifr2RKQKFUGiqKlLKjRQwnJiYwOjqKbDaL/fv3w2KxaDrWchEqAAv+OoxGo6LVHfGF9fl8GB0dlfjCkuewLKv671fMRF0sloFAAACEOZbi9UoqlpRyQ0WRUjbUimEmk8H4+DjGx8exZs0aHDhwAKdOnVLtHSqmmmuKlRTgxXxhV/KFTafTiMfj8Hg8iMfjeP3115FOp0v2hRX/C8zNsSTRZm5kyfO8pLiH+sJSSoWKIqVkyMSKUCiEmZkZbN68WVEMx8bGMD4+jnXr1uHgwYOCLygp0lDr5kLQKlQsy8LlciGdTgvrbXV1dar8SZdTVFoOzGYzGhsbhQhx8+bNAObEkkwbmZ6eRiwWK8kXlud5RU/X3MhSnMolaVgqlhQtUFGk6CZ3lmE2m0UsFsu76HAch9HRUUxOTqKjowOHDh3KuxiSiE+rKKqNFFmWhdvtxuzsLDZs2IBVq1YhHo9jcnJSqM4sxUVmJZO7pmg2m9Hc3JznPyu2upuamhLed4vFkieW4s9BNptV/DuoTcOKz5HOsqQUgooiRTNKg32NRqNEoNLpNEZHR+HxeNDR0YHDhw8rip7eNGix6I1Ehj6fD5s2bUJvb6/gntPa2ip5TWIXGZ/Ph3g8Dp7nhfaF2travAssRX3Rkl5f2EQiIaROtfjCiv8VHw+Yu1FLp9MA5iarsCyL9vZ2KpYUKooU9RSbck+EjWVZjIyMwOv1orOzE4cOHSp6MSu3KOaKYV9fnxBtyNmZKbnIiBvjw+Ew4vE4Tp48KfT6kQjHZrNJKjJXEqXcJKjxhY1EIvD5fPB6vWXxhRX/C8ytWZJ2EbFYEmhkubKgokgpSjExJHAch3A4jJMnT2LTpk04fPiw6ouVwWDQtV6XK6aFxFAP4t695uZmhMNh7N27V+j1i0ajee0L4lSgzWbTVUC0lKhE5Cz2hQ0EAmhpaUFLS4tsRJ9IJCS+sOT9V2t1R9L2hdYsqViuHKgoUhRRO9g3mUzC7XbD5/PBaDRqEkMCKbTRCokUyy2GSsciKPX6cRyHWCyGaDSKmZkZuFyuvCITLcU9S4FKp5Nzi2fU+MLOzMzI+sISE3Xx+ZayZknFcvmxPL6VlLKidpZhIpGAy+VCKBRCV1cXent7cfr0aV1iVEoT/tjYGBKJBLq6uioihgQ11acmkwmNjY15o6xYlhUqMicmJhCLxSSpQIvFItjfLbXinkqLopr3pBRf2EQigfr6ek2vo5hYptNpsCwrfGampqawYcMGKpZLACqKFAEihm+++Sa6u7thMplkv7DxeBxOpxORSASbN2/GwMCA8OWXW69Tg1ZRJJHh1NQU1q5di127di1qMbFYLEIKkCBOBQaDQcTjcZw+fRoAZFOBi/XiSXoFK7l/va9djS+s3+9HLBaDx+Mpiy8sAMkaOsdx8Hq96OjokIil+BypWC4eqCiucORmGYbDYdkLUSwWg9PpRCwWQ3d3N7Zt2yZ5TilfYrWimEqlhFRtV1cXNm7ciPr6+qoIYrn7FMWpQJvNhng8jp07d0qim0gkAo/HIzHyXmzFPYshUtSK2Bc2Go2isbERbW1tqn1htbz3ZM1S7jWIv39yYinusVRa96SUFyqKK5RCg32NRqMk4otGo3A4HEgmk+jp6UFra2vZv5jFCm1yxXDLli1gGAYul2tZONqIEV+AV69eLTy+WIt7qrmmWAnE/bHl9oUl2ypVXxcyUSdr5bmvX8kXlopleaCiuMJQM+XeYDAgk8kgEonA4XAgnU6ju7sbLS0tFfviKRXapFIpuFwu+P1+iRgW267YsZYipRT3ZDIZTX1+WliKkaLW/RfyhRVXwub6wurNYhSbOJLJZCTvezQahdlsRkNDg0QwqVhqh4riCkHtLENg7s52aGgIBoMBPT09ec4klSA3fZorhlu3bpU9V72tHHpYrBcXueIecuNDinvS6TTOnTsnce4hglmqc081IsWFFkUllAqriC9sNBqF3+9HKBTCiRMnSvKFBZTF0u/3o66uDjU1NQUjSyqWxaGiuMzRIobBYBAOhwPRaBTd3d3o7Oys2nkSUVQrhgTqRyoPwzCS4h6Px4N9+/ZJinui0ShmZ2fzWhe0FvdUI1JcapEo8YVtbGwUxKqvr6/svrCETCYjrD+KEUeWuRCRFE8doWJJRXHZonawLzB3l+lwOGAymdDX14fJyUnU1dVV9XzJGCmXy4XNmzcXFUOCwWDI6xOjKKPGuUdrcc9KSJ+WQiaTEfYv5wubG9WLfWGJ1Z2SL6z4GHIiqiYNy7Ks5HcrXSypKC4ziBhOTU0JpeVKVW9EDC0WC/r7+4W1qunpad2tFWTfar9AJDL0eDxoa2vDvn37NH35lkukuNAXnFKKe0ifX6VYDqJbaC03N6oniH1ho9Foni+sOLLkOE7TerFWsWQYBslkUoh6l7NYUlFcJuRGhsFgEAzDwGaz5T1vdnYWTqcTVqsVg4ODec8pZUah2mkX4jTp5s2bhchU6xesmvMUK81iFHc1xT2k0GRkZAQWi0VSBasU2WhhOUSKet4DNb6wpM/S5/MhEokIN8J614sLiaXD4UBXV5dwwzw2NoZHHnkEjz76qObXtpihorjEURrsm9tWwfM8ZmZm4HQ6YbPZsGPHDsUUae62WigmiqlUCk6nE4FAQJImnZiY0JUG1RopRqNR2O12RKNRoT9wKTTILzbEBSbxeBxtbW1oamqSpAHlnHtKuVhXisUqikqIfWHJpJezZ89i27ZtQuuIeL2Y53ndvrDiY5IIlbyWYDAouAMtJ6goLlGISbe4NFv8ISfCxvM8pqen4XK5sGrVKuzatQtWq7XgvklLhh5yx0cRcsWwv78/rw1ET8Sntvo0FosJvZbd3d2w2WxIJpOIRqN5a2jiiwdZQ6MoQ97/QmnAQsU94jTgQtyYVCMSrbTPLVlTtFgsBX1ho9GoxOpOS3EVx3GS1xEKhfJaVJYDVBSXGHITK+S+0AaDAX6/HyMjI2hubsaePXtQW1ur6hhGo1F38UquoCaTSbhcLkUxFG+n1xC80HbEki4WiwnGA6SwgXhl5q6hkYuHz+fDyMgI0um0UB2YTqcRDofLkhZcThQSFT3FPeLCkmw2W/F1xUruu1L9oWKKmZoX84WNRqPweDxIJpMAILG6s9lsqK2tzWuNCYfDaGpqqujrWgioKC4R1I5vymazmJqaEtYM9+3bh5qaGk3HMhqNwpdDK0TcxGLY3d2tKIbi7fSsqSmJYiKREPxZe3p60NbWJhy/0HGUmrSJobfX65WkBXPvtOvq6lZcClavYKkt7kkmkzh58qSkIZ6830thLFc1RLHU91+M2BdWfLMSj8cxNDQEnucxNDSE6enpkiLFS5cu4c477xT+2+l04otf/CI+9KEP4c4774Tb7UZXVxeefPLJqvRKE6goLnK0iOHExARGR0fR3t6O3t5eJJNJzYIIlFa8QhbkiT9qMTEk6HGmAfLFNJlMwul0IhQKoaenB4ODg2URKZIWtFgsGBgYACCflorH44vWo7RSlLuPMLe4JxAI4KqrrkI6nZZMu4hGo+A4TijuEadhF1MUr6bwrFTKWaQl9oUVc+LECXR2dmJkZAQvvvgizp07h9nZWfz0pz/FwMAA9uzZg09+8pOqj7N161acO3cOwNyNQ0dHB2677TYcPXoU119/Pe69914cPXoUR48exVe+8pWyvb5iUFFcpPA8j0gkUtTbMJPJYGJiAmNjY1i9ejWuuuoqWCwWzMzMIBaL6Tq2nkIbEhn6/X5s2rQJO3fu1HShLCV9SqrxyJpld3e3MLlDaRu9iC8+SmkpcaTj9/sxNjYmWH+JhVJrg3Y5zrlS+6+G4JvNZjQ1NUlSduK2hXIX95QLcZ9iJahG1TK58WloaMD27dvx8MMP4+jRo9i9ezduvvlmvPnmmxgbG9O9/2PHjqGnpwebNm3Cz372M7zwwgsAgLvvvhtHjhyhoriSEQ/2HR4eRkdHh2zqIJPJYGxsDBMTE1i7di0OHDggSSWZTKaqjHEiYhgMBrF582bwPI/m5uaqtVZkMhkEAgGcPn264JplNVFqYyCVmdFoVNKgXVtbC6vVKpTZa60MXGgWspWkUNuCmuIem81W8fOvdPq00jZ4gLw5QCgUQlNTE+rr67Fv3z7s27dP9/5/8IMf4AMf+ACAuT7pdevWAQDWrl2L6elp/SeuAyqKiwS5wb4mkylPKDiOw+joKCYnJ9HR0YGDBw/KRhulpEDVRIriNKVYjEKhkO40qNZ5im63G9PT0zCbzTh48GBVxLCUYyi5mSSTSQSDQczOzsLlcgmVgeL1s1JSsNWI5Bb6RiSXYsU9pOp4amoK8Xgcp06d0j0aqhiVFkWtjft6j5F7nQmHw3mer3pgWRY///nP8eCDD+b9biHMAagoLjCFptyLxSmdTmN0dBQejwcdHR04fPhwwS+C0WgU1iG1UkigcsUwN02p1JKh5phq7tjT6TTcbje8Xi82bdqE3bt3w+FwLLqLslrIxdtoNGJ6ehrbt28HMPe5kBtVRIpNxP2V1UjBFqKSolvuKC63uCSbzeLMmTPYs2ePpLhHnPIutbinGn2Qlf4MyIkiiRRL5ZlnnsHevXuxZs0aAMCaNWswNTWFdevWYWpqSlKEVQ2oKC4Q2Wy24PgmYC4FmkwmMTw8DK/Xi87OThw6dEjVXaFecSLb5kaKxcSQoLfHsVihDcdxGBkZgcfjwcaNG3H48GEYDAahOXm5YTAYFFOw4hL6aDSaZ/tls9kk62dLed5htSzeCqW8yc2J2MBbS3FPpSPFalS3KkWK5RDF73//+0LqFABuueUWPPHEE7j33nvxxBNP4NZbby35GFqgolhF1MwyJLAsi9nZWaG/joiAWsrhSgOoF0O5bfUeUwxJF09NTaGzszPvfai29+lCC7BSsUkqlRLWK30+n7B+RgqB0um0UI1cbpGp5JrWQvuS6inuEQtlXV1dxYW9WqKYe4xypE9jsRh+85vfSKzi7r33Xtxxxx341re+hU2bNuHJJ58s6RhaoaJYBbSIodgTtLGxEW1tbbpGOJUiiqR5f2hoCKFQqGg1p5hyiWImk5GsnSpFyMvJ+1QvYtuv3PWzeDyOSCSCTCaDS5cu5Zl5l6Pfr5I3CtUQFK2iq7W4JxaLYWhoKK8ZvlyvSy6KKzdyxyARcynU19fD5/NJHmttbcWxY8dK2m8pUFGsILmzDAFlMSQRGani3Lp1K7xeLyKRiK5j6xWLZDIpmamoVgzFx9UjxuR8SVXt+Ph4QTEkVDtSXEprl6TfrLa2FlNTU9i1axcAqZm3OCVIxhSJIx01glFJ4VpKZuBKxT0nTpxAV1eXoqWg+D3XU9xTjUgxd91yoTMmlYSKYgUQi+HFixexbds2RTEkzivhcDgvPVlKtKf1i5VIJOByuYTIMBgMCgvfWihlvmEikcDx48exbt06HDp0SNXdL40UtSM3LZ6kBEkK1u/3C32uuS0MuVHOclhTrCQMw8g6xxBLQTLpYnR0VFLcI745KRTJVyt9KjdAYCndJKqFimIZkZtYEQqFZD84sVhM8OTs7u6WdV4pRRTVkiuGRJQdDoeu/Wkt8BE78WSzWcUWEyUWYp5ipS7UlXodas5XnBIkkxcAaQtDOBzG1NQUkskkjEajcMFmWVaYHl9ullKkKAdZOpFDyVKwUHGPWChJcc9CtGQkk0nVXspLDSqKZaDQlPvci1E0GoXT6UQikUB3d7fEkzMXk8mku62iGOIIVU+aVAm16dNsNovJyUmMjIxg9erVOHDggOBvqfV4WiPFTCYDv9+Purq6sq7tlMJiOAc5xCk+ceZAnIIl68+5KVhSBVvKBXupR4p6ipCKFfeQgcOkuCebzQoCabPZKmL+kCuKwWCwLD2KixEqiiWgNMtQDGk1IKOLWJZFT08PWlpain7ZKxEp5ophubxBCcVEiud5TE1Nwe12o62tTbCl04uWc89msxgbG8PY2BiampowMTEhRD3i9KDNZlMUZxKZLlYRU6Lc5ytOwU5NTWHPnj1gGEZIwcZiMYyNjSEejyObzeYNv1U7ImqpR4rlSm0WKu65dOkSLBYLEomExLkn9z0v5QYwVxTL1bi/GKGiqAM1Yih+7tmzZ8HzPLq7uyUf6GKUGikSQTYYDBIxLKdRdi5KosjzPDweD1wuF1pbW7F///6qzSkUR6Vr167FwYMHJSbWuYUnTqcTHMdJvDPJHfhCpGtLpVrep0opWJ7nZacuqJldWWkLs6UuuuTa09LSIoksc9Pek5OTSCaTuot7VsosRYCKoia0iGEwGBSG2vb09OgqWik1UjQajYjFYhgdHdUlhnoiolxRFA85bmpq0jXKSi/iqLS9vV3whyWpKIJS4Ym4929mZkYYn3Pp0iU0NDQsqSHEC2nzJi400TK7sr6+vuCaXDmotFn3QjXWK42FEhf3+Hw+1cU9ueuW5XKzWYxQUVRBsSn3YgKBABwOB4xGI3p7ezEyMiJbtaWGUioryQX8tddeQ29vr+bIkBxb6xeaFNrwPA+v1wun04nGxkZNQ45LRXzs5uZmXVGpUu/f6dOnsX79eiQSCcmFXFwEQf5dLKbeizWyLTa7kly4o9EoTpw4IZldKY7cS6HSY52q1S6h9hh6i3uIWT1Zu1yuA4YBKooF4XleKEsn45vkLnQ8z8Pv98PhcMBisWDr1q2CXVQ12yoA6ZphbW0ttm/fnjcXTQ3kvLV+oRmGQTwex6uvvoqGhgbs3r0bVqtV8/H1wPM8fD4f7HY7GhoaKiLE5A4894Igbmcga2k8z0vWdWw2W0UcZZYbZHZlS0sL6urqEAwG0d3dLZldOT09LaRgyexKPb1+C+2YUw7KIbzFins8Hg8mJibgcDhw//33o7m5Ga2trXjyySexfft29PX1aao+DgaD+OhHP4qLFy+CYRh8+9vfxtatWxd0uDCBiqIM4sG+drsd7e3tkmhB/LzZ2Vlhyv3g4GCeAFWyglRMIpGQNN0PDg7i4sWLuiNNrVEqEaTh4WGkUikcOHBAd4SsB7/fD7vdjtraWuzcubPoscstTOILOUG8rhMKhTAxMSGYeucW9lSapSrERLTUzK7MTQeKhVLJOD2bzVbUDabS6VmgclMyxGvEZLj2wMAAbrrpJnz5y18Gy7Kw2+346U9/inA4jKefflr1vu+55x7ceOONeOqpp8CyLOLxOL785S8v6HBhAhVFEXJT7i0WS56o8TyPmZkZOJ1O1NfXY/v27Xm5e0KlRTFXDIlRAFC61ZvabUl0ZrVaMTAwgOHh4aoJYigUEtZN5W5Kyo2WQhuldgaxqTeZq8hxHFiWhcvl0lyhWYzFmj5VQ7F1bSUjb3H7wuTkpKw3qc1mWxZritWIduUm4Vx33XW45ZZbNO8vFArh97//PR5//HEAczeUFotlwYcLE6goQjrYF5BasYlFTVw0smrVKlURSSkjnAhyF4Z4PA6n0ykrhuJjVzJSDAQCsNvtsFgs2LZtG2w2m8TSTita2h0ikQjsdrtwodu9e7euY2qlHNWncqmqdDqNc+fOwWazCRMwSHpQHPHYbDbNTfJLsYWEoPeCTy60crMrxcVTwWAQRqMRs7OzZZldmUs1RBGobCZAbjRVKS0ZLpcL7e3t+PCHP4zz589j3759ePjhhxd8uDBhRYtiITEkmEwmsCwrVDFqLRopNVIkERv5UKoRQ4JeH1LxceUIBoMYHh6G2WzGwMCAJDorpTiItJAUuojEYjHY7XawLIve3l40Nzfj5Zdf1nysxSYUDMPAaDSivb09Lz1IIp6ZmRm4XC6hSV5pVNRyopx/J7E3KXmP7XY7mpqaUFNTIztLMbd4SmuqtVo2cpVEaZaiXlHkOA5nzpzB1772NRw8eBD33HMPjh49KnnOQgwXJqxIUSw02Df3eaFQSLiD2bt3r+Z2ApPJhFQqpftciTixLCuIYU9PT0ExFB+7lPFRuduGQiHY7XYYDAb09/fnpazIdnojqUIVryRNHIvF0NvbK+mDqybV7lOUqxbM9SkVj4oSRztLpV2kENVIDZpMJsVZiqQKlqS5M5mMpH+VjIdSOsdqDACu9OdRaWyU3iKYDRs2YMOGDTh48CAA4Pbbb8fRo0cXfLgwYUWJohYxJH6cpLeqv79f1zHLsab4xhtvCP2OasSQUGqkSCK+cDgMu90OnufR29tbMScLuSgzlUoJ00N6enrQ3t4u+/r1RBR6tlkMzfuFfEpJVCk2mCYm7VNTU2WxXqsmC9m8bzab0dzcLJuCzR0PBUCogiWCWVNTIwyAXojzLxflHjC8du1adHZ24tKlS9i6dSuOHTuGwcFBDA4OLuhwYcKKEUXSjA0UFsPx8XGMjY1h9erVuOqqq5BIJDA2Nqb7uHpFkaRJSdP9xo0bNV/ASy20iUajGB8fRyaTQW9vb8X7ksSiSIpOfD4furu70d/fXzBNrKX0nfz99YrbQouiEgaDQTbiCYVCcDgc4DgO4+PjiMVi4HlemH5RDhuwSrHYqkOVxkOR2ZW5lcbpdBo2mw0sy5ZldqXc+S+EOUC9yLdoAAAgAElEQVQsFiupqO1rX/sa/uRP/gQsy6K7uxuPPfYYstnsgg4XJqwYUVQy6AYgzPCbmJjA2rVrBecTAJJqVD1oFSbxmmFPTw8YhsGqVat0XayMRqPEuUUtpNCD53ls27atar1CJKKZmJjA9PQ0urq60NfXV/SiVc3obbGJhhpMJhNqamokw6rF1mtiGzDiAysWy0qn/wqxVGzYxAVRYi5duoS6ujowDKM4u7KUNeGFmJBBXIZKed92796NU6dO5T2+kMOFCStGFIH89S6O4wQxXL9+vezYIpPJpHs+INlejajG43FhzUycJg0EArpFWasgE9PyZDKJlpYW2Gy2qgkiKSg5d+4curq6cPjwYdVfumrOVFwM6dNyoGS9xnGcsFbp8XgQjUaFdTSxWJarXaQY1ZiSUUlR4XkejY2NeWvCxCEmd01Ya/S+EAOGCUvxBlENK0oUCRzHYWRkBFNTU9iwYQMOHz6s+MEqdU2w2PZEDOPxuOwoqWr0GorPobe3Fy0tLZiYmKiK6QCZXDE+Pg6GYbBjxw7NQrxchKqSr0HtBcxkMsk6m4hbGbxeLxKJhCCsLMsiEAjoahcpxlKJFJWQS8+KLQSVZleGQiFJ9J5rnE7e52oU8uQOGK7GOuZCsqJEkeM4OJ1OeL3eomJIKKWaElAWpmJiWGz7Uo5NEDf+9/T0SM5Bb+qVUOwOX25yxaVLl3TdfWqJFMm68ejoaF6q0GazFf08VEqAK3nXXer5yrUyAFei+0AggNnZWbjdbokPbKmpQXLuS3meopZIVM3sSnFbjsViEb7jkUikYgVUuenTSCQiW3m+XFhRouhwOGC1WjWl5kolV1TViiGhlEhVSRSTyaRkpqJS438p7RxKFzPxCKm2tjbJ+q3eNKgaoSITM1wuF9asWYO9e/cK3rbEWSYajSKbzUpSWDabTZLCWi5RaTkgNxUWiwV9fX0ACreLiD1K1TbIL4dIsVShUpriwrIsxsfHEYlEMDY2JhRQ6Z1dqcRKGhsFrDBR7O/vr9raUy6xWAxOpxPxeBw9PT1obW1V9UEtJWLLFTZxe0N3dzcGBgaKVnTqgbSCiC824skVSiOk9EblpOlfDmLJ53A40NTUJAw1ZlkWPM/LpgpJCisSiWBqagrJZFJo5CYVhkutraFS0Vbuvgu1i5D3LhAIYHx8XOIDK15HE7+v1TBZqLQbTCVEl7zPdXV1sFgsQhEVScGS2ZXk8yueXSk2TldDbjFPOBymorhc0PvhFw/r1UosFkMikcDFixc1iSGhHOnTVCoFl8sFv99ftL2hHMcVC6p4coXNZivoBlRI3IodT05M/X4/hoeHUV9fr3pah9h4WlyAQhq5w+EwZmZmMDExIUyUF1/UV9oUDLWipVSdSd7XaDSKiYkJwaOUROvxeFy4gVmK72u1R1OJU7BKsytnZmaEVDeZXSk2Ts8939xIMRgMLtuxUcAKE0W9kBSmFncQUsmZSCRgNptx1VVX6RLVUtKnZK3h9OnT6OrqwtatWzU1/pfqmxoIBDA8PKx6ckW50qfhcBiXL1+G0WgUPFmLbVMM0sjd0NCAtWvXorGxUWhrUJqCoRT9LATVihS1otQgT6J1r9eL8fFxuFwuScEJ+VnIdhG1VFLMOY5TdbNXbHal0k1JfX09UqlUXqRYKQOPxcDi/0SVEb0fTi2iKBZDEhmeOHFCd6SpJ2JLp9Nwu93wer1gGAaHDh3SfOxSIsVMJoMLFy7AYrFomlyhVxTJdrFYDMPDw+A4Dn19fRX/4orbGnKnYMhdaMRRZUNDQ9lMp4tRyTXQSkRw4mjd6/Wiu7sbdXV1koKT6elpOBwOwTEmt11kOVdHiil1zVJu5Jn4piQWiyGZTOLMmTMIh8N47LHH0NDQIPxt9FixdXV1oaGhAUajESaTCadOnYLf718UsxSBFSaKelETrcmJYe6kDT13tVrEKZ1OY2RkBNPT09i0aRMOHz6M48eP67pA6BEoMrkiGo2iv79fcLyv5DGBuQuD0+lEOp2uuC+qmghTLvopNFuRXGQymUzF+s4Wa6SoZf9KBSepVEoy+YK0i8gV9iw3KtGSIb4pAYCZmRkcOHAAyWQStbW1+MEPfgCXy4U//dM/hdfrRV9fH370ox9pOsbzzz8vcQQ6evToopilCKwwUSw1UpRDLIakx09u0obeFKiabUnfpcfjQWdnZ1mqa7WIsbjpv6+vD5OTk6pSOrloLbQhJumzs7PYuHGj4ABUSfRWnyqV25P0VSQSQSqVwtmzZ4UKQq2VmkostUhRTLEMi7jnT3yRzR0+PDIyIruGttQrias1mgoAamtrcfXVV+Oll17CzTffjDvvvBPAXDVqqSyWWYrAChNFvcgJkxoxLLS9WgqJE8dxGB0dlZgQlCttpCZqE/c5kgiNYRh4PB7da4NqthPfBHR1dWH9+vVoamqqWiFGOS6kPM9j8vgkZl6bARtjkQqnMOmaRLOtGalICslwEqlwCqlICul4GjAABosB5hozzFYzzLVmWKwWZLks6tfUY8utW9B9U3cZXp3211FpUdRrcag0fJhElWNjY4jH4zh58mReG045CqaqIbiVtnmT+/vmrilqXaZgGAbvfOc7wTAM/uqv/gp/+Zd/uWhmKQJUFFVhNpsFURNHRT09PQXFkFDuBnyxV2tHRwcOHTpU9i9GoXPOnVyR2+dYytpgIUu9bDaL0dFRjI+PSyLi4eHhqtq8lYPjR4/jlQdfyXt8AhO69jf0X0N41/ffhb4b+vLShNVsyajE/su5PiheQ+M4DqlUCrt371ZMbZcyT7Eazi+VdrRRGhtVynr9iy++iI6ODni9Xtxwww15E4gWcpYisMJEsZT0aTQaxYULFzSJoXh7vZGiOKWYyWQwPj6O8fFxRa9WOco1Ikk8uWLz5s2KrR2lFszkwvM8Jicn4Xa7sXbtWhw6dEjyuqttCF7qsYKuIF7936+W6YyuMH5qHOx6VkgTNjQ0CCnCSr0/izVSVLtvg8GgmNpOp9MScwcyT9FqtUqqYJWa46vlS1rJY5R7bBQAdHR0AABWr16N2267DSdOnFg0sxSBFSaKeiAVhMlkEtu3b9ckhoRS/VN5nsfo6CjGxsYESzS1d4dkLqLWL474NXIcB7fbLRTwFJtcUYooii/epOHf4XCgtbVVaLyXO9eltDb0yoOvIMuVN7KtaazBW/7iLWjYMJcuZFkWkUhEmK1IeizLPYR4qUWKYopFcmazuagP7PT0NBKJRF5zvM1mq5ooVjIaLbcokkrshoYGxGIx/PrXv8bnPvc53HLLLYtiliKwwkRRy5c3Go3C4XAglUqhvb0dLMvqrmrUK4rEHzQWi4FlWYklmlqIu4yeLyfP83C5XJicnNRUwCMeUKz1XMl2Pp8Pw8PDaGhowN69exUb/nO3U0spRg6lCPDsG7N444dvSB4buHMAjd2N8Aa96O7vhsVmgdlmhrneDIvNApPVBD7DI+6N48fv/bHMSQF3/upO2DpsCNgDaOxqhMViQWtrK1pbW7Fq1SrMzs6ip6dHMpmBFJ+U4lVazerTcqMnvVnMB1bsT8qyLDiOg8PhEASzFB/YQudUKeREMRQK6RbF6elp3HbbbcK+P/jBD+LGG2/EVVddtShmKQIrTBTVIBZD0loRDAYxMaFvrQfQbtWWzWYxNTUFt9uN1atXo76+Hr29vbqPrXU9k5hmx2IxANC8ZkmEWCsMwyCRSODkyZOwWCzYsWMH6uvrVW23VNKnL3/pZUC0eetgK971zXcBDHDmzBns3L9TcdvfffZ3eY/1vqcX7/rGu2AwG/D9678PzykPbOtt+ODzH4Rt3ZX+UIZhZIcQF/IqFUc+DQ0NsjdkS9VpBijvmp9cc3woFML4+DiampoU31vyb7V6VrUit2bJsmzBm9RCdHd34/z583mPt7a2LopZisAKFEWli5qcGBIqPT6KQEyr3W63JF04OztbleZ/8eSKNWvWoL6+Hps3b9Z8TD2iSN7/ZDKJPXv2aPJWrOY8xVLwnPbA/j92yWNvve+tMBiLn3/QFcS5R89JHtv7ib14+5ffDoZhcObrZ+A55QEARCejOP6V43jHv70DQOEqSCWv0tyWBrfbLQzHJRfzhoaGiq75VZpKW7Bls1lJxC5+XOwDOzY2BpZlhXYRsWCqmdpSSeQGDC93Vpwo5lJIDAmVFkXx5IiWlpY8s2wibJUSRfHxxWLs9Xp1RQLFqkjFJBIJ2O12xONxrFu3DpFIRLPZsNboLZVKCW4oDQ0NmqzY9ESKPM8j6Azi95/7veTxNXvXoOePelTt46UvvIQMe+XvaOuw4a33v1X427z0zy9Jnn/h2xdw/VevB2PQd9GUa2kgUSVZq5ydnUUkEgHHcchms5IUbLnnKlaCSq/HKe1fyQeWROyxWEzWCYmIJZnaUo0bQbn06UJXh1aaFSeK5KKmRgwJZrNZ9UVeDiVRzJ0cobR2RrbXc6EpJIq5EyRyxZhEYFrvptVEbizLwuFwCG0d7e3tiEQiCIfDmo4FzP1N1USmHMfB5XJhZmYGmzZtgtFolFyAcpvmiRVbLoVEked5RCej8JzxYPrM9Ny/Z6eRCqbynnv1P10tubgo7XfyxCQu/eSS5LGmzU04+W8nEZmIwP1rN9LR/M+n54wH6/avK1uKUxxVkkb5QCAAr9eLdevW5c37y7Vfq6urW1QX08U0SxFQtlyLx+N5Uy+MRqNgfRcKhTS3i6gld8AwiWiXMytOFKPRKIaHh8GyrNBaUYxS03O5wiQWo8bGxoKTI+S2L+XYBDK5otAECb2iWKjQhgiT1+vNa+soV9VqLqS/cWJiAhs3bsShQ4eQzWbBcZzk75+b1hodHZUUojQ0NCCdTueNvAKA+Ewcf/j8H+D+tRux6VjRc+68phMbj2wU/ltJLHiex+/v+33e4+MvjmP8xfGCx5g6OYV1+7VZ7WmF53nZ9TQ5+7V4PC5bpblQF9mlMEtR7K8rblPgOA5+vx/hcBgejwfRaBSZTAa1tbUSx55Sb0Ry+xSX+yxFYAWKYiAQQGdnpyoxJJR6d0siPZ7nMTs7C4fDAZvNpnqcUTlFMRAIwG63w2KxYPv27QULWci2eipec8Utk8lgdHS0YCVruUdHiQcL57ayyB1HLq0lLkSJRCIIBALC6Cjxhf3Ff3gRl398WdX5GmuMuOaBa1R9rhy/cGDy+KSq/ebS0HEl9Vnt5v1C9mtyU+TFF/NCvX/lZCmIohImk0m4udi6dSuA/HYRr9ebdyMiLuxRQ276dLlPyABWoChu3LhRt8DoxWg0IplM4sSJE6irq1M1RklMOWziwuEwhoeHYTAY0N/fn2d/JUcpkRt5j8XFO+vWrStYyVrK6CjxduTmw263SwYL6yG3EIU4+69evVrSq+Z501NwP5ZVFqzZvQZr9q7BwJ0DaN/WXvD5BMcvHLrO++p/uhq975mrWF5M3qdKUaVS718qlcL4+HhFRkVVI31aSbeZXNFV0y6itRUnVxRDoRAVRcoV9KzN+P1+2O12sCyLvXv3qmoxyKWUSDGdTmNiYgI1NTWaxynpPS4RRY/HA6fTiba2NlU9luWYpxgKhXD58mXU1NQUjMRLiUJIylA8sYG5l8HTf/Z03nPXvXcdOt/XidUDq9Gw6kpRj1pIM74W+t/fjwP/64Cq55qeNKHmCzVgxhnwG3ik/ikF7g71N2DlEFyli3k6ncbp06fBMIxkVFS5osqlHCkC8kUwcijdiIgLe0ZHR4UWLPG6OsuyeelTKorLDL0XQyIQau/8SJrSbDZjcHAQFy5c0CWIwFykqFWc4vE47HY7QqEQWlpasG3bNs3H1SNSPM8jFArB6/XCaDTmFe8UO57eCRSpVArnzp1DJpNRHQnrQenzs+XWLThy9AheuPcFyeP+5/245u+ugbXdKphQk6IeY8wIboJDciSJ2EgM4VQY/Q/2o2H9lXPff89+ZDNZOH7hgO8Nn6pzzKSlnxW5mznTkybU/EMNGD8DBnO/Y8YY1P5tLZJIahLGSqU5SVRObMEA5ahSPICY2NsV+65WO5JbTPtXasXJXVePRqM4f/48OI7DE088IXyv4vG4pmxX7nnv378fHR0dePrpp+FyuXDXXXfB5/Nh3759+O53v7ugY75WnCjqhVSgFvsSBYNB2O12GI3Gsl2cjUaj6vRpMpmEw+FAJBJBb2+vkObTe1wtYhwMBjE8PCxEUYODg5qOp0eEk8kkRkZGEA6HsWvXLk1rxXoo1JKx92/2gktxePGfXhQeS4VS+Nkf/wzv+Pd3IO6NY3ZoFrOvz2J2aBZsON/Q4YeXf4hr/79rhVYRa70VV3/ualz9uavB8zxefuBlvPqQsm9qa38r3vL/vkX4b9OTJnR+rhOmSRNgBJAB+BYeTIgBk5FZC0wwsHzegjPRMzDVmtB3Wx/MVuUIv9qOM0pRpXgAsbjwJHf6BWlnAOYuzpW8+FY6Eq3EhIzcdXW/34+rrroK8Xgc4XAYP//5z2G323H99dcjmUxi9+7deOyxxzQd4+GHH8bAwIBQaf6Zz3wGn/70p3HXXXfh4x//OL71rW/hr//6r8v6urSw4kSxEjMVgbm0gt1uB8Mw2LJli2yFlt4LiNFoLNoSIp5c0d3djcHBQTAMg9nZWd2pV7VN+JFIBMPDwwCA/v5+mM1mDA0NaT6elkIbcXvFmjVrUFtbW3FBVMOBTx8Al+Bw/Ohx4bFkICmbWpUj/HoY/AyPmCEmiYLIhWrb326D6zcueM97Jdvd+oNb0drfisauRqE30fSkCbV/WwsmMf+Zm/9TMv7Cn0HDuAG//dRvAQATr07ghodvUHxuJS/8Wr4vSgOIyfSL3HYGm82GZDIptPNUIqJbzJGiFkgF7E033YTh4WG84x3vwJ//+Z8jk8lgclJbEdj4+Dh+8Ytf4LOf/Sy++tWvgud5PPfcc/iv//ovAHOzFD//+c9TUVwKKIkiEQSe59Hb26uYb9eafs3dNplMyv4unU7D5XJhdnYWXV1deZMrSq1cLSRSJEVLhguTSfMsy+peiyyWPpVrrwiHwyXZ8GlBjXAf/sfD4JIcTv3bKV3H8L/kl6wJchwnpAs90x50/FUHvJ/wSuzixl4ZQ/dN3ZK/fc0Xaq4IogZCuDI0dvwP8m0f5O+02LxJxYgnyOe2M0SjUbhcLgSDQczMzCCbzRaMKvVQDVHUa7eml3A4jIGBAQBz14fOzk5N23/qU5/Cv/zLvyASiQCYaw1ramoSrosbNmyo2ndZCSqKKskVxUgkAofDAY7j0NvbW9QgtxRRlBNk8aDdTZs24dChQ7IXkEr0OBJHmFAohN7eXrS1tZVlnmIhctsrxFWsi837lGEYvO0Lb0MmmcHZb55VfJ6lwYLWgVYYa4wS8bE/bZeIoslkkkxr6O/vR+Z0Bhe+dUF4zpmvnUGmL4OWHS1C6rV3XLtfLgcOx3DFgzJgDyDqicLaasXUySmM/HYE7mNuTJ+bRtPmJhz+98OwbbYV2KN+KhWFkveT9P41NTUVjSrFa5Vqha7SNnKVFl2lAcN6zcCffvpprF69Gvv27cMLL7xQhjOsDCtOFEtNnxInHC3N/+Lt1RadiBGLE+n3m5iYUDW5Qq85N9lWLG7iqLS7uxsDAwNlnacoh7i9orm5Wba9YjF6nzIMgyNfOQJrmxUXv3sR5noz2gbb5n62zf27auMqMAyD6FQU/3fr/xW29Zz2wPFLR17kl8lkceqUB3V1Zhy+7y1w/NKB2NS8UUAWGPv3MRw6dghJbq4QhV3Dosaj7fOWQgoXcVHy2Ld3fRsGsyFvDTToDOLc/z6Ht33jbRrfHXVUY1Yj+e4oRZXpdLqg9Rr5qampyTvXatjIVbq6NXf/pUzIeOmll/Dzn/8cv/zlL5FMJhEOh3HPPfcgGAwKlbTj4+OSwqqFYMWJol6y2SzcbjeMRmNRWzg5Suk1JNuSmYrr16/H4cOHVX0hTCaTbsEgYpzJZDAyMoKpqSkhZVnoy16uyI0U7qhpr1hMkaL4uYf+4RAO/cMh2d9PvjqJ177zGlq3tmLN3jWYPjMt/O5nd/0MXTd04b0/fC+yAE6cmMSNN14Zp9Pb24wHPnUVLn3mBeEx/2U/Tj90Gtc8cA1aW1sxccsEOv+zE+aszHQL8ELVqZg65FcUcgkOSMi/xqnfTeFXt/4KzZuace1XrkVzX7PCu6GdavQRFtu/2WxGc3OzsDQAXLFei0ajCIVCmJiYQCqVgslkkghlJQphxKhtySjn/kuJFB988EE8+OCDAIAXXngBDz30EL73ve/h/e9/P5566incddddCz5LEViBoqj1zjMej8PpdCIQCKCpqQnbt2/XXSyjRxSJJZzX60VdXZ2mAcNAaZEiwzDw+/0YHx/H+vXrVY+QKvXuPhaL4fLly8hms6oqeKs5TxEoT29ewpfAU7c8NSc4AJp7RRddAD4Ar/7GjW9f932ccwYQzonS7PYAnnh+FB/90HZc/M6VyO7Uv59Cz809GPr+EF57/DVsx3Zcj+vRiEZBCEMIwQorapAfRYrXE9USGg4hNBzC8595Hu/7yfskv5t5bQa//NgvEZ+OY+O1G3Hg7w+gfbs644JqRopaEFuvrVmzRnicRJVkMHksFsOpU6dURZV6qMaaZTlnKSrxla98BXfddRfuu+8+7NmzB3/xF39R1v1rZcWJIqDubj+RSMDpdCISiaCnpwerV69GMBgsKf2qRZx4nsf09LRgFr5q1SpdMxX19DiS9Tu73Y7a2lrNQqyXZDKJRCKBixcvoq+vT3VqerFGioWYPjstCCIwt3bnA/AiAAcAwRb93HT+xvPU1Znx9n9+O0aeG0FkfK5wATzww3f9UHjOxfn/5bId2/EevAcWXElFs2Al64lGqxGGWgPSgbnK55qWGnRc0wHnT52y58NGpcIdn4njJ7f/REjxXnrqEi49dQlNPU04fO9hDNw5oPjagMURKWohN6qMRqPYt2+fsFapFFU2NDSgrq5Os8AthDlANBoti/fpkSNHcOTIEQBzMxZPnDhR8j7LxYoUxUIkk0k4nU6EQiH09PQIrQ2BQKBqMxXFZuH79u2DyWTCqVP6Khm1tDnk2qP19fUhEolUXBDFa5UmkwkHDhzQdPOhNVIkQ5TD4bAwHkntRPRS7/CjU1FMHp/ESw9IRz0lADwKQO0o6oYGC+6//y2oaazBOx95J3783h9rOg8ilDeYbkAD14AQQjiGYxIBzSQyGLhjAK39rWjd04qsLYvj9x9X2iW6b+oW/n82k8UzH33myprnPBswhi7HH/Dax17DMx/rxHt/+F5svnGz7Pu6WCNFLRDf0dyoUjzYmRg6AMiLKgsNH16I9Cxxc1rOUFGch/T5BQIB2SKSagwaVppcwfN8SSlQNQQCAQwPD8NqtQrH9vv9CIW0p9PUItdecfz4cc0XQrXRG7nhsNvtaG1tRUtLC+LxOEZGRhCPx4W0GBFKOVcUPZHiyPMjeP17r2Py1UmER+RHY3mgXhA3b27E+fNXUkztO/PTkVkA4/P77AYgd+m/iIto+1JbnguP5DlPzInkhqs3wPeGDwmfwuIigPiGOE6fPo36+nqMfncUI8+PSH6/AWP4EJ6AERlkYMR3cDd+eudP0TrYigOfPoCtf7wVBtOVM11qkaIW5MZEFRs+LJ5+oXeCjRZW4oBhYIWKovjCxrIsXC4XfD5f3igjMaWKotFoRCqVP1MPKD65opJ3y2Kj8MHBQcmEiFLaOQrB8zwmJyfhdrvzTMLJ36bckWIoFMKlS5dgtVqxd+9eWCwWsCwrO8EhEolIXFHI3TsZHaXl4nD5Z5dVNe6vA1ADQPwJsQHYBeClnOd+6UvXCP8/4AjgsT1XHEX4+ef/VvT8AQB3Khy3kCCKKTamCgDO/MWZgr/vghtGZGAADx4ZdMGNcXTCN+TDMx97Bi9/+WXc8r1bhDXHSkeKla4O1Uqx4cO5UWUikcDo6KhkrbKcKIniYpqJWQlWpCgCcx80t9stNL1v2bKl4B+7HIOGyYeZQASJYZiK+nXKEYvFBKPyLVu2yJoOlNrqkHtRU9NeQdK9Wu6AC0Vv8Xgcly9fRiaTwcDAgPAey70uJeNk8d273+9HKpUSUq9ELJVMqUdfGC16/hswhi64EUUXHsGVZugEgFzLBovFgLo6MwKBJOJv+oT1Qx7AEIAfyez/DQAbP9iLmV9NFIz0JMdZZZG1oSsFN7qQgRE8MsjCCDe6JL8PuUJ47n89hzt/NSfhlY7keJ6v2P6z2WzZxEMpqnz11VdRU1Mjmf0pjiq1LAvIkTtgOBaL6fZvrhbf/OY38c1vfhPA3I1wV1cXnn/+eU37WJGiODY2BrfbXbDpPZdSoyZxpBmNRmG328FxnObJFaUi9kbt6+sr2FpSymvOjfiCwSAuX74sSc/KoSctJCeKLMvC4XAgGAxiy5YtmltoxPsWrwn5fD5hJmckEhHm1ont2IhY1tfXY9N1myRN9rmIU4pvgxFv4G4cmxfGDIDTOc9n2Sxuu+0nAIBGzEWYLAD50pc51gA4dN9erPv3P4LjGQcufuciRo6NgM8qR7yH7z2MutV1uPDtC5h4+YrDiMFswOF/PIydH9mJb3R9o9Bblwf/lkP4dXItas8chxtdGEe+G0rIF8LQ0JCuqFwrlY5CK5naNBgMMJlMWLt2reRxlmWFzyVZFgCQN9hZjefrUhww/PGPfxwf//jHkU6ncd111+Hv//7vNe9jRYpiW1sb1q1bp+kOqhyDhpPJJC5cuIBkMone3l7NXp1600kMwyCVSsHtdsPn80kKiApRSqRILOLi8bhggyeO1Mp5TLEoinsqC6XDyXZ6EE8YEKdfiblDJBK50ui9Jov2t7Zj5qUZ2X2JU4pGZPCp7iCOOdVZZ4Xmf4rxQQB8hofRYsSWW7dgy61bMPHKhKRKNZdXH3oVHz79YQzcMQHygHwAACAASURBVIDZoVm88eQbSEfT2PmRnWgbbFPcTokPn/0wmnuaYf+FHT//gLLB+Fs+/RbhhsPn8yEej8Pn80nGRTU0NJRswVZpqhHlymGxWNDa2po3+ULPPMXcloxKtGNUinvuuQfXXXcd3vOe92jedkWKYl1dXVUHDSeTSbhcLvj9fuzcuTPPFk0NRGS03n1yHAeWZXHy5ElVaeLcY+p9n3iex9DQEBKJBLZs2SJpfi6EnvFRRBQnJibgdrs19VRqpVCqNteODZi7ILXf347/vvG/ZbfJTSlu+cgt2PWjGM7nGH7r4e0AjgBgAMyemUVLbQvi3jhi3hhO/tvJgtsm/Un8/r7f413feBfaBtvwts9LXWviM3EYbhyE6ZPXAmtXAZ4wuEeeB377JrLclZuamsYafODYB9DcM/f3T8eUlyA2v3Mztn1wGxiGkaS5Ozo6FMdFkaicROaLpTJyMZmBGwwGoXiMIJ6nSMQyN6okI84I4XB40UeKAPD4449jZGQEjzzyiK7tV6Qoltq4rXZ7lmXhdDrh9/uxceNGsCwrGXejBdL8r8V3cWxsDGNjYzAYDNi7d6/m+Wd6RJG0V0QiEaxfvx47duyoaHsFAMzOziIejyMSiciuU5YTrdWnDMPg/MPnJY+FAQQBtAIYRye+g7vRBfdcSvFzLtz+gb48UTyMuXTpFOYqVb2Q+IFL6MZcYY247OL5u5/H81BeW7nhazcgPhPHS1+8Utbz+vdex8AHBrDxmo3CY97zXpz5xhm86QdM998Mxjr/Xq9vhOn+m8EBwLNzE1IMJgPe8933oGXLXEYkYA/g2Kev9EGKsayy4B0PvyNvDZphmILjonIjc57nJW0NDQ0NCzKbbzGJohxK8xQzmYywhk6WWhKJBB544AGsXbsWmUwGQ0ND2LJli6ZWrWQyiWuuuQapVAocx+H222/HF77whbLPUjx9+jQeeugh/OEPf9Adqa9IUdSL2vUuIgwzMzPYvHkztm7dKvTG6UWtQIkrO9esWYNDhw7h4sWLutZmtAgA8WSdnJzEpk2bhJYHPe0VakUxHA7j0qVLsFgssFqt6O/v13SsSlc3AsClH1+C8xkngpgreBkCMDb/uwYA/8/+NRg/hSvra1kg871hbG+24GLgSqFLEMC7RPtNY04YzwLI7WD9IOS/2HYAxzA3VvFmAOvnH99w9QZs/9B2ZNNZvPmjNzE5P8y4BsBv7/kt/uzFP4Prty6c/cZZYX3R8vQnrgjiPIzVAtMXbwEHIPvsEK7/1+ux8cicoHIpDr/4yC/ARuSLd448eAQNHdLUerEhwEqRubgwihSgkFQhiSz1DshVSzVmKVaif5hE3w0NDZicnMTu3bvBMAweffRRPP7447h48SIeeOABXL58GS0tLfjNb36jar81NTV47rnnYLPZkE6ncfXVV+Omm27CV7/61bLOUnzkkUfg9/tx7bXXAgD279+P//zP/9S0jxUpinovhGazuWC0Jp5csXHjRolZd6mVnMVEked5eL1eOBwOtLS0SCImvWlQNe+TUnuF3+/X9XrVvE+JRALDw8NIpVJC5ezLL7+s6TgkAtGKlhuFodNT+MInf43zAOSmzkUAfOPUNL565wCGf/SmpPDlmgCL1zHXlnE9gJ3zj4/fOI43PvEGEmsTsHqs2Pn1Abzx7AaI65rHAGzOOdZrAH6CK9HldwB8HMDqhisRGmMy4OXeZnx/XhQbALQ6gvivzq+jmcuiFUAbgCZgLmUq9/6YDDB96VYw/3wrzph51IQ5bFmVxR/u/wO8565Ev6Ti1o0uGK+/Gtv+dFvevvTctMi1NYhThZFIBDMzM4jH44jH43jzzTcl62rlEprFHimqgQg7wzDo6upCR0cHenp68Dd/8zcAoKkan2EY4W+STqeRTqfBMEzZZylqHXgsx4oURb2YTCak0+m8fqBMJoOxsTGMj49jw4YNsutZ5SjUUeqTFDf97927N2/GWiX6DcXOOy0tLThw4ADM5isFFHpvAgptl06nhXS03MiqalBMFFMpDo8+eg5PPvkGLlyQL64R4wHwrSffwB9tXCVp7F8N4AMAugDBiG38xnGcv+88Mta5v2VifQKv3XcefQDOPbtB2PZVzIlwZP7ndZnjJgH8al09nv7J+9DSN5fe/NKXXsb3/8cuPIds7+by/x6f8YTRtF6+apoMOY5yDJ6fMcNzdlwyQkvaxG9C8CO3yP4dyxVtyaUKM5kMzp49i/Xr1yMSiWB6ehoOhwOZTEaYq0iiSj1epctBFAHpdSsUCkmqXcXfdzVkMhns27cPdrsdn/jEJ9DT07PoZikCVBQ1kStMJCU6NjYmREmVskSTE7ZQKITh4WGYTCbZpv9C25aCnPuN3DH1imKu6IidbzZt2qSpWKja/PVf/wpPPXVJ0zbHeCA6EsZbAYjjry05z3vjE28IgkjIWDOwfOINQCSKb87/FOONqRj8JgPaAfzwh2/goYdelfz+EMZwBG68gC4cz2mfePaR5/G++2+GxVp4/YfjGVxgpFGluOKWMWSw6vIZsLg2b9tKDzBW6kslXqVkeDXxKs0t6ikk2Asx1qnShMPhktrHjEYjzp07h2AwiNtuuw1vvqnmU1p9VqQoljpTUZwyXL16dV6UVAi9X3SxsbfWPke9ApVLNBpV3V6hdzqHOFLkeR4ejwdOpzNvsPBCUShSTKU4/Pd/X5b9XSGBAeaiu1cBdADYDWAf5qzZGrsaEXLPNV4k1so33nNrE2CgXHijhNVqhMmUwiuvjOGTn/x13vkewxOwIAMWRlyPuyXnfX6+mOb9X7wFRlORaC4n1TpXcWsCY8gCNRZwV18tu1kl1+WUagMKzVUk/X9iVxlx/19DQ4NwHaj0mqLegeVqkfuMl6slo6mpCddeey1eeeWVRTdLEVihoqgXo9GI2dlZYZFZa6VjKX6FRqMR8XgcFy9eRCwW0zRFQu/YKkIikYDD4UAsFlPdXqE3fUoKbXw+H4aHh9HY2FjxilItFBLF4eEAMhnp7zYBeB/GcHQ+XZiSERgxE/M/UQDXAoIgAoDVY0Vifb4wWj1W9AIYVvkaWgF09zTi7z6/Hx7PLD74wd8hlbpyA2OxGPCV6xlYnsnANG/JdgTuvHMmwnjH/TfDUChi9Ej9XpvuvgmRW25D/YWT4K6+GtmDB2U3q2SkqNXizWw2y7rKiPv/3G63MEicDCKOx+OKbkelnn+107OlRIozMzMwm81oampCIpHAb37zG3zmM5/Btddeu6hmKQIrVBS1fkCJPdnExITgnZm7bqcGksbU+mFmWRYzMzOIxWIYHBxEe3u7ptegN1JMp9NIpVI4e/Ysent7NR1Xryim02kMDw+jrq4OO3bsWJS2Ukqi+OabPsl/715Xj/fNJPAWzg3TfLqwRkFgcrkM5CUUB74+IFlTBABjwoiBrw/gAOb8ToO4Ut2ay/sBkHKWGz59ADtu3YH77vsdAgFpwURtrRH3PsPjt/P9k2kY8UKOJRsANNWaYH52CP8D4JpPXoum+YiQrCkCAJ9g5/oXyTn88v3ovHrutbM3XINCVDpSLHXfSv1/qVQKDocD6XRaaGkgbkckoiy1pzKTyVT0RlFpwLDafuNcpqamcPfddyOTySCbzeKOO+7Au9/9bgwODi6qWYrAChVFQH0VISliqaurw6ZNmwBAlyACV9Kvaj/MHMfB7XZjenoazc3NaG1tlaR01GIwGDRFiuL2CpPJhP3792v+AmoVxWQyCbvdDp/Ph46ODl2zI9VGFiT97fF4JFMx1BRUFPp9rii+7fZ+fOzv9mPswSZkHvs9eJ5DFkawMgKTS8bAADk2bBvm1w3F1acDXx8QHv/j+ec9BkA6n2KuGrVb9N+XZhO4/70/xnPP5T4TCIfTeAWduB53K6Z86wF8ck0tvjkSxeyzQzg+HzXuunEQ7yICOd/Qn53/3cfe/Bga1qv39620N2kl9s0wDGpra2G1WrFq1SrB8Yj0VEajUUxOTiIajQrRpLiop9CoKDG5vqTlRk4US0mf7ty5E2fPns17fLHNUgRWsCgWIxgMYnh4GGazGdu2bYPNZoPX6y1plJLaSRvialbS2jE7O6v72CaTCSxb3NyZuMKMjIwIrjBnz57V1eOoNjrlOE7o6ezt7UV9fb0ut39SoFPsguLz+XD58mU0NTVh8+bNSCQSiEQimJycRCqVgtlslhh919XVSfZZ6Gbq0iW/5L/7+1vR4L6InZ0+cP/6EC7+6+9wYaQZDejEPuT7morXHZEF/uajjXDzXTh9xozps3PDhjc8u0EQQSXWI18UXQAeAXAT5lK6f/svxxGLF/4sHkenYkR7DYDHR6KYzXn8/LNDsD77K3wEbqyFFU2rsmBuvxY9X/4T2Nba5HalCMdyCFwOYMI5AS7JYceHdmjavhDVWPMTR4JKPZWkqCd3VJS4qEfO1HshBgynUilFz+LlBBXFHMLhMOz2ubL03MkVlZ6pmM1mMTk5iZGREaGwhHwwS6kgLVb0Uqi9Qu9xi0WKxHGHCD8xZh8dHS1pLVLpQheNRnH58mUwDINdu3bBarWCZVnYbDaJSwoxVI5EIoJTjtjo22QyKZ5fbqR4VWYEdbd8BGBZZA0mXEj/KcbRCQZz4sQDIMOWxIUt3Pz0Q9N/ZtFssWDwF08j1vs+PHnTk/DlHEOO6zHX3J/b1B8D8BT5jyKCWIg1mEvvyllRXHkdHIwAEGHA/88xTLx/I2LXXAOLpRbf+97r+I//OA+GAT7wgW340Ie2w5jiMHtxFjMXZzBzcQazr89iZmgGfHruBqR+bf2SFkU5Cg0gJkU9xH6NGNOLR5gtxCzFxVr1XU5WrCjm3vGTis50Oo2+vj7ZNEE5RFFOYHiex/T0NJxOJ1pbW2WrWUs5diFhI+0VdXV12LNnT15quJR+Q7nmXvFrXb16NQ4ePCj58mlxtFFznizLwm63IxKJSIqEtBgqi+3EyJSMkydPStaIamqscDiCkn1tmbgIsCyYTAbIZIX5gcDcF+8WAGsB/BLAEbhhwVxhC4O5v5URQJpl8ey9j+OG5/4P/vhnf4wf3vhDhFyFMwYmAO8GsAPA/wB50VypTM//yHHldcyJPsPz4Nk0ov/xNJ54KobvHvNg0n8la3Hhwgt48LO/w59leRQyQIx5YojPxlHXVp6U4WIeYCz3GSSzPslUFp/Ph3A4jLq6OklUWS6jdLlIsRoOUIuBFSuKhHg8DofDgXg8XrSisxyDhnO3J/MFGxoaChbwlBIpym1LIicAecOFy3FcOZEKBAK4fPkybDYb9u3bJ5smVRLTYuTe5IinZXR3d2NgYKCkVhyS+mpra4Pdbse2bdskBtWTkwFwoib3hgYzcOQg8LAFmUxSdn4gMOdU8zyAF9AFdr6wJTMfKWaRRRpGPHS6Fv9685N44ol34/af345v7/p2wbFPhE2Yc615EcAfAKj9K+4FcB2ALIBnMWdNJ6ZQewl5HZiPFLNgwPEGfPa3Fvw35GdLhrM8XgRwW5Hzmr04K9jGlcpiiBS1kNtTmUgksHXrVvA8L9ywTU1NIZlMwmQy5RX1aH2tuX2QC9EXuVCsWFFMpVIYHh5GKBRS7Y5SzvQpWbO0WCzYuXNn0UXzcokiKWhR215RDlEU9zeS9Vk122mBiCLP85iamoLL5cL69eslVnvlxGAwSC5Sq1cnAVyxmuN5Hufr6mD5p39G7B9/oTg/sBbAnwD4NtOJ6/m5wpZZWHET7FiPCL6FPXPC89IEensfxdatLTh8Rz86fu2C2Z87ghh4+5ffDi7JYfrcNMZ+PwYEUziCuarT72LOjDyXtfO/H8Rcu4aYOwD8urcZL9sDAIr3Lx6fL9B5N9zohhUeJPCkQm+mGAfmI8ucxy1tFjRtbUL7tnaYmk1li1aWuuMM6VM0mUx5RunpdFq4YRsfH5cYpYujykLFcxzHSdYPl8qEjHKwYkVxZmYGLS0tmiKIcohiOBzG2bNnkc1m89Ysi21biiim02lcunQJPp9PU3tFqenT119/HdFoVHVfpZ7RUWS7QCCAkZGRivY2KhXaGAzS9zIa5fAfj07jrbbNmMbb8p7//7N33mFyFfT6/0zfvtmezZZstqZnE7IpEEoMCAImAhLEkiCIqIBY7s8bVLxyFUXFK3IBRUUpGkoQCBCJNCEhpJJOtvfed8pOnzm/P2bP2TO97GyIN3mfZ59sZue0OXPOe77tfeW44isL+dqty1i+8ikA/sUT6Cbiumr6OEmeRCr19SNSQ8/1GTqGR20YgDl0soQ2fvrDLrrzqrjhktnYry5nhbGZglMHaBjJZ89wDkd9tn07BExbmoGTwAE87hYi5GneYPOLYoNOKEGBMjxEKH1eAHMzWXBePjkLc8hemE2Pq4cFNQsk26gh4xBdB7ukSEhu5hztg8+ZnD6NBKFIV6PRkJGR4fXAKwqliyUA0VNRp9P5eSoqFAq/9Om/k5fiVHHWkmJRUVHUJBPrDRs86Y6uri7MZjNLliyJet4n1gF8l8tFd3c3IyMj5OXlRS2RFkuk6HQ66evrY2BggPnz50dkaCwiFhIeHx+XlDGme7YxGCkGOry/bj3FXyd+zwE+g0exxhdZG7IYM7bwxz8upPnW3WhwSRGTJsRc47ZRG+CJ3p6YiN6+7Vaxrncz9z4zzio6+eHE64tQ8TCbQbaefDzNMlo8tlSjTErEid2rvqnSI5TgQAUT/o9JlFBBYOGAQFdKBXApnmadp/Emxhm3LeXyW5ZI/x86PERiYiKpqakBIyGj0einLiPvHA6l+BLOgWOqOB3apNFasgUTShebekShdPkIl1qtRqVSodfrYx7c7+zsZNOmTfT396NQKPjqV7/KXXfdxcjICDfccANtbW2UlJTw/PPPxzwHGU+ctaR4ugrGNpuNlpYWxsbGyM3NJTU1NaYTHy0hy8crZs6cSWpqKoWFoVv5g203UpJyu910d3fT0dFBVlYWOTk5XgLCkSCaRhu73U5zczN6vZ6UlBTmzp37sQ37OwOIZssxCPxx4vdkYANQjqdL9I9f2sfwspm8s7uD8yZIRzkRKQYbnJcjWPQWLqrrBbaHWK9vqvReNpNEEVtl/o9uivgCHsGAt4E2n+VFQm2niMvwnpUswZsUd+/u5BYZKQaLtgJFQvJGlMHBQVpaWkKKe093JAfTe4+J9eFcDrlQujhPCZ6H2hMnTqBUKunp6eG2225Dr9ej0Wi47777qK6uZsmSJRQUFER0jGq1ml//+tcsW7YMo9HIeeedx2WXXcYTTzzBunXr2LJlC/fffz/3338/v/jFL6Z8XFPFWUuK0w2Hw0FbW5vkqTh37lwMBgOdncH0RuIDcbyiqalJ6mRVq9X09wfrFwyNSCJF+Tazs7NZsWIFVquV1tbWqLcXCQm73W7a29vp6emRPttYPSOjRaBI8ejRfq699sWI1zEObJ34XQW4xmwwMUS/jyLu5FPcwhGsqKkN2ZPpgbxJR06iwV6PFL6kuo429lJE18SPHEXAD+kkjTaepoRBkAjVhYqn2EyPzzK+9la7d3d51QyjqR+GE/fW6/WSuLdGo8HtdpOWluaVMjwHD9RqNUqlkoKCAnQ6Hbt27eLll19m7969VFRU8MEHH/DEE0/w/PPPR7S+/Px88vPzAUhNTWXevHl0d3ezfft23n33XcBjG3XJJZecI8WPE7FeBOGeMuVqMPL5O5h6TTIcxO7OYBZSsSAcSen1eurr6/3k7+x2e9yto+QC4XLvRnG5aEkxXn6Kv/jFPoaGAot1h4P8cWMVnWziGF/mCGpcqIA1dPBljvAXlvIUSwKmUfcFUZ8J9nqkGCLR0z2KgICKzhCkKreD+gwq9lCNDheqiRGTUtropohRoANPetZXYGBw0Exd3TDz5k1GLVMhq2Di3na7nYaGBkk4QkwZyiPKlJSUM7bb8nSMRvjWFE0mE2VlZWzcuJGNGzfGvN62tjaOHDnCypUr6e/vl8hy5syZMT+4xxtnLSnGimBSbfLUoagG43tRTYevIXgPpofr7owWKpUKm83m97rZbKahoQGXyxXQMWMqDTqByE0k/NTU1IBNNLHON8YC3/1LTQ3c0JONJ0V6BPD/BL0hpip1E6P7YpOKGgEVLm7jELdwmNu5kpPkBSTAYIQZLRmK+/NbdqJCAJT8kysCds+KkNtBCbhIBtwTtUcnKp6mhBfx+DOGQlub3osUpwNarVZKGYopWKfTyfj4uDTaIMqw+bpgnAnC9KejXukrmC6Wf6YCk8nEddddx4MPPujXyRqr6fd04KwlxanaR4kXhxi9tLa2kpOTE9JGaqqRom+UarFYaGpqwmKxUFFRMS1Fal9yE+t4Y2NjVFZWeg0Yh1ou1u2ZzWbq6+vDjnOcLlIM9L357/++kLq6YY4eHZDqaEMkko2FdymhjCL+63srueGX+wOs0eNkv5k2EiYmFAU8M4IKPJGkcuJHg5tH2YELJWrcAUciILxNVSS4hDYp0nMBSYSOhD12UCoEnICCFmbyFJvJpY1/BNiPQPuYlqZl2bLoatCxwjfbo1arSU9P92omETs2TSYTIyMjdHR0eHVsilGlrwvGdKfxPw6DYaPRGJMesQiHw8F1113HF77wBa699loA8vLy6O3tJT8/n97e3imTbrxw1pJirNBoNDgcDsk5o6mpifT09KDD6HJMNVIUlxfTPiMjI5SVlUXtmhHLNuXD8GIdL9Q2Y3XmEMktUvIVEUv6NJY0VKD06cyZKfzrX5/nH/ds5TMP34dW6s9USMT1+V8GXt/P6KSGNiwk4lZrwOnEjYJdLOUZZrKMPr7ChxNKN6CaiBw95OnfPBNujjBSvEsJNlRocSFMCA8oVAryl+cz+xOzmf2J2cw8bybvv/Y+qX2ptL7Zys53BviUawcKBD7NTh5kM/t8xlG8U8QeYv9G5bdJ/eRFfH59Oba6YY7vaGa0eZSTz53kVOYpknKTuPC/L2RmHAkzEusoecem2DAmumCI3a/9/f2SC4ZIktNdozxdpCiHXq+P+aFbEARuueUW5s2bx3e+8x3p9fXr1/Pkk0+yZcuWM8Y2Cs6RYtRQq9WMjo5SX19PQkICS5YsiVitfqoXilKplMSzS0pKoh6viJUE9Ho9+/btC5oWDravsT4AmEwmDh48SElJSVjyle/n6Wi0gcCRgEql5LqsfrRKNwr3ZOoz2CwfeAjiLhmB/TBlPXkaB/16JQq7iZPk8SeWc5iZPMI/UOIWFUVRTajd+DbPRDJHGAl865GLr7mYjbcupeaCYq/zkViQyNLPLKX6tmqemfcS9AioEAKOkgRKEStwcVPbMfY8ns4bD/tKpINtwMZw3TA7btrBLcfjZysUq6+p6IKRkJDg1bEZaGD+4MGDXmMiKSkpEZuRh8J0q8sEuk9MxUtxz549PP300yxatIjq6moAfvazn7FlyxY2btzI448/zuzZsyNu3JlunLWkGAtBGY1GBgcHUavVLFy4MOLB+6lCFAofGxsjNTU1Jgd6MeKLZjZrcHCQuro6lEplyLRwIERLUqImalNTE263mwsuuCCqY4wmXetyuWhpaaG3txe1Wu1lHxWuwSLU98a5Zg1arRa3zYZCcONEEbLr05fAFGN6nqfEL9L7E8s5SR6bOAbA86VODtfUYUixgukFOLgOWhYDkXWcRppeldcj973Uwh9eapH+duONpdxzzyrcbjcup4u/f38XT/Rk86UA21YBVcC3pJTsZIrYjYpmexFOQpcVbIZwVdnoEO+RDPmYSHZ2tpTuF+uUkYyJRIpor+NoEcxgONZIcc2aNUHvBW+//XZM65xOnLWkCJHfuM1mM42NjdjtdrKzs0lLSzsthCgIAgMDAzQ3N5OdnU1ubi4zZ86M6SkxmnSmwWCgvr4enU7H3Llz6e3tjfoJN5qLfGxsjPr6elJSUqiurqa2tjbqY4zkXIo+im1tbRQUFLBixQovpQ+xwUIQBL9BcFzQ/Fozja82MtQzhPtyN3M+OYfshR55wN5eE4+/6eREws0ssdYxRCK5WDhMCaeCEE8gAtvEMXQT8aBvpLeZo/y91Mn+C8Eino5UPVz4quf3lsV+ER7AFnbzLiWMUMRXFJ18S/k0KrcTt0rFrrK1/KA++macZ55p4ZlnWrjgghmUGU9gO27w23Y7RXwKjzB5EmCiRGq+EVByhKUcZ0nIBh4RF993cVT7Fw6nw8A4mAlxsDGRcHZRIj4O26ipDO//u+GsJsVwsFqttLS0YDAYKC8vJysri+7u7rg2ywRDoPGK2trauIqC+8JisdDY2IjNZqOyspL09HTGx8enrYFF3sEqNtGIztzRIlykKKa85RJwdrsdhULhN9/mdrulJ/zmvc00b2um/41+nMbJ8/7+vvfZfe/7DGcnciRVy/4OA06XAOTwus9s4dLqXH54zwVoNEqe/Mke/n6wD/BPUQJ8mSNSs40LpfS6GFXeUyMjRBEaB0mX7OJa7Voy64bJpIheikigkx0+UWfBXCPqBhcKwY3C6WZtw5u8PYXa4549Y+yR/V+MLvOB2yZeK9lQgn3IzsBxLU8ZJwf/IyFDgKueuIqqa6ui3rdQmE5SDEVaocZERGUZ0bJMXtOUZzGmO30aiBSNRuM5Ujyb4XA4aGlpYXh42M9hQa1WBxxRiBTBRjpEGI1GGhoaUCqVft2W02UfJR7vyMiInzj6dIyROBwOmpubGR0dpaKiwqs2E2sXabBIUSRet9vNwoULw46r9PeP8+Tjxzi1pxtXyyiKbhOpQCqQBugAJ5O6oL1DFggzo3jk6AAbN77MN7+5nEdfu57+q7fxvowYRTLawm7UuKWu0+0s5cZvbODiVC3H/zCIffQ9OlICn3+zYohymddiPvCpAPXF/60t4TqVCq3CDYKAQhBIUApcoeliny22TtVAEM9g9uJslt+33POQYTBi6VzEWIsNbYMZnm8JuY6kJknySgAAIABJREFUkiRuePEGMsrj31X9cZFiMASzixLrlPIxEVGJZnh4eFrGRAKRosvliks99N8BZzUphrIbmj17NhUVFX4Xjth9GiuCkaJ8vKKysjKg+G68jYblyjCzZ88O2LgT62hFIMiNhWfPnk1VVZXf9mJtmPHtPnU6nTQ3NzMyMhJR9yqAscvI1Rf/jfpBc9D3aPA0iNh9XpfX6Y4oitAkqTGNTxKYyyXwm98cpP3ZN/hVjY0XLsjlt/tUE9GlB77p1AdZwhfmZfGDzYvgB+djeucqcl/eSL/a5LdfCvMMnHhf0JoV56M4+j5Ouw0BBUMkeuTfXF/i+wUtXDWwB1xOFEoF3/vllZRlXsCmTa8Rj34lcQx78abFFBdP2j25ayaj8Ln/MZcTD56gc6u/ylNSXhLLHl02LYQI0zsAHy/CValUAcdExNrk6OgoHR0d2O12EhISQo6JRAPfSFR0nzlbcFaTogi3201XVxednZ0UFBSEbGSJh6einJzsdjstLS2Mjo6GtbCKl32U3F7JVxkmntsUIa+NBjIWlmOqSkOCINDV1UVHRwfFxcURd+gOnRpi29XbaA4T9QV6HJKPQThR8QSbWbjjTu796QHeeqvH631P9z6J9hUXSxRqqr/+IA/uV9Lw0TAKm5N9gr8CzYG73iIjI4H16ytI/MQa7kt/kDufvR2LazJboXSocO//BC1aJZV2zwNMIZ1UJLpxfuPraP73f1G43PyWnZLjxvruIr6mSuV/2YHK7SZhyxaueeUVfvazS7j77ncj+szDYRBo/kcz1V+tntxXWZ1t1qxZlD9YzkNbH/JbtnLeMAXPPkXfYDeqNWviaqArYrpIcTprfkqlUiJLUSQ9kjGRaHwVfSNFkRDPlOH66cZZT4o9PT20traSl5cX8mYtIl6eivKotKSkJGDUFGjZWKNUkdyGh4dpbGyM2F5pqpHi2NgYDQ0NJCUlxU16LhAUCgUGg4F9+/aRlZUV0bkUMXhikG2f3oZ1xEou0BN2Cfl24XphctBdgYs5Qhsv/+Yk5bOz2Zs8yPi455zJu00RnHz06It8ODHHd+dti0mv6+bw8QR+rS/GMWEi7HYL3HzzP9i27TOsXTubjedtRL3jH/yo60U6k6HIBEsOVvFqy2JaytOobBqblFx7z4VijwrcbpQIaH0ad2a4xvHYAAu4bTbsL79B+dgKvk4n6T7dqTrCq/L44hHgm++0c+HxAXIXBx7M3nv/Xr/XCunk6g/+inKXA+H552l+7DEay8qwWq1eDSmpqalnpG7p6WiEka8/kjERXzcReVTpmxb1JUWr1erlrfh/HWc1Kba0tGCxWKLy3otHpNjX18epU6fCRqWBlrVa/Y1lI4HT6aSxsZGkpKSo7JViveFYLBYpJRxIBi6eMJlMtLd7lDSXLl0a0dyoIAi43W76jvTx8rUvYx31fK7V+JPi+ecX0No6Rm/vuN96ntmyGtPPO6WuSjcqfkcJ27c3+b13iERUCAh4hvCHmLzRbH2hjj//eQHr0tJYd3ic733vQ1wTqVW73cWDNzzMgi8mM/OGT/HZy7/OpvWvI9jsWN0K1rEagAPNY6zBR3LN7UZQqUAQUGo0vGst8dofUdvU7lay+ZER1rKT3/Ikygkh7z+wmQGKmE2nV4PMrowE3hkN/118CMi9/Fnu6LoDpco7Suk/2s+h3x7yei0hM4GrrkxB+YwThdsNDgfFLS3MvOaaic/CvyFFpVL56ZZOtwNGKJwOL8VIHviC+SrKx0RaW1txOp0kJCRIn5/VavWqvZ9NnadwlpNieXl51AQXKymKKcSenh5mzJgR9dwfxJbKtFqtNDU1MTIyQn5+PhUVFVEtHy3kTTs6nY7q6uppm6mSW0fl5eWhVCrDEqKYCnK5XPQd7uOla17Cpp+MgRYBb6oUOGS1vpycJD76aMhvXcnJGpofO4KDIp6asFOqpYTtQboqs7HgwnPROYFl9EnjEvuGi+j9yMaFX1tCaakRq1XBD35wEPCkXV+zPon2Ty4UTz+K+ZXtmF95BdXu3dzypIl97Z7OWbcAx4G5E5JrKN0odDrav/UtsgSBYxkL2PefzdI6f8tOlAi4UXIXV/AqRVSzG+VE1Asu5tFGIkhi36Ljxdyqlbyzz/vxYcYMHWNj/vHkqXEnD2Y8yOq7V7PgiwtIK0rD5XDxxh1vIMg+55RZKWzev5nEuqPwwmO4BQGFVotzzRrpPYEaUpxOp5Q67O7uxmQyTZwf7/nT6Zztk8Plck2rRupUItFgYyJWqxWj0SiR5cDAAL29vbzwwguSN2WgBpxwuPnmm3nttdfIzc3l5MmTAGesj6KIs5oUY0Es8mUjIyM0NjaSkpJCUVERCQkJMXVyRUPITqeTlpYWhoaGKC8v9xPgjTfkdVmxlnfo0KFpGedwu910dHTQ3d0tSc4NDAxgNIaWmxajQ7fbTe/BXl78zIvYDd4tM+d/YymNfSZefHHSNnf79kAWujA+7kA/7iAJJDslJ55mnEBJbk8jjRoBFy6UXlJn69jM174P+3fo+f5jn+LOOy9EoUji+99/zyvt6rDZ6H3mObq+9AVSrrqKQnM3PHBc2sbbgJUihISbufFbuQhrL2JwRjq6wkL+/IN90vu81omHsAHepIQtqNDhAo0W65JVlBza5yX2XUIbu/Z5E//KlbN45ZXrePBX+/n5Awe8/rYV+DGw9+d72Xv/XkouLSE5N5nB44MUyiLQ5b+5E126DvfKlYxs24Zpxw6yr7sO98qVIc+rWq1mxowZXs1p8oiov7+f5uZmXC4XSUlJpKSk4HQ6sdvt00Je050+jff6FQoFiYmJJCYmkpubi8PhYObMmSQmJtLV1cXbb79NU1MTK1euRKPRsHjxYh599NGICPKmm27ijjvuYNOmTdJr999//xnpoyjiHClOIwKNV3R3d0+5LhgK8g5PuXVVT0/PlEZJgkHupZiTk+NVy4tV/zSSbeXm5nqlnkN1rcrJEOCx1w/w02MfMvZDNeljKta9bmfxERcrvr2Ci35yEbzd5kWKIffJ5/9qPOa5gZbeRxFb2Mz5tDGKnlsnNE3lQ/p/2dPNXxf+iasuLuJ7P7+EzZsX8u6TnV5dqS+cmEnTd+s52a7nuE9k5gJ2AUess/jn/bDg/g9ILk6itSSDJ3d3S+8Lpnwjn50suP4qvvjwzTTf/TSux95DmEgPt1GC7yNWd7eRxEQNd/9oDS9tb6SucdRvv1QTH1jbm22At92UW6XBlvk53BM2xPbzzqMvJ4fM+fMjOg++CBYRiUINTqeTjz76SBL4lkeUU23o+TiG66dj/Vqtlk9/+tPodDoyMjJ44IEHMJvN1NbWRrz9iy66iLa2Nq/XzlQfRRFnNSlOV4FeHIK3Wq1+4xVqtRqLJTbvvVCkKMqktbS0BOzwnI55Q71eT0NDAwkJCQGbaKaif+oLo9FIXV1d0G0Fmm/0JUOFQsFzdQ3c3XQYxwzPuddnKHj1szrKryrioi0XoVAoWLt2NoUFqXR1e0eeGRlajEYHTuckFQai4TL8SVEFXAsMU8QXKaKGTjZzNKAcmwN4+b1OXj7/6YlXfLpSP9QCA6E/L+CFiZ8FHWY+6vAeMwnltSjNTm41cOfWB/nGN5ahV93EfFcr/6SEDyjCtyLd22vC5XKjUin527MbOO+8J7z+3g/M8llGXvtU4EL9/vvYJ6JCcR5PuX8/6vffx7lmTdiIMRwUCgXJyckkJCTQ3d3N0qVLpc5NMXXY29srNfSIdcpoG3pOR03xdJKuwWCQ7mFJSUmcd955U1r/meqjKOKsJsWpINCcUyTjFWq1OmaiCJY+FdVvUlJSgrp1TIUUA1lWico3VVVVQVOzsXauilGfQqHAZrPR2NiI2WymqqoqaMFfPqcorxuK6xHPw38d2I/DJ3Pt0CrYmjHKf0mCBUruuqKUex4/hhUoBc4D7n7tKj517Rv090823AQixbnAm0waCCcBNwIW4LmJ16M1AA7lixhOy/Qjn/8rgGJgf4Rei48+ehgonPjxQO/zHpdLoK9vnIKCVMoDzBaqry4nt8vAwNFJMhftphRKAXxqh4IgkHryJEl33AETaU7zK69MmRjBm7TknZviiANMraHndLhYTGfHrW9369jY2LQ12pxJPooizmpSnKqnolgXdLlctLW10dfXx5w5c0KOV4gyTbHAl9hEc2EgrLnwVAf/xYhLrFOKSjShPsOpkKLD4aCrq4ve3l7KysrIy8sLuy1xyFjcV4VC4XfT6h73H3wH6DKN4zA70CR5zumF5+Xz3cc9Atwih7597Rvg8olGA6xrBvAZYC+QBawFDEwSoohYDYDliMUq6gY8xN0MPE/0oxbBYLd7jk6hUFBTk8/Bg73S3y6/dQlr186m/2g/J586Sd0LdXQbijn8xV+xpGQU14UXehGe2+0m7fBhsNtRuFwIdrtXJDkVRBLJTaWh5+OwdoonfPffaDQyZ86cuK3/TPVRFHFWk2KsEElRpVLR3d1NR0cHBQUFrF69OuzFFg+pNpvNRlNTEyaTicrKyog6t6ZKih0dHfT29lJUVCTVKSNZLlpSFAQBh8PBgQMHIv5MRdhsNmw2GyqVKugTaGFKCp0mf2JMH3Pz8o0vc83z1+AUnDirnPi2Qo33jyMfQlhFJxfTxmAAHc9FEz8i/o43IQbD5/F0pu7EQ6SBkJGkZnF5BsvOm0n+iwfQ6qOzihJVRMuArwB/A8Yi2LdQWLo0j9mzJ6OJbds+wzXXvMixYwNce20lF17o2ae86jzyqvNY9z/rpCg+UIVdEARMy5eDVotgt/tFklNBrOnNSBt6TCYTLpeL9PR0iSzj2dBzOtRl5NeOPH0aD5ypPooizpFiDFCpVPT399PT00NOTk5U4xVTIUWREA8dOkRZWRnz58+PONqNhRRFI2W9Xk9ycnJUA/EQPSmKbhkul4vly5eH1SkV99HtdpOQkIBSqeTw4cMApKSkSG4mqamp0pPvvatWcfu772KRnQONXWDd63bajrTx3PXPkfOZHPJm5FFzVw0Hf3vQa3vpeFKHk96AnuaTp9gcUuA6kk+hAg9Rrfnhar63IIu779vHmyf9R0E++PDLZGcncuutr7NXn8/tPg0zCToVq1bMon/ITG3tsN/ytYDYvpID3Ioniu0Is38//eliKisTSErSYrGoGRlx0tQ0yPz5ZVx66RyUysnvYmZmIu+99wWsVicJCdGrF7ndbqxLl2J+5ZW41RTl645XzS9QQ8/hw4cpKSnBYrFIUmzxauiZTnm6YJjKnOKNN97Iu+++y9DQEIWFhdx7771nrI+iiLOaFGP5co2MjDAyMoLL5YpJoSWWmqLb7ZYiUoVCEVX0JCJaUpTbR2VmZlJcXBx1x1uk3acWi4WGhgYcDgfz58+nsbEx7LZ864YajYYFCxYAns9LbJzo6enBaDR6alSpqaxOTeWXNTX88sQJukwmZpgUrH3FyuIjns+m5+0eet4OrmlTjIc8xJEGcZ6vhLaQpHgl8BKe6C8dyARKitPQDZtJHXeSBYhaJHt/upes+Vn8z88v4c/vtvOb33gPuH/yk8+Sma7j2Mkh/JpwKAKbi3d3d7J27Wx27foC69dvZe/eSXJ9C0+0KCbIkoGbgC486eBv7/sSs+dl84lPPMOHH/ZJy5WWlnDFFeWSUsrIyAjFxQ60WiNNTSe9GlNESbFghBgO4s3fvXJlXFKmckx3I4zb7SY1NdWLSOLV0DPd+x7MYDjWSPGZZ54J+PqZ6KMo4qwmRYhcgFocr1CpVOTm5pKTkxOTZFk0NUX5CEJ2djYrV67kwIEDMV0UkRKU1Wr16pxNT0/n5MmTMaVew3WfOp1OWltbGRwcpKKiQmp0CBdhhqsbKpXKgELKJpMJg8HAKp2Op+bOxe12M/bBGCeOnIj4mGYD7+M90gAqWinBDIwAoxM/4u8OoAhPs002nkYXADqCJUhh+NQwf9/wd+auyOeLN8zjr8/VSn/r7DTS2TnZGRusNvmvf7UzPu7gzjsrOXBgWFLIGQE+BFbI3qvEQ/hJOUkkWF2cePIEMy3e39OPfv48G06qUa69GM3Kleh0OiwWCwsXLvSqt4mSYmK3Z6QGznJMt9/h6XauD9XQI35uvg098qYe8bOY7nGMQPVQvV4f1/TpmY6znhTDIZDHoCiNFAsijU5FzdDExES/iDSWFEo4MpYTVHl5OTk5OdI2Ym2YCbacIAh0d3fT3t4esEYZzD4q0IhFpJ+DUqn08k3sO9nHP//fPxnYFXq0wRdi/6W8e/QEJbxFUciGlR48VlMb8EjJBcMYcAwwA4voJPfAbtI+qgBmhtwvBf5NP2lJahqePIH1VBfVLoEPZX97b2I/fCtd5kEzWy/ZCoD8kW8VnfzgxJPoTrjggQew7HgVYeFCL0s133qby+WS6m2BDJxDKc2IDzzTAZfLNe0ycJHuu1arJTMzk8zMTOm1UA09Op0Ol8s1beQYaL0Gg+GMUpyZbpz1pBgsUhQlxMbGxvzGK6aqfxoKcuPdQJqhItHEy5le7ioRrIkm1iadQKQ4MjJCQ0NDSKk73+WmQoa+sOqt/OM7/6D5+ebICn0+SATy8MzdRds9KgB78CdFUdWljhJ2AytpY4hEfstOT1fp+HscDtNVGijXUWp28sF/vQ/AJXgk4MSmlnGgF0/kGwzyrV1CG5qJmUKX3YZq926EBQtCngeVShXSwFmuNCMSpRgdCYIwrZHix6mNGg6hGnoGBwex2+0cO3bMS6EnXg09gQyMzWZzRHrC/1dw1pOiL5xOJ+3t7dJ4xdy5cwPOGsY6gB8Mch3PioqKoP5/IkHFQoq+GBoaorGxkczMzJDNQvGIFM1mM/X19QBhBcnl4xUiIfrOG8aCwcFBtl+/ndGDo2Hfq0nScHvH7bS+2cr2G7d7/a2YSb/AaOE7QSpXdVmDklsANW7cKFAiRNVV6ou5gBVP1JmGJ3XbK/t7uKJBOpACmPBXwPnVfi2Xz3iTgsZ9KK+9NuImmHBKM8PDw7S1tWGxWNBoNDidzrh3cJ7ppBgI4ufmdDpxOp1UVlZ6fW7xaugJZhv17/Z5TQXnSHECvp6KoZpZ4hEpijd5uYVUMBIOtO2p3CCMRiP19fVoNBqqq6vD2sLEGimqVCosFgt1dXWMjo5GbPYr1iJD1Q3B8xmO1I/w0daPOPm3k+jSdVR9pooL7rnA6zMcHx+noaEBtVodESECOMwOTL0mKtdX8pXjX2Hruq2YJ8yHl+Gpyfk+Juh0SrITVaSMOcgAMvA0tchHDq6//Tx45EMpOkxDL6m6qCeGNlSAc0Ko2wF+ijciivHULIMpvj4r+z1j4r1yZBIan7wxhYoD+3i4Od1PbIB/tvL//vljtLgQ/vRnTC9vR33h+WHWGBhi7TE5OZmZMz1p4o6ODgRBICEhweuGLzfTTU1NRafTRf2g9O9IiiLkDhnyz03EVBt6gqVlz7QB++nEOVIE+vr6aGlp8dPuDIZ4eCo6HA4GBgZob29n1qxZEXeUTmXe0O12c/LkybDqML6IJVJ0u92MjIwwMDBAZWVlRH6RgJQ26+joIDs7m/T0dC/SFgSB1n+28t497zF0yntcYbxvnL2/2EvhmkJKPlEiOXYMtg4y+LdBund1+24uJHZ8eQc3vnUjGWUZfPblz/LsFc9iN9rJB74KtAKZGQlseurTVMzLJjc3kZGuEXZ+bSfDHw1jG3bwfqIKh2XyfF2xvow8RRfzH/7phOuEEidKFLhx4Tn/btw4UPEHriANC0cpIYUi1uGJ7jLxjFQo8Qze/4LwmeBRvJVvDlFEKDOvQjpZ8fzTrBZcXIeSS9nEHlm6eAu7JwXFHXYe/uxvMNzu5nvfWxVzx6kcgiCg0+nIy8sjLy9Pes1qtUr1NlHPV6vVSjf71NTUsJHRdDfxTCeBhMsSTbWhx5cUp/t4zkSc9aTY2trK+Ph4UHm0QNBoNDGLeguCgMvl4uDBg2RnZ0dtIRXLSIfT6aStrQ2z2UxFRQW5ublRfdGjJeLBwUHJu7GgoIDCwsKwy8jrhsXFxej1ei8Xca1WS1paGkfvOUrbq20h17Xt09uk3xOyE7AORedBKTk3HCphz0+Kuei/LyKvOo9rnr+GFza8gMvuYiae1pcNj15O5SdKALDoLTTubiTriiyqv1ZN96iA8fZ/SutNSVbj0negbt+DUqzP4eZxzqODdCkavIQ2UinBQRF6PKo4gcbWq66touuDLor7xgn9ifgr31zPZhQh0rEltKF0OVAgoMHN92njfyjiEJ45Td906j8shex74ADd3SYee+yKCD/p4AhUU5S7Ochv+PLIqK+vD6vVilqt9iJKeWQ03aR4JuqeRtrQ43A40Gq1KJVK+vr6yMnJiZsX6s6dO7nrrrtwuVx85StfYcuWLXFZb7xx1pNiaWlp1CQTa6Qozv7Z7XYWLFjg5ZIdKaId6RC7PAsKCkhJSfHqKo0USqUyoocAk8lEfX09arWapUuXMj4+zvCw//C47z761g21Wi05OTl+N77+lv6whOiLWAjRyzvw19C+djaz186m+KJiNmzdwM6v78Q8aGbFd1ZQ8ekKBEGg/VQ7r254FWuvZ3sf8RG7fNZdNO7kXxv3YUDH/6BCM0EoT7HEq17oooirItjX9n+1Yx21UgZhSfGz2f1ohyaVb9bRhjEEKYq6pKIzRi8lXAhcALSnannfGFi79bXXmrDb7SiVSokcYiGJaMhFp9Oh0+m8rqdgkVFqaip2u52UlJRpIbDTIdYdi+1cIARq6BEdLVQqFc899xzvv/8+fX19bNy4kaVLl7J06VLWrVsX9T64XC5uv/123nzzTQoLC6mpqWH9+vXMj9EFZTpx1pNiLKmBaEnRV0C7s7Mz5nbqSKO24eFhGhoayMjIkKLRgYGBmNrRVSoVVmtwcrHb7TQ1NWE0GqmqqpIuMovFMqV5Qzl0Oh2FFYWkFadhCDHfFws0yRoc4x7S93KtnxjKf/7q57m97XaScpIo+1QZXz31VQRBQJuslR4Eht8ZlghRRLPPdoaABwEHRRwLQCha4HzgoiD7Kfce7KII66hne5V4fBQD4UffXM4Sl5vdT/V5RXY7KSEHj4pOoCugS2acLG4PPCnbOUY7c/A4frxHEftlyxmNdiwWN8nJkzOq4r/ycxzuOzhV5ZZgkZHRaKS9vZ2RkREGBwdRKBReXa/BRL4jxenwUgzXAzDV9c+YMYOsrCweeOABjh8/ziOPPMKPfvQjjhw5wjvvvMPatWujXu+BAwcoLy+ntNRjDfa5z32O7du3nyPF/yuIlJgcDgfNzc1+rhlTqUmGS5+KN2mVSsWSJUu8WqnF/Y72KS9YTdHtdtPe3k5PTw+lpaXMmzfP60YWak4xlhELpVrJ9duv5+BDBzn+l+Nh3x8pLvzxhShUCvY/sJ+2Hu8IqW0ipbn/1/tZe7/nZqBJ0jDWOcaxd4/hznSTp8rjnR+/47de38aWQdnv8nGOBGAVsBLPyEcgLLvAxhUfPIVKcHoiWJmsXB5wIR5RgZQEFVddUMTqGTrse7oYf+gQJ4AjzAwY2c0GLkPufzEJ0Tg5GLKAy4E6n2MdHLSSmelpqBIzAOL5DkSU4o+cjKYjilOr1WRkZDAyMiLd+F0uV8CZQHlTSjSiA/9uBsO+8K0pjo2NkZGRwdy5c5k7dy433nhjTOvt7u6mqGjyu1RYWMj+/ftDLPHx4awnxVieRsMtIyeL2bNn+zWZTFUUPNCycpFwebTmu2ysXaTy5QRBYGBggObmZvLy8rzMfn2Xi/e8YWZlJpc/fDmf/N9P0tvTy8k3T5LsSCbBmYBl0MJo8yh1L9RFtc7Wt1v57N8/y+KbFnPiyRO8cF8SuUOnvCKkgeMD0jG8/+j77L97P4Ir9FDDfGBfiL8n4YkMa/Af1fB7754DqHB6RbBywloHXAyorC4Ub7f5EfJC4M8yIpY33fyJIi5bNYtNKwpIGDLT+GojdqM9zB5NYlZBKqMy78nubiNVVR5SFIlN/v2QE6X4L3gTZby9P+WQE65KpfJTPwomOpCUlORVpwwmOnAm1hQjRSgvxbMFZz0pxhOCINDb20trayv5+fkhySLWi16lUmGzTWqnyG2rwomEx2MI32AwUFdXJ5mNhmpOEkcr4j1vKJobp6amsvYLa6XIt+uDLk49fyrq9bXubOWF775A2RfKSL80nZVXfouOlzro2vKe5zjUSpbcssQjWF5bz6H/PhSWEAE+iUdwezfeLhSpeGpzy/BXlAkG3xqfGMHKEepizgG+oVGyNzsJobeWt3zspt7cB/861MfNNy8mX6MMmFKVIyk7iVmrZrHi2ys4+bvDfCR7EOnpCWzPJSIcUY6MjDA2NkZBQYFUyw4WUcaCcMQVTHRAnAkcHBykpaVFGp6XE+XpJq3pXr/BYIiLl2JBQQGdnZ3S/7u6uigoKJjyeqcDZz0pTuXmLK97iDW8GTNmUFNTE3KOcKrpU7PZjCAI9PT00NbWFrHF0lQiRbvdzokTJ7BarcydOzeosbAcch/GSOqG4SBGwzabjXnz5kkuGk6bkz0/2cOB3xyIed2tv2+l9fetFFxcQHJZMppiDUv/uBSNQoNOq6PH3MOJH56g7S9tEa9TicecuBqP0W87UIDHUipcAlupUeJ2TEbZwWp84aBOVpN7QS4zVs0g7bw0rstJI/XhRnRbPRQrFwZwOt384Q9HyQFuY/LmcNOBm0gvSUebHPg7PWt7g9f/e3tDk2LA41Uqsdvt1NfX43a7Oe+889BqtWEjyliIMpZoTqlUSuMLomt8INEBq9WKUqmUmnrEWcp44XSQrnz9U3HIkKOmpobGxkZaW1spKCjg2WefZevWrVNe73TgrCfFWCESjOjwEKiGFwxqtdor2ot2u0ajkf3795Oenh6WgAPtczRwuVx0d3czNDTEokWLIh7nEFvqjUZsoJspAAAgAElEQVQjR48eJS0tjfT0dNLS0qKuaYrp6P7+fkpLS706aAVB4LWbXqPxlcao1hkM3e91e4RB4wgVsHjiR8TqLatZc88atl66le69/vOTboeb4puK6Xhi0tApXI3PF9e+cC2z185GPTE3KN7E7ddciuLvT+C2O3Ci4l2hxGu5QaABT/q34PwCchbk+K7aC7Nmebfsd3cHkxMIDLFLurOzk7KyMj/T2WARpTwVL2YklEplWKKMV4ozkOhAb28v4+PjJCcno9fr6erqwm63e6nMxCo6ANNPir4SkgaDgZKSkimvV61W8/DDD3P55Zfjcrm4+eabJVebMw3nSDFGKBQKTp48id1up7KyMqq8e6yRoslkoqmpCavVyvLly0NKpQVCNKQoTwXn5uaSkZEhDVGHW068WSmVSs4//3ysVisGg4Hh4WFJTD0pKUlKUaWlpQVMCYkuIS0tLcycOZMVK1b43cwsw5a4EeJpxUT2daRhJOhbKhZXoC/Ro2/TR736Tz/1aco+Veb1mnQTv/JKBn73DIa/vMKhhlQK+rOYgXeK9wgeUkyfm47dbg/54GWxeI/rRGOCazKZqKurIzU1lZqamrCpwVhrlHKinM66n6jCk5ubK5F7KJUZOVEmJiaGJcpA2qTxhnwf4hUpAlx55ZVceeWVcVnXdOKsJ8Von9acTictLS0YjUbKy8spLi6OybEimohNPvJQVFTE0NBQ1IQYzXZHR0elml1NTQ0Ax44dC7lMqLphUlISSUlJ0tO0GLHo9XqpYcflcknGwGlpaSgUCpqbm9HpdCxdujRoCkqbokWTosFhik1M4ePA7E/MZvmdyzEPmrEMB9fQffs7sXnOLdq8iKprq/xeN/YYqdtWx6nnTjFwbAAm6pKL8HSS/kH23iY83o8JpQmcOHECu91OYmKiZNyclpYmnZODB3u9trNkSfiHJ5fLRWtrK6Ojo1RVVUWUjg+GaIjS4XBgtVolEY141CjlkMuwiQilMiMS5cDAABaLJaTogO/xng4YjcZzjTZnIyLxVHS73XR2dtLZ2cns2bPJy8sjPT19Wucc5bqo4siD1Wqlr68v7LKBEI4ULRYL9fUe5/sFCxZM1uyczrjNG4J32mnWrFnApAvA6OgotbW1WCwWaSh7YGCAtLS0gK3x6gQ1G1/byN8u+VtEn0HpFaW07Gzxem390+uZc+UcTr57kv1378fUEH1NLBoUbSzCLJgxnIrvvCVA7pJcFty4AKfFiSbJk6ZufaOVAw8eoGNXR1AV8Hwm3T/A87bjwLevXkpmRSaCIGCxWDAajYyNjdHR0TEhYm9j1652r3UtX54fch+Hh4dpbGxk1qxZLF++fFpkxAIR5djYGPX19WRlZaHVav0iSnG5qRCly+WKuIao1WrJysry0gN2OBwSUfrKsYmNPNMV6Qa6B8YzUvx3wTlSDANBEOjv76elpYXc3FxWrVolNbvEKvUWjhTlqctZs2Z5dbFOtXM10D6L0e/Q0BCVlZV+SjvBthlPSyfw3LS6u7spKSlh5syZCIIgGQN3d3djNHpqVWKkIhLlrJpZfPnQl9m2fhumIJ2PBasLWP/0elLyU3jz229y9A9Hpb+13f8MSTsGyLnwEjY+v5HnrniO8b7xmI8jHDpe7iCjJoPat2vDvzlKDBwb4NkrnkWpVpJXnUdGRQanngndkatOUFN2ZRmftTl5ZMek5MCpNB0Z5R4fPYVCIUX8eXl5CILA739/mLvv/gCHrCEoIUGJzdZJc7NROkciSdhsNhoaGnC73VRXV8dk0h0LXC6XlGlZuHChV5bFN6IUf3c6nV5+opES5VQJS6PRhJRjs9vtfPjhhygUCq9ZyuTk5CmnVQPVK882L0U4R4pA8EhRTCOmpKT4jR+ItjaxIBQpin6DaWlpAZtopjrjKFemCWf2K8KX6OJNhsPDwzQ1NUlasOKFqVAo/FrjxWFrg8FAZ2cnRqNRstVZt3Ud733lPcaaxvy2UXp5KSn5nsh3yc1LJFIspJOrP3oS1UcuXM/+jafYzHiUFk3RomNnBwOHB7AORCdBFw3cTje9h3rpPdQb8O8KpYLii4uZf8N8KjZUoEvTsXrQzB/e+J1Ecj0GGw8+eJBNmxaRlTUpK2Ay2bnjjjfYts2f1DdtWszChfMxGAyMjY3R2dkpNZXZbDYKCwspLCw8bYQ4NDREU1MThYWFVFZWBk1F+qZeAb9xIvHBUBAEVCqVV61SxHQ0wohybOnp6fT19VFTU+MnOjA+Ph6xgXMwBDMYPpc+PQdMJhMNDZ42c3kaUY6pkpNv5CXaG0Fov8FIUr2htite8OIISTgvRTniPW8oGiorlUoWL14ckXxVsGFro9GIwWBg2cPL6Pzq8+R11HqNLjjMkxFyxtwMMhZlMHpiNKCsWzQdntFAnaQmMTcRY5sx7oToKwEXDCnFKSz/+nLmXT9PekgQkZOTxJVXlrF9+2Tj0g9/+B4//vFuLrtsDhs3zqOiIpNbb91Bba2/pu1dd9Vw770XotGoSExMJC8vD5PJxKlTp0hJSSEjI4Px8XFOnTqFzWYjISHBr0YZr1SqfLwjVE06EAJJ0YnXjfgwKK9TwiRRinXK6UA40YFwBs4iUQa71gM18YgC4WcTzpGiDHJVmMrKypBpg6mMVYgmuuBtLlxZWemVNgmEqVxwSqUSq9XK4cOHUSgUEY+QiIjXvKHT6ZSaLCoqKqacnlGpVJKwcUlvL2sGHkfA6iWHNqoZ5eDBgyiVSsbHx6n4UgUHvncgoqH4qSC1MJWyK8vIW5OHKccEbbDn63umvF6FRsHSR5cy/MYwrm3ve4uYyyTgfDF3w1xqvlkTdL2bNi32IkUAp9PN66838/rrvmquHqSmannssU+xYUOl9JrL5aKlpYWxsTHmzZvn57Qg2kCJDzNdXV0SUcrT49ESpVh6aG9vDzjeESsiIUq9Xs/Y2BiFhYU4HA6viNJ32VgQLgoNZODsKzogdn8nJiZ6NfRotVq/JqFYH77/3XGOFPF82ZqbmxkYGKC0tDSkKowItVrN+HjsdSdBEGhtbaWnpycic+GpwuFw0NXVxfDwMNXV1WHJV76fbreb5ORk9u3bR3JysjRzmJqaGlWqSH7DKioqory8PO7HrNy1Cxx2lAigcLMgd5D8z13LktuW0NDUgFarpaCgAGOakeKbihl6O4lnTbdS7GilS1PGqLaE5Il9NQ+Yw25vRtkM8qrzqP97vdfrizYtYtk3lpFWkUZTUxMOh4PqqmrU56k5+J8HsRuCy6jlLc2j5q4a3vr2W5Loty8Wb1qMcFig8+VOzo8i2k2fE7pp4vLL57B58yKefPJEmCP3YP78bLZu3UBFxeT3SUxZFhQUBG2kkdtA+Y4uGAwGL6LU6XRe4zvBiNJsNlNXV0diYiLLly+Pm5tEMIgkJ17LJpNJetAMFFHKZyljcRCJZRwjmOiA2DQlN3AWa6dDQ0MoFAopbTrVa3Tbtm38+Mc/pra2lgMHDrB8+XLpbz//+c95/PHHUalUPPTQQ1x++eVT2lY8oIjyaeD/5KNDU1MTbreboqKiiL+kIyMj9Pf3M2/evKi2JQgCfX19nDhxgvLycmbPnh31F/2DDz7g/PMjczkXu2a7urqYOXOmdOFGsp++dUNBEBgfH5duWkajEUEQpFGK9PT0oC4DY2NjNDQ0kJ6eTmlp6bTdsJT796O76iqw20GrxfjSS9RnZmKz2aisrPRLhdvtdgwGA3q9HoPBgNVqlQatd2/azfDx0NZXwXDNC9egXaCVxNLlggNvffctjvz+iNf753/OU9vb+bWduF1uClYVMGvlLD6474Ow2xLtrpQT0W6wSLHwskIue/QyMvIywn7nmppGeeGFWp57rpaGILOUn/vcfB566DKSJ5RubDYb9fWeh4PKysq41A19idJoNErnSN5s1d/fz8DAQFDd3+mC2ANQWFhIQUFBUALxJUmxFCG+PxIHEaPRSFdXV9T3nEggCAJdXV0YDAYSEhJ44IEH2Lt3LxaLhS9/+cssW7aMZcuWUVJSEjVJ1tbWolQque2223jggQckUjx16hQ33ngjBw4coKenh0svvVQSQpkmRLTj5yJFYM6cOVHXB2MxGpbP/yUlJTFnzpyYn8LCWesIgsDQ0BCNjY3k5OSwcuVKqc4Sbr3B6oaizU5qaqqkW+h2u6UUWGdnJyaTSXpfeno6Op2O7u5uacwjlvnKaOBeuRLbjh0o3nuP3qoqWjQaSie8GQN9XlqtluzsbK+OW6vVSk9tT8yECLDjazu46s2rqKmp8brIzYNmnBb/71r+inz23r8Xm96Tkm97q422t9oi2lYoCbiEjATm3zif0s+Wop6lZmBsgOZOTxpUPJdiXU9+My4vz2DLlvP5z/9czbFjA2zbVsu2bXV0dxtJTtbws59dwi23LJEelrq6uuju7qasrMxrFm+qkM/4BYooh4aG+Oijj6SIaGhoCJvNRlpaGgkJCdOWfXE4HJIdXCSdtPLoUIQvSYZzEJlONRtxGzNmzKCgoIDf/e539PX1ccstt7B69WoOHz7M008/ze9//3tp3jhSBCPx7du387nPfQ6dTsecOXMoLy/nwIEDrF69Oh6HFDPOkWKMiKbRRmwocbvdUuPO/v37Aw76RrrtUMsajUbq6+vRarUsW7ZMumBdLlfIcQ7xwhTJMJKoWalUBmx8GRsbo729Hb1ej0ajQafT0dXVJUWUkah3xIrB8nKagOzsbGpKSqK+kSQkJJBXkIc6QY3TGlszlWPEQWlpKUqV5zN02V28d897HHv8WEBSTCtMk+TYYkEgCbjyT5dz9V+uRpPoH5XLu3i7urokyyTfcRelUkl1dR7V1Xn85CcX09Y2RmZmIjNmeL5TRqORuro6SfN3utVWwHMD12g0jI2NYTabqampISUlxatG2dPT4xX1i8cUD6IcHBykqalJGh2KdX3RqvNYrVbJQSTeogPgSc/Kyd1oNJKbm8vVV1/N1VdfHddtgcdOatWqVdL/CwsL6e72lzw83ThHijEiElIUm2jGxsaorKz0GtIVl4+FFEX7KN9l5Y1Cc+fO9Ru6PR3zhmKE2tLSwqxZs6iurkapVOJwOKQU2MDAAGazGa1W66WJOtXuQ1GHFoi4mzUYknOTufQ3l7Lz6zvDvjdQ5+eSO5bgxo0Sz43r9a++Tm2AEQYRpj4Ty+9Yziv7X4l5n+VILUhlw183oFQHvnGG6uI1Go1eUb9caaikJF2KWsRGmrlz5/o10kwnxOH/goICKioqpO9MINUYeepVJErxexctUdrtdurqPG4gomB5vBGIKF0uFz09PXR1dVFZWRk3YXRfBHLIiFRp6NJLLw0oKnLfffexYcOGmPfp48A5UiS2QnIoUnS73XR0dEiD6IGaaOI50uF2u2lra6O3tzekfVQgX8R4zhsajUYaGhpITEz0u2loNBo/9Q75Dau7uxur1Sq16YtEGcmNR7TPGhoaory83GsbsUC5fz/KXbtIdxb7/a386nL6PuzDNOEEIdbz5J2fuVs2UHB9AUeOHPGIxp+0hCREgDfvepNPPvRJyj5VRnOQLk85ZpTNIPcTuTT+uTGgjdXF910clBCDQd7FK0I+7tLR0YHJZMLlcuFwOMjKyqKqqmraU+Ii7HY7DQ0NuFyuiIf/dTodORPpcxHy711vby8Wi0UiSjGqlGcyxD6Atra2uHa0RgKr1cqpU6dITEz00oYNJYwOsRFlIIPhSOuzb731VjSHBZy5dlLnSDFGBCIPufpNKPNdmJoyjUio4vaam5tD+jf67nO85w1FbVaLxUJlZWXEUYPvDUts0zcYDIyOjtLe3o7dbicpKUkiydTUVKlJRy4Ynp+fT01NzZRTSvJGnXJUFPIlKfpb8e0VXPzTiwH4c82fGT417DfneNUXUki7Zx0Apl4TB35zgA8f+TD8hgV44843qLmrho5dHTjGg9erU4pTWPjYQgpmFlCxuoIdN+/w+vvK/1jJvOvj04whJ0qr1UpDQwOCIJCfn4/VapWIUqlUepFKcnJy3NJ78R6zCEWURqORvr4+iSiTkpLQ6/WSf+jpmtkTreE6OzsDjmqFS73G4iDi29063YP769ev5/Of/zzf+c536OnpobGxkRUrVkzb9iLFOVJk6i3HMKmrGEj9JhCmGikaDAbq6+tJSkpi+fLlUQ0nR1s3DAaxs1XssIzUVioY5G36oiOHKB5uMBgYHByUBpJ1Oh1ms5mkpCSWLFkypVSpHMpdu8BuR+FyocDtNd7Q9nYb2pOHyeg5SuKpdKBImnNE4Uah1aK64hP0He7j6B+PUvt8bdQ1yYO/PUjSzKSQpFj1zSp0vzvMwOsP0GguAJ9aYrxrtfJGmvLycj8ZQPDcUMWIsq2tjfHxcclTUExTJicnR71vFouF2traaR+zCPSA1t7eTldXFzNmzMDlcvHhhx+i0Wik4wkm2D1ViMcsXtuRllim6iDi26eg1+unJNQu4qWXXuLOO+9kcHCQq666iurqav75z3+yYMECNm7cyPz581Gr1TzyyCOnpSYdDudGMvBcAHZ78LmxYPjggw+orq6WUjrRRElNTU2kpqZGZMckh9Vq5dChQyiVShYtWhTx9sQL4vjx45jNZpKTk6WaUiCh7XAQZ9Fyc3NjGiuJFU6nk6amJsbGxsjOzpYElMXREDGiDDYaEg5SpGiz43QrvMYbAqVKuyiKWE0mXgi2HyI2PLOByvWVIdYQOeSNNKWlpVGdZ6fTKUVfBoMhKqIUSxB9fX1UVVWdVv1Ns9lMbW0tKSkplJeXex2zOMIjHpPZbJaIUjyuWIlSlF0Ua4eRzhJHi0BEabfbOX78OEuXLkWlUqFUKnn44YcpLCzkpptumpb9+BhwbiRjOuFwOLDZbBw9ejSgiHY4iB2kkUK02hEdI3JyciIiRN+64aJFi6R5Q71eLwlti2MUYj0v2M1KlKNTq9WnVdRZnkIrLi6mqqrKT3MykCaqvKEikkhFHOmwPP86L/5+xItsgknCRWv+O1UE24/sBdms/O7KuBCiKGih1+tjbqRRq9V+4tbiQ4zBYKClpUVygZCfJ6fTKblZBPLQnC6I0WF/f3/QecdAIzyiBZS8iUz0ShSPKRxRyqPD6e7ilUeU8nppVVUVWq1Wutb++te/cuutt07bfpypOBcpTiBSyTb5MLwgCCxfvjwmYujs7EQQBIqL/Zs55BBrC21tbRQWFlJUVERnZycqlYrCwsKQy0VTNxQbKsQh9vHxcdRqtddF3dPTg8FgiNpUearQ6/WSSHo0g//ylJ5er8dsNnsdk29DhS+Ofvl3WP++k1bXbCkijGRIfrrhux//uPC7pN12OcWriuMy7iKmqcMNpMcLIlGOjY3R29uLzWYjOTmZGTNmREwqU4XJZKK2tpaMjAzPKM0UiVhOlPKI0pcoASk6PN0Rsd1up7a2FrVaTWVlJRqNBkEQePHFF/nVr37Ffffdx9VXXz3t5/80IqIDOUeKE7Db7SG1/gRBkAxxc3NzKSkp4eTJk1RUVMTUfSd2vZWWlgZ9z+joKPX19aSnp1NWViYV+bu6unC5XMyePTvovsrrhrF+qR0OB3q9nq6uLkZGRlCr1V4pSrkt0HTAbrfT2NiI1WqlqqoqoDB7tJCPhuj1esm30bdFf/ytt8jcuBGlw4Gg1dLz66cYzZuLct9+xp/7Byc6M70IMWtuFk6bE32rfsr7GAkK6WTtZQpyv7MR1+oV0jGJN+BYxg6sViv19fUolUoqKyun9dz6Qj5mUVhYKKVexR/RgDfSB5pI4Xa7aW1tZXh4OKBGazwh/+4ZjUZMJhM2m43ExESKiopCZmjiDfFeVl5eLtVRh4eH+e53v4tKpeLhhx+echf3GYhz6dN4Qa/XS00t8mH4qTTLhFpWPuwfyDEjmBh5vEcsTCYTLS0tZGRksHDhQtRqtdQdKjeaTUpK8hqjiGX2Ug632y01dsSjgUeOUKMher2ejo4OVAcOUPnMMyjtdhSCAA4HuYOnyLzpKvhUGcJ/3cje6j9B86RFVeWGStb8aA1up5tfp/86LvvqC22qFrvRU/vuoojW89cw86LVKMHvmMTal3zsIJiGqCAIUsNUsEaa6YI4ZuF0Or3S8YHOkzz66u/vx2Kx+NXzoiFKvV5PXV0deXl5LF++fNrTtOIxZWZm0tXVhdlsZtGiRSgUCoxGo186Wd7JG6/vv8PhkNxDxG5aQRDYuXMn9957Lz/4wQ/YuHHj/6XoMGqcixQn4HA4/NzlxWFwh8NBVVWV31NkXV0dOTk5MT1RjY6O0tvby/z586XXnE4nzc3NjIyM+A37yzE4OMjo6CiVlZ7aUbzJ0GKx0NjYiNvtprKyMqSThrw7VK/XYzQacblcXnqovhJioTAyMkJjYyPZ2dmUxKBGEytEIja9+SbLvvc9FHY7uN0IKHCh5tWVd5Gx5UryKvNIS0uj5dkW/vX//iUtr0vXkV+Tj6nHxNCpoai3r9QpKd5YzNq71/LqZ1+NaB2ZlZnccuSWiLchPtCIPzabTXrQSUtLo7KyMm5dvOEgr2WVlpZG3XAmQk7+YkQZLkoW66UGg4F58+adtjlL8FxbopWWbxOPCHndVYz849HJK8o+zpkzR5JqMxgMbNmyhdHRUR577LGoJdz+zXAuUowVcnKSpxd8EYv+qQh5o43Y8t7R0UFxcXFAM1Q5xBnHeM8bxjIEr1AoSE5OJjk5WVLhd7vdUtNLpI08IhELgjBlNZpoIYqVZ2Zmsmh4GIXDgcLtxo2CFkp5j0vo2p+C6nN7qPxCJUWfL8KxwIE6RY3T5In2bXpbxFqlciiUCnIvz+XCey6kZHEJhx46FDGpjjaPhtXAlUOuISp+x8fGxiguLsZut3Pq1CmvyF/8ifcYhNhUkpCQMOUxi2CNL4GiZPHhrL+/n8LCQi81nOmGfKxl7ty5IWvyGo0m4gYl3xplsPnDhoYG7HY7y5YtQ6fTIQgCu3btYsuWLXzrW99i8+bNp62h6UzHOVKUQYwWOjs7KS4upqKiIuQXRZRbiwXisuLTW1ZWFitXrowo9SimXuM1byh/ai8sLJzyELy861OEvJGntbVVauRJSUnBbrdLHpbxTN2JyjTuiy7CvXKl3+uOtDRGGhsZXbKEhdddR1JSEkP5C8l1K1Hixo3KQ4ji2IWtjbY/t/HuC2Ws+NYKFty4gGN/PBbz/mWuymTxdxaz9PKlqNVqLMMW9v5yb8TLl19VHtN2Re3OoqIivwcw0VbIYDAwPDxMa2srDodDsgzzFVCIBmKTWm9v77Q2lQQiStE43GKxkJCQQHd3N0NDQxFZUk0V4ohHampqzJ2l4YhSvKbkRJmamirV5YuLi8nPz0ehUGA2m/nRj35EY2Mjr7zyStDehLMV59KnE+jr66O2tpacnBzmzJkTETl1d3fjcDgoKSmJenujo6McPnyYrKyssClKOcSb1sGDB6Ua3owZM2KaNYTJzs7U1FRKS0tPq2JHb28vLS0tUmehKOAcj0YeXwsp244duFeunHzdagVBAKUSdDosr7zKC7/qpvWNVr+5w3BzgRBY/zQY5tw4h/RV6ay6YZVXSv7t/3ibw787LP1fk6Qhf0U+2hQt2hQtmhQN2mTPv1mVWZRfXR6ViHisjTTyFLn4Izq6y0kl1PfPaDRSW1tLVlYWc+bMOa1RSTABb7mAuNw2LF5EKa/VhosO4wWxQUmv19PT0yN18jY0NDA6OkpWVhaPPPLI/2/vzMObKvP2fydtupfuG12h6Qpl6QoOuKGDIoMO+MO6oaO8oNKBkRHB4UVgEAQcFLAOCMiwzKuouPCqiLiADC+lpSxi7U5buqV7m7Zp0iQn5/dHfY4nadImaXIK9PlcF9clsTQnbXK+z3e7bzz77LNYsmTJSMsOafnUEhwcHPSGaMzB0dERSqXSouchIuFyuRxOTk6YNGmSWf+O3zeUSCSYOnUqenp6jO4akqAyUN+BiIf39vYiISHBJpOd5qJQKFBSUgJnZ2ekpaVxN2e+JZDhIA8/UJpzYOEr07BqdV/GmJEB9cmTcO7thYhlwQIQ6XRg1Wrk3Ps6KjEdQH/HCVN7gQRzgibBdawrAhcEws3NDTKZDAqFAqNGjYKyVokre6/ofe3UV6YiY3mG0e9jCfybc0xMjMU9cFMlcuKt2dDQwPWg+eLhJOAT4fDExERB32fEKo1lWaMSbYYC4qZMjokeryUHtZ6eHhQWFsLLy0sw9xCg755ESsQREREICwsDwzDQarV45513UFBQABcXFxw5cgTFxcVYsmSJ3lwDhQZFDl9fX4tLoZZMn/JFwseMGYP4+Hjk5AxeKhuob0j88AgMw3Af6IqKCigUin5OFE5OTpxSiKH5rb3RarXcDdLYrqMp7zySpZAxcjLIw9dD7edVd/vtgJMT2F8zRdWUKSguKIBLVBSSyOM6HViRGFpGjCpEmbxuIuXG/roXaPi1gwVNPsoKJdiTLEThIqgCVVD5qFDvXo/Lf7sMnfa3Qa9REaOQ8kKKBT9d43R2dqK4uBi+vr42vTmLxWKj3pqkl1xfX4/29naoVCp4enoiJCSEO9TZOzuxVsB7MO9GuVyOmpqaAQMlOYDIZDKjTjX2RKfTcQfuCRMmcNWnwsJC/O1vf8PDDz+Mo0ePwtHREd3d3bhy5Yqg13ezQMunv0JOU5ZAPiTjx483+TVEtJpIoo0ZM4a7MZ07dw633XbbgP92qPuGarUacrkccrkcLS0tUCgUcHV1RXBwMBdU7KUnSTBUoxk9evSQAjH/5iuXy/WMjfmDPA55eb+ZDQcFcYcAh7w8iM+cATPKB7mvfIVrvYMr0gxUHuUv07NwQOWYO6Cc8//w9aEOqNpVg34fkVgEVqf/0Rq/bjy8p3nD1dVV7+ZrbnmbDNJ0dXUhPj5e8AyNrFnExsb2288jAZU/SWmrQKlSqVBcXN4tVtoAACAASURBVAyJRMItpNsavnA9eU29vb1wcnJCT0+P4JO8QN/hp6ioCMHBwYiIiIBIJIJGo8Fbb72FEydOYM+ePZgwYYJg13ODQpf3LUGn01k8SapQKFBWVmayBEp0I11cXBATE9OvNGsqKNpj37C0tBTOzs6Ijo6GTqfjAgrpD1m7QjEYRLjc09MT0dHRdgvA/CyZqNcAfWViLy8vSKXSfuVkRsPgLd+3+gUkS/BL8IP3GG9ESuqRqMyH7+ljAMOAlUjw3UOrkHPkt681t8wakh6Cx394HEDfTZ78njo7O6HVavtNhxqWk0lGbYsDiCUYrlmY2i/l21F1dnaiu7t7yCsHfN1Qa0rEQ4HIw8lkMgQHB3N9PbVabfWhxlyI+EBbWxsSExO59ZLi4mJkZWVhxowZWLNmjWCzAjc4tKdob0ytZPT29qKsrAw9PT2Ii4sbsETBH6m3dTDUaDRcthAbG6t3HW5ubtxOEukPGfYnyQfZy8vLYpktvp2UED1LBwcH+Pj4wMfHh/v59/b2Ijw8nOufktF8kiHrWnRDCoizD8xG9KxoHHvsGH44roQabbgLmr4yKsMisr4SV/wToWzp6zubW2a9a/Nd3M+auIaQ3xXRrTVWTnZzc0NraytcXFwEtTkC+tYsiouL4ezsPOiahTHfRr6CjSlNVFPvQb6At5D9O6DvYFxYWAgfHx9kZGToHSb5k7yGVmj8vqu1v6fu7m4UFhYiICAAKSkpnPnz7t27ceTIEezateuGsGK62aBBcQgY9hQZhuFOjNHR0QgKChowkIjFYq7HYst9Q51Ox52ao6Ki+olnG7sOU/1JuVyOa9euGe1PGhtKsqcazWDwnzs6Orpfv5TfH2ooa8DJP5y0+rlmvjMTsQ/F4vNHPuf2E/V7j2KcPeMIJX4bxOJspqAD6+AIRWIaPFo8OMNioG+4JjTDtNGqSCSCh4cHPDw8MHr0aAC/LaPX1tbCw8MDPT09uHz5cr9ysj16eSzLorq6GjKZbEjODqbEw0mgJELb/MV8T09PNDU1oaGhQbDpTgJfPDwhIcGoxZJIJIKbm5veAdTYyotWq+2XUQ50qCDP3dTUpCdNd/36dSxZsgSTJk3C2bNnBS3f3krQ8umvDMU+aurUqWhoaEBFRQVGjx6NyMhIs25A+fn5nHyaLYIh0KdfWF5ezo2+2/LUTPqTJFj29vbqSbwxDIOKigq7PPdg8BfwzXnuks9L8L+P/69Vz+Xs7Yw/1/4ZXy/+Gr/8zy96/4/fMwTQr39INEtHv/IYtzupUWrQVdMFZ29nuAdapq7CH6Thv+7BSpTWZP+GkPaAj4+PYL9vspjf0tICmUwGsVhsVI/XngcxfnZoC/Fww5WXrq6ufmVyshtKplqJlRc5WB88eBB79+7Fjh07cMcdd9jold5y0J6iJVgbFP/zn//AyckJHh4eiImJMbsUwrIsfvrpJ7i4uMDX13fIAy89PT0oKyuDSCRCTEyMIKdEcuptaWlBTU0N1Go1t+c10GSoLSGlUiLFZ+6+p6JJgT0Jeyw2AQYAt0A3JD6aiPwd+XqPixxEYJm+j8hA/cN7d9yLSQvNW8UxhTWDNIYi23zXEH72P1hAIYefjo4Owdd5jAl4G5Ovs2aNYjDMyQ5thbHdUKVSCYZhEBISAoVCgaioKKjVamRlZSEyMhJvvPGGXQXNbwFoULQUc+2jgN90UZuamjBlyhSz34z8vqFarUZbWxv3pucb5Zpb8tJqtaisrER7ezukUqndjEmNwTAMqqur0djYyAlJG06G8n0NyY3XFjZAg5VKzSH/7XycWnVq8C80g1GRo/D4D49DJBKhrbQNLtlvIfyrvRCzfXJxp3A3zmI6ApICMP/L+XDzNy94G8OWgzTGtEMNDzb8gNLW1obS0lKMHj0a4eHhggpH8wW8IyIiTH42DKdDydDLUHp5xFqKZORCLr3z9VJDQkLQ3d2Nd999F19++SUaGhowceJEzJ49G2lpacjIyBDM4/QmhAZFSzEnKJIg1NzcjJiYGJSXl/drsBvDHJ1SnU7HSaHJ5XI9T0MSKPnOBmTNITw8XBDfO/5raWlpwbVr17gR8IFePz9DIZOhTk5O3E2XvC5zsbRUagodo8PumN1QNCqs+vcEZ19npLyTAnGAmBvk8S8rQ/ATT3CKOt1HP4cybhLcg613PFAqlSgpKeH87+w1SGMq8+rt7YVYLLa7xZIhthDw5vfyyB9z5Ot0Oh2uX7+O5uZmxMfH2zU7NHbN9fX1qKmp0ZPFa25uxvLly+Hm5oY333wTTU1NyM/Px8WLF/Hyyy8P6LM6wqFB0VIG8lQkI98kCIWFhUEsFuPChQuYOHHigDeooewbEk9D0stTqVTD5mwA9PVTSktLIZFIEBMTY3VZir8QTW68g1lQWVsqHYj2a+3YN2Gf1f9e4i5B5olMBCcHc4M85DWx587BIz8fyowMOEybZpEiDx++ZuhQhlmsgaxZVFRUwN/fHyKRqN9qCFnjGaplmDFIZkp8Fm158BtMvk4ikaC2tnZYpOlUKhUnmh4TE8PNHXz11VfYsGED1q5di3nz5o1oiycroEHRUkwFRfLB9Pb27rdrd/nyZZM3aFuvWKhUKm7VwN/fnwssDMPoDVB4eHjY/APML9MaU6MZKqZuUGQqVqVSce4dtlbhubTrEr5/6XuL/51YIsa8T+YhakaUya8hr4u/a0jk0MzpuxIvTz8/P0GttAD9NYuYmBi99z1/NcTU67JWjxf4zdmht7cX8fHxgh38SPm/srISHR0dcHJy4iZ++RmlvX4P/F1Pvn1cR0cHVq5cCYVCgV27dllttTUQKpUKt99+O3p7e6HVavHwww9j/fr1qKysRGZmJlpbW5GSkoLDhw/frHuPNChaiqGnYk9PD0pKSgAAsbGxRss2P//8MyIjI/XKKrYOhmTVo6mpieuf8SEfZJJRkklDftnVnAEKY/A/pEKXaXU6HSca7ujoyLmB2Lo/yepYfDjrQ9T8p8bsfxMwPgB3bb4LkXdZ7jAwkCIPv49HVmHi4+MF9fzja6VakpnyXxeZouQHFHP75EQBKjIyknN2EAqy+0f8PMl0p+HrAtDvADDUg6harUZRURFXHpdIJGBZFqdOncLf/vY3vPTSS3jiiSfslrGSg46Hhwc0Gg2mTZuGHTt24M0338TcuXORmZmJ5557DhMnTsTzzz9vl2uwMzQoWgoJimTpnWRFA6ljFBUVISgoCL6+vjb3N2RZFk1NTaisrERISAjCw8PN/kCQHS+SoRDLHBIkzZl27ezsRGlpKTw8POyqRmMMYnmjVqv1MnFj/cmBBkPMpaOyAwemHICmu78Yg2eoJ4JTghGc3PcnaHIQXH1tm7nw90KbmprQ1dXFiVVbMhk6VGy9ZjHQaojhwYYv4B0XF2eTiVFz0el0nJeoOT1ThmH0AqWh1KCl8nVkeIrv39rd3Y01a9agqqoK+/btQ3j4wFKEtqSnpwfTpk3Drl278MADD6ChoQGOjo7IycnBunXr8M033wh2LTaEBkVLUavVuH79OmpqahAZGWlWVlRWVgYvLy8EBAQMWaeUT1dXF0pLS+Hq6gqpVDrkcgV/Io9klKbk3YiTR09PD2JjYwUdqiBTpfX19WYLlvP7eHK53GpnjYbLDTiz4Qy65F0ImhSEuLvjEJISAo9gYVYOSLnSyckJMTExAKB3AFCpVHqrBl5eXjYrYzEMw5XH4+Pj7fo7N1wNUSgU3EoUmWoV4gBAILZWAQEBZu8YG8Ma+TqNRoOSkhLodDrEx8dzv89z585hxYoVWLRoERYvXixYP5NhGKSkpKC8vBxLlizBihUrMGXKFJSXlwMAampqcP/996OgoECQ67ExNChaSlFRERiGMdtPEQAqKyvh4ODAWeoMNSCSgKRQKBAXF2fXmxNf3o2sT2i1Wmi1Wm7s3RblSXOx1VQpvz9JgqVOp9PLTgzLXb29vSgtLQXDMIiLixN0eIk4qDQ0NAxYrjQc5CGrBoa+hpYOvAznmgUR8BaLxQgMDOT6r7b01jQFPzu0l62VsQMAmSgXiURobm7G2LFjufuHSqXCa6+9hkuXLmHfvn2QSq0zkh4qHR0d+OMf/4gNGzbg6aefHlFBkcq88YiJiQHDMGZ9LSmTenl5oaysDNXV1XB3d+fKk5Y24/kZErGWsvfNiS/v5u7ujrKyMgQGBsLX1xddXV0oLy/nypP89QlbN9n5pdJx48YNuX9myv+PnOKrq6v1TvEMw6CjowMxMTFm2wzZCv4gTXp6+oAZAd/aiAxa8Ad5mpqaUF5ebvYgj0ajQWlpKdRqNSZOnCjoQWAwAW9Dyya+t6a5cmgDwc8OU1NT7ZaJGZOvUyqVKCoqgkqlgpubG3766SfMnz8f0dHRuHTpEubPn4/vvvvOLtO85uLt7Y277roLOTk56OjogFarhaOjI2prazmrsFsVminyMMc+it8z5JdKSZOan3URUe3BhkJaWlo4a6nIyEhBJwzJRCvDMIiNjTU6Rct3apDL5dBqtdwBYCjTeNaUSm1Ja2srSkpKIJFIIJFI+i2u2+MAQCB9a3sM0vAHrwz7XeTA1t3djaqqKowZM2ZQjV5bQwS83d3dIZVKzb7583cN+ZUNw0x5oPci31VCaDUe4LesPCIighsi0mg0eO2113DhwgUkJSWhuroaVVVVyMjIwJ49ewS7tubmZkgkEnh7e0OpVOL3v/89Vq5ciYMHD2LevHncoM2ECRPwwgsvCHZdNoSWTy1lsKBo6b6hqaEQcsN1dHREZWUlHB0djVpL2ROylMxXozEXUwcActM1R1fTVqVSa9BoNFwWHBcXp3dj5PddDcuTttjHI8NTFRUVgk5XarVadHV1obW1FXV1dXrZpFCDPHzxcFsJeJtaDeH38UimTLLDwMDAQQUnbA3DMCgvL4dCoUBiYiL3WS8sLERWVhbuu+8+rF69Wi/z7ezsFFQs4OrVq3jqqafAMAx0Oh3mz5+PV199FRUVFcjMzERbWxsmT56Mf//734IOQdkQGhQtxZSnoi1XLFQqFdra2lBTUwOFQqEXJO21Y2hIc3Oz2Wo05sKfniS9E6Jaw8+6+KVSU2su9oK/XhIVFYXg4OBBf4/G9gytkeMD+g/SCLnrZWzNgi/xZu9BHiKTJoR4uLGVl97eXrAsi/DwcAQEBNjNNcQYHR0dKC4uRlhYGDe8xzAMsrOz8cknn+Ddd99FSkqKINcywqFB0VIMg6Kt9w1JH4VMt4aEhHAn3Y6ODm4HylY7hoaQvUuJRAKpVGr3zJQ/FNLR0YGenh4wDIOAgACEhoYOWuqyJd3d3SgpKYGbmxukUumQ1kv4/Ulje6GjRo2Cq6sr9zszd5DGXpA1C+KsYOpnbo9BHktXHWwNcaQPCAiAj48PJ6OoUCjM9mu0Fp1Oh2vXrkEulyMxMZFrTVRUVGDJkiXIyMjA3//+d7t9DmtqarBgwQI0NjZCJBJh0aJFWLZsGdra2vDII49wh8OPPvqIk5C7xaFB0VJIUDTVNxwK7e3tKCsr407KA91Y+DuG5ATv6uqqV+qy5MZkbzWawSClUh8fHwQFBXH9rq6uLrAsq1d2tdRxfTCIo0N7e/ughs9Dge/9J5fLuf6ks7Mz2tvbERQUhOjoaMFLdkNdszClyGOsPGkIEfAmvXIhX7tOp+N+73xHej4ajYYLkmSXVyKR6AVKaw+kJBiTaoxIJIJOp8N7772HAwcOYOfOnZg+fbotXqpJZDIZZDIZkpOT0dXVhZSUFHz++ec4cOAAfH19sWrVKmzevBnt7e3YsmWLXa/lBoEGRUshJ2Vb7hsqlUqUlZVBp9OZHGQx57qUSqXesAt/xYCUXQ2vlV8uDAsLs7l25GCYUyolu118EXSJRNKv12UNpExsD93MwSD7Z11dXfDy8oJSqeSyLv6agb0yZXuuWRgO8hAnFPJ+9PDwgEwmQ1dXl9UC3kOBBGN+QDIXvmco0Rq2pKTMH+ThB+O6ujosWbIEUqkUb7zxhuA/EwB48MEHkZWVhaysLJw+fRohISGQyWS48847OeWuWxwaFC3lww8/xOeff47U1FSkp6cjKSnJ6p4KwzBc2UgqlQ6oimMNphw1SDBxcHBAZWUl3N3dER0dLWj/aqhTpeTGZCgWbu4yPnGTcHBwQGxsrKBDASzLorGxEZWVlf0GafhDIYaZsiUyaAOh0Wg4fVwhNUPJIE9DQwMaGhrg4ODQT+Dd3oM8JDskPo+2CDyGqyHGbKjIagiRiOOLAOh0Ohw5cgQ7d+7EP/7xD9x7773DIuJdVVWF22+/HQUFBYiIiEBHRwf3+nx8fLi/3+LQoGgpGo0GP/30E86fP4/c3FwUFBTA3d0dqampSEtLQ3p6+qD+dfybIslQhCobEX/GqqoqKJVKzvyYn3XZu4dnj6lSY8v4xoIJgGHt3VkzSEPkwvjrE9b0lPlTreYOEdkSrVaLsrIyKJVKJCQkwNXVlRvk4Tu82GuQZyjZoaUYVm5I6VWn02H06NFQKpUYM2YMVCoVli1bBm9vb2zfvn3Y+nbd3d244447sHr1asydOxfe3t56QdDHxwft7e3Dcm0CQ4PiUGFZFm1tbcjNzUVOTg5yc3Mhk8kglUqRlpaGtLQ0TJ48mRuq+OWXX6BUKuHh4YGxY8cKPl1IhnjITREA1w8iH2AAejdcWw0XCD1Vapgpk4zS09MTYWFh8PHx4bwn7Q1Zb2lqakJsbOyQb36mdGv5WRf/vUVUYYidl9AOBuYKeBua/8rl8n6ehpb2y0nPWC6XD0upVqFQoKioCF5eXggKCkJXVxf++c9/4uuvv0ZTUxMmT56MefPmIS0tbVCLOXug0Wgwe/ZszJw5E8uXLwcAxMXF0fLpQF9Eg6JlMAyDkpIS5OTkIC8vD5cuXYJGo4GDgwMYhsGuXbuQlJQk6FABf5BlsCEe/uqEMUFtS0/vLMuitrYWdXV1w7KAr1arUVpaCo1Gg+joaL2eEH9Aibw+W6uEdHR0oKSkBAEBAZyrgq0xnArlBxPS34uLi+vnnmJvyM+eYRjEx8dbVaa2xSAPEcsX8n3HX3FJSEjgBrja29vx0ksvQavV4q233kJDQwMuXLiA/Px8rFixAvHx8YJe41NPPQVfX19s376de3zFihXw8/PjBm3a2tqwdetWwa5rGKFB0d5oNBq8+eabeP/99/Hwww9DLBYjLy8PFRUVCA0N5bLJ1NRUeHl52fxDyzfdHUp2Rm645A9fsYYvFG7IcC7g8zNjYqdlbNDIsMxlK+9JIgCgVCoRHx9vE8NjS+jq6sIvv/wCJycnuLi4oLu7W2+Sl+i72iNQ8FsEY8eOtbm3n6lBHhIkPTw80NDQgM7OTr1VB6FQKpUoLCyEp6cnoqOj4eDgAJZl8d1332HNmjVYuXIlHnvssWE3AD579iymT5+ud0jftGkTMjIyMH/+fFRXVyMyMhIfffSR4K2GYYIGRXvDMAz279+PBQsW6J2SyV5aTk4Ozp8/jwsXLkCpVCIpKYkLlImJiVZnLfy9N2P+ikOFZdl+/ox8yTo3NzfU1NQMORhbS2dnJ0pKSuDl5YWxY8davDc3FO9JfkAYjt4dWbMgMmX8NQu+SwN/+MqWwy78Ui3x/BMCMsjT2NgImUzGDfLwh6+EUOSpr69HTU0N4uLiuDJ5V1cXVq9ejfr6euzdu9du2qDPPPMMvvzySwQGBnKC3CN459AaaFC8kejt7cXly5dx/vx5nD9/HkVFRfD29uaCZFpamlkalEQnNSgoSNDdL61WC7lcjpqaGrS1tUEikfQb4rH3DVKr1aK8vBzd3d2Ij4+3mW6lud6TPT09KC4uhouLSz8neiFob29HSUkJQkJCzB4mMbbzSl4bCZbmvA5+QDAm4G1vGIbBtWvXuDUPNze3AQd5jPVeh4JKpUJRURH3u3d0dATLsjh79ixefvllLFmyBAsXLrTr5/HMmTPw8PDAggULuKD48ssvj9SdQ2ugQfFGhkwLnj9/nutPtrS0IDY2lpt0nThxIrejV1FRgc7OTjg4OAiukwr0L5WKxWJOKJzclGxVmjRkoDUHe2DMe5JMF4aEhCA4ONhkSdke2HLNwvC1dXZ2Dirwbq2At60gMmmD7ZvaY5CHv+vLNxxXKpVYv349CgoK8N5772HMmDE2ea2DUVVVhdmzZ3NBcQQPzVgDDYo3G1qtFoWFhVzZ9erVqxCLxXBxcUFDQwP27t1rV5sbY1gyVWqqNMnXP7U0mCsUCpSUlAxbdsYfpPHz8+MmXvl9LvLa+NJutkCoNQtjAu8A4OnpyZUtExMTBS/LERHt7u5uJCYmWnUYsEQw3BC1Wo2ioiI4OjrqlYrz8/Px4osv4oknnsDSpUsF7aUbBkX+esUI2zm0BhoUb2ZYlsWHH36IjRs34s4774Svry/y8/Nx/fp1REVFcdlkcnKyXQYqyFRpbW2tyUEWc9BoNHpDPIaL+ERowBB+72w4pOlIdqZSqUwO0hAXFJJx2dJ7kvTuyA1Z6FF+uVyOwsJCODk5caVjW/cnB4KUiu2hRsTXrjU2yDNq1CgoFApUVFRAKpVyPXu1Wo3Nmzfj7Nmz2Lt3LxISEmx2TeYyUFAERtTOoTVQk+GbnaamJpw6dUrP1omIDOfk5OCLL77A+vXrodFoMHHiRK43GRcXN6TTK79Ump6ePqTvJZFI4O/vz70G/iI+McU11D9VKpW4du0aRo8eLXhmbImThjEDWVJS7ujowPXr1y32nuSvuAxH744v4D1u3Dg96yK+q0Z9fb1ef9JWfWV+dmgv42OxWMxdM4EccNrb23Ht2jXu9/btt99CrVZj9OjR2LRpE+bMmYNTp04JXrEwRVBQEGQyGVc+Fdok+1aEZoq3AD09Pbh48SKnxFNaWorAwEBOri41NRV+fn6DnrZJqbS3txdxcXGCTZWSqUni9ccwDDw8PODj48NlJkJItdljkMYS70lirzSYm4W9ICLW5gp4W9OfHAiSHfItloSkpaUFZWVl3JqJWq3GmTNnsG/fPly+fBkuLi6Ii4tDeno65s+fL+jOIcEwUxzBO4fWQMunIxUyKUh6k7m5udxOFym7jh8/nivJabVaVFRUoLW1dUilUmvhr5iQ7IivNymXy/XEtMnupK2CBl+RJi4uzu6lWmPek8TYlQwSCVkutaUqjE6n66fvyj8EjBo1qp8TCsMwKCsrQ09PDycRJyRarRalpaVQq9VISEjgDmBlZWXIysrCtGnTsG7dOjg5OaG6uhoXLlxAXFwckpKSBL3ORx99FKdPn0ZLSwuCgoKwfv16PPTQQyN159AaaFCk/IZGo8HVq1c5ubqCggK4ubkhIiICFy9exH/9139h8eLFwzJZWFJSAn9/f0RFRQ3o9Wcs4+L3uKyRrGtvb+cya6Htjcjzk0EeDw8PbpCH72NoT91a8vz2cNMgGDqhkP6kl5cXxGIxGhsbER4eLriTCfCbm0hERAQ31cwwDPbt24fDhw8jOzsbt912m6DXRLEbNChSTNPY2Ijly5ejpKQEv/vd71BUVDSgrqut4ZdqrVWEIX0gkpXwB13IH1MlUHMGaewJ//mNZUcDOWrYwnvSmIC3kBDDa4VCARcXF2g0Gj1JPnvvvfKz08TERG4quqamBi+88AISExOxZcsWQd8XJ06cwLJly8AwDBYuXIhVq1YJ9twjBBoUKab54IMP4OrqigcffFDPIZ7ouubm5uLSpUsAgOTkZK4/OVSjXFLara6uxtixYxEYGGjToMvfnZTL5Vx/ki991tTUhKqqKowZM8YswQRb09jYaNWaha28J0nvTIidT2OQ7Cw8PJxznSH9ScO9Vw8PDy5bttVuKNl75PcudTod/v3vf2PXrl3Ytm0bZsyYIbhSUWxsLL799luEhYUhLS0NH3zwARITEwW7hhEADYqUoUHk3vLz8zmBATIVmp6ebrGua1dXF4qLizFq1ChER0cLUqolPS65XI7W1la0trbCwcEBgYGB8PHxMduayRaoVCo9n0db9A0t8Z60hYD3UCDZKcmOBwvg/N8dWZ0YrD852Pe7du0a5HK5nmZqQ0MDli5disDAQLz11lt6U6lCkZOTg3Xr1uGbb74BALz++usAgFdeeUXwa7mFoUHRFqxbtw579+7ldpU2bdqEWbNmAeh747733ntwcHDAzp07MXPmzOG8VEEYSNeVZJOGuq5qtRqVlZXo7OxEfHy8nl6nUNdcVVWF5uZmbqrWUPqMlO7IzdaWAVvINQtT3pOOjo5QKBSIiIhARESE4L1TY707ayBDSvySuUQi6bc/aQiZrOX7LbIsi08//RRbt27Fpk2bMHv27GET8T569ChOnDiBffv2AQAOHz6M3NxcZGdnD8v13KLQPUVb8eKLL+Kll17Se6ywsBBHjhzBL7/8gvr6etxzzz0oLS0VfIxeaMRiMaKiohAVFYVHH30UQJ+u65UrV5CTk4Pt27dzawWpqakA+nolH330EVJTUwW/6ZBBkuDgYKSlpXHBwM/PjwtOfDcN4g1IVE9IoLRWIKG7u5vLjtPS0uz+/hCJRHB3d4e7uztCQkLQ29uLoqIi6HQ6hIeHo6urC3l5eXqL+GTlxR6/G352OGnSpCHLEzo4OMDHx0dPXYdvF1ZbW4ve3l7ukOPp6Yn29nZ0dHRg/Pjx3GRta2sr/vrXv8LBwaHfLjBlZEODopUcO3YMmZmZcHZ2xpgxYyCVSpGXl4epU6cO96UJjrOzMzIyMpCRkQGgL8jk5uZiyZIlYFkWISEhePTRR03qutoDvjzdYEvgIpEIbm5ucHNzQ0hICAB9E+OqqirOcYI/xDNQ+VGn06GyshKtra2Ij4/XW4IXAn7vNiYmpt9Nn680RBbxbe092drairKyMkRERCA+Pt5uByInJycEBARw1RxyyGlqasIvv/wCsVgMiUSCf/7zn9BqtfDw8MChQ4ewzFaWYwAAFo5JREFUZs0aPPLII8Nu8QQAoaGhqKmp4f5eW1trN7cNysDQoGgG2dnZOHToEFJTU7Ft2zb4+Pigrq4OU6ZM4b4mLCwMdXV1w3iVNw4sy2Lbtm3YuXMnfve73wHQ13U9dOgQrl69ColEgpSUFK7saouVCJZlIZPJcP369SEN0hhTPeH370hGQvp3/LUJsmYSFBQkuCIP0CdWXVRUBDc3N6SlpRkNbsaUhki23NLSgoqKCqsF3sneX29vr02yQ2toampCU1MTkpOT4enpCZ1OB61Wi507d6K8vBzOzs7YsWMHzp49i8WLFwu+c2hIWloaysrKUFlZidDQUBw5cgTvv//+sF7TSIX2FAHcc889aGho6Pf4xo0bMWXKFPj7+0MkEmHNmjWQyWTYv38/srKyMGXKFDzxxBMAgGeffRb3338/Hn74YaEv/6aEZVnI5XJcuHCBm3Ydqq6rQqFAcXEx3NzcIJVK7S7FxXeN5/cnxWIxIiIiEBAQYNXu5FCuh7jB8/3+rMUa70mSHUZGRgruNQn0vQcKCwvh4+ODsWPHQiwWg2VZ/Oc//8HKlSuxbNkyPP300xCLxeju7salS5cQERGBqKgoQa/TGMePH8df/vIXMAyDZ555BqtXrx7uS7rVoIM2toYvsWQ4HTZz5kysW7fO6vLptm3b8NJLL6G5uRn+/v5gWRbLli3D8ePH4ebmhgMHDiA5Odlmr+VGhEwHEs/JixcvQq1WY9KkSQPqupJSZUtLC+Lj44dlerCpqQnXrl1DeHg43N3d9dRqnJ2d9QKJPdRqhJKIM+U9SYQHAGDcuHGCT7byDwQJCQnce6Cnpwdr165FSUkJ3nvvPURGRgp6XZQbChoUbQER2wWAt956C7m5udyAzWOPPYa8vDzU19djxowZKCsrs+pmVFNTg4ULF6K4uBgXL16Ev78/jh8/jrfffhvHjx9Hbm4uli1bhtzcXFu/vBseU7quKSkpSE9PR3d3N44dO4bXXnttWKYqzVmz6O3t1csmNRqNnkHzUPbv+ALew9W7lMlkuHbtGtzd3cEwTL/dUHt7TyqVShQWFsLT0xPR0dHcZzAvLw/Lly/Hn/70JyxZskSw98bHH3+MdevWoaioCHl5edzAGTAyJ9ZvIOj0qS14+eWXceXKFYhEIkRFReHdd98F0Hcanj9/Prd+8M4771h9On/xxRexdetWPPjgg9xjx44dw4IFCyASiTBlyhR0dHToBeiRgpubG6ZPn47p06cD+G2A5OTJk1i7di06OzsREhKCDRs2GNV1tRcsy6Kurg61tbWQSqUDTi86OzsjMDCQczAg+5+kN9nd3a0nWWeuN2NnZyeKi4sREBAwLL1LjUaD0tJSaLVapKenc9khf7+wtrbWbt6T5L1QU1OjVy7u7e3Fpk2bkJubiw8++ABxcXFDfq2WMH78eHz66adYvHix3uMjdWL9ZoMGxUE4fPiwyf+3evXqIdf9jx07htDQUEycOFHv8bq6OoSHh3N/J4M8Iy0oGiISieDh4YGdO3di3bp1mDdvHrRaLa5evYrz58/j3Xff5XRdyQBPeno6p5xiCxQKBYqKiqxesyAL6J6enggLCwOg783Y2NjIlSWN2TIRAe+Ojg6MGzdOMDcTPmR1xdgwk1gsHvD1NTU19ZPkGzVqlEUHGZVKhaKiIri6uur9Dq5evYo///nPmDt3Ln744QfBtXwBmPRZpBPrNwc0KArAQIM8mzZtwsmTJ4fhqm5evLy8kJuby91EyRRrSkoKtwbS1taG3Nxcbtq1vr4eUqkUqampSEtLQ3JyssXZij3XLAy9GfmyZ/xpUCcnJ3R3dyM4OBjJycmCZxn87DA5Odns3uFA3pPt7e2oqqoyy3aK73cZGxvL7ZpqNBps374dX3/9Nd577z1MmDDBdi/aRtCJ9ZsDGhQF4LvvvjP6+M8//4zKykouS6ytrUVycjLy8vLo3tIgDJRViEQi+Pn5YdasWZz6EF/X9ZNPPuEyfKLrmpaWBqlUarIEKfSahUgkgqurK1xdXREcHMwtwXd1dSEsLAw9PT24cOECHBwc9HYn7bn+MFB2aA0uLi5wcXFBUFAQAH0nlPr6ej0nFOKCUlVVBYlEgtTUVC5zLikpQVZWFu666y6cOXNGENutgQ66/DYI5eaDBsVhJCkpCU1NTdzfo6KikJ+fD39/f8yZMwfZ2dnIzMxEbm4uvLy8LC6drlmzBseOHYNYLEZgYCAOHDiA0aNHj8jJVrFYjISEBCQkJOCZZ57hbsBE13XdunUoLy9HaGionq4rwzDYs2cPZs6ciaSkJMHdNAB9AW/DJXj+En5dXR23O8mXPRtqCVGj0aCkpAQMwyAlJcVuQYeUxj08PLgDIJF1q6urQ0lJCRwdHeHm5oZ169Zh3LhxqK2txRdffIFdu3YhPT3dLtdlDFMH3YEY6kH31Vdfha+vL/7yl78A6GvfBAYGYtmyZRZfC8U0dPr0BoIfFFmWRVZWFk6cOAE3Nzf861//0ptiM4fOzk6uxLdz504UFhZi9+7ddLLVBIa6ridPnkRLSwumTp2KGTNmICMjo5+uqz3hByNzBbz5S/h87VNrLadIdkjc6IWG/Ax0Oh3i4+Ph5OQElUqF/fv347PPPkNdXR1GjRqFhIQEpKenY/HixfDw8BD8Oo1x55134h//+Af3uR3qxHpVVRXmzp2LS5cuQafTISYmBnl5eXbV0r3FoNOnNxtVVVXcf4tEIrzzzjtD+n78npdCoeBuhHSy1ThE1zU0NBSffvopUlJSsHnzZtTV1fXTdSV7k2lpaXaxnyL2UpYGI2OSdXzLqYqKCigUCjg5OelJuhkGXH4wsmd2OBAkQ+b/DHQ6HT744AMcPnwY27dvx1133cWVxvPy8uwu2GAOn332Gf785z+jubkZDzzwACZNmoRvvvlmyBPrUVFR8PPzw+XLl9HY2IjJkyfTgGgHaKZ4i7N69WocOnQIXl5eOHXqFAICAjB79mysWrUK06ZNAwDMmDEDW7ZssTgTvZW5cuUKJk2a1O9xlmXR1NTEDfHk5uZy7hu20HXt7e1FcXExHBwcEBcXZ7ebPF+yTi6XQ61Wc0MuOp0OMpkM0dHRw5IdEpk4tVqNhIQELmDLZDJkZWUhPDwc27ZtE9xt5Ubgww8/xLlz59DQ0ICnnnqK65lTzIIu748EzG34v/7661CpVFi/fj0NijaGr+uam5uLn376CRKJBMnJyVygHEzXdTABb3vDsiw6OjpQWloKjUYDR0dHvd1CMuhib9k2YxZTLMvi448/xptvvonNmzfj/vvvF0w+bsWKFfjiiy/g5OSE6Oho/Otf/4K3tzeA4VnEV6vVSEpKgkajsVosZARDgyLlN6qrqzFr1iwUFBRg8eLFuPPOOznrp7i4OJw+fdqi8umNdrO4kRhM1zUtLQ0pKSmcrmt5eTna29vh6ekJqVQ6LLt1RKYuOjqaExkgQy4kmzTcLfTy8rJZJsswDMrKytDT04PExEQu025ubsby5cvh6uqKnTt36q10CMHJkydx9913w9HREStXrgQAbNmyBYWFhXj00Ue5/qCQi/jPPfccvL29sXnzZrs/1y0GDYojnbKyMsTExAAA3n77bfz44484evQovvrqK2RnZ3ODNkuXLkVeXp5F3/tGvFncyOh0OlRUVHBDPBcvXkRvby+8vLxQXV2Nbdu24Z577hH856RWq1FcXAyRSIS4uLhBe4dkt5AM8RArpqFIunV0dKC4uBhhYWEIDQ3lssOvvvoKGzZswNq1azFv3rxht3j67LPPcPToUfzP//yPzbWPzUWn0yE5ORkff/wx99mmmA0dtBnprFq1CiUlJRCLxYiMjMTu3bsBALNmzcLx48chlUq5yVZL+f3vf8/995QpU3D06FEAVLXDFGKxGFKpFFKpFE8++SSKioqwaNEihIWF4b777sP777+PtWvXcpJtZC3Ez8/PbsGADPPws8PBMNwtNCbpNpiTBoEIwMvlckyYMIFbd+no6MDKlSvR1dWF7777blj6msbYv38/HnnkEQDDs4hfWFiI2bNn449//CMNiHaEBsVbmE8++cTo47aYbOUz3DeLm5H/+7//w44dO/T2Q0lfMScnBzk5OdixYwfkcjkSExO5smtSUtKQJ0FJdigWi/WW4K3BmKQb30lDJpPpGRiTjLKnpwdFRUUIDg5GSkoKlx2ePn0ar7zyCv7617/iySefFETP1Zy+/MaNG+Ho6IjHH3/c7tdjisTERFRUVAzb848UaFCkmORmuVncjCxcuLDfYyKRCKGhoXj44Yc5X06NRoOff/4ZOTk52LNnD37++We4u7tbretKskOpVMo51dsaiUQCPz8/bl2ASNZ1dHSgqakJBQUF0Gq18PPzQ01NDRobGxEXF4f169ejsrISX331lZ7ur70ZbBH/wIED+PLLL/H9999zP2eqOHXrQoMixSTDebOg9jt9kCnW5ORkPV3XvLw8PV3X6OhoLps0puuqVqtRVFQEBweHIWeHlkIk6xiGQU1NDcLDwxEeHg6FQoFvv/0W+/fv5/Y/H3roIeTn58PJyemGKJueOHECW7duxY8//qinZjRnzhw89thjWL58Oerr61FWViaoog7FftBBG4pVnDhxAsuXL8ePP/6ol3HYymeyqKgIYrEYixcv1lMFoYM8/eHruubm5uLy5ctgWRbJyclISUlBQ0MDLly4gOzsbLtlhwPBsiyuX7+OpqYmJCQkcPuFKpUKGzduRH5+Pvbt2wcfHx/k5eXh/PnzGDNmDP70pz8Jfq2GSKVS9Pb2clnvlClTuN78xo0bsX//fjg6OmL79u24//77h/NSKYNDp08p9kOom4WhVNZwTf3dTBBd15MnT+LVV18FwzBwcXFBYGCgnq6rl5eX3Sc6FQoFCgsL4ePjg7Fjx3I9wsuXL2Pp0qXIzMzE8uXLR/ShhiIYdPqUYj/Ky8tN/j9b+Eyagg7yDI5IJIJYLMbmzZuxdetWzJo1i9N1PX/+PH744Qds2bIFPT09GD9+PCcwYEtdV5ZlUVNTA5lMhvj4eHh5eQHo65G+8cYb+P7773Hw4EGMHz/eJs9nDlQgn2IONChShg1qv2M/3NzccO7cOS7IEV3XqKgoZGZmAuiTlLty5Qqn61pcXAwvLy8uk0xPT7dK11WpVKKwsBCenp5ITU3lssDCwkJkZWXhvvvuw5kzZwTXKV2xYgU2bNgAoE8g/+9//zt2796Nr7/+GmVlZSgrK0Nubi6ef/55KpA/gqFBkTJsDIf9Dp8TJ05g2bJlYBgGCxcuxKpVq6z6Pjcqg2V9zs7OyMjIQEZGBoD+uq579uyxSNeVrJTU1NQgLi4OPj4+APrUat555x0cPXoUu3fvHjY5QSqQTzEHGhQpNxW2mvpjGAZLlizBt99+i7CwMKSlpWHOnDlITEy0w1XfHIhEIgQFBWHOnDmYM2cOgN90Xc+fP4/Dhw/jpZdeMqrrev36dZw+fRpTp05FWloalx1WVlbihRdeQHp6Os6ePWtXE2RzMBTIB/pK8vwVEFKSp0FxZEKDIuWGxF72O4S8vDxIpVKMHTsWAJCZmYljx46N6KBoDEdHR0yYMAETJkzAokWLwLIsOjs7uZWQo0eP4urVq9DpdJgzZw7GjBmD0NBQuLu7Y//+/di/fz927tyJ22+/XZDrHawkv3HjRmzcuBGvv/46srOzsX79ekGui3LzQKdPKSOSo0eP4sSJE9i3bx8A4PDhw8jNzUV2dvYwX9nNQ1tbGxYtWgRPT08sXboUBQUFnK7rtWvX8Ic//AE7d+68YUx/+dhaIJ9yU0CnTykUiv1wdnbGwoULcd999wEAJk+ejCeffBJAn0Gwr6+vIDJt5sIXyD927Bji4+MB9JXks7OzkZmZidzcXHh5edGAOIKhQZEyIqEyXUPH3d2dC4iGCO0HaQ72FMin3DrQ8inFrly4cAHPPvss8vLywDAM0tPT8eGHHwq6n2YMrVaL2NhYfP/99wgNDUVaWhref/99jBs3bkjf95lnnsGXX36JwMBAFBQUAOgrMz7yyCOoqqpCVFQUPvroI24yk0KhCIZZ5dMbp7ZBuSUhU53//d//jZdffhlPPPHEsAdEoG+AJDs7GzNnzkRCQgLmz58/5IAIAE8//TROnDih99jmzZs5ubsZM2ZQc1gbsm3bNohEIrS0tADoWwtZunQppFIpJkyYgEuXLg3zFVJuNmimSLE7arUaaWlpcHFxwblz5255Sa+qqirMnj2byxT5gxsymQx33nknSkpKhvkqb35qamqwcOFCFBcX4+LFi/D398fx48fx9ttvcwbay5Yto4v4FALNFCk3Bq2treju7kZXVxdUKtVwX47gNDY2coMbwcHBaGxsHOYrujV48cUXsXXrVj3FHVOL+BSKudCgSLE7ixcvxoYNG/D4449j5cqVw305w4pIJLK7CPdI4NixYwgNDcXEiRP1Hje1iE+hmAudPqXYlUOHDkEikeCxxx4DwzC47bbb8MMPP+Duu+8e7ksTjKCgIE42TCaTITAw0OrvVVNTgwULFqCxsREikQiLFi3CsmXLbslhnoEW8Tdt2oSTJ08Ow1VRbnVoT5FCsTGGPcUVK1bAz88Pq1atwubNm9HW1oatW7da9b1lMhlkMhmSk5PR1dWFlJQUfP755zhw4AB8fX2552hvb8eWLVts+bJuGH7++WfMmDGDM/2tra3F6NGjkZeXh7Vr19JFfIopzCvRsCxryR8KhTIAmZmZbHBwMOvo6MiGhoay+/btY1taWti7776blUql7IwZM9jW1labPd+cOXPYkydPsrGxsWx9fT3LsixbX1/PxsbG2uw5bnQiIyPZ5uZmlmVZ9ssvv2Tvu+8+VqfTsTk5OWxaWtowXx3lBsKsOEczRQrlJqWqqgq33347CgoKEBERgY6ODgB9B10fHx/u77c6UVFRyM/Ph7+/P1iWRVZWFk6cOMEt4g+XKwflhsOsTJEGRQrlJqS7uxt33HEHVq9ejblz58Lb21svCPr4+KC9vX0Yr5BCueGgKxkUyq2IRqPBvHnz8Pjjj2Pu3LkAfhvmATDkYR4KZSRDgyKFchPBsiyeffZZJCQkYPny5dzjc+bMwcGDBwEABw8exIMPPmj1c6hUKs5QeNy4cVi7di2APm/EjIwMSKVSPPLII1Cr1UN7MRTKDQgtn1IoNxFnz57F9OnTkZSUxDlQbNq0CRkZGZg/fz6qq6sRGRmJjz76CL6+vlY9B8uyUCgU8PDwgEajwbRp07Bjxw68+eabmDt3LjIzM/Hcc89h4sSJeP7552358igUe0J7ihQKZWj09PRg2rRp2LVrFx544AE0NDTA0dEROTk5WLduHb755pvhvkQKxVxoT5FCoVgHwzCYNGkSAgMDce+99yI6Ohre3t5wdOzT+6BKMZRbFUszRQqFMoIQiUTeAD4DsAbAAZZlpb8+Hg7ga5Zlh9/yhEKxITRTpFAoJmFZtgPAKQBTAXiLRCIiDRkGgKaKlFsOGhQpFIoeIpEo4NcMESKRyBXAvQCK0BccH/71y54CcGx4rpBCsR+0fEqhUPQQiUQTABwE4IC+g/NHLMv+XSQSjQVwBIAvgMsAnmBZtnf4rpRCsT00KFIoFAqF8iu0fEqhUCgUyq/QoEihUCgUyq/QoEihUCgUyq/QoEihUCgUyq/QoEihUCgUyq/QoEihUCgUyq/8f5kd2X+GmhfmAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "%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": { "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": { "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": { "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": { "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": { "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": { "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,flucts=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": { "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,flucts=gnm[1].getArray(),style={'cartoon':{}})" ] }, { "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": { "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,vecs=anm[0].getArray(),\n", " flucts=prody.calcSqFlucts(anm[0]),\n", " style={'cartoon':{}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Multiple Structures" ] }, { "cell_type": "code", "execution_count": 15, "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": { "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": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,prot,flucts=[gnm[1].getArray(),prody.calcSqFlucts(anm[0])],vecs=[anm[0].getArray(),[]])" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'1.10.8'" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.__version__" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.6" } }, "nbformat": 4, "nbformat_minor": 1 }