{ "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": [ { "name": "stderr", "output_type": "stream", "text": [ "/usr/local/lib/python2.7/dist-packages/prody/atomic/subset.py:35: FutureWarning: Conversion of the second argument of issubdtype from `int` to `np.signedinteger` is deprecated. In future, it will be treated as `np.int64 == np.dtype(int).type`.\n", " elif not np.issubdtype(indices.dtype, int):\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcUAAAE1CAYAAACWU/udAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzsvXmYHGd17/+t3mamu2dftIykGc2iZbRrtIIxXtjMBQfnZ7BZAuGSxRBuTJJ7L04guRByscx1SAwkMRCwTRLgOtyAjVmCg+0YG9larMWyLGmmp2dfemZ67+ruqq6q3x+jt6aqu6q7qnqZ7f08jx49Gk11VW/vt855zzlfRpIkUCgUCoVCAWxLfQEUCoVCoSwXqChSKBQKhXIdKooUCoVCoVyHiiKFQqFQKNehokihUCgUynWoKFIoFAqFch0qihQKhUKhXIeKIoVCoVAo16GiSKFQKBTKdRwmf5+Ov6FQKBTKSoQx8ks0UqRQKBQK5TpUFCkUCoVCuQ4VRQqFQqFQrkNFkUKhUCiU61BRpFAoFArlOlQUKRQKhUK5DhVFCoVCoVCuQ0WRQqFQKJTrUFGkUCgUCuU6VBQpFAqFQrkOFUUKhUKhUK5DRZFCoVAolOtQUaRQKBQK5TpUFCkUCoVCuQ4VRQqFQqFQrkNFkUKhUCiU61BRpFAoFArlOo6lvgAKpZJIkgRBEMAwDGw2GxjGkBk3hUJZI1BRpKwJiBhmMhmk02lIkgQAsNlssNlscDgcsNvt8r+pWFIoaxOGLA4GMfXLFMpSoxRDSZLAMAx4nofycy9JErK/BzabDXa7Xf5DxZJCWfEY+vJSUaSsSrTEkAgax3Hyz/SOJX8rvx+pVAo8z6OpqYmKJYWy8jD0JaXpU8qqQpIkZDIZCIIgC5/NZq6ejAhcttClUimEQiHU1dWB53nV/9HIkkJZHVBRpKwKiBhmMhkA0BXDZDKJkZER2O12eL1eeDweVFVVGRIv8jvZj0uiyUwmQ8WSQlnhUFGkrGi0xFBLcFiWxdDQEGKxGDZt2gRRFBEKhTA+Po50Og273Q6Px6P643K5VI/FMEzO3iP5ufJv5bUBC2LJcZzq/6lYUijLE7qnSFmREDEcGhrCli1bdMUwkUhgaGgILMti69ataG1tRSaTgSiKqt/PZDJgWRbxeByJRAKJRAIcx8HhcMgiKYoiWJbFjh07ir525R+GYcBxHEKhEDZs2EDFkkIpD3RPkbL6EEVRLqABgMnJSXR2dub8Xjweh8/nQyqVQnd3N5qbm/MKi8PhQF1dHerq6lQ/z2QySCQSiMfjiEQiiEajOHXqFBwOh5x+JX+cTqeh56Al4CRyXb9+PXie14wsHQ6HHGFSsaRQygMVRcqKIFsM9SLDWCwGn88HnufR1dWFpqamooTD4XCgvr4e9fX18Hq9mJ6exvbt28HzvBxRBgIBxONxZDIZOJ3OHLF0OIx/zbT2QUlEqVU1mx1V2u123deGQqEUhooiZVkjiqJcTQroi2EkEoHP54Moiuju7kZjY2NZr8vpdKKhoQENDQ2qn3McJ4vl9PQ04vE4BEGAy+WCx+NRCabdbjd0Lr3nLEmSfLOQTyyJYFKxpFAKQ0WRsiwxKoaZTAZnz54FwzDo7u5GfX192a5Jr9BGicvlgsvlUokyifKIWE5MTIBlWQiCgKqqKjma5DgOgiCURSyHh4fR2dlJxZJCKQAVRcqygaQJeZ6HKIoA9Bf+YDAIn88HjuOwb9++soohwYgo6h1XVVWFqqoqNDU1yT+XJAnpdBqJRAKhUAjJZBLnzp2DKIqorq6WI0qv1wu3222431LrNQsGg+jq6pLFMhsqlhTKAlQUKUuOUTGUJAnz8/MYGhpCVVUVduzYgUuXLqG2tnYpLrtoGIZBdXU1qqurUVNTg1Qqhd27d0OSJKRSKTmyDAaDYFkWoiiipqZGJZY1NTVFiSWgjiyzISJJinyoWFJWO1QUKUsGWYwzmQwuX76MHTt26Irh3NwchoaGUFNTg76+Pni9XgDWozcrVOpcDMOgpqYGNTU1aGlpkX8uSRKSyaQslnNzc2BZFgBUYunxeEoulhzHqX6fFPVQsaSsNqgoUiqOUgxJv2A4HNacFBMIBOD3++H1erFnzx643W7V71gRKqsLdyUFWO/8brcbbrcbra2t8s9FUVSJZSAQQDKZBAC43W54PB65D7Ompsbw8zcqlizLIhKJYOPGjVQsKSseKoqUiqElhnqR4czMDPx+P+rr67Fv3z7U1NRoPqbNZpNTrmsVm80mR4hKyLCBRCIBQRDg8/mQTCZhs9lksSR/qqurLYulIAhgWRYMw+SNLLXsuahYUpYbVBQpZUfPsUJLDKempjA8PIzGxkYcOHAA1dXVeR97NaZPS4XNZoPX64XX68Xo6Cj27NkDYFHEEokEIpEIJicnkUqlVOJK/hiZC5vvPSX/L4oi0um06ucMw2hO7zE7wJ1CKSVUFCllQ0sMtRY8URTBcRxOnjyJ5uZm9Pf3o6qqytA5rEaK+ayj9FhpoqiH3W5HbW1tToGSIAhyCtbMXNhCr6XeDRA5JxnIoPx9vaEEFEq5oaJIKTlmxHBiYgKjo6MQRRGHDh2Cy+Uyda7VIlQAlvx52O123VF3ZC7s/Pw8RkdHVXNhye9wHGf4/Ss0RF0plqFQCABkH0vlfiUVS0qpoaJIKRlGxVAQBIyPj2N8fBzr1q3DkSNHcObMGcOzQ5VUck+xnAK8nBd2vbmwPM+DZVlMT0+DZVm89tpr4Hm+6Lmwyr+BBR9LEm1mR5aSJKmKe+hcWEqxUFGkFA1xrIhEIpidncXWrVt1xXBsbAzj4+PYsGEDjh49Ks8FJUUaRqe5EMwKFcdx8Pv94Hle3m9zu92G5pOupqi0FDidTtTX18sR4tatWwEsiCVxG5mZmUEikShqLqwkSbozXbMjS2Uql6RhqVhSzEBFkWKZbC9DURSRSCRyFp1MJoPR0VFMTk6ivb0dx44dy1kMScRnVhSNRoocx2F4eBhzc3PYtGkT6urqwLIsJicn5erMYqbIrGWy9xSdTicaGxtz5s8qR91NTU3Jr7vL5coRS+XnQBRF3ffBaBpWeY3Uy5KSDyqKFNPoGfva7XaVQPE8j9HRUUxPT6O9vR3Hjx/XFT2radBC0RuJDOfn59HR0YGenh55ek5zc7PqOSmnyMzPz4NlWUiSJLcvVFdX5yywFONFS1bnwiaTSTl1amYurPJv5fmAhRs1nucBLDircByH1tZWKpYUKooU4xRyuSfCxnEcRkZGEAgEsHnzZhw7dqzgYlZqUcwWw97eXjna0BpnpjdFRtkYH41GwbIsTp8+Lff6kQjH6/WqKjLXEsXcJBiZCxuLxTA/P49AIFCSubDKv4GFPUvSLqIUSwKNLNcWVBQpBSkkhoRMJoNoNIrTp0+jo6MDx48fN7xY2Ww2S/t12WKaTwytoOzda2xsRDQaxcGDB+Vev3g8ntO+oEwFer1eSwVEK4lyRM7KubChUAhNTU1oamrSjOiTyaRqLix5/Y2OuiNp+3x7llQs1w5UFCm6GDX2TaVSGB4exvz8POx2uykxJJBCG7OQSLHUYqh3LoJer18mk0EikUA8Hsfs7Cz8fn9OkYmZ4p6VQLnTydnFM0bmws7OzmrOhSVD1JXXW8yeJRXL1cfq+FZSSopRL8NkMgm/349IJILOzk709PTg7NmzlsSomCb8sbExJJNJdHZ2lkUMCUaqTx0OB+rr63OsrDiOkysyJyYmkEgkVKlAl8slj79bacU95RZFI69JMXNhk8kkPB6PqedRSCx5ngfHcfJnZmpqCps2baJiuQKgokiRIWJ45coVdHV1weFwaH5hWZbF0NAQYrEYtm7dip07d8pffq39OiOYFUUSGU5NTWH9+vXYt2/fshYTl8slpwAJylRgOBwGy7I4e/YsAGimApfr4kl6Bcv5+Fafu5G5sMFgEIlEAtPT0yWZCwtAtYeeyWQQCATQ3t6uEkvlNVKxXD5QUVzjaHkZRqNRzYUokUhgaGgIiUQCXV1d2LVrl+p3ivkSGxXFdDotp2o7OzuxZcsWeDyeighiqfsUlalAr9cLlmWxd+9eVXQTi8UwPT2tGuS93Ip7lkOkaBblXNh4PI76+nq0tLQYngtr5rUne5Zaz0H5/dMSS2WPpd6+J6W0UFFco+Qz9rXb7aqILx6Pw+fzIZVKobu7G83NzSX/YhYqtMkWw23btoFhGPj9/lUx0UaJcgFua2uTf75ci3squadYDpT9saWeC0uO1au+zjdEneyVZz9/vbmwVCxLAxXFNYYRl3ubzQZBEBCLxeDz+cDzPLq6utDU1FS2L55eoU06nYbf70cwGFSJYaHjCp1rJVJMcY8gCKb6/MywEiNFs4+fby6sshI2ey6s1SxGIccRQRBUr3s8HofT6URtba1KMKlYmoeK4hrBqJchsHBne/nyZdhsNnR3d+dMJikH2enTbDHcvn275rVabeWwwnJdXLSKe8iNDynu4Xke58+fV03uIYJZ7OSeSkSKSy2KeugVVpG5sPF4HMFgEJFIBKdOnSpqLiygL5bBYBButxtVVVV5I0sqloWhorjKMSOG4XAYPp8P8XgcXV1d2Lx5c8Wuk4iiUTEk0Hmk2jAMoyrumZ6eRn9/v6q4Jx6PY25uLqd1wWxxTyUixZUWiZK5sPX19bJY9fb2lnwuLEEQBHn/UYkyssyGiKTSdYSKJRXFVYtRY19g4S7T5/PB4XCgt7cXk5OTcLvdFb1eYiPl9/uxdevWgmJIsNlsOX1iFH2MTO4xW9yzFtKnxSAIgvz4WnNhs6N65VxYMupOby6s8hxaImokDctxnOr/1rpYUlFcZRAxnJqakkvL9areiBi6XC7s2LFD3quamZmx3FpBHtvoF4hEhtPT02hpaUF/f7+pL99qiRSXesEppriH9PmVi9Uguvn2crOjeoJyLmw8Hs+ZC6uMLDOZjKn9YrNiyTAMUqmUHPWuZrGkorhKyI4Mw+EwGIaB1+vN+b25uTkMDQ2hpqYGfX19Ob9TjEehUbcLZZp069atcmRq9gtWST/FcrMcxd1IcQ8pNBkZGYHL5VJVwepFNmZYDZGildfAyFxY0mc5Pz+PWCwm3whb3S/OJ5Y+nw+dnZ3yDfPY2Bi+9rWv4etf/7rp57acoaK4wtEz9s1uq5AkCbOzsxgaGoLX68WePXt0U6TZx5qhkCim02kMDQ0hFAqp0qQTExOW0qBmI8V4PI7BwUHE43G5P3AlNMgvN5QFJizLoqWlBQ0NDao0oNbknmIW63KxXEVRD+VcWOL0cu7cOezatUtuHVHuF0uSZHkurPKcJEIlzyUcDsvTgVYTVBRXKGRIt7I0W/khJ8ImSRJmZmbg9/tRV1eHffv2oaamJu9jk5YMK2TbRxGyxXDHjh05bSBWIj6j1aeJRELutezq6oLX60UqlUI8Hs/ZQ1MuHmQPjaIPef3zpQHzFfco04BLcWNSiUi03HNuyZ6iy+XKOxc2Ho+rRt2ZKa7KZDKq5xGJRHJaVFYDVBRXGFqOFVpfaJvNhmAwiJGRETQ2NuLAgQOorq42dA673W65eCVbUFOpFPx+v64YKo+zOhA833FkJF0ikZAHD5DCBjIrM3sPjSwe8/PzGBkZAc/zcnUgz/OIRqMlSQuuJvKJipXiHmVhiSiKZd9XLOdjl6s/VEmhoeaF5sLG43FMT08jlUoBgGrUndfrRXV1dU5rTDQaRUNDQ1mf11JARXGFYNS+SRRFTE1NyXuG/f39qKqqMnUuu90ufznMQsRNKYZdXV26Yqg8zsqemp4oJpNJeT5rd3c3Wlpa5PPnO49ekzYZ6B0IBFRpwew7bbfbveZSsFYFy2hxTyqVwunTp1UN8eT1Xgm2XJUQxWJffyXKubDKmxWWZXH58mVIkoTLly9jZmamqEjx6tWruOuuu+R/Dw0N4S//8i/x4Q9/GHfddReGh4fR2dmJxx9/vCK90gQqisscM2I4MTGB0dFRtLa2oqenB6lUyrQgAsUVr5ANeTIftZAYEqxMpgFyxTSVSmFoaAiRSATd3d3o6+sriUiRtKDL5cLOnTsBaKelWJZdtjNKy0Wp+wizi3tCoRAOHz4MnudVbhfxeByZTEYu7lGmYZdTFG+k8KxYSlmkpZwLq+TUqVPYvHkzRkZG8MILL+D8+fOYm5vDj370I+zcuRMHDhzAJz/5ScPn2b59O86fPw9g4cahvb0dd9xxB06cOIFbb70V9913H06cOIETJ07ggQceKNnzKwQVxWWKJEmIxWIFZxsKgoCJiQmMjY2hra0Nhw8fhsvlwuzsLBKJhKVzWym0IZFhMBhER0cH9u7da2qhLCZ9SqrxyJ5lV1eX7Nyhd4xVlIuPXlpKGekEg0GMjY3Jo7+UQmm2QbsU11yux6+E4DudTjQ0NKhSdsq2hVIX95QKZZ9iOahE1TK58amtrcXu3bvx0EMP4cSJE9i/fz/e+c534sqVKxgbG7P8+L/85S/R3d2Njo4OPPHEE3juuecAAB/5yEdw0003UVFcyyiNfQcGBtDe3q6ZOhAEAWNjY5iYmMD69etx5MgRVSrJ4XBUxMaJiGE4HMbWrVshSRIaGxsr1lohCAJCoRDOnj2bd8+ykui1MZDKzHg8rmrQrq6uRk1NjVxmb7YycKlZylaSfG0LRop7vF5v2a+/3OnTco/BA7SHA0QiETQ0NMDj8aC/vx/9/f2WH//73/8+3v/+9wNY6JPesGEDAGD9+vWYmZmxfuEWoKK4TNAy9nU4HDlCkclkMDo6isnJSbS3t+Po0aOa0UYxKVAjkaIyTakUo0gkYjkNatZPcXh4GDMzM3A6nTh69GhFxLCYc+hNM0mlUgiHw5ibm4Pf75crA5X7Z8WkYCsRyS31jUg2hYp7SNXx1NQUWJbFmTNnLFtDFaLcomi2cd/qObLXmWg0mjPz1Qocx+HJJ5/E/fffn/N/SzEcgIriEpPP5V4pTjzPY3R0FNPT02hvb8fx48fzfhHsdru8D2mWfAKVLYbZaUq9lgwj5zRyx87zPIaHhxEIBNDR0YH9+/fD5/Mtu0XZKGTxttvtmJmZwe7duwEsfC60rIpIsYmyv7ISKdh8lFN0Sx3FZReXiKKIV155BQcOHFAV9yhT3sUW91SiD7LcnwEtUSSRYrH87Gc/w8GDB7Fu3ToAwLp16zA1NYUNGzZgampKVYRVCagoLhGiKOa1bwIWUqCpVAoDAwMIBALYvHkzjh07Zuiu0Ko4kWOzI8VCYkiw2uNYqNAmk8lgZGQE09PT2LJlC44fPw6bzSY3J682bDabbgpWWUIfj8dzxn55vV7V/tlK9jus1Ii3fClvcnOiHOBtprin3JFiJapb9SLFUoji9773PTl1CgC33347HnvsMdx333147LHH8Bu/8RtFn8MMVBQriBEvQwLHcZibm5P764gIGKUUU2kA42KodazVcyoh6eKpqSls3rw553Wo9OzTpRZgvWKTdDot71fOz8/L+2ekEIjnebkaudQiU849raWeS2qluEcplG63u+zCXilRzD5HKdKniUQCTz/9tGpU3H333Yf3ve99+Na3voWOjg48/vjjRZ3DLFQUK4AZMVTOBK2vr0dLS4slC6diRJE071++fBmRSKRgNaeSUomiIAiqvVO9CHk1zT61inLsV/b+GcuyiMViEAQBV69ezRnmXYp+v3LeKFRCUMyKrtninkQigcuXL+c0w5fqeWlFcaVG6xwkYi4Gj8eD+fl51c+am5vxy1/+sqjHLQYqimUk28sQ0BdDEpGRKs7t27cjEAggFotZOrdVsUilUipPRaNiqDyvFTEm10uqasfHx/OKIaHSkeJK2rsk/WbV1dWYmprCvn37AKiHeStTgsSmSBnpGBGMcgrXShoGrlfcc+rUKXR2duqOFFS+5laKeyoRKWbvWy51xqScUFEsA0oxvHTpEnbt2qUrhmTySjQazUlPFhPtmf1iJZNJ+P1+OTIMh8PyxrcZivE3TCaTeOmll7BhwwYcO3bM0N0vjRTNo+UWT1KCJAUbDAblPtfsFobsKGc17CmWE4ZhNCfHkJGCxOlidHRUVdyjvDnJF8lXKn2qZSCwkm4SjUJFsYRoOVZEIhHND04ikZBncnZ1dWlOXilGFI2SLYZElH0+n6XHM1vgo5zEI4qibouJHkvhp1iuhbpcz8PI9SpTgsR5AVC3MESjUUxNTSGVSsFut8sLNsdxsnt8qVlJkaIWZOtEC72RgvmKe5RCSYp7lqIlI5VKGZ6lvNKgolgC8rncZy9G8XgcQ0NDSCaT6OrqUs3kzMbhcFhuqyiEMkK1kibVw2j6VBRFTE5OYmRkBG1tbThy5Ig839Ls+cxGioIgIBgMwu12l3RvpxiWwzVooUzxKTMHyhQs2X/OTsGSKthiFuyVHilaKUIqVNxDDIdJcY8oirJAer3esgx/yBbFcDhckh7F5QgVxSLQ8zJUQloNiHURx3Ho7u5GU1NTwS97OSLFbDEs1WxQQiGRkiQJU1NTGB4eRktLizyWzipmrl0URYyNjWFsbAwNDQ2YmJiQox5letDr9eqKM4lMl6uI6VHq61WmYKempnDgwAEwDCOnYBOJBMbGxsCyLERRzDG/NWoRtdIjxVKlNvMV91y9ehUulwvJZFI1uSf7NS/mBjBbFEvVuL8coaJoASNiqPzdc+fOQZIkdHV1qT7QhSg2UiSCbLPZVGJYykHZ2eiJoiRJmJ6eht/vR3NzMw4dOlQxn0JlVLp+/XocPXpUNcQ6u/BkaGgImUxGNTuT3IEvRbq2WCo1+1QvBStJkqbrghHvynKPMFvpokvWnqamJlVkmZ32npycRCqVslzcs1a8FAEqiqYwI4bhcFg2te3u7rZUtFJspGi325FIJDA6OmpJDK1ERNmiqDQ5bmhosGRlZRVlVNra2irPhyWpKIJe4Ymy9292dla2z7l69Spqa2tXlAnxUo55UxaamPGu9Hg8effkSkG5h3UvVWO9ni2Usrhnfn7ecHFP9r5lqabZLEeoKBqgkMu9klAoBJ/PB7vdjp6eHoyMjGhWbRmhmMpKsoC/+uqr6OnpMR0ZknOb/UKTQhtJkhAIBDA0NIT6+npTJsfFojx3Y2OjpahUr/fv7Nmz2LhxI5LJpGohVxZBkL+Xy1Dv5RrZFvKuJAt3PB7HqVOnVN6Vysi9GMpt61Spdgmj57Ba3EOG1ZO9y9VqMAxQUcyLJElyWTqxb9Ja6CRJQjAYhM/ng8vlwvbt2+VxUZVsqwDUe4bV1dXYvXt3ji+aEch1m/1CMwwDlmXx8ssvo7a2Fvv370dNTY3p81tBkiTMz89jcHAQtbW1ZRFicgeevSAo2xnIXpokSap9Ha/XW5aJMqsN4l3Z1NQEt9uNcDiMrq4ulXflzMyMnIIl3pVWev2WemJOKSiF8BYq7pmensbExAR8Ph/+/M//HI2NjWhubsbjjz+O3bt3o7e311T1cTgcxu/8zu/g0qVLYBgG3/72t7F9+/YlNRcmUFHUQGnsOzg4iNbWVlW0oPy9ubk52eW+r68vR4DKWUGqJJlMqpru+/r6cOnSJcuRptkolQjSwMAA0uk0jhw5YjlCtkIwGMTg4CCqq6uxd+/egucutTApF3KCcl8nEolgYmJCHuqdXdhTblaqEBPRMuJdmZ0OVAql3uB0URTLOg2m3OlZoHwuGco9YmKuvXPnTtx222344he/CI7jMDg4iB/96EeIRqN46qmnDD/2vffei3e84x34wQ9+AI7jwLIsvvjFLy6puTCBiqICLZd7l8uVI2qSJGF2dhZDQ0PweDzYvXt3Tu6eUG5RzBZDMigAKH7Um9FjSXRWU1ODnTt3YmBgoGKCGIlE5H1TrZuSUmOm0EavnUE51Jv4KmYyGXAcB7/fb7pCsxDLNX1qhEL72nqDvJXtC5OTk5qzSb1e76rYU6xEtKvlhHPLLbfg9ttvN/14kUgEzz//PB599FEACzeULpdryc2FCVQUoTb2BdSj2JSipiwaqaurMxSRFGPhRNBaGFiWxdDQkKYYKs9dzkgxFAphcHAQLpcLu3btgtfrVY20M4uZdodYLIbBwUF5odu/f7+lc5qlFNWnWqkqnudx/vx5eL1e2QGDpAeVEY/X6zXdJL8SW0gIVhd8stBqeVcqi6fC4TDsdjvm5uZK4l2ZTSVEEShvJkDLmqqYlgy/34/W1lZ89KMfxYULF9Df34+HHnpoyc2FCWtaFPOJIcHhcIDjOLmK0WzRSLGRIonYyIfSiBgSrM4hVZ5Xi3A4jIGBATidTuzcuVMVnRVTHERaSPItIolEAoODg+A4Dj09PWhsbMSvf/1r0+dabkLBMAzsdjtaW1tz0oMk4pmdnYXf75eb5PWsolYTpXyflLNJyWs8ODiIhoYGVFVVaXopZhdPmU21VmqMXDnR81K0KoqZTAavvPIKvvrVr+Lo0aO49957ceLECdXvLIW5MGFNimI+Y9/s34tEIvIdzMGDB023EzgcDqTTacvXSsSJ4zhZDLu7u/OKofLcxdhHZR8biUQwODgIm82GHTt25KSsyHFWI6l8Fa8kTZxIJNDT06Pqg6skle5T1KoWzJ5TqrSKUkY7K6VdJB+VSA06HA5dL0VSBUvS3IIgqPpXiT2U3jVWwgC43J9HPdsoq0UwmzZtwqZNm3D06FEAwJ133okTJ04subkwYU2JohkxJPM4SW/Vjh07LJ2zFHuKr7/+utzvaEQMCcVGiiTii0ajGBwchCRJ6OnpKdskC60oM51Oy+4h3d3daG1t1Xz+ViIKK8csh+b9fHNKSVSpHDBNhrRPTU2VZPRaJVnK5n2n04nGxkbNFGy2PRQAuQqWCGZVVZVsAL0U118qSm0wvH79emzevBlXr17F9u3b8ctf/hJ9fX3o6+tbUnNhwpoRRdKMDeQXw/HxcYyNjaGtrQ2HDx9GMpnE2NiY5fNaFUWSJiVN91u2bDG9gBdbaBOPxzE+Pg5BENDT01P2viSlKJKik/n5eXR1dWHHjh1508RmSt/J+29V3JZaFPWw2WyaEU8kEoHP50Mmk8H4+DgSiQQkSZLdL0oxBqxcLLfqUD17KOJdmV1pzPM8vF4vOI4bv1lwAAAgAElEQVQriXel1vUvxXCARCJRVFHbV7/6VXzwgx8Ex3Ho6urCI488AlEUl9RcmLBmRFFvQDcA2cNvYmIC69evlyefAFBVo1rBrDAp9wy7u7vBMAzq6uosLVZ2u101ucUopNBDkiTs2rWrYr1CJKKZmJjAzMwMOjs70dvbW3DRqmT0ttxEwwgOhwNVVVUqs2rl6DXlGDAyB1YpluVO/+VjpYxhUxZEKbl69SrcbjcYhtH1rixmT3gpHDLIlKFiXrf9+/fjzJkzOT9fSnNhwpoRRSB3vyuTychiuHHjRk3bIofDYdkfkBxvRFRZlpX3zJRp0lAoZFmUzQoyGVqeSqXQ1NQEr9dbMUEkBSXnz59HZ2cnjh8/bvhLV0lPxeWQPi0FeqPXMpmMvFc5PT2NeDwu76MpxbJU7SKFqIRLRjlFRZIk1NfX5+wJkwkx2XvCZqP3pTAYJqzEG0QjrClRJGQyGYyMjGBqagqbNm3C8ePHdT9Yxe4JFjqeiCHLsppWUpXoNVReQ09PD5qamjAxMVGRoQPEuWJ8fBwMw2DPnj2mhXi1CFU5n4PRBczhcGhONlG2MgQCASSTSVlYOY5DKBSy1C5SiJUSKeqhlZ5VjhDU866MRCKq6D17cDp5nStRyJNtMFyJfcylZE2JYiaTwdDQEAKBQEExJBRTTQnoC1MhMSx0fDHnJigb/7u7u1XXYDX1Sih0h6/lXHH16lVLd59mIkWybzw6OpqTKvR6vQU/D+US4HLedRd7vVqtDMBidB8KhTA3N4fh4WHVHNhiU4Pk2leyn6KZSNSId6WyLcflcsnf8VgsVrYCquz0aSwW06w8Xy2sKVH0+XyoqakxlZorlmxRNSqGhGIiVT1RTKVSKk9Fvcb/Yto59BYzpYVUS0uLav/WahrUiFARxwy/349169bh4MGD8mxbMlkmHo9DFEVVCsvr9apSWKslKi0F5KbC5XKht7cXQP52EeWMUqMN8qshUixWqPRcXDiOw/j4OGKxGMbGxuQCKqvelXqsJdsoYI2J4o4dOyq295RNIpHA0NAQWJZFd3c3mpubDX1Qi4nYsoVN2d7Q1dWFnTt3FqzotAJpBVEuNkrnCj0LKatROWn614KM5PP5fGhoaJBNjTmOgyRJmqlCksKKxWKYmppCKpWSG7lJheFKa2soV7SV/dj52kXIaxcKhTA+Pq6aA6vcR1O+rpUYslDuaTDlEF3yOrvdbrhcLrmIiqRgiXcl+fwqvSuVg9ONkF3ME41GqSiuFqx++JVmvWZJJBJIJpO4dOmSKTEklCJ9mk6n4ff7EQwGC7Y3lOK8SkFVOld4vd6804DyiVuh82mJaTAYxMDAADwej2G3DuXgaWUBCmnkjkajmJ2dxcTEhOwor1zU15oLhlHR0qvOJK9rPB7HxMSEPKOUROssy8o3MCvxda20NZUyBavnXTk7Oyunuol3pXJwevb1ZkeK4XB41dpGAWtMFK1CUphmpoOQSs5kMgmn04nDhw9bEtVi0qdkr+Hs2bPo7OzE9u3bTTX+Fzs3NRQKYWBgwLBzRanSp9FoFNeuXYPdbpdnshY6phCkkbu2thbr169HfX293Nag54KhF/0sBZWKFM2i1yBPovVAIIDx8XH4/X5VwQn5s5TtIkYpp5hnMhlDN3uFvCv1bko8Hg/S6XROpFiuAR7LgeX/iSohVj+cZkRRKYYkMjx16pTlSNNKxMbzPIaHhxEIBMAwDI4dO2b63MVEioIg4OLFi3C5XKacK6yKIjkukUhgYGAAmUwGvb29Zf/iKtsasl0wtBYaZVRZW1tbsqHThSjnHmg5IjhltB4IBNDV1QW3260qOJmZmYHP55MnxmS3i6zm6kglxe5ZalmeKW9KEokEUqkUXnnlFUSjUTzyyCOora2V3xsro9g6OztRW1sLu90Oh8OBM2fOIBgMLgsvRWCNiaJVjERrWmKY7bRh5a7WjDjxPI+RkRHMzMygo6MDx48fx0svvWRpgbAiUMS5Ih6PY8eOHfLE+3KeE1hYGIaGhsDzfNnnohqJMLWin3zeimSREQShbH1nyzVSNPP4egUn6XRa5XxB2kW0CntWG+VoyVDelADA7Owsjhw5glQqherqanz/+9+H3+/Hhz70IQQCAfT29uJf//VfTZ3j2WefVU0EOnHixLLwUgTWmCgWGylqoRRD0uOn5bRhNQVq5FjSdzk9PY3NmzeXpLrWjBgrm/57e3sxOTlpKKWTjdlCGzIkfW5uDlu2bJEnAJUTq9WneuX2JH0Vi8WQTqdx7tw5uYLQbKWmHistUlRSKMOi7PlTLrLZ5sMjIyOae2grvZK4UtZUAFBdXY0bbrgBL774It75znfirrvuArBQjVosy8VLEVhjomgVLWEyIob5jjdKPnHKZDIYHR1VDSEoVdrISNSm7HMkERrDMJienra8N2jkOOVNQGdnJzZu3IiGhoaKFWKUYiGVJAnTFwXMD0jgk26k49WYHZXQVNsGLiFiKi4gHRfAsSL4VARgJNgcgMNlg7PaBmeVHa5qByQBqGlm0H2zCx03lLZx3ujzKLcoWh1xqGc+TKLKsbExsCyL06dP57ThlKJgqhKCW+4xb1rvb/aeotltCoZh8La3vQ0Mw+D3f//38Xu/93vLxksRoKJoCKfTKYuaMirq7u7OK4aEUjfgK2e1tre349ixYyX/YuS75mzniuw+x2L2BvON1BNFEaOjoxgfH1dFxAMDAxUd81YKzn47jTP/mMr6aQMCULbf2K7/0WPxRuvaT3ncdL+ArjfU5qQJK9mSUY7HL+X+oHIPLZPJIJ1OY//+/bqp7WL8FCsx+aXcE230bKOK2a9/4YUX0N7ejkAggLe+9a05DkRL6aUIrDFRLCZ9Go/HcfHiRVNiqDzeaqSoTCkKgoDx8XGMj4/rzmrVolQWSUrniq1bt+q2dhRbMJONJEmYnJzE8PAw1q9fj2PHjqmed6UHghd7ruiEgFceyRbE4pm8lESyflROE9bW1sopwnK9Pss1UjT62DabTTe1zfO8argD8VOsqalRVcHqNcdXai5pOc9RatsoAGhvbwcAtLW14Y477sCpU6eWjZcisMZE0QqkgjCVSmH37t2mxJBQ7PxUSZIwOjqKsbExeSSa0btD4oto9oujfI6ZTAbDw8NyAU8h54piRFG5eJOGf5/Ph+bmZrnxXutaV9Le0Jl/TEG0ljjQxeUFjrx3I7zrNgFYuIGJxWKytyLpsSy1CfFKixSVFIrknE5nwTmwMzMzSCaTOc3xXq+3YqJYzmi01KJIKrFra2uRSCTwi1/8An/xF3+B22+/fVl4KQJrTBTNfHnj8Th8Ph/S6TRaW1vBcZzlqkarokjmgyYSCXAcpxqJZhQyXcbKl1OSJPj9fkxOTpoq4FEaFJu9VnLc/Pw8BgYGUFtbi4MHD+o2/GcfZ5RiBjkUI8DBIQHX/l2dIu59hxN17TbMBCfQtW0znG4GzhoGTjfgdDNwVDGQRIANivjJvQmNiwLe83AtPG0MwmMC6jbY4HK50NzcjObmZtTV1WFubg7d3d0qZwZSfFLMrNJKVp+WGivpzUJzYJXzSTmOQyaTgc/nkwWzmDmw+a6pXGiJYiQSsSyKMzMzuOOOO+TH/sAHPoB3vOMdOHz48LLwUgTWmCgaQSmGpLUiHA5jYmLC8mOaHdUmiiKmpqYwPDyMtrY2eDwe9PT0WD632f1MMjQ7kVhYgM3uWRIhNgvDMEgmkzh9+jRcLhf27NkDj8dj6LiVkj49/Y0UoDi8qcuGmz/rBhgJr7wSRt8h/ff52ldyP0Nb3+zEzX/uhs0O/PB34wi8JsDTyuA3v1ULT+vi4sswjKYJcb5ZpcrIp7a2VvOGbKVOmgFKu+en1RwfiUQwPj6OhoYG3deW/F2pnlWzaO1ZchyX9yY1H11dXbhw4ULOz5ubm5eFlyKwBkVRb1HTEkNCue2jCGRo9fDwsCpdODc3V5Hmf6Vzxbp16+DxeLB161bT57QiiuT1T6VSOHDggKnZipX0UyyGwOUM/P+pjhIP/141bHYGophfaKMTAl79QVr1s713V+H4Hy4MK7/4/RQCry285olZCWe/ncKNn17oM8sn4nqzSrNbGoaHh2VzXLKY19bWlnXPr9yUewSbKIqqiF35c+Uc2LGxMXAcJ7eLKAXTiGtLOdEyGF7trDlRzCafGBLKLYpK54impqacYdlE2MolisrzK8U4EAhYigQKVZEqSSaTGBwcBMuy2LBhA2KxmOlhw2ajt3Q6LU9Dqa2tNTWKzUqkKEkSouMiXvo7dXFN6047Om80lg5/+eEURMVL6mljcPj3F907Tn9T/diXf8ThTf+jBozN2qKp1dJAokqyVzk3N4dYLIZMJgNRFFUp2FL7KpaDcu/H6T2+3hxYErEnEgnNSUhELIlrSyVuBLXSp0tdHVpu1pwokkXNiBgSnE6n4UVeCz1RzHaO0Ns7I8dbWWjyiWK2g0S2GJMIzOzdtJHIjeM4+Hw+ua2jtbUVsVgM0WjU1LmAhffUSGSayWTg9/sxOzuLjo4O2O121QKU3TRPRrFlk08UJUlCYlbC7OUMAq8LCLwuYPaKAC6We8zRe9Su6nqPO3MpA99/qD9/de02nP/nFBIzEkZP8uDZ3OMCrwtYt8tRshSnMqokjfKhUAiBQAAbNmzI8fvLHr/mdruX1WK6nLwUAf2RayzL5rhe2O12efRdJBIx3S5ilGyDYRLRrmbWnCjG43EMDAyA4zi5taIQxabnsoVJKUb19fV5nSO0ji/m3ATiXJHPQcKqKOYrtCHCFAgEcto6SlW1mg3pb5yYmMCWLVtw7NgxiKKITCajev+z01qjo6OqQpTa2lrwPJ9jeQUAyaCIl/8hhdGTPNj5wpHkxn4H2g+r20q0kCQJJ7+azPn51DkBU+fyfyYCry2IYjmRJElzP01r/BrLsppVmku1yK4EL0XlfF1lm0Imk0EwGEQ0GsX09DTi8TgEQUB1dbVqYk+xNyLZfYqr3UsRWIOiGAqFsHnzZkNiSCj27pZEepIkYW5uDj6fD16v17CdUSlFMRQKYXBwEC6XC7t3785byEKOtVLxmi1ugiBgdHQ0byVrqa2jlMbC2a0sWufRSmspC1FisRhCoZBsHaVc2E//jQO+/zCWYre7gOOfrDb0uRp+nsf0RWvvvadN7XNYDvSi0Hzj17Rc5JWLeb7ev1KyEkRRD4fDId9cbN++HUBuu0ggEMi5EVEW9hghO3262h0ygDUoilu2bLEsMFax2+1IpVI4deoU3G63IRslJaUYExeNRjEwMACbzYYdO3bkjL/SopjIjbzGyuKdDRs25K1kLcY6SnkcufkYHBxUGQtbIbsQhUz2b2trU/WqzQzVAdC/wXF5gJYdDrTttKP37S409xhbLId/Ze19P3JPNba+eeFmZjnNPtWLKvV6/9LpNMbHx8tiFVWJ9Gk5p81ki66RdhGzrTjZohiJRKgoUhaxsjcTDAYxODgIjuNw8OBBQy0G2RQTKfI8j4mJCVRVVZm2U7J6XiKK09PTGBoaQktLi6Eey1L4KUYiEVy7dg1VVVV5I/FiohCSMlQ6Njh/h8Mv/ix3Y6/1xhDW3xRHa1c1autIWis3/aqHd5356+x5mxMHP2KsZD78xFUEHjwJfioG54ZatP3342j4je2Gz1UKwdVbzHmex9mzZ8EwjMoqqlRR5UqOFAHtIhgt9G5ElIU9o6OjcguWcl+d47ic9CkVxVWG1cWQCITROz+SpnQ6nejr68PFixctCSKwECmaFSeWZTE4OIhIJIKmpibs2rXL9HmtiJQkSYhEIggEArDb7TnFO4XOZ9WBIp1O4/z58xAEwXAkbAW9z0/XzS684VMSfv236v2/yCuNeMNvbUBVKysPoV4o6gEcXC3E+Tqkp6uRnHYimmzBzg4R3rbFhXr/B6shCsDwr3iEhoy9F2JWTZjWzVz4iauY/svnIYQXq1b5yRgmP/MMAJgSxnKlOUlUTsaCAfpRpdKAmIy3K/RdrXQkt5weX68VJ3tfPR6P48KFC8hkMnjsscfk7xXLsqayXdnXfejQIbS3t+Opp56C3+/H3Xffjfn5efT39+Of/umfltTma82JolVIBWqhL1E4HMbg4CDsdnvJFme73W44fZpKpeDz+RCLxdDT0yOn+aye14wYh8NhDAwMyFFUX1+fqfNZEeFUKoWRkRFEo1Hs27fP1F6xFfK1ZOy9qwoCJ+Hlv18UGi4O/PxP0njzfV6wQQ+CvhYEfQKCQwI41dsiAGjEj0ZncfwLCdTWLkZBR+9Z+CNJEk5/I4VXHk1nn1qmcasNh353MUoMP3EVc196AeI0i9fsDCBIsDVUQYxxgKCxD5vMYOb/nMSk2Al7FYPuW5xwVOuLXqUnzuhFlUoDYmXhSbb7BWlnABYW53IuvuWORMvhkJG9rx4MBnH48GGwLItoNIonn3wSg4ODuPXWW5FKpbB//3488sgjps7x0EMPYefOnXKl+ac//Wn80R/9Ee6++27cc889+Na3voWPf/zjJX1eZlhzolgOT0VgIa0wODgIhmGwbds2zQotqwuI3W4v2BKidK7o6upCX18fGIbB3Nyc5dSr0Sb8WCyGgYEBAMCOHTvgdDpx+fJl0+czU2ijbK9Yt24dqquryy6IRjjwW9XIpCWc/daicKWjkmZqVYu4vwpM1I6ELaaKgshCte19Xoz+2oa5a+rX6R1f8qBxqw11G21yb2L4iauY/MwzkJLXP7fXRVAM64sqAPBTMbzwpYWId/qiC2++Tz8iKOfCb+b7omdATNwvstsZvF4vUqmU3M5TjohuOUeKZiAVsLfddhsGBgbwlre8Bb/9278NQRAwOTlp6rHGx8fxk5/8BJ/5zGfw5S9/GZIk4ZlnnsF3v/tdAAteip/73OeoKK4E9ESRCIIkSejp6dHNt5tNv2Yfm0ppuyrwPA+/34+5uTl0dnbmOFcUW7maT6RIipaYCxOneY7jLO9FFkqfarVXRKPRosbwmcGIcB/6WDWENHD+n/OLjx7hi14c/MhixWYmk5HThdMzU2i9ncfcX28CpMX3efJCGh03eFTvfeDBk4uCaIK0bTHNP3lO+3jyPi232aRKlA7y2e0M8Xgcfr8f4XAYs7OzEEUxb1RphUqIotVxa1aJRqPYuXMngIX1YfPmzaaO/9SnPoUvfelLiMViABZawxoaGuR1cdOmTRX7LutBRdEg2aIYi8Xg8/mQyWTQ09NTcEBuMaKoJchKo92Ojg4cO3ZMcwEpR48jmQgTiUTQ09ODlpaWkvgp5iO7vUJZxbrcZp8yDIOjn6hGhpNw6XH9mbdON9DUZYfdxWDylcX3d/h5XlUo43A4VG4NO3YAjC+Byz9czB5c/C6P9LpX0dCzOKWHn4qZfn4iGAx79sv/joyKSMyJqG5gELgkYOxlHmMvZTB7VUBduw39/x1wbzR9GmPXUqYolLyepPevoaGhYFSp3Ks0KnTlHiNXbtHVMxi2Ogz8qaeeQltbG/r7+/Hcc8+V4ArLw5oTxWLTp2QSjpnmf+XxRotOlCjFifT7TUxMGHKusDqcmxyrFDdlVNrV1YWdO3eW1E9RC2V7RWNjo2Z7xXKcfcowDN74qRrUNNhw5cccnDVAU7cdTV32hb97bKhdbwPDMEjMivin2xcn+QQuCxj+FY+OGxyq11eQJJxJpeG2MTh0TzWGf5UBO3ddoCUGM/+vE0e+4UKKX9hbQ0sVMGsuWhUYF2ZrulQ/+96dUdgcAJdl0BEdF/Had1w48qflaXGqhFcj+e7oRZU8z+cdvUb+VFVV5VxrJcbIlbu6Nfvxi3HIePHFF/Hkk0/ipz/9KVKpFKLRKO69916Ew2G5knZ8fFxVWLUUrDlRtIooihgeHobdbi84Fk6LYnoNybHEU3Hjxo04fvy4oS+Ew+GwLBhEjAVBwMjICKampuSUZb4ve6kiN1K4Y6S9YjlFisrf7f9oNfo/qp3imn41gytPcmjotKF1px2zry+Ky8//ZwKbjztw25c8EO3AqVQa75iYkv+/x+nE//6LRoz84WK0GB4RcfE7Eo5/cmEA9ZV3vxH8o/8Jm4Z5owRAS24cUq6IZtIAdLQ18Iod//knbjRsjOONf1SDhi2lW6Qr0UdY6PGdTicaGxvlrQFgcfRaPB5HJBLBxMQE0uk0HA6HSijLUQijxGhLRikfv5hI8f7778f9998PAHjuuefw4IMP4l/+5V/w3ve+Fz/4wQ9w9913L7mXIrAGRdHsnSfLshgaGkIoFEJDQwN2795tuVjGiiiSkXCBQABut9uUwTBQXKTIMAyCwSDGx8exceNGwxZSxd7dJxIJXLt2DaIoGqrgraSfIlCa3rxURMRT/y2+IDgA6rcsLs4SJMyuF/GiJ4WHL8dwzptGNMtFY5Dn8Uh7DJ94dx2u/HgxRXvhu2l0vsmJaz/j8PpT7Wj1HENn4jyqxAQWZFBC2uaBQ0zDgdzPo3I/0SixURtioxm8+DdJ/Je/UQ+5nhsQ8MznEmCDEjYdduDAh6sNDy6oZKRoBuXotXXr1sk/J1ElMSZPJBI4c+aMoajSCpXYsyyll6IeDzzwAO6++2589rOfxYEDB/Cxj32spI9vljUnioCxu/1kMomhoSHEYjF0d3ejra0N4XC4qPSrGXGSJAkzMzPysPC6ujpLnopWehzJ/t3g4CCqq6tNC7FVUqkUkskkLl26hN7eXsOp6eUaKeZj9nVBFkRgYe9udp2A596VwsBuHuEWxTl09N5ts+H4H9Zg/BSP+AxJowJP3LPY6zFb05WTDgWA1uQQemMvwY7Fz4YAu2o/0VYlwuYSkYktvPeuOhHrDzIYfU77O8Cz6tclGRTx0z+Oyynewad5DD7No36zDf0fq8a2t+dvh1gOkaIZsqPKeDyO/v5+ea9SL6qsra2F2+02LXBLMRwgHo+XZPbpTTfdhJtuugnAgsfiqVOnin7MUrEmRTEfqVQKQ0NDiEQi6O7ullsbQqFQxTwVlcPC+/v74XA4cObMGUvnNdPmkD0erbe3F7FYrOyCqNyrdDgcOHLkiKmbD7ORIjFRjkajsj2SUUf0oqPgWRHTFzM49Q11NTHrFvHQF6LgCo/CBQDUMgz+vKkRVS4Gb/4zN35yb6LwQQqIUHalz8PFJZCyeTDs2a8SUDFtw7a3V6NxK4OG7RmIVUm88vc2ANr74h03LE4sEgUJv/w8u7jneZ2N7FlsfuUkrlw9jmc+14/b/o8HW97o0Hxdl2ukaAYydzQ7qlQaO5OBDgByosp85sNLkZ4l05xWM1QUr0P6/EKhkGYRSSWMhvWcKyRJKioFaoRQKISBgQHU1NTI5w4Gg4hEIpbOawSt9oqXXnrJ9EJoNHojNxyDg4Nobm5GU1MTWJbFyMgIWJaV02JEKLWmoliJFMdP8bj6Uw7TFwXEprTFe7JDMCyIW50OXOhYLIVv6c1dpERGwmiPgHS1hN5LDtik3Nd0tqYLvX+6C7/621wXDsKVJxdSsxsO2BHyu5EK6z/3aP0VnD0rwOPxYOrnTRg/pY4EN7Jn8b7h98Mu8RAYJx7v/B5+9j/60dRlw/4PV6PnVidsjsXrXGmRohm0bKIKmQ8r3S+sOtiYYS0aDANrVBSVCxvHcfD7/Zifn8+xMlJSrCja7Xak09rVCoWcK8p5t6wcFN7X16dyiCimnSMfkiRhcnISw8PDOUPCyXtT6kgxEong6tWrqKmpwcGDB+FyucBxnKaDQywWU01FIXfvxDrKzOIw9Kz2TNRsNvkdqGaBlKJPvjbMoP9XLjz3bvXn5q+aFxfSyJiA771vsfVCgoTn/ksaP7t7Ueh2n3biw19R7/URssfS6VHIpgoAXr+/M+//b06chF3iYYMASVr496S7H8EhEc98jsWZb9rw9hMeec+x3JFiuatDzVLIfDg7qkwmkxgdHVXtVZYSPVFcTp6Y5WBNiiKw8EEbHh6Wm963bduW980uhdEw+TATiCAxDFPWeZ1aJBIJeVD5tm3bNIcOFNvqkL2oGWmvIOleM3fA+aI3lmVx7do1CIKAnTt3yq+x1vPSG5ysvHsPBoNIp9Ny6pWIpd5Q6vHThW+kNrJnsXn2JEIP9+Nv/njH4rV7JCQ96udVBcBtYxASBKQuSfL+oQQJrx7h8c//LTeNeukwjw23pxF6vjpvpKfE5cltwSiWMc9xCIwTkgSIjBNjnuOq/49OiPjVgyze8/Die1RO0ZIkqWyPL4piycRDL6p8+eWXUVVVpfL+VEaVZrYFtMg2GE4kEpbnN1eKhx9+GA8//DCAhRvhzs5OPPvss6YeY02K4tjYGIaHh/M2vWdTbNSkjDTj8TgGBweRyWRMO1cUi3I2am9vb97WkmKec3bEFw6Hce3aNVV6VgsraSEtUeQ4Dj6fD+FwGNu2bTPdQqN8bOWe0Pz8vOzJGYvFZN865Tg2IpYejwebjzhw+Yf6DfzKlOLxf3Hi0qFv4ukb9wIABCfw8i3qY9MA7picAQA0sjZsvNcOrlrCwG598d0wakf/b/HY8N/bMPICj9d/zGH85QykPPc7/R+rhrvJhss/SmPq/OJnwOZYmNjTd4cLj74jqv8AGkhvOIJn+x5H9ekXMOY5jkl3f87vRIMJXL48ZikqN0u5o9BypjZtNhscDgfWr1+v+jnHcfLnkmwLAMgxdjYy83UlGgzfc889uOeee8DzPG655Rb88R//senHWJOi2NLSgg0bNpi6gyqF0XAqlcLFixeRSqXQ09Njelan1XQSwzBIp9MYHh7G/Py8qoAoH8VEimREHMuy8hg8ZaRWynMqRVHZU5kvHU6Os4LSYUCZfiXDHWKx2GKjt1tC095NCF7UvgnITil+6vmLsigWItQqItRa+LX66INeSA8nYXcy6LrZha6bXZi6kFFVqWZz7rE07vp+LXrf7kLQJ2Dg3znwSQl976lCU7f5xf79jwEy7nwAACAASURBVNeifrMd/uffiH8f3q/7e0c+0oD2zR7EYjHMz8+DZVnMz8+r7KJqa2uLHsFWbioR5WrhcrnQ3Nyc43xhxU8xuyWjHO0Y5eLee+/FLbfcgne/+92mj12Touh2uytqNJxKpeD3+xEMBrF3796csWhGICJj9u4zk8mA4zicPn3aUJo4+5xWXydJknD58mUkk0ls27ZN1fycDyv2UUQUJyYmMDw8bKqn0iz5UrXZ49iAhQVpnSOBn/2B9uuYnVLc/qZbsK/KhQtp/ejSKG/5YTXe+m/VYMAgeMWGZo8INighGRRx/p/yT7pJRSS89LUUbv6sG03ddhz9hFrUk0ER7et57OrhUFMtIZli8NqgC1OzTihnBbi8wB3fXBBEAMgk9d/bLW9wYPs7q8Aw1ao0d3t7u65dFInKSWS+XCojl9MwcJvNJhePEZR+ikQss6PKBYuzxfcrGo0u+0gRAB599FGMjIzga1/7mqXj16QoFtu4bfR4juMwNDSEYDCILVu2gOM4ld2NGUjzv5m5i2NjYxgbG4PNZsPBgwdN+59ZEUXSXhGLxbBx40bs2bOnrO0VADA3NweWZRGLxTT3KUuJ2epThmFw+Xvq5x9pWIjwWqdsmEQ/Hu/8HjYnTi6kFL/Vh7vqU7iQFUy96adVaJy3YaJTwGSHgJl2AaLOR6H3VQc+/JAXVenF8578XzU4Cf1055v/tAbJoIRTX19sFbn6Ew7bbnOivX+x1WLuagYXH08j8QqLA31pOK5fg7tGwoG+NHAZmJhe+H2bHXjbFz1o7Fz4pfCYgOe/pF105PIAN37anbMHzTBMXruo7MhckiRVW0Ntbe2SePMtJ1HUQs9PURAEeQ+dbLUkk0l84QtfwPr16yEIAi5fvoxt27aZatVKpVK48cYbkU6nkclkcOedd+Lzn/98yb0Uz549iwcffBC/+tWvLEfqa1IUrWJ0v4sIw+zsLLZu3Yrt27fLvXFWMSpQysrOdevW4dixY7h06ZKlvRkzAkBmsk5OTqKjo0NuebDSXmFUFKPRKK5evQqXy4Wamhrs2LGj8EEKyl3dCAC+/+Ax8kIGoWYRrx7m8OphDiPbFt7HuiCDP/tJAyZ/0b+4vyYBjr+uwd7PcbjYvfg6hFtEvPt7izc1vFPC9CYBZ27kcPIt6qjvo3/thUPIfV7XdvP42XuTsAvAe77jxqbhha//hgN27Hi3C2IGGHiaw/TEwvVVpRk8/0ASdz7mwNhLPF59fHF/8W03cLIgEhx2oH/XwrVMTDvxpv9Zg02HFwRS4CT8x5+z4HUKcd9wb43KXBkobAKsF5krC6NIAQpJFZLI0qpBrlEq4aVYjv5hEn3X1tZicnIS+/fvB8Mw+PrXv45HH30Uly5dwhe+8AVcu3YNTU1NePrppw09blVVFZ555pmFQfU8jxtuuAG33XYbvvzlL5fUS/FrX/sagsEgbr75ZgDAoUOH8I//+I+mHmNNiqLVhdDpdOaN1pTOFVu2bFEN6y62krOQKEqShEAgAJ/Ph6amJlXEZDUNauR10muvCAaDlp6vkdcpmUxiYGAA6XRarpz99a9/beo8JAIxi5kbhdfn0vjfF0I49zkOY925r3+0ScJX3hTFVxxN8P9cXfjy1m948er9UdSFGbz9X2tw8MXrd8/hESBwEU6exebX3Nh0fg9ePdSAeMPiNQ1vy6DndafqXOeOcfj+xxOQrq/T3/x0HJ/6bB3akja8+XqExjgkPPvZNL7rXNhrrAsyaJ2247HvxtA8ZkeLZEdbmw2NszbUVGu/BjYbcGh3Gof2pGEfSIK9UAf3Pg9O/l0Sc1cXX4ON7Fk5OrbffBTb35UbHVi5adFqa1CmCmOxGGZnZ8GyLFiWxZUrV1T7aqUSmuUeKRqBCDvDMOjs7ER7ezu6u7vxiU98AgBMVeMzDCO/JzzPg+d5MAxTci9Fs4bHWqxJUbSKw+EAz/M5/UCCIGBsbAzj4+PYtGmT5n5WKQp19PoklU3/Bw8ezPFYK0e/oXLyTlNTE44cOQKnc3EhtnoTkO84nufldLSWZVUlKCSKaUnC18NRPB6L4yLHAe/J/3iTnQK+dTqG29d5VI396ybt+O0ve9H9ugMu7vpzDI8Ak6fBSNffS54FJs7grY/swQ//aNHZ4cW3pTGxVUCkUUS0QcTFY7mLV9Ir4cm/YPFU03o0dCx8Vv8qGJIFEVgQ7WhTBj6N6z72Awc2snoFTAt/ixEB0SfDmHs9g0uPL75m2U384bc/CYY5lvM4pYq2tFKFgiDg3Llz2LhxI2KxBUNnn88HQRBkX0USVVqZVboaRBFQr1uRSERV7ar8vhtBEAT09/djcHAQf/AHf4Du7u5l56UIUFE0RbYwkZTo2NiYHCWVaySalrBFIhEMDAzA4XBoNv3nO7YYtKbfaJ3Tqihmi45y8k1HR4epYqFK8/GZWfwgbq7B72fvTSFaJ+Gmn1SjPrQoAjsvZC06gYuLgngdRhJw6N+uqkTxtUM8XjtU+C7+tSYec1tENAP4v7E4Hgyppxcdu3ABN509g+f6D+GlfftU//fgQQF/ddIOt0aaVgUvIXM2DmDxs6msuGUYoO7yi+BuyxXFchsY6/WlklmlxLyazCrNLurJJ9hLYetUbqLRaFHtY3a7HefPn0c4HMYdd9yBK1eulPDqSseaFMViPRWVKcO2tracKCkfVr/oysHeZvscrQpUNvF43HB7hVV3DmWkKEkSpqenMTQ0lGMsvFTkixTTkoQf6ghiPoEBgBffnsaLb09js8+OQ89X4eizLtgkBnXtNkQnrr93OhtyzkQSjAg5NWqUGgCOaBQnq6vxycBczvX+8p7fhYvnwTmduPXhb6qu+6kuCYCAL71oh0NjhJzqPFmp1jHPcQg254J1VZULmRtu0DyunPtyerUB+XwVSf+fcqqMsv+vtrZWXgfKvado1bDcKFqf8VK1ZDQ0NODmm2/GyZMnl52XIrBGRdEqdrsdc3Nz8iaz2UrHYuYV2u12sCyLS5cuIZFImHKRsGpbRUgmk/D5fEgkEobbK6ymT0mhzfz8PAYGBlBfX1/2ilIz5BPFAY6H8jaAEYGtVx2484lL+MIPfg92kUO6KldglIx1CxjrZhGrF/G2H9YsCiIAON3awuh0Y/tFB67sN/Yet07a0OW045ObGUyHQviAZENasYC7ADzw2qtw8Twcoggpk8FNZ8/kXDMRxgeed8CZZ/1PptSiWXf3ccSPPgn3uReRueEGiEePah5XzkjR7Ig3p9OpOVVG2f83PDwsG4kTI2KWZXWnHRV7/ZVOzxYTKc7OzsLpdKKhoQHJZBJPP/00Pv3pT+Pmm29eVl6KwBoVRbMfUDKebGJiQp6dmb1vZwSSxjT7YeY4DrOzs0gkEujr60Nra6up52A1UuR5Hul0GufOnUNPT4+p81oVRZ7nMTAwALfbjT179izLsVJ6oniFU/cWHph24QMPeHB4+iwcIg8bRFSltQUmm9cP8HjbD7PS0m17ISn3FAFIjB1o24vf/LYHP707iVCLIFe3ZvOhr3iw9/TCzcWb/7QGO49X4bNzQYTC6rRpNSTct2s3/sPphJTJgHc48Fz/oZzHa0jZMOd34C8PZHDPoA0brgfJNoV9cUYAXhtcvKG5/e+82HjQAeAYuFtyU6ZKyh0pFvvYev1/6XQaPp8PPM/LLQ1k2hGJKIvtqRQEoaw3inoGw0b7jbOZmprCRz7yEQiCAFEU8b73vQ/vete70NfXt6y8FIE1KoqA8SpCUsTidrvR0dEBAJYEEVhMvxr9MGcyGQwPD2NmZgaNjY1obm5WpXSMYrPZTEWKyvYKh8OBQ4cOmf4CmhXFVCqFwcFBzM/Po7293ZJ3pNHIgqS/p6enVa4YRgoq8v3/FU69j/fGvmp86Ik6TH71zRD+9iuQxIUGfWnmjQWvUdCa7dyw8PmTAhcXIkanG2jbCzR0oCEEfOAfFm4g/uEzMfh3qN/vntcc6Hlt8ev++joefzoRwjPJ3IHgUTA4uW8/bn34m7opX2+YwZ98w42vfpDFS+0i/u++hff6XUMM/uQVOzYmGLmhn/QtfuiJupy2i3yUezZpOR6bYRhUV1ejpqYGdXV18sQj0lMZj8cxOTmJeDwuR5PKop58VlFKsueSlhotUSwmfbp3716cO3cu5+fLzUsRWMOiWIhwOIyBgQE4nU7s2rULXq8XgUCgKCslo04bympW0toxNzdn+dwOhwMcV3hCCpkKMzIyIk+FOXfunKUeR6PRaSaTkXs6e3p64PF4LE37JwU6hRaU+fl5XLt2DQ0NDdi6dSuSySRisRgmJyeRTqfhdDpVg77dbnVDeb6bqatZr/EOlwve10+jr/YUMl9+AFcencLlyFE0XT6Io8+wOXNNlfuOAPDRG69hrOY4Lozsx+yV69FfQ4csjnpsGrLniOLgrgwefCCK2//Zje4rDny+aQ6JPNNlAOClfft0I9pbn6jGN+9IItCufo+f6pIwFzuL3/3OebRfaUVXTRjdb3kTOu67AZ4WcyKU4URERoCZSQ6ZtISdt5fOBaISe37KSFCvp5IU9WRbRSmLerSGei+FwXA6ndadWbyaoKKYRTQaxeDgIADkOFeU21NRFEVMTk5iZGRELiwhH8xiKkgLFb3ka6+wet5CkSKZuEOEnwxmHx0dLWovUm+hi8fjuHbtGhiGwb59+1BTUwOO4+D1elVTUshA5VgsJk/KUQ76djgcuteXHSkevnAe7jvvBDgOot2Jy+3fw6S7HwyA3/iOGxIDnLp5QRiVhS0Zux1gGDgyGTS4XNj24x8jsf0wnvxEHCF/4dfmtsdrwFdJeOlWtejGGyR895PFW19sGLHj9QM8RntzPxfy8+A42CUJYBhIj/0tJm79DhI33ghXTQ3+JRbHNyMxMADeX+fFh+tq4YgA84OC/Cc4KGJ+aD2kjAiAhbuFWdGiqEU+A2JS1EPGr5HB9EoLs6XwUlyuVd+lZM2KYvYdP6no5Hkevb29mmmCUoiilsBIkoSZmRkMDQ2hublZs5q1mHPnEzbSXuF2u3HgwIGc1HAx/YZazb3K59rW1oajR4+qvnxmJtoYuU6O4zA4OIhYLKYqEjIzUFk5Toy4ZJw+fVq1R1Tl8cCX9Xy3vfQywHFgBAEQFv0DAcAhMLjz2x5sHLHjR7+dxE1nz8iFLcz1a7NLEniex89//u9465EjeNdDXjzx8bi6+EYDh8DgNx/14MCvq/D//msiJ5orlqkOAVM6/yc/D0mCBICRJEh8Buw/v4TvsNvxWJeEidpFMbo4F8SJwRB+94u1aJvKXuQXF2B2TkIyJKKmsTRCtpwNjLU+g8Trk7iyzM/PIxqNwu12q6LKUg1K14oUKzEBajmwZkWRwLIsfD4fWJYtWNFZCqPh7OOJv2BtbW3eAp5iIkWtY0nkBCDHXLgU59USqVAohGvXrsHr9aK/v18zTaonpoXIvslRumV0dXVh586dRbXikNRXS0sLBgcHsWvXLtWA6slEApm6xSKEWgbAG44DThdEgdP0DwSAgy9W4Rf/XwrP9R8Cd72wRbDZAIaBKAjgHQ48uHsv/mZiGo+tb8O7vurB9+6M5bV9Imy95sCnPluHZ9+VwjO3pyAY7LU+8qwL7/hBDQQ78OMPsrh4VP1+5GsvIc8DPA+7KEKEDRk48Wf9h/Bv+xgohY4QaZDw7LtTuOsb+Yuq5gcFbDq8MkSx1OnN7J7KZDKJ7du3Q5Ik+YZtamoKqVQKDocjp6jH7HPN7oNcir7IpWLNimI6ncbAwAAikYjh6SilTJ+SPUuXy4W9e/cW3DQvlSiSghaj7RWlEEVlfyPZnzVynBmIKEqShKmpKfj9fmzcuFE1aq+U2Gw21SLVJgiAf1T+f0mScMHthuvzXwL312O6/oHVKQb/9a+9+Pv/tU8ubJmrr8dtL76AjbOz+NZ77vj/2Xvv+DjqO///ObNVZdWbLcmSbTU3JNtyAQyYkhBqEuAcCAkhkIRcAiE/0rhLuXD5kgqXRhoJJAQOCBBaICFUU90b2JZVrN77Fm2fmd8f8qxmq3ZXK+Oc/Xo8/ADL2mk783nNu71e08TjdlPV2U2twcAZD5pY8k09ps7w+/X0L5mRPDDSLNG/xw92+MDTaZy208h9X7MzWRAeIS/s1HHaTgOn7TBSMBy88H3inkyerXDyVsm0puls84vb66fP47JX91Czs5i+umEe/cgattevjHk9W1b6UFAQQkjTkO0nq0Imv0qHzuJFUfQpiVb+1RVn1DlFvV4fJpTu8/kCL2y9vb1BQunaqDJW85zf7w+qH/6rOGSkAictKY6MjJCXl5dQBJEKUrTZbOzbtw9ZlsNqlrN9di6k6PP5aG5uZmxsLKHxirmmTw8dOoTD4Yh7rjIZ6yj1cxMTE3R1dc3rbGO0RhsxZDF3IPD7ojLOKVw6q+fhB+uNXL9gIY3HNvHaTZ/BdKxpp6GlmYNV1QHiafb5aMYH34NrH0tnVJCx5slUju9jVc8OvutvpH/lWq4pycBzCWx4az9lL77B0ZEzef3wMnafHVxn/OrXsyKkLWEqU+bARi/vXOBhLH/m+LVp3mjzi2qDjigR0clDkKH6oJ6W02aeJXuugv90WJlnJL9KJL9KR6/zEMtXLz1mG2VjxN5Hz66ZSEhr5pzoi8+JnD6NB7FI12AwkJubG/TCqwqlqyUA1VPRZDKFeSoKghCWPv1X8lKcK05aUiwvL0+YZJJdsGE63dHb24vT6aS+vj7heZ9kB/AlSaKvr4/x8XGKi4sTlkhLJlL0+/0MDg4yPDzM8uXL4zI0VpEMCU9NTQWUMeZ7tjEaKUY6vYdcDh5qcMCDUNwrsuXeDMo7wh+5rDN6mTzs4PdeL0f37Mbg8wUo1hCFeAD+d8v0IP/GAwe471vT0duXH5+O3r5TX8/GAwf4z7unf75M/BV3Z98LzGyntF1Hd5Ufo0cgZ1xkvFDi0Bofh9b6aF/mD2z7hn/OpEoPVKzDpzOC4kPBQM7Bs6jbZ+DI6vCUdyRCrNuv56K/pLOgV8cfvmYPIsb0/9JxbvZMxmRor0xaWhoWiyViJGS328PUZbSdw7EUX2Zz4Jgrjoc2aaKWbNGE0tWmHlUoXTvCpdfr0el0WK3WpAf3e3p6uO666xgaGkIQBD73uc9x6623Mj4+zsc+9jE6OzuprKzkscceS3oOMpU4aUnxeBWMPR4P7e3tTE5OUlRUhMViSeqLT5SQteMVJSUlWCwWysrKktpvvCQlyzJ9fX10d3eTn59PYWFhkIBwPEik0cbr9XL06FGsViuZmZnU1dW9b8P+/lm+m6EymV/+tx2ATKvAv/0+g9p39XhN8KeXChj4VB6v+p2sXduIz2BAPBYpRhuc1yJa9Bb0c8EXRq59SyQe+1wULyfCU6V3XvQgOa3n8GTZjP+jOLKWG/4Huqr8vPBvLo4u9wd9Xq099mSt4ZJH06g6PFPYXNpkCCLFN11ubsyeSdFFi7YiRULaRpSRkRHa29tjinvPdyQH87vGJPtyroVWKF2dp4Tpl9r33nsPURTp7+/npptuwmq1YjAYuPPOO2loaKC+vp7S0tK4zlGv13P33XezZs0a7HY7a9eu5QMf+AB/+tOfOP/887n99tv54Q9/yA9/+EN+9KMfzfm85oqTlhTnGz6fj87OzoCnYl1dHTabjZ6ennndrzpe0dbWFuhk1ev1DA0NJbW9eCJF7T4LCgpYv349brebjo6OhPcXDwnLskxXVxf9/f2Ba5usZ2SiiBQp7nd7uGJgMO5tOLIV/vjVaTcKnY/pBhgJEAS219dzy9e+wY1PP4XbaKJpyZJZt6dt0tGSaLSfxwstqeLzc8GO3ezO20h/+tqw+mhFm547vt1CWsbbPHh1A4OlcoBQJcHIYxUPM5gWvP+lTcHLz5suV1CHYyLdjrOJe1ut1oC4t8FgQJZlsrKyglKGpzANvV6PKIqUlpZiMpl44403ePrpp9m2bRvV1dW88847/OlPf+Kxxx6La3sLFixgwYIFAFgsFpYtW0ZfXx/PPPMMW7duBaZtozZv3nyKFN9PJPsQzPaWqVWD0c7fwdxrkrNB7e6MZiGVDGYjKavVSnNzc5j8ndfrTbl1lFYgXOvdqH4uUVJMlZ/ijyYmGZWSG3vQdoRuPHCA6577G5/+2zPofT50isKm/fv49LNP88fLP8KfL70sYhpVbW4J7QiN9vN4MZqdjSwI+AURRTHQbw7vnlWhtYO67NsGti+6ApPHhw4ZAR8Vju0MpK1lvFCmo9Yf+KPFiCRzxOtjmWmmFjwXsoom7u31emlpaQkIR6gpQ21EmZmZecJ2Wx6P0YjQmqLD4WDp0qVs2bKFLVu2JL3dzs5O9u3bx4YNGxgaGgqQZUlJSdIv7qnGSUuKySKaVJs2daiqwYQ+VPPhawjBg+mzdXcmCp1Oh8fjCfu50+mkpaUFSZIiOmbMpUEnErmphG+xWCI20SQ735gMQo/PEuUFqahPpPZdA7vO8eCeRZFLTVWavF5ERUEAFECvyOh8Mjf99XFufOYpvviN/+BgVXVEAoxGmImSoXo8P7/rx+gkBRB5ZcF/ReyeVaG1g1IUsNhEZGG69ijpDDzwsQae+LAVW17sF5dOvz+IFOcDRqMxkDJUU7B+v5+pqanAaIMqwxbqgnEiCNMfj3plqGC6Wv6ZCxwOB1deeSU/+9nPwjpZkzX9ng+ctKQ4V/so9eFQo5eOjg4KCwtj2kjNNVIMjVJdLhdtbW24XC6qq6vnpUgdSm5qHW9ycpKampqgAeNYn0t2f06nk+bm5lnHOY4XKUa6b/47P5cjXi/7Pd5AHW00O5sCq5Wtaxupea+R75hzuap+OOI2c0dEbvjffZi9PsRjQ+8y0xN9kiAgKgoi0003v/7B95F0InpJijgSAbPbVMWDzXt2Y/JOR3oSAmnSRMzf78k4HUkwHHthEOjIX84Dl/6JPHEbf9+8lu31K5im+ejHmCUKrJlnQlQRmu3R6/VkZ2cHNZOoHZsOh4Px8XG6u7uDOjbVqDLUBWO+0/jvh8Gw3W5PSo9Yhc/n48orr+Taa6/liiuuAKC4uJiBgQEWLFjAwMDAnEk3VThpSTFZGAwGfD5fwDmjra2N7OzsqMPoWsw1UlQ/r6Z9xsfHWbp0acKuGcnsUzsMr9bxYu0zWWcOldziJV8VyaRPk0lDRUqfluj1vFa2kL+/upWPfP5zGL0edIqCJIh4jdPEdfWyyOR015daaWzZiVtXhKw3gM+PjI63F1/FQx9fypojR/jM00+ilyQEQCdL6GQZESXiSMRsc4TxYuvaRjxGA0avH4Vp4QFBB0XLdZRvMFC2Xk/RMh1vv3iAnKlldG/bwKsv/hfn93wbAZkLe+7gx9/8PdvrbwjablCK2O/HazDwhfv/iGXj6VwrZeLfr3C4x4O1R6bphUrac2yk5Yls/KKZomWpW67isY7SdmyqDWOqC4ba/To0NBRwwVBJcr5rlMeLFLWwWq1Jv3QrisKNN97IsmXLuO222wI/v/zyy3nggQe4/fbbTxjbKDhFiglDr9czMTFBc3MzZrOZ+vr6uNXq5/qgiKIYEM+urKxMeLwiWRKwWq1s3749alo42rEm+wLgcDjYtWsXlZWVs5Kv9jiPR6MNRI4EdILAlfv3YvT7pqXNmE59Rpvlg2mC+PddKoEZ+dZXv0rx1BRDZCPIExysquYPV17F3ro6fvXD7yMqCn6dHlDQyXLE5pl45gjjQWg9smHBMq6q0rOuMCPo+zAX+lh1sYmVVxl55PZh+J2MTpYjjpJEShELHj+f/fp+dhZv4k1PqGOHAe+EzESHzMvfdvLxJ1I3PJ6sr6nqgmE2m4M6NiMNzO/atStoTCQzMzNuM/JYmG91mUjrxFy8FN9++20efPBBVq1aRUNDAwDf//73uf3229myZQv33XcfFRUVcTfuzDdOWlJMhqDsdjsjIyPo9XpWrlwZ9+D9XKEKhU9OTmKxWJJyoFcjvkRms0ZGRjhy5AiiKMZMC0dCoiSlaqK2tbUhyzJnnnlmQueYSLpWkiTa29sZGBhAr9cH2UfN1mAR677xb9qE0WhE9ngQZBm/KMbs+gwmMB+C3c5jaxvDIr0/XHkVB6uque65vwGgX1TAVfZushU/j+x/nv9Mk3ikZg0QX8dpvOlVbT1yOx5+axsB2wgA16Dw7axMZFlG8ss889cp/rR+NZ+8P3zfOh8s32fgK795F7PHPx3lAjICsmCg03g6/vCydRC8jtS+8KR6JEM7JlJQUBBI96t1ynjGROJFos9xoohmMJxspLhp06aoa8Err7yS1DbnEyctKUL8C7fT6aS1tRWv10tBQQFZWVnHhRAVRWF4eJijR49SUFBAUVERJSUlSb0lJpLOtNlsNDc3YzKZqKurY2BgIOE33EQe8snJSZqbm8nMzKShoYGmpqaEzzGe71L1Uezs7KS0tJT169cHKX2oDRaKooQNgguyjo43fHRs9TE6WIJ4lptFZxjIrxIRBIEBv5/7qmp473f3Ur9zF6PZ2RSNW9m/ZB1N5fVo62kqIhHYdc/9DdMxh4nQSO9Tzz2LqSwfMXM5gn76+lQ6Jvn9a38F4JGaNWERHsDt99/H1rWNWKVGbrr/MF/Y9Tl0sg9Zp+ONK6/imx/8UMLR5CMIPGKb4sz0TGrv68L/Xgbbvxy87x7LGj78gJGG7UYyHCIe55lIwj0oCijoOJi7hUM5V8Zs4FGx8Za5d1JrcTwMjKOZEEcbE5nNLkrF+2EbNZfh/X81nNSkOBvcbjft7e3YbDaqqqrIz8+nr68vpc0y0RBpvKKpqSmlouChcLlctLa24vF4qKmpITs7m6mpqXlrYNF2sKpNNKozd6KYLVJUU95aCTiv14sgCGHzbbIsB97wO/aN0/WSldGdFiSnuhBlsvOQmx2/dTG8XmLnR7xsK/PgF4CVp/GP73kAZQAAIABJREFUlacF7Xu1wcC3CnIxAA9tc/BY5bQCSyQC+/SzTweabSSdLvBzNarUra0CffCCmOH38fMdL2IaOpvip0TyRzbTV7GJtGVv89xtnwmKOovWHkC/w4cgSQiSxLmPPsIrTz7J+b+9N6k069tp6bx9PnD+zDltr6+ntF3Hrf81fU3LNvvxTYqMtKzlMWYG/+MhQ4ALvpdO1QWpbcCZT1KMRVqxxkRUZRnVskxb09RmMeY7fRqJFO12+ylSPJnh8/lob29nbGwszGFBr9dHHFGIF9FGOlTY7XZaWloQRTGs23K+7KPU8x0fHw8TR5+PMRKfz8fRo0eZmJiguro6qDaTbBdptEhRJV5Zllm5cuWs4ypDfj8PjNk50uVFeU+H7kgmWZMi2TmQpSiYXeA3wIGNXt7+gIe+xbNfm30+L1sGhvhSTja/PK+A/jaJt0Q3EJyivP3++wINNZIg8Mz6j/AJy3mckyvx7vr1eH9vQJdpjuAzAfm2cZbfOxPNl3bp+cBP3w2rL/5ibSNXGAwYZRkUBUFRMPt9fGjvnqQ7VSNBPray5FXDaf/uYmrKgd3mwD2cz8TAteh7MuDl2NswL/DwkZ8XkFOeegJ4v0gxGqLZRal1Su2YiKpEMzY2Ni9jIpFIUZKklNRD/xVwUpNiLLuhiooKqqurwx4ctfs0WUQjRe14RU1NTUTx3VQbDWuVYSoqKiI27iQ7WhEJWmPhiooKamtrw/aXbMNMaPep3+/n6NGjjI+Px9W9CuAYkrm8fZCmQh/kAecc+6OB0Q2CAp4QA3JtnW5fbT0Gs4JDQ18S8NNJK11vv81PDr3HEw1r+Hl17XR0eQyh6dSf3Xgp124Q+GZ2LnzkchzFf0N84heY5fD7rzczB79OQS/NbNC0/CyEnl/i93lRmB7G315fz+bf3Mt/vvB3LnnqSfD7EQSBr1ctYWlJEdcNDkdI9CaOgUXT99ryy9NYtCjYbV6NwpdcO0HzI0YGXgyvVaXlCyz7ah855cVh/5YKzOcAfKoIV6fTRRwTUWuTExMTdHd34/V6MZvNMcdEEkFoJKq6z5wsOKlJUYUsy/T29tLT00NpaWnMRpZUeCpqycnr9dLe3s7ExMSsFlapso/S2iuFKsOkcp8qtLXRSMbCWsxVaUhRFHp7e+nu7mbRokVxd+iOH5V49hYHrXfFfuHxRihtaccgJNHIw9UPU/PEKr7rmeJlDTFuPHCAB4/9Xr3RSMMjj/Kz6lpapnyIjsgKNDtHRsnViVyemUHa6RtRDD6UR3+NIM3cgy5Rx+2nfwizw8ey7dMvWwude1gs7cb/75/HcM89CLLMz+/6ccBx4/L6ej5fVc0vf/QDdLKM+fbb+ejy5Xy/to7/GB2P86rHxvACic63fKy8amZUSVtnW7gQqv5T4Q8vWsM+u7jwVUoffoXBoS50mzal1EBXxXyR4nzW/ERRDJClKpIez5hIIr6KoZGiSognynD9fOOkJ8X+/n46OjooLi6OuVirSJWnojYqraysjBg1RfpsslGqSm5jY2O0trbGba8010hxcnKSlpYW0tPTUyY9FwmCIGCz2di+fTv5+flxfZcqRlslnrvFgduqUNKro3dJ/C8BggJb/rZ3RtJM9lFh3c6zh9ZQtSKbbV47U8cWlc17dmP0+tArMni9HHrlVfZULAYD3FKQTsEjTvZkb+TuT9XjO7amysANg8M8vrCEc9PT8Deehf/vz6PvOYiQYUaZcvNM7mIeqVnD6bl+lm03zkiudfgQ9AJIEqKiYAxp3MmxWkGWEWQZ2evF+9wbLH9pFbeMNJGR8U5Qd6rJFR4dz4a7fmzjG1/JYmOLmYKayN/FnvvcYT9b6NzDRX+7CcHvRXnsMY7+7ne0Lp22kNI2pFgslhNSt/R4NMJotx/PmEiom4g2qgxNi4aSotvtDvJW/L+Ok5oU29vbcblcCXnvpSJSHBwc5PDhw7NGpZE+63aHLyLxwO/309raSnp6ekL2SskuOC6XK5ASjiQDl0o4HA66uroAWL16dVxzo4qiIMsyw00+/nGbG49tmrga3zDSuyR4Xu4Ms4kOn5+BCBHzY01FePadgyz8FhQfsmDgns+fxlNVVggpPY9mZ6NTZBRAJ8uMatJiD/sd3H9GD2dbLJwjZPF19Kh78wI/e+FFVrQeoWTzZrwXfhT95X9E9npx6/X8/LefAWBHvocLssyUj2gk12QRRacDRUE0GoNGNALapqKIV6/n00W1fPCJt7l7//WIihdJNPLA8v9lTGyk3L43qEHm5ZtcvLhp9nvxR3fbKL5J5IYXcxB1wffSSLOf/Q8HXyRztsAFG/YgPOBDkGXw+VjU3k7JRz86fS0iNKTodLow3dL5dsCIhePhpRjPC180X0XtmEhHRwd+vx+z2Ry4fm63O6j2fjJ1nsJJTopVVVUJE1yypKimEPv7+8nJyUl47g+SS2W63W7a2toYHx9nwYIFVFdXJ/T5RKFt2jGZTDQ0NMzbTJXWOqq4uBhRFGclRDUVJEkSw00Sf/+yC69j5t8bthl5/pOuQKQGUKjTccgbHqFnCAI99/rwp6/lscrprsrmog08dWlkOawCqxVJENArCn5RZM2RI4Fxie319Qy58jhr03KW2O24rTa+6Z+O0DceOMBzx9Kuwt1343zmGZzPPovuzTe5saaO7ctXACCLsO8ML7WD05JrggCYjHR9+cvkKwoHNm1i+8LywDZ/ftePEWUFWSdy61e/zjMXraK+8z7EfV50yCD7qB3eQU+GGBD7lgQDj1U+Qq2+nhcJJsUcQWRSCc8qvFft595NVho/Y6b2EiOWEhHJr7D1TheK5nbOKBTY8nAWaYfOgYfvQlYUBKMR/6ZNgd+J1JDi9/sDqcO+vj4cjukvNHT+dD5n+7SQJGleNVLnEolGGxNxu93Y7fYAWQ4PDzMwMMATTzwR8KaM1IAzG2644Qaee+45ioqKOHjwIMAJ66Oo4qQmxWSQjHzZ+Pg4ra2tZGZmUl5ejtlsTqqTKxFC9vv9tLe3Mzo6SlVVVZgAb6qhrcuqtbzdu3fPyziHLMt0d3fT19cXkJwbHh7GbrfH/JwaHcqyzNBBib/f6sQ7Ffw76y8zcZkljSedM9HiM1ORPQenFAWbWyYdMWCn5BcVjO7JiLXHrWsb8RqNKH4/kk7Hp599OkjD9HP19ex8YZJvnJPDLZW5CJNW/nN0PGjI3+fzMfD44/Recw2Zl1xCWVo6+Geu8Qv/5sKVvhz56Ye54tL9KJvPYiQnC1NZGfe7PGCbvkaBbSoyPkWgwDpd13vxjEa+8WcDJq8fDAY86zdRvvOtILHv8qltvP7aabBx5tw2mE08u7CEnw1Z+cHUZNB5//GrDn78yVx2/8HN7vvclG/Uk54nMtYqsdC5JxCB1n/9bEyZAvKGDYw//jiO55+n4MorkTdsiPm96vV6cnJygprTtBHR0NAQR48eRZIk0tPTyczMxO/34/V654W85jt9murtC4JAWloaaWlpFBUV4fP5KCkpIS0tjd7eXl555RXa2trYsGEDBoOB0047jV//+tdxEeT111/PzTffzHXXXRf42Q9/+MMT0kdRxSlSnEdEGq/o6+ubc10wFrQdnlrrqv7+/jmNkkSD1kuxsLAwqJaXrP5pPPsqKioKSj3H6lrVkiHA6OvvIh5+nas/5WDKnsne7RvoaK1m9SfNnH5zOj6nPogUY0EOyZLpJYElTZGd6LfX1/Ptq/7M6Qd3Mrq0m888/dcwObY/1Nl5oMfOpZ50vrYsh09lZYZ1pT615EO0tVRyMM3LvlovaCR3/UZ45aNudp9TxQsPn0b9bQbSiyUevNbFAw1TqH0/0ZRvtM0+peeewyfOO5uuH+iQfvILFAVkYVoHNWc8+MT7/H7SRJH/WJDL021TNBF8/pKooJMFUKBn2/SLndZuStYb8Bj+hsw0AXrXrmWwsJC85cvj+h5CES0iUoUa/H4/hw4dCgh8ayPKuTb0vB/D9fOxfaPRyGWXXYbJZCI3N5e77roLp9NJU1NT3Ps/++yz6ezsDPrZieqjqOKkJsX5KtCrQ/ButztsvEKv1+NyxbfghiIWKaoyae3t7RE7POdj3tBqtdLS0oLZbI7YRDMX/dNQ2O12jhw5EnVfkeYbQ8lQEAS8bQcxtv8Tfeb0wpyZ5eCMc19n8dkG6j63BkEQODc9jTKdnl4pOCrPRcGugF/riBDhFqo5qA8jRb0Xrv5dBhNKI1ffW8e6wwe47vlnI8qx+YzwlNHJU73HItSwrtRqwEEsWPMUHr55iodvhvrtBg6EHE8sr0Xt7OQtHV184abl2IX7Wda7g3+etZa3G5ZgsQbvf8AvISkKOkHgoUVFrO3uC/73comyruDlRms3JSigf+stvMeiQnUeT9yxA/1bb+HftGnWiHE2CIJARkYGZrOZvr4+Vq9eHejcVFOHAwMDgYYetU6ZaEPP8agpHk/StdlsgTUsPT2dtWvjE12IhhPVR1HFSU2Kc0GkOad4xiv0en3SRBEtfaqq32RmZkZ165gLKUayrFKVb2pra6OmZpPtXFWjPkEQ8Hg8tLa24nQ6qa2tjVrw184pauuG6nbU78G3fyt6ffA11Bv8lGe9gyBMP+w6QeC2lmy+WTKOK0Oh+j09G7aauO1mPxfp/QxprmMkUlyxx8DfP+bCfywzl2ETuP6nmTgzFP58qwPJkLgBcCxfxNm0TA9sDCZEUYKKVj07V8bntfhrqw0+tgJYEfjZRGHw9yoBg5JEqV5PVYTSgPzvAgV/1DHaPHPtVLspQQBCaoeKomA5eJD0m2+GY2lO57PPzpkYIZi0tJ2b6ogDzK2h53i4WMxnx21od+vk5OS8NdqcSD6KKk5qUpyrp6JaF5Qkic7OTgYHB1m8eHHM8QpVpikZhBKbai4MzGouPNfBfzXiUuuUqhJNrGs4F1L0+Xz09vYyMDDA0qVLKS4unnVf6pCxeqyCIIQvWk5bxM8rThs+t4LBPL2Ps/xpfOuW6YXA4Jv+2Tv/nwl+akUrKaNECAhyx3Rs+X0Gb1zkpmBQx4VPmJnMlwOEqCJZA2AtkrGK+uTPM1ixz0jLSh8P3TyFOyM1g9lezTzbOpOJXZp0/QfWpXPuOWmMNPs58qyXtpd8DIhrefdLf2VF9g6ks84KIjxZlsnauxe8XgRJQvF6gyLJuSCeSG4uDT3vh7VTKhF6/Ha7ncWLF6ds+yeqj6KKk5oUk4VKijqdjr6+Prq7uyktLeX000+f9WFLhVSbx+Ohra0Nh8NBTU1NXJ1bcyXF7u5uBgYGKC8vD9Qp4/lcoqSoKAo+n4+dO3fGfU1VeDwePB4POp0u6huokJGFMhVOjFP2TLZ+3c7Fd1nwyR48JW0YfAuCfsc5puAZU+DYKNjGAwc4c+A9xuUzwnQ8G7ZPC2GrePiLU0GEGA2fvisTv0Hhb9c6mSyITFZ5bpHTPEbWmE2U7NmbsFXU8n3TB1Jz0MAX/9vCH7/iYLxobrXf1SYjFZqU2+MLi/lo/yAHPF6uyMzgrLTpdHdhrZ7Cr+k562tqtuUsfJwVtj1FUXA0NoLRiOL1hkWSc0Gy6c14G3ocDgeSJJGdnR0gy1Q29BwPdRnts6NNn6YCJ6qPoopTpJgEdDodQ0ND9Pf3U1hYmNB4xVxIUSXE3bt3s3TpUpYvXx53tJsMKapGylarlYyMjIQG4iFxUlTdMiRJorGxcVadUvUYZVnGbDYjiiJ79+4FIDMzM+BmYrFYAm++xjXn4nnnedDUC/0+PXu3b6C71cdT/98QOZsGKckvo+ETZvY/FDx2kDsqMlkga7wBfcj8kscqH4kpcC3H8TXV7ddT856eDZ8xcpvXwu19k/yzNHwW8O26Ugp0Oj47OMw7a9fwxZCGGTMCGw0mhpBoitDUdbDRx6rd04t0cb+Om79r4cEvTdFRF/u+/H8eFzXOKdJ1OlwWC+MGI212O8srKrggIwNRcy/m6XS8Xl6KW5YxRyGgWPeuLMu4V6/G+eyzKasparedqppfpIaevXv3UllZicvlCkixpaqhZz7l6aJhLnOK11xzDVu3bmV0dJSysjLuuOOOE9ZHUcVJTYrJ3Fzj4+OMj48jSVJSCi3J1BRlWQ5EpIIgJBQ9qUiUFLX2UXl5eSxatCjhjrd4u09dLhctLS34fD6WL19Oa2vrrPsKrRsaDAZWrDg2ryfLgcaJ/v5+7Hb7dI3KYsFiySK7/lz0zTtQpmw4XZnsfnO6+xRgaJeBoV3lNAMQTkiVLXo66vwzrhWKDPgon9oWkxQ/8ud0/nLTFJP5MrmjIvnDIpU6A2mtMtmdegoHdBQOThP37nt95C0R+Z8v5/OHbBs/tQaPmnywp588RA5IvghNONM2VVt9bs5NM/PGolIub2ljm27mev5ji4vl+wzojumkZtpFbvp+Jt1VEooAX7g7m0X5es7rHWCPJgW6pKKCD2VmBJRSxsfHWWSzYmw+QtsxRwd10VclxaIR4mxQF395w4aUpEy1mO9GGFmWsVgsQUSSqoae+T72aAbDyUaKjzzySMSfn4g+iipOalKE+AWo1fEKnU5HUVERhYWFSUmWJVJT1I4gFBQUsGHDBnbu3JnUQxEvQbnd7qDO2ezsbA4ePJhU6nW27lO/309HRwcjIyNUV1cHGh1mizBnqxuKohhRSNnhcGCz2ejTZeJYtGGaPA/m0NFaGLqLqFjcrOc1NCMNPj8oBrozNjKVKTNWJDNeJDNWJAX+32dUWNSm5/qfZlI4ICJE9LkIxni7zPNfctKwysQnvgMPuWeIsUeS6GHmukarTb7mcjMly9zi87BTN6OQM7pAZsdmD2e8MnP/iopAZauetFyBtBGBpre9lBp17NFoPRx56E0+bD2IuPksDBs2YDKZcLlcrFy5MqjepkqKqd2e8Ro4azHffofH27k+VkOPet1CG3q0TT3qtZjvcYxI9VCr1ZrS9OmJjpOeFGdDJI9BVRopGcQbnaqaoWlpaWERaTIplNnIWEtQVVVVFBYWBvaRbMNMtM8pikJfXx9dXV0Ra5TR7KMijVjEex1EUQzyTRxudfHqT8YZ25/Yi01Fmw5BDu4ePVS6npfOWoo7PVzYWkXPUoltF3i46g/pNL4V3h2sYiJfZu8mD1OZCiu79pHdvpOc32+ATy6LeVzCsculbfzJ9gm0P+JBOpJF4yIdO87zBv7t5Y+6aXzThNEbfP1cEwpP3jDdQJK5UYBjpLjxwAFu/9bnMHl9cJcR13PPoqxcGWSpFlpvkyQpUG+LZOAcS2lGfeGZD0iSNO8ycPEeu9FoJC8vj7y8vMDPYjX0mEwmJEmaN3KMtF2bzXZCKc7MN056UowWKaoSYpOTk2HjFXPVP40FrfFuJM1QlWhS5UyvdZWI1kSTbJNOJFIcHx+npaUlptRd6OfmQoahcNslXv7RCF0v60FOPNJPc4qU9OgYqJAS7h6VdfD6Je4wUlRVXVoKNvDilU4aW3cxKmRz959+PN1V+oaB3afF7iqN1AVbvc3A7t97gSw+kC2z90wvvmO7dmQr9FVKLG6JvgRUtM7cY5v37Mbg9yIiI3m86N58E2XFipjfg06ni2ngrFWaUYlSjY4URZnXSPH91EadDbEaekZGRvB6vRw4cCBIoSdVDT2RDIydTmdcesL/V3DSk2Io/H4/XV1dgfGKurq6iLOGyQ7gR4NWx7O6ujqq/59KUMmQYihGR0dpbW0lLy8vZrNQKiJFp9NJc/N0pW42QXLteIVKiKHzhslgZGSEF77qwtY0uxi63gw3vphH9zYv//hG8KD64hY9AxXJdfKaXcHHr1V12TCu4xNfVdBLErIgIMoKeiX+rtJQrNhtwG1WmMqSyR4XKerT0adxAJHF2GWDnDERy4SAPVcJU8D5ycpVXLhtO6Uvv4R4xRVxN8HMpjQzNjZGZ2cnLpcLg8GA3+9PeQfniU6KkaBeN7/fj9/vp6amJui6paqhJ5pt1L/a9ZoLTpHiMYR6KsZqZklFpKgu8loLqWgkHGnfc1kg7HY7zc3NGAwGGhoaZrWFSTZS1Ol0uFwujhw5wsTERNxmv2otMua8IdPXcKJTovnvHo4878GUKbL0PCPrbwo2WJ2amqKlpQW93oCtaUHYdiLB74apUZkl55q49q96nvysFdf49AKx7nUj28/1BNzlVZgUhUKbQFabgbwRkfwhkX98zBWkg3rlqmkyUKPDLF9/QNVFL02/QOiOCYbLOhGfIoQp3qiobNYzUShhzYtMbg/cNiPsmj8kMlYc/GJTMBT7xWrzOe9R9thr/HJjfZjYAFMuvnZsPlL54x9xPPss+o0bY24vGtTaY0ZGBiUlJQB0d3ejKApmszlowdea6VosFkwmU8IvSv+KpKhC65ChvW4q5trQEy0te6IN2M8nTpEiMDg4SHt7e5h2ZzSkwlPR5/MxPDxMV1cXCxcujLujdC7zhrIsc/DgwVnVYUKRTKQoyzLj4+MMDw9TU1MTl18kEEibdXd3U1BQQHZ2dhBpK4pC1zs+tt3jZPxo8HVwjkrsvt/FwjV6ytcbA44doz0OJl9ZzOAeIAFf+Re/7eCK32eRU67jsp9n8dTnbfimFEq79HzpvywcXeYnV6fjmluzqLYYKRIExgemeO1+JxNt4LXrefUaJ15NY82FeWYWfngvtT84pvmJDr9Oh6CApJ+2eJJlGb/OwB/O+SZZ7gneXbSB3KbVfKh5+tjzhnWcttOAqAh4TArf/e0k0ixP8lixHKR8s6eunqyJ6N/HQuce1tx7DesUH1fcb+CC39zL22tm0sW3339fkEj5PX97Hlt1DV/PzUm641QLRVEwmUwUFxdTXFwc+Jnb7Q7U21Q9X6PRGFjsLRbLrJHRfDfxzCeBzJYlmmtDTygpzvf5nIg46Umxo6ODqampqPJokWAwGJIW9VYUBUmS2LVrFwUFBQlbSCUz0uH3++ns7MTpdFJdXU1RUVFCN3qiRDwyMhLwbiwtLaWsrGzWz2jrhosWLcJqtQa5iBuNRrKysmj6XQ49b8Re0J69eaZT05STj2eygETIMODcsOt0dv7uTE7/YgaFtXouucvCs1+yIftgYbeehd16LvpRJktyp+8bt91P+65xstfZWXnVQnoLTdiNM9c5U1aQd7QgHnkanXKsPifCfR+5gu4FCwLR4OY9u8l/92yU/vW4gZJxKNkffpxLzzcycMBHRYue9uWxX9JClW+23HEfAmdH/f3yqW2Ifh8CEgYJvvvf7/E/12xk2/keJgvksHTq39esZfuElT6/xO+K4+/mjYZINUWtm4N2wddGRoODg7jdbvR6fRBRaiOj+SbFE1H3NN6GHp/Ph9FoRBRFBgcHKSwsTJkX6gsvvMCtt96KJEl85jOf4fbbb0/JdlONk54UlyxZkjDJJBspqrN/Xq+XFStWBLlkx4tERzrULs/S0lIyMzODukrjhSiKcb0EOBwOmpub0ev1rF69mqmpKcbGxmY9xtC6odFopLCwMGzhG+62zUqIofBMJnau2hqfJBh47DeP0LPuXMrXGyhda+CiH1p49U4HrnGFNdeZWbzZiKIodLcM8eKXFbxjFsBCG/DqZZOwZWbblXuMbPtjOc+tPJcfm/6E4Rih/PnSy4LqheLwej7y2uyNDb27fHhsCrXvGWYlxav2ByvfXLB9N74YpKjqkqrOGEO6Mzj3OTPnPG/i6Jl+Xj99bUTt1uccU3hzsxFFMUAOyZBEIuRiMpkwmUxBz1O0yMhiseD1esnMzJwXAjseYt3J2M5FQqSGHtXRQqfT8Ze//IW33nqLwcFBtmzZwurVq1m9ejXnn39+wscgSRJf/OIXeemllygrK2PdunVcfvnlLE/SBWU+cdKTYjKpgURJMVRAu6enJ+l26nijtrGxMVpaWsjNzQ1Eo8PDw0m1o+t0Otzu6C7rXq+XtrY27HY7tbW1gYfM5XLNad5QC5PJROniAiwLJrEPpNajUZ8G/mN9U1rnBtU78Nmb13LDP3NJyxWpPMvIJ5/OBQUMaULgRcC6NxfvWHA6umVV8D0yskDiBz+14jMtY8+Z4YRicsHZ/zBz/tORu2K13oP96Wvx2Kaj32X7DfzjY5Ebv77ry2PtViNb/afjNfwmENn94+y1lKb7qDugjzg32a8xTlb3B9PzjNVvGah+y8BI8Vm8fd56dqyaGfC3KwouUSQDAvep+l/tdzzbPThX5ZZokZHdbqerq4vx8XFGRkYQBCGo6zWayHe8OB5eirP1AMx1+zk5OeTn53PXXXfx7rvv8qtf/YrvfOc77Nu3j1dffZVzzz034e3u3LmTqqoqlixZAsDVV1/NM888c4oU/68gXmLy+XwcPXo0zDVjLjXJ2dKn6iKt0+mor68PaqVWjzvRt7xoNUVZlunq6qK/v58lS5awbNmyoIUs1pxiMiMWol7gsl9ksf9/XRx+OnXekBu/kI4oCux5wEWPMzhC6sk4HYA9D7jY9OXphgaDWcA64KH17W6kzElK0mrY/vvw73O8MPh7GiqbuRbacY40h8Cmf5o480UT6c7IC/Kqqm2c/9yn0EnHIliNrFxJr47znjGz9VI3GZLApc4MztptQn5ewTmq0ISHAx9tCInsToPNDpY06bn40TQWtYcvBapxcjQUDum47JF0DjX6grRTRxDIM003gqkZAPX7jkSU6h8tGc1HFKfX68nNzWV8fDyw8EuSFHEmUNuUkojowL+awXAoQmuKk5OT5ObmUldXR11dHddcc01S2+3r66O8vDzw97KyMnbs2DHn450PnPSkmMzb6Gyf0ZJFRUVFWJPJXEXBI31WKxKujdZCP5tsF6n2c4qiMDw8zNGjRykuLg4y+w39XKrnDXMrdJz7n5ls/o8MBvoHOfzWABZKMJODewKsvRKtL3pn35AG3ds9XPbTHJZ/2MThZ8/mbz//CwW9bwdFSKMt/sA5bH94kL2/NICcB+ScFQHAAAAgAElEQVQBkb/LVbuMvHlRdPLOsAmc/Xczp79iwuyOfR3St+1BJwVHsFrC+tATaVzwtBmdHwQEHAS/jNTvMPKbDzQEiFjbdHPPHfV80JXO9UcsZBwRaH/di28qgRqsUc84M9e8z++n9hgpqsSmvT+0RKn+F4KJMtXen1poCVen04WpH0UTHUhPTw+qU0YTHTgRa4rxIpaX4smCk54UUwlFURgYGKCjo4MFCxbEJItkH3qdTodHo0epta2aTSQ8FUP4NpuNI0eOBMxGYzUnqaMVqZ43VM2NLRYL51yxKhD59u/30XJH4hFk99sST3+/iYoPechek0XDX+rpf20D/T+fXuhFHaz4qHlasLyphXfvXRyXwvclj6RR1KfjtcvdQZFU1rjA5ufNrN8arigTDaE1PjWC1ULvj76t4n4dX/5vC29+xovfvYOXQ+ymXqyv59XVTm7YnMXi/SK6qdjHlZYrULJKz+rr0thZ4OOgfYYU+2d54ZuNKMfHx5mcnKS0tDRQy44WUSaD2YgrmuiAOhM4MjJCe3t7YHheS5THm7Tme/s2my0lXoqlpaX09PQE/t7b20tpaemctzsfOOlJcS6Ls7buodbwcnJyWLduXcw5wrmmT51OJ4qi0N/fT2dnZ9wWS3OJFL1eL++99x5ut5u6urqoxsJaaH0Y46kbzgY1GvZ4PCxbtizgoiF5FXb8zsm+B6PXPWdD39NF9D0NJWtl0ko96ItGWfYNJwZdGmajmT5rH4d+kU7fc0vj3qaoCGx43UTjW0YObPDSUeOnvEPP6neMAY/GqJ/Vg6y5RaLV+GaD3qyQX+/FssJKRq2NS/PMWB7dg8nnQxdiN+UH7rXaKL5F5NZvZaE/Jhh+9SPZZC3UYUiLfMwLh4OXkYEk7m1RFPF6vTQ3NyPLMmvXrsVoNM4aUSZDlMlEc+IxwfPMzMyAa3wk0QG3240oioGmHnWWMlU4HqSr3f5cHDK0WLduHa2trXR0dFBaWsqjjz7Kww8/POftzgdOelJMFirBqA4PkWp40aDX64OivUT3a7fb2bFjB9nZ2bMScKRjTgSSJNHX18fo6CirVq2Ke5xDbam32+3s37+frKwssrOzycrKSrimqaajh4aGWLJkSVAHraIovPgtB+1bE0uZRsPgHhH2pAGhzQzJLww6SWDNOybWvDOzODbemMaGm9J58rNWBg6Ek4jshwUXjzLw95mOytlqfKG45G4LZesN6E0CsDCwiHvP2oTwu98i+3z4IwgDDJXJNK2etpdaUK8nf2nsZWJhSOTSlyApql3SPT09LF26NMx0NlpEqU3FqxkJURRnJcpUpTgjiQ4MDAwwNTVFRkYGVquV3t5evF5vkMpMsqIDMP+kGCohabPZqKysnPN29Xo999xzDxdeeCGSJHHDDTcEXG1ONJwixSQhCAIHDx7E6/VSU1OTUN492UjR4XDQ1taG2+2msbExplRaJCRCitpUcFFREbm5uYEh6tk+py5Woihyxhln4Ha7sdlsjI2NBcTU09PTAymqrKysiCkh1SWkvb2dkpIS1q9fH7aYua1KygjxuOJYyW6iM/r3UbWinKn9Hmz9iXfbfvDOTCrPCn5ZCiziF1/M6B+ew/7wVvaPrmPxWytpXuhionBmP7vO9rJqt5HMSh9erzfmi5dLCT6+RDxwHQ4HR44cwWKxsG7dullTg8nWKLVEOZ91P1WFp6ioKEDusVRmtESZlpY2K1FG0iZNNbTHkKpIEeDiiy/m4osvTsm25hMnPSkm+rbm9/tpb2/HbrdTVVXFokWLknKsSCRi0448lJeXMzo6mjAhJrLfiYmJQM1u3bp1ABw4cCDmZ2LVDdPT00lPTw+8TasRi9VqDTTsSJIUMAbOyspCEASOHj2KyWRi9erVUVNQhjQBQzr4nIlcifcX5RsM1H/cjGtCxm2NziBv/iQ5fd1ll5uouiCcxBzDEq0veWl5wcNocy1QC8Dq7VA4KPKL782IHjTX+7DmyJgWOHnvvS68Xi9paWkB4+asrKzAd7LLFZy2rjfPni6UJImOjg4mJiaora2NKx0fDYkQpc/nw+12B0Q0UlGj1EIrw6YilsqMSpTDw8O4XK6YogOh53s8YLfbTzXanIyIx1NRlmV6enro6emhoqKC4uJisrOz53XOUauLqo48uN1uBgcHE94nzE6KLpeL5uZmJElixYoVgZqd3+9P2bwhBKedFi5cCMy4AExMTNDU1ITL5QoMZQ8PD5OVlRWxNV5vErj8niz+eoMtrmtQcaaBrreDhQgu/H4mi87ScXhbD/t+acTZk7h7RiIoOc+K0+/HcTT182YFtTpqLzbh94Dh2Gl0bfOy70EXfXv8UYV9Sjt1LOjWMbBo+v5QRNh7ppcvnF9GTkUFiqLgcrmw2+1MTk7S3d09LWKfnsEb5uCSQeMsPqNjY2O0traycOFCGhsb50VGLBJRTk5O0tzcTH5+PkajMSyiVD83F6KUJCnuGqLRaCQ/Pz9ID9jn8wWIMlSOTW3kma9IN9IamMpI8V8Fp0hxFiiKwtDQEO3t7RQVFbFx48ZAs0uyUm+zkaI2dblw4cKgLta5dq5GOmY1+h0dHaWmpiZMaSfaPlNp6QTTi1ZfXx+VlZWUlJSgKMqMMXBfH3b7dCSjRioqUZasNHDNo9k8+yU7U8ORyXtBvZ4Lv28ho1Dk9R87OPjETE2355evYfrLC+SffR4f+Z+LePrf7ThHE8gBJoi+101kLRum9S0JKEnptkebJZ7+vA1RBwV1enIXiTT/I3Z6WWeCyk0m/g2FXzDjCfnuxV6yF00vvoIgBCL+4uJiFEXhtxNW/mNkFO0dZVYUPE2HOapJjask4fF4aGlpQZZlGhoakjLpTgaSJAUyLStXrgzKsoRGlOr/+/3+ID/ReIlyroRlMBhiyrF5vV727NmDIAhBs5QZGRlzTqtGqleebF6KcIoUgeiRoppGzMzMDBs/UG1tkkEsUlT9BrOysiI20cx1xlGrTDOb2a+KUKJLNRmOjY3R1tYW0IJVH0xBEMJa49Vha5vNRk9PD3a7PWCrc+b/y2L798zYesL3UXGmkYzC6XNb8VFzgBQXOvfwwVeunZZ1e+RBHqt8BGcCzSzJoP8dHWNNeXgm5rFhQoLhQ36GD0X+d0GE0kYDNR8ysnSzEWOmyFp/Or85ag2QXF+WxM8mJrkuO4t8zWLpkGVuHhzmcbsjbLvX5WSzsqIMm83G5OQkPT09gaYyj8dDWVkZZWVlx40QR0dHaWtro6ysjJqamqipyNDUKxA2TqS+GCqKgk6nC6pVqpiPRhhVji07O5vBwUHWrVsXJjowNTUVt4FzNEQzGD6VPj0FHA4HLS0tAEFpRC3mSk6hkZdqbwSx/QbjSfXG2q/6wKsjJLN5KWqR6nlD1VBZFEVOO+20uOSrog1b2+12bDYbK74yyOC3Oyhq3xM0uuBzzUSQ2ZUK2VU+rG2GiLJuiXR4JgK9GUx5MlP9YsoJMVQCLhrSS2QaPpZBzQfNgZcEFYV6HRdnZvCMY8Zy6lsjY3x3ZIwPZKSzJctCtdHAZweGafKGR5+35uZwR2E+hmPC3cXFxTgcDg4fPkxmZia5ublMTU1x+PBhPB4PZrM5rEaZqlSqdrwjVk06EiJJ0anPjfoyqK1TwgxRqnXK+cBsogOzGTirRBntWY/UxKMKhJ9MOEWKGmhVYWpqamKmDeYyVqGa6EKwuXBNTU1Q2iQS5vLAiaKI2+1m7969CIIQ9wiJilTNG/r9/kCTRXV19ZzTMzqdLiBsXDmwA9P+61Hc3iA5tHG5g127XIiiyNTUFIsvXsz+XxjiGoqfCzKLRSrPMlDQ4MWW1oYwUsyeH4a/ZCUKQS+z7Ks9TO4qgL8dCRYx10jAhaL6vHRWXxv9O78uOyuIFGFar+cfU07+MRW5m8kiCvyupJgPW2bOS5Ik2tvbmZycZNmyZWFOC6oNlPoy09vbGyBKbXo8UaJUSw9dXV0RxzuSRTxEabVamZycpKysDJ/PFxRRhn42GcwWhUYycA4VHVC7v9PS0oIaeoxGY1iTULIv3//qOEWKTN9sR48eZXh4mCVLlsRUhVGh1+uZmpqK+TuxoCgKHR0d9Pf3x2UuPFf4fD56e3sZGxujoaFhVvLVHqcsy2RkZLB9+3YyMjICM4cWiyWhVJF2wSovL6eqqirl5yy+8Qb4vIhMR+K15h0UfeJMll9TRktrM0ajkdLSUuxZ4yy42Mn47tN4yvK/lNl20J9zBpNZa0lneqxANRWOhexykcJaPW0vB0dOyy4zcdrVZjLL/bS1teH1+aivXYlBSOO9X07gjSGjVlinY/Un0nj9x1MB0e9QLL8sDbG9hsHXPaxLINrNLo39fV2Ykc6nsrN4wBpf49Jyo5GHS0uo1kQTasqytLQ0aiON1gYqdHTBZrMFEaXJZAoa34lGlE6nkyNHjpCWlkZjY2PK3CSiQSU59Vl2OByBF81IEaV2ljIZB5FkxjGiiQ6oTVNaA2e1djo6OoogCIG06Vyf0ccff5zvfve7NDU1sXPnThobZ+Zif/CDH3Dfffeh0+n4xS9+wYUXXjinfaUCp0iRackhg8EQtaYWCcmmTxVFYXBwMFADiCYFlyqoXbO9vb2BkYh4CDG0brhixQoURWFqagqbzcbAwAAtLS0oihIYpcjOzo7qMjA5OUlLSwvZ2dnzumDJZ58NRiOK14tgNFLz080053XS0ekJS4U3NExH6jZbI1ZrNQtsNnLdzccGrbPY9d0cJlpjLwjWHhlrT3gqcfG5Bpymftr294cJDtRebOK9x4PHGGoumq7tvfK9KSa7JY485+G0j5nZ9fvIYxmHntI0CsUZ7S7Y4GfBWW4kyRD1nhMEgV+XFHFbXg5P2Bz8xW6nxRu5oezqLAu/KC4k49j37fF4aG5uBkiqkUY7uhCNKPv6+nC73UFEmZmZydDQEMPDw1F1f+cLag9AaM0y9BkIJUmVINXfj8dBJNK4RzIIbZqC6evc29uLzWbDarVy1113sW3bNlwuF9/85jdZs2YNa9asobKyMmGSXLlyJU8++SQ33XRT0M8PHz7Mo48+yqFDh+jv7+eCCy4ICKG8nzhFisDixYsTJrhkjIa183/p6eksXrw46bew2ax1FEVhdHSU1tZWCgsL2bBhQ6DOMtt2o9UNVZsdi8US0C2UZTmQAuvp6cHhcAR+Lzs7G5PJRF9fX2DMI5n5ykQgb9iA5/nnEV5/nYHaWtoNBpYc82aMdL2MRiMFBQVBHbdut5vBNtushBgLL/33JOf9yse6deuCHnLXhIzfEx79law0sOt+F17H9L91b/fRvT2++yuWBJwpS6DmIiMV5/sR850Mj49xtGumi1dNVVoslqDFuMpo5PaCPL6Rn8sBj5fHbXYetzvo8/vJEAS+X1TAjdlZgRp3b28vfX19LF26NGgWb66YjShHR0c5dOhQICIaHR3F4/GQlZWF2Wyet+yLz+cL2MHF8wKgjQ5VhHa+zuYgMp9qNuo+cnJyKC0t5Te/+Q2Dg4PceOONnH766ezdu5cHH3yQ3/72t4GX63ixbNmyiD9/5plnuPrqqzGZTCxevJiqqip27tzJ6aentoSRKE6RYpJIJFJUG0pkWQ5EKzt27Ej6zU+1j4r2WbvdTnPzdKpwzZo1gQdWkqSY4xzqg6mSYTxRsyiKERtfJicn6erqwmq1YjAYMJlM9Pb2BiLKeNQ7ksVIVRVtQEFBAesqKxNeSMxmM0ULjehME0hJOlT5bXqWLF6CqJs+R8mnsO0eJ4eedOOPsM3MYhG9MfnrEUkCbvE5Bj74PQt6s7rdmQhK28Xb29sbsEwKHXcRRZEGs4kGs4nvFebT6fOTpxPJOXZN7XY7R44cCWj+Ho+3fEEQMBgMTE5O4nQ6WbduHZmZmUE1yv7+/kBEqT2nVBDlyMgIbW1tgdGhZLeXqDqP2+0OOIikWnQAptOzWnK32+0UFRVx6aWXcumll6Z0XzBtJ7Vx48bA38vKyujr60v5fhLFKVJMEvGQotpEMzk5SU1NTdCQrvr5ZEhRtY8K/ay2Uaiuri5s6PZ4zBuqEWp7ezsLFy6koaEBURTx+XyBFNjw8DBOpxOj0RikiTrX7kNVhxaIu5s1GtLzRM75egavfm/2unGkzs/lHwMZCfHYI/bKHY6YllZTozL115gZfC98zCEZZBaJfOgHFkR9dMeUaF28drs9KOrXKg1VHiNKbSNNXV1dWCPNfEId/i8tLaW6ujpwz0RSjdGmXlWiVO+7RInS6/Vy5MgRgIBgeaoRiSglSaK/v5/e3l5qampSJoweikgOGfEqDV1wwQURRUXuvPNOPvzhDyd9TO8HTpEiyRWSY5GiLMt0d3cHBtEjNdGkcqRDlmU6OzsZGBiIaR8VyRcxlfOGdrudlpYW0tLSwhYNg8EQpt4RqVaktumrRBnPwqPaZ42OjlJVVRW0j2Qg7tiB+MYbZJk2AsHO4IvPNjDc5GdqZDrNudC5J6zzM+8L9RSfN8m+fdP2Qt7OPFpfjH1Mr/9ois23Z1C5yUDnW7OnTbPLRfLWOuh8Ng0lgo3VGV9Kj0qI0aDt4lWhHXfp7u7G4XAgSRI+n4/8/Hxqa2vnPSWuwuv10tLSgiRJcdcsTSYThcfS5yq0993AwAAulytAlGpUqc1kqH0AnZ2dKe1ojQdut5vDhw+TlpYWpA0bSxgdkiPKSAbD8dZnX3755UROCzhx7aROkWKSiEQeWvWbWOa7MDdlGpVQ1f0dPXo0pn9j6DGnet5Q1WZ1uVzU1NTEHTWELlhqm77NZmNiYoKurmnNzfT09ABJWiyWQJOOVjB8wYIFrFu3bs4pJXHHDkyXXAJeL4sFAwvLZ8YbVn/SzBm3TBPAw1dPMNEuh805XrB5L5YvfBAoZmpUZu+fnbz3aBw5WAW2/mCKhk+Y0e/x4Y8he5pRorD0q0coW7CIqjWZvPSd4Gh2zfVpVH8wNXZFWqJ0u92B5qoFCxbgdrsDRCmKYhCpZGRkpCy9l+oxi1hEabfbGRwcDBBleno6Vqs14B96vGb2VGu4np6eiKNas6Vek3EQCe1une/B/csvv5yPf/zj3HbbbfT399Pa2sr69evnbX/x4hQpMveWY5jRVYykfhMJc40UbTYbzc3NpKen09jYmNBwcqJ1w2hQO1v7+6c7LOO1lYoGbZu+tivO6XRis9kYGRkJDCSbTCacTifp6enU19fPKVWqhfjGG+D1IkgSIsHjDT07fOx+9zWyO7aSNngmE+lrA52fAILRiPjBsxk+7OfgX920vOhJuCa5/yE3afnEJMXFW6xk/MXF+Iv30S5uhJBaYqortdpGmqqqqjAZQJheUNWIsrOzk6mpqYCnoJqmzMjISPj+cLlcNDU1zfuYRaQXtK6uLnp7e8nJyUGSJPbs2YPBYAicTzTB7rlCPWf12Y63xDJXB5HQPgWr1TonoXYVTz31FLfccgsjIyNccsklNDQ08M9//pMVK1awZcsWli9fjl6v51e/+tX73nkKICQ4oPl/cppTURS8ERQ6ZsM777xDQ0NDIKWTSJTU1taGxWKJy45JC7fbze7duxFFkVWrVsW9P/WBePfdd3E6nWRkZARqSpGEtmeDOotWVFRERUXFcbuZ/f7pub/JyUkKCgoCAsrqaIgaUUYbDZkN2kjRLwcPwkdKlfanr41bTSZViHYcKi76USZLzk1NpKhtpFmyZElC37Pf7w9EXzabLSGiVEsQg4OD1NbWHlf9TafTSVNTE5mZmVRVVQWd8/QIz8w5OZ3OAFGq55UsUaqyi2rtMN5Z4kQRiSi9Xi/vvvsuq1evRqfTIYoi99xzD2VlZVx//fXzchzvA+L6Uk5FiknC5/Ph8XjYv39/RBHt2aB2kMYL1WpHdYwoLCyMixBD64arVq0KzBtardaA0LY6RqHW86ItVqocnV6vP66iztoU2qJFi6itrQ3TnIykiaptqIgnUlFHOtxPbuX551cHkU00SbhEzX/nimjHkV+lY811aSkhRFXQwmq1Jt1Io9frw8St1ZcYm81Ge3t7wAVC+z35/f6Am0UkD835ghodDg0NRZ13jDTCo1pAaZvIVK9E9ZxmI0ptdDjfXbzaiFJbL62trcVoNAaetYceeojPfvaz83YcJypORYrHEK9km3YYXlEUGhsbkyKGnp4eFEVh0aJFMX9PrS10dnZSVlZGeXk5PT096HQ6ysrKYn4ukbqh2lBhtVoDb/V6vT7ooe7v78dmsyVsqjxXWK3WgEj6kiVL4k6haVN6VqsVp9MZdE6hDRWhOPiV1/A88wbd5tMDEeGWzmsQFR/yLHJq84nQ43jpw38g45OnUVafQXZ27HOKB2qauqysjNLS0nlVWoIZopycnGRgYACPx0NGRgY5OTlxk8pc4XA4aGpqIjc3lyVLlsyZiLVEqY0oQ4kSCESHxzsi9nq9NDU1odfrqampwWAwoCgKTz75JD/5yU+48847ufTSS+f9+z+OiOtETpHiMXi93phaf4qiBAxxi4qKqKys5ODBg1RXVyfVfad2vS1ZsiTq70xMTNDc3Ex2djZLly4NFPl7e3uRJImKioqox6qtGyZ7U/t8PqxWK729vYyPj6PX64NSlFpboPmA1+ultbUVt9tNbW1tRGH2RKEdDbFarQHfxtAW/amXXyZvyxZEnw/FaGLwnmeYLFuHsGMH7r++zmH7+iBCzF2sQ/Iq2Pqi+06mEgudezizfg+Fnz8X/8a1gXNSF+Bkxg7cbjfNzc2IokhNTc28freh0I5ZlJWVBVKv6h/VgDfeF5p4IcsyHR0djI2NRdRoTSW0957dbsfhcODxeEhLS6O8vDxmhibVUNeyqqqqQB11bGyMr3zlK+h0Ou655545d3GfgDhFiokgFilardZAU0tVVVUgMjx48GDgZk4UIyMjTExMUFNTE/Zv2mH/SC3vqkzc0qVLg36e6hGLiYkJWltbyc3NZfHixej1+kB3qBpRqt2h2jGKuUpRybIcaOxIRQPPbFA7D1VBZ93OndQ88ggF+/YhKAqKTofv29/G/7WvAaDICv/7b5NYe2YIsPGGNDZ8Ph3Zr/CbM8bn5TgNGQI+jWbqhs+n0XhDZHFvtfalJZVoGqKKogQapqI10swX1DELv99PXV1dzKxLaPTlcrnC6nmJEKXVauXIkSMUFxezaNGi45qm1TYuCYIQVKNU08naTt5U3f8+ny/gHlJXVxdIl77wwgvccccd/397Zx7eVJm+/ztp0r103+jepiurXcHBFceVwRnwQtwd5QsqHVAGhBkGAREEHBQQB1FgQGZcEBd+KiLuyFBaVrF2Sze6pXuTpmnWk/P7o76HkzRps54CPZ/r4hrJQHNCkvO8z3bfWLlyJebMmXMtZYds+KBoD3q9fpC7PFkG1+v1SE9PH3SKrKioQHh4uEMnqp6eHshkMmRlXd6DMxgMqKmpQXd396BlfzbmAdXVwVCtVkMqlcJoNCItLW1IJw32dKhCoYBSqQRFUSZ6qOYSYkPR3d0NqVSKsLAwJDqgRuMoJBD3ff01sp9/HgKdDjAaQUMISuCJo/e8jYCn0hGRNCCI3nDMC/979bJ+qae/AFHjRejrMKK7xv5VG6Enjehbe3HDgggcW2qw6WcEJQjx0Ie2l9vIgYb80mq1zEFnzJgxSEtLc9kU73Cwe1nJycl2D5wRLAX/4bJk0i/t7e1FZmYmZ3uWwMB3i1hpmQ/xENh9V3agdHaSl8g+JiUlMVJtvb29WLFiBXp6erBr1y67JdyuMvhBG0dhByd2ecEcR/RPCexBG3JybGhoQHx8vEUzVDZkx9HV+4aOLMELBAL4+fnBz8+PUeE3Go3M0IutgzwkENM07bQajb0QsfKQkBBM6OqCQK+HwGiEEUJc8puGkxHPoaUuB8J/AMl3U4i+vQOaGCU8fOJBqQduaro+2matUhOENELzezHlKT8kZCbhwrsadNfY1t9WNF1+322BrSFKPuNyuRzx8fHQ6XQoKyszyfzJL1evQZChEm9vb6fXLKwNvpgv5xO5N6FQiLa2NsTGxpqo4bgbdnaYkZExZE9eLBbbPKBk3qO0tn9YVVUFnU6H7OxseHl5gaZpHD9+HCtWrMCzzz6Lxx57jLNM+UqHzxR/Q6/Xw2AwoKmpCY2NjYiPj0dMTMyQH5T6+nqIRKIhB16s0d/fj8rKSsTFxUEqlSI0NBTJyck2lR57e3tx6dIlZGZmuiQYsk/tZLjC1V8Qa4M8/v7+0Ol0jIelK0t3RJnGeOONMBYUDHpcP2YMuqVS9EyahLGzZw8san9wEuHz/gCB0XSYhr120Rmei+se9kZfpxG/fuSgOCqAwPEqZD6ow6SbEwcyNrkR/5kth1Zp29cs+WZP3LnJ3+73nmh3xsXFDRqkIbZC7OxLr9czlmHmAgr2QIbUZDIZ50MlxDhcrVbD29sbOp3OZksqZyErHgEBAUhJSXFZ9YMdKJVK5aCVl4CAAKYvHx8fj+joaAgEAvT39+OFF16AVCrF7t27rc4mXIPw5VN7aG1tRXl5OcLDw5n+2XA0NzdDr9cjMTHR7ufr6enBuXPnEBoaOmyJkg25aZ0+fZrp4QUFBTm0awhcnuwMCAhAcnIyp4odMpkMtbW1zGQhEXB2xSAPe98Qnp7QfvEFjAUFlx/XaAZME4VCwMsL6v/3OT47mIWGk/pBe4fD7QUClvVPrRF3ez/8suTI/0OiSUn+py0qXPzgcklW5A1ETRBD7DPQTxT7CCD2Hfjf4EQPJN3oCZGX7TdxRwdp2CVy8os4urODylCfP6VSifLycoSGhiIpKYnTrMSagDdbQLy3t3eQJZWzgZLdqx0uO3QVZEBJoVCgpaWFmeStqqpCT08PQkND8cYbb+DJJ5/EwoULR1t2yJdP7cHDw8PEUcIWRCIR1OohpEcsQETCFQoFPD09MXnyZJv+HrtvKBaLMa/9DMMAACAASURBVHXqVPT391vcNSRBZai+AxEP12q1yMzMdMlkp62oVCpUVg74Fubl5TE3Z7YlkFwuR0NDwyCZN1sHedjKNLRON5AxFhRAd+wYvLTagQEaAAKjEbROhzMPfYWG8FQAgx0nrO0FEmwJmgTvsWoE3V4HX19fyGQyqFQqjBkzBtoOT5QeMvVYzJvni+xHnS8js2/OqampdvfArZXIibdma2sr04Nmi4eTgE+Ew7Oysjj9nBGrNJqmLUq0mQuIWzM5Jnq89hzU+vv7UVZWhsDAQM7cQ4CBexIpEcfHxyM2NhYURcFgMOCNN95AaWkpvL298f7776OiogILFy40mWvg4YMiQ0hIiN2ya/ZItbFFwpOSkpCRkYGioqJh/95QfUPih0egKIr5QtfW1kKlUg1yovD09GSUQszNb92NwWBgbpCWdh2teeeRLIWMkZNBHrYe6iCvOpbZMDw9oZkyBRWlpfBOTMQE8rjRCFoghMFo3ZgXGN7Ed7igyUbT4gPRzzkQR9LQBWnR4a9Ei2crynYGw0hdDhj+UQJMvN95YYTe3l5UVFQgJCTEpTdnoVBo0VuT9JJbWlrQ09MDjUaDgIAAREdHM4c6d2cnjgp4D+fdqFAo0NjYOGSgJAcQmUxm0anGnRiNRubAPXHiRKb6VFZWhr///e+47777cOjQIYhEIvT19eHChQucXt/VAl8+/Q1ymrIH8iUZP3681T9DRKuJJFpSUhJzYzp58iSuv/76If+us/uGOp0OCoUCCoUCnZ2dUKlU8PHxQVRUFBNU3KUnSTBXoxk7dqxTgZh981UoFCbGxuxBHo+Skstmw5GRzCHAo6QEwuPHQQWG4NyWZtSLpwxb8hyqPMpepqcFHrg0bi60f3oQ33wzDtpeetifIxACtNl6Y+q8DgRM6IGPj4/JzdfW8jYZpFEqlcjIyOA8QyNrFmlpaYP280hAZU9SuipQajQaVFRUQCwWMwvproYtXE9ek1arhaenJ/r7+zmf5AUGDj/l5eWIiopCfHw8BAIB9Ho9XnvtNRw9ehRvvfUWJk6cyNn1XKHwPUV7MBqNdk+SqlQqSKVSqyVQohvp7e2N1NTUQaVZa0HR1SsWZMjAy8sLKSkpMBqNJruGzqxQDAcRLidDBu4KwOwsmajXAANl4sDAQEgkkkHlZMpAY9e07kEByR6CkzwQGCtEnP4s0mQfIfib/wIUBVosxvePvIFzJ3/P/Flby6yR40WYvWdAiFmj0TDvU29vLwwGw6DpUPNyMsmoXXEAsQfzNQtr+6VsO6re3l709fU5vXLA1g11pETsDEQeTiaTISoqiunr6XQ6hw81tkLEB7q7u5GVlcWsl1RUVKCwsBDTp0/HqlWrOJsVuMLhe4ruxtpKhlarhVQqRX9/P9LT04csUbBH6l0dDPV6PZMtpKWlmVyHr68vs5NE+kPm/UnyRQ4MDLRbZottJ8VFz9LDwwPBwcEIDg5m/v21Wi3i4uKY/ikZzScZMnoDnAqIt7/kj8QbPHF0uRI/nZoIXcdxTNNRA2VUI5DQVI2yoNuhkQ+cJW0ts0579vK/NXENIe8V0a21VE729fVFV1cXvL29ObU5AgbWLCoqKuDl5TXsmoUl30a2go01TVRrn0G2gDeX/Ttg4GBcVlaG4OBgFBQUmBwm2ZO85lZo7L6ro+9TX18fysrKEB4ejpycHMb8+c0338T777+PnTt3XhFWTFcbfFB0AvOeIkVRzIkxJSUFkZGRQwYSoVDI9FhcuW9oNBqZU3NiYuIg8WxL12GtP6lQKFBTU2OxP2lpKIlrNRprz52SkjKoX8ruD7XVKfDTEgBwLCO+5e9+SL7VE18uUzL7iSa9R4hRXJEPje/l4gpjM0UDEInRP2ka/PQCxrAYAHLn+SBqovWAIhAI4O/vD39/f4wdOxbA5WX0pqYm+Pv7o7+/H+fPnx9UTnZHL4+maTQ0NEAmkznl7GBNPJwESiK0zV7MDwgIQHt7O1pbWzmb7iSwxcMzMzMtWiwJBAL4+vqaHEDZgbKrqwt1dXUwGAyDMsqhDhXkudvb202k6S5duoSFCxdi8uTJOHHiBKfl22sJvnz6G87YR02dOhWtra2ora3F2LFjkZCQYNMN6MyZMxg/fjxEIpFLgiEwoF9YXV3NjL678tRM+pMkWGq1WhOJN4qiUFtb65bnHg72Ar4tz13znRZHV/Q59FxeAQI8+U0wvntRhYovTPcU2T1DAIP6h0SzNGrRrczupEFDQ9lmhFeAAL4h9gUu9iAN+3UPV6J0JPs3h7QHiAwgF+83Wczv7OyETCaDUCi0qMfrzoMYOzt0hXi4+cqLUqkcVCYnu6FkqpVYeZGD9f79+/H2229j27ZtuOmmm1z0Sq85+J6iPTgaFH/66Sd4enrC398fqampNpdCaJrGzz//DG9vb4SEhDg98NLf3w+pVAqBQIDU1FROTonk1NvZ2YnGxkaTheihJkNdCSmVEik+W/c9+7uNeOfeHrtNgAHAJwRIv8sbF/5rukIh8ADo39TZhuof3rTcD+NnOzdZ6sggjbnINts1hJ39DxdQyOFHLpdzvs5jScDbknydI2sUw2FLdugqLO2GqtVqUBSF6OhoqFQqJCYmQqfTobCwEAkJCXjllVfcKmh+DcAHRXux1T4KuKyL2t7ejilTptht9kuMPbu7u5kPPdso19aSl8FgQF1dHXp6eiCRSNxmTGoJiqLQ0NCAtrY2RkjafDKU7WtIbryusAEarlRqCxfeVeN/W/udug6Cf7QA9+0JhEAgQM8lCt7b/4mYQy9DSFMwwgMnIv6K4vBChKZ64N43xsAnyPGDgisHaYYSDifvFzugdHd3o6qqCmPHjkVcXBynwtG2CnibT4eaC9c70ssj1lIkI+dy6Z2tlxodHY2+vj7s2rULn3/+OVpbWzFp0iTMmDEDeXl5KCgo4Mzj9CqED4r2YktQJEGoo6MDqampqK6uHtRgt4QtOqVGo5GRQlMoFCaehiRQsp0NyJqDJbkud0LTNDo7O1FTU8OMgA/1+tkZCpkM9fT0ZG665HXZir2lUmsYKRr7Z/Sgv8u5j7XnGCPGLW2FILCPGeQJk0oR9fDDjKKO6uPPoc7Kh2+o4yVytVqNyspKxv/OXYM01jIvrVYLoVDodoslc1wh4O2ofJ3RaMSlS5fQ0dGBjIwMt2aHlq65paUFjY2NJrJ4HR0dWLJkCXx9ffHqq6+ivb0dZ86cwdmzZ/H88887JDs5SuCDor0MZR9FRr5JEIqNjYVQKMTp06cxadKkIW9QzuwbEk9D0svTaDQj5mwADPRTqqqqIBaLkZqa6nBZir0QTW68w1lQOVoqHQp5I4X/zpY7/PdFPsCfdgYiIkvEDPKQ10SfPAn/M2egLiiAx7RpdinysGFrhjozzOIIZM2itrYWYWFhEAgEg1ZDyBqPs5ZhliCZKfFZdOXBbzj5OrFYjKamphGRptNoNIxoempqKjN38MUXX2DdunVYvXo1Zs+efa1aPLkLPijai7WgSL6YQUFBg3btzp8/b/UG7eoVC41Gw6wahIWFMYGFoiiTAQp/f3+Xf4HZZVpLajTOYu0GRaZiNRoN497hahWeix+o8dMW+8uoQhEw47UAxBUMfSAicnzkdRE5NFv6rsTLMzQ0lFMrLcB0zSI1NdXkc89eDbH2uhzV4wUuOztotVpkZGRwdvAj5f+6ujrI5XJ4enoyE7/sjNJd7wN715NtHyeXy7F8+XKoVCrs3LnTYautodBoNLjxxhuh1WphMBhw3333Ye3atairq8PcuXPR1dWFnJwcHDhw4Grde+SDor2YeyoSJwsASEtLs1i2+eWXX5CQkGBSVnF1MCSrHu3t7Uz/jA35IpOMkkwassuutgxQWIL9JeW6TGs0GhnRcJFIBIFA4Jb+JG2k8ekzvWg5Z7uiUajEA7971g9x+Y65RVhT5GH38cgqTEZGBqeef2ytVHsyU/brIlOU7IBia5+cKEAlJCQwzg5cQXb/iJ8nme40f10ABh0AnD2I6nQ6lJeXM+VxsVgMmqbx/fff4+9//zuWLl2Khx9+2G0ZKzno+Pv7Q6/XY9q0adi2bRteffVVzJo1C3PnzsVTTz2FSZMm4emnn3bLNbgZPijaCwmKZOmdZEVDqWOUl5cjMjISISEhLvc3pGka7e3tqKurQ3R0NOLi4mz+QpAdL5KhEMscEiRtmXbt7e1FVVUV/P393apGYwlieaPT6UwycUv9yaEGQ2xF0Uzhg4fk0FtIGP0jhIjIEiE80wORWSKEZ4jgHeh6ay3ymtrb26FUKhmxansmQ53F1WsWQ62GmB9s2ALe6enpLpkYtRWj0ch4idrSM6UoyiRQmksN2itfR4an2P6tfX19WLVqFerr67F7927ExcU5/Tptpb+/H9OmTcPOnTtxzz33oLW1FSKRCEVFRVizZg2++uorzq7FhfBB0V50Oh0uXbqExsZGJCQk2JQVSaVSBAYGIjw83GmdUjZKpRJVVVXw8fGBRCJxulzBnsgjGaU1eTfi5NHf34+0tDROhyrIVGlLS4vNguXsPp5CoXDYWaO93ICiXQooFX0IzxAjdWoIIrM84RfGTS+JlCs9PT2Rmjrg2ME+AGg0GpNVg8DAQJeVsSiKYsrjGRkZbn3PzVdDVCoVsxJFplq5OAAQiK1VeHi4zTvGlnBEvk6v16OyshJGoxEZGRnM+3ny5EksW7YM8+fPx4IFCzjrZ1IUhZycHFRXV2PhwoVYtmwZpkyZgurqagBAY2Mj7rrrLpSWlnJyPS6GD4r2Ul5eDoqibPZTBIC6ujp4eHgwljrOBkQSkFQqFdLT0916c2LLu5H1CYPBAIPBwIy9u6I8aSuumipl9ydJsDQajSbZiXm5S6vVoqqqChRFIT09ndPhJeKg0traOmS50nyQh6wamPsa2jvwMpJrFkTAWygUIiIigum/utJb0xrs7NBdtlaWDgBkolwgEKCjowPJycnM/UOj0eCll17CuXPnsHv3bkgkEpdfky3I5XL86U9/wrp16/D444+PqqDIy7yxSE1NBUVRNv1ZUiYNDAyEVCpFQ0MD/Pz8mPKkvc14doZErKXcfXNiy7v5+flBKpUiIiICISEhUCqVqK6uZsqT7PUJVzfZ2aXScePGOd0/s+b/R07xDQ0NJqd4iqIgl8uRmppqs82Qq2AP0uTn5w+ZEbCtjcigBXuQp729HdXV1TYP8uj1elRVVUGn02HSpEmcHgSGE/A2t2xie2vaKoc2FOzsMDc3122ZmCX5OrVajfLycmg0Gvj6+uLnn3/GnDlzkJKSgnPnzmHOnDn45ptv3DLNaytBQUG45ZZbUFRUBLlcDoPBAJFIhKamJsYq7FqFzxRZ2GIfxe4ZskulpEnNzrqIqPZwQyGdnZ2MtVRCQgKnE4ZkopWiKKSlpVmcomU7NSgUChgMBuYA4Mw0niOlUlfS1dWFyspKiMViiMXiQYvr7jgAEEjf2h2DNOzBK/N+Fzmw9fX1ob6+HklJScNq9LoaIuDt5+cHiURi882fvWvIrmyYZ8pDfRbZrhJcq/EAl7Py+Ph4ZohIr9fjpZdewunTpzFhwgQ0NDSgvr4eBQUFeOuttzi7to6ODojFYgQFBUGtVuP222/H8uXLsX//fsyePZsZtJk4cSKeeeYZzq7LhfDlU3sZLijau29obSiE3HBFIhHq6uogEoksWku5E7KUzFajsRVrBwBy07VFV9NVpVJH0Ov1TBacnp5ucmNk913Ny5Ou2Mcjw1O1tbWcTlcaDAYolUp0dXWhubnZJJvkapCHLR7uKgFva6sh7D4eyZRJdhgRETGs4ISroSgK1dXVUKlUyMrKYr7rZWVlKCwsxJ133omVK1eaZL69vb2cigVcvHgRjz32GCiKgtFoxJw5c/DCCy+gtrYWc+fORXd3N6677jr85z//4XQIyoXwQdFerHkqunLFQqPRoLu7G42NjVCpVCZB0l07huZ0dHTYrEZjK+zpSdI7Iao17KyLXSq1tubiLtjrJYmJiYiKihr2fbS0Z+iIHB8weJCGy10vS2sWbIk3dw/yEJk0LsTDLa28aLVa0DSNuLg4hIeHu801xBJyuRwVFRWIjY1lhvcoisKOHTvw0UcfYdeuXcjJGdrkmscl8EHRXsyDoqv3DUkfhUy3RkdHMydduVzO7EC5asfQHLJ3KRaLIZFI3J6ZsodC5HI5+vv7QVEUwsPDERMTM2ypy5X09fWhsrISvr6+kEgkTq2XsPuTlvZCx4wZAx8fH+Y9s3WQxl2QNQvirGDt39wdgzz2rjq4GuJIHx4ejuDgYEZGUaVS2ezX6ChGoxE1NTVQKBTIyspiWhO1tbVYuHAhCgoK8OKLL7rte9jY2IhHH30UbW1tEAgEmD9/PhYvXozu7m7cf//9zOHw4MGDjITcNQ4fFO2FBEVrfUNn6OnpgVQqZU7KQ91Y2DuG5ATv4+NjUuqy58bkbjWa4SCl0uDgYERGRjL9LqVSCZqmTcqu9jquDwdxdOjp6RnW8NkZ2N5/CoWC6U96eXmhp6cHkZGRSElJ4bxk5+yahTVFHkvlSXOIgDfplXP52o1GI/O+sx3p2ej1eiZIkl1esVhsEigdPZCSYEyqMQKBAEajEXv27MG+ffuwfft23HDDDa54qVaRyWSQyWTIzs6GUqlETk4OPv30U+zbtw8hISFYsWIFNm7ciJ6eHmzatMmt13KFwAdFeyEnZVfuG6rVakilUhiNRquDLLZcl1qtNhl2Ya8YkLKr+bWyy4WxsbEu144cDltKpWS3iy2CLhaLB/W6HIGUid2hmzkcZP9MqVQiMDAQarWaybrYawbuypTduWZhPshDnFDI59Hf3x8ymQxKpdJhAW9nIMGYHZBshe0ZSrSG7Skpswd52MG4ubkZCxcuhEQiwSuvvML5vwkA3HvvvSgsLERhYSF++OEHREdHQyaT4eabb2aUu65x+KBoLx988AE+/fRT5ObmIj8/HxMmTHC4p0JRFFM2kkgkQ6riOII1Rw0STDw8PFBXVwc/Pz+kpKRw2r9ydqqU3JjMxcJtXcYnbhIeHh5IS0vjdCiApmm0tbWhrq5u0CANeyjEPFO2RwZtKPR6PaOPy6VmKBnkaW1tRWtrKzw8PAYJvLt7kIdkh8Tn0RWBx3w1xJINFVkNIRJxbBEAo9GI999/H9u3b8c///lP/P73vx8REe/6+nrceOONKC0tRXx8PORyOfP6goODmd9f4/BB0V70ej1+/vlnnDp1CsXFxSgtLYWfnx9yc3ORl5eH/Pz8Yf3r2DdFkqFwVTYi/oz19fVQq9WM+TE763J3D88dU6WWlvEtBRMAI9q7c2SQhsiFsdcnHOkps6dabR0iciUGgwFSqRRqtRqZmZnw8fFhBnnYDi/uGuRxJju0F/PKDSm9Go1GjB07Fmq1GklJSdBoNFi8eDGCgoKwdevWEevb9fX14aabbsLKlSsxa9YsBAUFmQTB4OBg9PT0jMi1cQwfFJ2Fpml0d3ejuLgYRUVFKC4uhkwmg0QiQV5eHvLy8nDdddcxQxW//vor1Go1/P39kZyczPl0IRniITdFAEw/iHyBAZjccF01XMD1VKl5pkwyyoCAAMTGxiI4OJjxnnQ3ZL2lvb0daWlpTt/8rOnWsrMu9meLqMIQOy+uHQxsFfA2N/9VKBSDPA3t7ZeTnrFCoRiRUq1KpUJ5eTkCAwMRGRkJpVKJf/3rX/jyyy/R3t6O6667DrNnz0ZeXt6wFnPuQK/XY8aMGbjjjjuwZMkSAEB6ejpfPh3qD/FB0T4oikJlZSWKiopQUlKCc+fOQa/Xw8PDAxRFYefOnZgwYQKnQwXsQZbhhnjYqxOWBLXtPb3TNI2mpiY0NzePyAK+TqdDVVUV9Ho9UlJSTHpC7AEl8vpcrRIil8tRWVmJ8PBwxlXB1ZhPhbKDCenvpaenD3JPcTfk356iKGRkZDhUpnbFIA8Ry+fyc8deccnMzGQGuHp6erB06VIYDAa89tpraG1txenTp3HmzBksW7YMGRkZnF7jY489hpCQEGzdupV5fNmyZQgNDWUGbbq7u7F582bOrmsE4YOiu9Hr9Xj11Vfx7rvv4r777oNQKERJSQlqa2sRExPDZJO5ubkIDAx0+ZeWbbrrTHZGbrjkF1uxhi0Ubs5ILuCzM2Nip2Vp0Mi8zOUq70kiAKBWq5GRkeESw2N7UCqV+PXXX+Hp6Qlvb2/09fWZTPISfVd3BAp2iyA5Odnl3n7WBnlIkPT390drayt6e3tNVh24Qq1Wo6ysDAEBAUhJSYGHhwdomsY333yDVatWYfny5XjwwQdH3AD4xIkTuOGGG0wO6Rs2bEBBQQHmzJmDhoYGJCQk4ODBg5y3GkYIPii6G4qisHfvXjz66KMmp2Syl1ZUVIRTp07h9OnTUKvVmDBhAhMos7KyHM5a2HtvlvwVnYWm6UH+jGzJOl9fXzQ2NjodjB2lt7cXlZWVCAwMRHJyst17c854T7IDwkj07siaBZEpY69ZsF0a2MNXrhx2YZdqiecfF5BBnra2NshkMmaQhz18xYUiT0tLCxobG5Gens6UyZVKJVauXImWlha8/fbbbtMGfeKJJ/D5558jIiKCEeQexTuHjsAHxSsJrVaL8+fP49SpUzh16hTKy8sRFBTEBMm8vDybNCiJTmpkZCSnu18GgwEKhQKNjY3o7u6GWCweNMTj7hukwWBAdXU1+vr6kJGR4TLdSlu9J/v7+1FRUQFvb+9BTvRc0NPTg8rKSkRHR9s8TGJp55W8NhIsbXkd7IBgScDb3VAUhZqaGmbNw9fXd8hBHku9V2fQaDQoLy9n3nuRSASapnHixAk8//zzWLhwIebNm+fW7+Px48fh7++PRx99lAmKzz///GjdOXQEPiheyZBpwVOnTjH9yc7OTqSlpTGTrpMmTWJ29Gpra9Hb2wsPDw/OdVKBwaVSoVDICIWTm5KrSpPmDLXm4A4seU+S6cLo6GhERUVZLSm7A1euWZi/tt7e3mEF3h0V8HYVRCZtuH1TdwzysHd92YbjarUaa9euRWlpKfbs2YOkpCSXvNbhqK+vx4wZM5igOIqHZhyBD4pXGwaDAWVlZUzZ9eLFixAKhfD29kZrayvefvttt9rcWMKeqVJrpUm2/qm9wVylUqGysnLEsjP2IE1oaCgz8cruc5HXxpZ2cwVcrVlYEngHgICAAKZsmZWVxXlZjoho9/X1ISsry6HDgD2C4ebodDqUl5dDJBKZlIrPnDmD5557Dg8//DAWLVrEaS/dPCiy1ytG2c6hI/BB8WqGpml88MEHWL9+PW6++WaEhITgzJkzuHTpEhITE5lsMjs72y0DFWSqtKmpyeogiy3o9XqTIR7zRXwiNGAOu3c2EtJ0JDvTaDRWB2mICwrJuFzpPUl6d+SGzPUov0KhQFlZGTw9PZnSsav7k0NBSsXuUCNia9daGuQZM2YMVCoVamtrIZFImJ69TqfDxo0bceLECbz99tvIzMx02TXZylBBERhVO4eOwJsMX+20t7fj+++/N7F1IiLDRUVF+Oyzz7B27Vro9XpMmjSJ6U2mp6c7dXpll0rz8/Od+llisRhhYWHMa2Av4hNTXHP9U7VajZqaGowdO5bzzNgeJw1LBrKkpCyXy3Hp0iW7vSfZKy4j0btjC3iPGzfOxLqI7arR0tJi0p90VV+ZnR26y/hYKBQy10wgB5yenh7U1NQw79vXX38NnU6HsWPHYsOGDZg5cya+//57zisW1oiMjIRMJmPKp1ybZF+L8JniNUB/fz/Onj3LKPFUVVUhIiKCkavLzc1FaGjosKdtUirVarVIT0/nbKqUTE0Srz+KouDv74/g4GAmM+FCqs0dgzT2eE8Se6Xh3CzcBRGxtlXA25H+5FCQ7JBtscQlnZ2dkEqlzJqJTqfD8ePHsXv3bpw/fx7e3t5IT09Hfn4+5syZw+nOIcE8UxzFO4eOwJdPRytkUpD0JouLi5mdLlJ2HT9+PFOSMxgMqK2tRVdXl1OlUkdhr5iQ7IitN6lQKEzEtMnupKuCBluRJj093e2lWkvek8TYlQwScVkudaUqjNFoHKTvyj4EjBkzZpATCkVRkEql6O/vZyTiuMRgMKCqqgo6nQ6ZmZnMAUwqlaKwsBDTpk3DmjVr4OnpiYaGBpw+fRrp6emYMGECp9f5wAMP4IcffkBnZyciIyOxdu1a/PGPfxytO4eOwAdFnsvo9XpcvHiRkasrLS2Fr68v4uPjcfbsWfzf//0fFixYMCKThZWVlQgLC0NiYuKQXn+WMi52j8sRybqenh4ms+ba3og8Pxnk8ff3ZwZ52D6G7tStJc/vDjcNgrkTCulPBgYGQigUoq2tDXFxcZw7mQCX3UTi4+OZqWaKorB7924cOHAAO3bswPXXX8/pNfG4DT4o8linra0NS5YsQWVlJX73u9+hvLx8SF1XV8Mu1TqqCEP6QCQrYQ+6kF/WSqC2DNK4E/bzW8qOhnLUcIX3pCUBby4hhtcqlQre3t7Q6/Umknzu3ntlZ6dZWVnMVHRjYyOeeeYZZGVlYdOmTZx+Lo4ePYrFixeDoijMmzcPK1as4Oy5Rwl8UOSxznvvvQcfHx/ce++9Jg7xRNe1uLgY586dAwBkZ2cz/UlnjXJJabehoQHJycmIiIhwadBl704qFAqmP8mWPmtvb0d9fT2SkpJsEkxwNW1tbQ6tWbjKe5L0zrjY+bQEyc7i4uIY1xnSnzTfe/X392eyZVfthpK9R3bv0mg04j//+Q927tyJLVu2YPr06ZwrFaWlpeHrr79GbGws8vLy8N577yErK4uzaxgF8EGRxzmI3NuZM2cYgQEyFZqfn2+3rqtSqURFRQXGjBmDlJQUTkq1pMelUCjQ1dWFrq4ueHh4ICIiAsHBwTZbM7kCjUZj4vPoir6hPd6TrhDwdgaSnZLseLgAzn7vyOrExXpaKgAAGUpJREFUcP3J4X5eTU0NFAqFiWZqa2srFi1ahIiICLz22msmU6lcUVRUhDVr1uCrr74CALz88ssAgL/97W+cX8s1DB8UXcGaNWvw9ttvM7tKGzZswN133w1g4IO7Z88eeHh4YPv27bjjjjtG8lI5YShdV5JNmuu66nQ61NXVobe3FxkZGSZ6nVxdc319PTo6OpipWnPpM1K6IzdbVwZsLtcsrHlPikQiqFQqxMfHIz4+nvPeqaXenSOQISV2yVwsFg/anzSHTNay/RZpmsbHH3+MzZs3Y8OGDZgxY8aIiXgfOnQIR48exe7duwEABw4cQHFxMXbs2DEi13ONwu8puornnnsOS5cuNXmsrKwM77//Pn799Ve0tLTgtttuQ1VVFedj9FwjFAqRmJiIxMREPPDAAwAGdF0vXLiAoqIibN26lVkryM3NBTDQKzl48CByc3M5v+mQQZKoqCjk5eUxwSA0NJQJTmw3DeINSFRPSKB0VCChr6+PyY7z8vLc/vkQCATw8/ODn58foqOjodVqUV5eDqPRiLi4OCiVSpSUlJgs4pOVF3e8N+zscPLkyU7LE3p4eCA4ONhEXYdtF9bU1AStVssccgICAtDT0wO5XI7x48czk7VdXV3461//Cg8Pj0G7wDyjGz4oOsjhw4cxd+5ceHl5ISkpCRKJBCUlJZg6depIXxrneHl5oaCgAAUFBQAGgkxxcTEWLlwImqYRHR2NBx54wKquqztgy9MNtwQuEAjg6+sLX19fREdHAzA1Ma6vr2ccJ9hDPEOVH41GI+rq6tDV1YWMjAyTJXguYPduU1NTB9302UpDZBHf1d6TXV1dkEqliI+PR0ZGhtsORJ6enggPD2eqOeSQ097ejl9//RVCoRBisRj/+te/YDAY4O/vj3feeQerVq3C/fffP+IWTwAQExODxsZG5vdNTU1uc9vgGRo+KNrAjh078M477yA3NxdbtmxBcHAwmpubMWXKFObPxMbGorm5eQSv8sqBpmls2bIF27dvx+9+9zsAprqu77zzDi5evAixWIycnBym7OqKlQiapiGTyXDp0iWnBmksqZ6w+3ckIyH9O/baBFkziYyM5FyRBxgQqy4vL4evry/y8vIsBjdLSkMkW+7s7ERtba3DAu9k70+r1bokO3SE9vZ2tLe3Izs7GwEBATAajTAYDNi+fTuqq6vh5eWFbdu24cSJE1iwYAHnO4fm5OXlQSqVoq6uDjExMXj//ffx7rvvjug1jVb4niKA2267Da2trYMeX79+PaZMmYKwsDAIBAKsWrUKMpkMe/fuRWFhIaZMmYKHH34YAPDkk0/irrvuwn333cf15V+V0DQNhUKB06dPM9Ouzuq6qlQqVFRUwNfXFxKJxO1SXGzXeHZ/UigUIj4+HuHh4Q7tTjpzPcQNnu335yiOeE+S7DAhIYFzr0lg4DNQVlaG4OBgJCcnQygUgqZp/PTTT1i+fDkWL16Mxx9/HEKhEH19fTh37hzi4+ORmJjI6XVa4siRI3j22WdBURSeeOIJrFy5cqQv6VqDH7RxNWyJJfPpsDvuuANr1qxxuHy6ZcsWLF26FB0dHQgLCwNN01i8eDGOHDkCX19f7Nu3D9nZ2S57LVciZDqQeE6ePXsWOp0OkydPHlLXlZQqOzs7kZGRMSLTg+3t7aipqUFcXBz8/PxM1Gq8vLxMAok71Gq4koiz5j1JhAcAYNy4cZxPtrIPBJmZmcxnoL+/H6tXr0ZlZSX27NmDhIQETq+L54qCD4qugIjtAsBrr72G4uJiZsDmwQcfRElJCVpaWjB9+nRIpVKHbkaNjY2YN28eKioqcPbsWYSFheHIkSN4/fXXceTIERQXF2Px4sUoLi529cu74rGm65qTk4P8/Hz09fXh8OHDeOmll0ZkqtKWNQutVmuSTer1ehODZmf279gC3iPVu5TJZKipqYGfnx8oihq0G+pu70m1Wo2ysjIEBAQgJSWF+Q6WlJRgyZIl+POf/4yFCxdy9tn48MMPsWbNGpSXl6OkpIQZOANG58T6FQQ/feoKnn/+eVy4cAECgQCJiYnYtWsXgIHT8Jw5c5j1gzfeeMPh0/lzzz2HzZs3495772UeO3z4MB599FEIBAJMmTIFcrncJECPFnx9fXHDDTfghhtuAHB5gOTYsWNYvXo1ent7ER0djXXr1lnUdXUXNE2jubkZTU1NkEgkQ04venl5ISIignEwIPufpDfZ19dnIllnqzdjb28vKioqEB4ePiK9S71ej6qqKhgMBuTn5zPZIXu/sKmpyW3ek+Sz0NjYaFIu1mq12LBhA4qLi/Hee+8hPT3d6ddqD+PHj8fHH3+MBQsWmDw+WifWrzb4oDgMBw4csPr/rVy50um6/+HDhxETE4NJkyaZPN7c3Iy4uDjm92SQZ7QFRXMEAgH8/f2xfft2rFmzBrNnz4bBYMDFixdx6tQp7Nq1i9F1JQM8+fn5jHKKK1CpVCgvL3d4zYIsoAcEBCA2NhaAqTdjW1sbU5a0ZMtEBLzlcjnGjRvHmZsJG7K6YmmYSSgUDvn62tvbB0nyjRkzxq6DjEajQXl5OXx8fEzeg4sXL+Ivf/kLZs2ahe+++45zLV8AVn0W+Yn1qwM+KHLAUIM8GzZswLFjx0bgqq5eAgMDUVxczNxEyRRrTk4OswbS3d2N4uJiZtq1paUFEokEubm5yMvLQ3Z2tt3ZijvXLMy9GdmyZ+xpUE9PT/T19SEqKgrZ2dmcZxns7DA7O9vm3uFQ3pM9PT2or6+3yXaK7XeZlpbG7Jrq9Xps3boVX375Jfbs2YOJEye67kW7CH5i/eqAD4oc8M0331h8/JdffkFdXR2TJTY1NSE7OxslJSX83tIwDJVVCAQChIaG4u6772bUh9i6rh999BGT4RNd17y8PEgkEqslSK7XLAQCAXx8fODj44OoqChmCV6pVCI2Nhb9/f04ffo0PDw8THYn3bn+MFR26Aje3t7w9vZGZGQkAFMnlJaWFhMnFOKCUl9fD7FYjNzcXCZzrqysRGFhIW655RYcP36cE9utoQ667DYIz9UHHxRHkAkTJqC9vZ35fWJiIs6cOYOwsDDMnDkTO3bswNy5c1FcXIzAwEC7S6erVq3C4cOHIRQKERERgX379mHs2LGjcrJVKBQiMzMTmZmZeOKJJ5gbMNF1XbNmDaqrqxETE2Oi60pRFN566y3ccccdmDBhAuduGoCpgLf5Ejx7Cb+5uZnZnWTLnjlbQtTr9aisrARFUcjJyXFb0CGlcX9/f+YASGTdmpubUVlZCZFIBF9fX6xZswbjxo1DU1MTPvvsM+zcuRP5+fluuS5LWDvoDoWzB90XXngBISEhePbZZwEMtG8iIiKwePFiu6+Fxzr89OkVBDso0jSNwsJCHD16FL6+vvj3v/9tMsVmC729vUyJb/v27SgrK8Obb77JT7ZawVzX9dixY+js7MTUqVMxffp0FBQUDNJ1dSfsYGSrgDd7CZ+tfeqo5RTJDokbPdeQfwOj0YiMjAx4enpCo9Fg7969+OSTT9Dc3IwxY8YgMzMT+fn5WLBgAfz9/Tm/TkvcfPPN+Oc//8l8b52dWK+vr8esWbNw7tw5GI1GpKamoqSkxK1autcY/PTp1UZ9fT3z3wKBAG+88YZTP4/d81KpVMyNkJ9stQzRdY2JicHHH3+MnJwcbNy4Ec3NzYN0XcneZF5enlvsp4i9lL3ByJJkHdtyqra2FiqVCp6eniaSbuYBlx2M3JkdDgXJkNn/BkajEe+99x4OHDiArVu34pZbbmFK4yUlJW4XbLCFTz75BH/5y1/Q0dGBe+65B5MnT8ZXX33l9MR6YmIiQkNDcf78ebS1teG6667jA6Ib4DPFa5yVK1finXfeQWBgIL7//nuEh4djxowZWLFiBaZNmwYAmD59OjZt2mR3Jnotc+HCBUyePHnQ4zRNo729nRniKS4uZtw3XKHrqtVqUVFRAQ8PD6Snp7vtJs+WrFMoFNDpdMyQi9FohEwmQ0pKyohkh0QmTqfTITMzkwnYMpkMhYWFiIuLw5YtWzh3W7kS+OCDD3Dy5Em0trbiscceY3rmPDbBL++PBmxt+L/88svQaDRYu3YtHxRdDFvXtbi4GD///DPEYjGys7OZQDmcrutwAt7uhqZpyOVyVFVVQa/XQyQSmewWkkEXd8u2WbKYomkaH374IV599VVs3LgRd911F2fyccuWLcNnn30GT09PpKSk4N///jeCgoIAjMwivk6nw4QJE6DX6x0WCxnF8EGR5zINDQ24++67UVpaigULFuDmm29mrJ/S09Pxww8/2FU+vdJuFlcSw+m65uXlIScnh9F1ra6uRk9PDwICAiCRSEZkt47I1KWkpDAiA2TIhWST5ruFgYGBLstkKYqCVCpFf38/srKymEy7o6MDS5YsgY+PD7Zv326y0sEFx44dw6233gqRSITly5cDADZt2oSysjI88MADTH+Qy0X8p556CkFBQdi4caPbn+sagw+Kox2pVIrU1FQAwOuvv44ff/wRhw4dwhdffIEdO3YwgzaLFi1CSUmJXT/7SrxZXMkYjUbU1tYyQzxnz56FVqtFYGAgGhoasGXLFtx2222c/zvpdDpUVFRAIBAgPT192N4h2S0kQzzEiskZSTe5XI6KigrExsYiJiaGyQ6/+OILrFu3DqtXr8bs2bNH3OLpk08+waFDh/Df//7X5drHtmI0GpGdnY0PP/yQ+W7z2Aw/aDPaWbFiBSorKyEUCpGQkIA333wTAHD33XfjyJEjkEgkzGSrvdx+++3Mf0+ZMgWHDh0CwKt2WEMoFEIikUAikeCRRx5BeXk55s+fj9jYWNx555149913sXr1akayjayFhIaGui0YkGEednY4HOa7hZYk3YZz0iAQAXiFQoGJEycy6y5yuRzLly+HUqnEN998MyJ9TUvs3bsX999/P4CRWcQvKyvDjBkz8Kc//YkPiG6ED4rXMB999JHFx10x2cpmpG8WVyP/+9//sG3bNpP9UNJXLCoqQlFREbZt2waFQoGsrCym7DphwgSnJ0FJdigUCk2W4B3BkqQb20lDJpOZGBiTjLK/vx/l5eWIiopCTk4Okx3+8MMP+Nvf/oa//vWveOSRRzjRc7WlL79+/XqIRCI89NBDbr8ea2RlZaG2tnbEnn+0wAdFHqtcLTeLq5F58+YNekwgECAmJgb33Xcf48up1+vxyy+/oKioCG+99RZ++eUX+Pn5OazrSrJDiUTCONW7GrFYjNDQUGZdgEjWyeVytLe3o7S0FAaDAaGhoWhsbERbWxvS09Oxdu1a1NXV4YsvvjDR/XU3wy3i79u3D59//jm+/fZb5t+ZV5y6duGDIo9VRvJmwdvvDECmWLOzs010XUtKSkx0XVNSUphs0pKuq06nQ3l5OTw8PJzODu2FSNZRFIXGxkbExcUhLi4OKpUKX3/9Nfbu3cvsf/7xj3/EmTNn4OnpeUWUTY8ePYrNmzfjxx9/NFEzmjlzJh588EEsWbIELS0tkEqlnCrq8LgPftCGxyGOHj2KJUuW4McffzTJOFzlM1leXg6hUIgFCxaYqILwgzyDYeu6FhcX4/z586BpGtnZ2cjJyUFraytOnz6NHTt2uC07HAqapnHp0iW0t7cjMzOT2S/UaDRYv349zpw5g927dyM4OBglJSU4deoUkpKS8Oc//5nzazVHIpFAq9UyWe+UKVOY3vz69euxd+9eiEQibN26FXfddddIXirP8PDTpzzug6ubhblU1khN/V1NEF3XY8eO4YUXXgBFUfD29kZERISJrmtgYKDbJzpVKhXKysoQHByM5ORkpkd4/vx5LFq0CHPnzsWSJUtG9aGGhzP46VMe91FdXW31/3OFz6Q1+EGe4REIBBAKhdi4cSM2b96Mu+++m9F1PXXqFL777jts2rQJ/f39GD9+PCMw4EpdV5qm0djYCJlMhoyMDAQGBgIY6JG+8sor+Pbbb7F//36MHz/eJc9nC7xAPo8t8EGRZ8Tg7Xfch6+vL06ePMkEOaLrmpiYiLlz5wIYkJS7cOECo+taUVGBwMBAJpPMz893SNdVrVajrKwMAQEByM3NZbLAsrIyFBYW4s4778Tx48c51yldtmwZ1q1bB2BAIP/FF1/Em2++iS+//BJSqRRSqRTFxcV4+umneYH8UQwfFHlGjJGw32Fz9OhRLF68GBRFYd68eVixYoVDP+dKZbisz8vLCwUFBSgoKAAwWNf1rbfeskvXlayUNDY2Ij09HcHBwQAG1GreeOMNHDp0CG+++eaIyQnyAvk8tsAHRZ6rCldN/VEUhYULF+Lrr79GbGws8vLyMHPmTGRlZbnhqq8OBAIBIiMjMXPmTMycORPAZV3XU6dO4cCBA1i6dKlFXddLly7hhx9+wNSpU5GXl8dkh3V1dXjmmWeQn5+PEydOuNUE2RbMBfKBgZI8ewWElOT5oDg64YMizxWJu+x3CCUlJZBIJEhOTgYAzJ07F4cPHx7VQdESIpEIEydOxMSJEzF//nzQNI3e3l5mJeTQoUO4ePEijEYjZs6ciaSkJMTExMDPzw979+7F3r17sX37dtx4442cXO9wJfn169dj/fr1ePnll7Fjxw6sXbuWk+viuXrgp095RiWHDh3C0aNHsXv3bgDAgQMHUFxcjB07dozwlV09dHd3Y/78+QgICMCiRYtQWlrK6LrW1NTgD3/4A7Zv337FmP6ycbVAPs9VAT99ysPD4z68vLwwb9483HnnnQCA6667Do888giAAYPgkJAQTmTabIUtkH/48GFkZGQAGCjJ79ixA3PnzkVxcTECAwP5gDiK4YMiz6iEl+lyHj8/PyYgmsO1H6QtuFMgn+fagS+f8riV06dP48knn0RJSQkoikJ+fj4++OADTvfTLGEwGJCWloZvv/0WMTExyMvLw7vvvotx48Y59XOfeOIJfP7554iIiEBpaSmAgTLj/fffj/r6eiQmJuLgwYPMZCYPDw9n2FQ+vXJqGzzXJGSq8x//+Aeef/55PPzwwyMeEIGBAZIdO3bgjjvuQGZmJubMmeN0QASAxx9/HEePHjV5bOPGjYzc3fTp03lzWBeyZcsWCAQCdHZ2AhhYC1m0aBEkEgkmTpyIc+fOjfAV8lxt8Jkij9vR6XTIy8uDt7c3Tp48ec1LetXX12PGjBlMpsge3JDJZLj55ptRWVk5wld59dPY2Ih58+ahoqICZ8+eRVhYGI4cOYLXX3+dMdBevHgxv4jPQ+AzRZ4rg66uLvT19UGpVEKj0Yz05XBOW1sbM7gRFRWFtra2Eb6ia4PnnnsOmzdvNlHcsbaIz8NjK3xQ5HE7CxYswLp16/DQQw9h+fLlI305I4pAIHC7CPdo4PDhw4iJicGkSZNMHre2iM/DYyv89CmPW3nnnXcgFovx4IMPgqIoXH/99fjuu+9w6623jvSlcUZkZCQjGyaTyRAREeHwz2psbMSjjz6KtrY2CAQCzJ8/H4sXL74mh3mGWsTfsGEDjh07NgJXxXOtw/cUeXhcjHlPcdmyZQgNDcWKFSuwceNGdHd3Y/PmzQ79bJlMBplMhuzsbCiVSuTk5ODTTz/Fvn37EBISwjxHT08PNm3a5MqXdcXwyy+/YPr06Yzpb1NTE8aOHYuSkhKsXr2aX8TnsYZtJRqapu35xcPDMwRz586lo6KiaJFIRMfExNC7d++mOzs76VtvvZWWSCT09OnT6a6uLpc938yZM+ljx47RaWlpdEtLC03TNN3S0kKnpaW57DmudBISEuiOjg6apmn6888/p++8807aaDTSRUVFdF5e3ghfHc8VhE1xjs8UeXiuUurr63HjjTeitLQU8fHxkMvlAAYOusHBwczvr3USExNx5swZhIWFgaZpFBYW4ujRo8wi/ki5cvBccdiUKfJBkYfnKqSvrw833XQTVq5ciVmzZiEoKMgkCAYHB6Onp2cEr5CH54qDX8ng4bkW0ev1mD17Nh566CHMmjULwOVhHgBOD/Pw8Ixm+KDIw3MVQdM0nnzySWRmZmLJkiXM4zNnzsT+/fsBAPv378e9997r8HNoNBrGUHjcuHFYvXo1gAFvxIKCAkgkEtx///3Q6XTOvRgenisQvnzKw3MVceLECdxwww2YMGEC40CxYcMGFBQUYM6cOWhoaEBCQgIOHjyIkJAQh56DpmmoVCr4+/tDr9dj2rRp2LZtG1599VXMmjULc+fOxVNPPYVJkybh6aefduXL4+FxJ3xPkYeHxzn6+/sxbdo07Ny5E/fccw9aW1shEolQVFSENWvW4KuvvhrpS+ThsRW+p8jDw+MYFEVh8uTJiIiIwO9//3ukpKQgKCgIItGA3gevFMNzrWJvpsjDwzOKEAgEQQA+AbAKwD6apiW/PR4H4Euapkfe8oSHx4XwmSIPD49VaJqWA/gewFQAQQKBgEhDxgLgU0Weaw4+KPLw8JggEAjCf8sQIRAIfAD8HkA5BoLjfb/9sccAHB6ZK+ThcR98+ZSHh8cEgUAwEcB+AB4YODgfpGn6RYFAkAzgfQAhAM4DeJimae3IXSkPj+vhgyIPDw8PD89v8OVTHh4eHh6e3+CDIg8PDw8Pz2/wQZGHh4eHh+c3+KDIw8PDw8PzG3xQ5OHh4eHh+Q0+KPLw8PDw8PzG/we8sLgRtCOb3wAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "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" } ], "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" } ], "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" } ], "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" } ], "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" } ], "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" } ], "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" } ], "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" } ], "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" } ], "source": [ "prody.showProtein(prot,prot,flucts=[gnm[1].getArray(),prody.calcSqFlucts(anm[0])],vecs=[anm[0].getArray(),[]])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.12" } }, "nbformat": 4, "nbformat_minor": 1 }