acentral reconstruction

acentral reconstruction.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ancestral Reconstruction\n",
    "\n",
    "\n",
    "## Methods\n",
    "\n",
    "1. **Uniprot**에서 목적 단백질을 검색\n",
    "2. Sequnce 정보를 얻기 위해 **BLAST!**\n",
    "3. Clustral Omega 에서 *MSA* 실행\n",
    "4. http://fastml.tau.ac.il 에서 업로드해서 분석 시작\n",
    "5. 결과 파일로 Ancestral reconstruction sequence와 tree를 다운로드\n",
    "\n",
    "## Results\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "from Bio import Phylo\n",
    "\n",
    "def plot_tree(treedata, x, y):\n",
    "    #handle = StringIO(treedata)  # parse the newick string\n",
    "    tree = Phylo.read(treedata, \"newick\")\n",
    "    tree.ladderize()\n",
    "    matplotlib.rc('font', size=7)\n",
    "    # set the size of the figure\n",
    "    fig = plt.figure(figsize=(x,y), dpi=100)\n",
    "    # alternatively\n",
    "    # fig.set_size_inches(10, 20)\n",
    "    axes = fig.add_subplot(1, 1, 1)\n",
    "    axes.axis('off')\n",
    "    Phylo.draw(tree, axes=axes)\n",
    "        #plt.savefig(output_file)\n",
    "    return"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyIAAAEdCAYAAADw/tMPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xm8XdP9//HXWwxFoqaqr68hpqoICaIxpD+poKVKaYih\niJbqt1WdtLRUg68aSluK9osS1JyadTCmxiKjDBKiIjWVIJUIIfL+/bHWSXZOzj733Js7JDef5+Nx\nH3eftdde+7PXPjfZa6+19pZtQgghhBBCCKE9LdfRAYQQQgghhBCWPdEQCSGEEEIIIbS7aIiEEEII\nIYQQ2l00REIIIYQQQgjtLhoiIYQQQgghhHYXDZEQQgghhBBCu4uGSAghhBBCCKHdRUMkhBBCCCGE\n0O6iIRJCCCGEEEJod9EQCSGEEEIIIbS7aIiEEEIIIYQQ2l00REIIIYQQQgjtLhoiIYQQQgghhHYX\nDZEQQgghhBBCu4uGSAghhBBCCKHdRUMkhBBCCCGE0O6iIRJCCCGEEEJod9EQCSGEEEIIIbS7aIiE\nEEIIIYQQ2l00REIIIYQQQgjtLhoiIYQQQgghhHYXDZEQQlhKSeovyZK2qkobViPvcEmTJI2RNEHS\noTm9u6QRVXmHStqnzn7Xl3SLpOcljZB0s6RP1tp3Xt89L0+V1LU5ZUi6XNJVklQSy6n5eMblcjZu\npO7am6TBks7Ly3XrdzH2cbqkz7ZgO0k6TdIUSc9JulvShoX1DddxPq9/kvRPSSMl3SfpM03sf359\nFL8jIYTOb/mODiCEEEKLDQIeAw4Cft5A/oG2x0taFxgDXNfcHeYGwW3A72wfkNM+C3yitcuQ9Ctg\nDeAg265Rzs7A54Detj+UtD7wbjPiWM72vEbzL+lsn9rCTY8HtgN62n5f0mHArZL6ADvRYB0Xzusl\ntr+S03oDnwaebGFsIYROLHpEQghhKSSpC7AX8DVSQ6Q5ugIzW7jrAcC7tv9QSbD9sO3xrVmGpNOA\nHsAhtj8qKWddYLrtD3MZL9l+O2//eUmPSxot6Y+SVszpb0q6SNI44FOSphf2eZykIXl5uKTz8139\npyVtl3sKpkg6ruzAJG0q6WFJoyQ9mS/Ea9k7l/2MpF3ztvN7TfLnEbnHqruksZKuzT0Wv5P0ZUlP\nSBovafOcv7pnYYhSD9hTkv6rLGbgR8Dxtt/P9XgtMIt0nkrruIbdgZm2r6gk2B5j+4Yc03455tG5\nLlevU48L9dRJOk/S4Lz8S0mTc52cUqeM/SQ9m4//D4UeqW/mtLGSrpO0Qp26CSG0oegRqUPSqh0d\nQwghANiuvgu9GzDC9mRJ0yX1sj22iWKGSfoQ2BwY3MJQegCj6qzfXdKYwudPtaCM3YAtgM/Y/qBO\nvnuB0yRNzMvX2B4haW3SxfVutt+TdDpwDHAxsCbwF9vHAaj2iK+KWba3l3QycCOwQ05/BrioZJtX\ngd1tz5G0DXA+sEeNfP8N9AG2BG6TtEW9QHK+g4ApwPgcW19JxwLHAd+tsc1Ltnvn4z8aOKM6g6TV\ngFVsv1C1ahTpPF1JjTquE+OYknUADwF32Lak44FvA2fWyb8ISWuRegK7254n6eMl+VYGLgR2AV4D\n7gMqcd9k+/c5369I9Xptc+IIIbSOaIjUN6ujAwghhKz6inkQcHNevil/bqohUhmatTEwXNLdwCJD\nnrKydNVZB3Cf7YHzM1fNP2mwjMmkO/GfBe4py2R7pqRtSUOHBgD3SjoIWAnYBng8NzRWAu7Om71n\n++5a5dVwR/49jtTom5GPaaakNUp6BlYCLsqNkI8oH7J2Yx5uNlHSLFLDpJ7Jtifn/T9DurCuxLZ3\nyTa35t8jgX1L8pS1xAS4rI5t31u2zfwPaa5PD+BR28cAGwI3S/oksDLwRMm+6/lP/rlC0m3AXSX5\ntgAm2X4px/InYKO8rpekM4CP55/3WhBHCKEVxNCsGiStKqnef5IhhNBh8lCS/YBfSpoK/JhmDM/K\nd79fJd3BfpM0D6NoTWB69XbZRNJ8gsXRVBkvk47vinwRXMr2XNv32j6JdHd9P9L/bXfb7p1/trR9\nQt5kdnURheWVqtbNyb/nFZYrn7uUhPQ9YCqpIdSvRpm19lv5PJeF/18ublu9/2JsZbFU8nxUlsf2\nf4DZyg8UKNiW1PNTVse1TATmD0XLDdJvseD7dSFwru2tSfVUVjdQUhe255J6kv4EHAjcULJ9dQOr\n+PkPwDE5jnObiCOE0IaiR6Rp67Dof1whhNCR9gAesb1/JUHSQ3lycZPy0KVNgGm2Z0maIWln24/l\nychbAxNKNr8POEfSYNtDc3n9gBnNiL/JMmw/LekI0rClXW1PrXEcWwBzbT+v1PXRk3Qx/DhwoaSN\nbL+Yhx+tVWP4EcB/JG0EvALsA/y9GcdRy2rAlDz8aHCdfAdJuoE0kbtr3v+LwFH52HqQ7uq3h/OB\nC3JPxxxJh5CO4/46dVzLfcDZko60fVVOW7mwfjXg5VzOEU3E9DqwnqRupEbaHsBIpSdqrWL7TklP\nAo+WbD8J+LSk/yYNzdqf1DMEsCrwb6V5Q4eQvi8hhA4QDZGmza4xNjuEEDrSINLTiYpuBQ4mDVXZ\nW9JLhXX98u9hkt4HVgROtv1aTj8SuCRfsM8FjrVdc2hqvsD+MulC/2fA+6Q5C8cDazcSfKNl2H5A\n0k+AP0vqZ/utqqK6koZBrZY/jwR+m+eFHAP8KV9sziPdga/VEDkFeIB0sTqpkfibcEne7zEsGBpV\ny6ukOQurkO7OW9IjwPQ89GokuUeiHVxA6gWbkEcDTAH2zzHVrONahVSd1yGkOp0B/G/OchpwJ/AW\nqcG3Ua1yclkfSDoXGA1MIw1BA+gG3C6p0otxYsn270n6HvAgaSjXJOCdvHoIqe5fz+UD6fHHpCF4\nd5Qtl8UbQmgZ1Xgi4jIvT1Kv/CfcNRoiIYQQwtJF0qq231V6wtwtwGW2y+aUhBA6QMwRCSGEEEJn\n9K38BLfxpF6VRh9SEEJoJ9EjUkP0iIQQAkh6gkUn8g6w/eayHEuOZ2vgmqrkF4rzdpY0+THEB1Yl\nX2D7yhaUtcScD0kXkx7TW3Si7b+1dywhhOaJhkgN0RAJIYQQQgihbcXQrBBCCCGEEEK7i4ZICCGE\nEEIIod1FQySEEEIIIYTQ7qIhEkIIIYQQQmh30RAJIYQQQgghtLtoiIQQQghhiSSpvyRL2qoqbViN\nvMMlTZI0RtIESYfm9O6SRlTlHSppnzr7nZvLGSPpKUm9C+t2lTQy72uCpIFV5b4tacVC2l8kTS3E\nYkk/Lqzvm9MG589/zW+or6y/T9IXS+LsIukSSeMkjZf0mKSuZcfVkSQNkXRcXh4uqWcb7ONySZu2\nYLuVJF0s6XlJUyRdL2mNvK7hOs7nd3b+3oyV9JCkDVvhuBb5DhfW/VnSyou7j44SDZEQQgghLKkG\nAY8BBzWYf6Dt3sAA4FeLsd8Ztnvnss4GTgWQtB7p/TFH2f40sBdwlqRtC9u+BOyZ868JfKKq7OeA\nfQufDwLGFT4fD5wjaWVJBwLv2y57GeMgoBuwje2ewGDgw0YOUFKnuwa0fbTt51uw6dnA8sDmtjcj\nvQTziryuuXU8MX93egG3A99rQTwNs7237ffach9tafmODmApsIqkjo4hhBBC6NSq39klqQvpQv/z\nwG3Az5tRXFdgZiuFthowIy9/C7jS9tMAtqdJ+gXwQ+CrOc9NpMbFXcD+pNiPLpT3DjBd0ibAC0B/\n4IHKStvPSroVOBPYh1QHZdYFXnN+KZztZysrJB1OatSsCDxg+/uSugN3AhOA3pL2B66x3Sdvcx4w\n3vbQ3ItzbY5hVi7rl0B34Ie2b60VkKQdgV+TXng5EzjC9os1sn5d0m7APOCrtidIGgJMt31RLmu6\n7bUl9QdOAd4DegKXAW8B3yA1Cr5oe7qk4cBxtsdLmg4MJX1/Xgf2rfVeuPzuuCOA7rbn5eSzgCmS\nNqtXxw3oRv7uSFobuBLYKMc+2PZUSUNJ9bsDsGaO5bvAtsAttn+Sy1pR0k3A1sCTwNdtz83nqSew\nNqnhMwb4DPA0cLBtSzoN2Bv4GHCf7e+XBZzr7VpSY/pZ4LxcH+sCh9t+vBnH36RoiDTt9Y4OIIQQ\nQlgGVN/12w0YYXuypOmSetke20QZwyR9CGxOunPdUqtLGgOsAqwF7JzTewBXVeUdBXyn8Hk0sK+k\nlYCBwPdZuCECcDOpsfJ3YCQwt2r9GcC/gP9r4g7/MOARSZ8D7gOG2p4kaUtgP2CnfLF6dR7eNQHY\nEjjM9tO5YVLPFNu9JF0G/AbYndQQuQmo2RABJgL9bH8kaV9SA+KYGvmWz2V/HriY1CCrp3eOfTap\nAfcL29tJOgs4nNT4KVoL+KvtEyRdDRxA6s2qtikwzfb8hqvteZLGks53zTquE2eP/N1ZnfSd3j6n\nDwEetv0lSYOAC1nQM9bNdl9Jh5EaitsDrwKTJJ2f8/QEvmZ7hKTrSA3foVX73hI4BHgGeBDoBzwM\nXGD750p31odJ2sX2oyXxrwXcZvu7ku4BjgP+H7AHcBLpe9VqOl23XAghhBA6hUGkC3ZIF76DGthm\noO2tSRdk50jqBrgkb1k6LBia9SngWNKFMqQLy+rtaqXdAxxKugP9Uo3y7wC+RGqM3FRj/Q7A+8BW\nNdYtOAB7GrAFcDqwKvCE0nyaAcCOwIh8UbwjsFne7NlKj04D7si/xwGP2J5jezKwXp1t1gBukTSe\nNOSpR0m+6/Mx/A3YooGhYo/bfiP3arwE/KUQW/ca+WfZvi8vjyzJA7XP3/z0OnVcpjI0qzvwW1Id\nQGoU/DEv30Tqtago1vNztl+0/QEwBdggr5tiuzJP5IZcXrXJtifm3pvRLDjmAZKeBMYCu1B+TiDV\n24N5eTypN20e5fW8WKJHpGnrkFrfIYQQQmgHklYg3XndU9I5QBdgDvDTRra3/YKkV0kNkomki+Oi\nNYHpDYZzF3B1Xp4AbMeCC0dIQ2ieqdrmZtJwq9NK4psp6XXScJkTgC9U1klaHrgA+CLwK0lfsP3X\nsuDy/IA7gDvyHe+9gA+Ay2wvtP/cA1K8ppnLwjelV6oqfk7+Pa+wDIv2XhWdDtxt+9I8IX1oWeg1\n0urFU9x/MZ55pO9HtWL+j0ryQLrY30hSV9uzYP78mW3I57WkjieUlFd0F/C1knXF4y+r5+KxVddX\nrfpb5JglfYzUm9XH9qt5+F31eS4ro5F6XizREGna7FpjCkMIIYTQZvYg3YHfv5KQn0DUp5GN83j8\nTUhDbmZJmiFpZ9uPSVqfNM6+kQtJSMOy/pmXfwc8LmlYnoewPvAT4ODiBrbHSDqdfNe/xJnAhnkI\nUzH9O8CjtkdLOh64XtL9theZIC1pO+AV26/lxtungYfysd0s6SLbb0pah9oXka8D6xV6jvYg9R4s\njtWAl/Py4Dr5BpHqcndgUh4O9SLwuXxse5Dm+rQp2+9K+iNwrqTj8t3/E4GnbU+pU8eNKH53HiH1\nkp1HGrL3ZDND3VzSdrZHkXrS7m1wu4+Rzu2bkj4OfJnUMFkiREMkhBBCCEuaQaRJ3kW3ki747wL2\nllQc8lQZpjJM0vukCdon234tpx8JXCJpNdJd92Mrd79LVOaIKOf/BoDtl5Ues3tVnuT8EfCzfHG4\nENu/ASh71GseZlP9WOFPksbk75DzjJP0IGny8nk1ilkHuDxfIHcB/goMyxf1ZwL357v7c0iNgoVu\nrNr+QNK5pGE801j46V0tdS6pfs7I8ZSZl+vYLJjofwswWNI44H7gzVaIpxEnkuaYPJcbhU+xoCej\nZh3XKatH4bszkwXzg4YAQyUdQZ6s3swYxwMnStomx3ddIxvZniHpKlLP4CvAP5q530XkuT99bJ9a\nttxwWfkhAKEg/+NS+Qeqa/SIhBBCCCGE0LpisnoIIYQQQgih3cXQrBBCCCEskyQ9waITdwfYbq8h\nQQ2RtBZpqFLRHNt9OyIegPzY3XOqkh+1/e2OiKcRki4mPTWq6MT85K7mlLPEnY/mUnpXzcZVyYfb\nbo3heY3HEUOzFhVDs0IIIYQQQmhbMTQrhBBCCCGE0O6iIRJCCCGEEEJod9EQCSGEEEIIIbS7aIiE\nEEIIIYQQ2l00RJpB0tT89s3K5+MkDcnL+0gaL2mepJ4dFmQIIYQQQghLgWiINN8ukqofdwYwGRgI\nPNTO8YQQlmCS+kuypK2q0hZ5M6+k4ZImSRojaYKkQ3N6d0nVb2AeKmmfOvudWyjnTkmrV60fIum1\n/NblYtpLebtnJP2yKraeku7L66dJej0vj6kuP2/TRdIlksblGzWPlb1luqPlYz8uLw9vixtKki6X\ntGkLtltJ0sWSnpc0RdL1ktbI65pVx5K2lvSQpMmSnpXUKo9azd/pzzSQr03qNoSwdIqGSPNdAPy4\nOtH2c7YndUA8IYQl2yDgMeCgBvMPtN0bGAD8ajH2O8N2b9tbATOA6gvOgcALwK5V6Wfn/W8N7C5p\nm+JK27vn9acCV+d99LY9o0YMg4BuwDa2ewKDgQ8bCb7YQOosbB9t+/kWbHo26b1fm9veDBgPXJHX\nNVzHklYBbgeG2N4C6AscKWm/Gnm7NDPG/kCTDZEQQijqdP/Qt4OrgT0krdvRgYSwpJK06rL4U6Me\nugB7AV+j8YZIRVdgZiucDoBHgfULcW0DvAX8lnQhW8tKwArAO4ux33WB15xfWGX7WdtzcgyHS3pK\n0lhJv85p3fOd/RuAiZK2VKEnSNJ5kgbn5amSzszbPyppe0kPSPqnpP3LApK0o6THJY2S9HdJG5Vk\n/Xoue7Ryb1ax1yR/np5/91fqKbpT0guSfirpm3kfT0haO+eb3xsgaXo+nnGS7q/1/cn5VgWOAE6w\nPS8nnwX0krRZvTqu4VDgYdsP5LxvAycCP8r7Girpd5KeBE5R6jVZLq/bVNJTebmvpH/k+hkuaQPg\nm8BJSr1jvSVtktc9LekOSWvWOLZLJY1U6rU7oZD+TaXemkcl3VCs8xBC5xJvVm++D4CLgR8A0zo4\nlhCWVLOaztIpqerzbsAI25PzhWcv22ObKGOYpA+BzUl3txcvoNQY2oMFd9AhNT5uAu4Azpa0vO25\ned1Jko4GNgGusT11MXY/DHhE0ueA+4ChtidJ2hLYD9jJ9lxJV0v6IjAB2BI4zPbTkro3Uf4U270k\nXQb8Btgd6J6P7daSbSYC/Wx/JGlf4BTgmBr5ls9lf570b37/JmLpnWOfTepp+oXt7SSdBRwO/Loq\n/1rAX22fIOlq4ADgmhrlbgpMsz2/UWp7nqSxQA9K6rgkxh7AqKq0UTm9GFdf21YaRjYAuJfUGLpa\n0orAtcCXbY+XtKbttyT9Hphu+yIASXcBl9i+SdKJwBDg+Kp9n5S3XR54WNKNwDzgBGB7YG6O75GS\n4wkhLOWiR6RlLgUOBNbo6EBCWJLkngF3dBxLkEHAzXn5Jsp7H4oG2t6adFF7jqRuQFmd1qvr1SWN\nAf4NrAf8rbgP4E+2ZwFPkRpMFZWhWesC20vq10DMtYOzpwFbAKcDqwJP5N6FAcCOwIgc447AZnmz\nZ20/3eAu7si/xwGP2J5jezLpeMusAdwiaTxpyFOPknzX52P4G7CFmh4q9rjtN2y/C7wE/KUQW/ca\n+WfZvi8vjyzJA6lxW+s8K4VXWsfNKauYNqzSuwJcSRq6JdJ393rg08BU2+NJAbxVsq8dWPDdvwb4\nbI08h0gaTWpsbJHL3gG43/Z/cl3eVVJ+CKETiB6RFrD9rqQrgW8Af+joeEJYQq1Duju8TJK0Aumu\n/56SzgG6AHOAnzayve0XJL1KapBMZNEbH2sC0+sUMcN27zy0517gW8CFkrYHNgIeT9eXrEKaQ3JP\n1f5nSxoO7Mxi3JG2/R6pwXBHvqDdi9SzfJnt04p5cw9I8Tszl4VvmK1UVXxlCNK8wjIs2jNVdDpw\nt+1L8zCpoWWh10irF09x/8V45pHOfbVi/o9K8gBMATaS1DU3HCvzZ7YBnoHSOp5Qo6wJQHXDcrtK\nOVmx/oeThu/tBTxje7qkeo28omL9LdIAkrQJad7STrb/o/TwhpVY9NzVO5chhKVc9Ii03G+B+U+J\nkbS3pJeAnYD7JN1cumUIy4bZtt9dln6qjn8P0l36DW13t70B8IqkPo1UXp5XsAlpWM4sYIaknfO6\n9UmTyWtdbC4kx3U88MM8BGYQcGKOqTuwMfCFPOSmuP8upMnH/2wk3pJj2E55Pl1umH2aNKT1fmCQ\npLXyunUk/VeNIl4H1pPUTelJUHu0NJaC1YCX8/LgOvkG5dh2Bybl+RkvkoZgIWkP0jyeNpXP3x+B\ncwu9MicCT9ueUqeOa7kO2FVS/5x/ddJ8k/NK9m1Sr8blwFU5eRLQvTDXpTL3YyZp0nzFCOAreflQ\n4OGq4ruRhnC+k7/Pu+f0p4DdJK2mNLl+75JjCSF0AtEj0gz5P+3K8tsU/tG1/WcKk0FDCMu8QcBt\nVWm3AgeThptUbl5UVO5UD5P0PrAicLLt13L6kcAlklYj3Zk/tnKHvCm2R0gaR5owfyDwucK6d5Um\nhO+ZkypzRFYi3RG/pZF9lFgHuDxfIHcB/koa+jNP0pnA/fnieg6pUbBQY872B5LOBUaTLq7HLUYs\nFecCV0k6I8dTZl4eNmbgqzntFmBwrsv7gTdbIZ5GnEiaY/Jc7sV6ivQABCip41qF5F6u/YHfSlqH\n1Ntwke3q72nRDcB3gLtzGR9I+ipwhaSVSHWwG3An6bs7iHQujweulHQqqQF3ZFUsYyU9Q2rYTCX3\nutl+SdJvSA2ZN0i9Ne9Aevwx8Pv8fa65XL8aQwhLGi0YChoq8lCGyn/wXWvc6Qwh1BB/OyF0LkpP\nKetl+/vtuM9VcwN5ZdK7ub5muzUaoSGEJUz0iIQQQghhEZIuBf4f6eEC7el/81PAPkZ6V000QkLo\npKJHpIa4qxtCy8TfTvuT9ASLTuIeYLu9hg2R53rcX5U8x3bf9oqhWn7s7jlVyY/abpU3ibcFSRcD\nu1Qln5if3NWccpa48xFCCLVEQ6SGuJgKoWXibyeEEEIIjYqnZoUQQgghhBDaXTREQgghhBBCCO0u\nGiIhhBBCCCGEdhcNkRBCCCGEEEK7i4ZICCGETkdSf0mWtFVV2iIv+5M0XNIkSWMkTZB0aE7vnl/2\nWMw7VNI+dfY7N5czRtJTknoX1u0qaWTe1wRJA3P6UZKuqCqnj6Sn8vL0QjzOL5ys5BtWeVN6STzH\nSBon6WlJYyXtWFppHah4biQNkXRcG+zjm/mFiy3Z9tuSJuefhyVtU1jXcB1Lmprzjcl5dy/L28z4\nekvas+mcISxZ4j0iIYQQOqNBwGOkt8n/vIH8A22Pl7QuMAa4roX7nWG7N4CkrwCnAgdIWg+4BtjH\n9tOSNiS9Wf554FbgF5JWsP1hLucg4MYa5b8G/EDSFbbn1QtE0vrAD4A++QWBa7Hoo57rbd/F9keN\n5l/S2f59S7aTtB/pzfA72n5b0m7AbbmRuxbNr+Odbc+StD1wGbBdS+Kq0hvoCdyzuAV1tvMelmzR\nIxJCCKFTkdQF2Av4GumCvjm6AjNbKZTVgBl5+VvAlbafBrA9DfgF8EPbM4ARQPHu+FeAm2uU+TIw\nCjiwgf2vA7wDvJf3+abtV2B+j8vfcw/NnZLWzOlTJZ0taTSwW/7cNa/bR9LQvDxU0kWSnpD0nKSd\nJN2QewzOKgtI0mqSHpA0KvcKlPUI9MllPyvpoLztQj1axd4gSdMlXSDpGUm3S/qspEckTZG0U84z\nv6cl94KdI2mEpPHFnrMafkR6n8vbuR4fAB4GDqtXxw2Y//1Q8pscy/x6kbR1oa7GSFqnRp12AU4H\njsh59pa0o6TH87Z/l7RRdR1U6q1Qtw9I+jPwaL26DqE1RY9ICCEsJqX3p4QOUuN9NbsBI2xPzheo\nvWyPbaKYYZI+BDYHBi9GOKtLGgOsQrpbvnNO7wFcVZV3FPCdvHwjqdH0F0mfAf5t+8WSfZwNXE3t\nHpOisaT3+jwv6V7gBtsPSFoBOB/Y3/Zbkr4G/IR0wQ3wL9vbAkiqV343230lHQbcCWwPvApMknS+\n7ek1tnkP2M/2zNz79DegV418PUkvd1wdeFLSX5s41rWA22x/V9I9wHGkt8LvAZwE7Fdjmw9t98nH\n/wPg6yVl9yCdq6JROf1KatRxE7E+Jml5YCNSgxlSw3MzYBtgQ2C4pE8D3wB+Z/sySSsDi/RU2P5I\n0qlAT9snQGrwAf3yun2BU4Bjmohre2BL269EoyO0l2iIhBDC4pvVdJbQhqqvlgexoDfhpvy5qYZI\nZWjWxqSLwLuBsjf+1nsTcHFo1kDgYlJPh2psV0y7HTg3NxLKhmWlnac4/yVp73oHlC9Cdyc1hnYH\nrpV0CvAE6eL/gdzQWB6YUNi0Vk9MLXfk3+OA5yoNJ0lTgA2AWg0RkY6zH+miegtJK9bId4vtOcC/\nJY0kXaDXM8v2g3l5PDDZ9jxJ44DuJdvcmn+PJPVuNIcAl9Wx7T/U2bYyNOszwNWStgT6Adfl4XZT\nJT0LbAE8Dpyah3zdZPufDca3BnCNpE1Jo1/ebmCbR5vRmxNCq4ihWSGE0EKSVpVU76I0tLN8Ib8f\n8EtJU4Ef04zhWbZfIN3V3xJ4k3RBV7QmtS+wa7mLBT0iE1h0LsC2wDN5vzOBf5Du4B8ALDKpvspZ\npF6Mupw8avs0Uu/L/qT/+0fb7p1/etouTuKeXViey4Jrheq5D3Py73mF5crnLiUhHQasCmybG2yz\ngFoNkeq/K1fFUh1P9f6LsZXFUsnzUZ08ABNJ56qoeO5q1XGTbD9J6vH5RI3VlYbOdcA+OdZ7JTU6\nn+R04G7bPUnf/0pd1avDsvNenS+EVhM9IiGE0DrWYeH/yEPH2AN4xPb8i0FJD0nq08jGktYGNgGm\n5bvWMyTtbPsxpcnfW7Nw70E9OwOVO9i/Ax6XNCz3aKxPakgcXMh/I6mB8S/bL9cr2PY/8lCyvsBF\nJceyHrCO7TE5aWtgGjAJ2EDS9rZHSloJ2Nj2pBrFvAj0lvQw8GVqDA1qptVIw87mKj19bK2SfPtL\nOp90ob4Xn0csAAAgAElEQVQdqddldaBHHtZUGfZ2/mLG04jzgbMlfcH2DEm7koZ9/U+dOm6SpE8B\nK5AavI+Q5njcQOpN2gx4VtImwPO2f53z1xomBmleU7fC59VI84lg4aGGLwKfy/vfgzQnqpZpdExd\nh2VMNERCCKF1zK4xVyG0v0HAbVVpt5Iu+O8C9pb0UmFdv/x7mKT3SXfnT7b9Wk4/Ergkj7mfCxxr\nu95QvMocEeX83wCw/bKkwcBVeU7RR8DPbBcvKu8E/gCc2OCxnkX9pyStAPxG0ieBD4GpwNdtf6D0\nGNsLJHUj9QacQWqgVDsDuJw0qfoJFr7YbYlrgbvykKlHKL9on0CaEL4G8GPb7wDv5MnUE4HJwOjF\njKUhtm/NDccncw/oG6T5NbMlfYIaddxEkY9JqvTUHJWHd91C+i4+TfreHGP7/XyevpobnS+yYDhZ\ntQeBk5QeMnAycC7pu3YGUJxfcwswONf//aRGUK1jnlZW1zn96DyXpOZyE8cfwnyyY1RBtfyfROU/\nmq5xcRFCY5a1v51l7XhDCCGE1hRzREIIIYQQQgjtLoZmhRBCCM0k6QkWncA7wHbNoS7tEM+twMZV\nyYfbHtdB8axFGvpTNMd2346IpxGSjgK+W5V8s+0zW1BWm5wPSVuTXoxZ9EJxTlQIS5MYmlVDDLcI\noWWWtb+dZe14QwghhNYUQ7NCCKGVKb2N+o+Fz8dJGpKX91F6e/I8ST0LeSTp90pvgh6Rn/8fQggh\ndFrREAkhhLaxS345XrXJwEDgoar0fYC1bW8GDAHOadvwQgghhI4VDZEQQmgbF5BeprcQ28+VvK/h\nSywY+303sLPya69DCCGEzigmq4cQQg15/kdTVqmz7mrSewfWbXCX65FfQGbbkt4mvUis0bd4hxBC\nCEuVaIiEEEJt9V5a14gPgIuBH9DYm5Zr9X7E00RCCCF0WjE0K4TQ5pqYvP1DSc9IelrSrfkN1h1G\n0qr57cmt4VLgQNLboZvyMvDfOQblbd5qpTgWm6T+kixpq6q0YTXyDpc0SdIYSRMkHZrTu0saUZV3\nqKR96ux3bi6n8rOypMGSXq9KX71OPCtK+pWk5yWNzPH1LawfLGmOpNWbqIMtJT2U9zdJ0mn1a63j\n5L+5rrXqvJXKX0/StS3cdutcj5MlPSvp24V1DdexpCGSXip8B77fxH6n5999JP2yJbE3In8Xv9FW\n5YfQmURDJITQXsomb48EtrW9DfAMcEL7hlXXOkDXOj/r1Ns4P873SqCRi5K7gMPz8heBx71kPV99\nEPAYcFCD+Qfa7g0MAH61GPudYbt34ee9nH51VfqMOmWcA6wOfNr29sCRpGFvFYOAEUBT72K4ADgj\nH9dWwI2NHoSkLo3mXRrYfsX2Yc3dTtIqwO3AENtbAH2BIyXtl7M0t47PLnwHft1g7CNs/6i5sTfD\n6jT2Nx/CMi+GZjVtlZgvGkLD6s2ZqEze/p9iou3hhY9PkZ4etaSYXe/dIA3+2/BbYP5Fj6S9ST0l\nnwDuk/Sw7QNJDZF9JD0PzAAObnCeSqurPuZ8Eb0X8HngNuDnzSiuKzCz9aJrnlyHRwDdbX8IYPtF\n4MW8fk1gc1Lj5GekhmOZdYFXchkfARML+7iYdOG8HHCS7Xtzr9/GufzRkt4Aptu+KG833fbakvoD\npwDvAT2By0i9Yd8APgS+aLvmXKHcY7A38DHgPtu1egVWlHQTsDXwJPB123MlTQV62p6Ve6UG2h4s\naShpaOIOwJq5/r4LbAvcYvsnkroDw2z3kTQ4x7BGPt7f2T6/pA4PBR62/UCux7clnQicQWqg1Kzj\n5pD0edKT5z4GTAC+ZvuDwvr+wHG2B0r6HHAhMA/40HafkjK7AZcAvUhDJo+z/bCku4D/Ir3c8izb\n1wJnAj0kjSHV1+nNPYYQlhXREGna6x0dQAidRCOTt48Erm+neNqM7e6F5beBboXPfwbWr7HNPKru\norbiELHmqm5h7QaMsD1Z0nRJvWyPbaKMYZI+JF2ED16MWFbPF3TkGI7Oy0dI2j0vz7T92ZLtNwWm\n2S5rDH2F1Lh6DNhc0tplF/2kC9bHJD0E/AW4MvfQnAzclS/i1wYekbRl3mYzoL/tDyrDEUv0BrYE\nZgMvAL+wvZ2ks0g9ZWV3+y+w/fM8nG+YpF1sP1qVpyfpYnyEpOuArwJD68QC0M12X0mHAXcC2wOv\nApMk1WpkbAP0IV1XTJb02+LFf0EPYFRV2qicDuV1XOYkSZXvxNHAVFLDfzfb70k6HTiG1FCs5QfA\nD3LD8eN19vMz0vfocEnLseBv+gjbb+XG6FN5aODJwBZljZoQwgIxNKsG2+/ajm6QEFpXcfL2IiR9\nF1jOdsPDXTorte48ldYwCLg5L9+UPzdloO2tSRfX5+Q7ymXHVO9Yi0Ozji6kF4dmlTVCIDWq6pU/\nCLg5D4O7DTigLKPty0kX9XeS5v78La/aEzg1N5juA1YFPpnX3V5yQV7tcdtv5N6ol0gX4QDjgO51\nthsg6UlgLLALCy7oi6bYrswTuQHo10A8dxT2/5ztF/NxTAE2qJH/Xtuz8hC5V1hw/NXKzoehbh2X\nKQ7NGgHsSGoUPZ7Px4GkXpoyjwJnSzqe+j26uwO/yzHOs/2fnP59SWNJDdkN808IoUHRI1Jf144O\nIISlzCrU70W8FBgPXFFMlPQl0vCPXdsutKXWOqS75B1C0grAfsCeks4BugBzgJ82sr3tFyS9SmqQ\nTGTRiftr0raPKJ4CbCSpq+2FnoQmaR3SRfmNeZjdiqR5SpeWFWb7X8Clkq4A3sg9IMsBX8pDvorl\nw8Lnbi4L3wBcqbA8p7A8r/B5HqnOFyHpY8BvgD62X5V0XlWZ88Mu+VyMp3q74v6rY6sVTzHPR2Ux\nk4ZKVTeEtiPVewquRh3X6aWqthxwt+2jGsls+2xJfyENCf2HpJ1tv9zItnlY1y5AX9vv54cCrEQa\nThdCaEA0ROqoNzY8hLCopuZM2H5XUmXy9h/yNtsD5wEDqi8UA9DEPJV2sAfwiO35E7nzU40aGnaS\nL9Q3IQ1rmSVpRr7Ye0zS+qR5CxPaJHLmf+euAX4t6X/y3IgNSXfdu5PmM8yfVyHpn5I+afvfNY5l\nT+AB23NJQ74+Is3nuQc4Hvhhztfb9pjq7UnzUj6X8+zB4t/s+hipUfFmHlb0ZVLDpNrmkrazPYr0\nsIF7C/H0lvRw3vajxYynEdcBJ0vqb3u40pPKzgLOhrp13KjHgQslbWT7RaWn8K1l+4VamSVtmocZ\njpXUj9R7Uqshch9pftvJhaFZqwFv5kZIb9L8EUhzorrVKCOEUCWGZoUQ2ttvSU+VqTiH9B/6XfkR\nnGVjuUPHGEQaslR0K3BwXt47P0K18tM9pw/LQ2MeAk62/VpOP5I0FGZMLvfYFjZAj9DCj+/dtE48\nJwLvkuYujCO9wX56ybFVhgTVshcwIQ/FuZ40P2AuaaL1x5UeQT2R8ie/3ULqnRlHejLamy047vny\nMKirSD1NdwL/KMk6HjhR0jOkHo3rcvoZwOXAEzTvYr/FbM8mPZ3sDEmTSZPnr7NdOQ9lddxo+W+Q\n5oT8SdLTpO/fRnU2+b7SI6afJg0pe7wk3xlA93zuRpEa0H8FuuVzfjLpCYDYfhMYJWmcpFOVHnX8\nZ5j/2ONFlkNYVslL1NMhQwhLszxhs3JR2XVp7FVszjG01fGWlav0lKNHbH81fz4OWNv2EKWnHp1N\nmiOwje3xOc9nWPC0n/1t39UaMYYQQgiLK3pEQghh6VL2PpbJwEDSHeCiV0hPE1rqn0YWQgihc4k5\nIiGE0Dpa851DLXkfy3Ow6Dwd2y8BL0ma11rBtQVJT7DohOkBeZhLe8dyFOm9GUU32z6zvWOpyEMW\nd6lKPtF2U0+V6hCS1gLur0qeY7tvrfxNlNVu50PpHSTnVCU/avvbtfKHEBZPNERCCKF1tNc7hxp5\nH8tSpyUXqG3F9pXUf7Fhu1vaLoRzA7J3K5XVbucjN+yWyMZdCJ1RDM0KIYQW6qB3DtV9H0sIIYSw\ntIgekRBCWHyt/c6hFr2PJYQQQliaREMkhBAWU2s/Hawl72MJIYQQljYxNCuEEJZOC72PRdLekl4C\ndgLuk3RzTt8mpx8IDJVU9p6EEEIIoV3Fe0RCCK1mWXuPSGeOIYQQQmhr0SMSQgihXUjqL8mStqpK\nG1Yj73BJk/Jb0ydIOjSnd5c0oirv0PxCx7L9ri/pFknPSxoh6WZJn6za/skaZb4gafn8uaek4c2I\neYykz5bE01XSjZU3sUv6S1nsHa1Yt5KmSmrt+VBI+rOklVuw3RqSrpM0JZ/biyWtlNc1q44l9ZP0\nuKRnJI2WdH4Lj2VfSd9v5jZ1v78hdGbREAkhhNBeBgGPAQc1mH+g7d7AAOBXLdmh0oSb24C7bW9q\nuw9wIfCJvH5F4HPAypI2qdp8eeCQZu5yoO3e+efhkjzHA1Nsb2O7B+m9MI0eT5dmxrPEs7237fda\nsOkVwATbmwGbk87X2Xldw3UsaT3SCz+Pt70lsD3wbAviwfYdtn/dkm1DWBZFQySEEJZsq0hadWn7\nqT6IfAG9F/A1Gm+IVHQFZraw/gYA79qeP6nf9sO2x+ePXwD+TroQrY7r18CPcmOmNa1LeuN9JZ5x\nlWVJJ0p6Kt/JPyGn9Zf0gKQ/A49W98hIGiapf16eLumCfGf/dkmflfRI7jXYqSwgSftJeiL3Btwt\nafWSrKdIGpfLXC9vW+w16Sppal4enGO7P/emHClpSD62ewq9F1Pzdt0ljZV0VY7/xrK6l7QZ0As4\nK9fhPOAE4Ij8/Sut4xq+BVxh+6lKWbb/r1695GO+QNI/JD0nadfCMZ+Xlw/OxzFW0u119l88rvm9\nTpL2kTS0sL9LJI3MZe7aSHkhLOmiIRJCCEu210nzRZa2n2q7ASNsTwamS+rVwLEPkzSO9KjinzWQ\nv5YewKg66wcBNwE35uWiZ4HJwH7N2N8w5aFZdfIMBYbki/lTJa0PIGlPYH3gM8C2wN6SeuZttgeO\ntr1jE/tfC7gt39lfGTgO+H/At4GT6mz3ELCj7W1JL/Qre4HiG7a3Bq4BftFELJDq/0vAZ0nvv5lg\nexvgTWDvGvm3JL3ZvAfwSaBfnXLH5gYIALZnAtOATSmp4zpllX1H6tXLWqSHQxwLnFpj25OBfW33\nAo6os/9G/TfQB/gKcFkbNJBDaHfREAkhhCVMB70osa0NAm7Oyzex6EV/LQPzRe+WwDmSugFlT1gp\nS1fZOkkfI12k32P7eeAjSVtUZfsF8JMGYi3G3DsPKasdqD2KdLF8IbAJMFrSOsCewBeB0cBIYCPg\nU3mzR22/UqO4arNsP5iXxwMP5Iv1cUD3OtttCNybG37Hky7Oa7m+8LuskVB0v+3Ztv9FehnnHTm9\nLJ7Jtic6PUlndJ2Yy86rANep4+aUBfXr5dYc58iSOB8FLpV0dN7H4rrRyURSY/+/W6HMEDpUvEck\nhBCWXK0+MbgjSFqB1Kuwp6RzgC7AHOCnjWxv+wVJr5IaJBOBNaqyrAlML9l8IrB/ybov5rKezTeX\nVyMNzzqjsO/Rkt4mDfFqNbbfITXIbpJ0F6lBtBxwmu2rinnzsKvZhaS5LHwjcaXC8pzC8rzC53mk\nei9zIXCm7XvyMKvBZaEXfleWi/GsVJV/oXhsNxVPMf9HdWKeCPSWtFylVyQ3VDcA/gmldbzIQwZy\nWduxoJFUVK9eKrGWxfk/wI6kHqERkrZuYC5MvbqsbizFY0/DUi96REIIYQmVe0aWyp+qQ9kDeMT2\nhra7294AeEVSn0bqQdLapLva02zPAmZI2jmvWx/YGphQsvl9wGqSBhfK65eHPA0CDssxdScNe6nV\nU/ML0vyDViFpl8Jcg1XJxwbcAxwtaZW8rrukj9coYhrQQ9LySk//2rkVwloNeDkP96k3jGhQ4fcj\neflFoNIDdEArxNIk28+RelV+DCBpOeBc4Jr8ws+yOq7lEuDrkravlCXpmLyu0XqpZRPbj5OGaH1A\nGsrVlBdJDSwBX65ad5CSLUk3KRrpIQthiRY9IiGEENraINKTq4puBQ4G7iLNhXipsK4y5GeYpPeB\nFYGTbb+W048ELpG0GukO8rG5gbII25b0ZeBCST8D3icNWfox6WlZRxby/lPSR4V5GZX0hyS9WFV0\nWcyN2BT4vSST7qQPtf0kgKQewD/yheh/SPMBqo9pmtLE9YmkOSyjm7HvMqcBdwJvkSbvb1SS75N5\nmNJ/WDC5/3LgdkkH5DLay1HAxZKmkG6s/o0F82BK67ia7ZclHZbLWp3UW/PnvLrReqnlvDypXqRh\nXC81tQGpN+5yYAbwBNCtsO5VYASwCnBM/m73Ab5p++iy5WbEG0K7ixcahhBajTrBi/g6wzGEEDqX\n/PSsYbbv6uhYQmhNMTQrhBBCCCGE0O5iaFYIIYROQdITLDrBd4DtNzsonspjbotesF02eb7NSToK\n+G5V8s22z+yIeBoh6VZg46rkw5t4N0itcjr0fCxO3dse3CZBhdDBYmhWCKHVdIZhTZ3hGEIIIYSl\nQQzNCiGEEEIIIbS7aIiEEEIIIYQQ2l00REIIIYQQQgjtLhoiIYQQQgghhHYXDZEQQmiApKmS/lj4\nfJykIXn5B5KeljRG0j35bddLNUn9JVnSVlVpw2rkHS5pUj7+CZIOzendJY2oyjtU0j519mtJlxQ+\n/1d+yeCQqnwvSDq3Ku3jkq6U9E9JIyT9RdKnchyW9ONC3r45bXDhGHoW1tc81sL6XSQ9KWmspGck\nHVuWt6NJmp5/1z2mxSi/j6RftmC7/pJm5O/NGEnXtnZsrUnS9ySt2ESeDSXdLek5SVMknZ5fTln5\nLv8lH+skSZe1T+QhLLmiIRJCCI3bRVL1Y0QBLre9je3epLcw/7Sd42oLg4DHWPD27KYMzMc/APjV\nYuz3LWBHSV0q5QITihkk9QVeBw6oXORlQ4F/AZva7gP8CKg0Cp8D9i3kPQho1uNfq/yB9AjZXkAv\n4MFGNywcW6dge4TtH7Vw8/ts984/h1WvXMLq6ntAaUMkfxdvBa6zvTnQE9gWOD5nOQ24JR/rp4EL\n2jjeEJZ40RAJIYTGXQD8uDrR9juFj6sA7fZcdEmrLu5PjTK7AHsBX6PxhkhFV2DmYhySgYeBXfPn\n/YFbqvIMAv4P+CfQF0DS5qQGwRDn59LbHm/74bzNO8B0SZvkC8b+wAOLEedawBt5Px/YfjbH8QlJ\nt+QemcclbZvTh0r6naQngZOKPUOSukqampcHSxom6f7cC3ekpCG5x+0eSdXvSZlP0qWSRuZeqRNK\nsq0h6a+SJhd7MSq9Jnm52Ns3XNL5udynJW2X7/hPkXRczjO/pyXHermkh3LP1MHNrdi8z99IGgkc\nLumbkp7KvU/XSVoh59spH+toSRcXYhgq6SJJT+SeiZ0k3ZCP+azCfg4vlPvrnNY9f75KqafrRiXf\nBtYDHpN0R0noA4BZtq8FsP0+qRFSaaStC7xSyWx7fHPrJoTOJl5oGEJoK6ssfLN6qbFKnXVXA09K\nWrd6haSTgP8hvYNk1+r1bWhW01maVH2idgNG2J4sabqkXrbHNlHGMEkfApsDgxcznptIF6CTgA+A\n6cDaMP+u85eB/81xDwL+AWwJjLU9r065N5MaVn8HRgJzFyPGS4Apku4n9YJda/sj4DfAWbafyo2j\nP5IbS6TGS1/bljS0Ttk9gD45/zPAUbaHSLoe2Jt0172Wk2y/JWl54GFJN9r+V1WeHUl19QrwgKT+\ntoc3cayzbG8v6WTgRmCHnP4McFGN/JuSLso3BP4G3FCn7N0ljcnLv7f9+7z8oe3tASStWUmX9CvS\nObwWuJzUKzVKhWGTWTfbfSUdRjo/2wOvApMknQ98AtgP2Mn2XElXS/oiqfdtS+CQfHwPAv1sXyzp\nR8DOtsv+5noAo4oJtl/IDf7VgN8D10l6KtfLFbbfqlM3IXR60RAJIbSV1zs6gDbwAXAx8ANgWnGF\n7bOBsyX9APgO8PO2DEQLv3ixtQ0iXbRDahQMAppqiAy0PV5p6NpwSXdT3jPUVI/RY8BvgYOBYcDH\nCut2ASbnC+5bgVNynauBcu8A/koarnUT8IUmYiotz/bPJV0HfJ40ZGcP4HBgd2CrQiN8jcJmwyq9\nNU243/ZsYLakD3LckIaSda+z3SGSjga6AOsDnyYNVSt61PY0gNyD0A8Y3kQ8xf2PsD0jbz9T0ho1\n8t9l+0PgeUmrN1H2fbYH1ki/ubDcS9IZwMfzz3u53OVtVy78bwSOLIn5Odsv5pinABuQvkc7AiPy\nuVqF1DidQPp+Tcz5R5Pq/GGaVvc7aPvPuXG6F3AAcKyknrbnNFB2CJ1SNERCCK0mv4VcktptaFIH\nuBQYD1xRsv5a4H7auCFSZR1gdmsUlIe97AfsKekc0kXtHBqc95LvAL9Kuqs8kYUvxAHWJPVw1CvD\nkh4CTmLB3emKQcAOlaFMpGPvl/fVS9JyZb0itmdKep3Uq3ACCzdE3qyKtZE4JwOTc4NkamFVH9u1\neluK52guC4ZHVw+3Kl6YzitcqM4jnY9FSNoE+DbpDv9/ciOj1jCu6r9N10gvi2dedWwl8bTGhXWx\nrv4AfNH2M3k4WHcW7cWr/txUzMsBl9k+baFCpO5V+T+ipM5rmEAaSlgsb2Pg3crwTdtvkHpWr5Y0\nHtgaGFFdUAjLipgjEkJoC12X4p916h1YbmxdCXyjkpbvclbsC0xqsJ5ay2zb77b0p6qsPYBHbG9o\nu7vtDYBXJPVpJBBJawObANPyEJYZknbO69YnXXhNqFNExcXAibbfLJS9HGlY1hY5tu7AccAg28+R\n7n7/LA/fQtJWkvpVlXsm8JM8jKroIeCrhf0cQp274JL2ruwH2IYFPWQPkoboVfJtU1LEi0DvvHxA\n2X6aoRuph+ydXM+7l+TbRdL6efjWAcAjOf0/kjbKDdHSp5p1kFWBfys9seoQANtvA3MlVeqwVq9K\nPfcDgyStBSBpHUn/1cQ2M0n1XK/Mj0s6JJe5Emmo3nn58+ckrZyXP0kaevdyM+MOoVOJHpEQQqur\ncXG71GhwXstvWTABFeBESTuS7p7+C/hmG4TWXgYBt1Wl3UoaJnUXsLeklwrrKhf6wyS9T3qq0Mm2\nX8vpRwKX5DHyc4Fj64yxny83LJ6rSt4VmFJsnJCG4Jwh6bvAUaQLv+clzSY1Dr5bVe4Iat+B/j/g\nAkljSXfX76H+3IYjgV9Leg94jwXzYr4D/D4PkVoxx/d0je0vB26XdABpDsNisT1W0jOkRvBUFjQw\nqj2R970xcIftv+f0U0iT91+j/RvSTRlCOmevA6ML6d8ArpU0h2bO+bE9QdKZwP254TmHdA7r/dt1\nGfCgpGdt71u9Mvfk7Q/8TtLppJu91wMX5iw7kP4WPiR9x06y/WqtHUnal9SzdmrZcqPHGsKSTI0N\nVw0hhGVD1dyLrktqo2ppiTOEtiJp1cr3XtJFwCTbtSbPhxCWUDE0K4QQQghLo/2UXg44kTTM6Q8d\nHVAIoXmiRySEEAqWlp6GpSXOMpKeYNFJ0QOqhl11uPzI2gOrki+wfWVHxAOQnxZW/WLNw20vzgsa\n24ykrYFrqpJfsL1/rfxLg9b6/kr6PHBOVfKjtr+9OPGFsLSIhkgIIRQsLRf4S0ucIYQQQpkYmhVC\nCJ2I0tu4/1j4XHxL9gaS/p6Hs4yUtENpQSGEEEIbi4ZICCF0Prvk9xdU+zFwne3epKck/W/7hhVC\nCCEsEA2REELofC4gNTqqmQXvQVgNqPno0BBCCKE9xHtEQgihBfIcjY60Sp11VwNPSlq3Kv0XwL2S\nvke6EbVTWwUXQgghNCUaIiGE8P/bu/dwq6p6/+Pvj1iW4t3MzAtqpngBVPx5LfF60kyzQ+60VOqU\n1tFjWpal1SHLUvNGJZZWkpqlUhhqRyuMFPWoIHJVvCSRqXlQUfGCAt/fH2MsmMy95tprb2DtDXxe\nz7OevdZYY4455pgLnd85xpijazpclK8bvUlamfxLLF7xG+BY4PKIGC7pw6THnVatwG1mZrZceWiW\nmVknSFpL0orwuMErSI+dXb+Q9h/AjQARcSswoFEBkgZJCkk7ltJG1sk7VtIjeSL8NEnH5vQ+ksaX\n8o6QdHiD/c4ufb5Q0pDCfnaqVx9JQyUtKPYESbq8fL4kDZE0T9J6pbQFkt5fVY9SGb0kDZc0RdJU\nSfdI6l2Vvzvldjklv1+i/ZbhPn4maZsubLeGpMskPSHpcUm/lrR+/q7pNq74nY2UNKhLB9SNJB0h\n6fTurodZKzgQMTPruo2B3t302rhRxfLjfK8CTiwk/wM4EEDSnvlzI23APcDRHeSrGZwnwh8IXNzk\nNsvadGAwpAtZ4P8BL5TytAHjgfI6Fk8DZza5nzbSfJt+EbETMAR4q5kNJa10/++NiM9GxBNd2PQ8\n0uiMbSPifcBU4Bf5uy638YosIkZHxCXdXQ+zVljp/mNoZtZCr0XEq93xAl5ron4/AtYrfD4D+KKk\nSfm7E3MPT7v5Lvki/lDgMzQfiNT0Bl7p5DbLym+Bj+X3g4C7SJP0AZC0AbAtaTJ/W51t95W0eRP7\n2QR4NvJiXBHxaETMy/s4TtIDkiZJuiSn9cl39n8DTJfUt3gHv9TrM1PSuXn7uyXtJukOSX+TVLkI\noKQ9Jd0r6cH8mOYtK7L+Ry57Yq23q9hrkj/Pzn8HSfqzpJslPSnpLEmfz/u4T9JGOd+inhZJs/Px\nTJE0pmo+VU4/HjgjIhbm5O8D/SW9r1Ebd5akQ3Jv3VRJF0tSoa7DJD0s6feSPiBpXO6d2Svn2Sgf\n/+R8nH0a7Oc0STNy+16e047MbTVR0q21njilnsHhSo/SfljSfjl9iKQLu3KcZisaByJmZiuRiOgT\nEXPz+xcjYu2IGJo/T42IvSKif0TsHhEPkOa61JvvcgAwPiJmALMl9W9i9yMlTSHd1f7mUhzGevmi\n8fGv1PsAACAASURBVCFJD5EuVps1G3hT0qakAOqG0vf/DtxE6unZtnYhnc0Hfgh8uYn9jATaJI2X\ndJ6k7QEk9QWOBPaKiP7AhkrzcQD6At+LiO2B1zso//G8/XTgUlJQeCgwtME204F9I2JX4CLSI5rr\nWT2X/TXSXKKODCAFpDsBpwHvyPu4AziuTv4NgdsiYmfgnywODMu2AWZFxKKgNQckk4AdqGjjBnYo\n/W4OApD0TuBK4KNAP+D9LO4N2xC4KSL6Au8ETgE+CJxMah9IbX5XRPQDLif9Rqp8C9i10L4AdwJ7\nRsQuwO257Jr3AgNJv8srawGS2arCgYiZ2SpIHc91aSPPJyFdzJd7D+oZnC8++wLnS1qbQm9ESaN9\nz4mIAbUX6SlgjbYrp/0213d34N7Sd23Ajfku+020v0j+OXBEKUBpv8OIWcB2wDnAWsB9uXfhQGBP\nYHy+GN4TeF/e7NGImNyo3ILR+e8UYFxEzMtB4aYNtlkf+J2kqaQhTztU5Pt1Pobbge2aGCp2b0T8\nX+6Jewr4n0Ld+tTJPzci/pzfT6jIAyDqn0+l6lW2cZXppd9NrQ7bATMiYmYOdK4DPlCo61/y+6nA\nHTlP8dj2BWqLhN5AGu5X5X7gWqU5UrVhZFuQnlY3BTiVJc/L9ZFMJ90QeG+Dss1WOn5qlpmZLTHf\nRNLbSHf1D5F0PtALmAec1UxhEfGkpGdIAcl0lpwwD7ABqeeiK54vlVevrFF5v7+KiKjdZJa0Memi\n8vqc9nbgYdLE/lrd35B0BfDFjioSEa+TAobR+U72oaQnll0ZEd8u5s3DeYrD6eaz5M3ANUrF14Yg\nLSy8h3SRXuUc4NaIuCIPkxpRVfU6aY3qU9x/sT4LSb+NsmL+BRV5AB4HtpTUu9aLl4OifqTzUtXG\n0yrKa1YxAOrssUHjIPrDpCGBHwVOJwXDPwTOjYg/Kj2kYUiDslaEB2GYLTPuETEzs/J8k4NJd+G3\nyEO9NgeeljSwmcJyb8LWpGE3c4E5kvbO320G7EzXLybvBD6Vy1oNOIY0D2SRiJhNmnQ+vLTtYNLj\ni/vk16bAVpLeXco3PO+jHBwsImlX5adz5cBte9KjkseQhhNtmL/bWNJ76hTxHLCppLWVngR1cMeH\n3qF1SEOhYMmL3bK2XLeDgEdyD8DfyU9Rk3QwaZ7PcpV7WK4FLij0ypwJTI6Ixxu0cWfNAN4vacu8\nn09Q+s10YBzp0deQfkP318uUy948IsaQ5mNtoTTXah3gnzmQKg8zPFpJX1KbP92Jepmt8NwjYmZm\nZW2kYUtFo0gXcLcAh0l6qvDdvvnvSElvkHoazo6IZ3P6CcBwSeuQ7ryfVLsD3gU/BYYpTbgX8Efg\nN+VMEXFVnW3bSGP4i24mPeZ4UX0i4mVJ19G4B2hj4Gf5ArkXcBswMiIWSjoXGJMvTOeRgoJXS/V7\nU9IFwETSxfWUBvtq1gXALyV9J9enysI8bCzIQR3wO2BIHj40htTz1ApnApcAj+VeqgdI81Ggoo07\nu4OIeF3SicDvSdc9f6T977uRocAISceTnsA2pCJfL+BXeUiigHMiYoGkb5N+Zy8AfwWKDxF4hvQE\ntzWBzxV78AAkfT4fw0+q3nfiOMx6HOWHUZiZGYue5FO7KO2d79o2/X2rLG09qraXNJPUG1LrdTgF\n2CgihubhRSNIE3yfAj4RES8t3ZGYrZokjSAFr7d0d13MuouHZpmZWdk+kraqk34RaWjTzqQhNc2u\nuWFmZtaOh2aZma341uzCUz/XbPDdMNI6G18opfclDduB9OjWb9LkBPZ6JN1H+3kYB0ZEq4YFdSjP\n9RhTSp4XEXt0R30AJP0bcH4p+e6IOLle/p5A0mXAPqXkM/OTuzpTTo84H5LOJg3pKxpWMSSwrogY\nskwrZbYC8tAsM7OCFXRo1tIqD83alTQhd1/S5Nza0KzfAH+JiJ9KOhk4LyLWXkZ1MDOzVYyHZpmZ\nrYDyCuvLa/GzN0kL3X2plP5l0kT1B0mrXndLEGZmZisHD80yM1uxdfUxq2uSHiFb5QrSAm+/qCVE\nxD9J64vUHsP7oS7u28zMzIGImdmKrKtDwzqaUxIRr0q6CjiRtNp4bX2Q5/Oq5GdTWAjQzMysszw0\ny8zMqvwIWK/w+UBghqRHSYsg/rxbamVmZisFT1Y3MytYUSarL62V5TjMzGzF5R4RMzPrkKRBkkLS\njqW0ditdSxor6RFJD0maJunYnN5H0vhS3hGSDm+w3/m5nKmSbpS0Zk7fQtKtkh6T9Likc5THm0ka\nImmBpPcXypldKneIpHmS1iulPVeo9zV5Ve+quh2Z807O9TuyURt2l2K752O8cDns4whJp3dx28G5\n/R6RNEHSoMJ3Tbdx4Xc3SdLdkrbrSn2WJ0kzJXV1XpfZSseBiJmZNaMNuAc4usn8gyNiAGk418VL\nsd85ETEgInYiPc3r8zngGAVcFxHbAjsBuwCnFrZ7msYLLrYB44GjSulX53rvBKwLHFpv4xyg/Bg4\nJCL6AXsCk5s9KEm9ms27IoiI0RFxSWe3k7QL8H3gsIjYHvg0cLWk93SxjQdHRH/gl8APOlsfM2st\nByJmZrampLVqr/KX+aL5UOAzNB+I1PQGXlkGdQS4C3gfKbiZGxG/AoiIN0hByFcKeX8L7Ctp83Ih\nkjYAtiUt2thWsa/VgXcCcyq+XxsQ8FKuw9yIeDKXv42k2yWNl3SHpD45faykSyVNAI7Ln3fK3+0k\naWx+P1TSLySNk/SkpA9JulzSdEnXNmogSbfkXoWpkj5ZkW0rSXdKelTSaXm7JXqrJF0oaUh+P1PS\nuYWeht3ycf1N0lE5z6KeltzLNUzS/+Yeq/0aVPnLwPciYlZux8nAVcDJjdq4CXeSfitIOqTQq3Zx\noedsdj7OKZLG1H77kvbKPWITJV2m3Osnaet8ziZLGp1/R3VJ+rakB3LZ5QDtGzl9nKRNc/53Sfpd\n/s3cmwM0s5WeAxEzM3uONF+k9io7ABgfETOA2ZL6N1HmSElTSI8A/ubSVlDS6qRgaAqwA/Bg8ft8\ngbqWpHVy0nzgh6QL3bJ/B24i9fBsq/Q0sJrjJT0EPAO8GRF31qtPRLxAWuF7pqSraxfk2XDgpIgY\nCHyXJe/MvxURu0XEiA4OeUtgP+BTwEjSxfmOwNYdXKQeHxG7AXsAZ0sqr1wPsDtwBGnhys9L2qaD\nugA8nnsapgOXks7FocDQivwbAnsBJwHfalBuu3OZP+/QQRt35AhgiqR3AlcCHwX6Ae9ncS/YhsBt\nEbEz8E/gYzn9Z8BxEbELqVes5ofA8Nw7czfVxw5plfXd8z63kFRcVf7/8j6vAb6X0y4Fvp9/M8cD\nP+nEsZqtsByImJktA/mu8bWFz6dIGlrK819K8yy6fYx4JxdEbANuzO9voLoXoWhwvtjqC5wvaW2g\n6ukojZ6asl4ODMYDs0hP6lIH29T8HDiiFGhAPp78GOKbWHwBCouHZm0CvCnpU1WFR8QJwGHANOAC\nSd/J5/YDwE253pcA7y1sdmP7kur6Q0QsIAVer0TE/bm+U4E+DbY7XdIkUpC1RX6V3RYRcyJiLvAH\nUsDQkdH57xRgXETMy4HpphX5R+X6TuigvvXO5aK0em3cQT1H5nbfixSEbgfMiIiZEbEQuI50fiD1\nqv05v58A9FGaM7R6RNSCo+sLZe/O4vN3TaGceg6UdD8wCdiHFHDV/Lrwd9/8/iDgylz3G0m/P7OV\nntcRMTNbdvaRtFW94SOS3kW6oJrV+mo11DAoyuP0jwQOkXQ+0AuYB5zVTOER8aSkZ0gByXRg/VKW\nDYDZ7TZcbE4ODIp1mkZpboekrYBXI+LlPPKGiHhD0hXAFwv5NiZd/F2f870deJjSmigRMV/S7cDe\nQOVwqIiYCEyUNAYYQer9+Fe5zgWvFd7PZ/ENwXLPxbz8d2Hhfe1z3fklkvYnXfTukY99fC73rXK1\n63wu1qUz9akKZmt5FlTVN5tG6pkpzv3YhXROUuXat3GjHrbBETF1UeXaB6HFwKd4HLV6lo+n+DlK\n6XWDYUnvIPVwDIyIZ/KQtWJ7RuFvsYyBETG/XplmKyv3iJiZLTvDSPMO6vk+8N80dye/ZXLPyBKv\nUpaDSXfAt4iIPhGxOfC0pIHNlJ8vBLcGZuU78HMk7Z2/2wzYmXQx2hljgHUlHZPLWYN04VfvaVDD\nScObaheCg4HL87H0iYhNSXMm3l1n272Av1UcV29JHywk9cvH+DLwL0kfyfl6Kc8DqePvQC1g+VhF\nns5Yh7Tg5BuSBgBVQ+g+JGndPCfiQ8D/kobnbSpp7dyrc/AyqE8zLgbOyr8Fclt9Ghhe1cadLH8G\n8H5JW0paDfgEaa5RXRHxIjA/tx+k30vNeNKwPoBjG5TzDtK/8+clrUsaFlbUVvg7Lr//C/CFWgZJ\n/RodlNnKwj0iZmbV1lT7FcjXbJD/auB+SUsMq5C0J7BaRNxfp7yero00fKloFOmC7hbgMElPFb6r\nDTUZKekNUo/D2RHxbE4/gXSRuQ7pLvxJOUBpWkREni9wuaRzSDfVfk0aw1/O+7Kk61jcg9NG+zkL\nNwMfJ82POV7SQcDbSAHS5RXVEPD13OPyBvACaT4EpIvUn0j6bi7np6QhVWUXk3pmvgjc0eGBd+w2\n4AuSpue6T6jI9wDpmDchzXl4AkDSBcBE0sX+lGVQnw5FxIOSvgHcrvRQhNeAIRHxdB7OV9XGzZb/\nuqQTgd+Trnn+SPvfc9mJwK8kzSO1Ya2X4lTgKknfIgWRJ1Tsc46kX5J6AJ8mBXpF787zp15i8cMf\n/ov0m/ks6d/MaGCypM/nMn9S9b6JZjDrsbygoZlZgZZc6K8jixYClDST9MjXzwHvIV3MbQScQ7qD\nf2wepjET2KmzF99m1hqS1ir8u/4x8EhE/Libq2W2UvLQLDOzgk5O4q7nCtLd9dpciLVJAcq9OQjZ\nDJhWeLqTmfUstUUUp5OerPXz7q6Q2crKQ7PMzOqrmsS9Jmk8fV0R8aqkq0jDO34eES8B76p97x6R\n+iTdR/sJ0gdGxPPdUZ8iSZeRJoEXnRkRt3dHfaBnt1c9kv4NOL+UfHdEnNyFspbr+YiI60hP1+r2\nupit7Dw0y8ysE0pDt9oNzYqIuZLWJw3Nuigihpa2X5SvZZU2MzPrgRyImJl1QlUgYmZmZp3jOSJm\nZmZmZtZyDkTMzMzMzKzlHIiYmZmZmVnLORAxM7NVnqRBkkLSjqW0kXXyjpX0SH7E6zRJx+b0PpLG\nl/KOkHR4g/3OLn2+UNKQ/F6Svi3pcUmPSbpV0haFvPMLdbhZ0no5vZek4ZKmSJoq6Z68Wnpxm9rr\n+AZ1+1Yue4qk8ZK26qAZu4WkIZIuzO8btvdS7OMcSR/ownYdncOm27gz564nKf/GzYr8+F4zM7O0\n4vo9pJWu/7uJ/IMjYqqkTYCHaPJxr510KrAr6Slrb0j6JDBK0sBIT5qZExEDACRdA5wMnJuPZW2g\nX16F/v3AW7nMRds0ImlvYH9gQES8JWkzoOkHM0haLSIWNn+oPVtEfKuLm1aeQ2AvOtfGTZ07sxWJ\ne0TMzHoASWv51ZpXnbbvBRwKfIYUiHRGb+CVZfATqOcrwKkR8QZARPyK9MS2A+vkvZu0WCbAJsCz\nOVghIh6NiHmd3PcmwOyIeCuX8VREvAhpTRBJ90qaKOlaSW/P6c9L+rGkKcD7i3fCJZ0iaWh+P1bS\nRZImSJosaVelnoLHJZ1SVSFJ20i6S9KDku6XVHVRflgu+2FJ++VtF/Wa5M/jlXqw+kiaJOlXSj0W\nl0v6qKT7lHqTts35F/W0SJopaahSr8QDkt7ToB0bncPKNu4MSSfmut8j6Te1Nsz1rPWEHS5pRJ1j\n6a30SHEk7ZzbttbjsnHF/laXNEypF2eypKNz+hW53adJOqO0zWU5/VZ5MVcrcI+ImVnP4HVFWkel\nzwcA4yNihqTZkvpHxKQOyhgp6S1gW2DIUtRlPUkPFT5vCnw1X6ytGRFPlvI/COwA/LmWkAOpg4Ff\n1OoGjJO0f843IiIeqdjf6RHxlzr1+hPwbaXVxf8EXBMR4yVtRLq4PiAiXpd0DvA54DJgA+B/IqJ2\nIdzouOdGxG6SzgauB3bP6Q8DP67Y5hngoIiYJ6kfcFE+7rL3AgOBvsBNkrZrVJGc72jgcWBqrtse\nkk4CTgG+WGebpyJiQD7+zwLfKWdo4hxeRZ02blDPducOmAF8NR/vAmAiMK6D461yInB5RFwp6Z25\nvKp8GwD9I2Kh0rpJAF+LiBckrQ7cJen6iPgHaXX6OyLiZEnfB74EDO1iHW0l40DEzKwbacl1Sax7\ntAE35vc35M8dBSK1oVlbAWMl3QpULczVaMGuJYbbFO7aV13Fq1Be7cJ0M+Ax4HaAiJiVL74Pzq/7\nJO0dEdPK+6sSEa9I2oU0dOhA4E/5zvcaQD/g3hxorAHcmjd7PSJurVdeHaPz3ymkIHAOgKRXJK1f\n0TOwBvDjHIQsAN5VUfb1uTdouqS5pMCkkRkRMSPv/2EWB3lTgMMqthmV/04AjqjI0/AcVrVxRPyp\nYrt2507SR0kX+bX2G113y+bcC3xL0obADRHxt4p8BwEX14beFc7VMZI+C/Qi/Sa3B/4BzAN+l/P8\nGrh4KepoKxkHImZmPcfGwGvdXYlViaS3AUcCh0g6n3QRNQ84q5ntI+JJSc+Q7qpPB9YvZdkA6PRk\n3Yh4SdJrkvpExMzCV7sAt+T3c/Jd+bVId9T/E/hh3v510sX+aKWI4VBgWifrMD+X+6c8zOpI4I/A\nrRHx6TqblH+7xQBsjdJ3taFiCwvva597VVTpNGAm8Clgrfy+btXrfJ7PksPRi/Up779Yt6q61PIs\nqMrTzDmsaOOqQKRKVaBbPOY1OkqPiOsk3Q98JNfn4xHxYDMVkLQ1aY7SXvm4R9L+nHdUX1sFeY6I\nmVnP8VpEvOrX8n2V2vxgYFxEbBERfSJic+BppcnEHcpDlbYGZkXEXGCO0kRvlCYf70wnA4CCi4Bh\nktbI5R0DrAOMKWbKx3Qq8OU8fn9XpUn0tUBre2BWZ3YsaTtJ2+T3AnbKZdwL7C9py/zdOqp+0tNL\nkrbMdVgWT7JaB3g693YMaZDvaCV9SXN4ngb+DvTPdd4B6Gi41rJSeQ4btHFn3A8cIGndPB/kI4Xv\n/g4MyGV/tJye33+slpiDiSci4hJSwLlDxT7/DJwoabW83fqkhyPMBV7Ov/uDCvnXIAVYkIbAdXXo\nmK2E3CNiZmarsjbgplLaKOATpLvWh0l6qvDdvvnvSElvAG8Hzo6IZ3P6CcDwPD9gPnBSDlC6Yhip\nR2WapCDNYTgqX4gvIdL8jSmkC70XgJ/lAKAXcBtp3gi0n2fwy3zhWdabNAyqNrF4AvCjSPNCPgf8\nVmmS+kJST0V5HgTAN4A7gGeBR+p831nD834/x+KhUfU8A4wH1gQ+FxEhaRwwOw+9mkCai9IKlecw\nBw7t2rhBWXXPnaQfAA+Qet4mFL7/DvAzYA5wHylYIKf9XtLHgJsL+duAT+W5T3+nuo2vIAW3UyTN\nB74bETfmtn2E1FNVDDaeBw6SdG7+7hhIj0QmDcsbXfW+QVvYSkJ1/ntmZmYVSnM6ete5w96t5ZnZ\nqkvpyWSzI6Jqwr9Zj+KhWWZmZmZm1nLuETEz6wT3iFhXSLqP9pN3D4yI57ujPkU9rW6SdgauKSU/\nGRFHdUd9mpEfQ/zxUvKwiLiqC2X1iPOxIp4HW/E4EDEz6wQHImZmZsuGh2aZmZmZmVnLORAxMzMz\nM7OWcyBiZtYDSZop6drC51PyE3GQNETSc5Ieyq+2bquomZlZF3kdETOznmsfSVtFRL01Gq6OiDNa\nXiMzM7NlxD0iZmY91zDgq91dCTMzs+XBgYiZWdetKWmtpXmRVn+ucjVwsKRN6nx3jKTJkq6TtPHy\nOTxb1iQNkhSSdiyljayTd6ykR/Lwu2mSjs3pfSSNL+UdIenwBvsNScMLn98jaUFhuN9QSU8Vhvvd\nk9OLwwBnSDq9ap/FeuXtFkh6f+H72Q3q10vScElTJE2VdE9eebzHyW11Sn4/VtJOy2EfP5O0TRe2\nW0PSZZKekPS4pF9LWj9/13Qb53MZkr5aSNsjpw1pxbHYqsGBiJlZ1z1HevTu0ryea1D+m8BlwJdK\n6TcDW0dEP+D+nMdWDG3APcDRTeYfHBEDgAOBi5divy8Ae0rqVSsXmFbKc15EDMivvQvpV+c67At8\nQ9LmTe7zaeDMJvO2AWsD/SJiJ2AI8FYzG0pa6a5lIuKzEfFEFzY9jzTsftuIeB8wFfhF/q6zbfwY\ncETh89HAlM5URlKvpTgWWwWsdP94zcyWp4h4NSLUwl1eQVoobf1CHZ6PiHn545XA7i2sz3K1tD1M\nPelV59h6AYcCn6H5QKSmN/DKUjRtAHcB++XPRwG/61QBEf9HujjdtMlNfgvs22TgsgnwbOTFzSLi\n0dpvXNJxkh6QNEnSJTmtT76z/xtguqS+xV4iSRfW7twrPfjh3Lz93ZJ2k3SHpL9JqlycT9Keku6V\n9KCkv0rasiLrf+SyJyr3dBV7TfLn2fnvIEl/lnSzpCclnSXp83kf90naKOdb1NMiaXY+nimSxtT7\nbeV8awHHA2dExMKc/H2gv6T3NWrjCi8DsyVtLUnAIOCOwv62kXS7pPG5PfsU6n6ppAnAcaVjOS4f\nxxRJX8lpfXL7/VLSw5Kuz/tD0sDc9hNym22Q03+g1EM3SdI3GhyD9XAORMzMuqb3Mno1HFaVFzi8\nCjixlqYlh2p9lPZ3tldkS9vD1JNeZQcA4yNiBukCr38T7TFS0hTSne1vNpG/kRuAoyVtSuptKw+V\n+poWD826vrxxvtBcC5jU5P7mAz8EvtxE3pFAW76oPU/S9nmffYEjgb0ioj+woaQP5236At+LiO2B\n1zso//G8/XTgUlJAeCgwtME204F9I2JX4CKg6oJ39Vz212iud3IAKRjdCTgNeEfexx3AcXXybwjc\nFhE7A/8EPlZR7jbArIhYFLDmgGQSsAMVbdyBG0lB857ABNI5rRkOnBQRA4HvAj8ofPdWROwWESNq\nCZLeS2rv/Ug3T46RtFv+ui9wfq7nu0kB7NtI7X5UROwGjAK+LmlDUu9O39zuP2riOKyH8lOzzMy6\nYFmtgJ5v/HXkR8BXCp9PUxqbvwD4F3DSsqhLd9KSK8yvrNpIF3aQgoI2Or6oHxwRUyVtBYyVdCup\nd6OeqvSae0i/pU+QLkrfUfr+vIj4cZ3tjpf0b8D2wGkR8UaD/ZXTfk7qsfhuo4pFxCxJ2wEH59d9\nkvYG9iddBI/P/1bWJF0QTwMejYjJjcotGJ3/TgFm556AGTkoq7I+cI3S/IbVgBcr8v06H8PtSvNm\nOrrJe2/uXULSU8D/FOq2R538cyPiz/n9BKBPRbmi/jlRql79No6IRjcyRgO3kYKDG4AP5Xr3Bj4A\n3FTrvACK/028kfZ2B8ZExAu5jJGk4X6/B2ZExPScPjEf44tAf+COvI/VSef9pfz6haSbgFsa1N96\nOAciZmY9UET0Kbx/kTS2u/b5a6S7ryurjYHXursSy1K+u3skcIik84FewDzgrGa2j4gnJT1DunM8\nncJQvWwD2vdwlMsISXeSfjt9gWOarP7VEXGGpD2A2yX9NiKeBZ4v1aNdHSLiDUlXAF/saCcR8Trp\nwnd0HppzKKnn5sqI+HYxb+6dKf5G5rPkKI81SsXXhiAtLLyHdAFd5Rzg1oi4Ig8tGlFV9TppjepT\n3H+xPgtJv4uyYv4FFXkAHge2lNQ7IubCovkz/YCHobKNKwORiHhF0nPAYcAZ5EAkH9u/8tyhepr5\n91sMnOod42rAxIjYv92G0kDgEFJQ/SnSnCdbAXlolpmZ9TSv5bk4K/SrdEwHA+MiYouI6BMRmwNP\n5wuqDuW5A1uTht7MBebkHgMkbQbsTHND9C4DzoyI55s9GTURcR9wLXBqTroT+ISk2k3NT5HmoZQN\nz9+Vg4NFJO1aG3KYg7btgVnAGNJwog3zdxtLek+dIp4DNpW0dr5bf3Bnj6+OdUhDoSBN7K7Slut2\nEPBIHg71d9IQLCQdTBqGuVzl39y1wAWFXpkzgckR8XiDNu7IucDXI2JBYV8vA/+S9JFcXi91/PSw\n+4EDJa0vaQ3SPKVxDfI/AmxeG76l9ESw7fP5XTcibiY9yKMqGLIVgAMRMzOz5a8NuKmUNop0Rxfg\nMKXH59ZefXL6SEkPkS76z849EQAnAOfl724ijdXvcGhbRDwWEVdVfF2cI/KQpHfWyXMeMETSWhEx\nGngQeFDSJFLvyKV19vkycB2NL8Y3Bv6Q58NMAiYDI/OwoXOBMZImA7eSel7K+3gTuACYSLrj36mn\nO1W4ALhE0oOknpkqC/N5+AFQm6D+O1LvxBTgw6Teo1Y4k9TT8JikJ0hDmz6Tv6vbxh0VGBHjI6Le\ngw2OBf4rn/sppCe7NSrnaeDbpN/yBOD6iHiwQf43Sf9uhuV9TMjHszZwS077Iw2ezKb06OCBjd5b\n91J+eIKZmXWD0tyI3stq7smKxu1gZrbqcY+ImZmZmZm1nCerm5mZrSQk3Uf7uRgHdmVOyPKQ53qM\nKSXPi4h6T4tqifxEsPNLyXdHxMndUZ9mSLoM2KeUfGZE3N7Jcnrc+eiKZdUe1noemmVm1o08JClx\nO5iZrXo8NMvMzHospVWxry18PkXS0MLnM5RWWJ4u6fRuqaSZmXWJh2aZmVlPt4+krSLiyWJiflzq\n/sBOEfGWpIar1JuZWc/iHhEzM+vphgFfrZN+EvD9iHgLICKea2mtzMxsqbhHxMys51gzLXa8Slqz\nwXdXA/fXFmMr2BY4SNLFwBzglIh4dHlV0MzMli33iJiZ9RzPkSZsr4qvRr0Zb5JWBP9SKf1twJoR\n8f+Ai4Cqhfp6LEmDJIWkHUtp7RaakzRW0iN5scFpko7N6X0kjS/lHSHp8Ab7nV8o52ZJ65W+XwR4\nKgAAEiRJREFUHyrp2doK3ZJOLSx0OL/w/viK8o/M30+WNFXSkZ1rmdYotp2kIZIuXA77OKKr85ck\nDc7t94ikCZIGFb5ruo0Lv51Jkh6Q1OFq5JK+Wni/nqQTu3IMZo04EDEz60YR8WpErLLdIJ1wBfBx\n0urdNU+RVrAmP6Zz226o19JqA+4Bjm4y/+CIGEBaxfripdjvnIgYEBE7knqTyo+qHQw8CewHEBE/\nzPkHFLYdEBFXlwuW9Dbgx8AhEdEP2JO0indTJPXq2iH1TBExOiIu6ex2knYBvg8cFhHbA58Grpb0\nni628eCI6A/8lLQKfEeKwyHXA3pMILKy/UZWZQ5EzMx6ht6r+KvhRPP8ON+rWPJiaDRpsjqSdgdm\nNSqjp8kXU4cCn6H5QKSmN/DKMqrK3cBmhXr1A14AfkQKlDprbUDASwARMbf2oAFJ20i6XdJ4SXdI\n6pPTx0q6VNIE4Lj8eaf83U6Sxub3QyX9QtI4SU9K+pCky/NT065tX5XFJN2SexWmSvpkRbatJN0p\n6VFJp+XtluhxknShpCH5/UxJ5+aehrsl7ZaP62+Sjsp5FvW05J6qYZL+V9JjkvZrUOUvA9+LiFm5\nHSeT/g2c3KiNmzCOJc/3Fbldpkk6I6edC6yXe1x+ApwL7JA/f6vcaydpZE7rJenafD6mSPp0VSUk\nzc5tOUXSGKVHeCNpV0n3556eqyW9o9DW50maCBwu6epCWSdIurTJ47cexHNEzMx6gFV93Qw1Nzfm\nR8BXCp9/RrpDPBV4FfjcUtZhraXZviN1zvEBwPiImJEvyvpHxKQOihkp6S1S78+Qpa1TDoYOBn5R\nSG4DbiAFeudJWj0i5jdbZkS8IGkMMFPSn4BRETEqfz0cOCkiZko6gHRn/uP5u7ciYrdcryENdrEl\nqadmT+B2Ujv+J3C3pF0iYmLFdsfnuq0FPKA6w9+A3YF+wHxgvKSbgQUdHPLjEdFf0pXApcBBQB9S\nG46qk39DYC9SEP0tUu9WPTvQvufiQeCEDtq4Ix8mnduar+XyVgfuknR9RJwt6aTcA0YOGLeLiIH5\n86CKsgcAW0XEDjnfug3qsSFwW0SckYOKjwHXAL8EPhsR90m6nHRua71//4iIXXKP0KWS3h4Rb5J6\n8M5r8vitB3EgYmZmPVZE9Cm8f5F0J7j2eR5du2NfZW7HWZZKOdpqA27M72/InzsKRAZHxFRJWwFj\nJd0KVK1M3GjF4vUkPUS6M/4Y6YJ+0T6AQRExV9IDpAv9P3ZQryV3HHFCHlp0CHCBpF1Jq5d/ALgp\nB54iBZA1N7YrqL4/RMQCSVOAVyLifoAckPYBqgKR0yUdkd9vkV9vlfLcFhFzcnl/IAUM4zqoT+2i\nfgowO/8uZ0jatCL/qIiI3PvTp0G5ov05XJRWr40j4psNyhuZexfWAfoX0o+R9FmgF+n3sD3wjwbl\nNPI3YFOllc5/HxGNfjdzI+LP+f0EoE8OXNaIiPty+jWkmw+1QORGgPy47rHAIZLuJAVt93SxztaN\nPDTLzMxWaZLWktToon157PNtwJHADyTNJI3Hb3p4Vh6G8wzQF3ieJefOAGwAzG5QxJx8t3tL0oXt\nf+Z67ZbT7s31+gBdDPYiYmJEnA8cAxxFuub4V2F+Sf+I2LuwyWuF9/NZfI2yRqnoefnvwsL72ue6\ncwck7Q/sA+yR50k8UqdcaH/hH6W6dKY+Vd18tTwLquqbTQN2LaXtAjy8qHLt27iRwcBWpKfQDQOQ\ntDVpqNegPNfkDuq3S1ndNsk3C3YG/koK/BpN/i+2Va0tym1WDsaKv5EbSb1pHwFujoiW/hu2ZcOB\niJmZ2WIbs/zmwRQdDIyLiC0iok9EbA48LWlgM5WUtBGwNTArIuYCcyTtnb/bjHQxOK2jcvJwsVOB\nL+ehOW3AmblOfUgXrh+S9PZm6pX331vSBwtJ/XI9Xwb+JekjOV8v5XkgdfydNMwH0pCdpbUO8HxE\nvKH0xKj+Ffk+JGndPHzrQ8D/kp7otqmktSX1Jp27VrgYOCufT3JbfRoYXtXGHRWYL9bPAvaWtD2p\nh3Eu8HLez0GF7Au0eFL4KxR6I/O+dpC0uqR3A7Xf3kbAahFxAzCUxeewKbk3ap7SnC+AY4G7KrL/\niRQof5Lme9Osh/HQLDMzs8Vea9F8nTbgplLaKOATwC3AYZKeKny3b/47UtIbwNuBsyPi2Zx+AukC\ndR3S3eqTcoDSoYgYn4c5HU26w7x/4btXlSZqH5Lr1QwBX5d0BfAGaeL7Sfm7Y4GfSPou6fHLPwWm\n1injYuB6SV8k3aVfWrcBX5A0nRSgTajI9wBwM7AJMDwingCQdAFpyNcs0hCs5S4iHpT0DeD2HBC8\nBgyJiKclrU11G3dU7mu5p+JLEXGipIdJPUQzWXIY2i+BKZLujIjPS3ow/05ujIhz8tC16cAMFg+H\ney8wQumxz/OB07pw6EOAy/MwsoeAyyuO4y1JfyUFhouGZUl6qDC3pe576znkniwzM+tu+Q507cK5\ndysn71ftOw9NGhcRn8qfTwE2ioihki5h8QX7usBLvsgxM+sc94iYmZlV20fSVuVHo0bEogXq8t39\npp8qZWZmiQMRMzOzasNIE8m/0CDPx4GPtqY6nSPpPtpPPj4wIp5fRuVfRpoEXnRmXmCyWyzvY17W\nJP0b6YliRXdHRHmRyWbK6onnYxRprlHRcRHRkiFu1rN5aJaZmXW70vCojVny6TjL25qkCcnQfmjW\nrsD9pDkag8lDs2ob5onPV+cnDpmZWSe4R8TMzHqa5zrO0jJvApcBX6L+U4mOBq5vaY3MzFYSfnyv\nmZl1u4h4NSKaWl69G1xBGn5VXqsDUiByQ2urY2a2cnCPiJmZ9STl9TZaoTg0q538CNurgBOBn9fS\n8+J/L0XEY8u/imZmKx8HImZm1mO08rG9NVJTHTE/Ar5SSnNviJnZUvBkdTMzW6V15xomZmarMs8R\nMTMz66EkDZIUknYspY2sk3espEckPSRpmqRjc3qfvDp6Me8ISYc32G9IGl74/B5JCyQNLaSdLGlG\nft0lqV/hu9ml8obmBSGRdGSu42RJUyUdWajT3/J3D9U7xkJ5+0i6X9IkSQ9LampV8e5Qa4uq87YM\nyh8o6Qdd3HY/SRPy72aapMGF75pu43zuXpT09kLa/+Qnz5lV8tAsMzOznqsNuIc0DOy/m8g/OCKm\nStoEeAi4rov7fQHYU1KviFhAenTxtNqXOXg4AdgzIl6UdABwk6QdI+L1qkIlvQ34MbBbRDwnqTfw\nrkKWUyPilibq93PgyIiYkS9++zR7YIVjWilExHhgfIcZSyRtClwDHB4RkyVtAYyR9ERETKTzbfwU\ncAhwi6QNWPK8mtXlHhEzM7PF1pS0Vne8yhWR1As4FPgMKRDpjN7AK0vRDgHcBeyXPx8F/K7w/VdI\nC+W9CBARd+T8n+yg3LUBAS/l7eaWV61v0obA/+Uy3oyIRwEkvUvS7ySNl3SvpF1y+ghJl0u6H/ha\nsUdIUu/anXtJQySNlDRG0kxJJ+TenMmS/iipvFDiIpKuyL0L0ySdUZFtfUm35V6kHxS2nV14f0qt\n50mpl+uiXO5kSbtKulXS44UepkU9LbmuP5N0Z+5d+kSDNvxP4KqImJzbcRbwPeDLjdq4gRtY/Ds9\nCripcExL9AblNh4kqZekayVNlzRF0qerCpc0W9Kw3Dvze0kfkDQut8VeheM/pbhN/ruppLtz785k\n5d47SZ/JnydJurDQ5jvl9ztJGtvBcdtScCBiZma22HOk+SLd8So7ABgfETOA2ZL6N1H/kZKmAFOB\nb3bu0Nu5ATg63zl/EygOt9oBeLCU/8GcXikiXgDGADMlXS3pqFKWH2rx0KxLGhQ1HHhc0o2Sjs9B\nG8ClwPcjYiBwPPCTwjYbAntExLmN6piP4SPAB0hryEzLC1Y+DxzWYLuvRcRuQH/g3yVtXifPnqSn\nr+0I7CFpUAd1AZiby70+vz4JDATOrsi/DXAgcDDw3QbldnQOq9q4ykSgbw7WBgPNDEMbAGwVETtE\nxM4sGeyWbQjcFBF9gXcCpwAfBE4GvtbBfo4BxkZEf9IipY9L2hk4HfhgTv9eE/W1ZcyBiJmZrdJ6\n8BombcCN+f0N+XNHBucLur7A+ZLWJvVu1NPR02ruAf4f8Amau6hUB2UGQEScQLqgnwZcIOk7hTyn\nRsSA/Dq9sqCI/wb2IvXCnAaMyF8dBFwp6SFS221S2GxkNPeEnjER8VpE/IMUgI3O6VNoPDzpGEkT\nSRfz2wHb18lzd0TMioj5pDbdt4n6FPc/PiLmRMQc4BVJ9da2uSUi3oqIJ4D1GpRb73wtSmvQxo38\nETgWeAdpqFZH/gZsKukySYdExEsN8s6NiL/k91OBOyJiIR2fF4AHgGMlfRvoGxGvAfsD1+e2rAXJ\n1mKeI2JmZpZ0xxomdSnNpTgSOETS+UAvYB5wVjPbR8STkp4hBSTTab8Y4wYs2cNRr4yQdCfpbnNf\n0l3lmunALsDYQtoupItWgLckrZYvFGv7e7hQ9kRgoqQxpAvcTvfe5J6iGZKuA2YWvhqYL/TLXiu8\nn8/im7Hl4VbzCu8XRkTt80LSeWhH0takO/N7RcRLeRhSvWFc5Qv/qJNeVZ+F5bpV1GdenbR6ppF6\nB0YX0nZhyfNU1cZVbgTuAL5dSi+2N+RjzPOLdiYFpqfnYKRqWFv52Oudl6r93ClpX+Bw4NeSGv07\navTbsGXMPSJmZmYs6hnptlepOgcD4yJii4joExGbA09LGtjMsUjaCNgamBURc4E5kvbO320G7Exh\n8nkDl5HmgjxfSr8IOE/SernM/UjDZGqT4+8l9aSQ8xwC/G+ej/HBQjn9gFnNHFPp+A6TFi0AUyzj\nL8AXCvn6lbfN/k4aFgTwsc7uv461ScPrXs7te1BFvn0kbSZp9bzfcTn9JUlb5gC08mlmy9jlwGcK\n8yE2A74OXJw/V7VxpYh4CDgH+HXpq1nADpJWl/RuoPZb3AhYLSJuAIay+Jx01aLzKulg8s0FSVsC\n/4qIn5Im6PcjBUxtktbNeTYol8Gy+W1YA+4RMTMz63naKEz2zUaRLu5vAQ6TVBz6UhviM1LSG8Db\ngbMj4tmcfgIwXNI6pDu+J+UApaG8any7leMjYlS+cL1fUpAmNR+Vh7wAfBG4QtJX8+fvRsTf81Cx\nr0u6AniD9HSu4mNhfyipNq/hxYjYv6JqJwCXSHodeB0YktP/C/iJpM/mNhgNTK6z/c+A30v6GHBz\n41boWERMkvQw8Aip52BcRdb78r63AkZHxF9z+jdIF8bP5jKWu4j4p6QhwC+VHpawAPhmRNTmjVS1\ncUflXgrpIQCFtFmS/kDqSZtBmk8C8F5ghKTVSL/L05bysH4HDMnzpMaQ5vUADAK+IuktYA5wTEQ8\nK2kYcLek+aRhZV8lBWLXS/oi6ZyQj+cIUm/bt6reL2XdV0le0NDMzMzMzFrOQ7PMzMzMzKzlPDTL\nzMxsFSXpPtpPyD2wzpyQbiHpbODjpeRhEXFVd9QHQNIo0tCqouMiYkp31KcjeTL4NaXkJyOi/Ojk\nZspa7udjRWtfWzoemmVmZmZmZi3noVlmZmZmZtZyDkTMzMzMzKzlHIiYmZmZmVnLORAxMzMzM7OW\ncyBiZmZmZmYt50DEzMzMzMxazoGImZmZmZm1nAMRMzMzMzNrOQciZmZmZmbWcg5EzMzMzMys5RyI\nmJmZmZlZyzkQMTMzMzOzlnMgYmZmZmZmLedAxMzMzMzMWs6BiJmZmZmZtZwDETMzMzMzazkHImZm\nZmZm1nIORMzMzMzMrOUciJiZmZmZWcs5EDEzMzMzs5ZzIGJmZmZmZi3nQMTMzMzMzFrOgYiZmZmZ\nmbWcAxEzMzMzM2s5ByJmZmZmZtZyDkTMzMzMzKzlHIiYmZmZmVnLORAxMzMzM7OWcyBiZmZmZmYt\n50DEzMzMzMxazoGImZmZmZm1nAMRMzMzMzNrOQciZmZmZmbWcg5EzMzMzMys5RyImJmZmZlZyzkQ\nMTMzMzOzlnMgYmZmZmZmLedAxMzMzMzMWs6BiJmZmZmZtZwDETMzMzMzazkHImZmZmZm1nIORMzM\nzMzMrOUciJiZmZmZWcs5EDEzMzMzs5ZzIGJmZmZmZi3nQMTMzMzMzFrOgYiZmZmZmbWcAxEzMzMz\nM2s5ByJmZmZmZtZyDkTMzMzMzKzlHIiYmZmZmVnL/X8sJCqhXNjMSgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7c56630>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Albumin phyrogenetic tree\n",
    "plot_tree('tree.newick.txt',8,3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    ">N1\n",
    "MKWVTLISLLFLFSSATSRNLQRFRRDAEAHKSEIAHRYNDLGEEHFKGLVLITFAQYLQKCPYEELAKLVKEVTDLAQACVADESAADCSKPLHTIFLDKICAVPKLRDTYGAMADCCAKADPERNECFLSHKDSQPDLVPPYQRPEPDVLCQAYQDNKESFLGHYIYEVARRHPFLYAPAILSFAQKFKAVLTECCEEADKGACLTTKLTALREKALIVSVKQRLSCGILQKFGDRVFQAWQLVRLSQKYPKAPFAEVSKLVTDLTKVHKECCHGDMLECMDDRADLTKHMCEHQDTISSKLKECCEKPIVERSHCIVELENDEMPADLPSLVEKFVEDKEVCKSFEEAKDVFLAEFLYEYSRRHPEFSVQLLLRIAKGYESTLEKCCETDNPHECYANAQDELNQLIKEPQDLVKQNCELLQKLGEYNFQNALLIRYTKKMPQVSTPTLVEISKSMTKVGSKCCKLPEAQRMPCAEGYLSVVINELCVLQETTPINENVTKCCSQSYANRRPCFTALGVDETYVPPEFNADTFTFHEDLCTLPEEERKIKKQTLLVNLVKHKPHVTEEQLKTIAGEFTAMVDKCCAAEDKEACFAEEGPKLIEQSKATLGLGA\n",
    ">N2\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEIAHRFNDLGEEHFKGLVLITFSQYLQKCPYEEHAKLVKEVTDLAKACVADESAANCDKSLHTIFGDKICAVPSLRDTYGDMADCCEKQEPERNECFLQHKDDKPDLVPPFARPEPDVLCKAFHDNEEAFLGHYLYEVARRHPYFYAPELLYYAQKYKAVLTECCEAADKGACLTPKLDALREKALISSAKQRLRCASLQKFGDRAFKAWALVRLSQKFPKADFAEISKLVTDLTKVHKECCHGDLLECADDRADLAKYMCEHQDTISSKLKECCDKPILEKSHCIAELENDEMPADLPALAEEFVEDKDVCKNYEEAKDVFLGKFLYEYSRRHPDYSVSLLLRLAKAYEATLEKCCATDDPHACYAKVLDEFKPLVEEPQNLVKQNCELFEKLGEYNFQNALLVRYTKKVPQVSTPTLVEISRSLGKVGSKCCKHPEAERMPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCSESLVNRRPCFSALGVDETYVPKEFNAETFTFHADICTLPETERKIKKQTALVELVKHKPHATEEQLKTVVGEFTALVDKCCAAEDKEACFAEEGPKLVESSKATLGLGA\n",
    ">N3\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEIAHRFNDLGEEHFKGLVLIAFSQYLQQCPFEEHVKLVNEVTEFAKTCVADESAANCDKSLHTLFGDKLCTVASLRETYGEMADCCEKQEPERNECFLQHKDDNPDLVPPLVRPEPDAMCTAFHDNEETFLGKYLYEVARRHPYFYAPELLYYAEKYKAVFTECCQAADKAACLTPKLDALREKVLASSAKQRLKCASLQKFGERAFKAWAVARLSQKFPKADFAEISKLVTDLTKVHKECCHGDLLECADDRADLAKYMCENQDSISSKLKECCDKPLLEKSHCIAEVENDEMPADLPALAADFVEDKDVCKNYQEAKDVFLGTFLYEYSRRHPDYSVSLLLRLAKAYEATLEKCCATDDPHACYAKVFDEFKPLVEEPQNLVKQNCELFEKLGEYGFQNALLVRYTKKVPQVSTPTLVEVSRSLGKVGSKCCKHPEAERMPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADICTLPETEKQIKKQTALVELVKHKPKATEEQLKTVMGDFAAFVDKCCAAEDKEACFAEEGPKLVASSQAALALGA\n",
    ">N4\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEIAHRFNDLGEEHFKGLVLIAFSQYLQQCPFEEHVKLVNEVTEFAKTCVADESAENCDKSLHTLFGDKLCTVATLRETYGEMADCCEKQEPERNECFLQHKDDNPNLVPPLVRPEPDAMCTAFHDNEETFLGKYLYEVARRHPYFYAPELLYYAEKYKAVFTECCQAADKAACLTPKLDALREKVLASSAKQRLKCASLQKFGERAFKAWAVARLSQKFPKADFAEVSKLVTDLTKVHKECCHGDLLECADDRADLAKYMCENQDSISSKLKECCDKPLLEKSHCIAEVENDEMPADLPALAADFVEDKDVCKNYAEAKDVFLGTFLYEYSRRHPDYSVSLLLRLAKAYEATLEKCCATADPHACYAKVFDEFKPLVEEPQNLVKQNCELFEKLGEYGFQNALLVRYTKKVPQVSTPTLVEVSRSLGKVGSKCCKHPEAERMPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADICTLPEKEKQIKKQTALVELVKHKPKATEEQLKTVMGDFAAFVDKCCKAEDKEACFAEEGPKLVASSQAALALGA\n",
    ">N5\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEIAHRFNDLGEKHFKGLVLIAFSQYLQQCPFEEHVKLVNEVTEFAKTCVADESAENCDKSLHTLFGDKLCTVATLRETYGEMADCCEKQEPERNECFLQHKDDNPNLVPPLVRPEPDAMCTAFQENPETFLGKYLYEVARRHPYFYAPELLYYAEKYKAVFTECCQAADKAACLTPKLDALKEKVLVSSAKQRLKCSSLQKFGERAFKAWAVARLSQKFPKADFAEVSKLVTDLTKVHKECCHGDLLECADDRADLAKYMCENQDSISSKLKACCDKPLLQKSHCIAEVENDDMPADLPALAADFVEDKDVCKNYAEAKDVFLGTFLYEYSRRHPDYSVSLLLRLAKTYEATLEKCCAEADPHACYATVFDEFKPLVEEPQNLVKQNCELFEKLGEYGFQNALLVRYTKKAPQVSTPTLVEVSRSLGKVGSKCCKLPEAERLPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCTESLVERRPCFSALEVDETYVPKEFKAETFTFHADICTLPEKEKQIKKQTALAELVKHKPKATEEQLKTVMGDFAAFVDKCCKAEDKEACFAEEGPKLVASSQAALALGA\n",
    ">N6\n",
    "MKWVTFLLLLFVSGSAFSRGVQRFRRDAEAHKSEIAHRYKDLGEKHFKGLVLIAFSQYLQKCPYEEHVKLVQEVTDFAKTCVADESAENCDKSLHTLFGDKLCAIPNLRENYGEMADCCAKQEPERNECFLQHKDDNPNLVPPFQRPEPDAMCTAFQENPETFMGHYLHEVARRHPYFYAPELLYYAEKYNAVLTECCAAADKAACLTPKLDALKEKALVSAVRQRLKCSSMQKFGERAFKAWAVARMSQTFPNADFAEITKLATDLTKVNKECCHGDLLECADDRAELAKYMCENQASISSKLQACCDKPLLQKSHCLAEVEHDDMPADLPALAADFVEDKDVCKNYAEAKDVFLGTFLYEYSRRHPDYSVSLLLRLAKKYEATLEKCCAEADPHACYGTVFDEFKPLVEEPQNLVKTNCELYEKLGEYGFQNAVLVRYTKKAPQVSTPTLVEAARSLGRVGTKCCTLPEAQRLPCVEDYLSAILNRVCVLHEKTPVSEKVTKCCSGSLVERRPCFSALTVDETYVPKEFKAETFTFHADICTLPEKEKQIKKQTALAELVKHKPKATEEQLKTVMGDFAEFVDKCCKAEDKEACFSTEGPKLVARSQEALALGA\n",
    ">N7\n",
    "MKWVTFLLLLFVSGSAFSRGVQRFRREAEAHKSEIAHRYKDLGEQHFKGLVLIAFSQYLQKCPYEEHVKLVQEVTDFAKTCVADESAENCDKSLHTLFGDKLCAIPNLRENYGELADCCAKQEPERNECFLQHKDDNPNLVPPFQRPEAEAMCTSFQENPTTFMGHYLHEVARRHPYFYAPELLYYAEKYNEVLTQCCAEADKAACLTPKLDAVKEKALVSAVRQRMKCSSMQKFGERAFKAWAVARMSQTFPNADFAEITKLATDLTKVNKECCHGDLLECADDRAELAKYMCENQATISSKLQACCDKPLLQKSHCLAEVEHDNMPADLPAIAADFVEDKEVCKNYAEAKDVFLGTFLYEYSRRHPDYSVSLLLRLAKKYEATLEKCCAEADPPACYGTVLAEFQPLVEEPKNLVKTNCELYEKLGEYGFQNAVLVRYTQKAPQVSTPTLVEAARNLGRVGTKCCTLPEAQRLPCVEDYLSAILNRVCVLHEKTPVSEKVTKCCSGSLVERRPCFSALTVDETYVPKEFKAETFTFHSDICTLPEKEKQIKKQTALAELVKHKPKATEEQLKTVMGDFAQFVDKCCKAADKDTCFSTEGPNLVARSKEALALGA\n",
    ">N8\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEVAHRFKDLGEEHFKGLVLIAFSQYLQQCPFEEHVKLVNEVTEFAKTCVADESAENCDKSLHTLFGDKLCTVATLRETYGEMADCCAKQEPERNECFLQHKDDNPNLVPPLVRPEVDVMCTAFHDNEETFLKKYLYEVARRHPYFYAPELLFFAARYKAAFTECCQAADKAACLLPKLDELRDEGKASSAKQRLKCASLQKFGERAFKAWAVARLSQKFPKAEFAEVSKLVTDLTKVHTECCHGDLLECADDRADLAKYMCENQDSISSKLKECCDKPLLEKSHCIAEVENDEMPADLPSLAADFVESKDVCKNYAEAKDVFLGMFLYEYARRHPDYSVVLLLRLAKAYEATLEKCCAAADPHECYAKVFDEFKPLVEEPQNLVKQNCELFEQLGEYKFQNALLVRYTKKVPQVSTPTLVEVSRNLGKVGSKCCKHPEAKRMPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADICTLSEKEKQIKKQTALVELVKHKPKATKEQLKTVMDDFAAFVEKCCKADDKEACFAEEGPKLVAASQAALALGA\n",
    ">N9\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEVAHRFKDLGEENFKALVLIAFAQYLQQCPFEDHVKLVNEVTEFAKTCVADESAENCDKSLHTLFGDKLCTVATLRETYGEMADCCAKQEPERNECFLQHKDDNPNLVPRLVRPEVDVMCTAFHDNEETFLKKYLYEIARRHPYFYAPELLFFAKRYKAAFTECCQAADKAACLLPKLDELRDEGKASSAKQRLKCASLQKFGERAFKAWAVARLSQRFPKAEFAEVSKLVTDLTKVHTECCHGDLLECADDRADLAKYICENQDSISSKLKECCEKPLLEKSHCIAEVENDEMPADLPSLAADFVESKDVCKNYAEAKDVFLGMFLYEYARRHPDYSVVLLLRLAKTYETTLEKCCAAADPHECYAKVFDEFKPLVEEPQNLIKQNCELFEQLGEYKFQNALLVRYTKKVPQVSTPTLVEVSRNLGKVGSKCCKHPEAKRMPCAEDYLSVVLNQLCVLHEKTPVSERVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADICTLSEKERQIKKQTALVELVKHKPKATKEQLKTVMDDFAAFVEKCCKADDKETCFAEEGKKLVAASQAALGLGA\n",
    ">N10\n",
    "MKWVTFISLLFLFSSAYSRGVQRFRRDAEAHKSEIAHRFNDLGEEHFKGLVLIAFSQYLQQCPFEEHVKLVNEVTEFAKTCVADESAANCDKSLHTLFGDKLCTVASLRETYGDMADCCEKQEPERNECFLQHKDDNPDLVPPLVRPEPDAMCTAFHDNEQRFLGKYLYEIARRHPYFYAPELLYYAEKYKGVFTECCQAADKAACLTPKIDALREKVLASSAKQRLKCASLQKFGERAFKAWSVARLSQKFPKAEFAEISKLVTDLTKVHKECCHGDLLECADDRADLAKYMCENQDSISSKLKECCDKPLLEKSHCIAEVEKDEMPADLPPLAADFVEDKDVCKNYQEAKDVFLGTFLYEYSRRHPEYSVSLLLRLAKEYEATLEKCCATDDPHACYAKVFDEFKPLVEEPQNLVKQNCELFEKLGEYGFQNALLVRYTKKVPQVSTPTLVEVSRSLGKVGSKCCKHPEAERMPCAEDYLSVVLNRLCVLHEKTPVSEKVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADICTLPETEKQIKKQTALVELLKHKPKATEEQLKTVMGDFAAFVDKCCAAEDKEACFAEEGPKLVASSQAALALGA\n",
    ">N11\n",
    "MKWVTFISLLFLFSSAYSRGVQRVRREAEAHKSEIAHRFNDLGEEHFRGLVLVAFSQYLQQCPFEDHVKLVNEVTEFAKACVADESAANCDKSLHTLFGDKLCTVASLRDKYGDMADCCEKQEPERNECFLQHKDDNPGFVPPLVTPEPDAMCTAFHDNEQRFLGKYLYEIARRHPYFYAPELLYYAEKYKGVFTECCQAADKAACLTPKIDALREKVLASSAKERLKCASLQKFGERAFKAWSVARLSQKFPKAEFAEISKLVTDLTKVHKECCHGDLLECADDRADLAKYMCENQDSISTKLKECCDKPVLEKSHCIAEVERDELPADLPPLAADFVEDKEVCKNYQEAKDVFLGTFLYEYSRRHPEYSVSLLLRLAKEYEATLEKCCATDDPPACYAKVFDEFKPLVEEPQNLVKTNCELFEKLGEYGFQNALLVRYTKKVPQVSTPTLVEVSRSLGKVGSKCCKHPEAERMSCAEDYLSVVLNRLCVLHEKTPVSERVTKCCTESLVNRRPCFSALEVDETYVPKEFNAETFTFHADLCTLPEAEKQIKKQTALVELLKHKPKATEEQLKTVMGDFGAFVDKCCAAEDKEACFAEEGPKLVAAAQAALALGA"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [conda root]",
   "language": "python",
   "name": "conda-root-py"
  },
  "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
}