{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2018.03.4\n" ] } ], "source": [ "%matplotlib inline\n", "from rdkit import rdBase\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from rdkit import Chem\n", "from rdkit.Chem import AllChem\n", "from rdkit.Chem import DataStructs\n", "from rdkit.Chem import Draw\n", "from rdkit.Chem.Draw import IPythonConsole\n", "from mhfp.encoder import MHFPEncoder\n", "mhfp_encoder = MHFPEncoder()\n", "print(rdBase.rdkitVersion)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "mols = [mol for mol in Chem.SDMolSupplier('cdk2.sdf') if mol != None]\n", "for mol in mols:\n", " AllChem.Compute2DCoords(mol)\n", "nmols = len(mols)\n", "#Calc morgan fp\n", "mg2fps = [AllChem.GetMorganFingerprintAsBitVect(mol, 3) for mol in mols]\n", "#Calc min hash fp\n", "mhfps = [mhfp_encoder.encode_mol(mol) for mol in mols]" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "tanimoto_sim = []\n", "mg_sim_mols = []\n", "for i in range(nmols):\n", " for j in range(i):\n", " tc = DataStructs.TanimotoSimilarity(mg2fps[i], mg2fps[j])\n", " if tc >= 0.5:\n", " mg_sim_mols.append([mols[i], mols[j]])\n", " tanimoto_sim.append(tc)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "mhfps_sim = []\n", "mhfps_sim_mols = []\n", "for i in range(nmols):\n", " for j in range(i):\n", " jaccard = 1. - MHFPEncoder.distance(mhfps[i], mhfps[j])\n", " if jaccard >= 0.5:\n", " mhfps_sim_mols .append([mols[i], mols[j]]) \n", " mhfps_sim.append(jaccard)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "1.033917242502858 -0.031604772419224866\n" ] } ], "source": [ "a, b = np.polyfit(tanimoto_sim, mhfps_sim, 1)\n", "y2 = np.int64(a) * tanimoto_sim + np.int64(b)\n", "print(a, b)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'mhfp sim')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzt3Xt4VOW1+PHvShggQSBB4oVIABVB1AI1Rf2htlQUWhXQKoLa6mkFUal3jlQ4XmhBlBbktHhB1FZFUcEiKG3UYit4RAgGVEAUUAIBBIEoQoBc1u+PuTCZzD2zMzPJ+jxPnmbP7JlZmeJe+72tV1QVY4wxBiAj2QEYY4xJHZYUjDHG+FhSMMYY42NJwRhjjI8lBWOMMT6WFIwxxvg4mhREZKCIrBeRDSIyNsjzBSLyroiUiMjHIvJzJ+MxxhgTnji1TkFEMoHPgQuBrcAKYLiqrvU7ZyZQoqqPi0gPYJGqdnYkIGOMMRE52VLoA2xQ1U2qehiYAwwOOEeBNp7f2wLbHIzHGGNMBM0cfO98YIvf8VbgrIBzHgDeEpHfAq2A/pHetH379tq5c+cEhWiMMU3DypUrv1HVvEjnOZkUJMhjgX1Vw4G/quqfROQc4HkROV1Va2q9kchIYCRAQUEBxcXFjgRsjDGNlYhsjuY8J7uPtgId/Y5PoG730G+AVwBU9QOgJdA+8I1UdaaqFqpqYV5exERnjDEmTk4mhRVAVxHpIiLNgWHAgoBzSoELAETkVNxJYZeDMRljjAnDsaSgqlXAaKAIWAe8oqprRGSCiAzynHYXMEJEVgMvAderlW01xpikcXJMAVVdBCwKeOw+v9/XAn2djMEYY0z0bEWzMcYYH0sKxhhjfCwpGGOM8bGkYIwxxseSgjHGGB9HZx8ZY0xTN7+kjClF69lWXkGHnCzGDOjGkN75yQ4rJEsKxhjjkPklZfzutU+oqKwGoKy8gt+99glAyiYG6z4yxhiHTCla70sIXhWV1UwpWp+kiCKzloIxpklIRjfOtvKKmB5PBdZSMMY0et5unLLyCpQj3TjzS8oc/dwOOVkxPe5VXl5OQUEB48ePdyKssCwpGGMavWR144wZ0I0sV2atx7JcmYwZ0C3ka+666y5yc3PZsmULL730kqPxBWNJwRjT6CWrG2dI73weuvwM8nOyECA/J4uHLj8jaLfVihUrEBGmTp0KwL333svGjRsdjS8YG1MwxjR6HXKyKAuSACJ14yTCkN75YccuDh06RPfu3fnqq68AaNmyJV9//TVt2rQJ+RonWUvBGNPoxdON0xCmT59Oy5YtfQnh7bffpqKiImkJAaylYIxpArx36qmyiGzjxo2cfPLJvuPhw4cze/ZsRILtYtywLCkYY5qESN04DaGmpoYBAwbwzjvv+B4rKyujQ4cOSYyqNus+MsaYBjBv3jwyMzN9CeG5555DVVMqIYDDLQURGQhMBzKBWao6OeD5aUA/z2E2cIyq5jgZkzHGBOPU4rZvvvmGvLw833GfPn14//33adYsNTtqHItKRDKBGcCFwFZghYgs8GzBCYCq3uF3/m+B3k7FY4wxoThVo+jGG29k5syZvuO1a9dy6qmnRhVPssY/nOw+6gNsUNVNqnoYmAMMDnP+cKDhV2oYY5q8RC9uW7p0KSLiSwgTJ05EVaNOCMlYfe3lZPslH9jid7wVOCvYiSLSCegCLA7x/EhgJEBBQUFiozTGNHmJWtx24MABOnXqxDfffANA+/bt2bx5M9nZ2VG/R7gE1RCtBSdbCsHmVmmIc4cBc1W1OtiTqjpTVQtVtdC/b84YYxIh3hpF/h566CFatWrlSwjvvfceu3btiikhQPKL6DmZFLYCHf2OTwC2hTh3GNZ1ZIxJkvosblu3bh0iwr333gvAiBEjUFXOO++8uGJJRIKqDyeTwgqgq4h0EZHmuC/8CwJPEpFuQC7wgYOxGGNMSLHUKPKqrq7m7LPPpkePHr7Hdu7cWWtgOR7JXn3t2JiCqlaJyGigCPeU1GdUdY2ITACKVdWbIIYDc1Q1VNeSMcY4LpbFbS+88AK//OUvfcevvPIKV155ZcLigOStvpZ0uxYXFhZqcXFxssMwxjRB27dvr7XY7IILLuCtt94iIyP11wGLyEpVLYx0Xur/JcYYk2SqyjXXXFMrIWzYsIF33nknLRJCLBrXX2OMMQnmvfC/+OKLAEybNg1V5aSTTkpyZM5IzXXWxhiTZPv27eOYY47h4MGDAHTu3Jl169bRsmXLJEfmLGspGGNMgPHjx9OmTRtfQvjwww/58ssvG31CAGspGGOMz+rVq+nVq5fv+Pbbb2fatGlAcusRNSRLCsaYJq+yspJevXqxdq2vXid79uwhNzcXcK5gXiqy7iNjTFqbX1JG38mL6TL2TfpOXhxz4binnnqK5s2b+xLCwoULUVVfQoDEF8xLZdZSMMakrfrcwZeWltKpUyff8aBBg5g/f37QLTGTXY+oIVlLwRiTtuK5g1dVBg8eXCshbN68mddffz3kHsn1qUdU35ZMQ7OkYIxJW7Hewb/xxhtkZGSwYIG7ys6TTz6JqkYsyR9vPaJk740QD+s+MsakrQ45WZQFSQCBd/B79+6lXbt2vuPTTz+djz76CJfLFdXnxFuPKNl7I8TDkoIxJm2NGdCt1pgC1L2Dv+OOO3j00Ud9x6tWraJnz54xf1YsBfO80nEswrqPjDFpK1zJ6+XLlyMivoQwbtw4VDWuhBCvZO+NEA9rKRhj0lrgHfzBgwfp1KkTpaWlAGRnZ7Njxw5at27d4LFF05JJNdZSMMY0GtOnTycrK8uXEN5++23279+flIQA8W3ek2zWUjDGpL2NGzdy8skn+46vvvpqXnjhhZBTTBtSPGMRyeRoS0FEBorIehHZICJjQ5wzVETWisgaEXnRyXiMMY1LTU0NF1xwQa2EsG3bNmbPnp0SCSEdOZYURCQTmAH8DOgBDBeRHgHndAV+B/RV1dOA252KxxjTuLz66qtkZmayePFiAJ5//nlUleOPPz7JkaU3J7uP+gAbVHUTgIjMAQYDa/3OGQHMUNW9AKq608F4jDGNwK5duzjmmGN8x2effTZLly4lMzMzzKtMtJzsPsoHtvgdb/U85u8U4BQReV9ElonIQAfjMcakuREjRtRKCGvXruWDDz6whJBATiaFYB16GnDcDOgK/AQYDswSkZw6byQyUkSKRaR4165dCQ/UGJPali5diogwa9YsACZOnIiqcuqppyY5ssbHye6jrUBHv+MTgG1BzlmmqpXAlyKyHneSWOF/kqrOBGYCFBYWBiYWY0yKSdSGNAcOHKCgoIDdu3cDkJeXx1dffUV2dnaiQzYeTrYUVgBdRaSLiDQHhgELAs6ZD/QDEJH2uLuTNjkYkzHGYYkqAjdp0iRatWrlSwhLlixh586dlhAc5lhSUNUqYDRQBKwDXlHVNSIyQUQGeU4rAnaLyFrgXWCMqu52KiZjjPPquyHNunXrEBHGjRsHuMcRVJVzzz034bGauhxdvKaqi4BFAY/d5/e7And6fowxjUC8ReCqqqro27cvy5cv9z22c+dO8vLyEhqfCc/KXBhjEiqeInDPP/88LpfLlxBeffVVVNUSQhJYmQtjTELFUgRu+/btdOjQwXfcv39/ioqKyMiw+9VksaRgjEmoaDakUVWuvvpq5syZ43tsw4YNnHTSSSHfN1Ezmkx4lhSMMQkXqgjc/JIyxv1lNmufucf32KOPPsptt90W9v28M5q8rQ/vjCbvZ5nEsaRgTJzszjU2Ly1dzzX9foBWHQYgs+2xnHTTTDqdf2bE16bjtpbpypKCMXGwO9fYjB8/nokTJ/qOj/vVVFocfwqHFO56ZTUQ/ntLx20t05UlBWPiYHeu0Vm1ahW9e/f2HbcuHEy7C0bUOqdaNWJC7ZCTRVmQBJDK21qmKxviNyYOducaXmVlJaeddlqthNDxtjl1EoJXpMVtYwZ0I8tVu+hdqm9rma6spWBMHNLpzrWhxz6efPJJRo0a5Tt+4403mPRJ8O/LX7iEGs2MJpMYlhSMiUO6bMjekGMfpaWldOrUyXc8ZMgQXnvtNUSE0UvejPj6SAk13ba1TFeWFIyJINyddqrfuTbE2IeqMnjwYBYuXOh7rLS0lI4djxRJDtWy8krFhNpUWVIwJoxId9qplgQCOT32sXDhQgYNGuQ7njlzJiNG1B03CNayEtwbrOTXI6HatODEs6RgTBjpPsvIqbGPvXv30q5dO9/x6aefzkcffcSbn+6k7+TFdS7STrSsbFqwMywpGBNGus8ycmLs44477uDRRx/1Ha9atYqePXtG1apK5MU63RN2qrIpqcaEEU/Fz1QypHc+D11+Bvk5WQjurpqHLj8jrovm8uXLERFfQhg/fjyqSs+ePYHo91GYX1JG38mL6TL2TfpOXhzz5jte6Z6wU5W1FIwJI11mGYVT3zv0gwcP0q1bN0pLSwHIzs5mx44dtG7dutZ50VykE9nlk07TgtOJtRSMCSORd9rpaNq0aWRlZfkSwjvvvMP+/fvrJASIrlVV313Z/NmCNmc42lIQkYHAdCATmKWqkwOevx6YAnjbj39R1VlOxmRMrBLVF57KM2UCY7v21ObcPOQ83/PXXnstzz33HCIS8j2iaVUlsssnXaYFpxvHkoKIZAIzgAuBrcAKEVmgqmsDTn1ZVUc7FYcxqSCR3SbzS8p4cOEa9h6oBCAny8UDg06L+2LoH5tqDSsfv4P/K/3Y9/y2bds4/vjjI76HtxWQKUK1atCpponu8kmHacHpxsnuoz7ABlXdpKqHgTnAYAc/z5iUlahuk/klZYyZu9qXEADKKyoZ8+rquAdsvbHt/2wppY8M4pAnIXQd+jtUNaqE8LvXPvFd7KtVfS2EwAu2dfmkPieTQj6wxe94q+exQL8QkY9FZK6IdAzyPCIyUkSKRaR4165dTsRqjKMS1W0ypWg9ldVa5/HKGo2rXx5gy7YdbH74Er553d2726JDdwrGvE5ll75RxxRtwmvqYzTpwMkxhWCdj4H/mhcCL6nqIREZBfwN+GmdF6nOBGYCFBYW1v0vwpgUl6huk3BJJJ5++REjRrBl1pFhvA43PI7r6I4xxRaqfEVZeQXzS8rqXPCtyye1OdlS2Ar43/mfAGzzP0FVd6vqIc/hU0DkLZiMSUOJ6jYJd6FWiHre/5IlSxARZnkSQl6//6LTPW/4EkK0sc0vKQt69+f1u9c+ibtbK9EStT6isXMyKawAuopIFxFpDgwDFvifICL+nZWDgHUOxmNM0iSq22TMgG64MkNfhr0D2KEuePv37+foo4/m/PPPB+CYY45h//79zPzT7+OKbUrR+jrNf3/xTjdNNP9xDyXy99SUOdZ9pKpVIjIaKMI9JfUZVV0jIhOAYlVdANwqIoOAKmAPcL1T8RiTbInoNvG+3n/2UaBQpR4mTpzI+PHjfcdLly6lb9++nvfNjiu2aLqsUmGFsZXEiJ6j6xRUdRGwKOCx+/x+/x3wOydjMKax8U8uXca+GfRO3f9CvHbtWk477TTf8Y033sgTTzyRkFgilcT2npNsVhIjerai2Zg0lpPtCvl4VVUVP/rRj2olhF27diUsIUDwsRJ/qTLdNN1rWDUkSwrGpDEN0aG/u+RtXC4XxcXFAMydOxdVpX379gn9/MCxktxsFzlZrpSbbmrrI6JnBfGMSWPfVtQeV6jat5uyx67zHed0LaTtZfcx9YtWZAaZHpoI6TDF1EpiRM+SgjFpzNunr6p8s+ARDny2xPfcibc8Q/VRxwC2AQ2kR/JKBdZ9ZEwaGzOgG5n7vqb0kUt9CeGYi0bR68EiX0LwSpXpoSa1RUwKInKJiJSIyB4R+U5E9onIdw0RnDEmtEOHDvHJG8+w9embAWiWcxxn//4fPDl5POUhpqvabBsTSTTdR48ClwOfqIYa1jIm/aVyaetA//rXv7j55pv5/PPPGTp0KFOnTiU//0isU4rW2wY0Ji7RdB9tAT61hGAas3RZ8bpjxw6uueYa+vfvT3V1NUVFRbz88su1EgKkx2wbKzuRmqJpKfw3sEhE/gN46xShqlMdi8qYBpbqK16rq6t54oknGDduHBUVFdx///2MHTuWli1bBj0/1WfbJHJ/CZNY0SSFicD3QEugubPhGJMcqbzitbi4mFGjRrFy5Ur69+/PjBkzOOWUUyK+zunZNvXpbkv1JNyURZMU2qnqRY5HYkwSJaK0daLHJMrLyxk/fjyPPfYYxx13HHPmzGHo0KFht8RsKPW900/lJNzURTOm8I6IWFIwSedkH3R9++ATOSahqrz44ot0796dxx9/nN/+9resW7eOq666KiUSAtR/JzkrO5G6okkKtwD/FJEKm5JqksXpgeD6lrYOdZG865XYtslcv349F154Iddccw0FBQWsWLGC6dOn07Zt21j+HMfV904/HQbCm6qI3Ueq2rohAjEmnIbog65PH3yoSqHVqlF1q1RUVDBp0iQeeeQRsrKyeOyxxxg5ciSZmaGLzSVTfbvbUn0gvCkLmRREpLuqfiYiPwz2vKp+5FxYxtSWyn3Q3t3HQs3ZjpS8/vGPfzB69Gg2bdrEtddeyx//+EeOPfZYx+JNhDEDutUaU4DY7/St7ERqCtdSuBMYCfwpyHNKkL2UjXFKovY4dkKk3ccgePLaunUrt99+O/PmzaN79+4sXryYfv36ORNkgtmdfuMVMimo6kjP/8b9r1REBgLTce+8NktVJ4c47wrgVeBHqloc7+eZxisRd6ZOiaa14p+8qqqq+POf/8x9991HVVUVEydO5O6776Z58/Sa8W13+o1TNLWPrhSR1p7fx4vIayLSO4rXZQIzgJ8BPYDhItIjyHmtgVuBD2MN3jQdidrj2AmRWiv+yeuDDz6gsLCQO++8k/PPP5+1a9dy7733pl1CMI1XNOsU/kdVXxWRc4EBwB+BJ4CzIryuD7BBVTcBiMgcYDCwNuC83wOPAHfHErhpehr6zjTadQfBWjHeMYZ8z+vO75TFyJEjeeqppzjhhBOYN28el112WcpMMTXGK5qk4P2XfjHwuKq+LiIPRPG6fNx1k7y2EpBIPC2Ojqr6hohYUjApI5bFWeH611WVv/3tb3S7aAx79+5l8C9vZMdJl3LXcvjT5+9aP7xJOdEkhTIReRLoDzwsIi2Ibn1DsFsg33iciGQA04DrI76RyEjcg94UFBRE8dHG1E+sU2CDtWLWrFnDTTfdxJIlS+jes5Cjr/w9q9p0BM8QhNX7Makomov7UKAIGKiq5UA7YEwUr9sKdPQ7PgHY5nfcGjgd+LeIfAWcDSwQkcLAN1LVmapaqKqFeXl5UXy0MfGbX1IWct1BNIPK+/fv55577qFXr16sWbOGW+6bgl7yIAfbdKxzrm18Y1JNxKSgqgdU9TVV/cJzvF1V34rivVcAXUWki4g0B4YBC/ze91tVba+qnVW1M7AMGGSzj0wyebuNQok0qPz666/To0cPHnnkEX71q1+xfv16SrJ+yMGq0JNWU2GthTFejm3HqapVwGjcrYx1wCuqukZEJojIIKc+15j6CNZt5BVuCuzmzZsZNGgQQ4YMoU2bNixZsoSnn36a9u3bR7zop8JaC2O8ohlTiJuqLgIWBTx2X4hzf+JkLMZEI9wFPNgU2MOHDzNt2jQefPBBMjIymDJlCrfddhsul8t3TqiFd+D8Wot02k3OpIaoWgoicpyIDBKRS0XkOKeDMiZZQt215+dk1bmYvvfee/Tu3ZuxY8cycOBA1q1bx913310rIUDw4m8AOVkuR9dapMtucia1RLN47QZgOe59mq8AlonIr50OzJhkiKZ6565du/jpoKH8+Mc/5vOy3XT/1R/41f/8mY4d6w4kQ/CFd49e1YtV91/k6F17fctbm6Ypmu6jMUBvVd0NICJHA/8HPONkYMYkQ7g1BzU1NcyaNYu7xvw333+/nzZnX0nb/3cVFa6WEaeWJqMkRCoXETSpK5qksBXY53e8j9qL0kwT0hT6qINdwFetWsVNN93EsmXLaHNiTzpcOQpX+yMtg1TcSjKViwia1BXNmEIZ8KGIPCAi9+OeOrpBRO4UkTudDc+kkqbYR71v3z7uuOMOzjzzTDZu3Mhzzz1H7hV/qJUQvFLtDtw2sjHxiKalsNHz4/W6539t850mxumNbhLdCqnP+6kq8+bN47bbbmP79u3ceOONTJo0idzcXJ6YvDgt7sCtvLWJR7hNdp5X1V8C5ao6vQFjMinKyT7q+m4En8j327hxI6NHj+af//wnvXv35rXXXuOss46U7erXPY/Zy0pr7aGQqnfgVt7axCpc99GZItIJ+LWI5IpIO/+fhgrQpA4nN1tP9EyZeN7v0KFDTJgwgdNOO43333+f6dOns3z58loJYX5JGfNWltVKCAL84kxnLr7zS8roO3kxXca+Sd/Jixt1V51JDeG6j54A/gmcCKykdoE79TxumhAnN7pJdCsk1vf717/+xc0338znn3/OVVddxdSpU+nQoUOd84IlGwXe/WxXXHGGk+jWkzHRCNlSUNX/VdVTgWdU9URV7eL3YwmhCXJyo5tEt0JCvU6h1h33jh07uPrqq+nfvz81NTUUFRUxZ86cWgnB/269PoXyYmXrDEwyRBxoVtWbPLuoHet/vqqWOhmYSU1O9VEnuhUS7P28ysorGDt3FYvmPMvLT0zh4MGD3H///YwdO5aWLVvWOjfwbj2U+nShhRoQt3UGJhkiJgURGQ08AHwN1HgeVuAHzoVlmppEz5Txf7/Au/tD279g+1szWL9jA21PPpO8fiN5u+XJ9Fq3u87nhSuQ51Wf5BWui8jWGZhkENXQJX0BRGQDcJZ3RXOyFRYWanGxVdc20esy9k0UqDn4PXvfe57vSxaReVQuuT+9gezu5/m2xMxyZdbpDvO+NhiBeievviGmt3q38QzWekqVvalNehGRlapaZ7+aQNGsU9gCfFv/kExjlA4rnNu0bEbZyrfZ++7T1Bz4jtZnXkLOeb8ko0V2rfOCrbkIdbeen5PF+2N/Wu/YwnUR2ToDkwzh1il4Vytvwr072pvAIe/zqjrV4dhMikuH2TEz/v4eX/z1Hio2r6b58afQ7ooHaHHcySHPD7xIxzLWEU+CjNRFZOsMTEML11Lwrlgu9fw09/wYA8S3wrmhWhYVFRVMmjSJiQ89jDRrTruLbuaongOQDHfZhwyBmiD9QoH99dHercebIJ2c5mtMPEImBVV9sL5vLiIDgelAJjBLVScHPD8KuAWoBr4HRqrq2vp+rgkvURfmcF0fwT4DaJCWxaJFixg9ejRffvklrU7rR26/X5PZKrfWOTXqvvhGczGO5m493hIg1kVkUk00A82FwDigE7WnpIadfeSZxvo5cCHuSqsrgOH+F30RaaOq33l+HwTcrKoDw72vDTTXT7AplvEOXoYaJM3NdnGwsqbOZ7RolkF5RWWd8+vbP+9NQKWlWziw5Gn2fLqE7t2789hjjzH+Qw07kJuoi3GoAWkBvpx8cVzvaUwiJXKgeTbuPRU+4ciU1Gj0ATao6iZPQHOAwYAvKXgTgkcrCDnRwyRIIovaher6UCXoZ4Sa2llWXkHfyYvjujjPLylj7NxV7Pzg75S//yLUVJPX73omTL6Pfn26MCYneBL0fkai7shDjQ1kiDC/pMzu/E3aiKZ09i5VXaCqX6rqZu9PFK/Lp/a+C1s9j9UiIreIyEbgEeDWqKI2cUvkgqhQK5y/DdIaiMS/HPeYV1dHXePnvpmvsemp0ex992ladjyd43/zGNl9ruDRxV+GjTHRF+lQW25Wqzb68uKmcYmmpXC/iMwC/kXt2UevRXidBHmsTktAVWcAM0TkamA8cF2dNxIZCYwEKCgoiCJkE0qiF0QFu9sOtmAMoFXzTPYfDr8QDKCyRnlgwZqwF+7du3czduxYPpk1i8zW7cm77F6yup7jW3Pgn+QaYgaP9/3vemU11QFdsqm4AY8xoUTTUvgvoBcwELjU83NJFK/bCvjvRHICsC3M+XOAIcGeUNWZqlqoqoV5eXlRfLQJpSE2Xgn1Ga7MaP65uQUbewD3PgfPPvss3bt359lnn6XDeUPpcMPjZJ/y/3wJAaBtlqvBq4sO6Z1PTYgxOitNYdJFNC2Fnqp6RhzvvQLoKiJdcO/eNgy42v8EEemqql94Di8GvsA4qiFmu4T6jNtfXlWv9/3000+56aabWLp0KX379uXKWx/gr5/VsPdA7QTiyhD2H67yJZay8grGzF3NAwvW8G1FpaMzfKw0hUl30SSFZSLSI9apoqpa5ambVIR7SuozqrpGRCYAxaq6ABgtIv2BSmAvQbqOTOI1VHeK/2fMLylDiH4mQW62y/f7/v37mTBhAlOnTqVt27Y8/fTT5PS8kHHz19QZvM7JciFCnURRWa21koRTi+xs3YFJd9EkhXOB60TkS9xjCgJopCmpuE9aBCwKeOw+v99viy1ck66mFK0PmRACF5K5MoX7Lz0NgNdff51bb72V0tJSfvOb3zB58mTat29P38mLg85matWiWVRdNU7189u6A5PuokkKYdcNGBOK/wK2cC2EGoVMEapVfesHeuVWMmjQIBYuXMjpp5/u6zLyCjeLKlQXTrBznWClKUw6i2Y/hWimn5pGJBErnqPdh8CrWpUsVya3/7QLnxU9z9UTJpCRkcGUKVO47bbbePPTnbXWMuRku+p0EcGRqqVO74FgTGMVTUvBNCGJKnIXzT4EgfZuXMV1l46iYudmLrvsMqZPn07Hjh2DxuTKEFyZQmX1kTaI/6I0bwzbyivIDjIV1vr5jQnOkkIjFs8df6JWPMfSNVO9v5y9/36G/Z8uplnbY3njjTe4+OIjpSGCxVRZo+RkuXxjCIF/n7cLx5tQ/AnwizOti8eYYCwpNFLx3vEnasVzNP36qjV8v7qI8v/8jZrDB2lzzlBOHXgdF1/886g++9uKSlbdf1HYzwiWUBR497Ndkf8IY5ogSwqNVLx3/ImaZx+pX//w15vYXTSDw9vX06LgDI6+8GbaHN+ZsZee4WvhlJVXkCkScpA6mphsn2NjYmNJoZGK92IYzzz7cN1U/o/3657HO6s3s+6NWXy3ciEZZOpZAAAXY0lEQVSt2uTQaehYDnfuS35udtDy2oElI6KNycsWkxkTG0sKjYz3Ah3L3XXgRf0XZ+bz7me7ohqLiNRN5X2dqjJ37lyenXk7+7Zv56ZRo5g4cSK5ubX3OQi1/sBffgwzomwxmTGxsaTQiESaBhrsYhjsoj5vZVnUlUSj6abasGEDo0ePpqioiN69e/P3v/+dPn36BH2/SC0ZgZj2XrDFZMbExpJCIxJuGmiou+v6zjYK10116NAhHn74YSZNmkTz5s2ZPn06N998M82ahf5nF2mAOp5uH1tMZkz0LCmkscBun1AX03B315HGHiJNaw31uZnbP+GMM+7giy++4KqrrmLq1Kl06NAh4t8UaYB6/6Eq27TGGAdZUkhTwbp9QvG/uw68yLfNcgUtU90hJyuqaa1jBnRjzNzVvkVkVd/vYe/iWRxY9x7Hd+zMW2+9xYUXXhj13+Xf3VNWXoEI+I81l1dUOlbMzhgTxR7Nqcb2aHYLtT9yoMwMoXWLZnxbUUnbLBf7D1fVWgUcjHfP5lCb5QTuqdzrwbfYu/8g+0oWUf7e82h1JW3PvpJTB/ySD/6nfqWzQv2d9d3X2ZimJpF7NJsUFO08++qaIyWjQ21c4y8328X9l57GkN753BFi/4PAz965aQ2733qMwzs20LJzb9pdOApXu3x27A/eBRTLSutQf6d/okhErSZjjJslhTQVbSXQWKm6u27ueHkVGZ7KpcE+G6C8vJxx48ax/fnHyWyVS/tB95Dd/VzfDmihpr/GstI61N8pnvcCElKryRjjFv3+iCalhNoovr7KKyop85S6DpYQslyZ3H3RKcyePZvu3bvzxBNPcMmwX3PiTTNpdep5voQQai1AuNlOwYwZ0C3kZt9TitbH/H7GmPAsKaSpIb3zeejyM8jPyUJw7zjmlEwRBHc//i29s/jLmOu49tpr6dSpE8XFxSx8cRaPDD/bF0t+TlbIdQ6xrrQe0js/5EK8beUVVsbCmARztPtIRAYC03FvxzlLVScHPH8ncANQBewCfm37N0QvcP5957FvOvI5Naqsvf+nTJw4kTuHPUKrVq14/PHHGTFiBJmZmUFjmV9SVmv/A28/fzxlJ/IjvMbKWBiTOI61FEQkE5gB/AzoAQwXkR4Bp5UAhZ6tPecCjzgVT1OQH+FCGOn5UA5sXEFu/klMnDiRYcOG8dlnnzFq1ChfQgjkHTfwdkN5+/nnl5QF7faKVHYi3GvieT9jTGhOthT6ABtUdROAiMwBBgNrvSeo6rt+5y8DrnUwnkZvzIBu3PHyqqDdLZkibPNUHQ1VZC5Q1XffsPdfMznw+f/hOrojHa99mMvvuIZjjz027OvC9fN7p5HGMlsomlIVNvvImMRwMinkA1v8jrcCZ4U5/zfAPxyMp9Eb0jufV4tLeX/jnjrPeRNBsITgyhRQ98Y1AFpTzb7iBZQvnQ2q5Jz/K9r0uQzJdNUqfxFqKmikfv54yk6Ee42VsTAmcZxMCqEmjdQ9UeRaoBD4cYjnRwIjAQoKChIVX6P01e7oBlgzRahR9V3MAR5cuIbtn3/MnrdmULnrK7JO+hG5/W/ElXOc73X+5S+CTQUt3rwn4lRWY0zqcjIpbAU6+h2fAGwLPElE+gPjgB+r6qFgb6SqM4GZ4F7RnPhQG49oZ93UqPLl5CNbXu7evZvS16ex56N/ktm6PXmXjSOr69m+KaZebT2znEJ1Ec1eVho081s/vzHpwckpqSuAriLSRUSaA8OABf4niEhv4ElgkKrudDCWJiPau3HveTU1NTz77LN069aNPavepk2fy+lww+Nkn3JOnYQA4H0oVPIJNZ4RbSluY0xyOVr7SER+DjyKe0rqM6o6UUQmAMWqukBE3gHOALZ7XlKqqoPCvafVPgpvfklZyMFmL8F98c45uIPD/5nJulXLObVXH8p/+Cua53WO+Bn5OVkcOFzF3gORy2Z4P8+/VWKMaXgpUftIVRcBiwIeu8/v9/5Ofn46iqWOT7BzIcTAjZ/qwwf59v0X2Vz8OpktWnHLfX9kqZxG84PhdzzzKiuvwJUhuDKlVnE9b7IJZGMJxqQPq32UQmKpCxTs3DFzV1NdEz4lHPhiGXvefpLqfbs46gcXkfOT63njUBsguoTgVVmj5GS5aNWiWa09mOetLLOtL41JY5YUGkC0d/+x7IIW7NxwJbGrvv2aPe88ScWG5bjyOtN+0BhanhC4ljA231ZUsur+i2o9VtipXdC/1SqZGpMeLCk4LJa7/1jq+EQ7y0irK/lu+d/59v9eBhFy+/2a1mcOQjKj+79eBNq2DL0RT6BgawZirYxqjEkeK4jnsFiqeIbqew98fH5JGRlBZgYFOlj6CdufvZXy954j68Qz6XDD47Tpc3nUCcGVIUwb2osHBp1Wr1ISVsnUmPRhLQWHxXL3H2x/4sCLr/euO1ypiur95ez99zPs/3QxzdoeyzFX3E/WST+KOfar+nQMWkoiJ9uFKtzx8iqmFK2P2BVklUyNSR+WFBwWS1XQSDV+5peUcdcrq0MmBNUavl9dRPl//kbN4YO0OWcobc8ZSoarZVyxz1tZRmGndr4uIe/YQKxdQfFURjXGJIft0eywwIuol/+2l6Fe98CCNb6+/GxXBpU1GnIw+fDXG9ld9BiHt6+nRcEPOPrCm3C17xj03FjkZrsoue/IYHI8eyYH+w68+0DbmIIxDSMl1imYI3fP/hd4gL0HKn132FC7ddCvex4vL9/iK1AHcKCyJuj71xw6QPmSF9j30RtkZLXh6EvuolWPnwRdjRyPvQcqmV9S5vs7Qm0BGq4rKJoqp8aY1GAthQYS6g47J8vFoaqaWnfRoRaB+VNVDny2lL2Ln6L6+70c1ftn5Jz/KzJbHpXYwDnSCgi3WjpcS8EYk3zWUkgxoe6kg031jJQQKvduY8/bT3Dwy49ofuxJ5F02jhYdnFsg5o19StH6oLEJ2AI1YxoJSwoNJNRgayy06jDffjiPbz94BclsRm7/G2nd++dIRvAd0BLFOyAcrgiedQUZ0zhYUmgg/brn8cKy0rhfX/FlCXvefpyqvdvIPvV8cn96A82OapfACINzZYivFRAqscW7zacxJvXY4rUG8u5nu+J6XdX3e9i14BF2vvI/ABwz9PfkDfrvhCQEwX1Bv/bsAt+FPcNvfDony8WUK3v6WgG2H7IxjZ+1FBIs1i0qQ9GaavaVLKL8vefR6krannsNbc/6BdKseULi9B8Ynl9Sxruf7UKA49uGnhlks4iMafxs9lEChZuPP6VofdRjCoe2f86eohkc/nojLTv3pt1FN+HK7ZDQWEVA1b0O4fuDVbWmv9oaAmMan2hnH1n3UQKFq/HTr3texNfXHPye3W89xo7n7qJ6/17aD7qHY4ZOqHdCyHZlELhswXsvsPdAZa2E4B+zMabpcbT7SEQGAtNx77w2S1UnBzx/Pu6d2X4ADFPVuU7Gk2iBXUWhWgJl5RXMDjPIrKrsX/tv9i5+mpqK72h95qXknHctGS2yExJnZY0Sa4PQ6hIZ0zQ5lhREJBOYAVwIbAVWiMgCVV3rd1opcD1wt1NxOCGwBAWEXunrFeqaXLl7C7vfepxDpR/T/PhuHD30QZofe1ICow2/z0IobbNcCY3BGJMenGwp9AE2qOomABGZAwwGfElBVb/yPBe8hkMKClXLKFbVFd+x9X+vBiCjRSvaDbiFo3oOQCQ1evQSVCXDGJNmnEwK+cAWv+OtwFkOfl6DCDZuEKs9bz/Ovo/e9B13GPEkma1y6htaQpUfqLvS2hjT+DmZFILda8Y11UlERgIjAQoKCuoTU73Vp6/9UNk6drwwxnfctu/V5Jx7dSLCSjgra21M0+RkUtgK+NduPgHYFs8bqepMYCa4p6TWP7TIQq03iKdcRU3lIbbNHEH193sAkBatOOHmv5LRPDUvvLYgzZimy8mksALoKiJdgDJgGJCat8UBwm0k0/no2JLCtx++Rvm/n/EdHzt8Ei0LfpDYgGMQqQJrvi1IM6ZJcywpqGqViIwGinBPSX1GVdeIyASgWFUXiMiPgL8DucClIvKgqp7mVEzRCrXe4PaXV0X9HpV7ytj21I2+41an9+fon9+WsH0O4iHANWcX8NKHW4Lu3mblr40xjq5TUNVFwKKAx+7z+30F7m6llFKfcQOtqebrl37Hoa1HZt6eMPp5MlvlJiK0iHKzXVz8g+OZt7Kszh4N15xdwB+GnBFyzYStTTDGWO2jAPNLysgQCbkPcjj71/6HbxZO8R23v3QMrXr8OJHhRZTdvBl/GHIGgK9FkCnC8LM6+h63PZONMaGkxqT4FDG/pIwxr66OOSFU79/L5ocv8SWEFh1Pp2DM6w2eEMB9tz9+/ifMXlbq+zuqVZm3soz5JWWAVTs1xoRmLQU/DyxYU6cOUDiqyu5F09n/6Tu+xzqMeBJXu+QN0rbNcjF7WWmdwWRvPaMhvfOt2qkxJiRLCn6CbY0ZysHSj/n6pXt9xzk/+S/anvULJ8KKiUjo2UX+XUb+ycEYY7wsKcSo5nAFW2dchx4+AEDmUe3oMPIpMlwtkhyZ294wK5EFdxeZJQNjTCiWFPxkCITrPSpf+iLfvv+i7/i4a6fQIv/UBogsMRR8XUjGGBOMJQWOrF4OlRAO7/qK7c+M9h23/uHFtLvwpgaKLrFs2qkxJpwmnxTCVT3Vmmp2PHcnh7/e6HvshFtfJDOrTUOGmFA27dQYE06TTwqhqp5+//Fb7P7H//qO8y4fT3bXsxsytHrJyXJxqKqmztagNu3UGBNOk08KgYu4qr77hrLHr/cdZ51YSN4V9ye1PEU8Lul5PIWd2tm0U2NMTJp0Uhg//xPf76rKN69P5sD6932P5Y96mmZtj01GaPX27me7+MOQMywJGGNi0qSTwksfuvcAqti0kp2v3u97vN1FN9O698+TFVZC2ICyMSYeTTopVB78ni3Th4O6dwNt1i6fDr/+C5KZ/vsT24CyMSYeTTYpvPTSS2x59Mj2Dsdd9ygtjjs5iRGFluXKoF2rFmwrryAn24UqfFtRGXZfBBtQNsbEo8kmhSVLlgDQ5qwryP3J9ckNJoKKypqgg8V9Jy8OWu00J8tlYwnGmLg0ySqp1zz1AYvaXEyne95I+YTgpRzZAS5StdMHBiV9nyJjTJpqcknhmqc+4P2Ne5IdRty81U7BXdTuocvPID8nC8G9c9pDl9uMI2NM/BztPhKRgcB03NtxzlLVyQHPtwCeA84EdgNXqepXTsaUzgnBa5tVOzXGOMSxloKIZAIzgJ8BPYDhItIj4LTfAHtV9WRgGvCwU/E0JjazyBjjFCe7j/oAG1R1k6oeBuYAgwPOGQz8zfP7XOACSbelww3MSlUYY5zkZFLIB7b4HW/1PBb0HFWtAr4Fjg58IxEZKSLFIlK8a9euegXV96R29Xp9ouVmR78mwsYMjDFOczIpBLvjD5xaH805qOpMVS1U1cK8vLx6BTV7xDkpkxhys12U3HcRj17Vq9Zg8bVnF9Q6fvSqXnw1+WLeH/tTSwjGGEc5OdC8Fejod3wCsC3EOVtFpBnQFnB8JHj2iHPqPNbQs5JcmcL9l7qnjtpgsTEmVTjZUlgBdBWRLiLSHBgGLAg4ZwFwnef3K4DFqhpuoa5jnGpBZLsyePSqXnVaA1Ou6GmJwBiTchxrKahqlYiMBopwT0l9RlXXiMgEoFhVFwBPA8+LyAbcLYRhTsUTjWAtiESyJGCMSXWOrlNQ1UXAooDH7vP7/SBwpZMxGGOMiV6TW9FsjDEmNEsKxhhjfCwpGGOM8bGkYIwxxseSgjHGGB9LCsYYY3wsKRhjjPGxpGCMMcZHklRVIm4isgvYHOG09sA3DRBOfaVDnBZj4qRDnOkQI6RHnKkWYydVjVhRNO2SQjREpFhVC5MdRyTpEKfFmDjpEGc6xAjpEWc6xBiMdR8ZY4zxsaRgjDHGp7EmhZnJDiBK6RCnxZg46RBnOsQI6RFnOsRYR6McUzDGGBOfxtpSMMYYE4e0TgoiMlBE1ovIBhEZG+T5FiLysuf5D0WkcwrGeL6IfCQiVSJyRUPH5xdHpDjvFJG1IvKxiPxLRDqlYIyjROQTEVklIktFpEdDxxhNnH7nXSEiKiINPkMliu/yehHZ5fkuV4nIDakWo+ecoZ5/l2tE5MWGjtETQ6Tvcprf9/i5iJQnI86oqWpa/uDezW0jcCLQHFgN9Ag452bgCc/vw4CXUzDGzsAPgOeAK1L4u+wHZHt+vylFv8s2fr8PAv6Zit+l57zWwHvAMqAw1WIErgf+kox/jzHE2BUoAXI9x8ekYpwB5/8W9y6USfleo/lJ55ZCH2CDqm5S1cPAHGBwwDmDgb95fp8LXCAikkoxqupXqvoxUNOAcQWKJs53VfWA53AZcEIKxvid32ErIBkDZtH8uwT4PfAIcLAhg/OINsZkiibGEcAMVd0LoKo7GzhGiP27HA681CCRxSmdk0I+sMXveKvnsaDnqGoV8C1wdINEF/D5HsFiTAWxxvkb4B+ORlRXVDGKyC0ishH3BffWBorNX8Q4RaQ30FFV32jIwPxE+//3LzzdhXNFpGPDhOYTTYynAKeIyPsiskxEBjZYdEdE/d+Op8u1C7C4AeKKWzonhWB3/IF3htGc46Rkf360oo5TRK4FCoEpjkYU5KODPFYnRlWdoaonAfcA4x2Pqq6wcYpIBjANuKvBIqormu9yIdBZVX8AvMORFndDiSbGZri7kH6C+w58lojkOBxXoFj+Gx8GzFXVagfjqbd0TgpbAf+7lxOAbaHOEZFmQFtgT4NEF/D5HsFiTAVRxSki/YFxwCBVPdRAsXnF+l3OAYY4GlFwkeJsDZwO/FtEvgLOBhY08GBzxO9SVXf7/X/8FHBmA8XmFe1/36+raqWqfgmsx50kGlIs/y6HkeJdR0BaDzQ3Azbhbo55B3hOCzjnFmoPNL+SajH6nftXkjfQHM132Rv3gFrXFI6xq9/vlwLFqRhnwPn/puEHmqP5Lo/3+/0yYFkKxjgQ+Jvn9/a4u3GOTrU4Ped1A77CszYslX+SHkA9/w/5OfC552I1zvPYBNx3sgAtgVeBDcBy4MQUjPFHuO829gO7gTUp+l2+A3wNrPL8LEjBGKcDazzxvRvuYpzMOAPObfCkEOV3+ZDnu1zt+S67p2CMAkwF1gKfAMNS9f9v4AFgcjLii/XHVjQbY4zxSecxBWOMMQlmScEYY4yPJQVjjDE+lhSMMcb4WFIwxhjjY0nBNDkikiMiN9fj9RM8C/nqG0dnEbm6vu9jTCLZlFTT5HhKqL+hqqcnOY6fAHer6iXJjMMYf9ZSME3RZOAkT337aZ79IT7y7MUwGHx38etE5ClPrf63RCTL89xfvXtfiMhXIjJJRD4QkWIR+aGIFInIRhEZ5TlHRGSKiHzq+Yyr/OI4zxPHHSLSUkSe9ZxTIiL9kvDdmCauWbIDMCYJxgKnq2ovT02sbFX9TkTaA8tEZIHnvK7AcFUdISKvAL8AXgjyfltU9RwRmYa7XElf3Kvp1wBPAJcDvYCeuMsxrBCR9zxx+FoKInIXgKqeISLdgbdE5BRVTUZ5bdNEWVIwTZ0Ak0TkfNx7WuQDx3qe+1JVV3l+X4l7Q6RgvEnkE+AoVd0H7BORg56qnecCL6m7OubXIvIf3OVNvgt4n3OBPwOo6mcishl3eeiP6/k3GhM1SwqmqbsGyAPOVNVKT+XSlp7n/CvBVgNZId7De15NwGtqcP83Fu3GTg25AZQxQdmYgmmK9uEuYQ3ucuo7PQmhH+DE3tPvAVeJSKaI5AHn4y7Q6B+H97xrAETkFKAAdzloYxqMtRRMk6Oquz27dX0KrAC6i0gx7uqqnznwkX8HzsFdcVSB/1bVHSKyG6gSkdW4xyIeA54QkU+AKuB6bfh9K0wTZ1NSjTHG+Fj3kTHGGB9LCsYYY3wsKRhjjPGxpGCMMcbHkoIxxhgfSwrGGGN8LCkYY4zxsaRgjDHG5/8DmtF1mmHDI4oAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.scatter(tanimoto_sim, mhfps_sim)\n", "plt.plot(tanimoto_sim, y2, color='black')\n", "plt.xlabel('tanimoto')\n", "plt.ylabel('mhfp sim')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### MorganFP / Tanimoto" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAyCUlEQVR4nO3deVxU5f4H8M/MsC9iISqKYmqigJqiuIC5Zyq5IV0zyTUsu879paldrUu7lMsF0wxbAL2CUmZC7opXERXEpSsupYCZioqA7Pt8f38cGkZAhJlzZsD5vl/+MXM48zwPepzvOc/yfWREBMYYY8xYyQ3dAMYYY8yQOBAyxhgzahwIGWOMGTUOhIwxxowaB0LGGGNGjQMhY4wxo8aBkDHGmFHjQMgYY8yocSBkjDFm1DgQMsYYM2ocCBljjBk1DoSMMcaMGgdCxhhjRo0DIWOMMaPGgZAxxphR40DIGGPMqHEgZIwxZtQ4EDLGGDNqHAgZY4wZNQ6EjDHGjBoHQsYYY0aNAyFjjDGjxoGQMcaYUeNAyBhjzKhxIGSMMWbUOBAyxhgzahwIGWOMGTUOhIwxxowaB0LGGGNGjQMhY4wxo8aBkDHGmFHjQMgYY8yocSBkjDFm1DgQMsYYM2ocCBljjBk1DoSMMcaMGgdCxhhjRo0DIWOMMaPGgZA1JxcvYtOmqtc5OYiKwuXLOHCg6khmJrZuNVTTGGPNFQdC1pycPo358xEdDQD37mHDBpw7h+3bq356+zbWrTNg6xhjzZI2gbD2XTnAN+ZMT8aMwdKlyM01dDsYY08KbQJh7btygG/MmZ64u2PGDCxfXn1kxw64u8PdHb6+hmsWY6zZ0rJrlO/KmQG99x7i4pCYWPV20iScOYMzZ6o6JxhjrFFMtPuYuzv69cPy5VAqqw/u2FH13VRSAnt7MVrHWF0sLPDVV/DzQ/fuAKBQwNwcAMzMDNsuxlizpP1kmRp35eAbc6YXMhkADB+OCROqXgvKy6te3L2LZcvw3nsoLDRA8xhjzY6WT4SodVcOvjFn0ps1q/r1999XvXBxQbducHbGwYNITMS8eVi6FIcP4+ef8eqrhmglY6xZ0fKJ8FF35Zr4xpzpR/v2uHoVSUlQqQDg229hZYXYWAwYYOiWMcaaAxkRiVVWWRk2bsTJk4iKgkxWfWPeogXfmDNpOTvjxg1cuoQePXDkCDZtQlAQnJ0N3SzGWHMg5oJ6U1OsXo3t23H+PMA35kyPhGtMGLFeuRIAAgNx8KAhm8QYay7EDIQyGXx8ACA2FgCOHMGSJdi4EV27ilgJY3XQDIQHDiAqCuHhGD3asI1ijDUPIqdYe+kl4K9AyDfmTG80AyFjjDWKmGOEAEpL0aoVCgtx4wacnEQsmLH6FBfDzg5EePAA1taGbg1jrFkR+YnQ3BwjR4IIe/aIWzBj9bG0hLs7Kipw7pyhm8IYa27E331Cs3eUMb3h3lHGmHYkCYRyOQ4fRlGR6GUz9kgcCBlj2hE/ELZujf79UVyMQ4dEL5uxR+JAyBjTjiQb83LvKNM/FxcMG7amqMgtIyPD0G1hjDUn0gZCIeUVY3ogl8PEZN/9+5eSkpIM3RbGWHMiSSDs1QuTJh2ytByTnMxfSUx/BgwYACCRu0cZY40hSSAE4OS06/r1A7/88otE5TNWGwdCxpgWpAqEL730EoCYmBiJymesNiEQnj59urKy0tBtYYw1GyJnllErKytzcHDIy8tLT0/v1KmTFFUwVtszzzxz/fr1lJQUNzc3Q7eFMdY8SPVEaGZm9sILLwDYwzlmmB5x7yhjrLGkCoT4q3c0lldRMD3iQMgYaywJA+H48eMVCsWRI0fy8/Olq4UxTZ6enuBAyBhrDAkDob29/aBBg0pLSw/yPkxMX/r27WtqapqSklJQUGDotjDGmgcJAyG4d5TpnaWlZefOnSsrK99///3bt28bujmMsWZAqlmjgsuXL7u6ujo4OGRkZCgUCukqYkywZcuWOXPmtGrV6s6dOwBcXV39/PymT5/erVs3QzeNMdZESRsIATg6OmZmZg4fPnz8+PFeXl59+vQxMTGRtEZmtEJCQhYtWqRSqcaOHWtubr5///7i4mLhR3369PH19Z0yZUqPHj0M20jGWFMjbSD8+uuvFyxYYGdn9+DBA+GIlZVVnz59vL29vby8hgwZ0rJlS+lqZ8aDiJYtW7Zq1SqZTLZq1arFixcDKCws3Lt3744dO/bs2ZOXlyec2aNHjylTpvj6TuvTx92gTWaMNRUSBsLIyEh/f38i+uijj5ycnI4fP56QkHDlyhX1CSYmJr179/b29n7++RGDB09o21aihrAnXFlZ2ezZsyMjI83MzCIiIqZNm1bjhNLS0vj4+NjY2O3bt9+9exfA0KHLrl8PmjgRfn4YPBhyacfKGWNNmlSBMDY21tfXt7y8PCgoaNmyZerjeXl5SUlJhw4dOn78eHJycmlpKYDnnht2/vwRR0d4e8PLC97e6NsXMpkU7WJPmsLCwqlTp+7bt8/GxubHH38cM2ZMPSdXVFT897///emnn/73v3kJCX2Fg+3aYfJkTJmCoUPBA9mMGSFJAuGRI0fGjRtXUlKyYsWKTz755FGnFRUVJSUlxcfH377daetWf83Vhg4OGDwYQ4bAywseHjA1xcWLSEhAQAAA5ORg3z688oroDWfNTFZWlo+Pz6lTp9q0abNnz56+fatiW2ZmZqtWrWSPvplSqXDuHGJjERmJq1erDj79NMaPh58fXngB167x9caY0SCxJSYm2traAliwYEHDP1VRQWfP0rp19Le/Ufv2BFT/sbKioUNpzRoCaPt2IqIrV8jLS/SGs2YmLS1NmAv6zDPP/P777+rjKSkpHTp0WL58eQPLSU6m5cvJxaX6kmvRgqZM4euNMWMh8thISkrKuHHj8vPzZ8yY8eWXXwoHCwsL9+3bV/8HFQr06YOFC7FtG27exK1biI6GUgkPD5SU4NQp2NpizBgsXYrcXHGbzJqllBRMm9Y2J8fMw8Pj1KlTzz77rHA8Pj7e29v7zz//PHHiRFlZWUOK8vDAp5/iyhWkpiI4GF5eyM+HhQVfb4wZCzG7RlNTU4cMGZKRkTFhwoQdO3YIyyRKS0snTJhw6NCh8PBwf39/LYrNzsalS7h2DSkpsLBATg6USsydi+PHxWr4EyUlJeXYsWNnz54NCAgQ8o09eeLiMHky8vLw2mv31683F3ogAMTExEybNq24uHjSpEmRkZGWlpbalZ+Whr17kZ7O1xtjxkGsR8ubN2n8+C8BjBw5sqSkRDhYUVHx8ssvA3BwcLh8+bIu5YeF0eLFVFxM3btTRAR3VdXtwIEDdnZ2wtiYXC4PCAjIysoydKNE9tNPZGFBAE2ZQsXF1cfDwsKEe685c+aUl5frWAtfb4wZD3ECYWYm9ehBpqY0Y0Zkfn6+cFClUs2dOxeAnZ3d2bNndaxC+GIiorg4srfnL6Y6fP/996ampgAmT568YsUKc3NzAE899VRwcHBlZaWhWyeO9etJLieAFi4kzd/pww8/Em7sPvzwQ1Eq4uuNMeMhQiDMzaV+/Qignj1J8/FDWNRsZWUVHx+vey1hYfTOO1WvX3mF+vfXvcgnh0qlCgwMFCKBUqkUwt7vv/+uXkvg4eGRmJho6GbqRKWiwEACSCajwMDq4xUV9OabNHToKYVCsXHjRrGq07zeZs8mb2+xCmaMNTm6BsKiIho6lADq2pVu364+/t577wEwMzPbu3evjlXUsHEjyeX097+LW2ozVlJSMn36dAAmJiZfffVVjZ/GxMR06NBB6Cn19/e/f/++QRqpo4oKCggggBQK2rSp+nhxMfn6EkCWlrR79xUpqj54kLy96Z//lKJsxliToFMgLCuj8eMJoPbtKT29+nhICA0detLExPSHH37QsX21JSURQB07kkoletnNT1ZW1vPPPw/AxsZm9+7ddZ5TUFAQGBhoZmYGoG3bthEREapm9XdXUlIV7ays6Jdfqo/n59Po0QRQy5Z09KhUtR88SAANGiRV+Ywxg9M+EFZW0rRpBJCDA126VH38m29IJiO5nKKjb4jQwFpUqqqFhr/+KkXxzcm1a9dcXFwAtG/f/ty5c/WffOXKlZEjRwo9pc8///yFCxf00kYRLFlCANnb08mT1QczMui55wggR0c6f17C2nNzSS4nc3P6awYYY+xJo2UgVKmquqpatKDk5OrjP/1EJiYE0Nq14rSvTq+/TgB98omEVTR9J06ccHBwANCrV68///yzIR9RqVQRERFt2rQR+lGVSmVeXp7U7WyslBQKDa16nZ1NkZGUn09TppDmpOPUVOralQDq3p3++EPyJvXoQQCdPi15RYwxg9AyEAo36ZaWD3VJHThA5uYE0KefitM4tdxcio6ufu6MiSGABgwQuZZm5IcffhAWyb344ouNDWY5OTlKpVLYHrJdu3YRERESNVI7YWGPT+myYQMBNHAgSTTimZhIISFUVlb1dtYsAmj9eknqau5ycnL27dsXExNTWFho6LYwpiVtAmF6OtnakpkZ7dtXffDECbK2JoD+8Q+x2lZt0SICaOnSqrdFRWRlRXI5ZWSIX1fTFxwcLJfLAQQEBGi9YO7s2bODBg0SekqHDx9+SbN326DCwmjMGHJ2pgcP6sttFh5O0n3xdu9OQHVXx8aNBJC/v1TVNV+JiYldunQxNTWVy+Xdu3fXfZUUYwbRoEBYu7cqKYl++qn6hNxcsrcngAICJJnDcuQIAdSjR/URHx8C6Ntvxa+rKSsvL3/jjTcAKBSKoKAgHUsTekqF/lVTU1OlUqleA2pAwgK+FStowQKDJfmcOZMA2rCh6u3ZswRQt24GaEmTVVlZ+emnnwoZDHr27CkkfbWwsAgJCWleU7EYowYGwob0VsXE0KxZJNG67YqKqkCrTq0cGkoATZwoSXVaKC8v//333xcvXhweHn716lUpqsjLyxs7diwAa2vrn3/+Waxis7KylEql8Ijp5ORk8J7SppDS5auvCKDXXqt6W15O1tYkk0nVE9vs3LlzR1iiKpPJlEplaWlpcXGxUqkUUhqNHj36tuZSKsaavIYGwob0Vklq+vSH5uDcvk0yGVlZSdg/1nBCDp0WLVqoE9e1adPGx8cnKCgoPj6+tLRU9ypu3rz53HPPCesfTkswbeP06dP9+/cXGu/j45OWliZ6FQ3UFFK6nDlDALm4VB/x9ibgobEAo3XgwIG2bdsKeRNrrNj5+eef7e3tAbRu3XrPnj2GaiFjjdXQQGjw3qqoKAJo+PDqI0I6G82FZYayaNEiAJaWlgsXLpwyZYowLVPNxsZm1KhRH3zwwaFDhwoKCrQo//z58+3btwfg5uZ2/fp10dsvqKioCAkJsbOzE36Xjz/+uEw9XUSPNFO6/O1v1K+f/ptAZWVkZUUyGWVnVx1ZvJgAEil9W3NVVlYWGBgodB6MHDmyzse+2g+L+m8nY43ViEBo2N6qBw/IzIxMTKq/mz76iLp3L/rXvww8q33FihVCDp19Gs8Lt27dio6OViqVHh4emtvDKhQKV1fXgICAiIiIPxo28X/nzp1WVlbCV8+DBw8k+z2qZGRk+Pv7y2QyFxcXPVRXj82bCaBp0wxTu5cXAbR/f9Xb6GgCaNw4wzSmKUhLSxs4cKCw8CYwMLCe7LUqlSo4OFhI4NCvXz/NrSIZa5oaEQjJ0AmIR4wggCIjq96ePXteWEtuwMH5kJAQIbzVk0Pn7t27MTExy5Yt8/LyEhJhqzk6Ovr5+QUHBycnJ9f5W6gniM6ePVufz2eTJk0C8Nlnn+mtxtouXCCAnnnGMLULE5U/+qjq7Y0bpUOGRPTr91rTmQmSlZWlt1zqERERNjY2AJydnRMSEhrykaSkpK5duwKwtbUNVc+1Y6xJamggbAoJiL/55vrAgV/Mnl2dZtTZ2RlAsuaSfj0KCwuTyWQymey7775r4EcKCwvj4+ODgoJ8fHyeeuopzaDYunVr9bBiSUlJRUXFwoULhS6mQM0k03rRu3dvAKJkS9daZSW1aEEA3bljgNp/+unesGH75s1bpz7Srl07ANeuXTNAax6WlpY2dOjQESNGODg4+Pv7x8TESHeTlJeXN2PGDOESnTp1ak5OTsM/m5ubq/nZbHVnDmNNTOPWEe7cSR4e0maNqUdqaqqwqZP6v/2bb74JQP9xgoh27NghrEn/97//rV0JFRUVZ86cCQkJ8fPzc3R01AyKVlZWQoy3sLCIiooSteGPl5ubq1AozM3NizW3+zOE4cMJoNhYA1R9/fp1APb29upHwIkTJwLYunWrAVqjITk5uXXr1gA0J2e1atVq7ty5u3fvFndMLjk5+dlnnxXGjIODg7UrRPNp8vjx4yI2jzGxNC4QbttGAA0bJlFjHs/V1RVAXFyc8HbPnj0A+vbtq+dmHDhwQOjkXLlypVhlag4ryuVyR0dHc3PzQ4cOiVV+wwl/q95NYOehd98lgN57zzC1C3MjU1NThbefffYZAKVSaZjWEBHR4cOHhfg3cuTIvLy8lJSUwMDAHj16aN5C+fj4RERE6Jg8T3Ocz83NTcfMtOnp6UL2BmF8saKiQpfSGBNd4wJhbi6ZmZFCYbAFVcuWLQPw9ttvC29LSkpsbW1lMlkDk22KouTYMWcnJwDLli2TqIrMzEx3d3cAmzdvlqiKegjTf9599139V13Dzp0E0OjRhql9woQJACL/GpQ+fPgwgAGGy+y3Y8cOCwsLAK+++mqNvtDU1NTg4GAvLy/1zCxLS0shImox4+nevXvjxo0TuuUDAgKKiop0b3x5ebl6xunw4cNv3rype5mMiaXRKdZGjiSADNU/dPz4cQBdunRRH5k8eTIAEXdkfYzTp6lFi8udO7+/eLGk8ya+++47Q33tDh06FECsQXokH3b7NgFkZydVoob6ffrppwD+8VfOwLy8PKHHuMQQ+1CsW7dOiCLqjZfrlJ6eXiMiKhQKLy+v4ODgjIYlJDx8+LAwGtqqVauYmBjxfgMiokOHDgmFt2zZcruQoYOxJqDRgTA42JCT2isrK4VVeurcmN9//z2AcfqZ2P7bb9S6NQE0fbrU381FRUXC2uSkpCRJK6qhrKzMyspKJpNlZWXps95H6dCBgIe2ntCbQ4cOARg4cKD6iJubG4DExER9NkOlUgUGBjZ22tSff/4ZGhrq4+NjampaIyI+6mlM86Ft2LBhEj203bt3z8fHR2iSv78/p+pmTUGjA2FqatVNuiEWWxMRzZw5E8Dnn38uvL17965wny55nsw//qCOHQmgl17Szy+/ZMkSALNmzdJDXWonTpwA4O7urs9K6zF1KgEUHm6AqvPy8uRyueYj4Jw5cwCsW7eu/g+KqKKi4vXXXxdG177VKrXu/fv3IyIifHx8hAE/AHK53MPDIzAw8LffflOfdv36dS8vL/0M46lUqtDQUGH7lB49evzKO4syQ9Nm9wk3NwLorwkr+vbjjz/a2tpq3hoL4/Aipt+sw61b1LlzVW4bfc2lvH79uhDj7969q58aieiLL74A8MYbb+itxvqtWkUAvfmmYWoXJmepH8pDQ0OFITo9VV9Y+NOCBQCsra337t2rY2HZ2dmbN2+eOHGiEIEEffv2/eSTT9atWycs5unQoYPe1sykpKT07NlTmBodHBzcdBZoMiOkTSAU5vL9NWFF38rKymrMERem87344ov3JZrDk5lZFfw9PUm/O9kKnUi6bzTRcMIigS1btuitxvodO0YyGb34omE60BISEq5cuaL+jj5//jyArl276qPu+/dp4EACPhk3TtzO2KKiopiYGH9/f/UCDKET3tfXV89L/YqKipRKpTBFmQMhMyBtAmFCAgHUubPojdHSxx9/rF6c3rlzZ39//9DQ0JSUFHH+a+XmVmU1dXfX/2TZvXv3AujYsaN+ZpyrVCphV6b09HQ9VNcQBQVl9vZdTE1NRZm7qKOKigpra2uZTLZr1y5p/0Vu3aKePQkgZ2e6ckWiSoqLi2NiYkaPHi1k1JOolsfq0qULAE7SzQxIm0BYWUlt2hBATWEz161bt8rlcplM5urqKuTkVGvXrt3LL78cEhJy9uxZLb+2iopo2DACqEsXMsTOMiqVStjpbdeuXXqo7tKlSwDat2+vh7oaTth2oymsxY6IiFAoFMJ6dnt7eyGri/h5pS9erJoj5OZG0q8LysvLs7CwkMvlBtk7KT09XSaT2djYGDx7AzNm2gRCIpo1i+RyCg7OFbc1jRUTEyNMihM6D8vKyk6dOrVmzZpJkyYJTzZqtra2GbNm0YcfUlxcI7ZuOnKEzMzIyYkM94S0du1aAC+88IIe6tq0aROAaYaaE/wI8+fPB7BmzRrDNmPlypXCmgRvb28hi6agZcuW/v7+O3fuFOeZ9dQpatWKABo6lBqTz0wXL730EoCvv/5aP9VpWr16dRO85Jix0TIQxsRcaN26l2GTj8TFxQnri5cvX17nCampqREREQEBAa6urpampiprawIIIIWCXF0pIIAiIujGjcdUs2+fYZ98c3JynJx6Dxu258oVyRfTCTNyv/zyS6krahRhhczLL79sqAaoVKp33nlHWMCgjse//vprYGCgMN1DYG1tPXXq1MjIyNxcbW8Qd+0iS8uqLaf12BUcHh4OYLQhMhcMHjwYQHR0tP6rZkxNy0CYn59vYWGhUCgyMzPFbVADJSYm2traAliwYEFDzs+6c4d+/JH+7/+oXz8yMamKiMIfFxeaM4fCwujmTcrMpDVr6M03KTiYtP46E1tAgAqgvxZ2S0gYrTl37pzkNTXGxYsXATg7Oxuk9tLS0ldeeQWAmZnZtm3bap9Qew17xqBBNGoUBQc3Ll94WFjVlTl7NpWXi/YLNEBOTo6ZmZmJicm9e/f0WW9GRoZcLrewsNAxIRxjOtIyEBLRiy++CCAiIkLE1jTQhQsXhHluM2bM0GYnmoICio+noCDy8aGWLasj4vr11K0bTZ1K4eE0fjy5uBhsseTDUlIIIFtbaaesZmRkAGjRokVTSwVZWVkpbBfcwNwoIsrPzxe2mbWxsdmv3pzwEf7444/g4OCXRo2qeqoDyMSERo2ir756/ABzaGjVRz74QLTWN4bw31m7pYpaW79+PYBJkybps1LGatM+EG7YsAHA1KlTRWxNQ1y7dk3Yq2HChAnlut84l5fT6dMUHEy+vhQc/FBey379yBA5r+s0ZAgBJGkiue3btwurUCSsQ1sjR47U24whtTt37vTt2xdA27Ztz54924hP3r9PERHk40NmZlXhTS4nDw8KDKRH7VJ76xZ16SLtP3C9vvnmGwBjx47VZ6UjRoyAgRLqMqZJ+0B448YNYbqXPlMv3rx585lnngEwYsQI8aeZzZ5NmsNjgYH0iNFH/RP2/XB1JelWWwkruj7++GOpKtDB8uXL6xkMlkJaWpqwA1Hnzp2vXr2qZSnZ2RQeThMmkIVFdceDhwf9+msdnfAGXR9y//59ExMTU1NTvS0l1H+NjD2KHNrq0KFDr169CgoKjh49qnUhjXL//v0XXnghPT19wIABu3btEmbKiKljR2RnV79NS4O7u8hVaGvKFLRvj0uXcOyYVFUICc29vb2lqkAHnp6eABITE/VT3YULF7y9va9evdqvX7+TJ09qzhFtnKeewsyZ2LUL2dmIiYG/P1q0wK+/wsYGXl44eRIDBuDgQXh6orwcGglf9M/e3n7YsGHl5eWxsbH6qVFYizlixIgaO1QzZgC6RNH3338fwN///vfHn6qz3Nzcfv36AejZs6dU+aDPn6cePUjYf+7OHWrXjhq/hY10AgMJID8/SQrPz88Xbs+bZhLku3fvArC1tdVmSLiR4uLihJQrI0aM0H7+56MUF1N8PG3b1gQ74Tdu3AhgwoQJ+qlO2Olp06ZN+qmOsXroFAiTkpIAdOzYUer0SEVFRcLeQF27dpV22e/ixWRiQgMGUMeOhsn0/Gi3b5OpKZmYSLLGev/+/QAGDRokftEiad26tUKheOWVV/bv318m2SQm9Z5/06dPF3+lvFqT7IS/c+fO2sGDb3p76yGPoHoVv/4nQDFWm/ZdowD69evXvn37GzduXLhwQbfn0vqUl5f7+fkdPXrUycnp4MGDwkwZqaxejdxc/Oc/SEvDzJkSVtR4jo6YPBkVFfjmG/ELT0hIQFPtFwWwe/fuzMxMU1PTqKioMWPGtGnT5uWXX968eXNBQYGItaxfv97Pz6+kpESpVG7ZskW9XYP4mmQnfJs2bd42NW1//Dh++UXqumJjY0tKSoYMGdK2bVup62Ls8XQMpPPmzQPw4YcfihKWa6uoqJg2bRoABweHS00hpZtBHT1KALVtS6I/qwjz96TdwUNbd+7cETahXLJkSWBgoIeHh/rqVe/DrmMfpnZ7/mmvyXbCr1tHAPn6Sl2Pr68vgJCQEKkrYqwhdA2EmzZtMjc3VygUHh4eSqUyOjpaxCX2KpVK2IytRYsWycnJYhXbrPXqRdbWJO5mveXl5TY2NjKZLDEx8YsvvhCzaJ1VVlaOHj0awAsvvKDugU9LS6uxht3CwsLHxyc0NFSLLat03/NPG02zEz4jg+RysrKiggLpKiksLBRyl//xxx/S1cJYw+kaCMePHw9A2NVaIJPJ3Nzc5s+fv2XLFh03MRB2prWysjp69KiO7Xxi/O9/omWgLCgoSE9PF7KzAujevXv79u0BrF+/XpwKxBAUFASgdevWdY4NC2vYR40aZWJiIlx+6n3Yb9261ZDyCwsLhb2urK2td+/eLXbzH1M3Xb1KTSx9AQ0eTAD98IN0NezYsQPAgAEDpKuCsUbRKRDu2rVLeFz7/fffDx48GBgYOGrUKMuHZ4E7Ojr6+PgEBQXFx8c3ao7DBx98AMDMzIz3ZyGilBQKDa16nZ1NkZF06RKpU53cu0f/+c9D5xcW0h9/UFIS7d5N4eG0ejWtXJkwc+bM8ePHe3p6duzYscZOHb169frggw+ioqJkMplCoWgifaTJyclmZmYymSw2Nrb+MzMzM2vvw+7l5RUUFFTPKsDs7GxhWPTpp59OSEgQu/nN05o1BJCUWbBfffVVAJ9//rl0VTDWKNoHwqKiImFt+7p16zSPl5eXJycnBwcH+/n51dgCwsbGxsvLa9myZQcPHqw/Vf+XX34p3N1v375d6xY+ScLCCCDhL+PKFfLyoq1bac6cqp+eP0+enrRhAw0YQJ06kTq7uOafYcMSagwPW1lZdezY0dPTc9SoUd98841QlHD/YWVlJe5msFrIz88XtqBatGiR5vGTJ0+Gh4c/aglNdnZ2RESEn5+ftbW1+jd1dXUNDAy8fPmy5pm3bt3q1asXAGdn5xo/Mmo3bpBMRjY2Ei3wLysrExYO/vbbb1KUz5gWtA+EK1asANCzZ8/685xpbgGhHtERxmPUw4o1Uv1u3rxZ2GJQ/e3MwsJozBhydqYHDx4ZCIWFhsIfCwvq0IE8PGjsWHrtNVq0iDZsuBYWFhYbG3vq1Kn09PSCggIiKi8vf/PNN9u1a3fjr104VCqVsAeFo6OjYYdw/P39hQusRgoh4fhju0ALCwuFfdiF5OzqiLhs2bL4+PiLFy926NABgJub25/S7/nXzPTvTwDt3ClF2bt37wbQu3dvKQpnTDtaBsKrV69aWFjIZLITJ040/FO3b9+Ojo5WKpV9+/ZVKBSaw4o9evSYN29eeHj4xo0bhfGe1atXa9e2J1JYGC1eTCtW0IIF1YHQzo7c3MjNjbp0IU9Pun6dTp6ka9coP7+hxZaWlgrzRfv06ZP/18fKysqE3J5ubm4PDDSbUUh8am1tXftZbfPmzZqDgnK5fMiQIcHBwY8K20VFRTt37pwxY0bLli01OycAPP/88zn62vOvOQkKIoBmzJCi7Llz50LKeeaMaUHLQDh27FgAc+fO1Tz43Xffde7c2d/fPzQ0NCUlpf5V9gUFBfHx8UFBQT4+PsLeAgLhNf8/qUEIhMXF1L07RURUBcKZM6mkhEpKKCmJPD21LDkrK8vFxQXA2LFj1Q/3ubm57u7uAMaMGSNCZvNGSk2lkSMDAHz33XePOqfhXaBqFRUV8fHxSqXSzs7OwcHBw8ODd0WvW1oaAWRnR2KnEa6oqBCGSy5cuCBuyYzpQptAKNytP/300zW6NAMCAjSHoBwdHadOnRocHJycnFz/l2l5eXliYuLatWuHDBkCwMXFRYtWPdmEQEhEcXFkb19316jW0tLSWrduDeD1119XH0xPTxdW79W43ZFaWRkNGEAAvf324YacX38XaJ0fWblyJYClS5eK2vAny3PPkVxOj/gL1FpcXByAZ599VtxiGdNRowNhXl6eMMk+VD2L8S8VFRUpKSmhoaH+/v7CAIyatbW1ME0mJiamns6oS5cucSCsU1gYvfNO1evZs8nbW8xASERJSUnCPNK1a9eqD548edLS0vL5Dh3yH54PJal33yWAOnSgxiaULS4ujomJCQgI0Jyi1alTJ6VSGR8fr9k/8cUXXwBYsmSJyE1/kpw+TQ1bf9Iob731FoB//vOfopfMmC4aHQiXLw8E4Onp+dj0x7du3YqOjq49TUahULi6ugYEBERERNQY18nMzBSeNRvbKqa76OhouVwul8t37NihPhi3c2dl69Ykl0u6sEztv/8lhYJMTOj4ce0LKS0t3bt37+uvv64ZEfft26c+QQiE76jvLFhttXeJ0plKpXJycgJw+vRpUQpkTCyNC4QpKfTUU3nDhi1qbJ6XO3fu/PTTT2+//banp6d6moPAT2M/BZVKZWJiIpPJpMuqzOoh9BlaWlqePHmy+ugXX1TNQ5V4pV12NnXsSAB99JE4BVZUVBw5cmThwoWurq6au2auWrWKA2F9Hjygbt1o6lQKD6fx48nFhcT4/3jixAkAHTp0kDpHP2ON9VBMqh8RFi5ETo6tm9sajXSPDdKmTZvJkydPnjwZQFFR0dmzZxMSEo4fP56QkNCpUyf1aTKZzN7e/u7du1lZWZyNV//efffd1NTUb7/9dtKkSVcSE1s6OwPAkiX44w9s2ICJE3HiBJ59VoqqiTB7Nm7cwJAhWL5cnDIVCsWwYcOGDRtW47jQP0FE4lTz5Nm3D87O+OEHAJg5E/3749gxrFuHlJR6PpTdpUv/1FT126KiotLSUs0TSkpKAPj6+mr2DzHWFDQiEP7nPzhyBK1b4+OPdarSysrK29vb29t72bJllZWVRUVFmj91cHC4e/fuvXv3OBAaxMaNGzNu3lxZVNRy7FgkJEDYNDUkBDduIDYWEyYgIQFPPy16vRs2YNcutGyJLVugsbKGGcL+/Zgwofrt+PGIi8OtW0hLq+dDChubtHpPALBhwwYhKSNjTUpDA2FeHpYtA4BVqyDihtIKhUJzsh8AYfqiMFjI9M/ExCR2+3aZtzcuX8bUqdi3D6amUCiwdSuefx7nz2PSJBw8CHNzESu9eBFLlwLA119DeAqVFD8RPkbtXaLGjsWCBXj4Ca8GC5ksVeOv1MrKyrzWRWJtbS3h5laMaauhgXDFCmRkwNsb/v6StgfCBId79+5JWw17NFmLFti7FwMHIi4OAQEICwMAW1vs3o2BAxEfj1mzEBkJkTq4SkowfTqKizF/Pv72N1GKZLqZPBmvvIIZM9C5M+7exeHD+PJLaCz2rZM50Fk/zWNMbA3amPfCBXz9NUxMsH69WN9+j8RPhE1C+/bYtQs2NggPxyefVB1s1w579sDODtu24cMPxaoqJweWlnB1xdq1YhX5GPxE+Bi9e2PcOLi4YOBAeHris88eGwUZa9YeHwhVKsyfj4oKKJXo3VvyBglPhBwIDa9vX2zfDoUC//oXtmypOujujshImJggPBx5eY0q7+JFbNpU9TonB1FRAHD5Mi5cQHw89u9HYSG2bhXxF2A6WL0aubn4z3+QloaZMw3dGsak9fhA+P33OHkSbdviX//SQ3u4a7QpGTcOGzaACPPm4ciR6oORkTh5Ei1aNKqw06cxfz6iowHg3j1s2AAA585h+3aYmsLJCbdvY906cX+ButnaDh4yZMvTT/OsjXpZWaFrV565xIzBYwJhdnbVXPbgYD31jrRvP6lfv0yV6gt9VMYea/58KJUoK8PUqfjtt6qDfn5wdNSisDFjsHQpcnPFbKAWCgoGxMfPyMoaaeB2MMaahsdMlrl7F23bolcv/c1iaNmybXKyuHMSmW7WrsX164iJgY8Pzp+HRpLrxnJ3R79+WL4cSmX1wR07kJgIACUlsLfXubUNIIxz8xAhY0xQ/URYewjn8mX8+SfOnkVUFDIz9TR+07o1APAQYROiUGDbNgwejIULq6Lg/ftYuxYLFiAkpLEjhe+9h7i4qsgnmDQJZ87gzJmqUUPGGNOz6kBYewhHGL8xMYGDg/7Gb4T0kDxE2LRYWuLo0arnuNxceHnh5EkMGICDB+HpifLyhpdkYYGvvsKiRdVHFAqYm8PcHHpbYMZPhIwxTQ+NETaFIZynnoKpKXJzUVZmyGawmtQZYtX5t2bOxC+/wNYWx441sAwhAg0fjgkTaq7DKSwUsa2MMdYID40R1h7CMcj4TatWyMjA/fto104fNbLGqTP/1sjHTzyZNav69fffV72YPh0TJmD+fJw7h4SEh7pMpcNPhIwxTTVnjdYYwjHI+A33jjZptfNvubvrUl5FBQ4cwOnTCAzUtWmMMaaFmoGwxhCO/sdvwPNlmrjJk7FtW1X+ZSH/1rhxupQnJNo2McHnn+PgQXHaWD9+ImSMaXooED5qCKeyUq9DOPxE2KRJkH/L2xvvvw+VCv7+uHNHlFYyxlhDVY8RPmoIp7AQ06cjO7vhUyJ0xU+ETd3q1fjoI9y+jWeeESvzyHvv4ehRxMVh9mzs2SNtSlt+ImSMaXp8irXcXJw7h+PHMWeOnr47hCdCDoRNmtj5t+RybN6MVq2wbx/+/W9xyiSCSlXHcQsLPPUULC3FqYUx1tw9PhC2a4e9e2Fnh8hIPaUb7d8fc+eif3991MWajvbtEREBmQzvvqvT9FGVCseP49138eyz+PnnOk6YOxfZ2Vi1SvsqGGNPkgZtw+TmVrWy/pNPEBoqTsWZmVi5Ev/7X83jFy/i+nV8+y0mTareo4AZiXHj8Pe/o7wcr77a2JQ1qKxEXBzeegtOThgyBJ9/jtRUHD780Dl17oDBGDN21GDffksAmZrSgQMN/1AdkpMpIIAsLQmguXNr/jQsjADavp2I6MoV8vLSqS7W7JSU0HPPEUD/93+ZDTm/rKxs//797713wcGBgKo/XbrQ0qWUlEQq1UMn89XFGKutQU+EgrlzsWQJyssxdSouXGh0xC0pQXg4+vdHv37YtAllZZgwAdOn13FmU0hwwwzF3BxRURg9OmH9+g6bN29+1GmVlZXHjx//xz/+4eTkNGbMmJiYjZmZeOYZKJWIj8fVq/j8c/Tv/9Ckm4oKqFR8dTHGapJRYybAEGHGDERGolOnqk0KG+L69etRUb+uWTMxKwsAWrXC3Ll44w106lTzzKIiREcjJQUWFsjJgVKJuXNx/HjDG8ieEN9///3cuXOtra3PnDnj4uKiPl5SUrJ///4ff/wxNjY2969o5u7u7uf3tylTVri71zHZtLISJ0/ihx+wfTteeQUKBV9djLGHNfYRsqiIBg0igPr3p8LCivpPjo+P9/PzMzExUSjMnJzKPTwoNJQKC2ueVlxMERH03HM0cyaFhdHixVRcTN27U0QEd14Zr1dffRVA3759S0tL1QfT09Nlfz3lubq6BgYGXrp0qc6PFxdTTAy99hq1bFndZTpuHF9djLGaGh0IiSgzk7p1o6FDj/j4+FRU1BELHzx4EBwc3K1bN+ELy9zc/LXXXktOvlH7zKtXadGi6q8qZ2f69ltavJiIKC6O7O35q8p45efnC5fQokWLNI/PnDnz888/v3btWp2fEuKfvz/Z2VXHP1dXCgykixerbrOIry7GmAZtAiER/fbbn08//TSAt99+u8aPgoKCbGxshBDYsWPHzz777N69ezXOqaykgwfJz48UiqqvKuFhsaiIwsLonXeqTps9m7y9tWsgexKcPn3azMxMJpPFxMTUf2Z+fv727dtnzy62tq6Ofx4etHIl/f579Wl8dTHGatMyEBLRsWPHzM3NAaxbt07zeGhoKAAvL6/o6Ojy8vIan3rwgEJDqUePqq8qc3Py86MTJ7RuBXvCBQUFAXBwcLh161btnxYWFsbExPj7+wv3XkOGXFI///32m/4byxhrlrQPhEQUFRUlk8kUCsXPP/+sPlhQUHD58uXaJ589e3bRonXCqgmAOnemVasoK0uX+tmTr7KycvTo0QCGDh2q7ofPysoKCwsbP368cCsGQC6Xe3t7h4buv1FHBzxjjNWncbNGa/vwww8/+OADKyurI0eOeHp61j6hrKxs165dmzZtOnToEAA3twJHR+uAAEyZImJ+LvYku3fvXu/eve/cufP+++937dr1hx9+OHDgQFlZGQC5XD5o0CA/Pz8/P792vH0lY0wrugZCIpo9e3ZERISjo+OpU6c6duyo/tHNmzdDQ0O/+eabu3fvAmjZsuWcOXMWLHi7SxcnXVvNjMy+ffvGjRtnZmZWWloKwNTUdPjw4b6+vpMnT3YQUtMyxpi2dA2EAMrLy8eOHXv48GE3N7eEhAQ7O7szZ86EhIRs27atvLwcQPfu3d9444158+ZZW1uL0WZmjKKioq5evZqUlOTr6ztx4kRhrhZjjOlOhEAIIC8vz8vLKyUlxc3NjYguXboEwNTU1NfX96233vL29ta9CsYYY0wK4gRCAOnp6QMHDrS2tk5PT2/btu3MmTPfeuutDh06iFI4Y4wxJhHRAiGAGzduPHjw4MqVK5MnTzY1NRWrWMYYY0w6YgZCxhhjrNlpxO4TjDHG2JOHAyFjjDGjxoGQMcaYUeNAyBhjzKhxIGSMMWbUOBAyxhgzahwIGWOMGTUOhIwxxowaB0LGGGNGjQMhY4wxo8aBkDHGmFHjQMgYY8yocSBkjDFm1DgQMsYYM2ocCBljjBk1DoSMMcaMGgdCxhhjRo0DIWOMMaPGgZAxxphR40DIGGPMqHEgZIwxZtQ4EDLGGDNqHAgZY4wZNQ6EjDHGjBoHQsYYY0aNAyFjjDGjxoGQMcaYUeNAyBhjzKhxIGSMMWbUOBAyxhgzahwIGWOMGTUOhIwxxowaB0LGGGNGjQMhY4wxo8aBkDHGmFHjQMgYY8yocSBkjDFm1DgQMsYYM2ocCBljjBk1DoSMMcaM2v8DBah28iruJOkAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Draw.MolsToGridImage(mg_sim_mols[0])" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAwJElEQVR4nO3deVxV1fo/8M8ZGBUFRVBLDVMTnEhUvII5myZgDihpZE5U9xr3plcx+3ax3+vbN7UyrprlNVMoTVAkQU3FGQdwTMXhmhOaIjMKyHjO8/tj0eEEynTO2Qfcz/vly5dnc85eC1juZ6+113qWgojAGGOMyZXS3BVgjDHGzIkDIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWjBYIL1zAypXl/87Nxfr1uHQJO3aUH8nMxPffG6soJmtVWxrAjY0xVn9GC4TnziE4GD/+CADZ2fjmG5w/j+jo8q9mZFRcvBgzRNWWBnBjY4zVnzGHRkePxqJFyMkx4ikZewJuaYwxIzJmIHR1xcyZWLCg4kh0NLp2RdeuGDPGiOUwuava0sCNjTFWX0aeLPPhhzh+HEePlr8cNw7JyUhOxvbtxi2HyV2llgZubIyx+lIb93SWllizBuPGwcUFAJRKqNUAyv9mzFgqtTRwY2OM1Zfxl094e2P8eCgUAMr/1klNxT//iYUL8eiR0YtlsqPf0gAoFCgsrGha3NgYY7WkICLTnT09HQsXIi8PW7YAwOzZCAnBwYOwsMDbb5uuWCZHS5fi44/xySf48EOAGxtjrNZMu6Dezg4bNyImBnl5ALB2LaysEBMDLy+TFsvkqF07lJYiKan8JTc2xlgtmTYQ2tjg5Zeh0ZRfng4exPz5+OYbdO5s0mKZHHl6AqgIhNzYGGO1ZPIUa+J+/NgxAPjsMwD417+wZ4+pi2Wy07EjHB3x4AHu3AG4sTHGas3kE+y8vLB8eXkg3LvX1KUx+VIo0K8fdu1CUhLat+fGxhirLYl6hCdOoKzM1EUxuas0OsoYY7Vh8kDo7IxOnZCfj+RkUxfF5I4DIWOsHqTYhkl0CvWTgDBmCp6eUChw5gxKS81dFcZY4yFdIBSPCRkzHXt7dO6MwkIefmCM1YF0gTAhQYKimNzx6ChjrK6kCISurmjZEvfu4e5dCUpjstavHwCcOmXCfEmMsWeMFIFQocDkyZe8vL5NTOQlXcy0Bgy4/9JLASdP9jZ3RRhjjYYUgRBAu3Zxx469d/hwnDTFMdnq0aPVnTuxly6dz83NNXddGGONg0SB0MvLC8AxnjDDTMzCwsLd3Z2ITp8+be66MMYaB4kCYb9+/aysrC5evPjw4UNpSmSy1a9fPwBJPGGGMVY7EgVCKyur3r17azSakydPSlMiky1PT08A3NIYY7UkUSAEj44yqYhAmJiYaO6KMMYaB6kD4VFOMMNMrGPHjk5OTunp6SkpKeauC2OsEZAuEHp7eysUisTExDJOv81MrG/fvuDHhIyx2pEuEDo6Onbq1KmgoODChQuSFcrkSYyOciBkjNWGdIEQgLe3N3h0lJkeB0LGWO1JGggHDBgAIDIyMisrS8pymdz07dtXqVSePHly48aNeXl55q4OY6xBky4QajSa/fv3Ozg4HD9+3NnZuU+fPosXL758+bJkFWAyodFoFi5cSERqtfrNN99s1arViBEj/v3vf6elpZm7aoyxhkhBJEV64pKSkjfffHPLli3W1tY9e/b89ddfS0pKxJdcXV39/Pz8/Pz69++vVEraQ2XPnuLi4qlTp0ZHR9va2k6dOvXq1avHjh3TarUA1Gr14MGDJ0yY8Prrr7du3drcNWWMNRRSBMLHjx9PmDBh9+7d9vb2cXFx3t7ejx492r17d2xs7C+//JKdnS3e5uTk5OPj4+vrO3LkKFtba1PXij17cnJy/Pz8jh496uDgEBcXJ1bsPHjwYPv27dHR0YcOHSotLQWgVCoHDBgwfvz4CRMmt2/f1ty1ZoyZmckDYW5urq+v79GjR52cnPbs2ePu7q7/VY1Gc+LEiR07dmzfvv3q1aviYOfONzt0cPHxwcSJeO45k9aOPTvu378/evToCxcutG3bdvfu3T169Kj0hpycnLi4uB07duzataugoACAp+e2vLxx/v6YPBmuruaoNGOsATBtIExLSxs1atSvv/7avn37+Pj4Ll26VPPmS5cuxcbGJiVdjov7QasFAKUSffpg7Fj4+qLKZY2xCleuXBk1atSdO3fc3Nx2797drl27at6cl5e3a9eu7dtj9+1bm5FhKw52747x4zF+PHr1kqTGjLGGg0wmJSVFRL6XXnrpzp07uuNarbb6D2ZkUHg4+fuTnR0B5X86dKCgIIqNpeJiOn+eVqwof3NODn3/vem+CdYIJCYmOjo6Aujfv39mZqbuuP6/n6iwkGJjKSiInJwqWtoLL1BwMCUkkFbLLY0xWTBVILxy5Yq4K+/du3daWpru+IEDBzw8PPSPVKOggLZvp5kzydm54jplb0/z5xNAP/xARHTjBvXta6JvgjUCsbGxtra2APz8/B4/fqw7vm3bNjs7u927d9fmJCUltGdP5YjYrh399a/c0hh79plklubZs5g9e93du3dfeeWVAwcOODk5iePR0dGjR48+c+bMunXranMeW1v4+eG773D/Pk6fRmgoPDzw8CHatMHo0Vi0CDk5pqg+azQiIvDFF66FhUVvv/12dHS0jY2NOL5y5cqJEyfm5eUdOXKkNuexsMDIkVizBqmpSEhASAg6dcLdu2jShFsaYzJg9NB6+DA1a0YODtpp0/5X/w49IiJCrVYDmDNnjkajqff5796lDRto7lxavJhmzeL7dPn63/8lhYIAWr78vO6gVqsNDQ0FoFAoQkND631yrZaSkmjdOm5pjD37jBwId+wgGxsCKCCASkoqjq9cuVKsEQwJCTG8FBEIi4vJzY3Cw/nyJDtaLc2bRwCpVLR6dcXxsrKyoKAgACqVau3atYYXxC2NMTlQG7FzuXkz3noLpaV45x2sXg3d4vilS5cuXLhQoVAsX778gw8+MFZxlpZYswbjxsHFxVinZI1ASQmmTcPmzbCyQkQEJk0qP66/lD4qKmrMmDHGKpFbGmPPNqM9I/z2W0yditJShITgm2/Ko6BWi4ULHy1btk6tVkdERBgxCgre3hg3zrinZA1afj58fbF5M5o2RVxcRRTMycGoURa3bnVycHDYu3evEaOg4O2N8eOhUBj3rIyxhsEo/colS0ihIIWCli6tOFhaSlOnEkDu7o9iY2ONUlAl/+//kZ0drV9vinM3emlpaWfPnt2zZ09xcbG562IcDx5Q794EUOvWdPZsxfE7d8jNjQDq1k1z5cp/TVH0pk00cCCFh5vi3IwxMzO0R0iEkBAsXAilEt9+iwULyo8XF2PyZGzciCZNsGyZna+vr6ER+0maNkVeHo4dM8W5G7fTp09369Zt1qxZr776qoODg6+v73/+859GnXX61i0MHIizZ9GxIxIS8PLL5cevXIG3Ny5fhpsbfvlF2bVrdUkb6i0jAwkJqN0UVMZYY2NIFNVoaPZsAsjSkjZvrjiel0fDhxNADg504oShsboaSUkEkKurCYtojPbt29e0aVMA3bt379mzp+53rVarhwwZsnz58uvXr5u7jnVTWkpduhBAfftSenrF8cREcnQkgPr3p5pWzxskMZEA6t7dhEUwxszFoECo1dK775KtLe3aVXEwO5v69y8fv7pwwdD6Va+khGxtSaGgjAzTFtSIxMbGWltbA5g6dWppaSkR3b59e82aNT4+PpaWlrqg2LFjx+Dg4Pj4ePGehqZqSpdDh2jsWMrLq3hPbCzZ2hJAfn6kt07HJIqKyMqKlEp6+NC0BTHGpFe3QFj18qTR0KVLf3rPwIEEUKdOdPOmkeqop6SE4uLoiy8qjgweTACZ5hFk4/PDDz+IxZp/+9vfqi7WzM/Pj42NDQwMtLe310VER0fHwMDAqKioR48emaXOT7RhQw0pXbTa8l/97NlUVmb8Cmg0dPbsn8Y5+vUjgA4cMH5ZzwCtVhsfH3/u3Ll0/Q47Y41E3QJhjZcnIjp3jgYMoN9/N0btqigqImtrUigoK6v8yEcfEUDGWJ3Y6K1ataqWizXLysoSEhJCQkJeeuklXUS0sbEZPnx4WFjYvXv3pKlwNTZsoNGjqV07ys5+akvLzqaVK01VgYICUqvJwqKir/n++wTQZ5+ZqsTGKyMjQ0wCcHR0dHJy2qU/QMRYY1DnQFjj5cnUvL0JqBiM3bWLAPL2NkNNGpQlS5YAUCgUX+j3l2vhxo0bYWFhXl5eul2RlUqlh4dHaGjo6dOnTVTbGjWE5EG9ehFAx46Vv/zxRwLo9dfNUJOG7MCBA8899xwAe3t7sfWVUqlcsGBBiX5CDcYatjoHQrNfnhYuJIA++qj8ZW4uqVRkZUVFRWaoTEOg1Wrnzp0LQKVSrVu3rt7nSU9PDw8P9/f3FxNthBdeeCEoKCg2Nlbi61pDSOkSFCTyt5W//O03AqhNGzPUpGEqLS0NDQ1VqVQAPD09b968qdVqw8LCxKNoDw+Pa9eumbuOjNVKfQKheS9PcXEE0KBBFUd69PjTnbuslJWVzZgxA4ClpeWWLVuMcs78/Pxt27ZNnz69VatWuojYokWL7Oxso5y/NkRLI6KEBHJ0NE9LW7eOAJo8ufylVkstWxJAeluKyVdKSoq3t7fo/wUHB+vfJ508ebJTp04A7Ozs1qxZY8ZKMlZL9QmEZNbLU3Y2KZVkY0O6ZeLvvksALVtmhsqYV3Fx8YQJEwA0adJkz549Rj+/RqM5ffp0aGho79697e3tO3ToINnafF1LI6I336Q+faQp9k+Sk8u3J9QZNYoAMtL9RiO2bdu2Fi1aAHj++ecPHz5c9Q0PHz588803xS3UxIkTpbyFYqwe6hkIjxyhnj3Jzc0kdaqRSCOSmFj+8ocfCKBx48xTmSeSIGDk5+ePHDkSgIODwzHTd4e7du0K4IRJl4U+yeTJBFBcnMTFEhFpNNSsGQH04EH5kdBQAmj+fDNUpoEoLCwMDg4WEW7s2LHVb30cFRUl5id36NDh6NGjklWSsbqqW2aZadPw5ZcAEBeHCxcwYUKdPm00Xl4AKhLKeHvn9+r19zNnXInIPBXSU1xcPHv27A8++MDFxeWdd96Ji4srKSkxeik5OTkjRozYu3dv69atDx06NGDAAKMXUckrr7wC4JjkWXxEnuukJImLBQClEn36AMDJk+VHvL2zBw2K/f33JWaozVOYonU9zZUrVzw9PVesWGFtbR0WFhYTE9OyZctq3u/v73/u3Lm//OUvKSkpgwcPXrx4sUajkay2jNVB/eLnzz8TQEOHGjUo19rmzfdeeWXTrFmf6o48//zzAP77X5Pkmay9vLy8YcOGAdCfb9KiRYupU6dGRkY+NNJi7NTUVJEv5oUXXvjtt9+Mcs4ahYeHAxgneb87JoYAGjFC4mLLffZZ2uDBx5Yt2yFeZmVlKRQKW1vbhpCFIDMzc/LkyZMmTXJzcwsJCUlISDBpceHh4ba2tgC6du167ty52n9QzKkRc5KHDBnyu4lWVjFmgHoGwowMUiioSRMyywXht99+A+Ds7Kw74u/vD8CQOZOGy87OFj0zZ2fns2fPiqdrHh4euoioUqm8vLyWLFliyGy6W7duiZkIrq6ud+/eNWL9q3f9+nUArVq10mq1khVKRPfvE0DNm5MBeznX388//wxg+PDhuiPih1+nSGAKt2/fFoPVIjgJnTt3DgkJSUpKMu7vKDc3NyAgQBQRGBiYn59fj5Ps37+/bdu2AOzt7SMjI41YPcYMV/8Uay+9RACZa6VZmzZtAOj6Q2FhYQBmzJhhntoQpaam9urVS/TSKsW5mzdvhoWFDR8+3MLCQnfN0t3F1+madfnyZbFmq0+fPhmSp5UTP3Pp58S3a0cAXbkicbFERKmpqQCaN2+uS9MzdepUAN9++60ZavOH5OTkdu3aAejWrdvt27cTEhKCg4PFb0do166dWPRieM81MTGxY8eOAJo1a7Zx40ZDTpWenu7j46MLqAUFBQbWjTFjqX8gnDGDAAoLM2Jl6mD8+PEANmzYIF6eOnUKwEsvvWSWyty6datz58419tKysrKioqICAwObN2+uu2Y5OTkFBgbGxsYW1bQQ8tSpU46OjgAGDx5srFHWOhEzVNdLvuvVxIkE0B+/aqm1b98ewOXLl8XLFStWAJg+fbp5akN04sQJ8WRu0KBBubm5uuMajUZERPGYQGjZsqVoXfWYvSUWBYq7t759+xolUbtWq12zZo3oxbq6uv7666+Gn5Mxw9U/EIpVVv7+RqxMHXz55ZcAZs+eLV6WlpY2bdpUoVBIn+rw8uXL4tLj4eFRy15aaWmpuGaJ+3rB1tbWx8dnzZo1qampVT9y8OBBOzs7AL6+vo9NnWH6KZYvXw5g1qxZEpf7+ecE0HvvSVxsOTHqrgv/SUlJokNvlsps377dxsYGwOuvv15NM0hOTg4NDdXPn+fg4ODv7x8eHl7Lgc0HDx6IOckKhSI4ONi4s6CTk5NFDhox6UbiwXbGqqp/ILx6lQBq29aIlamDxMTEStejoUOHAti+fbuU1dD10gYNGlS/XlpycvKSJUu8vLwUf2x/rlKpRIYzXS9Ef0MJM2auEjHAVfJdr44cIYB695a42HKff/45gHfffVe8LC4utrKyUiqV+r0xaaxfv15kVJ8xY0YtxzxFRNR/UG1jY+Pj4xMeHl5Nc92zZ4+zs7MYrvjll1+M9x1UePz4ce2XYTBmavUPhFottWpFAN26Zbzq1FpJSYmtra1Cocj6I/32xx9/DGDBggWS1eH+kSOil+bn51dYWGjg2VJSUlatWjVy5Ej9zZJcXV19fHzE5S84ONi89866n7nE3W5d/muzPFQ6cuQIgJdffll3xNPTE8C+ffukrIbIJYtaZFR/IvGgWv9+y9raWoxApKWl6d5WVFQUHBws3jN8+PD79+8b7zt4At3CfGdn5927d5u0LMaqYdB+hGPHVmxGIb1BgwYB0P3/2b17N4ABAwZIVHxcHNnYfDp48JQpU4zbSysoKIiNjQ0KChJ35ba2tmq1un6XP6MbPHgwgFjJd71ydyeAzLImu6CgQK1Wq9Vq3aCi6Mp8+umn1X/QaMrKtv3rX2Ko4D//+Y+BJ0tJSRFTt8TdFf6YzBwWFnbkyJGXX34ZgFqtDg0NrbqNlyncuXNHLFE1xRgsY7VkUCBctowA+mPQSGrnz5/Xn5ny8OFDlUplZWVleOesZhERpFYTQHPmaE12vSgtLd2/f3+zZs0ApKSkmKiUOvnoo4/q3SkxxDvviJlZ5lm9t379+oSEBN3tzsaNG8UwgBRlFxbS+PFaheKdIUNiYmKMeOL79++vXr1aPyKKOVydO3eWeNeRsrKyxYsXq1QqGxsbydbFMqbPoEB4/HjxX/6SEhAQY6TKGEqs8ZozZ45p/zt9/TUplZLtgih2ejNw5rqx7Nq1C4C35Ltebdp0rG1bz0mTJklc7hOJJZUODg76g4omkZdHI0YQQPb2dOSIiQrJzMz8/vvvhwwZYmtra29vb679mUU/+/333zdL6UzmDAqERUVF1tbWSqWyISTV3bhxo1qtdnJyEre3HTt2DA4Ojo+PN3ISkCVLCCCFgj7/3Jinra7AJQD++te/SlNc9XJzc6Xrduu5dOkSgA4dOkhZ6NOEh4fb2Ng0a9ZMqVSKDAlGWVpQWWpq+YhwmzZk+mUGGo1GrEQ015KGgQMHAvjpp5/MUjqTOYMCIRF5eXkBMNHUstpbvXq1yOE0fvz4KVOmODg46C/UmzFjRkxMjKELeLVamjePAFKpaO1aI1W8ZgkJCQB69eolWYnVExPfJUjzrU+j0YiBO1NP36jRkiVLxFwSV1dXKysr0cYUCoWnp+fSpUuNFhFv3KBOnQigrl1JqlHx9957D8DHH38sTXH6Hjx4IG6wzLJAljFDA2FISAiA//mf/zFKbepHtzn7sj+2YiotLT1w4MA//vEPkRRDf+L4/ogIqsfFtKyMZs4kgCwtJd6GR9ftln6+/hO9++67AJZJvuuVSOIq8fIYfVqtdv78+aKlffHFF/THtKbAwEAxeVhwc3MLDQ016DHbqVPk5EQA9e1LEk7Q3bdvH8yUleKbb74B4OvrK33RjJHhgTA2NhbA4MGDjVKbutJqtf/85z+rn1B348YN/WlyZwcNIoDc3Cg0lE6fptosSCguLs9u0qQJmWOSt0hh2kDml//www9i7ZfE5S5atAjAokWLJC5XKC4ufuONNwBYWlpWHb4rKCiIjo6eOnWqfs6gtDFj6MMP65yEcN8+srMjgIYPJ2kf15WVlYmtmJOTk6Usl4hGjBgBcyQtYkwwNBBmZmYqFAobGxvp5z2XlZXNmjVLXJuioqJqfP+DBw/Wrl2bN20a2dgQUP7HxYX+/nfav7+69OFpadS5Mzk4kOS78QmiI2KWMauqbt68CcDR0VHiRY0i//WwYcOkLFTIz88fNWoUgKZNm1a/AXJRUdHOnTtnzpzp3aULqVTlbax9ewoOpvj4mlPUR0aShQUBNG2aWfLZi/9QixcvlrLQnJwcS0tLtVotffpcxgRDAyH9sWXryZMnDT9V7RUXF4vcV7a2tnV+Qvn4McXHU3AwtW1bERGbNCEfHwoPp5wcSkujzz6jmTNp6VISOS9u36YLF0zxjdSGiAFDzbXrVRUipdzVq1elLDQtLQ2AnZ2dNOvbdDIzM/v37y8WfZ85c6a2Hysro4QECg6m556raGOOjhQYSLGx9LS7xgsXyN6egoPNs9cG0S+//AKgR48eUha6fv16ACPMtdUWY0YJhDNnzgTw1VdfGX6qWiooKHj11VcB2NvbG7TztUZDp09TaCh5eFRcrWbPpk6dKCCANm2iCRPIxYXMlNtTJyMjQ6FQNGnSpCFsg0fm2/SqUv5rCdy8ebNLly4AXFxc6rnthkZDx47RvHnk4lLRxlq2pOnTKSXlCbdcEm6tVVVJSYlI9XJFws0+xAKhb775RrISGaukbjvUP5GYOCrZ3uW5ubkjR44U6RAPHjwoSq8npRIeHli8GKdP4/p1LF+OwYPRqhVcXPDTT3jjDWzdirZtceCA8apfH46Ojl26dCkoKPj111/NWxNB4t+4jshtliTVdvXJyckDBw68du2ah4dHYmKi2GCkzpRKDBiAL77AzZtITkZoKNzckJWFH35AaSm8vHD+PIYNw8mT6NsXhYXQ2ztCehYWFn5+fgC2bdsmTYn5+fnx8fFKpXLs2LHSlMhYVUYIhN7e3gCOHj1q+KlqlJaWNnjw4GPHjnXo0CEhIcHd3d1op37xRXzwAQ4exP378PWtOP7aazh40Gil1Je5Ys8TySEQHjp0yNvb+969e0OHDj1w4IBufapBunXD4sW4dAmXL+P773HqVEO75QIgNtuKjo6WpridO3cWFRV5eXnp76fImMSMEAg7d+7cunXrBw8enD9/3vCzVSMlJWXgwIHnz5/v2rVrQkJCPe/Qa/TCC8jKqnh54wZ69jRJQXXRoAKhu7t706ZNr127lp6eLmW5ffv2BbBt27bo6OjHjx+brqCYmJjRo0c/fPhw/PjxO3fuFFnujMnVFYGB2LOnAd5yjRgx4tWuXRc3a6a5eVOC4kTEFduLMmYuRgiEROTg4NCuXTsPDw9vb++lS5deu3bN8NNWcvXqVW9v799++83Dw+PIkSP6O/kZ2bhxiIzE9esAkJ6OPXswZoypyqq1BhUI1Wp1v379iOjEiRNSlrtv3z5ra+v09PSJEye2bNnS19c3IiLi4cOHxi3l66+/njhxYlFR0fvvv79lyxaxAZZJNMhbLisrq90eHr6HDqlM3yksKioSufJ5XJSZmeGPGVeuXAnA2tpal70XQK9evT7++ONTp04ZZZL96dOnxQqnV155RYrcE/Pnk1pNfftSu3b03XcmL64WtFqtGJ27ZZZdr6oQm17Nnz9fshKTkpIsLCyUSuXs2bP79++v207I0tLytdde++677wyffK/VakNDQwEoFIrQ0FBj1Lpa589T164k8uKmpdFzz1ED2ZZv2zYCyNPT1OWI6dB9+vQxdUGMVc/QQHj58mWxZXZUVFR2dnZUVFRgYKD+suJWrVoFBgZGRUXVcmvsqg4dOiTGpsaMGSPd5uwFBXTtmlnWcj2NuGv+wVy7Xv2ZmEwh2aZXubm5Li4uABYuXCiOpKenh4eH+/j4WFhYiJam207o999/r0cRZWVlQUFBMNJuR7XV8G65iIgKC8nOjhQKun3bpOW89dZbAP7v//7PpKUwViODAmFJSYl4bDNjxgz942VlZQkJCcHBwfoDmCLD2Zo1a+qULjIuLk4EWqNv+9foLFu2DHpbpZtRRkZG7969ra2tJUt2M2XKFNF1qJq3ITMzU0RE/Q2NRZ6z//73v7U8f1FRkZgkYmtru2PHDmNXv1oN75aLiGjSJAJo+XLTlWCWpRqMPZFBgXDhwoUAXFxcqhmuTE5OXrJkif7W2Eql0sPDozb5GDdu3Cju99977z2Jl1E3QOIBocSLnau6f/++yLvt4uJikl0Xqli3bh2Apk2bVh/YcnJyxIBE06ZNK0XE6ltadna2mPns4OBg0LLUZ0lkJAHk5WW6EsTTwe7du5uuCMZqqf6B8OjRoyqVSqlUHjp0qDbvT0tLE3fuupz90NssqWpvT7ehRAPZnN3sGsKmVzdv3nzxxRdFgKnfCGRd/fbbbyKldaUx4RMnTqQ/JSH148ePRS5s/SF60dISEhIqPbS+d+9ez549AbRt2/aC+ZIHNTgFBdSkCSmVdO+eiUoQA9FSPItlrCb1DIR5eXliF9yPPvpI//jx48eriW2CyNkfFBTUunVr3XWqRYsW/v7+4eHhonNZdUMJRkSi47Jr1y6zlJ6cnNy2bVsAffv2zZRkWkdJSUm/fv0ABAYG6h/XarUdOnSocTvAoqKi+Pj44OBgZ2dnXUvr0KGDiIgajeby5csiW42bm9udO3ck+I4ak3HjCKCVK01xbo1GI/77nz9/3hTnZ6xO6hkIp02bBqB3796Vntn861//0l1xHB0dp02btnXr1ry8vCeepKys7NChQ/PmzRMxVbC2thZJrVQq1XcNZ/pAwyA2vap08yGNkydPtmzZEsCQIUMk28R87ty5ojNXaew9Ozv7tddeq/12gKWlpfHx8e+9957+qm1nZ2cxiDpw4MCGsLN0g/PjjwSQaTaWOXTokBhdN8XJGaur+gTCmJgYEbEuXrxY6UtlZWWnT58ODQ11c3PTj23Dhw8PCwu7+/Q8ivqbJdnZ2dnb29dmQwm5MdemVwcOHBDjk35+fpLtTb9vX2rTpnYWFhZJSUlPfEM9tgPUaDSifYpsDB06dBg0aJB0U5Ebl0ePyNqalEpKTTX6uYODgwEsWLDA6GdmrB7qHAjv36fBg/c6O7dZWdOYSaWNAHXXqZCQkKqPanQuXLgg7tbrWjE5yMrKUiqVEm969fPPP4tF5W+++aZkWb8fPCBnZ3Jxufrvf0fU+OYnPhR0cXF54kNBnaFDhwLg+63q+PoSQN9+a9yzarVaMSKdmJho3DMzVj91C4RaLY0aRQCNH59d+5XymZmZYjqffqqq9u3bBwUFxcbGFhUV6b85IyNDDKvWqWLy4erqCuBpnSSji4iIEPcxc+bMkWzirlZLPj4E0KBBVFZWhw/q1u3oP35u3769eGhdKYqL9RJbt241cu2fJQcOUEQE5eYa96wiW+xzzz0n8ZaWjD1N3QLhihXlu6rVZSlghcePH+/YsSMoKEj/UU2zZs30Z47l5uYCaN68eX0KkAGxdepyUy7w0lm5cqVZJu5++WV5M6v3jMWysrL9+/f/7W9/E7N7hDZt2uTqXdMnTpwIYMuWLcap9DOp6i5RxrBgwQIAwcHBxjohYwaqQ67RK1cQEgIAq1ejfpnibWxsxowZI9bUJycnh4aGenh45OXliSXzguh/lJWV1acAGRBJR3fu3GnqH9HSpUvff/99Ilq+fLmYxCuNixfx0UdQKLBuHfSiWN2oVKqhQ4euWrXq3r17oqV16dKlbdu2+mOnrAa5uU/YJcoYRGY1TrTNGg51zW8BAJSVYdo0FBZi+nT4+xuh4G7dunXr1m3x4sW3bt3SXwEtVtCXlpYaoYxn0cCBA0eOHLlv3z4nJ6fhw4f7+PiMHTvWuNd3Ipo7d25YWJhKpVq7du306dONePLqFRRg0iQUFSE4GH5+xjmnrqVlZmbqHxcZHojIOMU8e3bvLt8lCsAbb8DbGwcO4PffsX9/9Z8LUCi0T/+pPnz48Nq1a46OjgbtJMqYcdWy5/jhhwSQiwuZOuW1RqMBoFQqTVtMY3bmzJmuXbvqfoOWlpYjR45ctWpVSkqK4ScvKysTkc/Kyio6OtrwE9bJ9OkEUPfuJMFETn9/f/BkmWq8/TatWFHx8tNPad48mj2bgOr/KP/IIfU0o0aNioioeQ4UY5KpVY/w2DEsWwalEuvXw+hbs1WiVCqVSqVWq9VqteIBFaukd+/eV65cuXXrVmxs7I4dOw4fPrx37969e/fOmTPHzc3N19fXx8dHP6dd7RUXF0+ZMmXbtm1NmjSJiYkZMWKEKer/NFu3Yv16WFtj0yboDZabCvcIa1B1l6hBg9CtG2pqFZFA9T9TFxeXPn36GKGGjBlLjaEyL486dSKAJFvGLRZKV5pNyp4mKyur6qYfTk5OgYGBVSflViMvL09EPgcHhxMnTpi0zlXduUMtWhBAa9ZIVOKkSZMAREZGSlReo9Ngd4lizNhqDoTTphFAvXuTZEvXmjRpAuBp+WjY05SWllbd9MPW1lZs+pFa7bLo7Ozs/v37A2jdurVZUm6GhYllOdKVOHnyZACbN2+WrshGp2HuEsWYsSmo2qEhIixejK++QlISXF1N0iWtysHBITc3Nycnx97eXqIinzmXLl3asWNHXFzc8ePHxa9YpVK5u7v7+PhMnjzZ9c+/y9TU1FdfffXixYsuLi579+7Vz3gnpbg4eHmhRQuJinvjjTc2b978008/BQQESFRkY/T4Me7dg4sL1LWdWMdYo1PxEO7CBaxcWf7v3FysX49Ll7BzJz75BLduoVUrfP+9RHVq335gu3aDSku1EpX3LOrWrVtISMjRo0dv3769Zs0aHx8ftVp95syZTz75xM3N7cUXX/z73/++b9++srKyW7duDRw48OLFi25ubgkJCaaOgk9sZjt2AICvL7Ra6ZqZnd2L7doNsrBwkqi8RsrWFp07cxRkzzhd33DDBgJI7HVz4wb17UsbN9Lbb5d/9fJlcneXqJfapg0B9Vyzz57m4cOHkZGRU6dOdXBw0P32HR0dRbqf/v37S5N4uuE0szfeIIA2bZKoOMZYg/WnaZmjR2PRIuTkSB2MK7GwAABeSWhczZo1mzRp0o8//piRkZGQkBASEtK1a9fMzMyhQ4cOGTJk7969+gHSpBpIM2OMMeFPIx6urvD0xIIF+PDD8iPR0ThxAgBKSiBZUg4xDMO5ZUxEpVJ5e3t7e3svWbLk0qVLzz//vK2trchjII0G0szE6hJePcEYq7xQ78MPcfw4jh4tfzluHJKTkZyM7dulqxP3CCXTrVu35s2bSxkFhYbQzBhjTKj8DNzSEmvWYNw4uLgAgFJZ3j+T8mE59wifeQ2hmXGPkDEmPCF1i7c3xo8vv0yIv7OzJa0T9wjloFIzy8/H9Ok4e9bc1WKMyU/FHfi0aRVH16wp/8eUKdixA9Om4euvce6cVHXiHuGz62nN7PPPsWED9u/Hr79KVBPuETLGhJqTed67h+xszJ6Nq1clqA/APUJZ+uADDBiAu3cxa5a5q8IYk5maA+E77yAwEPn5mDQJjx9LUKXyHiEHQllRq7F5M1q0QEwMvvlGihK5R8gYE2q1vcPq1XB1xcWL+Mc/TFwdAH/0CHloVG7atSsfLJ07F+fPm7s2jDHZqFUgbNoUUVGwtcXatYiIMHWVeGhUviZOxKxZKCrClCkmH37gHiFjTKgh6ba+deswaxaaNMHJk3BzM7TgCxcQH495857wpWPHUFgId3c4OhpaCmt0iorg6YkLFxAUVDGbpn7y8rBzJ+Li8P33sLKq/NVLl5Caim7d0KaNQaUwxhq7OgRCAG+/jfBwdOuGkydha1uf8lJS8PPPiIgonyh/4QJ69Kj46oULOHwY778PALm5iInB9On1KYU1apcuoW9fFBZi0ya88UadP56bi/h4xMVh2zYUFABAbCx8fSvewM2MMfYndcpMmp9Pbm4E0JQpdUtpeucOff45eXgQUP6nVSt67z26fv1Pb6uakZnJ09dfE0D29nT7traWH8nMzAwPL3rtNbK0LG9jSiW98gqtWEEPHvzpndzMGGP6avWMUKdJE2zZgiZNsH07Nm8+WeP7s7KyIiIiRowY4eeXO38+zpxB8+YIDERsLO7dw+rVePHFyh/hjMwMwF//inHj0KPHlcDAUaXVPi4WbczX17dt27Zr197ZtQsaDby8EBaGu3fLe37OzpU/xc2MMaZTt6FRYdOmx6GhU+7fj09KSurevXvVN+Tk5Gzbtm3z5s0HDx7UaDQAhg6NaNUqMCAAo0c/4WmNRoPDh7FjB3r2xMWLaNYMv/+ODz9EQABO1hxt2bMpO/uhu3uPu3fvLlq06NNPP6301fv378fExGzdujUhIUG0MQsLi7feWtW/f9Drrz/56XJODmJjsWsXRo7E5cvczBhjf6hfR3L69OkAunTp8ujRI/3jhw8f9vX1tbS0FCe3srLy8/PbtGlTfn7+E8+TnEwhIeUbEAK0cCHNnUvFxeTmRuHhPGYld0eOHFGpVEqlMj4+Xv/4ihUrlEqlro35+Phs2LDhafspZmVReDj5+FQMmc6dy82MMVahnoGwsLDQ3d0dQEBAgP7xqKgoACqVysvLKywsLCMj44kfP3uWFiygDh0qHhl26UIff0xffklz5xIRJSSQoyNfoRgtXrwYgLOz8wO9B32JiYnW1tY+Pj7h4eG5ublP/GBqKq1eTcOGkUpV3sbUahoxgr79llat4mbGGKtQz0BIRNeuXRObm69du1Z3sKCgYPXq1enp6U/8yM2btGQJde1aEf+ef56CgykhgbRaIqING8qvUEQUFET9+tW7duwZodFohg4dCmD06NFabfnEGa1W+7Qxhnv3clesoFdeIaWyvI1ZWtJrr9G6dZSZWf4ebmaMMX31eUaoExkZGRAQYG1tffz48Zdffvlpb7t9+3ZkZOT27fvPndtbVAQAzs7w90dAAAYMKF/XzNjT3Lt3z93dPTMz86uvvvrHU5Ib3blzJyYmZsuWLWfPXlars/LyFNbWGD4cvr4YNw6tWklbY8ZYo2JQIAQQFBS0du3aTp06nTlzRnQQdVJTU6OioiIjIxMTE0Upw4bdfuGFDgEBGDIEKpVB9WaysmvXLh8fH7VanZCQ4OnpqTt+48aNrVu3RkdHnz59WrQxW1vbmTMvDhjQccwY2NmZr8aMscbD0EBYVFTk5eV19uzZSZMmRUZGAsjJyYmLi9uyZcvu3bvLysoA2NjYDBs27K233ho7dqxuHg1jdTJnzpyvv/5a3HJlZWXFxsZu2bLl+PHjogGLNubv7z9u3Dg7DoCMsbowNBACuH79ep8+fR4+fPj222+np6fHx8eLhV/W1tajR48OCAjw9fW1sbExRm2ZfBUVFfXv3//8+fMtWrTI/mOraHt7e19f3wkTJrz66qvW1tbmrSFjrJEyQiAEEBUVJR4WFhYWqlSq/v37v/XWW5MnT27evLnhJ2dMuHLlyrx5806dOqXRaHx8fPz9/UeOHGlVdV0qY4zVhXECIYCLFy/u3bu3adOmEydObNmypVHOyVhVV65c6dy5s1rsWskYYwYzWiBkjDHGGqO65RpljDHGnjEcCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmv/H9F7PUBMswf1AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Draw.MolsToGridImage(mg_sim_mols[1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### MHFP6" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAzOklEQVR4nO3deVyU1f4H8M8wLM4IKiCKCiogoqDiWiKgqCxlgP1QTG+CmkWWhmbZkObFbnUdrVvkUhdbFJGrQYoBbow74I6WCojK4hKKIpCA7HN+fzzTgIDKMvMMON/3y1cv5pmH53yHYL5zznPO9wgYYyCEEEK0lY6mAyCEEEI0iRIhIYQQrUaJkBBCiFajREgIIUSrUSIkhBCi1SgREkII0WqUCAkhhGg1SoSEEEK0GiVCQgghWo0SISGEEK1GiZAQQohWo0RICCFEq1EiJIQQotUoERJCCNFqlAgJIYRoNUqEhBBCtBolQkIIIVqNEiEhhBCtRomQEEKIVqNESAghRKtRIiSEEKLVKBESQgjRapQICSGEaDVKhIQQQrQaJUJCCCFajRIhIYQQrUaJkBBCiFajREgIIUSrUSIkhBCi1SgREkII0WqUCAkhhGg1SoSEEEK0mmoSYUYGtmxRfF1UhF27ACAzE4cOKQ4WFCAmRiVNEUIIIaqkmkR46hTmzUNsLADk5WH9egA4dw7btytOuHkTYWEqaYoQQghRJZUNjXp6YulSlJSo6nqEEEIIH3RVdSFHR4wYgZUr8dZbdQd37cLZswBQXg4zM1U1RQghhKiMKifLhIZi3z6kptYd8fHByZM4eRKbN6uwHUIIIURlVJkIRSJs3IilS+uO6OlBLIZYDJFIhe0QQgghKqOaRCgQKL5wd8eUKXUPlV9w7t3D8uVYtQqPHqmkWUIIIaStBIwx9V399GkEB8PWFtu2AcCbb2LZMiQmwswMM2eqr1lCCCGkudS7oL5rV5w5gyNHFA83bYKhIfbuxZgxam2WEEIIaS719ggZQ8+euH8fOTno3x8pKdiwAVIp+vVTX5uEEEJIC6i3RygQwMkJAFJSAGDVKgAIDa2rOEMIIYRolsrWET6JszPi4pCSgtdfh0ym7tYIIYSQllF70W1nZ+DvHiEhhBDS3qj3HiGAykp064aqKjx4gG7d1NoUIYQQ0mJq7xEaGGDUKMjlOHVK3U0RQgghLcbHfoQ0OkoIIaTdokRICCFEq6n9HiGABw9gZoZOnVBcDH19dbdGCCGEtAAfPUJTU7z66s9du47+44+zPDRHCCGENB8fiRCAicmJu3dTU2h4lBBCSDvDUyJ0dnYGQImQNJNcLi8uLvbz81u2bNkPP/xw9OjRO3fuaDooQsjzSe2VZTguLi4AkpOT+WmOdHQxMTHvvPNOUVFR/YNdunSxtbUdOHCgnZ3dQAU7IyNDTQVJCHk+8DFZhtOrV6+7d+9mZWVZW1vz0yLpoORy+fDhwy9duvTOO+9YWlpevXo1MzPz6tWrDx48aHDmhAknLl50sraGvT0cHGBtDe7r+htBZ2Tg9GnMnQsARUU4cgR+fvy9FkJI+8dTjxDA2LFjd+/enZKSQomQPF1MTMylS5f69esXFhamX2+ecWFhoTIpcmpr+xUVITUVqal1366ri/79MXAg7OwwcCD09PDmm+jaFf/3f8jLw/r1lAgJIY/hLxE6OztziTAgIIC3RkmHI5fLv/jiCwCffPKJ/uOrbUxMTMaOHTt27Nj6B4uKkJaG9HRkZyM7G2lpyMzE9eu4fh1798LEBF99BU9PLF0Kd3deXwghpKPgNRGC5suQZ+G6g3379g0MDGzO+cbGcHGBi0vdkaoqZGUhMxNXr6KmBgAcHTFiBFauxFtvqSdoQkhHxl8iHD16tFgsTktLKywsNDEx4a1d0oEwxhp3B6Oiot577z1ra2t7e3sHBwdra2tra2sHB4dOnTo1eRF9fQwejMGDFQ83bwaA0FAMH46RI3l4ER1TYSGiopCVBTs7BAaic2dNB0QIf/hLhHp6eqNGjUpKSjp16tSUKVN4a5d0IMru4Jw5c5QHr127VlRUlJqamlrvTqCenp6n56cCwcfcjUBbW9jZoXfvJ15ZJMLGjZg5E0OHqvUVdEx//QVnZwwaBG9v/PILvv8eFy5AKNR0WITwhL9ECMDOzi4pKemTTz4BMGnSpCd9oifaSdkdXLFiRf27g6tWrQoKCkpPT8/+W1paWmZmZmmp7bFjSEiou4KREWxt8dJL+OKLuoMCgeILFxdMmIDHV2QQAMCBA7C0RGwsAMyfD0dHnDgBV1dNh0UIT/hbPnHmzBlPT0+BQFBcXAxAJBJNnjzZx8fHx8enV69e/MRA2rOYmJgZM2b07dv32rVr+s8qSltZWXn9ekVmZterV8H9y8xEQQEAvPYaduxoeH5aGhwdYWeHtDT1RN+hvfEGRo7EokWKhytWQF8foaEajYkQ/vDUIzx69OjUqVMfPnzo5ubm4uKyb9++8+fPJyQkJCQkvPvuu05OTj4+PlOnTrWzs+MnHtLeyOVYt+6/AJYvX/7MLAjAwMDAwcHAweGxg0VFuHoVTQ402NgAwLVrqKqiyu+NWFqiuLjuYU4OLTEh2oWpX1xcHDcK+o9//KOqqoo7eOvWre++++6ll14yMDBQBjNw4ECJ5OPjx1lNDQ9xkXYkJoaJxWWvvrqhsrJSTU3Y2jKAXb6spst3ZBcuMAcHlp3NGGN5eax3b/bwoaZjIoQ/aq81GhUVNW3atIqKinfffTcyMlJPT487bmFh8c477+zbt6+wsDAuLi4oKMjc3Pzq1auJiafHj0ePHpgxA1u34uFDdQdINI8xfP45Hj0Sv/TSwuZ0B1uH6z7S0GhD588jNxdOTrCzg5MTXnwRK1bU3VklRBuoNc1u2LBBR0cHgEQieebJNTU1x44dW736OPfJnfvXqRN7+WX2/ffs9m21Rko0KSaGAczSklVUqLGVFSsYwEJD1dhEhzR/PgPYjz+ysjJ27RrbtYsBzNdX02ERwh819gjXrFmzaNEixthXX30llUq5g+Xl5Tk5OU2eLxQKx48fHxLievUqsrIQFgZ3d9TUYN8+vPMOLCzg4ICQECQnIz0dW7YovquoCLt2qe9FELXjuoMAli9HvWFy1bO3B6hH2FhpKQB07gyxGAMGoLxc8ZAQraGWRMgY+/DDD0NCQoRC4Q8//PDBBx9wx0tLS319fcePH5+bm/v0K1hbY/FiyGS4exfR0QgIQJcuSE/HmjXw88PJk5g3TzHZm6seSTquXbvwxx+wtMS8eeptiIZGm1ZWBtTLfA0eEqIFVD9rtLYWCxYIrlyx19fXj4qKmj59One8oKDgpZdeSk1NtbCwKOc+dTaDqSn8/eHvj/JyHD6MuDh07QodHaoe+ZxgDJ99BgAff6ze7iCAQYOgq4tr11BZqfa2OpImE6EhbW5FtIiKe4SVlfD3x48/4vLlNw4cOK3Mgnfu3Jk0aVJqaqqVldXRo0cHK+tfNZtIhFdeQXg41q4FAEdHvPYaVq5UbfiEb7Gxiu7gG2+ovS0DA1hbo6YG166pva2OpEEiVI6UEqI1VJkIy8rg64vYWBgbY88euLkN547n5OS4urpeunTJwcEhOTnZhlvS1Wahodi377H9d0jHwhj+/W8ACAnhqYtGo6NNoKFRovVUlgiLiuDpicREmJvjyBGMG6c4fvFiupOTU1ZW1rhx45KSkno/pRxkC3HVI5cuVdX1CN9270ZqKnr14qM7yOESYUaGnKf2OgRKhETrqeYe4d278PLCxYvo3x+JibC1VRw/cwb/+IedoeEL9valv/32m5GRkUqa49Y4yeUwNsagQbTkqaPauhUAVqxouhaMOgwfntK79weXLvUBdvLUZPtHiZBoPRUkwtxceHjg+nUMHozERFhYKI4fOIBp01BWJpw3b8d//6urwoXSc+cCQHU1xo9HebmiwiRRYoxt2LChuLj42rVrvr6+Xl5eqvoIoloxMdixA/7+/LVoa2uUl3c6PX0Qf022f9xNQeXsGEqERAu1cR1iejqzsGAAGz2a3b9fd3z3btapEwNYQACrrm5jI080fjwDWEKCuq7fEVVXVwcEBABQJj9dXV1nZ+ewsLCbN29qKqr0dLZ5s+LrwkK2cye7coUdPKg4cv8+i47mI4yKigpdXV1dXd0KtS7d70Bqa5lAwAQCJpcrjnh7M4DFxWk0LEJ41aZ7hI8eYfJk3L6NSZNw+DC6d1ccj4zE9OmoqMB77yEiArpqq+zt7AwAtOm9UmVl5cyZMyMjIzt37rxhw4awsDB3d3eBQJCSkrJkyZK+ffs6ODiEhIQkJyczvnYd4Zw61XDp57lz2L5d8ezNmwgL4yMMAwMDGxubmpqaq1ev8tFeuyd/9KimVy+5uXndDQbqERLt07JEmJHxWEmX/fuxYQP8/LBnD5Rjb+vXY84c1NRAIsG6deq9gUeJsL6ysjIfH5+dO3caGxvLZLLAwMDFixfLZLI7d+5ER0cHBAR06dIlPT19zZo1rq6u5ubmgYGB8fHxlZWV/ITHLf0sKeGntSdycHAAkEYzRwEA98vK9PLyetXWKo/M1def1LdvRrscSydETVqWCBt/rvfzw86ddZMdqqqwbRsEAmzYgL+rqqnS/v144QVIJIqH48ZBRwdnz6KqSvVtdSxFRUUeHh4ymczc3Pzo0aNOTk7Kp0xNTf39/bdu3Zqfny+TyYKDgy0sLO7duxcZGenr62tiYuLj47Np06b8/Hy1Rth46eeuXXB0hKMjZs5Ua8uPoURYX1lZGYDO9fp/Z2/dOnLzplws1lxQhPCtxUOjT/9cr6+PvXuxcycWLmxrZE3S1cXZszh2TPHQ2BiDB6O8HOfPq6W5juLu3btubm4nT57s379/UlLSsGHDmjytU6dO7u7u33777a1bty5fviyVSp2dncvLyxMSEt5+++0+ffq4uLisWbMmIyNDTXE2WPrp44OTJ3HyJDZvVlODTbC3tweQnp7OX5PtWONE2PgIIc+9FifCZ5Z0MTXFq6+2JaSnGTsWuro4fx6PHimO0OgoV6/g4sWL9vb2ycnJAwYMaM53OTg4SCSS5OTknJyc8PBwb29voVCYkpISEhJib29vY2OzePHigwcP1tTUqDDUBks/9fQgFkMshkikwkaegXqE9ZWWlgIwrFdQjRIh0UKtmSyjwZIuhoYYNgzV1Th7VnFEyxNhenq6q6vr9evXR48efezYsT59+rT0Cv369QsKCoqPj793796OHTtmzZplbGycnZ29bt06Dw8PCwuLN998c+/evW2MU3mr2N0dU6YoHnL/LStDcDB/G0/a2dnp6upmZWXxdnO0PaMeISFAC5dP/PwzW7aMMcZkMmZqytzc1DGR9Rnee48B7IsvFA+vX2cAMzOrm/6tPc6cOdO9e3cAbm5uD1W3pXhNTU1SUpJEIrGzs+N+SVxdXY8cOaKq6zfw+ecMYEZGLClJTS00NGjQIAB//PEHT+21Y7t37wbg+/fug7W1tQKBQEdHR66Ff05Ei7WsR1j/c72nJ/idga/QoAtoY4NevXD/Pq5f10AwGnT06NHJkycXFBT4+Pjs27dPhUvmhUKhi4uLVCq9cuXKxYsXg4KCkpKSFi9erKrrNxASgjlzUFICLy8cPqymRurU1tZyI4Hblas3NKWwEOvXY8kSfP+9YtEC7xr0/8rKyhhjYrFYQOWaiFZpXf5cuJDp6bHfflNpUm6e27cZwLp1Y7W1iiOBgRUvvHB3x47LGoimkXv37mVnZ2/atCkvL099rfz222+dOnUCMHv27Gr1FSxgjDFWUVHRqVMnHR2d4uJiNTVRW8vmzWMAE4uZTKamRphcLo+Li+NmEpmbmwPw9/e/d++eutp7uuJiNmgQe/VV9uOPzMODDR3Kamr4j2LTpk0A3nzzTe7hnTt3APTs2ZP/SAjRoFYuqDcxQXW1Zu7M9emDfv1gbFx15UoRd2T48O/OnDE/eDBMA9E87ubNm66urlOmTAkKCurdu7eDg8OqVatSU1OZSvvO27ZtmzZtWkVFxcKFCyMiInTVV7AAAGBgYDBy5Ei5XH7q1Ck1NaGjg59+wrvv4tEjeHsjIUH1Tezdu3fUqFG+vr4XL160sbGZNGmSWCyOiYmxt7ffytU85dmBA7C0RGws5s9HYiIEApw4wX8UjXuEeHzuDCHaoJWJkBufTE5WZSjNN2nSezk5BklJ0X8H4wwgRdMTZjIzM11cXDIzM6urq728vEQiUXp6+qeffjp69GgbG5slS5YcPny4urq6ja1s3Lhxzpw5NTU1Eolkw4YNOjoq3lGySS4uLlDzT5hbe7poESorMW0a4uJUduUTJ05MnDjxlVdeuXDhgoWFRVhYWFpaWlRU1OXLlz09PQsKCubMmePm5sZ3rZn9++HrW/fQ25uPceFGGiRCbhIpzZQhWqd1Hcm//mJCITMwYOXlKu2gNs/GjRsBBAQEcA+rqqq4uxoaG+ZiLDU1tUePHgBefPHFBw8eMMYePXrErV6vv/NU586dvb29IyIiioqKWtGKVCoFIBAI/vOf/6j4BTzVb7/9BmDSpEnqbkguZ4sXM4Dp67PY2LZe7fTp097e3txPvnv37lKp9NGjRw3OiY6ONjMzAyASiUJDQ6uqqtra6tNVV7Mff2RffMH++U/22Wd1x2fNYjEx7Pff1dt6Ix9//DGAkJCQ7777LiEhgZs74+TkxHMYhGhW64tuDxvGAJacrMJgmuv3338HYG1trTzi5uYGIE5DlYKPHz/etWtXAJMnTy4pKWnwbG1t7blz50JDQ0eNGqXMiEKh0NnZWSqVZmZmNqcJuVy+dOlS7ht/+uknNbyIp7l//75AIOjcubPa8wRjcjlbsoQBrE+fmt2797fuImlpaf7+/tyMDyMjI4lE8tdffz3p5MLCwqCgIO7/i6Oj45kzZ1ob+1PV1rL//Y/Z2jKAGRiw/fuZgwPLzmaMsbw81rs3i4hgAgELCHiser2a3bhxIzg4mBsLFYvFRkZGc+bM0dTfESGa0vpE+M47DGBr1qgwmOaqra3lEs+ff/7JHVmxYgUAiUTCfzB79uwRiUQAZs6c+cw8oVy9rqenp0yK9vb2EokkKSnpSXPWa2pq3njjDQD6+vq//vqrGl7Es3FLKc6ePctPc8uXV9jY+Ojp6cXExLToG3NycoKCgoRCIdf/lkgkhYWF9U+oecKclCNHjtja2gLQ1dUNDg5u/IGmTWQyNmIEAxjArKxYeDirqWHvv8/09NjYsaxvXxYZyTZsUOzY0rMn27FDla035eHDh59++in3dyQQCKZMmeLu7s79Qjo7O6elpak7gCalpaX973//00jTRJu1PhFu28YA9vcCJL55eXkBUL5Lciu+XVxceA5j+/btXEp7++23a5XTWJuhsLCQK4TNvRNxzMzMAgICoqOjS0tLlWdWVFRMmzaNe1s/cOCAGl5Es3CZOCwsjLcWV65cyfWAIyMjm3P+rVu3goODDQwMuE8MQUFBDSbu1tbWRkdH29raPmlN5KNHjyQSCZdErays9u9vZX/0MUlJit3CAGZpycLDH9uWrKyMXbtWN1/0+nXm7q44ecoUlpurggAaqaysDA8P79mzJ/db5+7urvx8ExcXZ2FhAUBPT08ikfC5WVVubi73CaZz5853797lrV1CWFsSYW4uA1j37ppZyf6vf/0LwJIlS7iHxcXFQqHQwMCAzz/d77//npurIpFIWr0AmVu9Hhwc3LdvX2VGFIlE3t7e4eHh169f9/DwAGBsbHzixAnVxt8iP/30EwB/f38+G+XuiQqFwi1btjzltPv370skEq5frqOj4+/vn5WVVf8EuVy+c+dOrrgagMDAwKdc7cKFC6NHj+bObNP6ipMn2eTJiqxmZsak0mbdUZfLWUQEMzVVrCaRSlW4rKKqqio8PFxZfsjJyenw4cMNzikuLg4ODuZ+sQcMGND4BJW7ffv2ggULuA+U+vr6ixYtus/j4DAhrC2JkDHFlrwZGaoKpgUOHToEYMyYMcojQ4cOBZCSksJTBP/+9xwXF4FAsEZ1o8Pnz59ftWrVyJEjlcuZufej3r17X76s4VWSmZmZAMzNzXlul8uFOjo6m5Ub+9bz8OFDqVTapUsXbnzP39+/8T1XmUw2ZswY7ufZt2/f8PDwZ668rK6uDgsL4+6cGRsbh4eHt+iDzoULF+IXLlSkQBMTJpWysrLmfztjjN29ywICFFcYMYKdO9eyb2+E6w0r69AOHTo0+qlbIScnJ3PVyQUCQUBAADf/S+UKCgoafIK5fv26Ohoi5OnalAhnzGAA+/FHVQXTAqWlpZ6enl8oK60xtmDBAgBr165Ve9tyOVu2jAHV+vr7IiLU0UJ+fn5ERIS3t7eRkZFQKDx//rw6WmkRuVzOTYvN5uZ38Gjt2rXcO/LGjRuVB0tLS6VSqbGxsXJ8r/FPiVs4oRx2lkqlLRowyMrK4rrjACZMmNCcaU0ZGRkBAQE6OjpmpqYPe/ZkEglr1fRghbg4ZmnJAKanly+Vlrd2irZMJnN0dOReyODBg6Ojo5uT16uqqqRSKTfUbG5u/vTE2VIlJSVSqVR5h9Lb25sq3hENalMi/PZbBrB339XEEopGIiMjAXh5eam3mZoa9tZbign+6p/R4OrqCmDPnj3qbqg5pk6dCqCZd+xU66uvvuLeMdevX8/d4uJKw3AzO44fP97g/IsXL/r7+3MnmJqaNrlwopmio6O5gq4ikUgqlT6pN5mTkzN37lzu/qJIJPrggw/+Usn4XlkZk0iYSORjY2NjYyNrYd0dmUymHObt169feHj4k+YKPcnVq1e5KdkAvL29b9682aJvb6ysrCwsLIz7UMV9guEqThCiQW1KhOfP3zA3HzVw4EBVRdMWP/zwg5GRkb6+vpWVVXBwsEwmU/1c/8pKRS9YLGZ796r44k3hlnktX76ch7aeieuZLViwQCOtf/PNNwKBQCAQmJiYcO+h48aNazztJT09XblwwtDQUCKRtL0y3IMHD56yvuLevXsSiYSreKenpxcUFHT79u02ttjAjdOnhwwZwn0UmD9/foN5sE1KSUlR9ob79OkTFhbW6tvncrk8PDycG3/u2rVrWFhYi+aFKXGfYHr16qX8BHP06NHWhUSIarUpEdbU1BgZGQkEgvz8fFUF1DpbtmzhKo3Vrz1tamoaGBgYExOjmp0ZysrYyy8zgHXtyttGCQkJCdy4HD/NPR1XWWbo0KGaCuD77783Njbu1KnTkCFDGo/vKacdAhCLxcHBwar9tdy7d2+/fv1Qb33FgwcPJBKJWCxW3uK6du2aClusjxuo5NJtz549I548Jq/C3nB9eXl53OxltHx9BXeH0tramvv2F154gZYqknalTYmQMTZ58mQAsW2vAtIGykpjEolEuXqdu9XP0dXVdXZ2DgsLa/3ATlERc3FhAOvRg/F4x66wsFBHR0ckElVWVvLW6JMoq2+raerEMxUVFXGfdRqkwNu3bysXTnB9MuUCU9UqKSlZsmQJl2t79OjBzaYRCATTp0/nZ+HdtWvXuL84AFOmTLlx40b9Z9XRG26gpesr5HJ5dHS0cj8vBweHZt6hJIRPbU2E//znPwF8+OGHKommFZRVx7766qsGT125cmXNmjUuLi7cOxcAA13d6vHj2apVLUtm+fls+HAGsL592ZUrqoy+GbiMfurUKZ7bbRJXdHQvL8PCjd2+fRtA79696x9MTEzk+km6urrz5s3LVc/au/ouXLgwbNgwMzMzgUBQfxEeP+RyeUREhKmpKdfxlUqlNTU16u4N11d/fYWtre1T1lfIZLKRI0dyf3r9+/dvxR1KQvjR1kR44MABaKg4Yf2qYz8+derqvXv3Nm/e7Ofn99GLLyqmpHN9u4AAFhfHnnnvhFvjbGfH2jxToBW4u1M8Fxd9EolEAmDFihUaaZ1bwtHgnnRJSUnPnj39/f2v8PgZRSaT4fHVOzyrP1BpaWnJLcIzMDB477337ty5w0MASUlJgwcPxhPWVyQlJY0fP14ZXnOWrBCiQW1NhA8fPhQKhXp6emUtXSnVNvWrjrWgCld5OZPJWHCwYgkk908sZt7eLDycKetZPHjA1q1jixez775jpaUsM5O99BLT0H3QiIgIAH5+fhppvYG4uDgAbm5uGmk9NTUVwIgRIxocb10F87bg7t1OmTKF53YbiI+P79u378SJE5ssI6Bu9ddX9OrVi/szPHXqlHLwlit03upVH4Twpq2JkDE2fPhwAI2nsKtPZWVlW6uOyeXs7Fn2ySfM0bEuI+rqsokT2d277WHHVKXr169zd6Q0GIPSgwcPNHjP8vjx49BEIT3G2KFDh7Zv3668J/fLL7+A9zo7Tfrrr7+KiorUN0nnmS5dujR27Fgu81lZWXFfmJiYrF69un6lQELaMxXsZsfdN0rma3PCsrIyHx+fnTt3GhsbJyYmenp6tuYqAgFGj8Znn+H335Gbi/BweHtDRwcPH+LYsfawY6qSjY1Nr1697t27x2VEzTIxMbGzsysvL+c2AOFZg83z+LRhw4ZZs2adO3dO45E00KVLl27duilLxvBvyJAhJ06cCA8PF4lEjDGRSCSRSLKyskJCQtrDz4eQ5lBBIuRzX9yioiIPD4/ExERzc/OjR4+OGzdOBRft1w9BQYiPx717iIhoJzum1se9TI3vPMzR4DbIGkw/tJP7UwgEgqCgIF9f39zcXKlUKpVKu3XrpumgCGkBFSRCrvrJiRMn5HJ526/2FHfv3nVzczt58mT//v2TkpKGDRum4ga6doWDAywtUVxcdzAnB38Xa9YUDeaexjQYDLd/ukbSD+3k/kzcn7+y4g8hHYgKEmGfPn169OhRUlLy2muvxcbGcm8ZKpebm+vq6nrx4sXBgwcnJSWpcSzo//4PO3YgJwcA7tzBsWPw8lJXW83D5R7eBp+fToPBaLBH2CAHt5+h0faDfiak41JBIpTJZPfv39fV1f3111/9/Py6d+/u4eHx7bff/vnnn22/OCc9Pd3FxeX69eujR48+fvw4t6RXXYYPh6cn7Ozg5ISxY7FmDepVq9GIESNGiMXiK1euFBQUaDYSALa2tubm5vn5+dnZ2Tw33d6GRulNvz76mZCOq62JsLi4eP78+YyxRYsWhYWFubu719TUHDx4cMmSJRYWFg4ODqtWreKK6ra6iXPnzk2YMOHPP/90c3M7dOgQVwFZvb7+GsXFiIxEdjZmz1Z7c8+ip6fHLVk7deqUpmMBgBdffBGa6BQKBJaDB882NR3Ic7ugRNgMGhy4JqSN2poIFyxYcOvWLWdnZ6lUunjxYplMdufOnYiICH9/fyMjo/T09E8//XT06NFWVlZvv/12fHx8VVVVi65/9OjRSZMmFRQU+Pj47N27l6v8ywexGAMG4O+SNBrHTc1tD7cJS0tLs7KyRo4cqaybxZv8/JkZGZFC4Vs8twtKhM1APxPScbUpEW7duvWXX34xNDTcsmWLsoxZ9+7dAwMDo6Oj8/PzZTJZcHBw7969b9y4sWnTJl9fXxMTEx8fn61btxbXn5DyBPHx8S+//HJJScnrr7++c+dObgNP7dRObhMWFBRMnDjx8uXL+fn5/E/Z5+4+a+SddtSoHa6ukcp3+c6dRw4ePLtLl14aCKW9okRIOrBWr0C8desWty3qli1bnn6mshD2qFGjlO0KhUKuH/mkyliRkZHchhILFy5s3bYvz5OioiIdHR0DA4NWb6bTdnl5eUOHDgVgZWWlkRXcgYEMYM/6dVO98nIGMAODuiMeHgxgravl8LzitscqKCjQdCCEtFgrE2FtbS23Xeerr75a//j58+dXr179lEr8OTk54eHh3t7e+vr6yqRobW0dHByclJSkTHj1N5RoXYTPH25HuhMnTmik9ezsbBsbGwD29vYq32+vmfz8GMB+/ZXvdgsKGMBMTOqOjBvHAJaczHck7RlXa40KqpGOqJWJcPXq1QB69OjRoMg9VwWb6zQEBQXFxcU9aXfcBw8ebNu2bcaMGfVv+/Xu3fvtt9+eP38+nrChhDZbsGABgLVr1/Lf9OXLl/v06QNgzJgxGvzI7+XFALZvH9/t3rjBAGZpWXdk2DAGsN9/5zuSdqu6upob5tF0IIS0RmsS4fnz5/X19QUCwZ49exo8JZPJ5s6da2ZmpsxtJiYms2fPjo6O/uuvv5q8Wk1NTVJSkkQisbW15b6lZ8+eQqHwmSOu2iYyMhLA1KlTeW73zJkz3KY/EydOVM0Wx63F7QjJY1FbhfR0BrBBg+qO2NgwgGmuwGe7w93y79Kli6YDIaQ1WpwIy8uZv/8n3K27J53zlJuCT98d9/z58z4+PgB8fX1bGthzj1u31717dz73NT18+LCRkRH3f0Tjo14jRjCApaby3e7Zswxgo0bVHTE3ZwDLy+M7knaLWzTcq1cvTQdCSGu0OBEuXswA5ucX28x9l7Kysrj1hdyWaRx7e3uJRJKUlNT4Pf3nn38GMG/evJYGpg248cnMzEx+mtu9eze35+3s2bPbw35yAwcygP+tkdnRowxg48fXHTEyYgB7whiHNrp69SqAAQMGaDoQQlqjZcsnDh7EunXQ1YVE8qpYLG7Ot1hbW3PrC/Py8iIiIqZNm2ZoaJienr5mzRpXV9e+ffu+++67ly5dUp7PzRStqalpUWBags/q25GRkdOnT6+oqFi0aFFERAT3/0WzNLV8gmtXuVKcMTx6pJlI2i1aO0E6tBYkwuJivPEGGMOnn+KFF1rcEre+8Ndffy0qKkpKSgoODra0tLx9+/b333+fm5urPI3rOHL33kkDvBW83rBhw9y5c2tqaiQSyfr167kZvBqnqUQ4YAA+/xyzZikeVlSgthadOrWfcguaR4mQdGgt+Ji/YAFu3YKzMySStjWpq+vi4uLi4hIWFnb+/Pn4+Hh3d/f6z4IS4RM4OzsbGBhs3749Pz/f39/fx8eHW8qpWmvWrAkJCREIBF9//fX777+v8uu3mqYS4cCBWLFC82G0Z5QISYfW3EQYGYlffoGhIbZsUdkHYYFAMGrUqPoTavB3j5CGRps0YsSIvXv3Tp48OSEhISEhQU9Pz9XV1dfX19fXV7k5eFswxj744INvvvlGKBT+8MMP8+bNa/s1VSg7GxUVqLcAVTMoETZGhUZJh9asIa/bt7F4MQCsXw9119WiodGnEAqFkyZNUhYlEAgEhw8fXrJkibW1tY2NzeLFi5OTk1u9K2Rtbe38+fO/+eYbAwOD6OjodpUFMzKwZQssLDBgAIqKsGuXJoMpLQXq3TIkoB4h6eCenQjlcgQEoKgIr76KuXPVHhBNlnmm/v37BwUFxcfHFxYWxsXFBQQEdOvWLTs7e926da6urj179gwMDIyJieE+pDdTZWXljBkzNm/e3Llz5/j4eD8/P/XF3wqnTmHePMTGAkBeHtavb3gCY0hIQBv2OGmBQYNw4wYOHOCjrY6CEiHp0J6dCL/8EkePokcPhIfzEA/1CFugc+fOXAXzgoICrijBwIEDCwoKIiMjZ8yY0aNHD25jyLy8vKdfp6yszMfHZ9euXcbGxjKZzMPDg5/4W8TTE0uXoqSk6Wc3b4aPD9zckJmpxhi4jqlQiL590bmzhjum7QolQtKhPSMRyuVITIRAgC1b0KMHHwFRj7AVhEKhi4uLVCrNzMy8fPmyVCp1dnauqKjgNoa0tLQcPXo0tzFk4+8tKipyd3eXyWTm5ubHjh1zcnLiP/7mcHTEa69h5cqmn+3SBWZmOH4cI0ZgzRqo6dfnmR1TrUWJkHRoz0iEOjpITMS+fXj5ZX7ioR5hWzk4OEgkkuTk5Pz8fG5jSJFIlJqaym0MyS3rPHjwIPcTvnPnzoQJE06dOmVlZZWUlMRtLtFuhYZi3z5w2XzhQuzbV/fU9OnIzERQECoqEBKCUaNw5oxaYnh6x1RrUSIkHVpdIuSGfTjcfITMTBw6BKEQXl4oKEBMDB8B6euPtLSsMjI6ykdjzzUzMzNuY8i7d+/u3Llzzpw53bt3z8nJWbdunYeHR69evfz8/EaPHn3p0qWhQ4empKTwv79gS4lE2LgRS5eisBDffYcpU/D667h3T/GssTHCw3HgAKytcfEixo3DRx9Vcm/QKpGfj5qaZ3RMtRYlQtKh1SXCxsM+585h+3bFszdvIiyMj4CEQv1bt/Tu3TPgozHtYGho6Ofnt2XLlvz8fK4GrL29/YMHD2JjY8Vi8ZgxY44cOdKrV7veY1YgUHzh7o4pU2BigrAwGBrif/+DnR02baqbJuPhgcuXERoKHR2cPfvFkCFD9u/f38bWS0qwZg0GDkRSEvB4x5RwKBGSDu2xodH2MOzDVSSlkVF10NHRGTVq1KpVq9LS0q5cufLll18eOXLkyJEj3OYS7dncuVi7VvH11q04cgSLF+PiRXh5obgYb7/92DQZkQirVuH06drycllubu7LL788e/bs+/fvt6Ld0lJ8/jn69kVICEpKcOuW4vpcx5Qo0TpC0qE9lggbD/vs2gVHRzg6YuZMngLiSlrSXBl1s7Oz+/DDDy0sLDrup3grK+zfj+joumkyq1ahqkrx7IgRwuTkpLCwMENDw6ioKDs7u02bNrFmL7CoqsKmTbC1xcqVKC6GszOOHcOcOYpnuY6psp9KDA0NjY2Nu3btqulACGkVZfntn39my5axR4/YwIEsIoK5ubFt21hgICsrY2VlLDmZjRvHRxXw7GwGMCsrPtoiz4fCQhYUxAQCBrBhw9jp0489m52d7eXlxf22T5gw4cqzdq+oqmIREczKigEMYGPHskOH1Bg8IUTjGs4abTDso6cHsRgiEUQinhIz1yOkoVHSfNw0mf37YWWFixcxbdrtZcskymkyVlZW+/fvj46ONjMzO3bs2Llz5550HbkcMTFwcMCcOcjJwZAhiI7GyZOYNImvV0II0YS6RNhgPgL3UCBAZiZefBEZGTwFxN0jpKFR0lKenrh8GR98AAuL9776au3QoUMP1Kv+4u/vn56e/uWXX77++utNfvvBgwdffHHchx9WX7uGQYMQEYE//oC/P1/RE0I0R8Ceddfkww/xn//A2hqpqejWTe0BPXiA7t1haoqCArW3RZ5Lf/zxx1tvvXX27FkA/v7+GzduNDMze8r5iYmJn3zyCXe+j8+P06bNnz2btlgiRIs8u8Ta6tVwckJ2NgID+ajlSJNlSBs5OjqeOHGCmyYTExPzlGkyJ0+enDx5speX19mzZ83MzKRSaXT063PmUBYkRLs8u0cI4OZNjByJBw8QFqbYhkJ9Hj1C584Qi6G6ldBES+Xk5CxYsCAxMRHAhAkTNm3aNHDgQO6py5cv/+tf/4qJiQFgamq6bNmy9957TywWazJcQoiGNCsRAkhIgK8vdHVx9CjGjVNjQFVVMDCAvj4qK9XYCtEeMTExCxcuvH//vkgk+uijj6ZPn7527dqoqCi5XG5oaLhw4cKPP/6Y5v0Tos2amwgBLF2Kb76BpSUuXEDbV2DfuoWcHIwf3/A4Y1i7FgIBPvqorU0Qwrl///77778fFRUFQEdHRy6Xi8XihQsXSiSS9l9MgBCibi1IhNXVmDgRKSl45RXEx7dyNXFhIRISEBmJw4fRpw9u3HjsOhkZOH1asethURGOHEE72xePdGD79++fNWuWXC739vb+8ssve/furemICCHtQrN2qOfo6WH7dnTvjj178PXXLWumuBibN8PLCz16YM4cHDyITp3g5ISHDx87jba5Ierz0ksviUSihw8fUhYkhNTXgkQIwNISmzdDIMDKlTh9OveZ55eVlW3fvn3atNcGDGBvvIHERAiF8PFBVBTy8/HLL2h8a6Y91DslzysqiUkIaawFQ6NKoaFlsbEfFhfvuXDhQpO3WCorKxMTE2NiYmJjY7m3HheX64zZ+PvjH/9A4zVdjOHECeTlobQUGRnQ0UFFBd56C4sW4ciRVr0sQpqiq6tbW1tbXV3N7f9MCCEAWvN2sHKlweHDly9duhUQEJCQkKCjU9etTElJ+fHHH3fv3l1cXAxAIBA4OzvPnDnT39+oZ88mLpWWhpgYbNuGrCyYmWH1agAIDcXw4Rg5slUviJAnqKioqK2tNTAwoCxICKmvNe8Iurq627dvHzFixL59+9auXRsSEqJ8SiaTbdmyBYC9vb2/v39AQICNjU3jK2RkYMcO7NiBq1cVR/r1w2uvKbYO4OqdzpyJ9r1fOulgaM88QkiTWvnR2MLCYseOHZ6enitXrnR2dnZ1deWOz549WygUvvbaa8qVy/XdvInYWMTEICVFccTUFNOmISAAzs4QCLBli+I4V+/09u3WRUdIEygREkKa1Jp7hErLly9fvXq1ubn5hQsXzM3Nn3RaXl5edHT0yZO1MTEfcK2ZmMDPDzNnws2NylkRnmRkZNjb2w8aNCiDtxLyhJCOoE03Sz777LMzZ84cOnTo9ddfT0xMFD6e04qKiuLj42NiYvbv319TU6Ojo9uv3+IhQ3QDAzF1KvT12xY4IS1EPUJCSJPalAiFQmFUVNSIESMOHz68evXqTz75BMDDhw937969Y8eOgwcPVldXA+jUqZOvr+/MmTNfeaVGLKZ5CkQzaO0EIaRJbU1LPXv2jIqK8vDwCA0NraqqSktL27t3b0VFBQA9Pb0pU6bMnDlz6tSpXbp0UUW0hLQe9QgJIU1q0z1CpdDQ0M8++0wgEMjlch0dHScnJ39//1mzZvXo0aPtFydEJX77Lfebb+6OHFnz9dcumo6FENKOqGag8p///KeHh8cvv/wyYMAAf39/ql9F2qHCwv7HjvW3stJ0HISQdkY1iVAoFLq4uLi40Adt0n5xO1zSyCghpIGW1RolpOOiREgIaRIlQqItKBESQppEiZBoC0qEhJAmUSIk2oISISGkSZQIibYoLQUAWk9PCGmAEiHRFtQjJIQ0iRIh0RaUCAkhTaJESLQFJUJCSJOoBDbRFtbWePQI3bppOg5CSDtDPULyHMrIqNvkuagIu3YhMxNz5+LCBTg4oKAAMTGaDI8Q0q5QIiTPoVOnMG8eYmMBIC8P69fj3Dls36549uZNhIVpLjhCSDtDiZA8nzw9sXQpSko0HQchpN2je4Tk+eToiBEjsHIl3npLcWTXLpw9CwDl5TAz02BohJD2hXqE5LkVGop9+5Caqnjo44OTJ3HyJDZv1mhYhJB2hhIheW6JRNi4EUuXKh7q6UEshlgMkUijYRFC2hlKhOQ5JBAovnB3x5QpiofKg0r37mH5cqxahUePeA2PENKuCBhjmo6BEM14800sW4bERJiZYeZMTUdDCNEQ6hES7bVpEwwNsXcvxozRdCiEEM2hWaNEe508iQ0b8N//ol8/TYdCCNEcGhol2svDA927w8AAAQGYPFnT0RBCNIQSISGEEK1G9wgJIYRoNUqEhBBCtBolQkIIIVqNEiEhhBCtRomQEEKIVqNESAghRKtRIiSEEKLVKBESQgjRapQICSGEaDVKhIQQQrQaJUJCCCFajRIhIYQQrUaJkBBCiFajREgIIUSrUSIkhBCi1SgREkII0WqUCAkhhGg1SoSEEEK0GiVCQgghWo0SISGEEK1GiZAQQohWo0RICCFEq1EiJIQQotUoERJCCNFqlAgJIYRoNUqEhBBCtBolQkIIIVqNEiEhhBCt9v87MuykVzzxKwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Draw.MolsToGridImage(mhfps_sim_mols[0])" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAADICAIAAAC7/QjhAAAwJElEQVR4nO3deVxV1fo/8M8ZGBUFRVBLDVMTnEhUvII5myZgDihpZE5U9xr3plcx+3ax3+vbN7UyrprlNVMoTVAkQU3FGQdwTMXhmhOaIjMKyHjO8/tj0eEEynTO2Qfcz/vly5dnc85eC1juZ6+113qWgojAGGOMyZXS3BVgjDHGzIkDIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWOBAyxhiTNQ6EjDHGZI0DIWOMMVnjQMgYY0zWjBYIL1zAypXl/87Nxfr1uHQJO3aUH8nMxPffG6soJmtVWxrAjY0xVn9GC4TnziE4GD/+CADZ2fjmG5w/j+jo8q9mZFRcvBgzRNWWBnBjY4zVnzGHRkePxqJFyMkx4ikZewJuaYwxIzJmIHR1xcyZWLCg4kh0NLp2RdeuGDPGiOUwuava0sCNjTFWX0aeLPPhhzh+HEePlr8cNw7JyUhOxvbtxi2HyV2llgZubIyx+lIb93SWllizBuPGwcUFAJRKqNUAyv9mzFgqtTRwY2OM1Zfxl094e2P8eCgUAMr/1klNxT//iYUL8eiR0YtlsqPf0gAoFCgsrGha3NgYY7WkICLTnT09HQsXIi8PW7YAwOzZCAnBwYOwsMDbb5uuWCZHS5fi44/xySf48EOAGxtjrNZMu6Dezg4bNyImBnl5ALB2LaysEBMDLy+TFsvkqF07lJYiKan8JTc2xlgtmTYQ2tjg5Zeh0ZRfng4exPz5+OYbdO5s0mKZHHl6AqgIhNzYGGO1ZPIUa+J+/NgxAPjsMwD417+wZ4+pi2Wy07EjHB3x4AHu3AG4sTHGas3kE+y8vLB8eXkg3LvX1KUx+VIo0K8fdu1CUhLat+fGxhirLYl6hCdOoKzM1EUxuas0OsoYY7Vh8kDo7IxOnZCfj+RkUxfF5I4DIWOsHqTYhkl0CvWTgDBmCp6eUChw5gxKS81dFcZY4yFdIBSPCRkzHXt7dO6MwkIefmCM1YF0gTAhQYKimNzx6ChjrK6kCISurmjZEvfu4e5dCUpjstavHwCcOmXCfEmMsWeMFIFQocDkyZe8vL5NTOQlXcy0Bgy4/9JLASdP9jZ3RRhjjYYUgRBAu3Zxx469d/hwnDTFMdnq0aPVnTuxly6dz83NNXddGGONg0SB0MvLC8AxnjDDTMzCwsLd3Z2ITp8+be66MMYaB4kCYb9+/aysrC5evPjw4UNpSmSy1a9fPwBJPGGGMVY7EgVCKyur3r17azSakydPSlMiky1PT08A3NIYY7UkUSAEj44yqYhAmJiYaO6KMMYaB6kD4VFOMMNMrGPHjk5OTunp6SkpKeauC2OsEZAuEHp7eysUisTExDJOv81MrG/fvuDHhIyx2pEuEDo6Onbq1KmgoODChQuSFcrkSYyOciBkjNWGdIEQgLe3N3h0lJkeB0LGWO1JGggHDBgAIDIyMisrS8pymdz07dtXqVSePHly48aNeXl55q4OY6xBky4QajSa/fv3Ozg4HD9+3NnZuU+fPosXL758+bJkFWAyodFoFi5cSERqtfrNN99s1arViBEj/v3vf6elpZm7aoyxhkhBJEV64pKSkjfffHPLli3W1tY9e/b89ddfS0pKxJdcXV39/Pz8/Pz69++vVEraQ2XPnuLi4qlTp0ZHR9va2k6dOvXq1avHjh3TarUA1Gr14MGDJ0yY8Prrr7du3drcNWWMNRRSBMLHjx9PmDBh9+7d9vb2cXFx3t7ejx492r17d2xs7C+//JKdnS3e5uTk5OPj4+vrO3LkKFtba1PXij17cnJy/Pz8jh496uDgEBcXJ1bsPHjwYPv27dHR0YcOHSotLQWgVCoHDBgwfvz4CRMmt2/f1ty1ZoyZmckDYW5urq+v79GjR52cnPbs2ePu7q7/VY1Gc+LEiR07dmzfvv3q1aviYOfONzt0cPHxwcSJeO45k9aOPTvu378/evToCxcutG3bdvfu3T169Kj0hpycnLi4uB07duzataugoACAp+e2vLxx/v6YPBmuruaoNGOsATBtIExLSxs1atSvv/7avn37+Pj4Ll26VPPmS5cuxcbGJiVdjov7QasFAKUSffpg7Fj4+qLKZY2xCleuXBk1atSdO3fc3Nx2797drl27at6cl5e3a9eu7dtj9+1bm5FhKw52747x4zF+PHr1kqTGjLGGg0wmJSVFRL6XXnrpzp07uuNarbb6D2ZkUHg4+fuTnR0B5X86dKCgIIqNpeJiOn+eVqwof3NODn3/vem+CdYIJCYmOjo6Aujfv39mZqbuuP6/n6iwkGJjKSiInJwqWtoLL1BwMCUkkFbLLY0xWTBVILxy5Yq4K+/du3daWpru+IEDBzw8PPSPVKOggLZvp5kzydm54jplb0/z5xNAP/xARHTjBvXta6JvgjUCsbGxtra2APz8/B4/fqw7vm3bNjs7u927d9fmJCUltGdP5YjYrh399a/c0hh79plklubZs5g9e93du3dfeeWVAwcOODk5iePR0dGjR48+c+bMunXranMeW1v4+eG773D/Pk6fRmgoPDzw8CHatMHo0Vi0CDk5pqg+azQiIvDFF66FhUVvv/12dHS0jY2NOL5y5cqJEyfm5eUdOXKkNuexsMDIkVizBqmpSEhASAg6dcLdu2jShFsaYzJg9NB6+DA1a0YODtpp0/5X/w49IiJCrVYDmDNnjkajqff5796lDRto7lxavJhmzeL7dPn63/8lhYIAWr78vO6gVqsNDQ0FoFAoQkND631yrZaSkmjdOm5pjD37jBwId+wgGxsCKCCASkoqjq9cuVKsEQwJCTG8FBEIi4vJzY3Cw/nyJDtaLc2bRwCpVLR6dcXxsrKyoKAgACqVau3atYYXxC2NMTlQG7FzuXkz3noLpaV45x2sXg3d4vilS5cuXLhQoVAsX778gw8+MFZxlpZYswbjxsHFxVinZI1ASQmmTcPmzbCyQkQEJk0qP66/lD4qKmrMmDHGKpFbGmPPNqM9I/z2W0yditJShITgm2/Ko6BWi4ULHy1btk6tVkdERBgxCgre3hg3zrinZA1afj58fbF5M5o2RVxcRRTMycGoURa3bnVycHDYu3evEaOg4O2N8eOhUBj3rIyxhsEo/colS0ihIIWCli6tOFhaSlOnEkDu7o9iY2ONUlAl/+//kZ0drV9vinM3emlpaWfPnt2zZ09xcbG562IcDx5Q794EUOvWdPZsxfE7d8jNjQDq1k1z5cp/TVH0pk00cCCFh5vi3IwxMzO0R0iEkBAsXAilEt9+iwULyo8XF2PyZGzciCZNsGyZna+vr6ER+0maNkVeHo4dM8W5G7fTp09369Zt1qxZr776qoODg6+v73/+859GnXX61i0MHIizZ9GxIxIS8PLL5cevXIG3Ny5fhpsbfvlF2bVrdUkb6i0jAwkJqN0UVMZYY2NIFNVoaPZsAsjSkjZvrjiel0fDhxNADg504oShsboaSUkEkKurCYtojPbt29e0aVMA3bt379mzp+53rVarhwwZsnz58uvXr5u7jnVTWkpduhBAfftSenrF8cREcnQkgPr3p5pWzxskMZEA6t7dhEUwxszFoECo1dK775KtLe3aVXEwO5v69y8fv7pwwdD6Va+khGxtSaGgjAzTFtSIxMbGWltbA5g6dWppaSkR3b59e82aNT4+PpaWlrqg2LFjx+Dg4Pj4ePGehqZqSpdDh2jsWMrLq3hPbCzZ2hJAfn6kt07HJIqKyMqKlEp6+NC0BTHGpFe3QFj18qTR0KVLf3rPwIEEUKdOdPOmkeqop6SE4uLoiy8qjgweTACZ5hFk4/PDDz+IxZp/+9vfqi7WzM/Pj42NDQwMtLe310VER0fHwMDAqKioR48emaXOT7RhQw0pXbTa8l/97NlUVmb8Cmg0dPbsn8Y5+vUjgA4cMH5ZzwCtVhsfH3/u3Ll0/Q47Y41E3QJhjZcnIjp3jgYMoN9/N0btqigqImtrUigoK6v8yEcfEUDGWJ3Y6K1ataqWizXLysoSEhJCQkJeeuklXUS0sbEZPnx4WFjYvXv3pKlwNTZsoNGjqV07ys5+akvLzqaVK01VgYICUqvJwqKir/n++wTQZ5+ZqsTGKyMjQ0wCcHR0dHJy2qU/QMRYY1DnQFjj5cnUvL0JqBiM3bWLAPL2NkNNGpQlS5YAUCgUX+j3l2vhxo0bYWFhXl5eul2RlUqlh4dHaGjo6dOnTVTbGjWE5EG9ehFAx46Vv/zxRwLo9dfNUJOG7MCBA8899xwAe3t7sfWVUqlcsGBBiX5CDcYatjoHQrNfnhYuJIA++qj8ZW4uqVRkZUVFRWaoTEOg1Wrnzp0LQKVSrVu3rt7nSU9PDw8P9/f3FxNthBdeeCEoKCg2Nlbi61pDSOkSFCTyt5W//O03AqhNGzPUpGEqLS0NDQ1VqVQAPD09b968qdVqw8LCxKNoDw+Pa9eumbuOjNVKfQKheS9PcXEE0KBBFUd69PjTnbuslJWVzZgxA4ClpeWWLVuMcs78/Pxt27ZNnz69VatWuojYokWL7Oxso5y/NkRLI6KEBHJ0NE9LW7eOAJo8ufylVkstWxJAeluKyVdKSoq3t7fo/wUHB+vfJ508ebJTp04A7Ozs1qxZY8ZKMlZL9QmEZNbLU3Y2KZVkY0O6ZeLvvksALVtmhsqYV3Fx8YQJEwA0adJkz549Rj+/RqM5ffp0aGho79697e3tO3ToINnafF1LI6I336Q+faQp9k+Sk8u3J9QZNYoAMtL9RiO2bdu2Fi1aAHj++ecPHz5c9Q0PHz588803xS3UxIkTpbyFYqwe6hkIjxyhnj3Jzc0kdaqRSCOSmFj+8ocfCKBx48xTmSeSIGDk5+ePHDkSgIODwzHTd4e7du0K4IRJl4U+yeTJBFBcnMTFEhFpNNSsGQH04EH5kdBQAmj+fDNUpoEoLCwMDg4WEW7s2LHVb30cFRUl5id36NDh6NGjklWSsbqqW2aZadPw5ZcAEBeHCxcwYUKdPm00Xl4AKhLKeHvn9+r19zNnXInIPBXSU1xcPHv27A8++MDFxeWdd96Ji4srKSkxeik5OTkjRozYu3dv69atDx06NGDAAKMXUckrr7wC4JjkWXxEnuukJImLBQClEn36AMDJk+VHvL2zBw2K/f33JWaozVOYonU9zZUrVzw9PVesWGFtbR0WFhYTE9OyZctq3u/v73/u3Lm//OUvKSkpgwcPXrx4sUajkay2jNVB/eLnzz8TQEOHGjUo19rmzfdeeWXTrFmf6o48//zzAP77X5Pkmay9vLy8YcOGAdCfb9KiRYupU6dGRkY+NNJi7NTUVJEv5oUXXvjtt9+Mcs4ahYeHAxgneb87JoYAGjFC4mLLffZZ2uDBx5Yt2yFeZmVlKRQKW1vbhpCFIDMzc/LkyZMmTXJzcwsJCUlISDBpceHh4ba2tgC6du167ty52n9QzKkRc5KHDBnyu4lWVjFmgHoGwowMUiioSRMyywXht99+A+Ds7Kw74u/vD8CQOZOGy87OFj0zZ2fns2fPiqdrHh4euoioUqm8vLyWLFliyGy6W7duiZkIrq6ud+/eNWL9q3f9+nUArVq10mq1khVKRPfvE0DNm5MBeznX388//wxg+PDhuiPih1+nSGAKt2/fFoPVIjgJnTt3DgkJSUpKMu7vKDc3NyAgQBQRGBiYn59fj5Ps37+/bdu2AOzt7SMjI41YPcYMV/8Uay+9RACZa6VZmzZtAOj6Q2FhYQBmzJhhntoQpaam9urVS/TSKsW5mzdvhoWFDR8+3MLCQnfN0t3F1+madfnyZbFmq0+fPhmSp5UTP3Pp58S3a0cAXbkicbFERKmpqQCaN2+uS9MzdepUAN9++60ZavOH5OTkdu3aAejWrdvt27cTEhKCg4PFb0do166dWPRieM81MTGxY8eOAJo1a7Zx40ZDTpWenu7j46MLqAUFBQbWjTFjqX8gnDGDAAoLM2Jl6mD8+PEANmzYIF6eOnUKwEsvvWSWyty6datz58419tKysrKioqICAwObN2+uu2Y5OTkFBgbGxsYW1bQQ8tSpU46OjgAGDx5srFHWOhEzVNdLvuvVxIkE0B+/aqm1b98ewOXLl8XLFStWAJg+fbp5akN04sQJ8WRu0KBBubm5uuMajUZERPGYQGjZsqVoXfWYvSUWBYq7t759+xolUbtWq12zZo3oxbq6uv7666+Gn5Mxw9U/EIpVVv7+RqxMHXz55ZcAZs+eLV6WlpY2bdpUoVBIn+rw8uXL4tLj4eFRy15aaWmpuGaJ+3rB1tbWx8dnzZo1qampVT9y8OBBOzs7AL6+vo9NnWH6KZYvXw5g1qxZEpf7+ecE0HvvSVxsOTHqrgv/SUlJokNvlsps377dxsYGwOuvv15NM0hOTg4NDdXPn+fg4ODv7x8eHl7Lgc0HDx6IOckKhSI4ONi4s6CTk5NFDhox6UbiwXbGqqp/ILx6lQBq29aIlamDxMTEStejoUOHAti+fbuU1dD10gYNGlS/XlpycvKSJUu8vLwUf2x/rlKpRIYzXS9Ef0MJM2auEjHAVfJdr44cIYB695a42HKff/45gHfffVe8LC4utrKyUiqV+r0xaaxfv15kVJ8xY0YtxzxFRNR/UG1jY+Pj4xMeHl5Nc92zZ4+zs7MYrvjll1+M9x1UePz4ce2XYTBmavUPhFottWpFAN26Zbzq1FpJSYmtra1Cocj6I/32xx9/DGDBggWS1eH+kSOil+bn51dYWGjg2VJSUlatWjVy5Ej9zZJcXV19fHzE5S84ONi89866n7nE3W5d/muzPFQ6cuQIgJdffll3xNPTE8C+ffukrIbIJYtaZFR/IvGgWv9+y9raWoxApKWl6d5WVFQUHBws3jN8+PD79+8b7zt4At3CfGdn5927d5u0LMaqYdB+hGPHVmxGIb1BgwYB0P3/2b17N4ABAwZIVHxcHNnYfDp48JQpU4zbSysoKIiNjQ0KChJ35ba2tmq1un6XP6MbPHgwgFjJd71ydyeAzLImu6CgQK1Wq9Vq3aCi6Mp8+umn1X/QaMrKtv3rX2Ko4D//+Y+BJ0tJSRFTt8TdFf6YzBwWFnbkyJGXX34ZgFqtDg0NrbqNlyncuXNHLFE1xRgsY7VkUCBctowA+mPQSGrnz5/Xn5ny8OFDlUplZWVleOesZhERpFYTQHPmaE12vSgtLd2/f3+zZs0ApKSkmKiUOvnoo4/q3SkxxDvviJlZ5lm9t379+oSEBN3tzsaNG8UwgBRlFxbS+PFaheKdIUNiYmKMeOL79++vXr1aPyKKOVydO3eWeNeRsrKyxYsXq1QqGxsbydbFMqbPoEB4/HjxX/6SEhAQY6TKGEqs8ZozZ45p/zt9/TUplZLtgih2ejNw5rqx7Nq1C4C35Ltebdp0rG1bz0mTJklc7hOJJZUODg76g4omkZdHI0YQQPb2dOSIiQrJzMz8/vvvhwwZYmtra29vb679mUU/+/333zdL6UzmDAqERUVF1tbWSqWyISTV3bhxo1qtdnJyEre3HTt2DA4Ojo+PN3ISkCVLCCCFgj7/3Jinra7AJQD++te/SlNc9XJzc6Xrduu5dOkSgA4dOkhZ6NOEh4fb2Ng0a9ZMqVSKDAlGWVpQWWpq+YhwmzZk+mUGGo1GrEQ015KGgQMHAvjpp5/MUjqTOYMCIRF5eXkBMNHUstpbvXq1yOE0fvz4KVOmODg46C/UmzFjRkxMjKELeLVamjePAFKpaO1aI1W8ZgkJCQB69eolWYnVExPfJUjzrU+j0YiBO1NP36jRkiVLxFwSV1dXKysr0cYUCoWnp+fSpUuNFhFv3KBOnQigrl1JqlHx9957D8DHH38sTXH6Hjx4IG6wzLJAljFDA2FISAiA//mf/zFKbepHtzn7sj+2YiotLT1w4MA//vEPkRRDf+L4/ogIqsfFtKyMZs4kgCwtJd6GR9ftln6+/hO9++67AJZJvuuVSOIq8fIYfVqtdv78+aKlffHFF/THtKbAwEAxeVhwc3MLDQ016DHbqVPk5EQA9e1LEk7Q3bdvH8yUleKbb74B4OvrK33RjJHhgTA2NhbA4MGDjVKbutJqtf/85z+rn1B348YN/WlyZwcNIoDc3Cg0lE6fptosSCguLs9u0qQJmWOSt0hh2kDml//www9i7ZfE5S5atAjAokWLJC5XKC4ufuONNwBYWlpWHb4rKCiIjo6eOnWqfs6gtDFj6MMP65yEcN8+srMjgIYPJ2kf15WVlYmtmJOTk6Usl4hGjBgBcyQtYkwwNBBmZmYqFAobGxvp5z2XlZXNmjVLXJuioqJqfP+DBw/Wrl2bN20a2dgQUP7HxYX+/nfav7+69OFpadS5Mzk4kOS78QmiI2KWMauqbt68CcDR0VHiRY0i//WwYcOkLFTIz88fNWoUgKZNm1a/AXJRUdHOnTtnzpzp3aULqVTlbax9ewoOpvj4mlPUR0aShQUBNG2aWfLZi/9QixcvlrLQnJwcS0tLtVotffpcxgRDAyH9sWXryZMnDT9V7RUXF4vcV7a2tnV+Qvn4McXHU3AwtW1bERGbNCEfHwoPp5wcSkujzz6jmTNp6VISOS9u36YLF0zxjdSGiAFDzbXrVRUipdzVq1elLDQtLQ2AnZ2dNOvbdDIzM/v37y8WfZ85c6a2Hysro4QECg6m556raGOOjhQYSLGx9LS7xgsXyN6egoPNs9cG0S+//AKgR48eUha6fv16ACPMtdUWY0YJhDNnzgTw1VdfGX6qWiooKHj11VcB2NvbG7TztUZDp09TaCh5eFRcrWbPpk6dKCCANm2iCRPIxYXMlNtTJyMjQ6FQNGnSpCFsg0fm2/SqUv5rCdy8ebNLly4AXFxc6rnthkZDx47RvHnk4lLRxlq2pOnTKSXlCbdcEm6tVVVJSYlI9XJFws0+xAKhb775RrISGaukbjvUP5GYOCrZ3uW5ubkjR44U6RAPHjwoSq8npRIeHli8GKdP4/p1LF+OwYPRqhVcXPDTT3jjDWzdirZtceCA8apfH46Ojl26dCkoKPj111/NWxNB4t+4jshtliTVdvXJyckDBw68du2ah4dHYmKi2GCkzpRKDBiAL77AzZtITkZoKNzckJWFH35AaSm8vHD+PIYNw8mT6NsXhYXQ2ztCehYWFn5+fgC2bdsmTYn5+fnx8fFKpXLs2LHSlMhYVUYIhN7e3gCOHj1q+KlqlJaWNnjw4GPHjnXo0CEhIcHd3d1op37xRXzwAQ4exP378PWtOP7aazh40Gil1Je5Ys8TySEQHjp0yNvb+969e0OHDj1w4IBufapBunXD4sW4dAmXL+P773HqVEO75QIgNtuKjo6WpridO3cWFRV5eXnp76fImMSMEAg7d+7cunXrBw8enD9/3vCzVSMlJWXgwIHnz5/v2rVrQkJCPe/Qa/TCC8jKqnh54wZ69jRJQXXRoAKhu7t706ZNr127lp6eLmW5ffv2BbBt27bo6OjHjx+brqCYmJjRo0c/fPhw/PjxO3fuFFnujMnVFYGB2LOnAd5yjRgx4tWuXRc3a6a5eVOC4kTEFduLMmYuRgiEROTg4NCuXTsPDw9vb++lS5deu3bN8NNWcvXqVW9v799++83Dw+PIkSP6O/kZ2bhxiIzE9esAkJ6OPXswZoypyqq1BhUI1Wp1v379iOjEiRNSlrtv3z5ra+v09PSJEye2bNnS19c3IiLi4cOHxi3l66+/njhxYlFR0fvvv79lyxaxAZZJNMhbLisrq90eHr6HDqlM3yksKioSufJ5XJSZmeGPGVeuXAnA2tpal70XQK9evT7++ONTp04ZZZL96dOnxQqnV155RYrcE/Pnk1pNfftSu3b03XcmL64WtFqtGJ27ZZZdr6oQm17Nnz9fshKTkpIsLCyUSuXs2bP79++v207I0tLytdde++677wyffK/VakNDQwEoFIrQ0FBj1Lpa589T164k8uKmpdFzz1ED2ZZv2zYCyNPT1OWI6dB9+vQxdUGMVc/QQHj58mWxZXZUVFR2dnZUVFRgYKD+suJWrVoFBgZGRUXVcmvsqg4dOiTGpsaMGSPd5uwFBXTtmlnWcj2NuGv+wVy7Xv2ZmEwh2aZXubm5Li4uABYuXCiOpKenh4eH+/j4WFhYiJam207o999/r0cRZWVlQUFBMNJuR7XV8G65iIgKC8nOjhQKun3bpOW89dZbAP7v//7PpKUwViODAmFJSYl4bDNjxgz942VlZQkJCcHBwfoDmCLD2Zo1a+qULjIuLk4EWqNv+9foLFu2DHpbpZtRRkZG7969ra2tJUt2M2XKFNF1qJq3ITMzU0RE/Q2NRZ6z//73v7U8f1FRkZgkYmtru2PHDmNXv1oN75aLiGjSJAJo+XLTlWCWpRqMPZFBgXDhwoUAXFxcqhmuTE5OXrJkif7W2Eql0sPDozb5GDdu3Cju99977z2Jl1E3QOIBocSLnau6f/++yLvt4uJikl0Xqli3bh2Apk2bVh/YcnJyxIBE06ZNK0XE6ltadna2mPns4OBg0LLUZ0lkJAHk5WW6EsTTwe7du5uuCMZqqf6B8OjRoyqVSqlUHjp0qDbvT0tLE3fuupz90NssqWpvT7ehRAPZnN3sGsKmVzdv3nzxxRdFgKnfCGRd/fbbbyKldaUx4RMnTqQ/JSH148ePRS5s/SF60dISEhIqPbS+d+9ez549AbRt2/aC+ZIHNTgFBdSkCSmVdO+eiUoQA9FSPItlrCb1DIR5eXliF9yPPvpI//jx48eriW2CyNkfFBTUunVr3XWqRYsW/v7+4eHhonNZdUMJRkSi47Jr1y6zlJ6cnNy2bVsAffv2zZRkWkdJSUm/fv0ABAYG6h/XarUdOnSocTvAoqKi+Pj44OBgZ2dnXUvr0KGDiIgajeby5csiW42bm9udO3ck+I4ak3HjCKCVK01xbo1GI/77nz9/3hTnZ6xO6hkIp02bBqB3796Vntn861//0l1xHB0dp02btnXr1ry8vCeepKys7NChQ/PmzRMxVbC2thZJrVQq1XcNZ/pAwyA2vap08yGNkydPtmzZEsCQIUMk28R87ty5ojNXaew9Ozv7tddeq/12gKWlpfHx8e+9957+qm1nZ2cxiDpw4MCGsLN0g/PjjwSQaTaWOXTokBhdN8XJGaur+gTCmJgYEbEuXrxY6UtlZWWnT58ODQ11c3PTj23Dhw8PCwu7+/Q8ivqbJdnZ2dnb29dmQwm5MdemVwcOHBDjk35+fpLtTb9vX2rTpnYWFhZJSUlPfEM9tgPUaDSifYpsDB06dBg0aJB0U5Ebl0ePyNqalEpKTTX6uYODgwEsWLDA6GdmrB7qHAjv36fBg/c6O7dZWdOYSaWNAHXXqZCQkKqPanQuXLgg7tbrWjE5yMrKUiqVEm969fPPP4tF5W+++aZkWb8fPCBnZ3Jxufrvf0fU+OYnPhR0cXF54kNBnaFDhwLg+63q+PoSQN9+a9yzarVaMSKdmJho3DMzVj91C4RaLY0aRQCNH59d+5XymZmZYjqffqqq9u3bBwUFxcbGFhUV6b85IyNDDKvWqWLy4erqCuBpnSSji4iIEPcxc+bMkWzirlZLPj4E0KBBVFZWhw/q1u3oP35u3769eGhdKYqL9RJbt241cu2fJQcOUEQE5eYa96wiW+xzzz0n8ZaWjD1N3QLhihXlu6rVZSlghcePH+/YsSMoKEj/UU2zZs30Z47l5uYCaN68eX0KkAGxdepyUy7w0lm5cqVZJu5++WV5M6v3jMWysrL9+/f/7W9/E7N7hDZt2uTqXdMnTpwIYMuWLcap9DOp6i5RxrBgwQIAwcHBxjohYwaqQ67RK1cQEgIAq1ejfpnibWxsxowZI9bUJycnh4aGenh45OXliSXzguh/lJWV1acAGRBJR3fu3GnqH9HSpUvff/99Ilq+fLmYxCuNixfx0UdQKLBuHfSiWN2oVKqhQ4euWrXq3r17oqV16dKlbdu2+mOnrAa5uU/YJcoYRGY1TrTNGg51zW8BAJSVYdo0FBZi+nT4+xuh4G7dunXr1m3x4sW3bt3SXwEtVtCXlpYaoYxn0cCBA0eOHLlv3z4nJ6fhw4f7+PiMHTvWuNd3Ipo7d25YWJhKpVq7du306dONePLqFRRg0iQUFSE4GH5+xjmnrqVlZmbqHxcZHojIOMU8e3bvLt8lCsAbb8DbGwcO4PffsX9/9Z8LUCi0T/+pPnz48Nq1a46OjgbtJMqYcdWy5/jhhwSQiwuZOuW1RqMBoFQqTVtMY3bmzJmuXbvqfoOWlpYjR45ctWpVSkqK4ScvKysTkc/Kyio6OtrwE9bJ9OkEUPfuJMFETn9/f/BkmWq8/TatWFHx8tNPad48mj2bgOr/KP/IIfU0o0aNioioeQ4UY5KpVY/w2DEsWwalEuvXw+hbs1WiVCqVSqVWq9VqteIBFaukd+/eV65cuXXrVmxs7I4dOw4fPrx37969e/fOmTPHzc3N19fXx8dHP6dd7RUXF0+ZMmXbtm1NmjSJiYkZMWKEKer/NFu3Yv16WFtj0yboDZabCvcIa1B1l6hBg9CtG2pqFZFA9T9TFxeXPn36GKGGjBlLjaEyL486dSKAJFvGLRZKV5pNyp4mKyur6qYfTk5OgYGBVSflViMvL09EPgcHhxMnTpi0zlXduUMtWhBAa9ZIVOKkSZMAREZGSlReo9Ngd4lizNhqDoTTphFAvXuTZEvXmjRpAuBp+WjY05SWllbd9MPW1lZs+pFa7bLo7Ozs/v37A2jdurVZUm6GhYllOdKVOHnyZACbN2+WrshGp2HuEsWYsSmo2qEhIixejK++QlISXF1N0iWtysHBITc3Nycnx97eXqIinzmXLl3asWNHXFzc8ePHxa9YpVK5u7v7+PhMnjzZ9c+/y9TU1FdfffXixYsuLi579+7Vz3gnpbg4eHmhRQuJinvjjTc2b978008/BQQESFRkY/T4Me7dg4sL1LWdWMdYo1PxEO7CBaxcWf7v3FysX49Ll7BzJz75BLduoVUrfP+9RHVq335gu3aDSku1EpX3LOrWrVtISMjRo0dv3769Zs0aHx8ftVp95syZTz75xM3N7cUXX/z73/++b9++srKyW7duDRw48OLFi25ubgkJCaaOgk9sZjt2AICvL7Ra6ZqZnd2L7doNsrBwkqi8RsrWFp07cxRkzzhd33DDBgJI7HVz4wb17UsbN9Lbb5d/9fJlcneXqJfapg0B9Vyzz57m4cOHkZGRU6dOdXBw0P32HR0dRbqf/v37S5N4uuE0szfeIIA2bZKoOMZYg/WnaZmjR2PRIuTkSB2MK7GwAABeSWhczZo1mzRp0o8//piRkZGQkBASEtK1a9fMzMyhQ4cOGTJk7969+gHSpBpIM2OMMeFPIx6urvD0xIIF+PDD8iPR0ThxAgBKSiBZUg4xDMO5ZUxEpVJ5e3t7e3svWbLk0qVLzz//vK2trchjII0G0szE6hJePcEYq7xQ78MPcfw4jh4tfzluHJKTkZyM7dulqxP3CCXTrVu35s2bSxkFhYbQzBhjTKj8DNzSEmvWYNw4uLgAgFJZ3j+T8mE59wifeQ2hmXGPkDEmPCF1i7c3xo8vv0yIv7OzJa0T9wjloFIzy8/H9Ok4e9bc1WKMyU/FHfi0aRVH16wp/8eUKdixA9Om4euvce6cVHXiHuGz62nN7PPPsWED9u/Hr79KVBPuETLGhJqTed67h+xszJ6Nq1clqA/APUJZ+uADDBiAu3cxa5a5q8IYk5maA+E77yAwEPn5mDQJjx9LUKXyHiEHQllRq7F5M1q0QEwMvvlGihK5R8gYE2q1vcPq1XB1xcWL+Mc/TFwdAH/0CHloVG7atSsfLJ07F+fPm7s2jDHZqFUgbNoUUVGwtcXatYiIMHWVeGhUviZOxKxZKCrClCkmH37gHiFjTKgh6ba+deswaxaaNMHJk3BzM7TgCxcQH495857wpWPHUFgId3c4OhpaCmt0iorg6YkLFxAUVDGbpn7y8rBzJ+Li8P33sLKq/NVLl5Caim7d0KaNQaUwxhq7OgRCAG+/jfBwdOuGkydha1uf8lJS8PPPiIgonyh/4QJ69Kj46oULOHwY778PALm5iInB9On1KYU1apcuoW9fFBZi0ya88UadP56bi/h4xMVh2zYUFABAbCx8fSvewM2MMfYndcpMmp9Pbm4E0JQpdUtpeucOff45eXgQUP6nVSt67z26fv1Pb6uakZnJ09dfE0D29nT7traWH8nMzAwPL3rtNbK0LG9jSiW98gqtWEEPHvzpndzMGGP6avWMUKdJE2zZgiZNsH07Nm8+WeP7s7KyIiIiRowY4eeXO38+zpxB8+YIDERsLO7dw+rVePHFyh/hjMwMwF//inHj0KPHlcDAUaXVPi4WbczX17dt27Zr197ZtQsaDby8EBaGu3fLe37OzpU/xc2MMaZTt6FRYdOmx6GhU+7fj09KSurevXvVN+Tk5Gzbtm3z5s0HDx7UaDQAhg6NaNUqMCAAo0c/4WmNRoPDh7FjB3r2xMWLaNYMv/+ODz9EQABO1hxt2bMpO/uhu3uPu3fvLlq06NNPP6301fv378fExGzdujUhIUG0MQsLi7feWtW/f9Drrz/56XJODmJjsWsXRo7E5cvczBhjf6hfR3L69OkAunTp8ujRI/3jhw8f9vX1tbS0FCe3srLy8/PbtGlTfn7+E8+TnEwhIeUbEAK0cCHNnUvFxeTmRuHhPGYld0eOHFGpVEqlMj4+Xv/4ihUrlEqlro35+Phs2LDhafspZmVReDj5+FQMmc6dy82MMVahnoGwsLDQ3d0dQEBAgP7xqKgoACqVysvLKywsLCMj44kfP3uWFiygDh0qHhl26UIff0xffklz5xIRJSSQoyNfoRgtXrwYgLOz8wO9B32JiYnW1tY+Pj7h4eG5ublP/GBqKq1eTcOGkUpV3sbUahoxgr79llat4mbGGKtQz0BIRNeuXRObm69du1Z3sKCgYPXq1enp6U/8yM2btGQJde1aEf+ef56CgykhgbRaIqING8qvUEQUFET9+tW7duwZodFohg4dCmD06NFabfnEGa1W+7Qxhnv3clesoFdeIaWyvI1ZWtJrr9G6dZSZWf4ebmaMMX31eUaoExkZGRAQYG1tffz48Zdffvlpb7t9+3ZkZOT27fvPndtbVAQAzs7w90dAAAYMKF/XzNjT3Lt3z93dPTMz86uvvvrHU5Ib3blzJyYmZsuWLWfPXlars/LyFNbWGD4cvr4YNw6tWklbY8ZYo2JQIAQQFBS0du3aTp06nTlzRnQQdVJTU6OioiIjIxMTE0Upw4bdfuGFDgEBGDIEKpVB9WaysmvXLh8fH7VanZCQ4OnpqTt+48aNrVu3RkdHnz59WrQxW1vbmTMvDhjQccwY2NmZr8aMscbD0EBYVFTk5eV19uzZSZMmRUZGAsjJyYmLi9uyZcvu3bvLysoA2NjYDBs27K233ho7dqxuHg1jdTJnzpyvv/5a3HJlZWXFxsZu2bLl+PHjogGLNubv7z9u3Dg7DoCMsbowNBACuH79ep8+fR4+fPj222+np6fHx8eLhV/W1tajR48OCAjw9fW1sbExRm2ZfBUVFfXv3//8+fMtWrTI/mOraHt7e19f3wkTJrz66qvW1tbmrSFjrJEyQiAEEBUVJR4WFhYWqlSq/v37v/XWW5MnT27evLnhJ2dMuHLlyrx5806dOqXRaHx8fPz9/UeOHGlVdV0qY4zVhXECIYCLFy/u3bu3adOmEydObNmypVHOyVhVV65c6dy5s1rsWskYYwYzWiBkjDHGGqO65RpljDHGnjEcCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmscCBljjMkaB0LGGGOyxoGQMcaYrHEgZIwxJmv/H9F7PUBMswf1AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Draw.MolsToGridImage(mhfps_sim_mols[1])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.6" } }, "nbformat": 4, "nbformat_minor": 2 }