rdkit_edit_molecules.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": [],
      "authorship_tag": "ABX9TyNSqlwZyPIoSFT20mgIOKp0",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/alihojatnia/rdkit/blob/main/rdkit_edit_molecules.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 1,
      "metadata": {
        "id": "cBnYUjgHFm_U"
      },
      "outputs": [],
      "source": [
        "%%capture\n",
        "!pip install rdkit"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "XGO94Ym4b7I4"
      },
      "source": [
        "from rdkit.Chem.Draw import IPythonConsole\n",
        "from rdkit import Chem\n",
        "from rdkit.Chem import AllChem"
      ],
      "execution_count": 3,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "mol = Chem.MolFromSmiles('ClCCCCCCCO')\n",
        "mol"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "SlkSQkRUOMgy",
        "outputId": "bd6d8eb8-b853-458e-92f8-4605edf9666e"
      },
      "execution_count": 4,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c82a20b0>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAQ1UlEQVR4nO3de0xUVx4H8N8MICKggO8H2hYFbdUioFifKcpqWxuzcQmru6Tdbhc1JmxSTWmziWy2myyxiY7+06VdjbP+0fpKt2O1rdTWB1IEB5BS5aEgCoIgb+Qx48xv/zh0ekVAypmZOzN8P/EPnbnc+xM433PvOfee0TAzAQDAcGnVLgAAwL0hRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApCBGAQCkIEYBAKQgRgEApDgjRm+bbv+1+q8H6g844VgAMOJ0dFBTE1mtah3fbjFa1lO2v35/yt2Ut6vefrfmXUOroYd7xFu15toD9QeOtxy317EAAOjKFUpMpJAQCgyk8ePJx4diY+nAAerpeWyz//6XQkLoT3/qfyf//CeFhNDu3TKFeMt8sdBiadl6Z+vx5uNMbHvxw/sfzvGdc+SZI7H+sfKHAAB4zAcfUFoaMdOMGbR6NY0ZQxUVlJdHubl06BCdOUPTpvVu2d1Nzc3U3t7/fjo7qbmZOjtlapGN0S5rV1x5XEFnQZBX0DuT3okfGx/iHVLVU3Wo8dBnzZ9trNh484WbkocAAHjM4cO0ezf5+tJHH9Gbb5JG0/t6aSlt2UL5+fTb39Lly+Rth9PEoZA9zN9r/17QWTDJe1JWRNYc3znixXDf8Pix8XGBcZFjIgO0AdJFAgD8rKuLdu0iIvr3v+nNNx97KyKCzp6l+fMpN5f0evrzn51TkdTYaKe1M+NBBhHtmb7HlqE2f5nwl8VjFsvsHwCgr//9jxobKSyM3nijn3fHj6d33iEiOnjQaRVJxWjOw5xWS6u/1j8xONFeBQEADCYri4ho3bpfruX72LCBiMhopK4u51QkdVFf1FVERAv9Fo7WjrZTPSNLbm5uc3PzunXr1C7EtdTU1BQWFkZFRU2dOlXtWlzLN998ExISsnjxyL7Iq6ggIoqIGHCD8HDSaslkoupqmvPzVXJ5Oe3b18/GeXnyFUnFaOOjRiKa4jNFvo6RpqysbNOmTcXFxUQUFBT08ccfJyQkqF2U+iwWS3Jy8uHDh61Wq1arTUxM1Ov1Pj4+atelvqNHj27durW1tZWIwsPDv/jii7lz56pdlEra2oiIAgMH3MDLiwICqK2NWlt/ebGoqPdi3wGkLuotZCEiL/Ia+pfUmmsPNR6ykmo3yqrObDbv27cvMjJSZKhGo2lpaUlMTNy2bVtDQ4Pa1akpKysrMjLy0KFDIkOtVuunn34aExOTnZ2tdmlqqq+v37p16+bNm1tbWzUaDRGVlZVFRUXt37/fbDarXZ0a/PyIqO/NoX10dxMR+fv/8sqKFXTmTD9/Eu0xIMkS0uvSyUhrytYMvll2RzYZaUXpCmZ+4/YbZKSoG1EX2y/KHNpNZWZmvvDCC+I7P23aNIPBUFlZGRMTM2rUKCIKCgpKT0/v7u5Wu0xnq66uTkpKEhkxbty47du3P3r0aOfOncHBweJ7tWHDhsrKSrXLdDaTyaTT6caNG0dE3t7eixYtun379unTp0NDQ8W3JTw8/Msvv1S7TKfbsoWJ+G9/G3CDujomYiJ+8ICZOSODiXjTpv43fv99JuKdO2UqkorRk80nyUhTiqYMvpkyRj9r+iz0x1Ayksao+X3F7+/03JEpwI2UlZXZLtvnzJlz7Ngx5bulpaUbxLg4UXh4+KlTp9Sq08kePnyYnp4eEBBARH5+fqmpqe3t7f2+O2bMmD7vejZlj7t27dri4uI+7z7//PMDvevh9u5lIo6LG3CDkyeZiMPCev/p4jFaa6r1yvciI13uuDzIZsoYZeZOS2d6XXpAQQAZya/AL7U6td3iyW2jubk5NTXV19eXiAICAtLS0gY63xy85Xgeg8HwzDPPPPV88+7du7Zz1enTp+v1eqvV6txKnWqIfaryXNXHxyclJaWlpcXJpTpVZ2fvXyoqWKtlLy8uKel/y/h4JuLU1N5/uniMMvPGmxtFRJqspoG26ROjQrWpOqkySWPUkJGmF03XN+qt7Gltw2Kx6PX6yZMnE5FWq01KSqqrqxv8S55sG83Nzc6p1pmMRuPKlStFUkRFRV28+PQRnpycnKVLl4ovWbJkyQ8//OCEOp1M2eMOcYTnwYMHKSkpXl5eRDR+/HidTvfo0SPnVOtUBgPPmsUnT/b+8w9/YCJesoSf7Dl0OibisWO5urr3FdeP0fLucv8CfzJSfHl8aXep7fVGc6Puvu5UyykeIEaFKx1XlpYsJSORkZaULPmhw3Paxvnz5yMjI0Wzj42NzcnJGfrXKttGSEiIJ7UNmWY/jG7JXYj/2qRJk4b3X8vPz1+1apX4ZVu0aNGFCxccV6qz5eXxsmW9Y52/+13vi42NPG8eE3FoKO/dy9nZfO0anzjBGzcyEXt7s3LQzPVjlJkvtF8IuRYiojC8OHxl6cqw4jDvfG8y0m/Kf8ODxigzW9mqb9RPKZoiBkyTKpNqTbXyValIeRE6Y8aMYV+EFhQU2NrGvHnzvv76a7uX6kz2ughtb29PS0sTp2z+/v6DDJK4i++///7FF18UP+jVq1cXFhYObz8Gg+HZZ5+1DZJUVFTYt05nu3ePk5PZy4uJeMIE1ulY2eM+eMBbtrBG05uwtj/h4ZyZ+dh+3CJGmbnR3PiPe/+ILYmdeG2id773pGuTVpWu+qD2gxpTDTMXdRZF34h+u+rtQfbQYelIu5fmm+9LRvIv8E+7l9Ztdb+28fDhw7S0ND8/PzElkpaW1mkb0Bkuz2gbfaZEfvrpJ8kdKqfsZs+e3WfKzl3cuXPHLj2uTWdnZ3p6emBgoG3Krq2tzV7VOo/JxDodjx3LROzjwykp3Nra/5ZVVfyf//Du3fzuu/zhh3zpEj95cXPrFh89ytnZ/e+hqIiPHuXhdl2CfWLUXsq7yxMqEsSJ7ezi2cea3KZtWK3WY8eOzZo1i4g0Gk1CQsLt27fttXNl2xg1alRKSoobtY2SkpJXX31V5F1ERMTp06ftuPNvv/12/vz5Yudr1qwpKiqy484dSvS4o0ePtvW4XV1d9tq58gayadOmudmknMHAYWG9p5Zr1/L162oX9HSuFaPCubZzC64vEGEaVxZX1OnqbePq1avLly8XjTk6OvrSpUuOOEpNTU1ycrJWqxVtIyMjw2KxOOJA9tLU1JSamipuiQ0ODk5PT+/p6bH7Ucxmc0ZGxoQJE8TNlcnJyQ0NDXY/ih2JHnfmzJmO6HGVcnNzX3rpJfFruXjx4uyBTsdcR0kJv/JKb4BGRLBde1yHcsUYZWaz1ZzRkDHx2kQykne+d3JVcr25Xu2i+nHv3j1btE2dOtUJ0Zabm7ts2TLRNmJiYi5fHuxWM7WICZOJEyfaJkzu37/v0CM2NjampKR4e3vbJuXMZrNDjzg8eXl5yh43KyvLoYezWq16vX7KlCkispOSkmprXXLioamJU1LY25uJODiYdTp2yR/fQFw0RoWmR00pd1PEbFVwYbDuvs5sdZVvbk9Pj06nGzt2rG3CpHWg4Rt7E6cz4lEWcTpTVVXlnEMPxXfffbdw4UKRFC+//PK1a9ecdugbN26sX79eHHru3LlfffWV0w79VM7vcW06OjpsAwhiUs6OAwiyzGbOyOCJE3tn2JOTud4VT5gG59IxKtzouvFK+SviGn/uT3PPtJ5RuyI2GAzPPfecbdrn5s2bzq/BBdtGVVVVUlKS+LaEhobq9XpVyujz07l165YqZdgoe1wxtO20HlepvLzcNik3c+ZMtX46jzl3jhcs6L2Kj4tjJ/a49uUGMSoYWgxhxWFkpKjDUa+//roqycVPnO+cOaNypoupXtWTy9Uy3UWSi10v08+dO6fWtcJjyss5IaE3QGfPZve80cLGbWKUmbut3f+q/deCmAVE5Ovr+9577zlzwtqVR99UvI7uM/qWkJBw546rrJOg4nU0M1+/ft01RxjEpJxy5LreidfRbW1t1vffZ19fJuLAQE5PZze/7ZfdK0YF0TbEYzCibTj6CZ8n54Kd+Ws3RE8+BuPoWR12k7lg59xHoSR6XOVDaK7T49o45z4KJVuPW7xqFWs0nJTErjnf9eu5X4wKV69eXbFihWgbUVFRjmsb7nVnotPahnvdmejQu3qV3O7uK4fe1auUlZUVExMjDvTH9es5L89BB1KFu8YoP942xMCTfduG+z4no2wbdl+S0n2fk3HEM2ZKfXrcH3/80Y47d6jMzMx58+aJyu3yjJmSssf11AW63DhGhT5twy5LUooJE3d/atvuz1+yRzyZaq8VD5QGX0zWLYgVD5Q38Mkvu+e+Pe6v5fYxKtirbSjXEHLp25WHxo5LUnrYGkLnz59XrgZSUFAwvP0o10kZfDFZt9DQ0GCvZff6LCbrjj3u0HlIjAo5OTmxsbHiJ/dr16Zj5itXrnjkipaSS1J66oqWksvuefCqfcNYDVbJw3rcofCoGOXh/nKPhPXVlW1jiL/cI2F99WGslMxPdNge0+MqDfGzCZQ8tcd9Kk+LUWHol1oj7dN+ht42RtQnmig/t2PwwU1HDK26rME/KUtpJPS4g/DMGBWeOvA/jP7WA4iBf2Xb6DPwX1pa+tprr9km+kfO5+sN3nOMtB7X5qlT7Y6YzHQvnhyjQmZm5pO3oRw5csT24jBGfzyAsm1Mnjx5+/btFoulsrIyPj5+JH/as8lkst346ePjs3z58ps3b1qt1p07d06fPn2k9bhKypmDiIiIw4cPM3NJSYmyxx2Jn/bMzCMhRpnZZDLt3bs3KCiIiLy8vMQdGCI+Dh486OKrdjrUpUuXoqOjxXdj9OjRIlW1Wu22bdtc/L5xh6qvr7c9RarVasVaAeTCKxM6h8Vi+eSTT8RTpEQUGBgohkGDg4N1Op3JNOCHWnq8ERGjQmNj444dO+hn0dHRNTU1ahelPqvV+tZbb4nIECehJ06cULsol/D555+HhISIb4tWq928efMImTAZXF1dnW3FW3FToBMeO3ZxGmamkeT48eOFhYVxcXFr1qxRuxYXUlNT89FHHwUFBe3atUvtWlzLnj172traduzYMXXqVLVrcSFnz569ePHiokWLNm3apHYt6htxMQoAYF9atQsAAHBviFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAApiFEAACmIUQAAKYhRAAAp/wcWDQHt3J0JOQAAAJh6VFh0cmRraXRQS0wgcmRraXQgMjAyMi4wOS40AAB4nHu/b+09BiDgZYAARiDmBGIOIG5gFGRQAImxMSQAKSYmUmkOMM3IyA00lpGJgYmZgZmFgYWVgZWNgY2dgZ2DQYRBvAxqKRhwPnRbdmD1qlV2IE5QVuWB0NCr+0Dsh25qcHGgmv0wcRBAErdHUu+AZI4DQnwZXFwMACBcLMIR/09RAAAA9npUWHRNT0wgcmRraXQgMjAyMi4wOS40AAB4nI1SW2rEMAz89ynmAjGSHT/0uUmWUsom0Ka9Q2E/e38qb0m1C22I7AFZjIQ0skOz1+nl8wu/FibnANq5IoKPSETuguZgOD89zxjX07BFxuV9Xt8gqJqh55F5WpfLFmGMV3TJs2ROAeTjzdShm1luwIgu+iqFuUdHvuecc/mDGRsz+CSVSt4r2Tci+yBCse6WTMq06P8Vs/IOFSxKPNJiVd6hoQULjsh4nqeHBfysZFjmyVbSTjDd9YFo4rKiNwlZkUwnVmSTgxXFhmZFtdlYIfed3ffR3tunUt99A+QKgUKhWPbnAAAAbXpUWHRTTUlMRVMgcmRraXQgMjAyMi4wOS40AAB4nPN3hoIchRoNUz1DSzNDUx0DPWMw0LE21rOwNDc01NE10DMxNDMzM9exNtIztbQwsEBSZKhnZGlpYIKsyABJWheLvC6mKbpY7NLFcJFmDQD38yRXmgpBRAAAAABJRU5ErkJggg==\n"
          },
          "metadata": {},
          "execution_count": 4
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "changed_mol = Chem.DeleteSubstructs(mol,Chem.MolFromSmiles('Cl'))\n",
        "changed_mol"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "R6fmZ6_v6C-V",
        "outputId": "55678dfe-e71a-4f6b-da35-7acc6a811809"
      },
      "execution_count": 6,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c82a2430>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAP4ElEQVR4nO3de2yUVf7H8dMLpRdASotQbAGhIKUooeUOIpcSrvvHLo7rH9tsiDAxu2vXpAmzP6IpRrIMYd3UGIyzv2xcdgnJFo1akkUYEFoV5H6Rm3JruQgtttILaK9n/zizD+MA7dAzM888M+9X+gcdnj5zOsz5MM853+ecGCmlAAD0VKzZDQAAayNGAUALMQoAWohRANBCjAKAFmIUALQQowCghRgFAC3EKABoIUYBQAsxCgBaiFEA0EKMAoAWYhQAtBCjAKCFGAUALcQoAGghRgFACzEKAFqIUQDQQowCgBZiFAC0EKMAoIUYBQAtxCgAaCFGAUALMQoAWohRANBCjAKAFmIUALQQowCghRgFAC3EKABoIUYBQAsxCgBaiFEA0EKMAoAWYhQAtBCjAKCFGAUALcQoAGghRgFACzEKAFqIUQDQQowCgBZiFAC0EKMAoIUYBQAtxCgAaCFGAUALMQoAWohRANBCjAKAlgiP0cbGxuLi4lWrVtXW1prdFgCRKUZKaXYbgqKjo+PVV19977332tvbhRBxcXHFxcXr1683u10A9DQ0iD17xJUrorlZ9OkjMjPF7NliwIAHHHn+vDhzRsTHiyVLHno2t1vcvSuys0Vubs+bJCPR9u3bn3nmGfUL9u7dOzExUf15wYIF1dXVZrcOQI9UVcnCQhkfL4X42VdcnPz1r+WFC77Hr1snhZBJSV2dMytLCiGLi3XaFWkX9YcPH543b96iRYtOnjyZnp6+fPnypqamxsbGF198MTExcceOHTk5OWvWrGltbTW7pRazc+fON9988+OPPza7IYhWBw+KSZPEv/4l2ttFTo5YuVI4HGLlSpGbKzo6xL//LSZNEpWV5rRNJ4PDSlVVld1uj42NFUIMGDDA6XT++OOP3gfcuHGjsLBQ/dajR4/etWuXWU21ltOnT0+YMMF4w+Tl5d24ccPsRiHK1NTIQYOkEDItTW7b5vu3//mPTE+XQsjUVHnlyr3HQ/VpNBJi9NatWw6Ho3fv3kKIhISEoqKiH3744WEH79mzJycnRwgRExNTWFhYW1sbyqZaS2Nj4+uvv56QkKACtG/fvsYfPvroI7Nbh2iyfLkUQvbqJY8cefABR4/KhAQphHzhhXsPEqP+aG5udjqd/fr1E0LExsbabLbLly93+1Otra1Op1MNmKamppaWlnZ0dAS/sVbS2trqcrkGDx6scvOJJ5745JNPpJRbtmzp37+/enDp0qX+vNrw1tbW9vbbb2/cuLGtrc3stljHrVsyMVEKIX/3u64O++MfpRAyPl5eu+Z5hBjtmk8/LygoOHHixCOd4cKFCwsWLFA/PmPGjJMnTwapqdbS2dlZVlaWnZ2tXpmpU6dWVFT4HLBp06a0tDQhRFJSUklJyU8//WRWa63F6XQmJyerFzYlJWXz5s1mt8gitm71TCXt39/VYceOeQ77xz88jxCjXSgvLx81apR6O06ZMmXv3r06p8rMzBRCxMfHFxUVNTY2BrCdlrNr1678/Hz1wo4ZM6asrKyzs/OBR968ebOwsDAmJkYIMWrUKLfbHeKmWsuRI0dmz56tXti4uDg1gi+EePnll6P8LeeXVaukEDIhQXb9H3Z7u0xJkULI3//e8wgx+kBffvnljBkz1Fvwqaee6qKf++/27dtFRUVxcXHq6nXTpk0Baaq1nDp1ymazGZfwLpfLn6vOvXv3jh071hhorqmpCUFTraW6utqY+UxOTl68eHF9fX1NTc3UqVPVWy4jIyM633KP4Le/lULIrKzujxw9Wgohly3zfKtiNC5O/uIXD/1KSoqiGD19+rTRz4cMGeJnP/ff0aNHp0yZEoWjfleuXLHb7apL9+nTp6Sk5O7du/7/eGtra2lpaUpKihCif//+DDQbvv/+e++ZT7vd7jOfeeLEiWnTpqm33Ny5c8+dO2dWU8Pdr34lhZBjxnR/ZH6+FELOn+/5VsWoP18RH6M+/dzhcATpOqijo8MY9UtOTi4pKWlpaQnGE4WJuro6h8OhptpUP+/xx8mLFy8uWrRIJUJ+fv6hQ4cC21RruXPnjtPpfOyxx4yZz0uXLj3wSDXQnJ6ezkBzV37zGymEHDas+yNzcqQQ8pe/9HyrYjQhQW7d+tAvVSkVwTHq3c979eplt9tv3rwZ7CdV5aVq1C9Sy0tVP1dz7qqfX7x4Uf+05eXlWVlZ0TzQ3NbW5nK5MjIyjJnP48ePd/tTdXV1drtdveWys7N37NgRgqZaSVGRFEKmpMhuL3TS0qQQcsUKz7dRPjba0tJSWlqq+nlMTIzNZrtw/51ewRSp5aXqE/eQIUOMfn706NEAnr+5udnhcKhLhyFDhkTVqJ/b7R43bpx6YSdPnvzZZ5890o9XVFTk/u+2bpvNFoJPDJbx9797Lr27HveorvYc9s47nkeiNkY7OjrKysqGDx9u9PMjDyu4DbLIKy91u91PP/20emEnTZq0e/fuID3R8ePHp06dqp5oyZIlD7ukjRj79u179tln1e87evToHs98qoHmPn36GAPN7e3tAW+t9Zw/L2NipBBy3bquDvvrXz0x+vXXnkeiM0bdbrexpMjEiRPD4YL6woULCxcutHp56f79+2fNmqV+i2HDhm3atEm/wqFrUVJeeubMGWPmMz093el06o+nX7p0afHixeqceXl5Bw8eDEhTLaa6WhYXS+PFnDNHCiEHD5b19Q8+vrFRZmZKIeSMGfcejLYY/eqrr5577jn11hk6dKjL5Qqrj37WHfU7e/aszWZT426qn4cyziK4vPTq1avGzGdKSorD4WhoaAjg+cvLy4cOHaoGr+12e2BPHtbq66XD4alDMi7P9++XcXFSCDlvnry/9zU3ywULpBAyNlZWVt57PHpi9Ny5c0Y/T0tLC3E/9593eaklRv2uXbtmt9vj4+OD1M/9Z5SXqlE/q5eX1tfXOxyOpKSkYM983rlzp6SkRK1pEBXlpS0tsrRUpqZKIWRMjLTZfrb23V/+4rm0z8iQa9fKigp56pSsrJTr1nk+hwoh//znn50wGmL0+vXrPv389u3bJrbHH8eOHTNG/cK2vLSpqamkpMS7n5u+JlNklJe2tLS4XK6BAwcaM5/nz58P9pOeOHFi+vTp6i03Z86cs2fPBvsZTdDRIcvK5PDhnjQsKHjwEiT//KenPun+r7Q0+f77vsdHdow2NTU5nU61YpDq5999950pLekB7/JSNeoXPuWlqp8//vjjRj//9ttvzW7UPdYtL1Uzn08++aQx83n48OGQPXuEl5e63XL8eE8ajhv3gHXwvDU0yL/9Tb7wgpwyRY4YISdNkjabdLnkAxd1275drlwp//CHrk74pz/JlSvl1q06v0GoY1QtKTJo0CDjA11Y9XP/+ZSXmj7qp/r5iBEj1As7ffr0L774wtwmPYwx6meVgWa32z1+/Hj1wubm5paVlZnSjAgsLz1wQM6e7QnQoUOly9V9ZWhYCl2MqqWDRo4cqd6O06ZN+/zzz0P27EGyd+9e7/JSs0b93G63sbLy2LFjzern/lPlpWo8J5wHmg8ePDhnzhzvmU/TK5AqKyuN8tKlS5deMxaFs5Zz56TN5hnrTEuTTmc3y46EtxDFqNvtzsvLU//2OTk54d/P/WdueemhQ4fmzp2rXtisrKxw6Of+O378uHFTebiVl37zzTc+M58+mymYyNrlpbW1sqjIs59ScrJ0OGTYz4h0K+gxeujQoXnz5qmukpmZaa1+7j/v8tKJEyeGYNSvqqrKGFV44KYplhCG5aVqMwU1P56cnBy2M59Xr15dtmyZUV564MABs1vUnaYm6XTKvn09iyvb7dI6MyJdC2KMem+OlJqa6nQ6H2npICsKTXmp96Ypqp93sWmKJXiP+o0aNWrnzp2mNMNnM4XCwsLr16+b0hL/lZeXDxs2LNzLS1tbpcvl2UxJTcQbNxpFhKDEqM8SYUVFRfUPu/cg4gT1pvKebZpiFRUVFWaVl/rMfBYUFFjodrWwLi/t7JRlZXLkSE+ATpv2s/L4SBHgGFX93HuJsEjq5/4LeHmp/qYplnB/eWmwh4C63TTFKsKxvNTtlnl5ngDNyZERNCPiI2Axqvr5oy4RFsHuLy/t2aif6ueB2jTFEi5evGjcVJ6fnx+8m8rdbrfPpilBeqLQ8C4vTUxMNHGg+fDhw802mydAs7Lk++/LSJwRMQQmRr03R5o8efKePXsCctoIoFleGoxNU6zCp7w0sKN+Pds0xRLMLS81Nk35/5kzZWqqdDplpM+ISP0Y3bdv38yZM6Ozn/uvB+Wlwd40xRK8y0sDNeqn+rkavO7bt++jbppiFT7lpVevXg32M9bW1r7yyitqiDYxMfH/HI6HrsYUcbRi9Pnnn1f/ThkZGdHZz/3nf3lpyDZNsQrv8lKdPYsCuGmKJYSsvNT/TVMilVaMrl69mn7+SLzLS++/qdyUTVMsQfOm8iBtmmIJ165dM8pLJ0yYENjy0p5tmhJ5tGK0sbHx1q1bgWpK9DDKS2NjY/Py8m7cuFFXV7dw4cLk5GRh0qYpltCD8tL7N005duxYCJoaboJRXuqzaUo0z4iYv95odGpoaFixYoV6C8bGxqqbFKK5n/vP//JSn01THnVzpAgTwPLSQG2aEjGIUTNt2bJFDV0JIZKTkzds2GB2i6yh21G/0G+aYhWa5aXB2DQlAhCjJmtra9u4ceM777xjxWWMzeVdXtqvX7/NmzdLKbdt22asFhj6TVMsQQ00q8WnVXmpP6sxeG+aomZEwvTGUzMQo7C2rVu3Gp/o1ULgQojevXu/8cYbTU1NZrcufNXV1RUVFanRpOzs7E8//fRhR4Zs0xTrIkZheTU1Ncatt0KI3NzcryNr5YvgqaysNKaJ7i8vNWXTFCuKkVIKwPo++OCD06dPz5o1y1hoGf5oa2t79913X3vttebm5sTExKVLl27ZsiUuLu7DDz90OByXL18WQhQUFKxfv95YMhg+iFEAoqqqasWKFbt37xZCJCUlJSQkNDQ0CCEmTJiwfv36+fPnm93AsEaMAvBYvXr1hg0b2tvbhRADBw5cu3btSy+9pKaV0AViFMA9NTU1b731Vmdn55o1a4y5O3SNGAUALbFmNwAArI0YBQAtxCgAaCFGAUALMQoAWohRANBCjAKAFmIUALQQowCghRgFAC3EKABoIUYBQAsxCgBaiFEA0EKMAoAWYhQAtBCjAKCFGAUALcQoAGghRgFACzEKAFqIUQDQQowCgBZiFAC0EKMAoIUYBQAtxCgAaCFGAUALMQoAWohRANBCjAKAFmIUALQQowCghRgFAC3EKABoIUYBQAsxCgBaiFEA0EKMAoAWYhQAtBCjAKCFGAUALcQoAGghRgFACzEKAFqIUQDQQowCgBZiFAC0EKMAoIUYBQAtxCgAaCFGAUALMQoAWv4LFgpuYRkjbl8AAACkelRYdHJka2l0UEtMIHJka2l0IDIwMjIuMDkuNAAAeJx7v2/tPQYg4GWAAEYg5gBidiBuYGRjSACJMUNoJiZiaQ6IPkZuoHmMTAxMzAzMLAwsrAysbAxs7AwiDOJZUKvAgMP26MQD0h9q94E4hRL+B6YGvbMDscNX/NpvLbYELN7Tobw/w+SkHZRtD2Tvg6qxB6qxg+p1AOoFiwPNdACaCRYXAwBRKianxvaxhQAAAPx6VFh0TU9MIHJka2l0IDIwMjIuMDkuNAAAeJx9kk2KwzAMhfc+hS4QI8t/8rJJSimlCXQyc4fZz/0ZqSV2C8FyHtjOZyvvEQNaj/n2+we1aDYGADtPKQV+PCKaO+gExvPlusC0ncZ9Z1q/l+0LGLKckPFJnrb1vu84mGAINsboKcOAlkIunAAtPqudJSW9JR85ybU2pBhyOAC9gs6WyJRIr/SEgemADEqiTZ4Z9b0vVDIfgFHAyg0dMAlYW3c6Z+Gql6FjhmGFGk8nnfMyf+T6Snpcl7klrYNanLIA30JzotCScaLY/DtRai6dKDczTsTvn/LeWNf7zyFz8w/q8XcMacrCuwAAAHl6VFh0U01JTEVTIHJka2l0IDIwMjIuMDkuNAAAeJxFjkkOgCAMRa/iUhPa1A7QhiUH8EIcXjQB/vLlT639eo5+gqKZiSQgZC3hOVUQZDEviVCzadGBbgxzzp9NmNR5MMIs7sTDJ8FRPNWFYLMZ3cnZD3tg3lgvrv4CLckh52rptbkAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 6
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "mol = Chem.MolFromSmiles('c1ccccc1Cl')\n",
        "mol"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "UQLrVpDpOMZ_",
        "outputId": "2e186c53-6689-4643-c9b8-fe9802e80ece"
      },
      "execution_count": 7,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c82a25f0>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAZSElEQVR4nO3de1SUZR4H8N87DIPSgAgSkqgpGgi4KKitKScvpIWdXDL0aKQiXtYjwUAsCHgjUPDCRZH1QpiZVAtRHDdhRVwjb+XKxWBK4pKBF0gHFOUy9/3jbae2y8sMzDzPMPP7HP+Y7Mfpmwe/vDPv+zwPo1arASGEUH/xaAdACKHBDWsUIYQGBGsUIYQGBGsUIYQGhE87AOpDeXn5nTt32NcuLi62trZ08+hdU1NTb28v+3rBggUODg508yCkK6xRoxYQEFBSUkI7BTlCobCkpGT27Nm0gyCkAwYfeDJaRUVFgYGBADB06FD2d5ydnYVCIdVQ+tfc3CyVSgFAJpMplUoXF5eWlhbaoRDSAdao8RoxYoREIvH29q6urqadhYSKiorp06er1erU1NTY2FjacRDSFt5iMlIpKSkSiYRhmI8//ph2FkJ8fX3nzZsHAElJSQqFgnYchLSFNWqMZDLZu+++CwBr1qyZMGEC7TjkFBQU2NnZdXV15ebm0s6CkLawRo1RWlpafX29h4fH4cOHaWchavjw4UePHgWA+Pj49vZ22nEQ0grWqNFpa2tLTU0FgIyMDD7f7B6lCAoKmjNnTnt7e3JyMu0sCGkFbzEZndWrV7/33nuBgYGffPIJ7Sx0XL9+3dfXl2GY6upqT09P2nEQ6gNejRqXysrK999/XyAQ7N69m3YWary9vdeuXatQKEQiEe0sCPUNa9SIqNXqiIgIlUoVFRU1ceJE2nFoSkpKsrOzKysrKy4upp0FoT5gjRqRvLy8ixcvOjk5xcXF0c5CmaOj49atWwEgIiKCfTgfIaOFNWosuru7ExISACA1NdX0Fs73Q3h4uKenZ0NDQ3Z2Nu0sCHHBGjUWKSkpzc3NPj4+K1eu5BiTyWSvv/761atXiQUzhJ6enuTk5Ly8PI4ZPp+fkZEBAImJia2traSiIaQ7NTICzc3N1tbWDMNcuHCBezItLQ0AJk2apFQqyWQzhPz8fAAYOXLkw4cPuScDAgIAYP369WSCIdQPWKNGISgoCABef/117rG2trZhw4YBQHFxMZlgBqJSqfz8/AAgJiaGe7K+vl4gEPB4vGvXrpHJhpCusEbpu3jxIsMwQ4cO/eGHH7gn161bBwCLFi0iE8ygKisreTyeQCCoq6vjnoyMjASAWbNmqVQqMtkQ0gnWKGVKpdLX1xcA3n77be7JqqoqCwsLS0vLGzdukMlmaCEhIQCwePFi7rGHDx+OHDkSAPLz88kEQ0gnWKOUHTlyBABGjx7d1dXFPfn8888DwFtvvUUmGAGtra3sZxT/+te/uCfZvQW0+VNCiDysUZq0v8766KOPAODJJ5/s6Oggk40MdveASZMmyWQyjjHNNXtSUhKxbAhpCWuUpqioKG0+9evu7n766acB4OjRo8SykSGVSp955hkAyMrK4p68cOECwzDW1tZ9foKMEGFYo9TU19dbWVnxeLz//Oc/3JOJiYkAMGXKFIVCQSYbSUVFRQAwfPjwe/fucU++9tprABAcHEwmGEJawhqlRssnIm/duvXEE08AwOeff04mGHkLFy4EgLCwMO4x7Z+uRYgkrFE6zp49CwC2trZ3797lnlyxYgUALF26lEwwKsRisaWlpYWFxddff809uWXLFgDw8fEZ1KsPkInBGqVALpez22impaVxT16+fJl9pPT7778nEo2asLAwAJg3bx73WFdX15gxYwDg+PHjZIIh1CesUQrS09MBYMKECb29vRxjSqVyxowZALBt2zZi2Whpb28fMWIEABQVFXFPnjhxAgCcnJz6XEiKEBlYo6RJJBJ7e3sA+Oyzz7gnjx07BgCjRo16/PgxmWx0ZWVlAcD48eO5f7qoVKrZs2cDQFxcHLFsCHHAGiVtw4YNAODv78891tnZ6ezsDAAnT54kE4w6hUIxefJkAEhNTeWevHbtGruQ9LvvviOTDSEOWKNE1dbW8vl8Pp9fW1vLPRkbGwsAM2fONKuF5OfOnQMAGxubO3fucE+uWrUKAAIDA8kEQ4gD1ihRL7zwAgCIRCLuscbGRvaR0q+++opMMOOxePFiAAgJCeEea21tZTe3PnPmDJlgCP0RrFFyCgoKAMDe3v7+/fvck1pWiUnS/kfIzp07AcDDw0Mul5PJhtDvwholpKenZ9y4cQBw6NAh7kn2ja1QKOzzja2piomJ0eYDDalUyh78l52dTSwbQr+FNUpIUlISAHh6enJfOmlus6SkpBDLZmw0t9fy8vK4JwsLC7W8wEfIcLBGSbh165ZQKASA0tJS7knNQz89PT1kshmn3NxcLR/2WrBgAQCEh4eTCYbQb2GNkvDGG28AwJIlS7jHNI+gf/rpp2SCGS3tlx6IxWL24Yeamhoy2RD6FaxRg7ty5QrDMFZWVvX19dyTWi6INBPaL4TduHEjAMyfP59ILoR+DWvUsFQq1bPPPgsACQkJ3JPffPONlttzmI/ly5cDwLJly7jHJBKJg4MDAJw6dYpMMIR+CWvUsI4fP67lScLsZnGbNm0iE2xQaGlp0XKTwMzMTABwdXXlXkiKkCFgjRrQo0ePnnrqKQA4ceIE96T2Wxebmx07dmizZbVcLvfy8gKAvXv3EsuGEAtr1IDi4+MBwNfXl3tzTM1BGgcOHCCWbbDQ/gCVsrIyLReSIqRfWKOG0tTUNGTIEG22atfyWDezpf1xfi+//DIArF27lkwwhFhYo4by6quvAsDKlSu5x9ra2rQ8ZNicaXm4dENDA7uQ9OrVq2SCIaTGGjWQf//73+yCztu3b3NPrlmzBgBeeeUVMsEGqaqqKgsLC4FAcOPGDe7J6OhoAHjuuefMamcsRBfWqP4pFIo//elPALBz507uycrKSnbfzLq6OjLZBq+1a9cCwMsvv8w9pllI+tFHH5EJhhDWqP79/e9/B4Bx48ZxL+hUqVR+fn4AEBMTQyzb4KX59KO4uJh7MicnBwBcXFzM5NQARB3WqJ5pFnR+/PHH3JN5eXnsnZMHDx6QyTbY7d27FwDc3d2578Uplcpp06YBwI4dO4hlQ+YMa1TPIiIiAGDu3LncY93d3WPHjgWA3NxcMsFMgFQqdXNzA4DMzEzuyUuXLrELSW/evEkmGzJnWKP69O2337ILOqurq7knt2/fDgBTp07F89Z18s9//lPLdQrLli0DgOXLl5MJhswZ1qg+vfTSSwDw17/+lXtMs8axvLycTDBT8uKLLwLAxo0bucdaWlqsra0Zhvniiy/IBENmC2tUb9gLJTs7ux9//JF7Ei+UBkKzh0ufl/zbtm3DS35EANaofmg+tsvIyOCexI/tBg4/gEZGBWtUP/AmMkmaxyEKCwu5J/FxCEQA1qge4CON5OHDuch4YI3qwbp16wBg0aJF3GOaBTYffvghmWAmDJeKIeOBNTpQ7HJvS0tLXO5NmPYbF4SEhADA4sWLyQRD5gZrdKBw8yGKtNxGq7W1FbfRQoaDNToguBUmXZpNXb/88kvuSdzUFRkO1mj/4cbsxiAuLk6nIwaysrKIZUNmAmu0/xITE3U6JmjPnj3EspkPPPBKv8bWjBVWCR2uO9AOMpgwarUakO5u377t5ubW1dX1+eefsx+P/pH9+/eLRCJXV1exWGxlZUUsofl47733Vq9ePXLkyLq6OltbW47JF1988cyZM2FhYVlZWcTiGQlxr/jC4wut8laJQmJjYePEd5r5xMzpT0xngPnl2IivR0gUEgEjkE6V0oo6+NDu8cFqxYoVALB06VLuMTxCnQCVSvXss88CQEJCAvekWCxmF5J+/fXXZLJRJ1PJjtw78nTN01ABv/01pmbMCckJlfrnR0ccrjtABQgqBRQzDzpYo/1x+fJldkHn999/zz25ceNGAJg/fz6RXObrypUrDMNYWVnV19dzT4aFhQHAvHnzyASj60f5j3++8effLdBf/lrSuESm+unOG9ZoP2CN6kypVM6YMQMAtm3bxj0pFov5fD6fz6+pqSGTzZwFBwcDwGuvvcY9pllIWlRURCYYLZ2KTnexu6YrF9YvLOgouCO7o1ApOhWdVx5fCW8Jt6y0ZP9t7K1Y9quwRvsBa1Rnx44dA4BRo0b1uaBzwYIFABAeHk4mmJm7deuWUCgEgNLSUu5J9oPR8ePH9/b2kslGxeqbq9mKtKi0OHb/2O/OVHVXPXn9yYCGgAeKn/YcwBrtB6xR3WgWdJ48eZJ7srCwEADs7e3v379PJhtKSkoCAE9PT7lczjGmUCgmT54MAKmpqcSyEVbbU6u5Dt3Xuo9j8qb0pkL186MmWKP9gDWqm9jYWACYOXMm94JOqVQ6ceJEAMjOziaWDfX09IwbNw4ADh06xD157tw5MOknecOaw9gO9RB7/LIl+4Q12g9Yozo4f/48j8fj8XhfffUV9+SuXbsAwMPDg/uyCOldQUEB+yZAIpFwTy5evBgAQkJCyAQjTHNrPq0tTacvxBrtB6xRHbBv5ydOnMg91trayj69eObMGTLB0C/5+/sDQGRkJPdYY2Mju8tBnz8UB502eZvmHb24R6zT12KN9gOP1POpg15aWtrdu3cZhsnPz+eejI2N7ezsDAwMZG8xIcIyMzP5fH5WVpZYLOYYGz9+fEREhEqlEolEatNahHKj9wb7YghvyDNWz9ANYw74tAMMGunp6QDg7+8/ZcoUjrHHjx9fuHDBysqK3Q8fkefp6RkaGpqTk1NWVubp6ckxmZCQkJubW1VVxeMZ/HpiytUp1RbVhv6vyKbKLBnLDmUH+48OFg58Bv+OGxz+EeuZUCgUi8VXrlxxdXWlncV8JSUlhYaGTp8+XZthhmH6HhpUOpWd7AuhhZBuEjOBNaqtqKio6OjosrKyqqqqqVOnckwOGTJk7ty5xIKh33J0dHR0dOxzLDk5WSKRzJw5kz1nkEAwMgSMgH0hV8vpJjET+Nmott566y1nZ2f1/45HRoNdY2PjgQMHeDze/v37TalDAcDOwo59obksRQaFNaqDDz74wMLCoqGh4erVq7SzoIGKioqSSqWrV6/W8r3/IOJk6cS+kCgkms9JkeFgjepgzpw5f/vb39RqdUREhInd2zU3586dO3XqlI2NTXJyMu0s+uc5xHMobygAqEFd0V1BO47pwxrVTUJCwlNPPfXll1+yB6CjwUihUERGRgLAli1b2GeBTYwlYznNehr7uqCjgG4Yc4A1qhuhULhz504A2Lx5c1dXF+04qD8OHTpUU1Pj6uoaERFBO4uhBNsHsy8+7PjwlvwW3TAmD2tUZ6tWrZoxY8bt27d3795NOwvSWUdHB3v6S3p6ugkfRrDCfoUD3wEAHikfhf4QKlXjVvYGhDWqM4ZhMjMzGYbZu3fvzZs3acdButm6datEIpk3b94rr7xCO4sBCXnCQ6MPsa9LO0sXNSxqkjb9dkyulv+j4x/v3H+HbDpTg2cx9VNwcHBeXl5QUFCfa0OR8fjmm2+8vb3VanVVVRW7V55p23x78+62n94zWTFWC20XzrWZ62zpbMVY3Zbfvtp19XTnaYlC4mLp8sPkH3jAAzyLqV+wRvtJc6Td+fPn58yZQzsO0srChQtLS0vffPPNAwcO0M5CyKF7hzbf2cz9AOkwi2HX3K9NsJoAWKP9gm/q+2nUqFHs3qMikUipVNKOg/pWVFRUWlo6fPjw7du3085CzkbHjY2ejVudt/pY+/zqEFBbC9u5NnP3u+xv8mxiOxQAhDyhkCe0sbChEXawwqvR/uvt7Z00adLNmzePHDmyfv162nEQF5lM5uXlVV9ff/DgwU2bNtGOQ0eXquuu/O5D5UMLxsLBwsFF4PKrYkX9gzU6IPn5+cuWLXN0dPzuu+/s7Oxox0F/KDU1NS4uzsPDo7q62tLSknYcZFLwTf2ALF269Pnnn7937x57ChAyTm1tbSkpKQCQnp6OHYr0Dq9GB6q6unratGk8Hq+mpsbNzY12HPQ7QkJCjh8//pe//OXTTz+lnQWZILwaHagpU6aEhobK5fKoqCjaWdDvqKysPHHihEAg2LNnD+0syDRhjepBcnKynZ1dcXFxSUkJ7Szo/6jVapFIpFKpIiMj2bNaEdI7rFE9cHR03LJlCwCEh4fLZDLacdDPPvjggwsXLjg5OcXFxdHOgkwW1qh+hIeHu7m5NTQ0ZGdn086CftLT0xMfHw8AKSkpw4YNox0HmSysUf2wtLRMS0sDgB07drS2ttKOgwAAUlJSmpubfXx8Vq1aRTsLMmVYo3qzaNGil156qbOzc8eOHbSzIGhpaUlLS2P3kSFw8CcyZ/jtpU/sY4k5OTkVFbjlOGXR0dHd3d3Lly/38/OjnQWZOKxRfXJ3d9+0aZNKpRKJRPhALkWXLl0qKCgYOnTorl27aGdBpg9rVM+2b9/u6Oh48eLFwsJC2lnMlObH2ObNm8eOHUs7DjJ9WKN6Zmdn9/bbbwNAVFRUd3c37TjmKDc399q1a6NHj46OjqadBZkFrFH9W7dunbe3d0tLS0ZGBu0sZufRo0fbtm0DgH379llbW9OOg8wC1qj+WVhYsAW6a9eu5uZm2nHMS2JiYmtr63PPPRcUFEQ7CzIXWKMGMXfu3CVLlnR3d7OrmxAZDQ0NBw8e5PF4+/fvZxjcSRMRgjVqKOnp6dbW1idPnrx48SLtLOZCJBJJpdLQ0NBp06bRzoLMCNaooYwZMyYyMlKzNQbtOKavrKzs9OnTNjY27C0+hIjBGjWg+Pj40aNHV1RUnDx5knYWE6dQKCIjIwFg+/btI0eOpB0HmResUQOytrZOTk4GgJiYmM5OrqMZ0QAdPHiwtrbW1dU1LCyMdhZkdrBGDeuNN96YNWtWW1vb7t27aWcxWe3t7eyPq8zMTCsrK9pxkNnBQ0QMrqKiYsaMGXw+v7a2FncONoSNGzcePnzY39//7NmztLMgc4RXowbn6+sbHBwsk8k2b95MO4sJEovF77zzDp/Px8UOiBasURL27Nlja2v7ySef4OWS3kVGRioUik2bNnl5edHOgswU1igJTk5OMTEx8L+/87TjmI7CwsKzZ8/a29tv3bqVdhZkvrBGCYmOjp4wYYJYLM7JyaGdxURIpVL2c5Lk5GQHBwfacZD5wholxMrKKjU1FQC2bNkikUhoxzEF+/bta2ho8PT0XLduHe0syKxhjZKzZMmSF154ob29PSkpiXaWQa+trY09dz4jI4PP59OOg8waPvBElFgsnjJlCgBUVVXhLZGBWLVq1YkTJ1599VXcHhtRh1ejRHl6eq5du1azchH1D7u+ViAQsJ+TIEQX1ihpO3fudHBwKCsr++yzz2hnGZTUanVERIRKpYqOjsblDMgYYI2SZm9vz25Cyu7qRjvO4PP+++9funTJyckpNjaWdhaEALBGqQgLC/Py8mpsbDx48CDtLIOMZidsdkUD7TgIAWCNUqFZucieeEE7zmCya9eulpYWdn0t7SwI/QRrlA5/f/9FixZpzl9D2mhubs7IyGAYJjMzk8fDb11kLPB7kRp2Vzf2NGDaWQYH9szq4ODg2bNn086C0M+wRqmZMGFCWFiYSqWKiIjAp3f7dP78+cLCQs1O2AgZD6xRmtgTLy5fvlxQUEA7i1FTKpXsk7bx8fFjxoyhHQeh/4M1SpPm/LXo6Oju7m7acYxXTk7O9evXR48ejcsWkBHCGqWMPQ24paVl3759tLMYqQcPHrA34tgzq2nHQejXcE09fZcuXfLz8xsyZMi33347duxY2nGMTmRkZGZm5uzZs7/44guGYWjHQejX8GqUvlmzZgUFBfX09MTHx9POYnRu3LiRnZ3N4/EyMzOxQ5FxwqtRo9DS0uLu7t7T01NeXu7n59fnfGJiYlNTE4FghhMZGcludsUtICCgpKRkw4YNhw8fJpAKof5QI+PAHoPh4+OjVCr7HJ4+fTrtb5yBOnXqVJ//m+zuLba2tnfv3tXHnzFCBoFXo8aip6fH3d29ubn52LFjISEh3MPFxcX37t0jE8xA5s+f7+LiwjEgl8snT55cV1eXnp6ON+iRMcMaNSJ5eXnBwcFOTk51dXXDhg2jHYeytLQ0zQFWAoGAdhyE/hDeYjIiK1as8PPza2trS0lJoZ2Fsnv37rGrlQ4cOIAdiowcXo0al8rKyunTp/P5/NraWnPek3jDhg1Hjx4NCAg4ffo07SwI9QGvRo2Lj4/PypUrZTIZe669eaqurs7NzbW0tExPT6edBaG+YY0andTUVFtb26KiojNnztDOQodIJFIqlW+++aabmxvtLAj1DWvU6Dg5OcXFxQFAVFSUXC6nHYe0/Pz88vJyR0dH9gkwhIwffjZqjGQymZeXV319/Zo1a3Jzc2nHIaejo2P8+PEPHjw4cuTI+vXracdBSCt4NWqMBALBmjVrAODdd99taGigHYecoKCgBw8eCIXC0NBQ2lkQ0hZejRqvESNGSCQSb2/v6upq2llIqKiomDZtGgBkZGSIRCLacRDSFtao8SoqKgoMDASAoUOHsr/j7OwsFAqphtK/5uZm9qBpmUymVCpdXFxaWlpoh0JIB1ijRo3dmIN2CnKEQmFJSQketYQGF6xRY1deXn7nzh32tYuLi+kdzt7U1NTb28u+XrBggYODA908COkKaxQhhAYE79QjhNCAYI0ihNCAYI0ihNCAYI0ihNCA/Bdfnta0SJB5rQAAAJh6VFh0cmRraXRQS0wgcmRraXQgMjAyMi4wOS40AAB4nHu/b+09BiDgZYAARiBmh+IGRjaHDCDNzEwkQwPEYBFkUAAZxM3AmMHEyJTAxJzBxMySwMKawcTKxsDKkCDCCKRYmJkYxeOgFoIBu+ekJXYP3ZbtB3Hybk/bj2B/A9INqkji9iA2VD2YPUmlE0g3LIWIqzjADBUDALtqJMCNhFzYAAAA4HpUWHRNT0wgcmRraXQgMjAyMi4wOS40AAB4nH2RXQrDIAzH3z1FLjCJH631sdUyxqiFrdsdBnvc/VmsdGkZmBjQ+FOTvwKy3eL19YGf6SgEAFaG9x6eBhHFBHkCw3i+JAhLP2yZMD/ScgdHjtmPZL/M05ZREACl0cpqDycltfdoOkrhanxUE0j7ynWNU1XSFNLrleT9P9Aerqzc2OyLrHAtcUqiK81UXnYQ3qBlU8gKOKZ4kKsIOMwpsoCWQrNMlsKwFjYHd5y94b5oAS2XrygcF2lzal/K/uG83v6c5uILV+5vM9XVIYUAAABmelRYdFNNSUxFUyByZGtpdCAyMDIyLjA5LjQAAHicc85JNkwGAUOFGg0jPVNzQxNjHQMda0M9AxjTQM/YyNDEyFLHUM/I0tLARMda11DP0NzC1BxFxNIIJGKAJKsLl4YbARfSrAEAe9kZ0KCtudMAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 7
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "combinemol = Chem.CombineMols(changed_mol,mol)\n",
        "combinemol"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "ZZsbYyP0OfN6",
        "outputId": "18f35f07-2379-4199-9969-4d291374a513"
      },
      "execution_count": 8,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c82a2580>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAbUUlEQVR4nO3daVxTV94H8H/CvgoIomwuIKuiiFsBdaqodR21xbaPpVatodUWgfZjsLVitaNpRcCxOGKHVkfbadFhRh11NFhbxV1wwwoCLmARQSAGCYQs53lxNY0xIGY7Ifl/P30xnoTLD0d/3nvPuScsQggghBDSFJt2AIQQ6t6wRhFCSCtYowghpBWsUYQQ0grWqN4JBAKZTEY7BUJIX7BG9Wjv3r1eXl6urq7Ozs5///vfacdB6MW1tEBjI+B6nk5hjepFdXX1ggULZs+efe/ePQAQiUSLFy9OSEhoamqiHQ2hLmhogI8+Al9fcHSEnj3Bzg4mT4Zjx556T79+MHq0mq9dtAgcHKCuzjBJjQHWqI41NTWlpqYGBQVt377d0tIyPDz81KlTY8aMsbKy2rZtW0BAwLZt20xvra5IJHrjjTd69erl7e29d+9e2nGQdu7ehVGjIDMThg+HLVtgxw5ITobiYpgwAbKz/3hbSwu0tqr58rY2EInM6wSWIB0Ri8U5OTkeHh4AwGKx4uLiKioqFK/euHFj4sSJzO/52LFjS0pKKEbVIblcnpeXx/zUCh9++GFLSwvtaEhTsbGExSI//PDUYG0tCQwklpbk0qXHI+7uJDxczZf/3/8RAFJbq/ecRgNrVAdkMlleXl7//v2ZEomNjS0qKlL7zry8PE9PTwCwsrJKTExsbm42cFTdOnnyZHR0NPNT29vbv//+++Hh4SwWCwC8vb13795NO6BeHDlyJCoqaty4cSdOnKCdRQ+KiwkAmT1bzUt8PgEg77zz+JdYo09gjWqLz+cPGTKEqZKwsLC8vLzO39/U1JSYmGhhYQEA/fv3P3DggGFy6ta1a9fi4uKYn9rDwyMrK0sikTAvFRcXjxw5knlp+vTpt2/fphtVh+rr6+Pj45XPuxcuXNjW1kY7l0599RUBILm5al6Sy0nPnsTH5/EvsUafwBrV3NmzZ//0pz8xf538/PxycnKkUmkXv7aoqGjEiBGKrrlz545eo+pQVVUVh8Nh/hlwdHTkcrlCoVDlPTKZLCcnx9nZmTlLTUtLE4vFVNLqilAoXLVqlaOjI3PHJiAgwNfXV/Fv56+//ko7oO4kJBAAUlCg/tWoKAJAmH853N1Jr17kk09U/wsLwxpFz1dWVhYXF8dcvfbs2ZPH47W2tr7oQSQSSVZWVjfqmoaGBi6Xa2try9yU4HA4tZ3+Vbl3757i3C0oKOjo0aMGi6pD7e3tOTk5vXv3VtyxuXDhAvPSgQMHQkNDmWKNj4+/f/8+3ai6ER9PAMjFi+pfnTiRAJCGBkIIcXcnFhbEzU31P2trrFHUmfr6ei6Xa21tzXQfl8sVCATaHLCmpkbRNeHh4SdPntRVVB0Si8VZWVkuLi5qZ8869/PPPwcHByu6pq6uTq9RdYiZPQsICGD+3xk9evSzZ53t7e08Ho/5p8XV1TUrK0smk1FJqzNLlxIAcuiQ+ldHjiQsFmFu4OBF/RNYo1316NEjHo/HnDyy2ez4+Pjff/9dVwc/evRoUFCQomvq6+t1dWQtMbNn/fr1e+7sWSfEYnG365rCwkLF7FlwcHBeXp5cLu/ozZWVlVOmTGHeHBkZef78eUNG1bHNmwkA2bJFzUsSCXFyIoGBj3+JNfoE1ujzMZd1zAw7UyVXrlzR+XcRiURpaWk2NjYA4ObmlpOT08nfW8Pg8/nh4eHMTz18+PCCju6XdU15efnkyZOZo8XExFy9elVXOXWrpKREMXvm7e2dk5OjmD3r3L59+5gbppaWlomJic/eMjZqNTXk+HFCCCkvJywWefllNe/Zs4cAkKSkx7/EGn0Ca7Qzz17WHWf+qOlNeXn5pEmTmG83ZswYWstLz5w5M27cOOXZM12dP+7bt8/Hx0fRNUa15Etl9iwtLU0kEr3QER49esTlcpkjeHl57dixQ09Rdam5mfB4xMmJ9OlDHj0ihJDXXycAJD39qbeVlxMfH+LoSG7dejyCNfoE1miH+Hx+ZGSk8mWdwb41xeWlpaWlKrNnOl/QIxAIFEu+vL29Dfkb2xHl2TNra2sOh6PNfNHFixdHP3lKctq0aTdv3tRhVF0Si8mmTcTDgwAQFovMmUPu3SOEkAcPSGQkASDR0WTNGpKZSd59lzg4EBsb8q9//fHlWKNPYI2qofFlnQ4pLy/18fHZs2ePvr9jXV1dYmKipaUlADg4OGg/e9a54uLiUaNGKZZ83VKc4xhWS0sLj8djZs/YbHZcXFxlZaX2h5XL5Tt27OjZsycA2NnZpaWlGdfyUrmc5OURf38CQADI6NFE5TJLJCJffkkiIoitLQEgXl4kPp6oXBvFxJA5c9QcfPlyEhpKHjzQY34jgzX6lDt37igu65ycnDS4rNMtleWlelrK3tzczOPxnJycmGttDodTU1Ojj2+kglle2qNHDypLvmQy2Y4dO7y8vBS3vIuLi3X7LWpra+Pj45lT+4EDB/L5fN0eX0N8Phk27HGBhoQQI7ga6O6wRh/T7WWdDul1KbvK7Nn06dPLysp0dfAuUl5eGhgYqOVEVhepzJ7pdU3rL7/8wiwvBYC4uDiaf65KSkhc3OMC9fEhOTmkyw+MoE5gjerrsk63dL68lJk98/f3Z4750ksv0X1C/NixY4ZZXnr69OmxY8cyP3Xfvn137NhhgBUR7e3tWVlZDg4OAODi4kJhydedO4TDIRYWBIC4uhIej1C9zDIxZl2jz17WXezo4Q3joLKUXePlpXw+f9iwYcxPHRISYgyTPET/S9mvX7+umD1zd3fXx+xZ5+gsL33wgHC5j29xWlsTDod0nycgugvzrVE+nz948GDmz/SIESN+/vln2om6RMvlpVevXp02bRrzU/v4+LzQPgCGUVFRofPlpXfv3uVwOMqzZw8fPtT+sJox3PLSlhbC45EePQgAYbNJXBwx2jUD3Zw51iiVyzrd0mB5KTN7xmazmXM9Ho9Hd/ascyrLSzXumubm5rS0NDs7O3iyD8A9Zk0PVfpeXiqRSLZu3XrrlVce3wadNo3o4YERpGBeNcpc1jHtQ+WyTre6uLz0wYMHXC6XOYG1trZOTExsbGw0cFQNaLm8lNlFu1evXvBkH4AbN27oKapmLl26pI/lpYrLrDE+PrLoaPLLLzo5LOqEGdXonj17FCuZPv/8c6N6fkZjystLvb29VZaXMvsAMCuKmNkz410K3gENlpcy+wAMGDCA+aqoqKjCwkL9J9WEbpeXnj59esyYMYo1D53vA4B0yIxqVCgUent7G8llnW4VFRUpdkqOioo6ffq0SCT66KOP+vTpo5g9u6T47IfuhpkJdHNz68qSLz6fHxERwfzUoaGhRjJ71jntl5f+9ttvKpdZRr7pookxoxolhJjwBwRJpdKsrCxmCT1zS5H5HwbYB8Awnru89Ny5c+PHj2fe4Ovra4SzZ53TbHlpdXW14oER6rNnZsu8atTkXb16tW/fvsxfRWtr67Vr15rYZd2xY8dCQkKY251z584tLS0lhJw7d05xNufm5qbZLtrG4IWWlzY2NnK5XGObPTNPWKMmaM+ePdu2beumVfJcra2taWlpzPJSFouluHFhb2+fmpra1NREO6C2nru89NnPoC0vL6cSFTGwRlG3dP36dWZmhhEeHk5rcxM92bdvn5+fn2LJF/N5LSqfQTthwgTFJ5ogiliEEECoe1qxYkVdXd2sWbNmzJhBO4vutbS0rF27duPGjVKp1MLCYvz48fX19ZcuXQKAsLCwtLQ0xbQSogtrFCGjdvny5fHjxzc2NjK/9PX1Xbly5aJFi5hpJWQMsEYRMnZSqfTtt98uLS0NDQ3Nzc1lnqRAxgNrFCGEtMKmHQAhhLo3rFGEENIK1ihCCGkFaxQhhLSCNYoQQlrBGkUIIa1gjSKEkFawRhFCSCtYowghpBWsUYQQ0grWKEIIaQVrFCGEtII1ihBCWsEaRQghrWCNIoSQVrBGEUJIK1ijCCGkFaxRhBDSCtYoQghpBWsUIYS0gjWKEEJawRpFCCGtYI0ihJBWsEYRQkgrWKMIIaQVrFGEENKKJe0AyHi1tLQIBAILCwsnJycqAVpbW8VisbW1tYeHB5UACHUF1ihSb/PmzampqSKRiHYQAIDXXntt9+7dtFMgpB6LEEI7AzI6QqHQ3d1dIpEAAJvNtrW1pRJDIpFIpVLmj2hubu7ChQupxECoc1ijSI2kpKRNmzbZ2NgIBAJaHaoQHR196tQpPz+/W7dusdl4Nx8ZHfxDiVTdv3//u+++A4AtW7ZQ71AAyM/P79WrV1VV1a5du2hnQUgNrFGkisvlCoXC2bNnG8lFtKen54YNGwBg+fLlQqGQdhyEVOFFPXpKcXHxiBEjLC0tS0pKBg4cSDvOY4SQl1566ezZs59++ukXX3xBOw5CT8GzUfQHQsiyZcvkcnlKSorxdCgAsFis7OxsNpudnp5eUVFBOw5CT8EaRX/4/vvvCwsLPT09V6xYoTyenZ1dUFBgyCTV1dUff/wxs1SAERkZOW/ePLFYnJqaasgkCD0fQYgQQkhLS4ufnx8AfPfdd8rjFRUVNjY2bDa7rKzMMEnkcvmQIUMAYNOmTcrj9+7dc3Z2BoAjR44YJglCXYE1ih5buXIlAAwbNkwmkymP//nPfwaABQsWGDLMvn37AMDV1bW+vl55nLkxGhYWJpFIDJkHoU5gjSJCCKmqqrK3t2exWCdOnFAeP3r0KAA4OjrW1NQYONIrr7wCAEuWLFEebGtrCwgIAIC//e1vBs6DUEewRhEhhMTFxQHAvHnzlAelUungwYMBYP369YaP9Ntvv1lZWVlYWFy+fFl5nHkq1M3N7cGDB4ZPhdCzsEYRKSwsZLFYdnZ2d+7cUR7fvHkzAAwYMKC1tZVKsMTERAB4+eWXVcYnTpwIAElJSVRSIaQCa9TcyWSyyMhIAFizZo3yeGNjo7u7OwD8+9//ppVNkSE/P195vKSkxNLS0tLS8urVq7SyIaSANWrucnJyAMDX17elpUV5/IMPPgCA8ePH0wrGyM7OVntG/N577wFAbGwsrWAIKWCNmrWHDx/27t0bAPLy8pTHFfclr1y5QisbQyqVhoeHA8C6deuUxxsaGtzc3ADgv//9L61sCDGwRs1aSkoKAERHR8vlcuXxyZMnA8DSpUtpBVOmWC3w+++/K49nZGQAQEBAQFtbG61sCBGsUXNWXl7OrKs/f/688vh//vMftWs2KZo9ezYAzJ8/X3lQIpGEhYUBQHp6OqVcCBGCNWrOpk6dCgAcDkd5UCwWBwYGAsBf//pXWsGeVVlZaWtry2Kxzp49qzzO5/MBwMnJ6d69e7SyIYQ1aqaYAnJ2dlYpoC+//BIAQkJC2tvbaWVTi3mUfvTo0Sr3H6ZNmwYAixcvphUMIaxRc6S4HN64caPy+P3793v06AEA//vf/2hl60hzc7OXlxcA7Ny5U3m8o1sTCBkM1qg56mhyhtmneebMmbSCdY7Zk9/b27u5uVl5vKOJMoQMA2vU7HS0VKi4uJjNZltbWxtsJ6cXJZfLR44cCQArV65UHlcs2/rpp59oZUPmDGvU7CQkJDy7cF0ul48ZMwYAli9fTitYV5w+fZrFYtna2t66dUt5vKOHCBAyAKxR86J4jLKkpER5/IcffgCAXr16CQQCWtm6aN68eQAQFxenPCiTyYYPH/7sI60IGQDWqHlRu6mHSCTq27cvAOTm5tIK1nV37951cHAAgGPHjimPd7TBCkL6hjVqRjraYi4tLQ0AIiIiVDZsNlpr1qwBgCFDhkilUuVxtdv9IaRvWKPmorW1tX///s9ueFxdXc2c3P3666+0sr2o1tbWfv36AUBOTo7yuGLz6ePHj9PKhswQ1qi5WLt2rdqP33j99dcB4M0336QVTDM//fQTAHh4eDQ1NSmPf/bZZ2o/CgUh/cEaNQt37951dHR89sPgTp48ydxPvH37Nq1sGhs3bhwApKSkKA+KRCK1H8yHkP5gjZqF+Ph4AHj11VeVB2Uy2YgRIwBg9erVtIJp4+LFixYWFlZWVqWlpcrju3btAgBPT8+HDx/SyobMCtao6WPWWtrY2JSXlyuPf/PNNwDg4+Pz6NEjWtm0xOFwAGDq1KnKg3K5PCYmBgBSU1NpBUNmBWvUxMnl8lGjRgHAp59+qjwuFAr79OkDAP/85z9pZdNeXV2di4sLABw8eFB5vKioiHki68aNG7SyIfOBNWritm/fDgC9e/dWucL9+OOPASAqKqq7P4eenp4OAMHBwSpbUr3zzjsAMGvWLFrBuuj0o9MjS0fuF+ynHQRpDmvUlCl2RfrHP/6hPF5RUcHsinTu3Dla2XSlvb09KCgIADIyMpTHa2trnZ2djWS3Kolcsq1+W0xZjPMlZ4tiC+8r3gl3Eu623yWEHHp4CIrg2wff0s6INMcGZLrWrVtXU1MTGRnJPECpkJycLBaLFyxYwEwxdWtWVlYbN24EgNWrV9+/f18x7unpuWLFCgBISUmRSqXU8gGIiXhG5QxOFae6vXq+2/y03mnjnMZtb9gecT2iQlxBMRjSGdo9jvTl5s2bzI7xZ86cUR4vKCgAACcnp5qaGlrZdG7KlCkA8N577ykPisXigQMHAsDXX39NKxghZMXvK6AIXq18tU3+x7aEl0WXk6uTxXIxno2aAKxRkxUbGwsAb7/9tvKgRCIZNGgQAHz11Ve0gunD9evXrays2Gz2hQsXlMfz8/MBwNXVVeX5V4NpkbU4X3J2u+z2UKp+9RXWqAnAi3rTlJmZWVBQYGVltX79euXx7OzskpISf3//xMREWtn0ITg4eOnSpXK5nNl1RTE+e/bsSZMmNTU1ff7551SCnRedF8qEU5ynOFs4UwmADMCSdgCkFzweDwDGjh3LTDExCCHff/89AGRkZNjY2FALpx+rVq3atWtXcXFxWVlZcHCwYjw9PT0iIiI3N3fnzp0sFutFDzt0ydBLcZc0i3Rr0C3m7ucAmwGaHQF1C1ijpok5I/Pw8FAeZLFYJ06cyM/PnzlzJqVceuTq6vrjjz8GBQX5+Pgoj/fo0YPZYlUgEGhw2EfiR02yJs0iyUHeImsBAA9Lj+e+GXVfWKOm6a233srMzMzPzxcIBMwCdYaNjc2bb75JMZheTZgw4dnB5cuXi8XimTNnbt26VYNjsqxYxI48/33quFi42LPtAaBOUqfZEVC3gDVqmjIyMvh8fklJydq1a5n1QObp1KlTeXl5dnZ2GzZsYD6BysD8bfwB4Hb7bcN/a2QwOMVksnbu3GlhYbF58+aysjLaWeiQy+XLli0jhHC5XGZ7f8Mb6TDSycLpsPCwSC6iEgAZANaoyRo6dOiiRYskEsmHH35IOwsdubm5Fy5c8PHxYZ58pcKB7cBx59RL65dUL5ERmfJLApkm92qREWIprw5BJqa+vj4wMFAgEBw8eJBZoG4+mpubAwMDa2trf/zxR2ZralpEctGkikknH50Msw17w+2N3pa9ayQ1Bc0FVe1VNwfdPCI8MqViyrd9v13QcwHFkEgbeDZqyjw8PJjd4BMTE9vb22nHMag1a9bU1tZGRUXNnTuXbhJ7tn1BQAHPmycH+Wc1ny2uWpxely6Siz7w+EBGZNYsazdLN1u2Ld2QSBt4NmriJBLJ4MGDy8rKMjIykpOTaccxkMrKSubjUs6cOWNU+wZIiKRN3uZk4UQ7CNIlPBs1cco7d9TW1tKOYyBJSUlisXjhwoVG1aEAYMWywg41PXg2ahamTp166NChhIQEzdZOdi8FBQUTJ050cnIqKytjtqZGSK+wRs1CaWlpeHi4TCY7d+5cZGQk7Th6JJVKIyIiSkpKNmzYQHGCHpkVvKg3Cx3t3GF6FHuvmO0yL2R4eDZqLgQCQWBgYH19/e7du1977TXacfSisbExMDCwoaFh3759M2bMoB0HmQs8GzUXLi4ua9asAYCUlBSRyDSfqFm5cmVDQ8OECROwQ5Eh4dmoGZHJZJGRkZcvX167du3KlStpx9Gxa9euDR06FAAuXrzIbE2NkGHg2agZsbCw+Prrr1ks1vr166uqqmjH0bHk5GSpVLp06VLsUGRgWKPmJSYmZs6cOSKRyMTORvPz8/l8vpubG/PUFkKGhBf1ZqeqqiokJKS1tfX48eMxMTG04+hAe3v7oEGDysvLs7OzlyxZQjsOMjt4Nmp2/Pz8kpOTCSFJSUlyuZx2HB1IT08vLy8PDQ3lcDi0syBzhGej5kgkEoWEhFRVVW3fvn3+/Pm042jl/v37gYGBQqHw8OHDkyZNoh0HmSM8GzVH9vb2X3zxBQBwuVyhUEg7jlaWL18uFArnzJmDHYpowbNRM0UIGTNmzMmTJz/55JO//OUvtONoqKioaOTIkZaWliUlJQMHDqQdB5kpPBs1UywWa9OmTWw2m7mxSDuOJhS3dz/66CPsUEQR1qj5ioyMfOutt9rb21NTU2ln0cSuXbsKCws9PT27aX5kMvCi3qx13/kZU5olQ90dno2aNU9PTy6XC08eAaId5wWsW7euqqoqMjIyPj6edhZk7vBs1Ny1t7eHhYVVVFRs2bLl/fffpx2nS6qrq4ODg03pCQLUreHZqLmztrbm8XjwZHsk2nG6hNmkat68edihyBjg2SgCAJg0aRKfz1+2bFlWVhbtLM9RWFg4duxYOzu769ev+/n50Y6DEJ6NIgAAyMzMtLS0ZLaOp52lM4oN/FesWIEdiowE1igCAAgLC1u8eLFUKk1KSqKdpTPffPNNUVGRr69vSkoK7SwIPYYX9egxxSdw7N+/f/r06Sqv7t+//8CBAwYL4+XltWrVKpVBoVAYFBRUW1trwp+DgrolgtATmZmZAODv79/W1qbykoH38QwNDX02XnJyMgBER0fL5XKD/H4g1CV4Nor+0MmnE1+4cOHChQsGS+Lm5jZ37lzlkYqKirCwMKlUavKfEY26HaxR9JSCgoKJEyc6OTmVlZX16dOHdpw/TJ069dChQwkJCVu3bqWdBaGn4BQTekpsbOz06dObm5ufvTVJ0cGDBw8dOuTs7Lx69WraWRBShWejSFVlZWVYWJhEIjl8+HBsbCztOCASiYYMGVJRUZGRkcHcHkXIqGCNIjXefffd3NxcZ2fnpqYmNpvyJcusWbP27t3bu3fv27dv29jY0A2D0LOwRpEa1dXV/fr1k8vl1tbWdnZ2ffv2pRKjoaGhsbGxtbUVANLS0vCKHhknrFGkXkpKyqZNm4zkM++GDh168eJF2ikQUg9rFHXoypUr165ds7Gx8ff3pxKgvr6+vr7excVlypQpVAIg1BVYowghpBVc8IQQQlrBGkUIIa1gjSKEkFawRhFCSCv/D9+n9y9GSQgPAAAA/HpUWHRyZGtpdFBLTCByZGtpdCAyMDIyLjA5LjQAAHice79v7T0GIOBlgABGIOYHYj4gbmBkY0gAiTFDaCYmYmkOiD5GNocMIM3MTCRDA8RgEWRQAGnmBrqFkYmBiRkoxMDCysDKxsDGzsDBmcHEyZXAxZ3BxM2TwMObwcTLx8DLkSDCyMbBy8PNxSm+D+oLMOC3PTrxwCE+RQcQp1DC/4DYDj8wO3zFr/1/smTA7J4O5f0BIl4wtr3ZcXGYGvvuAleYXodCZWEwG2imw9EzDmC256Qldg/dlu0HsfNuT9uPYH8D0g2qSOL2SOrB7EkqnUC6YSlEXMUB5mgxACoYQ5MD8BlwAAABaHpUWHRNT0wgcmRraXQgMjAyMi4wOS40AAB4nH2T3W6DMAyF73kKv8Ci+I8kl22ppmkqSFu3d5i0y72/ZsNKUpWRYBSiLwYfHzrw8Ta8fv3AOmjoOoC4c5VS4JNjjN0FfAHH8/PLCKfr4XjbOU0f4/UdUAHFzti8Zw/X6XLbQTjBkwRVZUpAQbFX6iGGOI96lBzkQKy5T2ALKiVvgewghqKZerKMIhJt+xEUB2PoOedIlhE1SdkC1cCVo8AJlWWD641bX8whZilli0vGrZVQsEIi5w0uwwSrNJYPScoWV+bvY0Kh4qX/mxDjog2mrAn3UfyTkWa0Ao8k3SXdy8nth+6BMisZ01LR3sutN99umwXdIc/jcGe8xYrHaRyqFX1SNZw9AFdboYVU86CFVougRV+dgBapNhwtcu2rWM+w6Z94axCbPol3AKlph7jSyI3s4oqiNPKiS4fa6Ch+K60Ibcn+fPtzbd39AiTRwmxMknFoAAAAqnpUWHRTTUlMRVMgcmRraXQgMjAyMi4wOS40AAB4nFWOyw0DMQhEW8kxkTDiuzby0QVsEXtNCVt8cCRbCRfQ8GZgjG+dON4XX7P4cT+LoburgqDz4Qq9KIp6q5BdIiIVxvAmRzJmRpMhPLQ1koTYqyW0FUGtPIOWS5GazZwVLJi5ZNDX7SRYHHr+UNkUaJppjYQqbBLAy5cPcW1e/5SQqdDPtuz1jtjS6/4A3h84HcKM1YwAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 8
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "mol = combinemol\n",
        "for atom in mol.GetAtoms():\n",
        "    atom.SetAtomMapNum(atom.GetIdx())\n",
        "mol"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "PANVt3Ht8X96",
        "outputId": "923798c5-5077-43ed-8b30-a1bb4772f85d"
      },
      "execution_count": 9,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c82a2580>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3de1yMWR8A8N/MNE0Xle5I91BYUpZCii5eKyHyukRYl7W7lbyyCFlsYUlh3S/LrhDrvlpbWkmx6YouoiQRukz3yzQz5/3j2R0jk0pTT+n3/fTxac5z5szv0fTrmXPOcw6DEAIIIYQ+FpPuABBCqHPDNIoQQq2CaRQhhFoF0yhCCLUKplGEEGoVTKNtSCAQ5ObmZmRkVFVVfbjahysghDoyTKNtorKy0tvbW0NDw9DQsH///urq6h4eHiUlJe9X27VrV9++fZWUlO7fv09LqAihVpKhO4BPUH19/bhx4+7evevt7T116lQlJaWYmJgNGzYoKyvv3r1bVO3GjRvTpk3jcDhGRkY5OTkCgYDGmJuppqYmJyeHz+cbGBioqKg0Vu358+dcLldPT6979+7tGR5C9CBI2vbs2QMAO3bsEC98+vRpXV2deMmbN29iY2MFAsHp06cBICkpqX3DbJnCwkJ3d3c5OTnqbcNms93d3UtKShpUS0xMtLCwoOowGIwpU6YUFRXREjBC7QY/1Evf6dOn1dTUPD09xQsNDAxkZWUBoLa2lsfjAYCmpuaIESOYzE7wI6iqqrK1tT179uzq1atTUlIyMjJ27doVHh4eFBQkXi0vL8/Ozs7ExCQ5OTkvL+/QoUO///77V199RVfYzSEQCB4+fHj79u2srKzG6uS869mzZ+0ZIeoE6M7jnyBlZeWRI0c2dlRDQ+OLL74QL+n4V6MbN24EgBMnTogXPnv2TCAQiJeUlZWFhITU19eLSiZPnqyiotJOUbaQQCAIDAxUU1MT/S6YmpreunWrQTUul9vgV8bMzIyWgFGHhX2jUiYQCMrLy3V0dBqrMH/+fAMDg3aMSArCwsL09fXd3d3FC/X09KhvuFwum83u1q2bsrKyl5eXqAIh5NmzZ9ra2u0aa7MtX748JCRk+vTp//vf/3r37p2dne3v77906dLk5GQ2my2q9urVKwAIDg4eNWoUVdKtWzd6Iv5olZVw/z5UV0PPnmBmBhI/AFVXw6tXEsq1tKCuDt68gb59gcV651BGBqirg5ZWm8TcudCdxz9BHA7Hycmp+fU7+NVofX09k8mcMGFCYxW6devm6uoqXpKTkxMdHT137lwOh3PlypW2j7HF0tLSGAyGq6urUCgUFdbW1hYXFzeoeePGDQBISEho3wClpKKCLFxIZGUJwD9fOjrkl18k1Lx8+W0d8a/Tp8mWLQSAFBS8U5/HIwDE17d9zqODw6tR6TMwMMjJyaE7CqmpqKgQCoW9evVqrIKnp2e/fv3ES4YMGVJWVgYAXl5ejo6ObR5iy509e5YQ4uvry2AwRIUcDofD4QBAaWlpZWVl7969AeDly5cA0LNnT7pC/XgCATg7w927EBgI//0vqKtDRgasXw9z5kBtLSxc+E5le3vIzn6n5OBB2LEDhg+H3Nx2DLpT6gTjG52Ok5PTkydPEhIS6A5EOhQUFADg/S5CkYCAAA8PD/GS0tLSgoKC0NDQY8eOzZgxo81DbLmHDx8CwMCBAyUe9fX11dXVra+vh38/1AcGBtrY2NjY2Kxfv768vLw9Q/14v/4K0dGwdy8sXw46OiAnB0OGwMWLMGoUrFgBFRXvVFZQACOjt1/6+nDuHLi6QmfrgKIFplHp8/Hx6dat2+zZszMzM0WF9+7dKygoAID9+/dfvHiRvuhajMPh6OjoZDe4VGlKjx49Zs6c+e233168ePHp06dtFNtHKy8vV1JSaqyXc8KECWvWrGGxWABgZmZmZ2enrKw8bdq0/v37b9u2beTIkZ3jrrPTp0FNDebOfaeQxYJly6CsDK5dAwDIzYX8fAnP/e03yM4GH5/2iPMTQHevwqfpzz//1NDQYLFYQ4cOtbe3NzQ0BIC9e/eSd0fqx4wZY2RkpKWlBQA6OjpGRkaHDx+mNXDJFixYwGAwHj582GTNBmP31IyoDtixOHnyZAaDUVtb29Innjx5UvSj7Oh69SLDh0soz8ggAGTNGkIIMTQko0ZJqGNl9bac6hsNCyMREW+//vgD+0ZFsG9UOvLz8zdt2mRhYbFkyRIAcHR0zM7OvnDhwv379wUCwYQJE2xsbIYOHQoAN2/epD4mA8CMGTMafFgeNGgQAAiFwg41n3TlypWnTp1yc3M7f/68qakpABBCbty4MXjwYE1Nze3btxsaGk6dOjUlJcXNze3YsWPUoHZJScnPP/+sqqpqZmZG9xk0ZGJiQgjJyclpaWwuLi4AkJKS0jZxSVV5OfToIaGc6uctLQUAWLEClJQaVoiJgbt34cKFdwqnT2+LGD8RdOfxTq+4uHjFihXU7T1Uh1orGywoKDAyMjpw4ID4IDLtwsPDNTQ0mEzmwIEDbWxsNDU1GQzGqVOniNhIfWlpqYODA5PJHDJkyOjRo1VUVBQUFH777Te6Y5eAGn9fvXp1S59IDR6uWrWqLaKSMm1tInH+8uPHBICsXNnoEydNIkZGhM//5yF1Nfr0KamufvtVXo5XoyKYRj9eXV1dcHCwqqoqADAYDDc3tydPnrS+2Q0bNlB/4caOHZuZmdn6Bj9OdXX11q1bxaejl5WVnTx5cu3atX5+fvv373/06BFV/uTJk5cvX1LfC4XCP//8c926dcuWLQsODs7Ly6Mh9GYQCoVjxoxhs9kHDx4UdUSkpaVFR0cTQi5evOjn5ycQCPh8/owZMwICAmpqagghZWVlLi4uLBYrOTmZzuibydaW9OxJ3v9jTH0eP3ZM8rOysgiTSX766W0JTnhqCqbRjyEQCMLCwqgeTwBwcHBITEyUYvthYWFUhymbzfby8qqsrJRi400SP7vPP/9cKhfFUVFR8+bNe/PmTeubkpaioiJ7e3sAUFdXt7S01NXVBQA3NzdCyMKFCwGAx+MJBILvvvuOurlgwIABioqK8vLyR48epTv2xvH55OhRcvEiIYTs2EEAyPXrDevMmEFkZcmrV5JbWLqUqKkR8bccptGmYBptsYiIiMGDB1MJdODAgW00vZzL5Xp5eVE9pEZGRr///ntbvMr7xM9uwIABUjk7oVBobm4OAGpqaocOHaKxs+LKlSubN28WL7l79+6OHTvWrl27b9++27dvU4VcLvf58+eiOtT6AJs2bTp8+HBBg2zSoUREkMGDCQAxMCC1taSighgZkV69SGzsPxV4PLJpE2EwiKhT4ocfyJ49b1soLiaKisTP751mMY02BdNoC/z99992dnZUitHT0ztw4ECDgWmpS0hIoAamAMDZ2fnZs2dt91rvnx1f1DvWatnZ2ePHj6cat7S0vHfvnrRabqY7d+6MHj0aAFgsVnp6eusbLC0tHTVqVHh4eOubkoK//yZ2dv/cd6SnRw4c+KdnMzubmJsTAGJkRCwtiYoKYbGIj8/bfs8GI/UbNxIOh/zbRfMPTKNNwTTaLJmZmW5ubtQdL+rq6lu2bKE6y9pBfX19cHCwsrIyACgqKm7ZsqX1o1gNtNvZXb58mfrsLCMj4+XlVV5e3hav0kAbnZ2oC3vGjBkvG+Sd9vToEXFzIwwGASDq6mTLFtLg7AQCEhVFtm0ja9eSfftITs47R7lcUlb29mFkJLlxo+FLpKeTkycbNisUkpMnSafoI257mEab8ObNGy8vLxkZGQBQUFD47rvvSktL2z+M58+fu7q6Ur+3Q4YMiY+Pl0qz7X92lZWV3333HTWzvVevXsePH2+712rTs6P+vCkpKQGAiopKcHCwFC/em6OmoIAsWUJkZAjAP5/ExRMiakeYRhtVUVGxZcsW6vdERkZm8eLFdF50EEIIuXr1KjXyc/Lkyda3dubMGUVFRerslixZ0p5nl5ycbGVlJeqsePr0qXTbb7efXX5+/pw5c6gTMTc3v3v3blu8SgPU2elqatb36kWYTDJnDnnxoh1eFzUG06gEPB7vwIEDohXeHBwc7t+/T3dQ/6iqqjp48KBUmsrOzpaVlXVwcHjw4IFUGmwRgUDw008/UbuMKCgoBAYGSqWjmZaf3eXLl6nFD5lM5pw5c9puwf+6urqQkBBNTU1qjt0FLy9C35Q4JIJp9B1CoTAsLMzY2Jj6JbSysnp/Hd8OJSUlZeHChebm5kZGRqNHj961axePx2tQRygUHjlyZOzYsX369LGzs9u8ebNoO5M2HbNqjlevXs2ZM4fBYIwfP14qDW7ZsoX62Y0ePfrOnTtSabM5qqqq/P39qQ0OevTocfz4celOSKDemSYmJp3lndmlYBp9KyIiQrSPkJmZWVhYGN0RNeHXX3+VkZHp2bOnj49PYGDg7Nmz2Wz2+wuDzp8/n8VizZs3b8eOHbNmzWIwGNOnT6cl4MZERkZK63N9WVmZnZ0dXYucPnjwQLS6s6OjY25urlSa7XTvzK4G0yghhDx8+NDZ2Zl6m/bu3Vu6c33aSG5urry8vLm5OZfLFRVGR0ffvHlTvFpZWdno0aN/EVupd968eQCQn5/ffrE2m1AoPHv2rLOzs6GhoY6OzpgxY5rsBf7+++8tLS0tLS3F/x9oJBQKjx8/rqmpqaSkJJX/5DVr1ojemUeOHOn478wuCNMoIYQMHz4cAFRVVbdu3VpdXU13OM1CTbj5iE92J06cAIDIyMi2iKqVqIVdhg4dunnz5p07d1LrgPj7+zdWPz4+nsViffbZZwBQWFjYjpE2oaioKCIiQipNpaSkUPO0Oss7swvCNEoIITdv3lyxYsX7G0h0ZOPGjZOVlW1sWGb16tVTpkyReGjnzp0A0P4T4Jt07tw5AFiwYIF4r+LZs2cbu6arr68fMmTIqFGjtm3b1tHSqLiampqgoCArKystLS1tbW17e/urV69KrJmenu7u7q6vr6+hoWFhYbFv3z6qvKqqqh3jRS2GC+UBANja2tra2tIdRcsUFBTo6Og0tp7ey5cvJW5kQggJDQ01MjIaMmRIGwfYYvv371dUVAwODhbf2GPatGmN1d+6dWtaWlpycvI1agXiDqm2ttbJyen27dsuLi6zZ88WCoVXrlxxdnY+efLkrFmzxGvGxsaOGzeuV69ec+fONTIyevLkCTVgBf9uQIA6LrrzOPpI5ubmBgYGLX0WdSl64cKFtgiplZSUlGxsbBo76uLi4unpKXr46NEjOTk5Pz8/QsiPP/4IHfVqdN26dQCwe/du8cLQ0NAGn9D5fL6xsbGBgQEtN3egVupAawOjFunRo0d+fr5AIGj+U8LCwnx9fX18fCZPntx2gX2c6urqiooK0abN7ystLa34d/sgoVC4cOHCHj16rF69ur0C/EiHDx8eOHDgN998I144c+ZMeXl58ZKbN29mZ2f7+vqqqKi0b4BICjCNdlYjRozg8/kxMTHNqUwI2bx586xZszw9PamNPToaqnfiA38VoqOjjx07Rn1/4MCBmJiYPXv2UHdhdVh5eXkFBQWjRo0S76YQ4XK55ubmwcHBAHDnzh0AGDVqVFxc3O7du3fv3t3Sza8QjTCNdlbz5s1TUFDw8fEppXaD+FdNTQ282zdaXl4+bdq0zZs379u3r2PmUACQk5NTVVV99uxZkzXz8/NXrVo1YsSIHj16JCYmJiYmvnjxAgBSU1Nfv37d9pG2QGFhIQBQq7F8GLXd4eLFi21tbX/66Sc/Pz8zM7Pjx4+3eYhIGnCIqbPS1dU9dOiQh4fHgAEDZs+eraen9/z58ytXrjg5OQUHBy9atCg+Pr6wsDA/P9/Ozi47O3vMmDHZ2dmrVq2inu7g4ODg4EDvKTQwcuTIyMjI4uJidXX1D1S7du1aeXl5XFycaAlBioODw7Zt23x9fds4zBagLrGpjZrfp6qqKtrTiboMt7GxiYyM7NatW1lZ2RdffLF06VJnZ+cP/2+gjgCvRjuxWbNmJSQk/Oc//wkPDw8ODo6Lixs3btzixYsBwMPDY+3atQBQXFxsaGjo4ODAYrESxbx8+ZLu8Bvy9PSsra39+uuv+Xy+qLC2traurg4AMjIyqOtrDw+Pkndt3LgRAB4/fuzp6UlX8BL16tULAJpziU3VXL58ObXns4qKiqenZ01NTVxcXFsHiVoPr0Y7t8GDBx85cuT98un/7uM4ePDgiIiI9g3qIzk5Oa1ZsyYgICA1NXXy5MlKSkrZ2dmXLl3atWvX7NmzHR0dTU1NIyMjORwOh8MRfyI1XNO9e3dqY8GOQ1tbu2/fvlFRUXw+n1qvrzHUB//MzMye1Lad/17JdqgNYlFjWKLVZxGinb29vY2NTWFhYVxcXFpaGpvNnjt37rRp0xQUFJSVlUePHi1xP+QXL15wudyZM2d2tDQKACwWi5re5OjoKBpoKioqUlBQ4PP58fHxQqFQRUVFW1t79+7dpaWl1ArTfD5/+fLlL1++DAoKwkmjHR+DEEJ3DAh9soRC4Zw5c0JDQwcNGuTk5CQrK/vgwYPr16/fvn1bX19fW1vb19eXugtr48aN/v7+VlZWVlZWMTExSUlJO3fu9Pb2pvsMUNMwjSLU5i5fvnzmzJns7Gw2m21gYDBhwgRXV1ehUPjzzz+bm5uLVrA+e/ZsaGjoixcvDAwMFi5c6OTkRG/YqJkwjSKEUKtgBzZCCLUKplGEEGoVTKMIIdQqmEYRQqhVMI0ihFCrYBpFCKFWwTSKEEKtgmkUIYRaBdMoQgi1CqZRhBBqFUyjCCHUKphGEUKoVTCNIoRQq2AaRQihVsE0+qnh8XgHDx784osvBg4caG5uvnjx4szMTIk179y5M2PGDA8PjwblOTk53t7e48ePp7arbDfp6enffPPN8OHDTU1N7e3t9+/fL74pk0hlZeWuXbscHByuXLkiXn7v3r05c+aYm5uPGDHC09OT2i4UofZA0CekoqJi5MiRAGBvb79+/foVK1YYGxvLycnFxsaKVwsPD7e2tgYAWVlZNTU1UXlaWtq0adNYLJasrCwApKWltVvkZ86cYbPZWlpaS5cu/f7776dPn85isSZOnChep6KiYuXKld27d6fCCwgIEB06f/48m822srLauHGjj4+PoqKinp4el8ttt/hRV4Zb2n1S/Pz8YmNjDx48uGjRIqpk8+bNISEh5ubm4tUKCwsnTZp09OjR4ODgs2fPisqrqqqMjIyio6OfPHkyb968dgs7Pz9/wYIFpqamN2/eVFNTowpv3LhB3l1TnMPhlJWVHT582NraWkdHR/xQbGzstGnTfv31V2oPOGtr6+nTp4eFhVH7pCLUtujO40hqqqurFRUVbW1tm/+UJUuWiF+Nipw4cQLa8Wp006ZNABAVFdXM+mVlZfDu1WgDJSUlALBs2TIpBYjQh+DV6KcjNTW1qqpq7NixEo8+fvzYy8vr66+/njhxYktbFgqFNTU1rQ5QAhaLJScnFxcXx2azR48eLbHOpk2bUlNTz5071/xmuVwuACgrK0snSoQ+CNPop4MaEdLX15d4tKqqKikpqbCw8CNajouLs7GxaVVwjRg0aFBqampBQUHv3r1ZLJbEOjk5OampqS1q9tSpUwwGY9KkSdKIEaEmYBr9dIi2QZfI3Nz89evXH9cyi8Vqo93S5eXloanIjx071qI209PTAwMD586da2Fh0argEGoeTKOfjh49egDAs2fPpN6ytbV1VVWV1JsV0dbWfvjwoUAgaOyCtPmePXvm7Ozcp0+fnTt3SiU2hJqE80Y/HYMHD1ZQUIiOjqY7kBaztraur6+Pi4trZTu3bt2ysrJSUVGJjIxUVVWVSmwINQnT6KdDXl5+/vz5UVFRp06dEi+vq6sDAB6Pl5ubW1FRQVN0HzJv3jx5eXkfH58G4VGRFxYW5uXlNdnIwYMHHRwcbG1tb9++ra6u3laxIvQeBnl3ah7q1MrLy8eOHZuUlDRlyhRbW1uBQJCSknLlypWsrKynT58OGzZs9+7d3377bWlp6ZYtWwDg+vXrmZmZ3t7eAODl5dWrV68tW7aUlpampaVdvXr1yy+/1NDQcHFxGTFiRFtHfuLEiQULFujo6Li7u+vp6eXn51+9etXOzm7nzp1TpkyJioqiJjldvnw5Li6urq4uODjY1tbWysrKxMRk4cKFs2bNOnXqVK9evdzd3UU9rTo6Op6enm0dOULYN/pJUVZWjomJ2bt37/nz57du3aqsrKyvr79t2zZFRUUdHR1/f/9hw4YBQG1tbWRkJACwWKwBAwZQ38+fPx8Abt++/erVKwCwtLRMSUkBgCFDhrRD5HPnzu3fv39ISMhvv/1WXl6up6dnZWVF3QIwY8aM4cOHU9UePXpERWtpaVlZWRkZGVleXg4A+fn5lpaWAHDjxg1Rm2ZmZu0QOUJ4NYoAAMrKyiorKxvcGoQQag7sG0UAABs3buzXr9/JkyfpDgShzgfTKIInT57s2bOnpqbG1NSU7lgaIoSIf05HqAPCNIrAy8uLx+MtWrSI6l6k5Obm0hfRW05OTg4ODlR/KEIdE6bRru7atWvh4eHKysobNmwQFd68edPExGT58uX0xfUPe3t7APD29q6vr6c7FoQkwzTapdXX11O5csOGDdRNUAAgEAiWLVsmEAhomX1ZW1t7//590cPly5f36dMnPT394MGD7R8MQs1C8wpTiFY7duwAABMTk9raWlHhvn37AEBPT6+qqqqd48nLyzM0NOzZs2d5ebmo8MKFCwCgqqpaVFTUzvEg1Bx4Ndp1FRYWUgt9hoSEcDgcqrC0tHT9+vUAEBQU1EbLkXxA7969dXR0CgoKAgICRIWTJ092cnLicrnff/99O8eDULPQnccRbaiV4R0dHcULly1bBgCjRo0SCoW0RJWYmMhkMmVlZbOyskSFaWlpMjIyLBbr/v37tESF0AdgGu2ikpOTWSyWjIzMw4cPRYUZGRlsNpvJZCYkJNAYG3Xz0uTJk8ULv/nmGwAYO3YsXVEh1BhMo12Ura0tACxfvly8cPz48QDw1Vdf0RUV5dWrV9TC9devXxcVlpSUUENely5dojE2hN6HabQrCgsLAwA1NbXi4mJR4dWrVwGge/fub968oTE2SmBgIAD079+/vr5eVLhr1y4AMDY2Fh8QQ4h2mEa7nJqaGgMDAwDYv3+/qJDH4/Xr1w8Adu7cSWNsInV1dX369AGAPXv2iArr6+s/++wzANi6dSuNsSHUAKbRLmfjxo0AMHjwYD6fLyrcvn07AJiamvJ4PBpjE3f+/Hl4b54TdTuTkpLSy5cvaYwNIXGYRruW/Px8RUVFAPjrr79Eha9fv+7evTsAXLt2jb7QJHBycgIAT09P8UIXFxcAWLBgAV1RIdQAptGuZfbs2QDg5uYmXrho0SIAmDBhAl1RNUbiPKcnT55wOBwmkxkfH09jbAiJYBrtQu7cucNgMOTk5J4+fSoqpGY+sdnszMxM+kJr1Ndffw0A9vb24oW+vr4AYG1tTdfkVoTEYRrtKoRCIbX0/dq1a8XLJc586jiKi4vfn+dUXl7es2dPAAgNDaUxNoQomEa7Cmq3dx0dnYqKClHhmTNnAEBTU5PL5dIY24eFhIS8P8/p8OHDANC7d+/KykoaY0OIYBrtIioqKnr16gUAv/zyi6iwurqamvl04MABGmNrksR5TgKB4PPPPwcAf39/+kJDiBBMo13EqlWrAMDKykq8M5Ga+WRubi4+86ljkjjPKS4ujsFgyMvL5+bm0hgbQrjC06dPIBBEREQwGIzg4GDR5sMvXrzYunUrAAQHB7NYLFoDbJq9vf3EiRMrKirWrVsnKrS2tv7vf/9bU1ND/ZFAiC64M2iXwOfzb9686eDgICpxd3c/efKkm5sbdWNox5ednT1gwID6+vq7d+9SH+cBID8/39TUtKqqKjo6evTo0fRG2IAQhADAxJXRuwCW+NYRqIOIiIjYt2/fiRMnIiIiKisrTU1NmUwJv40FBQUhISFZWVkWFhbi5Tk5OUFBQUePHr1161a3bt309PSYTKaRkZGowp07d3x8fOTk5C5dukRNvO/4ZGRk4uPjnzx5cuHCBWqZkp49eyorK/N4vJiYGGNjY1Ea/euvv/bv36+lpSVazx8ASkpK9u7de/DgwaioKC6Xa2ZmJvG/9OOUCcr+rv47rSathtRos7UZwACAYZnDrpRdmaU2S+JTMmozCBBFpqLEo5m1mc94z6qF1WoyatIKErUhunsV0DvKy8upW3eMjY0nTZpkbW3NYrGGDx/eYCQ9NTV18eLF8vLyADBgwADxQ7du3VJQUNDS0po6deqAAQMYDMb7d6BTG4f4+fm1+flISUxMjIaGBovFkpWVBQBNTU0mk/nDDz8QQqqqqqh5+DweLywsbPjw4dQbe/fu3aKnJyUlaWlp6evrT5061c7OjsFgjB49Wiq3vVYIKhbkLpBJkoFEoL7M0sziq+IJIUMzhk540vCOBgERXC697JDlAIlgn2UvqUmSWp3KTmIzE5njH49vfYSoHWAa7Vg8PDwAICQkRDQWlJCQ4O7uXlpaKl7tm2++mTNnzrlz5xwdHcXTKJ/P79Onj76+fmFhIfXQ1dVVRkbm8ePHDV7o6tWrnWWq0Js3bzQ0NHR1dR8+fCia57R+/fp9+/Y1qGZpabl+/forV640SKPu7u6TJk0S5c2ffvoJAM6dO9fKwAREMDZrLDORue7luqd1TysEFXGVcVaZVpOeTCKNpFGHLAfVFNVZT2fpPtCVmEYFRGCdaT0sc9hn6Z9hGu0sMI12IM+fP2exWFOnTm3+U1xcXMTT6K1btwBgy5YtopJ79+4BwPfffy/NQNsXtaHI1atXCSECgcDCwkJOTu7GjRuN1X/+/HmDNCoUCqurq0UPExMTASAoKKiVgZ3jnoNEWP9yvXhhpaCSJ+SRRtJoTl0OdXRY5jCJaXT7q+0ySTLJ1cnm6eaYRjsL7P/uQKKjowUCgaurq8SjDx488PPze/z48QdayMjIAABHR0dRydChQ1VVVdPT06Ubanu6ceOGvLw8taQ0k8n8+eefMzIyxo4dSx09evTotm3bPsZgxXYAAA0cSURBVNwCNS+K+p4QcuLECRaLJT7g9nEulF5gAtNb01u8UJGpyGawxUvKBeVptWmVwkoAMJQ1bHBUXC4vd0PBhuVay83lzVsZG2pPMnQHgN569uwZABgaGko8mp6eHhAQMGrUKGohTony8/MBgLpRUkRbW5sq76Ty8vL09fVFI0LUVHyRkydPPnr0aOXKlU228+eff4aGhiYmJhYVFZ07d45qx8rKqrS0VGL9z0I+e6D7QOIhc3nz04an02vTtdnaTY4CXSu/NvPpzHCT8P8o/+fDNZfkLVGTUVvXc92Hq6GOBtNoB8Ln8wGgsf04XVxc8vPzNTQ0PtBCbW0tAFDjMCJycnLV1dXSC7O91dfXKykpNXY0LCysvr6+mU2x2eyePXs+fvz4jz/+mDBhApvNzsrK4nK5Eitr87Uf1T6SeKg7qzsAVAgqtGW0m3xRc3nzbTrbTOVMP1ztWPGxP8v/vGx8uRuzW5Ntog4F02gHQq3B8fz588GDB79/VF5eXkdH58MtUNehRUVFVFOUN2/eNJgR1bmoq6vn5eV94Ggz23FycqJmQVy5csXFxcXMzMzb2/vvv/8WCAQS6zM1mMJuQomH5JnyAKDIVCziFzX5uqZypk3m0CJ+0coXK91U3SaqTGyyQdTRYBrtQKhkl5SU5Ozs/HEt6OnpAUB6ejq1IwgAFBcXv3nzhirvpCwsLJKSkvLy8qR1Fs7OzkpKStHR0d7e3h/oIWlSX7m+50vPVwmrGpv+2Xy+L3yL+EV1wroleUuokuf1z4sFxUvylmzT2abCUmll+6hN4RBTBzJ8+PB+/frt2bOnoKDg41qws7PjcDjUYk6UU6dO8fl8anymk5o7dy4hhFrf7+NaSElJyc7OFj0sLi6urq6mNh9tjQnKEwREcKz4WNNVm6Inq+em6sZhcrgCLvUlIIJ6Us8VcKm7oVCHRvNMAfSu6OhoeXl5XV3dn3766c6dOxERERs3bnRxcSGEXLhwQU1NLTIykhBSXFwcFhYWFhY2dOhQXV1d6vuamhry7yokS5cujY2N3bt3b7du3UaOHNnZlzdeunQpADg7O1+4cCE+Pv7cuXNz5syhtuSbMmXKwIEDqWrx8fFhYWH79u0DgAULFoSFhcXGxhJCxowZo66ufvjw4aSkpD/++GPEiBEsFuvWrVutjKpWWDsofZB8svyJ4hP1wnpCiIAIwsvC71TeIWITnhKrEr/J+ya9Jl38uY1NeBLBCU+dCKbRDichIWHcuHGiYSJ9ff01a9bU19fHxMQ4OzsnJCQQQu7evfv+X8SCggJCCJ/PX7t2raqqKgDIysrOmDGDmorfqQmFwt27d4s+gHM4nBEjRvz++++EED8/v3nz5lHVqJsXxLm6uhJCiouLFy5cSG1CBQCmpqbS2uz+Be+FzSMbSATlFGWThyZyyXLMRCY1k1SURk+VnIJECC8LJ4RcLL0out9J9CWxZUyjnQguTdJB8Xi8kpISVVVVDofzEU8nhJSUlCgpKYnSMSFEtLxT51VZWVlZWamtrf0R58Lj8d68eSMvL9/8UalmSqxOTKhOqBZW92b3tlSwNOIYAcBfFX9xmJwRiiP4hF8jrFFgKrAYrJf1L2MrYxs83U3V7f02I8oj5JhyNt1spBsqaguYRruE0NDQgICAqKgoLS0tumNB6FODQ0xdQmhoaFpa2tq1a+kORJpKS0vr6urojgIhTKNdQ3BwMIfDOXLkSEJCAt2xSI23t/eAAQPu3LlDdyCoq8M02iWYmJh8++23QqHQ29v70+jGSUxM/PXXX58/f66pqUl3LKirwzTaVfj7+/fo0SMuLq6zLHf/AYQQb29voVC4YsUKExMTusNBXR0OMXUhhw4dWrx4ce/evR89etTYnfudwokTJzw8PLS1tbOyslo/ix6hVsKr0S7kyy+/HDp0aH5+/vbt2+mO5eNVVlauXr0aALZt24Y5FHUEeDXatcTGxtrY2MjJyWVkZOjr69Mdzsfw8/MLCAiwtLSMj4+X4n5KCH00fBd2LSNHjnRzc6upqVmzZg3dsXyMp0+fBgUFUZtFYw5FHQS+Ebuc7du3KygonDp1KiYmhu5YWszX17e2ttbd3X3UqFF0x4LQPzCNdjm6urorVqwQDXbTHU4L/PXXX7/99puCgsLmzZvpjgWhtzCNdkWrVq3S09NLTk4+fvw43bE0l0Ag8PHxAYA1a9Z06uVT0acHh5i6qJMnT7q7u2traz969EhFpROsCrxv376vv/5aT08vIyOjU8/WQp8evBrtombNmmVjY/P69evAwEC6Y2kal8tdv349AAQFBWEORR0NXo12XUlJSZ9//rmMjMyDBw/69u1LdzgfsmzZspCQkDFjxkRFRdEdC0IN4dVo12VhYeHh4cHj8ZqzOzGNMjMz9+7dy2Kxdu7cSXcsCEmAabRLCwwMVFFRuXTp0vXr1+mOpVHLly+vr69ftGiRxA1TEaIdptEuTVtbm7qxkkpVdIcjwdWrV8PDw7t3775x40a6Y0FIMkyjXZ2Pj0+fPn3S09MPHDhAdywN8Xi8FStWAIC/vz8uiIc6LBxiQnDx4sUpU6aoqqpmZWVpaGjQHc5b27dv9/X1NTU1vX//PpvNpjschCTDq1EEkydPdnBw4HK5R44coTuWt3g8XnBwMAAEBQVhDkUdGV6NIgCAtLS0e/fuzZ07t7H1PlavXv33339L/XWNjY0PHTrU2NEXL16cOXNm+fLlUn9dhKQI02gXkpube+vWrcLCQnV1dXt7e11dXYnVYmNjExMTORzOmDFjRPNJx48f/8cff0g9pEGDBqWmpooepqamJiQklJaW9u7d28nJSVVV9f2nCIXCqKioR48eKSoqDhs2rH///lKPCqEWwTTaJdTV1S1duvT48eMyMjLa2tqFhYV8Pt/Pz2/Dhg3i1aqqqqZOnXr9+nVNTc3S0lKhULhp0yZqKD81NbW4uFjqgSkpKX3++ecA8OrVq5kzZ968eVNeXl5DQ6OgoEBBQeHAgQMzZswQr19YWOjg4JCRkTFo0KDS0tKcnBwvLy/qsz9CtCGoC/Dw8ACAVatWlZeXE0Kqq6t/+OEHa2tr6qGIl5cXg8E4evQoIaSmpsbd3Z3BYMTGxrZ1eDwez8LCgtq7lMfjEUKKiorc3d3nzp3boOZXX33FZrOTkpKoh+vWrQOAqKioto4QoQ/ANPrpy8zMZDAYM2fObFAuEAgalKirq48ZM0b0sKioiMPhLFiwoK0jDA0NBYDt27d/ODxCyMiRI42NjUUPk5OTASA4OLitI0ToA3Ck/tMXHh5OCPnyyy8blFOjSZmZmefPn6+urhYIBCUlJeJr0Kmrq/fr1+/evXttHeG1a9cAYMGCBe+HBwB37969evUq9b2JiUlubu7ly5eph9Qm9VS3AEJ0kaE7ANTmHj9+DACNbUR85syZDRs25OTkGBoa9u3b9/bt23V1dRwOBwBKSkr4fH5RUREA/PjjjykpKVKPTV9fPyAg4PHjxxoaGhIHlAAgMDDw1q1bXC4XADZs2BAdHT158mRHR8exY8du3bo1KChoxIgRUg8MoebDNPrpq66uBgB1dXWJR2fOnDlkyBBtbW0A8Pf3nz17trW19aRJk168eHH58uXKykp5eXkAiIqKaqOR+oCAgOrq6sbCA4C1a9eKRrcUFRX79OnDYDAKCwtXrVqloqIiI4PvYUQzfAt++qhdiAsKCvr06fP+0b59+4pmNc2cObN79+579uy5ePFiv379zpw54+PjQ2XhlStXuru7Sz026gpUWVn5wYMHjdUR/8zu6ur68uXLxMREVVXVtLS0tWvXenl58Xi8//3vf1KPDaHmortzFrW5gwcPAsClS5da+sTa2lpFRUU3N7e2iErckiVLACA/P//D1TIzMwEgICBAVCIUCs3MzD777LM2DhChD8Ehpk/fxIkT2Wx2SEhISzewO378ODWTtI0CE3F1dQWAJpcT5fF4on8pDAaDzWYLBII2DQ+hJtCdx1F7oKbZu7q63r9/n8/nv379+vTp0zt37iSE/PLLL46Ojq9evSKEVFVVhYeHV1ZW5ufn//zzz8rKytbW1nw+vx0inDhxIoPB8PX1zcnJEQqFz54927Nnz2+//UYI8ff3d3V1JYTw+XxjY2N1dfWYmBhCiEAgCAoKAgB/f/92iBChxmAa7RKEQuGPP/5IdZJS2Gz2/PnzCSE7duzo2bNnXl4eIeTcuXOiERsmkzlt2rSSkpL2ibC6upqaWi+KUElJiZpJOm/ePDMzM6paWlrayJEjAUBFRUVOTo7JZC5ZsqS+vr59gkRIIrwZtAupq6tLTk4uLi7W0tIyMTGROMHo9evXmZmZhJD+/ftraWm1c4Tl5eUpKSkVFRW6urrGxsaKiooSq2VlZeXk5ACAhYVF+weJUAOYRhFCqFVwiAkhhFoF0yhCCLUKplGEEGoVTKMIIdQqmEYRQqhVMI0ihFCrYBpFCKFW+T8t4V3xEPio8wAAAR96VFh0cmRraXRQS0wgcmRraXQgMjAyMi4wOS40AAB4nHu/b+09BiDgZYAARiDmB2I+IG5gZGNIAIkxs3GAaCYmCUY4iwnOQsiywFmscBYbB5jFyCjBzuaRAWQxMzNKcCCYnAgmF4LJjWDysHlogJgsEryCHAogoyT4uIHuZGRiYGIGCjOwsDKwsjGwsTNwcGYwcXIlcHFnMHHzJPDwZjDx8jHwciSIAF3Ny8PNxSm+D+pDMOC3PTrxwCE+RQcQp1DC/4DYDj8wO3zFr/1/smTA7J4O5f0BIl4wtr3ZcXGYGvvuAleYXodCZWEwG2imw9EzDmC256Qldg/dlu0HsfNuT9uPYH8D0g2qSOL2SOrB7EkqnUC6YSlEXMUB5mgxALWcRbzeXoRXAAABgnpUWHRNT0wgcmRraXQgMjAyMi4wOS40AAB4nH2TW06FMBCG31nFbMCmc+nt8QAnxpgDiR7dg4mP7j9OwV5QYqGk03xMmf8fBsjjZX7++II6aB4GAPvPnVKCd7bWDjfICxivj08LTPfLWHam9W25vwI6QNF39Dqyl/t6KzsIEzyIcc4xBSDj0DvyYI3dRnuVMsiG2EUfQBeUUuxALCBnEE1ykTxpRhGxul1BKqBk0BrPMVrSjOiCpA7kAjoFK0eGAzqWxknhvHL1YDY2Skod5woXlKuVkNFCLMfG+cJFWKFKo/mQJHVcKFzavo8JhVIu/XfCWMW2uzYYogt4hqaK4o+MtKHNj6p3dQbpkPVvUqzeIPefekJWc1A2MW3Yizo5v9qjfTZ95tbZ2RO0OnRd5kP37f04rsvc+jFf1LpOg9YJrCE3w0VDab5qv+v0rR1QZ2iuo87YzBU1rqmYcraDVlt+7DXJp6naXemUd7irMAdZPHSdkpIfqRehLznH5ffV9fANeXzDmGBvkDcAAAEjelRYdFNNSUxFUyByZGtpdCAyMDIyLjA5LjQAAHicbZI7b8MgFIX/SsdWIoj7MnAzVV689LFHHtKoWyxHUbvlx/dSYqsVXjjwcTi8bn/oB1QYfwWrUBWuIlW68fA2aBz94TRoGqFIHksLoQpUsZCTApmjPyvw+HB73LEXESKHXqATcvsdeSRJ0ZliztkI+CwJO/Mwcyie4DtKKaCZQCKbaSXoKUIJWlaRD4lLzhKM3nIDu/2ytzkApWQQAmN2tuPdYj2ISeJ/lLGg8Gd6nV0zVgI+RGAqdrvkvfvkjl/z9H6dLwp+ms/PNno5Xl6/p4/PqwfFFtrrt5CUW8gqLRTtWthpbGHU1MKkuYXZfnjj9MF+fAODVcAGthrbuBiQVcgG5tsPXsa6Ow6gwpMAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 9
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "editable = Chem.EditableMol(combinemol)\n",
        "editable.AddBond(2,11,order=Chem.rdchem.BondType.SINGLE)\n",
        "editable.GetMol()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "VaYEBXSvOMW6",
        "outputId": "00cfb980-c0ff-45d8-901a-67c33c0b47a0"
      },
      "execution_count": 15,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c79fc820>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dd1iTV9sA8DuEEAhhb5AhIggIjqg40DpfRNQqBQcVnIVWKQ4cuIra+oktbQFHC67iFrfWYqkgAiqKVAEZimwQGbLCzni+P06bxjBEyYTzu3q9F5zn5OQOvtw85zmLRBAEYBiGYR9LTtIBYBiGyTacRjEMw3oFp1EMw7BewWkUwzCsV3AaxYSPy+VyuVxJR4FhYoLTKCZ8165ds7S03L9/f01NjaRjwTCRw2kUE75Lly7l5eUFBASYmZmtWbMmJydH0hFhmAiR8LxRTOi4XG5cXFxoaOitW7fQ/8EmTJiwdu3a+fPny8vLSzo6DBMynEYxEXr58uWhQ4eOHz/e2NgIAObm5t7e3qtWrdLS0pJ0aBgmNDiNYiLX0NBw/vz5kJCQ7OxsAFBUVHR3d9+0aZOdnZ2kQ8MwIcBpFBMT3NPH+iqcRjFxy83NPXbsWHh4eF1dHQAYGhp+8cUXvr6+2trakg4Nwz4GTqOYZDCZzHPnzoWGhmZlZQEAlUpdsGCBv7//sGHDJB0ahn0YnEYxSSIIIjY2NiIi4sqVKxwOBwAYDIafn5+Hhwfu6WOyAqdRTCrk5eUdOXIkIiKitrYWAAwMDLy9vXFPH5MJOI1iUqSxsfHs2bNhYWGZmZkAQKVS586du2HDhrFjx3b/wvb2dgUFBbHEiGGC8ComTIrQ6XRvb+/nz58nJia6u7uz2eyLFy+OGzdu1KhRJ0+eZLFYAvWzsrIWLVqkpaVFpVIVFRWdnJxSUlI6NvvmzZsdO3YYGRn5+vqK5XNg/Qu+G8WkV35+fkRExJEjR9DafH19/aVLl3799ddGRkYAkJycPGPGDBqN5ufnZ2NjU11dHRERkZmZmZOTY2Jiwmtk9erVx48f19bWrqurmz59+rVr1yT2ebA+CqdRTNo1NjaeOnXq4MGDvDH9qKioOXPm2NjYVFZWpqammpmZoZrt7e1JSUlTp07lf/n58+cHDRo0atQoW1tbS0tLnEYxocNpFJMZSUlJYWFhf/zxR15eXn5+/vjx43ft2hUYGNjDl9vY2OA0iokCfjaKyQxHR8eoqKiSkhI9Pb0nT54AwJgxYzqt+dtvvxkYGKSlpYk3QKyfwmkUkzEaGhoA8PbtWwAYMGBAp3W0tLTs7OyUlJTEGhnWX+EZzphMolAoANDc3Nzp1Tlz5syZM0e8EWH9F74bxWSSsbExABQVFUk6EAzDaRSTTZMmTZKTk7tx44akA8EwnEYx2WRmZubm5nbu3LnTp0/zCrlcbm5uLgAUFBRcv34d7SCFYaKGJzxhsqqmpsbFxSU5OZnBYDAYjIaGhkePHjU1Nb1+/frXX3/19fV9/Pjx6NGj4+Li9u3bBwAPHz6k0WjDhg1TUlLCt7GYEOEhJkxWaWpqJiUlXbx4MSYmpri4WEdHZ9myZR4eHmQyedasWVeuXBk8eDAA0Ol0c3NzAED/CwCKioqSjBvrc/DdKNZHEARBIpEkHQXWH+Fno5jMKy4uNjExsbGxkXQgWD+F70Yxmdfc3KysrEylUltaWvANKSZ++G4Uk3k0Gk1bW7utra2iokLSsWD9EU6jWF9gamoKAMXFxZIOBOuPcBrF+gKURvGiJkwicBrF+gKcRjEJwmkU6wtwGsUkCKdRrC9Ap4bgNIpJBE6jWF+A70YxCcJpFOsLUBotLCyUdCBYf4TTKNYXaGlpqaioNDQ01NfXSzoW7D3YBPtF64unzU+r2dW8wkkvJ3kXe3f1krfst2yC3dXVBk5DLae2idsk5EB7DKdRrI9AGznjqaPSrI1o21y2WTtde0jWkJE5I/XS9Wa+mlnYXggAb9lv6zmd/AlMakxyzXfVTdd1y3frtM3M1kyddB3NNM1lhctEGXt3cBrF+gj8eFTKEUB8lv/ZDxU/eGh6JFslv7R9eW7guVdtrzaUbujqJc6vnKfmTmVymJrymgR0smydC1yfYh97JfuBCgNFGft74I3ysD4Cp1Epd63u2q36W1v0tgQZBaGSwdTBk1Umq8ipdPWSPYZ7zlHPqZPVh2YN7bTCgcoDyU3Jj60eLyxYKJKgewanUayPwGlUyl2ovSBPkt+st5m/UFdeV6BaO9HeyGlUJavKk+RH00Z302Bxe/HO8p1rddaOpI0UfrgfAnfqsT4CTx2VchktGQMoAzTlNbuvdqrmlFa61oOmB+9t0LfEV52svttwt5AC/Hj4bhTrI/DdqJSr49SZKZi9t5qdot163fUDKAO6r3a65vTN+pvXBl2jy9GFE18v4DSK9RE4jUo5JTmlWk7te6uNUR4zRnlM93Xest/6l/rPV5//qdqnQoquV3AaxfoIQ0NDBQWFioqK1tZW2TpticPhVFdXc7lcXV1dMpncVTUul8vlcuXlBX9nWSxWdXW1jo5Ox0tSZRB1UFJjEptgy5N6G+fGso2V7EplOeWAsgBU8pb9NqM1I6AsYLfhbiqJ2utgPwx+Nor1EXJyckZGRgRBlJSUSDqWnnr79u1XX32lq6urr69vaGior68fGBjIYrEEqrW0tISHh9va2tLpdP6VWtXV1UuWLKHT6YaGhjQabdasWWVlZWL9AD3ABW4ztxkAnFSdmrnNV+uu9r5NFsFi0BjZrdl3mHfQf83c5kpW5R3mHQ7B6X37H0qq/3xh2AcxNTUtKCgoLi5GZ4JKOSaTOWnSpLy8vLVr186cOZNMJkdHR+/fv19dXX39+vW8ardv3/b09ORwOCYmJm1tbbwky2Kxpk6dWltbe/ToURsbm8ePH/v7+7u6uiYnJ0vPSSrV7GqvQi95kvz1QddXaq38ueLnNSVr9Ch6k+iTUIXU5lS6HN1K0Yr3krSWtN/rf/fS9DJWMO6q2dNmpwVKBmcOHq40/KL5RVF8ivcjMKyvWLp0KQAcO3ZM0oH0yJYtWwDg/Pnz/IVpaWkcDoe/pLi4ODY2lsViHTt2DABevnyJyhsaGrZv356RkcGruWHDBgAoLCwUQ/A9kcBMMEo3glTQSdMpaCsgCCK9Od3iuQWkgsVzi0kvJplkmEAqfPP6G4IgbDJtFuQvIAjiaPVRSIV7zHsEQVyvu27+3Nz8ubnC3wq0pzTz5+a2mbadvpfFcwu3PDfxfbZ34btRrO+QrVGmM2fODB06dOHCd+aN29vbC1QzNjZG61wFqKiofPfdd/wl6GzUoqIi9HOQIAKIsMqwTWWbWARrjPKYCwMvoDF6OyW7TJvM2w23/27+u41oc6e4j1UeO4o2CgB2GOzQJGsCgLu6+0T6RGOKMQDYKNps0dvC33JXz1W/NfxWi6wl6s/VFZxGsb5DhqaOVldXl5aWTpkypdOrXC6XQqGsWLHiyJEjPW8zNzcXAMzMzIQS4UdDHfnohmgSkPx0/YKNgikkCu+qAklhrtrcuWpzBV61WGMx+kKVrKpKVkVfW1AtLKgWPXnTRRqLhBH7R8JpFOs7ZOhutLa2FgCMjIw6vUoikTw8PMaOHdvzBplMZmRkpJOTE/pbIimJjYmLCxaXscq05bUjTSNnqc2SYDBig9Mo1neYmpqSyeSOI91SiEqlAgCTyez0KolEOnXqVM9b43K5q1atqq+vDwkJEU58H66rjnx/gNMo1ndYWFi0trZK+fRJRF9fX1FRsaCgoPdNsdlsHx+fa9euXbt2bciQIb1v8CNUVVUFRAZETovkEtwA/YBvDb7t/eRQGdKPPirWx1RXV58+fTo1NbW5udnIyGjRokXjx4/vmEMfPHhw/vz5169fDxgwwNnZ2cnJSSLRClBQUJg0aVJ8fHxVVZWOjs5Ht1NZWblo0aKUlJRLly45OzsLMcKeS0xMXLx4cVlZ2awzs3xn+zqrSiYMCcLT7zGZFBcXZ2FhERAQUFlZSaFQYmNjHR0dt23bJlBt165dEyZMePz4MZ1Ov3PnzsyZM/ft2yeRgDvaunVrS0vL4sWL3759yytMTk5ubW0lCCIiIiIhIaH7Fp48eTJmzJiSkpKHDx/OmTNHxPF2giCI0NDQadOmlZWVjRkz5tD4Q/0whwLgeaOYDCovL9fQ0Bg4cGBubi4q4XK5P//8c2RkJH+1qqoqXV3doKAg9C2LxXJ0dFRUVGxqahJ3xF04fPiwgoKCsrLyjBkz5s6dO2jQIACIjo5ms9kAsGLFClTNycmJwWCgAbShQ4cyGIzLly8/efJEQUEBAAYPHszgc+LECfEEX1VVNWvWLAAgkUh+fn7t7e3ieV8phDv1mOw5evRobW3tuXPnLCz+mQ1DIpHWrVsnUE1bW7uwsJC3vl5eXt7NzS0pKenFixcjRowQa8Rd+Oqrr2bNmnXp0qWcnBx5efmxY8fOnDkTxXb//n09PT1UbeLEiQKDUUZGRnQ6nX+xEw+a8JSTkyPS56SPHz9euHBhYWGhtrZ2ZGQkyqf9l6TzOIZ9sGnTptFoNIHVPjwHDhzYs2dPp5d++OEHAHj27Jkoo/tgDx484N0yC0VWVhaJRGIwGOHh4S0tLUJsmSAILpcbEhJCoVAAYMyYMQUFBcJtXxbhNIrJHktLyyFDhnR1ddKkSWZmZl1dMjAwYLPZIgvtw3C53ODgYAqFQiKR4uLihNXs1atXtbT+WdKjr6+/a9eu+vp6YTX+6aefAgCJRAoICGCxWMJqVqbhNIrJHnNz8+HDh3d1lclkdpo1zp8/DwACz08lqK6uzs3NDf59ttjW1ibExltbWyMjI4cNGwYA6urqjY2Nwmo5LCxMW1v71q1bwmqwDyARRCfn7WGYNBs7duyrV6+qq6vfX/Vfd+/enT179rx5806fPi0NGyA9efJkwYIFBQUFampqR48eRflUFO7du1dSUrJkyRL+wtevXzOZTB0dHU3N7o70YDKZKiqC582x2ey0tDQDAwNDQ0PhhyujJJ3HMfEpKCiIj4+/f/9+Q0PDeyuXl5fn5eWVlpaKIbAPtXr1agDIy8vrYf2zZ8/SaDRXV1dp6ISiZ4tokH3UqFE9/xRCeeuwsDDealESiTR58uTMzMyO1W7fvv2///2PRCIFBgbyXwoNDeXNcjUwMBDbrAAph9Nov5CQkMC/dZCiouKWLVu6eURYVFSkoqJCIpHGjBkjzjh76NGjRyQSyd3dvatRJh4Wi4X2o/P29paGHCrSjvx7+fv7A8Ds2bP//PPPtLS0S5cu2dnZjRgxgv/H2NTUNHToUCqVOnPmTADYvn0771JwcLCSklJISEhWVtbjx4+dnZ3l5OQSEhLE+RGkE06jfd/Dhw+pVKqFhcWVK1fKy8vz8/O/++47JSWl+Pj4rl4yd+7cwYMHjx8/XjrTKEEQGzduBIApU6acPXs2Li7ut99+c3V1DQ0NJQhixowZVlZWBEFUVlY6OjoCwIwZM8L5pKWlsdns3bt3i3mUOSUlZeDAgQCgpqZ28eJF0b1RQUHB8ePHBQrT0tJIJNL8+fO5XC6vsK6u7s2bNwI1Y2JiGhoaOByOQBp98ODBpUuXeN+WlpYCwLZt20TwCWQMTqN939ixY1VVVYuLi/kLKysru6p/5swZALhz587UqVOlNo0SBHHixInhw4ejB53q6uozZsyIiYkhCGL37t1ff/01QRDR0dEanTl48ODly5cBQE5Obvr06Tdu3ODPLKLw3o68cG9LN23alJiYKFC4detWAPj77787fUl7e7vA1KiOaVTAo0ePACA4OLj3Acs6nEb7OHTL4O3t3enVtra2vLy8uro6Xkl1dbWurq6XlxdBEFKbRl+8eMH7msVifcQw9PPnz5cuXYq2WQIAW1vbX375RYjD2fy678i3trb6+flNnjxZiNOw+H8+PC4uLmQyuau1Rp6eniQSib+kqzTKZDJfvHhx9epVGxsbW1vb2tpaYYUtu3Aa7eNu374NACEhIZ1eTUtLA4Aff/yRV7JkyRJNTc2KigpCKtNofX39woUL6XR6VlZW71urrKwMCgriDbmoqqp6e3sLpWWe7jvyr169GjlyJABQqdTk5GQhvm9HEyZM0NfX7+rqb7/9tnr1av6SrtLo6dP/nIOkqanZ8Z63f8JptI+LiooCgAsXLnR6taKiYvfu3Q8fPkTfRkdHA99ZRtKWRlNSUszNzVFKun37trCabW9vv3DhAnqKinr6mzdvFlbjaKDGwcGh4xFJV69e1dDQAABTU1PeP4HoTJ8+XUlJqeePL7pKo21tbVVVVX///feiRYsoFMqdO3eEHanswWm0j/v9998B4NChQ++t2dTUZG5uPnHiRN5vmlSl0fDwcNQHHzlyJG9HEuF6+vSpt7c3jUY7efKksNosLy/fuXOnQFcadeRR1p43b15NTY2w3q4bX3zxBQCUl5f3sP57n422t7draGi4uLgIKUAZhjfK6+MsLS0BIC8v77019+zZk5+fr6io+OWXX/r4+Pj4+OTk5BQWFvr4+Lx+/Vr0kXapoaFh0aJFPj4+bW1tnp6eSUlJvB1JhGv48OHh4eFFRUULFizgFSYmJvr5+bm4uMyePXvr1q38x8Tze/PmTWBgoK+vr8De+6qqqoaGhl9++SXvn6CoqGjy5MlhYWEUCiUoKOjKlSvonlTUZsyYAQDnzp3rTSME32odCoWira1dVVXV28j6AEnncSlVWFi4adMmR0dHOzu7qVOnhoaGdjqWyuVyL1686Orqam9vP2HChE2bNqGnitKDy+UOHjx4wIAB790d7ptvvpn+Lk1NTVVV1enTp0vwzN7U1FSUNFVVVc+dOyfOt+ZyuV999RUAmJmZLV261NPT08DAgEajCXRjs7Kyli9fzhut4v2cq6qqtm/frq2tjcrRLAIxd+T5tbe329nZ0en033//nVdYWFiINmqJjo7+6aef+Ot3vBvdvn37ggULeB/wzz//lJOTQ5Mi+jmcRjsRGxtLp9NVVVVXrFixc+dODw8PBQWFiRMndpy/vXz5cgCYOnXq1q1bt23bNmXKFCkcuERryZ2dncvKylBJTU0NWn9SXl6+adOmrmZQS7xTHxkZqaSkBKLsyHcjIiICAL744gtel7ypqWn9+vUCs01PnTq1cuXKGzduoE2jeVkmIyNj3rx5J06cQI+nb9265efnh6Znffrpp+LpyAvIz8+3s7MDAEtLyxkzZtjZ2cnJyS1evJh4d6T++++/Nzc3R4+h1dXVzc3NXV1dCYK4ceOGmpqarq7urFmzxo0bJycnN3To0KqqKvF/EGmD06ig+vp6HR0dExOToqIiXuGjR4/Onz8vUBP1j7rak02qhIWF0Wg0MplsZWVlYWEhLy+vrq5eWlracaSeH38avXLlyoQJE6KiosSzFgiNyKP7OE9Pz+bmZjG8qYAhQ4YMGDAA7UXfE3v27OFPozyxsbEAcOHCBX19fSqVGhYWJuxIPwCbzb5x48bWrVvXrFnzf//3f3/88QeaaFVcXJyWlobqxMfHB72L97C4vLw8JCRkzZo1a9euPX36dM9/OH0bTqOC0MngZ8+efW9NtOua0PdzFJHKysqzZ8/u3bs3NDT05s2b6LedzWZXVFR0laTi4+N5HVgXFxeU1MzMzL7//vu3b9+KLlReR15FRUXMHXmeiooK4Nt/XkBdXd2qVauuXLnCX9h9Go2JiUlISEhJSRFVxMLD4XBycnLwKs+ew2lU0LJlywCgq0U+YWFhLi4uHA6nvb1dUVHR3d29ubk5KSnp1KlTqampYg61J549e+bi4tL77SYbGhrCw8NtbGxQMqVSqZ6enqLY/5jXkR8xYoT4O/I8qampALBr165Or6Ixt4CAAP7C96ZRUcXaO0wmMzQ0dN++fbySpqYmAFBUVBT14q4+A6dRQU5OTjQaraur69at09LSYrPZZWVlADB+/Hg9PT15eXk0wjBv3jwxbzbRvfDwcHSExpYtW4TSIJfL/euvv9zd3clkMsqnDAYjMjJSWD39b775BjW7Zs0ayXYYHz9+DAB79+7tqkJ7e7vAxigymkbr6+sBQFlZmb8QjYx1XG6PdQqnUUHdp1EedMK4oaHhxYsXm5qa2Gx2UFAQAHz//fdiCPK9mEymh4cH79mi0A9xe/Xq1ZYtW3gzdQwMDAIDA3s/2lBUVGRmZiapjjw/9O/75Zdf9vwlMppGCYJQU1MDgOrqal4Jg8EAgEePHkkwKhmC06igpUuXAsB7M0JzczMArF+/nlfC5XINDQ0nTZok4gDfLzMz09bWFgDodPrp06dF90YNDQ0HDx7kHZ2mqKgYERHBu8rlctPT03///fc7d+508yy1uro6Li4uNjb29evXBEFIyQGTXC5XV1d3xIgRPX+J7KZRtIki/1MpV1dXAIiKipJgVDIET78XNHbsWAC4d+9e99WUlJS0tLRevXrFKyGRSHQ6XWD2tfidPHly9OjRmZmZ1tbWycnJn3/+uejeS0VFZc2aNdnZ2YmJie7u7iwWC82nAYBbt24NGjTI3t5+9uzZ06dP19fX9/f3R1MReQiC2LNnj7Gx8dSpU6dNm2ZsbPzll1/yHhdIFolEWrFixdOnTyMjIzutwGazuVyumKMSEXR0c1FRUTclWHckncelTk1NjYaGhqWlpcANKXroWVNTU1JSgkoWLVpEpVJ586LQXsIbNmwQc8A8TCaTlzRF0ZF/L96PIiYmhkwm29vbx8bGNjU1lZeXBwcHq6qqCuy+gaZhnT59urm5mclkogejhw8fFnPYXWlsbBwxYgSZTF6+fHlUVNTVq1d3795taWmZl5fHP8TE5XLRNqZz584FgAMHDqClUARBXLt2LTw8fMOGDQCwbt268PBw6ewm+/r6wrv71/z8888A4OvrK8GoZAhOo524fv06lUrV19ffunVrRETEt99+6+jo+PnnnxMEsXLlSgBAIyrZ2dkqKirGxsY7d+7csWOHhoaGgYEB6pkeOnQoMjJSnIMkYuvI94Stra2Ojo7A3yH+7fiQyspK/hVBHA5HQUEB/ZylBJPJ3Llzp6WlJZlMVlRUHDJkiL+/f3l5eV1d3cqVKy9fvkwQBIfDMe8gNjaWIIjPPvtMoFxKHp0LQOdO8z+hunLlCgDMmTNHglHJEJxGO5eRkbFs2TJLS0sdHZ3hw4d7eXklJSURBHH16tXt27fzhmjT09M/++wzExMTCwuLlStXohvVpqYmdFKYrq7uli1b+Kfxi0hkZCSNRgMAa2vrjIwMUb9d93JyctDNV6dXGxsb09PTO13rlZWVBQA7duwQcYAf40On/hw5cmTDhg0d/3JIJ7TOCi1VQtB8L3t7ewlGJUNwGhV0+/btXk6Xa2trO378ONpHEgAoFMrChQtRFhY6iXfkO0I3MuHh4Z1evXv3LgDwH4VWU1Nz+/btgwcPmpubjx49WqQT+0XE3d3d1taWf3Iuekbc1Vbz0iY5ORkAGAwGrwSduqquri7BqGQITqP/aWlpQduXBQUFCaXBJ0+eeHp6UigUlOZGjBgRHh4uxEyXlZUlPR15HjQmc/369U6vlpSUBAcHP3/+nFeSlJSEJk4pKSl1lXylnLW1NQDw9wP8/f3d3NzS09MlGFXPlZeXA4C2tjZ/ITpaWVZuqCULp9F/8J4t0mi03377TYgtl5eXBwUFGRkZoWSqrq7u5+cnlD2TduzYAQB2dnbZ2dm9b01Y0DFHR44c+aBXtbS0oGGNH374QUSBiY6TkxMA8O+cJFu4XC5aOcZ/jApasSYrfwkkC6dRgiCIkydPKisrA4CVlRVvgwbham1tjYyMHDVqFEqm8vLybm5uDx486E2bLBbr+++/l8i2Hd3IyMiAj103xWAw0KGessXb2xt6tje21Bo8eDAA8B+g4uzsDAA3b96UYFSyor/PG21tbV27dq2Xl1dTU5Onp+eTJ0/4z3MXIiqV6uXllZKS8uTJE29vb3l5+UuXLv3555/oKovFioiImDlzppWVlY2NzeLFi1NSUjo2cvDgwRl8nJ2dHRwc0H2E9LCxsTE3Nz937hxaodAN4t/l2zw0Gq2xsVGU0YlEH5hluWXLlkOHDvF2R4U+8aHEpl+n0RcvXjg4OISFhSkqKoaEhJw8eZJOp4v6TRkMRnh4eH5+fmBgILqLaW1t/d///ufj48PhcDw9PRcsWJCdnT1u3Dh0MhK/xMTE7OxsBp8BAwaIOuAPJScnt3fv3uLi4vnz5/P2ii8rKzt8+DAA5Obmrl69Go1pfPvtt6NHj37x4gWqc/Xq1QcPHkyfPl1CgX88lHGKi4slHcjHW7ly5erVq3V0dHglOI1+AEnfDksML2mKriPfQ2i739DQUF4Jm80+cOBAx/N+HR0d586dK97oPtKvv/6KxiiMjY319PQAwNDQsKKign+kPjc3197enkQimZubo2fHDg4O/Cu7ZUViYiIAjBs3TtKBCNOZM2cAYMGCBZIORAaQCL7DVfqL5mbCz2/148e/ZmR4enr+8ssv6MGoRKC126ampmimXvcGDRo0Y8aMX3/9VQyB9V5dXV1cXFxxcbGamtrAgQMnTJhAoVC4XG5rayuVSkWLPjkcTnx8fGZmJkEQ9vb2n3zyiZyc7PWQSkpKTExMDA0N0b5ffcP9+/cdHR0dHBxQ1wHrRv9Lo9nZ4O4OmZlsc/OzO3Z4LV8u2XBycnKsra39/f2Dg4M7Xi0rK1uyZMmqVavQ5FBlZeVVq1ZNmDChoqLCzMzMyclJQUFB7CFjgjgcjpKSEofDaW5u5h3KJOtKS0uNjY319fXRdCisG/KSDkC8Tp2C1auhsRGsrOSjorxEM5r0Qd68eQMAJiYmnV5lsVj5+fm1tbUA0NTUxOVyjx49eu/evdbW1pcvX1paWkZHRw8cOFCsEWMdkMlkQ0PDoqKi0tLSQYMGSToc4TA0NFRQUKioqGhtbUW71mJdkb0O1EdqbYW1a8HLCxobYckSePIEpCCHAgA646yrvYLMzMyKiorQzhHKysolJSXV1dXPnj3Lycm5d+9ecXGxj4+PWMPFutD3BmTk5OSMjIwIgigpKZF0LNKuf6TRFy/AwZ8QVuYAAA+kSURBVAHCwkBREUJC4NQpEP2IfA8ZGBhAjwd5tbW1edObJk6c6OrqGhsb29bWJsL4sJ7pe2kU+sQMBPGQRKeezYarV8HGBmxt3ylPTYXiYpg//53Cp08hNRUaGkBfH6ZMAQODThp88wYSEzspd3AAExP+jjxERUnJTSiPhYWFrq7u3bt3CYJAd6Y9hw5+YLPZfeZ5nOzqw2m0j30okZDA7ICGBgKA2LlTsHzlSoL/QJjSUsLRkQAgNDUJc3OCSiUoFMLfn2CzBV8YHU0AdPLf2bPEhg3/fO3pSTCZov1cHyswMBAAfvrpJ/5CtPS+vb09OzsbzQGqra09cuQIb9uUmpoaIyOjkSNHij9grCN0oOzy5ct5JfX19SdPnhTuwmIxQzvAfvPNN5IORNpJaxptbiZsbAhtbeKPPwiUOBobia1bCQBi48b3v8WuXYSiIvHmDXHrFkGjEXz70Uqh1tbWadOmAcC0adP27du3f/9+Dw8PZWXlV69evXz5EgC+++474t9fVAaD8c0332zdutXExIRGo92/f1/S4WMEQRAxMTEAMHXqVF5JXl4eAJiamkouqN46duwYACxdulTSgUg7aU2jBw4QAETH9byLFxMUClFa2l37ra2Evj7h4/PPtxUVvQ9Z1NBiUCcnJ0tLSzs7u7lz5x49erSpqammpmbfvn28pff37t3z8PCwtbW1s7Pz9vZ++fKlZMPGeNAuq4MGDeKVtLW1ycnJycvLC+vYVPG7c+cOAEyePFnSgUg7aU2j06cTOjpEx30/Y2MJAAJtp/bXX8S9e520Hx5OkEgE3yYLGCZqzc3NJBJJQUGB/9RltDRLDPt2C8WrV6+WL1++efNmXgnqDA0cOFCCUckEyc0bffQIfv75nZLnz//7OicHzM2h45DL4MEAANnZAABr1oCmJjx8+E4FgoCQEHBxAWtr4ceMYV1QUlLS0dGprKx88+aNoaEhKjQ1NS0rKysqKupqXrBUaW9vP3HihKWl5f79+1GJiYmJiYmJmZmZROOSAZJLoykpkJf3TklV1X9fNzeDllYnr0I70KBtgX75Bf7dEfk/t25BdjYcPizMUDGsB0xNTSsrK4uKinhp1MTE5MGDB0VFRRMnTpRsbD1hamq6atUqtGMegk5slGBIskJyadTXF/bseadk1So4f/6fr9XUoNMlaK9fAwBoaAAATJ3aSYUff4SRI2HyZCFGimE9YWpqmpKSUlRUNG7cOF4JyM6EIRqNhoYxeRoaGsLDw+/cuVNZWUmn06dMmbJ+/Xp0VAG/169f//jjj0+ePGlra7Oyslq/fv3w4cPFGLjkSev0e3t7yM+H9nbB8qysf652KjUV4uPB31+0sWFYZzpOVpetNCqgoKBg+PDh27Zto9Pps2bNsra2PnDgwPDhw9HyZZ6qqqqxY8deuHDB2dnZw8MjOzvbwcHhocCjtr5OWtfUL1wI16/DyZOwatV/hQQBYWFAp8Ps2Z2/6scfwcgI3N3FEyOG8euYNGU6jXp5eZWXl8fFxfGeSAQGBl64cAFte8hz7dq1kpKSuLi4KVOmAMDy5cv19fUPHz7MuyXvD6T1bnThQpg8Gb7+Gn75BdB26IWFsGQJ3LkDQUGgpgYAsHIlrFv330tKS+HSJVi3rpMHphgmemgcqWMalcXFlE+fPk1KSlq9ejX/U10jI6MNGzYIrLV7+/YtAKiqqqJv6XQ6jUarq6sTZ7QSJ61pVE4Obt6EBQvA1xdUVEBJCQYOhD//hIgIWLPmnzq5uZCf/99LQkJAUfGdu1cME6Nu7kYJ6d6OsrCwkMPh8Jfcv38fALo6ieD27dtz585FxxY4OTmRyWQvL68HDx4AwNmzZ2tqatCxDv2IBCZZcThEenons+JLS4nMTMHCykoiOpqIiiLu3yfa2rprNjOTePFCmHFi2IeoqakBABUVFf5CTU1NAKisrJRUVO91+fJldXX1Xbt28ReiExkyO/4+EgRBECdOnNDQ0EhJSUHfHj9+HCUTKysrCoUi0+tfP44kno3KyYGdXSflRkbw7ynE/9HRgZkze9SsjU1vA8OwXtDQ0FBVVW1oaKitreUNZ69bt45EIqGt/qUNOs8xIiICAHJzc4keb46zbNmyZcuWoa9LS0uDgoIcHR19fHwuXbpUUFCwceNGU1PTyf1qtoyk8ziG9R1//vnn48eP29vbJR3I+7148QJNS0LnOQpcDQ0NBYDo6Oj3tjN//nwtLa36+npesxYWFnp6em3d9x37Fml9NophYsfhcFJSUk6dOhUZGZmcnEx0/UCztLT0r7/+amlpESgfOHBgTk7O+fPnecedSqcrV644ODg8e/Zs8ODBDx8+XLt2rUAFNM6ekJDw3qYePnzIYDB4Q0yWlpbLli2rqKjgHQrbL0g6j2OYVEhJSbG2tgYAKpVKo9EAgMFg5OfnC1R7+vSpp6cnhUIBgIcPH/Jf2rhxI4lEUlZWptPpJBLp66+/5nbcFELSWlpa/Pz80O++q6trbW0t/1UOh5OQkIC+ZjAYdDo9IyODvwK60W5qaiorK0Nf29vbDxgwAO3riLi5uVEolLq6OpF/GKmB0yiGEQUFBWpqakZGRtHR0Ww2myCIpKSkIUOGRERE8Fe7deuWkpLSnDlzVq1aJZBGz549CwDr1q1ra2trb2/fvHkzABw/flzcn6Rb3Xfkq6qqnJ2d5eTkYmJiCIJ4/vy5rq6usrLymjVrjh49GhwcPH/+fDs7O4IgfvnlFwBAG4+dOXOGRCKNHz/+1KlT165dQz+ZjT3ZzbIPwWkUw4iVK1eSSCTe0DPS8eleW1tbY2MjQRDXrl0TSKNjx47V0dHhPRXlcDiGhoZStan25cuXVVRUAMDS0vLp06cCV+Pj49FWADo6On/99RcqLCsrW79+/dChQ7W0tKytrefOnXvx4kUOh/Po0aOtW7cWFxfzXjt79uyBAwdqa2s7ODiEh4fzb3PVH+A0imGErq7u6NGju7p69+7d9PR0/hKBNMrhcCgUyrJly/jr+Pj4yMnJSc9Iy927d8lkcseOPJfLDQkJQY8pJk2aVNr9Zr5YZ/AQE9bf1dbWVlZW8u9sJGD27Nm7d+/upoXKykoWi2VsbMxfaGxszOVye3jIOyH6+fmTJ09+9OgRmiXKK6yurnZxcVm3bh2bzfbz87tz545Rx0mH2PvgNIr1d01NTQCgo6PTVYVjx451HMvmhybe80arEZStqqur3xtAUlKSlZXV/v37Rb2GksFg8H+bkJAwfPjw6OhobW3tW7duhYaGUvBC6o+C0yjW36EnhhUVFV1VWLhwYfcbhqIULJAE0WJzPT09LpfLYrG6efnp06dzc3MDAgLMzMzWr1//6tWrD4r/IxAEERoaOn369LKyskmTJj179szZ2VnUb9qH4TSK9XdojB4dpvRxtLW1O+5wXFxcTCaT9fX1r1+/bmJiEhAQUFZW1unLDx8+fOPGjenTpzc0NISEhFhZWc2YMePmzZsi6unjjrzwSfjZLIZJAV9fXwBITEzsYf2OI/WffPKJhoZGa2sr+pbFYunq6o4fP54giBUrVqDfNSqV6uXlJTAfgF92drafnx+atQoAFhYWQUFBNTU1vfhkgu7du4eSpra29h9//CHElvsznEYxjHj9+rW2trauru6lS5dQKszPz9+7dy+a0/P555//8MMPqGZ9fX1NTc2ZM2cAICYmpqampqWlhSCI69evA8CiRYuYTGZTU9NXX30FAJcvX0avSkxMdHd3l5f/ZwsLBoMRHh6OXthRdXV1UFAQ2h0KAFRVVdetW9f7z8g/Ij9x4kQ8Ii9EOI1iGEEQxPPnz0eNGgUAZDIZ3Q/q6enFxsYSBKGvr79kyRJUbejQoQL9ub1796JL+/btk5eXJ5PJ8vLy8vLy3333ncBblJWVBQYGav17yJient6WLVtKSko6jYfD4fz111+zZ88mkUiurq69/HRoaj0AkEgkPz8/mVj1L0NIhHTvhIhh4pSZmYmWw5uYmNjb2ysoKAhUyMrKElhKb2RkpK+vj76uqKh4/PgxAIwePZpXKKCtre3ChQs//vhjeno6ACgoKHz66afe3t5dbe6ZkZFBJpNterGBWUJCgoeHR1lZmba2dmRk5KxZsz66KaxTOI1imGSkpqaGhoaeO3eOzWYDAIPB8Pb29vT0VFJSEtZbEAQRFha2adMmFos1ZsyYCxcu4NOSRQGnUQyTpPLy8vDw8EOHDqEZprq6usuXL1+zZo3AZP6PExYWtnbtWhKJFBAQsGfPHt7DWUy4cBrFMMlra2u7cePGTz/9lJycDAAUCmXevHnd9PR7qLW1dc6cOf7+/jN7uPc59lFwGsUwKZKYmHjgwIGrV6+inn5wcLD/vweG19XVxcfHl5SUKCkpOTg42HV6hARAbm5uQkJCa2vrqFGjHBwcxBd6fybJ8S0MwzpTXl4eGBior6/P2/A0NDQUbWOqr69Pp9MBwM3Nrbm5WeCF27ZtI5FIdDpdWVkZABYsWIAH5cUAp1EMk1Jo51Pi3zPjXFxc0CGjbDY7Kipq8ODBz58/569/8eJFAPDx8WltbWWz2T/99BMABAcHSyD0fgZ36jFMqrHZbDMzMyUlpYyMDEVFRV45l8uVk3tnMfe8efPu3r1bUVHBqzZs2LD29vbs7GyxRtz/4DX1GCbVnj17VlZWtmTJEv4cCgAoh1ZXV1+/fr20tBQAKisrdXV1+atNmDDhxYsXDQ0NYo65v8FpFMOk2suXLwHAwsKi06vPnj2bN29efHw8AAwZMqSoqCg/Px9dYrPZ6ECkbjavwoQCzyPDMKnW3NwMAJqamp1eHTZs2JUrV9Ay1o0bN0ZFRU2cOHHZsmVtbW23b99Gx3NyuVwxxtsf4btRDJNqampqANDVLvo6Ojrz589Hc/VtbGwePnz4ySef3Lx5Mzs7e/369UuXLgUA3ip+TETw3SiGSTVbW1sA6OF2qHZ2duiMUuTUqVOGhoba2tqiCg4DAHw3imFSztra2tra+uTJk7W1tR/0wszMzISEBFdXVxEFhvHgNIphUo1EIv3888/V1dXTpk2Lj49vb29vbGy8d+/e5s2bAeDZs2czZ86Mi4tDlWNiYmpraxsaGmJiYtzc3DQ1Nbdt2ybR8PsF3KnHMGnn5OR09erVr7/+esqUKbzCcePG1dbWMpnMZ8+eoRvVgoICNzc3JpOJKowaNSoqKsrAwEAyQfcnePo9hskGLpebnp5eUlKirq5uZmbW6RZQjY2NaWlpTCbTwsKiqzlSmNDhNIphGNYr+NkohmFYr+A0imEY1is4jWIYhvUKTqMYhmG9gtMohmFYr+A0imEY1is4jWIYhvXK/wPwMTPOT2js5gAAAUV6VFh0cmRraXRQS0wgcmRraXQgMjAyMi4wOS40AAB4nHu/b+09BiDgZYAARiDmh+IGRjaGBJAYMxsHiGZikmCEs5jgLIQsC5zFCmexcYBZjIwS7GweGUAWMzOjBAeCyYlgciGY3AgmD5uHBojJIsEryKEAMkqCjxvoTkYmBiZmoDADCysDKxsDGzsDB2cGDydXBg8XdwYPN08GDw9vBg8vHwMvRwYPEzeDCNDtvDzcXJzi+6D+BAP+BNvz++uPqTiAOGrepfusNf7ZQ9l7gWw7ENvTfuH+gr3nQBoZFr3UOKC72hYsLrWl9AA3h/l+EJvn/7IDsr8YwGo8Us4eWB/RDBZPa2lyMFuyCMwOsbNwuPiO6wCIvfbFfvv7Faf2Q8xcYa+7eu9eiF2KDgV7/cBucOusdNh8Ph5s133NxQ6Jxv1gcTEAO0ZMW8D/z5wAAAG8elRYdE1PTCByZGtpdCAyMDIyLjA5LjQAAHicfZNBbtwwDEX3PoUuEIGfEiVymZkJiiLIDJBOc4cAXfb+KOmBLM2msmVY8hMp/U9vKdrn5f37bzoaX7YtJfrPbWbpqxDR9pHiJZ3efvy8pvP99TRmzrff1/uvBImb4npmX++3jzGDdE4vyA1GjRJn6bXFqkx7m0s5QMpcTGpNyKbSqEwQAywPkMRUm3+v1iphgjzA+kjNYt0j+pJKpatNsgxSguTcirFvkjJUmJfcdYAtwJK1GNAjZIenb5OUQfYgJTO0lRIkfB9SJ9kGqemWXlqu0GqIDRM31IXsgzSPWTN1bY33k3VqZUmuh+jukWutXZVqnAydlBfSDjL8gWtorBoxRaFYkuMwCLyjBdZ8FA64Wj450cMilD2/wKpHpazUpetCHh4hTHJBzYgiq+slq544PPKiOf9J4lnddvhGAO2yVBIOl96ul6cafFTl6Xa9zKqMi2ft+WCWQ/FhmaZ7NXqXWSzw3mZJwHufvsO7TnOr2zZFtIj2JNUeH6skkc11XY7OMVOWE8YgtIMsQtZ42KJXnGLGDVVWDWI8/mp/3/4BWtzH+Ct23uIAAAFTelRYdFNNSUxFUyByZGtpdCAyMDIyLjA5LjQAAHicbZK9bkMhDIVfpWMiEcQxxth0qrJk6c9+xZBG3RIlitotD19zM967YPFxOMcI76f9gRq6l0Z9M+9ynws/S3kW6dPnodW+nU4N6JhOh4bUR7F51T5O8vA4N/AQuoI6Xh6bHaLAkgSKpbJQeN2lSNkKc0A0LTKTVExVEFJkE05wiEjFKgc/5ZSrmjOKko3EZdBClB3lqNlAQ1bhl4dfiQSVPBjcpLAziQxlC+6bSOAIMcOequSmgAzGZjSrikK9D4paVYNHoyb1RI6pqseMBmsSN/IWzLwnDy0jv7ifv8bfB2iV4VFgPCSaaqkatuH4e7183a+3hni5nt989368ffxdvn/uEY2W0D9nCXPjJeRWlrA0WUJpdQlr05WW0GyFJh+EJTYfkyVVH4kVi+yTs4LZB2kF0+MfGzDARADIsjoAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 15
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "editable = Chem.EditableMol(combinemol)\n",
        "editable.AddBond(0,7,order=Chem.rdchem.BondType.DOUBLE)\n",
        "editable.GetMol()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "nm8zIUIg9kgt",
        "outputId": "915d7993-61e0-4bd5-841a-e415c0c74a65"
      },
      "execution_count": 34,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c74f6eb0>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO2dd1hT1xvH3ySQEDYKskSGyHACKi7EAe5dFwpuwYFbq1b9ldo6oGpddaA4i0JFca/iQEWoynIgqOwRtuwRMs7vj9OmMQyRkRvgfJ48PuS9JzffC/HNuee8g4YQAgKBQCA0FDrVAggEAqFlQ9wogUAgNAriRgkEAqFREDdKIBAIjYK4UQKBQGgUclQLaKPw+XwOh1NWVqarq6uurv7V8dnZ2QUFBdra2hoaGlKQRyAQ6g+ZjUobDoezePFiTU1NQ0PDrl27tm/f3s7OLiwsrLbxPj4+ZmZmOjo6lpaW7dq1s7a2jomJkaZgAoFQN2Q2KlUSExPt7e0zMjIcHBwmTZrUrl272NhYb2/voUOHPn/+vE+fPhLjN27cuGfPHktLy4MHD+rp6eXk5Dx//tzU1JQS8QQCoUZoJPxemgwdOvTJkyc7duzYunWryMjhcHx9fTds2ECnf3Fz8OzZsyFDhvTt2/fx48eKiopSF0sgEOoFcaPS4/Xr11ZWVlZWVpGRkTQa7avjZ82a5e/v/+zZMzs7OynIIxAIDYOsjUqPx48fA8DkyZNr86EVFRW//vrrxYsXRePZbPbAgQMBoLCwsLS0VGpSCQRC/SFuVHokJiYCgJGRUW0DSktLN23adOTIEQAoLCzMzs42NTU9deqUubm5hoaGiopKr169/vrrL6kJJhAI9YFsMUmPsrIyANDR0altgIqKio+Pj7a2NgAUFxcDwPv37zdt2jRnzpx169YlJiYeO3Zs3LhxwcHBgwYNkppsAoFQN8SNSg9lZWUAqKioqG2AgoLCokWL8M/4xt/IyOjvv//W1NTERgcHh1GjRu3atev27dvNr5dAINQLclMvPUxMTAAgNja2PoPbtWtHo9FYLJbIhwKAo6OjoqJidHR0c0kkEAjfDnGj0mPIkCEAcPny5fpERygpKRkaGiYlJZWXl4uMfD5fIBAoKCg0o0oCgfCNEDcqPaysrAYPHhwZGfnbb79JHOLxeACAEMrKyvr8+TM2Tpw4saKi4ty5c6Jhf/zxB5fLxe6YQCDICCRutHnhcrkBAQGFhYUrVqwAgE+fPg0ePDg7O3vixIlTp07V0dFJS0t7+PBheHh4XFxcfn5+hw4dBg4c+Pz5cwDIyMjo2bNnRUXFihUrrK2tX79+vX//fjabHR4ebmpqGhYW9uDBg8WLF+vq6lJ9lQRC2wYRmof09HQPDw8tLS0AUFNTKy0txfbk5GQnJyfxG3N1dfXFixeXlpYWFhba29u7u7uLThIdHW1raysa2b9//+joaHxo+vTpAMBgMMaPHx8UFCQUCim4SAKBgBCZjTYxQqHw/v37R48evXPnjlAoBABra+tly5bNmTNH3HWWl5cnJiZyuVx9ff06QqAwmZmZHA5HYuTjx49///33Gzdu8Pl8AOjWrRt+F1VV1ea5MgKBUDPEjTYZhYWF586dO3ToEA6zZzKZkyZNcnNzc3R0bL43zczMPH/+/JEjR9LS0gBAWVl59uzZy5cv79WrV/O9KYFAEIe40SYgIiLixIkTvr6+eFddX19/8eLF7u7u+I5eCvB4vGvXrh09ejQ4OBhb7OzsDq1daz1+PDCZ0tFAILRZiBttOJWVlTdv3jxw4EBoaCgA0On04cOHu7m5TZkyRU6OmryGjx8/nj59+uTJk58/fy61tFTKz4cFC2DpUqg9A5VAIDSSVu5Gc3Nzg4ODMzIyGAyGhYXF0KFD5eXlqw/Lz89/9OiRhLFr167dunWr8bSfPn06deqUj49Pfn4+AKirq8+dO3f16tU4wJ5ySktLgy5dmnLwILx5AwDAYMD48bBsGYwYAfRaQtwEAnj5El6/huJi0NSEfv2glmsHAHj0CPLza7AzmTBpEoSEQFYWTJ4MEt8leXkQHAwWFtC9e4MvjUCQRSjd4GpGeDzeunXrmF/e0hobGz979qz64Fu3blX/zRw+fFhimEAgCAoKGj9+vKhEU+/evb29vcvKyqRyTd9OeDhyc0NsNgJAAKhzZ+TpiXJzJYc9fYrMzf8ZI3rY26OEhJpPO3Cg5GD8aN8eIYRGjUIAqLBQ8lVPniAAtHFjU18kgUAxrdaNOjs7A4CNjc3t27ezs7M5HM758+f19PSsrKz4fL7EYB8fHwA4f/58ghglJSWiAVlZWZ6enoaGhth7KigoTJ8+PTQ0VLrX1FCys5GnJzI2/sfZsVho+nQUFPTP0dBQpKCA1NXRkSMoIwNVVaGPH9HmzYjBQDo6KDOzhhNmZKCEhC8eoaGIwUATJiBE3CihzdE63WhQUBAAdOvWTRStiUlMTMzPz68+fvv27QDw7t276ofevXs3Y8YM0VKAubn5gQMHCgoKmkt688Hno+vX0ejRiE5HAEheHmVnI6EQde2K5ORQWJjk+MOHEQBycanXybduRQDoyROEmsyNPnz4cPny5cOHDx84cKCTk9PFixcFAkFtg+Pi4pYvX+7q6lpcXCxxqKSk5NChQzNnznz48GE937qZePTokbu7O76imTNn+vr61nFFCQkJq1evXrRoUXZ2tri9qqrq4sWLM2bMsLOzmz59upeXV/VLJkiZ1ulGXVxc4N/s9fqwdOlSAKjRw0ZGRgIAnU53dHS8ceNGa4hy//QJbdiAlixBCKGQEASAZsyoYZhQiExMEJNZg0OUoKwMaWqi3r3/edpoN1pRUYGTC+Tk5Pr06TN06NAOHToAwLBhw8RvETDPnj2bPn06g8HA33McDkd0KDMz08PDo127dvjQrl27vvrWzURlZaWTk5P4FeFyiPb29tWdYHh4+Jw5c0S7lO/fvxcdKi4uHjZsGF6eGj16NF6L79atW+FX/0aE5qR1ulHc9C0rK6u2AX5+fvb29jdv3sRPJ02axGKxPn36dOHChdOnT4eGhoq7y99//z09Pb3ZRVOClxcCQCdO1Hx0yRIEgP76CyGE5s1D9vao2noIQgj9/jsCQH5+/zxttBvF32oODg4ZGRnYwufzPT09bWxsJP4Qubm5CgoKgwcP3rt378iRIyXc6JQpU4yMjFavXv3LL79Q60ZxKvDQoUPT0tKwhc/n79mzx8rKSmTBlJaWqqqq9u/ff/fu3VOmTJFwox4eHgDw22+/4c8nn8+fO3cuAHh4eEjxagiStE43qqSkpKCgUMfM8ddffwWAE/+6j759++JvftHeUf/+/Vut6xRn1SoEgO7fr/no7t0IAJ05gxBCXbsiAFRVJTlGIECmpqhjx/8OYTc6cCAaPPiLR8+e9XGjSUlJdDpdX1+/+jSt+qI2Qkg0P8W3IOJuVHQIL/JQ5UZTU1MZDIaurm5RUZHEobqvCDtfcTdaWVl569Yt8cEfP34EADs7u6ZWTfgGWmeFJx6Pp6amVkfbOHd398zMTPwfDwD++OMPX19fDocjFArT09NXrFjx999/z5w5U1p6qYPHAwDJyCQReEW4qgoA4OlTyMyE6uFi165BfDysXSt5SEEB2OwvHixWfRThlZMFCxaoqKhIHMJ37gKBIC8vr6SkBBtxMewaqeOQNLl586ZAIJg3b171PF18RUKhMC8vD/c7gDpls1iscePGiVvw4gBurECgitZZ/V5dXb2oqEggEIiWzCRQVFQUb1lsbm5ubm6Of9bX1z906FB0dHRISMjbt2979OghDcVUoaEBAPBvaT5JcHBo+/b//VudfftARQX+Ldr/H4GBoKb2heXpU6hHib+3b98CQG0RuwDw4cOHbt26TZo06dq1a189W218/vyZw+E0+OX1RFFR0cTE5KtXlJqaamxsPGLEiAY02sJVwM3MzBqjk9BIWqcb7d69+6NHj969e9ew1HIajda3b9+QkJDExMRW7kZ79gQAiIiAadNqOPrqFQBAHb/DV68gNBTWr5f0mI0Al1vt1KlTbQMUFRUHDBhgYWHRmHe5cOHCqlWrGnOG+mBra/vixYuvXpGCgsKAAQPq8LN1gBsgzp49u8EiCY2ndbrR8ePHP3r06MSJE/hD1gAyMjIAAG8Qt2ZGjgRFRbh4EX78EdjsLw7Fx0NwMPTsCaamtb78119BTg6a1B/hHeo6smmNjIxw9m1jaN++fcPc1jfRuXNnqMcV6ejoNOyKrl275uvr6+DgMGHChAaLJDSe1ulGFy5cuG/fPm9vbysrK1dXV5Gdw+FUVlbi+6wXL14MHjzY3Ny8sLDQxcVl48aN9vb2eNjTp0+vXr1qZGQk2npqtWhowJo1sGsXLFoEp07950mzssDJCfh82L79H8uVK1BQAIsWgWjFOSkJrl6FGTOg9nlWA8D1XNLS0sQLrTY5s2fPltoMTnRFTXvakJCQOXPmdOnSxdfXt45tAIIUaJ1bTGpqaoGBgRoaGm5ubt27d1+yZAkO5DY0NPzpp58A4N69e66urk+fPgWAhISEsLCwoUOHDhgwYPHixWPGjBk+fDiLxfL19aWqwohU8fCAqVPBzw+6dAFXV9i8GWbPBnNziIiAn3+GyZP/Gfbjj+DqCnz+fy/cvx8EAlizpmnl9OnTBwBEpapaAc1xRX/++eeoUaO0tbUfPXr01Xq1hOamdbpRALC1tX3//v3WrVvZbPaNGzfu3bsnEAhWr169detWABg6dKiXl1e/fv0AoHfv3nFxcVu3bhUIBHfu3ElPT1+2bNm7d+/aSi94JhMCAuDSJejZE27fhiNHIDQUxo6F58/hf//7b9j69eDlBaItOy4X3rwBZ2eoPme0tIT+/WvY/VdVhf794d+E2tqYPHmykpLSuXPnUlJSGnVdMsPEiROVlZX/+OOPpKSkxp+Nz+dv3rzZycnJxsYmJCREX1+/8eckNBaqI64IBElwyz8TE5P79+/jyEoul3vz5s1ff/0VIZSamjpmzJiff/5Z4lXV40ZFUBs3ihA6dOgQABgZGd29e1d0Rbdv3/b09EQIZWdnjxkzZuvWrRKvqh43mpOTM3z4cABwc3Orqh7DS6AI4kbR2rVrN27cmFljDY42i1CI/P3RxImospKS99+9ezfuuaKoqGhgYIAD1/T19Xk8XkxMDABMmjQJj/yf+JT5X6ZNm4YQio+Pr3HqUD2jVAp4eXmx2WyJK9LV1eVyuXiWOmLECDxy79691TU7OjoihHAmaI1Uj+0nSI1WXm/0q+Tl5RkYGFRVVcXHxxsbG1MtR5bo0wciIuD8eZgzh5L353A4N27ciIuLEwgEHTt2tLCwGDlyJJvN5vF4KSkpysrKeE3w2bNn1be5zc3NJ0+eXFhY6O3tXf3Ma9euZVLRFCAzMxNfEZ/P79ixo7m5+ciRIxUVFfl8fnJysqKiop6eHgC8evWqevVbIyMjXM0Ex5BUZ82aNaz6JTgQmpy27ka9vLw2b948bty4GkuOtiEQgqIiUFf/z+LjA66uMGAANDq6SPo8fPjQ1ta2eh4UgdActNotpvogFArxbGX58uVUa6GUv/8GCwvJTCRnZ2jXDsLCICKCIlkNJD09fdKkSebm5vk1lugnEJqaNu1G7969m5SUZGhoOGrUKKq1UIqRESQlwY0bkJ7+n5HNhrlzAQBqui+WZTZu3FhWVmZvb9++tgRWGaOkpCQrK4tqFYSG06bd6LFjxwBg+fLltaXetxV0dGDKFODzwcfnC7u7O9BocOECFBRQpOybCQsL8/f3Z7PZnp6e4naZXbx6+vSpmZlZW78fauG0XTeakpJy7949Fos1f/58qrXIAMuWAQCcOPFPzSeMqSk4OEB5OZw/T5Wub0IoFK5ZswYh9P333xuJNUNFCDk6Ov7yyy8VFRXUqasZc3Pz8vLyq1evNqAuCUFWoDROgEo2bdoEAHPnzqVaiMzQvTsCQJcufWG8cgUBoC5dUEso+3/q1CkA0NfXl2gec+7cOQDQ1taWzaignTt3AkDXrl15PB7VWggNoY260crKSlx2JKx6D6I2C26+NHToF0YeD3XsiADQgwcUyaovxcXFurq6AHDhwgVxe0lJCQ4kOnfuHFXaRJSXl//8889nz54VN3K53C5dugDAkSNHqBJGaAxt1I2eP38eAKysrKgWIksUFyNVVQSA3r79wr59OwJAU6dSJKu+4NuLAQMGSHQ92LJlCwD07t27jv5xUiMwMBAAtLS0JLoiXrlyBQDatWuXl5dHlTZCg2mjbnTAgAEAcPLkSaqFyBhLlyIAtGLFF8bMTCQvj+Tk0Jddg2SKhIQEFotFp9NfvHghbk9MTFRQUKDRaM+ePaNKmwRDhgwBgHXr1knYcS+pVatWUaKK0BjaohuNjo4GAHV1dYkVNAKKiUE0GlJRQRJriDNmIAAkw33TJk2aBAALFiyQsH/33XeytgIeFRXFYDDk5eXj4uLE7TExMXJycnJycm8l7gYIMk9bdKO4Aine0iVIMngwAkDHjn1hfPwYASBd3Rpa2skADx8+BAAVFRWJuiQ4pVJRUTE1NZUqbTXi5uYGYkn0IpYtWwYADg4OlKgiNJg250YLCwuVlJRoNJrEXIDwDxcvIgDUrZukvcZ9fBmAz+fjRi+4WpK4HbeQ2bFjB1XaaiMnJ0ddXR0A7ty5I27Pz8/HKQM3btygShuhAbQ5N3rgwAH4t14OoQa4XKStjQCQxGIi3scfNowiWbVy+PBhADAxMan8shjV0aNHAcDY2LiiooIqbXWwb98+ADA1NeVyueJ2/Pns3LlzJUW1tQgNoM25UdyB58qVK1QLkWG2bkUAaPbsL4y17eNTyufPnzU1NQHg6tWrNdovX75Mlba6qaqqws1of/vtN3E7j8fr3r07AODiqoQWQdtyow8ePAAAPT09UvK2LlJTEYOBmEyUlfWFHe/jr1xJkawawIWNhw8fLmFfvXo1AAyTvbmzOLiomKqqqkStW/wprb7US5BZ2pYbnTp1KgBs376daiEyz6RJCADt3PmF8fVrvq7uvVGjKCl7XJ3379/Ly8szGIw3b96I22NjY7E9OjqaKm31ZMyYMQCwZMkSCfv48eMBYPHixZSoInwrbciNcjgceXl5OTm59PR0qrXIPPfuIYASOzvc8ULE0MGDAeD48eNU6RIH1+Vyd3eXsGPftHTpUkpUfRPY49Pp9PDwcHF7fHw8DoN9+fIlVdoI9acNuVEPDw8AmDFjBtVCWgJC4c4xYwDg+vXr4uYLFy4AQM+ePanSJeLatWsAoKGhkZubK26/efMmDgrOycmhSts3sWbNGgCws7OTSL7asGEDAAwcOFDYEqoZtHHaihvl8Xi4h+Ljx4+p1tIywB2BRo8eLW7kcrna2toAEBISQpUwLMPMzAwADh06JGHH+zb79++nStu3UlBQgBvZBwQEiNtFJQL8/Pyo0kaoJ23Fjf75558AYGlpSb7b60lBQYGioiKNRvvw4YO4HaeoOzs7UyUMIYRriVpaWkpsFe7ZswcALCwsWtYWIq57a2BgUFZWJm4/efIkAHTs2JGk28k4bcWNDh06FAAOHz5MtZCWxIIFCwBg/fr14saUlBQGg8FkMrMk9vGlRXZ2tpqaGgDcu3evRrtETLvsU1umgEAg6NOnDwD89NNPVGkj1Ic24Ubfv39Po9GUlZULCwup1tKSiIyMxOuMErOkiRMnAnVt3xcuXAgAEydOlLDjHN9x48ZRoqqRPH78GOetpqSkiNufP39Oo9HYbHZycjJV2ghfpU24URxd2CK2bmWNvn37AsDp06fFjXfv3gWATp06SezjS4HIyEg6nc5kMiWWGmqr99GCwNF4c+bMkbDPnDkTAGbNmkWJKkJ9aP1utKSkBN/rRUZGUq2l5XHmzBkAsLa2FjcKhUK8wyOxj9/cCIXCwYMHA8DGjRslDuHqcxLrDy2LlJQUvBgtUdMvLS1NSUkJAJ48eUKVNkLdtH43GhUVZWRkNGjQIKqFtEjKy8txsQyJAMYa9/Gbm9LS0lmzZuno6Ej0AvH39weADh06SNRCbnFs3bq1xgrTOFbP2tpaFipPE6rT+t0oQkggEJC8ugaDAxglSnnWto8vBfLz88WflpeX4+51J06ckLKSJqesrMzAwKB6v5Py8nJDQ0MAOHXqFFXavglOFSetSnaLfDc5NCSrjWfFycnJ8fPze/XqVUFBQfv27QcPHuzs7KyoqFj3q2JiYg4dOgQANjY2S5YskYrSVkhCQoKZmRmLxUpLSxPv/L5w4cIzZ85s2LABhxmJCAkJuXnzZnx8PI/HMzExmTZtmp2dXW0nz8jIOHbsmEAg+PHHH9lstvihtLQ0Hx+f169fMxiM3r17u7q64vhKCX7++WcPDw8rK6vw8PCW3ig7Ly9v7dq1vr6+LBZr8uTJQ4cOdXFxUVZWBgA/P7/Zs2c7OTn5+fnhwTweLyAgICQkxMXFZeDAgeLnSU5O9vHxiY2NVVBQMDc3d3Nz09HRaQ7BFcKKyPLIPEGeOkPdim2lxlDDdps4myxeFqcHp7YXfuZ/5vA45grm8jT5GgcgQFHlUQgQm87uqtC1OcQ3JVT78a/j5+eHP0mampo2Nja4FZ2xsXHdmwl8Pr9fv374hdOmTZOa2lbJ6NGjAWDv3r3iRryPr6GhIdrHLysrmzZtGgDQ6XQzM7MePXqwWCwAcHNzq343GhUV5ebmpqCggD+HEslIoaGhqqqqANCzZ09LS0sA0NLSevfuncRJSktLNTQ0ACA4OLipL1raXLlyBV+ynJwcAOD10E6dOuGrFgqFQUFBeGRRUdGBAwc6deqEf3W//PKL+HkuXLjAYrGUlJRsbW3xr05RUVEiOKzxlAhK3FPdFaIUIALwgxnJdE91LxeUI4SsY6113+jW+MKEyoRVaauUopQgAi4X1Fp/62TuSXzanu+pT5n7KrLuRp8+fcpgMNTU1K5evYr/KwqFwsDAwA4dOgQGBtbxwv379wPAb7/9Rtxo47l+/ToAdO7cWcIb4n38M2fO4KfOzs4AMGLECFHUTkFBwYIFC/r06SPRqS0zM5NOp/fo0WPbtm14riruRisrK42MjOTl5UX/+f38/Oh0upWVVXV3/OnTJwn/3hL5+++/5eTkVFRULl269PLlSxyN4O3tra2tXT2L6bvvvtPQ0HBxcdm8ebOEG01OTmYymZMmTSouLsaW+/fvy8nJGRsbN6HaSmFlv7h+EAGDPgy6kH/hRemL24W3ZyXOgghYl7YO1e5GZyfNpkfQlaKUzGPM63CjmVWZGtEa/eL6tX/dnrjRJgD3nrt586aEvby8vI5XJScnKysr29vb5+fnEzfaeAQCAV5/lJjUnD59GgBsbW0RQlFRUQBgZmYmUSZZKBTWWDhZtFqNA3rE3eilS5cAYP78+eLjx40b1zpmnTUybNgwEMsHnTNnDgBMnTq1xs95Tk4O7miPu6RIzEb/+usvifjonj17AoDEvlxj2JG5AyJgUvwkvvCLiLcrBVfKBGWodjd6MvfkzcKbFYKKIzlH6nCj3yV8JxcpF10erfVaq0W4UXpzLRY0BRkZGWFhYRYWFrhumDiidTSEkFAolDi6ZMmSqqqq48eP02g0aQht7dDpdBzcjtMWRTg5ObVv3/7ly5fh4eGXL18GAHd3d9F9OoZGo4ksOBMX/4wTxmskKCgIALB7FTF9+nQAwLU4Wxn5+fnBwcHGxsY4dBQAPD09HRwcvv/++xo/51paWvjGv0ZGjBiBI/wwlZWVaWlpmpqaKioqTSXYO8+bQWP83ul3Bu2Lxejv1L9TpNe1Y7FYc/F4tfEKdIU6xtwsuhlYGLhRe2Mvdq+mkdv8yLQbxS08cZ5cbVhbWzMYjIqKCpHl/Pnz9+/f/+GHH/DCEKFJcHV1ZbFYt27dSk5OFhnZbPb8+fMBYPfu3V/9Y4WGhjIYDJxgWjcfPnwAgC5duogbcaQqPtTKwHVRe/bsKfrW19PTe/DgQb9+/URj+vTpw2AwSkpK6n9aoVD44cOH2bNnFxYW7tq1i0ajCYVCTp1kFWVxeJw6HnzET6lKSatK66rQtaN8x/rI2JSxaVnqsnpqLhIULUtdZsQ02qKzpf5XSjm1fqfJAviWHIeA1IaxsXFlZSWd/s/3QV5e3oYNG7p06YKXjQhNhZaW1tSpUx88ePDx40d8g4+ZN2/evn37rl+/jrO/Rfse1VFUVDQ3N69jEiqiqKgIv6O4EW8tFhYWNvQKZJd6fs7LysrqH4rg4eHh5eXF5XIZDMbp06fxt11GRkYdfyAAGHJ6yJNeT+oYEGkRWSosBQAjplE9lfgX+Gfxso51Ovb1oQAbMjZk8DJudb6lRFeq5/llAZl2o9g54t3e2rh69ar405UrV+bm5vr5+UncWhIaz/79+9XU1CT+HC9fvsQ/4JlUHX8sKyuruLi4+rwRPhX6MhRPIBDAvx+JVkZ9Pud4zaT+jB07VldXl8Ph+Pn5ubq6CoXChQsXMhiMur/GVOVUdeXrGiBPk+ciLgC0l2tfxzBxAowDqlBVfUY+KX1yKu/ULI1Z49TG1fPkMoJMu1E8H+Fwao0+k+DOnTv+/v5z5851cHBoTl1tFDwflABv4vfr1w+3kMvIyNDT02vkG+EYppycHPE1vtzcXNGhVsa3fs7rQ79+/fCawObNm62trd3d3ceOHaunp9f4dykvKweAYkFxPcfbKtnWZ1iFsGJxymJVhurejnsbLo4iZPq73draWpRiXJ/x//vf/wDg/PnztH9p164dAFy+fJlGoz15UtfdCqEB5OXlBQUF0el0X1/f3r17A8CzZ88af9quXbsCQGxsrLgRP8WHWhlWVlZ0Oj0kJKT6ZmnjUVRUHDt2bGVlpei+oZF0UehCA9qbijdNcjYRO7N2xnPj9eT1tmduX5K6BD9KhCXpvPQlqUtCSkOa9u2aFpl2ox06dBg+fHh8fPyVK1fqM97Z2XnTl+AOkebm5ps2bap7VYjQAHx8fCorK8eOHWtsbDxjxgwajXb48OGysrJGnhY3WcKVtkXgKCicCNDKUFNTGzNmTFpa2sWLFxt/ttTUVAlLaWkpAMjL10EUD1kAABgbSURBVJwv9K1oMDT6KvWN58aHlYU1yQkxSVVJGgyNLF5WQEGA6MEVcosERQEFAYlViU34Xk0PldFW9eDVq1fy8vJKSkrnzp3DsXIIobS0tF9++aWyshIhtGXLFmdn59qqnX/+/BlI3GjzIBAIjI2NAeD27dvY4ubmBgB2dnbv37/HFpx7c+HCBYTQx48fnZ2dvb29Jc5TPW6Ux+NZWFgwGAwcRykUCr29vWk0mr29vTQujApev37NYrHYbLaPj4/ow5yRkbFjxw6cJLZt2zZnZ2eJCNzqcaNRUVFsNnvjxo2ikQ8fPmSz2erq6qKA/MYTUBAAEWAZY8mp+qJURQ4vhyfkoS/jRl+UvnhW8kziDHXHjYpoKXGjsu5GEUI3btzA9+aqqqo9e/bELZVoNBqOxMYRNrVF4xM32nzg5nEmJiaizCIulzt//ny8QWRoaNijRw+cjNutWzeE0PPnzwFg3rx5ePD27ds1NDQ0NDSYTCYAqKura2hoiFoKv3nzBjd9MjIywn9xY2Pj1l26+N69e3h9WUVFpWfPnh07/hNOdP/+fYSQjY0NAGBXmJKSgn91OBSUzWbjp0VFRTweb9GiRQCgrKxsZWXVuXNn/POtW7eaVu2qtFUQAWrRaktTlh7KObQzc+d3Cd8xI5m3Cm+hL91op7edmJFM/HNcRdym9E2b0jeN/DQSImBywuRN6Zt2Zu6s7V1aihttGaVJCgsLAwICXr16VVhYqKOjY2RkNGbMGBwWGhMTU1ZW1rdv3xoj7SsqKrZv396jRw+cp0hoQsaNG3fnzp1ff/31+++/F7dHRUVdu3YtMTERAAwMDLp16zZx4kQVFZXS0tL3799raWnhOey9e/eCg4MlztmrV69Zs2bhn/Pz88+fP//69WsajWZra+vi4tKEAeSySXFxcUBAwMuXLwsKCrS1tQ0NDUePHt29e3f493Pep08fOp1eUFDg5eVV/eU//fQTDlAJCwu7evVqQkKCgoJC9+7dXVxc6o6mahiXCy8fzjkcWhbKR3wmjWnANHBUcdyovdGEZeLwySGbl/2u6zsA2JSxqVhQjAOeHpQ8mJE4Q+I8OvI677u+r/EtzGLMtOW1n5k1wYJ7s9Iy3Ggj4fF4qamp+JuZ0CSkpKR07txZTk4uLS2txsJLzQeHw2l8MAChqUCACgWFagw1umxvtDQrrf/K4+LijIyMxo8f3xa+MKTG0aNHBQKBk5OTNH0on8+fNWtWly5dqu+iEKiCBjQNhkZb9qHQFtyoqakpnU6Pi4vD6/GExsPlcs+ePQsAy5bVN8mvScCJ5OXl5Zs2bZLm+1JLbm5ucwRCEZqQ1u9G5eTkaiyrQWgwf/75Z05OjpWVlXjSt3TYs2ePkpKSv7//06dPpfzWlIAQmjx5sq2tbXx8PNVaCLXS+t0oALi5ucnLy1+/fj09PZ1qLa0B/IWE+61KmY4dO+KmJmvWrMHpoa2bixcvhoaGpqen15hCRpAR2oQb1dHRmTJlCp/P9/HxoVpLiyc6Ovrvv/9WV1d3cnKiRMCmTZsMDQ2joqJw19JWTEVFxZYtWwBg9+7duDA+QTZpE24U/l3FO3HiBI/Ho1pLy+bIkSMAsGDBAtzlQvqw2WxPT08A2Lp1K64F1VrZvXt3amqqjY3NvHnzqNZCqBNqw1alCY6/u3TpEtVCWjAFBQVKSko0Gq3uRlhSwN7eHgC+//57amU0H6mpqbj36tOnT6nWQvgKbWU2CgC4OejRo0epFtKCOXv2bFlZmaOjo7m5ObVKDhw4QKfTDx48+PHjR2qVNBPr168vLy+fPXv24MGDqdZC+BpU+3HpUVxcjBeY3r59S7WWFolQKMTes+5mglIDZz1OmDCBaiFNT0hICI1GY7PZouaABFmmDc1GVVRUZs+eDQDe3t5Ua2mRPHz48MOHD3p6etVbY1HCrl271NTUbt68ee/ePaq1NCVCoXDNmjUIoR9++IGUJWsRtCE3CgArVqyg0Wjnzp0rLq5v0VmCCLwesnTp0qYqudZIOnTosHXrVgBYt25da9o59PHxCQ8PNzAwWL9+PdVaCPWD6umwtMErTceOHaNaSAsjIyNDXl5eXl4+IyODai3/weVycau7gwcPUq2laSgqKtLR0QGAP//8k2othPrStmaj8G/k0++//061kBbG8ePHeTzed999J1NlQZhM5p49ewDgp59+ysvLo1pOE7B9+/asrKxBgwbhhtKElgHVflzacLlcXMgS9yYh1IeqqirsPXGNV1kDl8Rfvnw51UIay6dPn1gsFp1Of/XqFdVaCN9Am5uNMpnMxYsXA0mx/xYCAwM5HE7Xrl1xtKas8dtvv8nLy3t7e79508QNgqTM6tWruVzu4sWLcbdqQouBaj9OAampqQwGg8lkZmVlUa2lZTBkyBAAOHLkCNVCamXVqlUAMGzYMKqFNJygoCAAUFVVzczMpFoL4dtoc7NRADAwMBg/fnxVVdWpU6eo1tICiI2Nffr0qbKysouLC9VaauWnn37S1NR8/Pjx1atXqdbSEPh8/po1awDAw8MDbzERWhBt0Y3CvxtN3t7ebaFKUCPBk9C5c+fKcnUMDQ2N7du3A8CGDRsqKyuplvPNHD58OCYmxtTU1N3dnWothG+H6ukwNQiFQhwoc/36daq1yDQlJSXYe75584ZqLV+Bz+f37NkTAHbt2kW1lm8jPz8fN21s8sZzBOnQRmejNBoNdwMmG01188cffxQXF9vb2/fo0YNqLV+BwWDs378fAHbt2sXhcKiW8w1s2bLl8+fPjo6O48aNo1oLoUFQ7ccpo6CgAFfQ+fDhA9VaZBc8v/Pz86NaSH2ZMmUKiLVxln3evXsnJycnJyf37t07qrUQGkgbnY0CgLq6+syZMxFCJ0+epFqLjPLs2bM3b95oaWlh39Qi2Lt3r4KCwvnz51++fEm1lnqxYsUKPp+/cuXKbt26Ua2F0EDarhsFgJUrVwLAqVOnysvLqdYii+AVjyVLlrBYLKq11BcTExNc12P16tVI5nvBXr58OTg4uF27dtu2baNaC6HhtGk3am1t3bdv34KCgkuXLlGtRebIzc0NDAxkMBg4W6EFsWXLFl1dXQUFhYKCAqq1fIW8vDwlJaWdO3fiLSZCC4Um+9/YzcqZM2cWLlxoa2v74sULqrXIFjt37ty2bdvkyZNbYiRmWlqagYFB3WOsra1jY2ObVYatre1XO5hyOBxtbW0Gg9GsSgjNSlt3oxUVFQYGBvn5+a9evSIZeCIEAoGpqWlycvL9+/dHjhxJtRxJXrx4ce/evZSUFACwtLScNm2asbFxjSOLioouXLgQFRWFEOrdu7eLi4uKigo+1LVrVym4UfGv56SkpGvXrn348KG8vLxTp04TJkyorUM1QujmzZuhoaHZ2dn6+vojR46UzTRcwj9QuL0lI+zevfvHH38kiaHi8Hi8M2fOTJs2TSgUUq3lC4qLiydOnIg/utra2jjhh8ViHT58uPrg6Ojojh07AkD79u0VFRUBoFOnTvHx8fgol8utaGa4XC5+L6FQuG3bNlynVU1NrWPHjnJycgAwf/58Ho8nIbukpASXczQwMLCzs8OtlZ2cnAQCQbP+bgkNhrhRQotBKBTiYk4TJkxISEjAxvfv39vZ2Q0aNKiyslJ8MI/HMzU1pdPpZ86cEQqFPB7v7NmzdDp9wIAB0le+c+dOALCwsHj27Bn+ZsrNzV2yZIment7Hjx8lBu/YsQMAli9fzufzEUIVFRW43UALCjtra7RmN5qRkbF+/foePXqoqqpqaGj079//6NGj+KMpgVAoPHPmTP/+/TU0NPT09EaMGHHv3j3pC6YQLpd75MgRe3v7Dh06KCkpmZmZrV27Nicnp8bBfD7/0qVLdnZ2tra2BQUFUhN5+/ZtAOjXr5/EDK60tFTChyKEcGeRWbNmiRtHjRoFAOHh4c2uVYy8vDwFBQU1NbX09HSJQ58/f64+HoeXRUZGiiyBgYEAsH79+uYVSmgorXanPjQ0tHv37vv27WMymXPmzHFycioqKlq+fPmkSZP4fL7E4AULFuDG64sWLZo4cWJERMTo0aPbTg/RwsLCIUOGuLu7x8fHjx492tXV1cjI6ODBg7169YqPjxcfWVFRceDAAVNT0xkzZoSEhLx8+VKaCex//vknAGzcuBHfEYtQUlKqHpKVkJAAABKLjxMmTACAr277NC3Xr1+vrKycO3euvr6+xCENDY3q49u3bw8A0dHRIgtOyiJ9mWQXqv14s1BQUKCrq8tgMM6dOycyCgQCDw+PvXv3SgyOjY1VVlbev3+/yJKYmMhisXR1daUkl2pmzZoFAM7OzuXl5SJjcHDwzJkzi4qKxEeWlJR069Zt5cqVDx48wOXZpVnVrXv37gCQlJRU24Dz58+rqanhhiK+vr5QrZE9ns+uXLkSSXFtdMWKFQBw5syZ2mQnJSWpqalNmzYNPw0NDZWTk5OXl9+6dWtGRkZqamrHjh0tLS0l/hYE2aF1utF9+/bBt5RDz8nJkdhL6dy5M4PBqHEFoJURHx9Pp9ONjIwqKiq+6YW4zao03SgODKqqqqptgK+vr5aW1u+//44Qys3NZbPZ6urqolv46OjoAQMGAMCCBQsQQpaWls09R7G1tUUIzZw5EwAePXpUm+zk5GQtLS0nJyeRxc/Pj0ajAQCDwVBTU7OwsMjLy2uqXyOhyZH72iehRYLXxebPn1/bgPT09I8fP5qZmeGdXC0tLfGjCQkJycnJNjY2bSGa76+//hIKhbNmzVJQUKhxQEVFRVhYmKamJs6vpxY8TavtqLOzs7OzM/5ZU1PT29t70aJFtra2RkZGVVVV6enpOKYNXymLxWru7Czx87PZ7NqGGRoa5uTkiJ7Gx8dv2bLF0NDwf//73+PHjwMDA+Pi4qZMmXLlyhWJDypBRmidbhTHA3bp0qW2AVevXl21atWhQ4dwPigmNTU1JSXlzZs3e/fubdeu3fHjx6WhlWrw78rU1LS2ARkZGQ4ODmPGjLlz544UddWApqZmdnZ2QUFBjUuK1ZkzZ461tbWvr29qaqq+vv6UKVPev38fHh6Om3FFRUU1s95/0NTUBIDMzMz6DEYIzZw5k8PhxMTEdO7ceeHChUeOHFm7du3p06ddXV2vXbvWzGIJDaF1utGioiJlZWV1dfXaBvTo0WPp0qUSxd88PT1xFjmTyTx69Ki1tXWzC5UBioqKAADPymtETU1t6dKleF2SWqysrGJiYp4/f44DgOpD9+7dPT09RU9x1q+NjU2z6KuFXr16AcDz58/rU+ElNTU1MjJy5MiRnTt3xhZVVdWTJ08+ePDg7t27AoGgLdwhtTyoXlVoFrS0tNhs9re+Kjc3Ny4uLigoaOrUqQCwbt265tAma+C6q3fv3v3WF0p/bRSnpQ4fPrxhSQElJSU6OjpqamriO2lSICsri8lkampq1md98+PHjwAwePBgCXunTp2YTGb1WH2CLNA6A5709fUrKirqeRslQlNT09zc3NHR8dKlS7179z5w4EBqaiqfz587d+6dO3eEQmEzqaWE9+/fP3jwAABwFE5iYiLVir7OhAkTBg0a9OjRo4ULFxYWFmIjQujJkyfBwcEAEBYWtmrVqocPH+JDhYWFotpdubm5s2bNysrK2rZtWx3LlM2Btrb2unXr8vLyxo4d++nTJ5H948ePfn5+AJCfn79q1Sq8iNS5c2djY+PQ0NBbt26JRp4+fTo1NdXR0VEi0osgK1Dtx5uF5cuXA8CpU6caeYZ79+4FBATgX5S+vr6Hh0dtEektBT6ff+PGDUdHRxqNZmhoyOfzcUPKCRMmfOuppD8bRQhlZWXh3XYWi2VrayvKlRwzZgz6t3Ssl5eXSKGKikq/fv369OnDZDJpNNqKFSsoSank8Xiurq4AQKPRunXrNmzYMBMTExqNpq6uzuVycXDu6NGj8eCnT5/ixi2DBg1ycXHBoa9GRkZ1RHoRqKV1liaJiIjo27evoaFhREREfUqQ+fv7T5s2Tfyr3sHB4dGjRxERESYmJufOnTt06BCer7FYrIkTJ7q5uTk6OjbjBTQD6enpJ06c8PHxwZN0FRUVFxcXLy8vRUXFLl26JCcn//XXX990Uc7OzhcvXszMzJRyJ0uhUHjjxo27d++mpqbKy8sbGhr26tVr+vTpOE3o3bt35ubmuFJJTEzMuXPncLalmZnZ9OnTbW1tpSlVgpcvXwYEBMTFxfH5fH19fQsLi2nTpuFQsydPnmhpafXu3RuP5HA4Z8+efffuHYfD0dDQGDp06KJFi5SVlSkUT6gLqv14c7F69WoA6NKly8WLF9PS0jIzM4ODg93d3XFs9unTp/X09E6fPo0QwveDDg4OoaGhRUVFSUlJGzZsAAAbGxvRGpxAIAgKCpo+fbrI1VpaWh44cKC0tJTKi6wfz549mz59uihOyNzc3NPTUzwNMSgoSF5eXlFRcceOHTExMXl5eVFRUXv27MHdAZKSkvT09ObMmYMHZ2RkhIeHh4eH4/T2+/fvh4eHizLcCYQ2SKt1o0KhcOfOnRJf4HJycjit5dixY4qKiseOHcODT5w4IdE9uF+/fikpKdVPm56e7uHhIQrfU1VVdXNze/v2rVSvrX4UFRV5e3uLdtgZDMb48eODgoJq3J8JCgqqHh9mY2NTWlqakJCgqKg4ZcoUPLLGOu0zZsyQ7sURCDJE67ypF1FaWvr8+fOUlBQFBYVOnTr16NEDJyzXOPLx48fJyclsNrtXr159+vTBaSQ1UlVVdf369RMnTuBdGgAYNGjQ6tWrJ0+eXEdwuNSIjIz09va+cOFCWVkZAOjq6s6dO9fd3b3uSsZCoTA8PDw2NraqqkpPT8/CwkIUcyNOQkJCUlKShFFbW1v2W4cSCM1EK3ejzU1sbOzx48dPnz5dWloK/zqs5cuXU1JFgsvl3rhxo7pznzJlCtnhJRCaD+JGm4Di4mJ/f//Dhw+/e/cOABgMxpgxY1avXu3g4FDHlLYJycjIOHny5NGjR3NzcwFAVVXVyclp1apVpNkkgSAFiBttSiIiIg4ePOjv78/j8QDAzMxs4cKFrq6uzdSwTCgUPnr06MSJE4GBgQKBAABsbGyWLFni7OyspKTUHO9IIBCqQ9xo04PnhidPnsRlIpWVlZ2dnX/88Uc9PT3RmPLy8ri4uMLCQm1tbUtLSzq91jyIqqqqT58+ff782cTERFSwEiG0Z8+e48eP42VKBQWFGTNmLF++vLbePgQCoRmhcn+rVSMe6M5kMrOzs7E9MzNzzpw54uWUdHV1T548WeNJDh8+LB6V2bVr14iICHxo2LBhAGBiYuLp6dnSkwIIhBYNmY02O7GxsS9fvpw3bx4A4CSc5OTk0aNHz549W1dXNzEx8dixY9HR0f7+/rgwpQgvL6/NmzcPGTJkxYoVWlpaQUFBXl5eGhoacXFx7dq1CwkJKSkpGTVqVB0zWQKBIA2o9uNtixkzZkC1kuxcLtfb27t6ieiNGzeOGTNGvBoFrqN+5MgRaWglEAj1g8xGpUd2dnbHjh21tbUTExOZTGYDzhAYGDh16tQ1a9bs37+/yeURCISGQe4HpUdoaCifzx85cmQdPvSHH37Yvn17bUcrKioAQEVFpVn0EQiEBkHcqPRITk4GAENDwzrGHD582Nvbu7aj9+/fBwB7e/umlkYgEBoOSW6RHrgXMS7sVhuXL1+uLeMoKirK39/f2trawcGhWfQRCIQGQdyo9MA34zhttDZw2aTqZGdnT58+ncVinTlzRjqZUQQCoZ6Qm3rpgUsovX379ltfmJGRMWzYMA6Hc/36ddzYh0AgyA5kp156FBcX6+rqslislJSU+m8TPXnyxMnJqby8/NKlS6NGjWpWhQQCoQGQ2aj0UFVVXbx4cUFBgZubG5/PFz+E23MCwNu3b2NiYkT2EydOjBgxQkVFJSwsjPhQAkE2IbNRqVJaWurg4PDy5UtLS0tnZ+dOnTplZWWFhobev38/JSVFS0tLWVlZVVWVw+EghGbOnBkQEKCgoLBu3TrxqtLa2trz58+n7iIIBMIXEDcqbcrLy/fs2XP27Fkc/0Sj0QwMDCZOnOjh4aGpqTl37lxFRcXjx4/zeLzawkutra0jIyOlKppAINQOcaOUUVBQgDunNyyjiUAgyAjEjRIIBEKjIFtMBAKB0CiIGyUQCIRGQdwogUAgNAriRgkEAqFREDdKIBAIjYK4UQKBQGgUxI0SCARCo/g/0IyfT95wcR0AAAETelRYdHJka2l0UEtMIHJka2l0IDIwMjIuMDkuNAAAeJxNj71OwzAQx8/nxHFiuxUSoh4ZiMTGyoKSCVV9hS4RUxZ4gS5siAdpBySEmNiSF0FIvADsTNwlqZ2TTv75f9+/3esXkDkYTZAvJ38UChrWpNL8InoRCAPFaBIoDaT0QEL4TK1bIimF1xHziEVEE9Gq9SVj4t2JPudWfmFoT4GAkmRIUkgVqAx03mJeNIVp0djGuhbdAhyNzzSeooZMpYlEOsFZU+Srbjp3sOXV519FJ1/w5+nuvSLvma+f3/YU6ye9P+qs7bYv5UyvmHfbTTn2Cn0Gth+H+vv2MNSuHu5nfFPT3HKmz/Mn7jhnP/JPfVz67B/KQ0gsE1q9jgAAAVl6VFh0TU9MIHJka2l0IDIwMjIuMDkuNAAAeJx9k01uhDAMhfecwhcAxfmBeDnAaDSqBqSW9g6Vuuz9VZtp4iAGAkFJ9OnF9jMVyHgf375/IQ87VhWAOXmJCL6cMaZ6gCygv97uEwzLpU8nw/w5LR+AQV4jz5a9LPMjnSAMgA0Fit6Dacw6ikXi7Mq5GCIS1HmVQUygY7DOx/Ve2ybSC3mm6RIZnuRhlD6B7UZyrxgS2DGose+U28RFmOFEr0scsV5orFllOAJLZFxUMOZ6sz3gmu6QpEyKNfY4aczeoC0195KYzUFXxvmCzOagX8n/2+sX12d3uMmGH+gafxhp9uc6jZvue/ZjP0+j9qM8VtuON9oIjrdO7eYUeAZtE+TZajMgz04tR55RnfXsmhaRRG1Tq1Ufy5rIbVzsInUrJ67IUDZSPAxFJb18qKiXpBXLopQlkH36nXld/QGhl8aFXcYJ7wAAAPl6VFh0U01JTEVTIHJka2l0IDIwMjIuMDkuNAAAeJxtkruOwjAQRX9lS1ZyLI8fsT0RBUpDA0sfuQBERxSEoOPjmcgkRDtucpJz71iKxi2suz/0qWu3GusMl2EzTIbOgASyO28xJBgR0/gElQEZVD0jGGq0VwSbfl4rkNHFYIUSDUgTXNBiYqPElDYVC6vF5JxWi9m50PDUSa3GiIyOUVHHSP/P6O/5Uzhn0/xCzGUv7ef1VxwfQ3+4DzcE2Q/XDX3tjrf9sz9d7tKj5rJGw6VDy6VFx6XBmkuNnkvAwGXAyGWkLRbmFW21oIG2XNB0Pwo/BoZuQUHb1xvYRLVkm2x6jQAAAABJRU5ErkJggg==\n"
          },
          "metadata": {},
          "execution_count": 34
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "m = Chem.MolFromSmiles('CCCC')\n",
        "m = Chem.EditableMol(m)\n",
        "m.AddAtom(Chem.Atom(8))\n",
        "m.AddBond(0,4,order=Chem.rdchem.BondType.DOUBLE)\n",
        "m.GetMol()\n"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "eB7yfZu7_e38",
        "outputId": "e889c0c9-cebc-41ac-bb8f-212367a53a25"
      },
      "execution_count": 51,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c742dd60>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAASY0lEQVR4nO3da1BU5xkH8Ie74A0XFFDwgmAgKJdoMIkaL5CASj9kLJP6YTP1Q+lMx2Fsmglpp8lqZzJDTW23maYt5kNnU5uqM2EyGzVWUIgagyKQaAQFNALiBRQwIArIvv3wrrvLgsDue86ec3b/v/HDkeye80bO/vdc3uc5fowxAgAAd/krPQAAAG1DjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYwCAAhBjAIACEGMAgAIQYx6wpEjR9LT09PT0z/99FOlxwJeoqmpKTMzMyMj45133lF6LL7OjzGm9Bi8WUtLS35+fnV1te0nCQkJpaWly5YtU3BUoGmPHj3atm3bgQMHbB/e2NjYw4cPp6amKjswn4WjUbkMDQ1t27YtPj6eZ2hCQsKKFSuIqLm5OS0t7dVXX+3t7VV6jKA9u3fvDg8P379/P2NMp9Pl5ub6+/vfuHEjIyPjjTfe6OjoUHqAPomBDCoqKhITE/m/sG2nZ4wdO3YsJiaG/zwmJubgwYPKjhM0pKmp6bXXXuM7T1BQ0Ntvv81/fuPGjaKiouDgYCKaNWtWcXHxwMCAskP1NYhRibW1ten1ej8/PyKaPn369u3bLRaL02uKi4t1Oh3/PGRlZV24cEGRoYJW9PX1GQyGkJAQIgoJCcnOzv7xxx+dXnP58uVNmzbxnWrJkiWHDh1SZKi+CTEqmQcPHhgMhilTphBRWFiYwWDo7+9/2ouHhoZKSkoiIyOJKDAwsKCgoLOz05OjBU2wWCwmkyk6OpqI/Pz89Hr9rVu3xnl9WVnZs88+y8M0Ozv70qVLHhuqL0OMSsBisRw8eHDBggV8X8/Pz79+/fpk3njv3r3CwsLAwEAi0ul0RqNxaGhI7tGCVpw9e/aFF17gmZiZmXnmzJnJvGtwcNBoNM6cOZOf+xcWFvb09Mg9VB+HGBV1/vz5VatW8X19+fLlp0+fdnUNDQ0Nubm5fA1JSUlffvmlHOMEDblx44bt0tC8efNMJtPoS0Pju3v3bmFhYUBAABFFREQYjcbHjx/LNFpAjLrv5s2bBQUF/v7+/H5RSUnJ8PCw22szm83x8fE8TPPy8pqbmyUcKmhFf39/cXHxtGnTiCg0NLSoqKi3t9fttdXU1KxZs4bvVBkZGV999ZWEQwUbxKg7BgYGjEbjjBkziCg4OLiwsPD+/fvSrpafjkmyWtAKs9m8cOFC21fpDz/8oObVgg1i1GVOh41Xr16Vdv3SHuSCJjgdNp48eVLa9Y8+yB19rx/chhh1gScvYjpdcj116pR82wIFefIipvglVxgTYnRSFLml7vYEANAEpW6pnz179sUXX+Tf0M8///wkJwDAOBCjE1B8giefjhoaGjqZ6aigFcpO8HR1OiqMDzE6nvLy8qVLl9rKjS5evKjUSByLo2JjY3E6pl2O5UbPPPPM4cOHlRqJY3HU1KlTDQbDw4cPlRqMpiFGx9bY2Jifn8/39YSEBJUUv1dWVqalpfFRrV27tq6uTukRgQu6urpUWPze1NSkwl1dWxCjznp7e52+oh89eqT0oOyGh4dNJlNUVBQR+fv76/X627dvKz0omAD/rc2ZM8f2W7tz547Sgxrh+PHjtuaNGzZsQJ8HlyBG7TSUUN3d3UVFRTzrw8PDi4uLVZX14OjEiRO2TqDr16//7rvvlB7R2PhtgNmzZ9tuA3R0dCg9KG1AjFpVVVU51i9/8803So9oYleuXMnLy+NjTkxMxOmY2rS2tur1ev4LiouLM5lMSo9oYl1dXbZJKbNmzUKfh8lAjI64e6PFyXRlZWUpKSm2e77ff/+90iMC690b3u5Li3dvGhoaNm7cqIb7YJrg0zH64MEDW2lHWFiYYP2ygvgMxPDwcD4DEW33FMRn+8bFxdlm+7a0tMi3uU8++WTv3r0y1bmZzebFixfbvqHr6+vl2IoX8N0Y9b5CY14jwOtheI0Amvp42Llz51566SXbzPavv/5a1s319PTwGc3ytR3h39Do8zA+X4zRmpqa1atX8339ueeek7x+WVn19fU5OTn8/y45Ofno0aNKj8gntLe32zohzJ0712OdEMxm86JFi2xHA9euXZNjK7zPA/+GjoyMxDe0E9+K0c7OTtvxmnfvDU6fLsn7p4AN7/oxffp0etLuy8NdPxwHIGvbEe8+/hDhKzHqg+cmvO2egh9vX+CZg8HJcGw7MnfuXJnulFosln379sXGxhLRMxERQz//OWtrk3wrmuMTMVpWVpacnOybV8qVOtn0enV1dS+//DLfqdLT0ysrK5UeEWOMnTt3zgNtR/r6+t59991LGzYwIhYWxnbtYr7d58HLY1Q99cvKcrz1sWLFCrlvfXg3x9Z2KryV59R2JD8/v7W1VZYttbUxvZ75+TEiFhvLTCamqZmCEvLaGHWaRayS+mUF8Yk48+fP98xEHK80urVdd3e30oMam+cmrlZVsZUrGREjYpmZTAt1K5LzwhhFTds4tD4tXEFaLHNobm62tR2RsYxqeJiZTCwqihExf3+m1zO1VlHLxNti1KnDgmrrl5WlxSJFBTkW3S5ZsuSLL75QekSuOX78uCeK+nt7mcHAQkIYEZs6lRkMzGf6PHhPjDr2+1q8eDEKzCfk1DLj22+/VXpEquM1LWB42x1+isbb7sh1itbYyPLzref4iYnMNz6G3hCjOFF1m/obuCll9L+Matt9TZ7nGp6Wl7OlS61hmpXFlOt37hnajlE8C0ESjp8ufszl47fjKioqHNtje9lxuoemrwwNsZISFhnJiFhgICsoYN7b50HDMeqZKXK+4/Lly5s3b7ZdATx06JDSI1IAv2rsCw9rcZpMLdfDoO7dY4WFLCCAETGdjhmNzBvb7mkyRvGcWPko+6g1BfFHB/JLQ/zRgV5/achzjyatr2e5udZz/KQkJueTyRWhsRj1WPmwL1Pqwb9KGT2j1qceZO1YTRARESFjNYHZzOLjrWGal8e8qM+DlmLUqbWdgvXLvsBzny5FVVdXr1q1iu9Uy5cvP336tNIjUkZNTc2aNWtkbzsyMMCMRjZjBiNiwcGssJB5RWsLbcRobW2trX5ZvtaKMFptba3t0+Vl//K8+RvvNhATE4NuA8xjTXhv3mQFBczfnxGxmBhWUsI0/i+v9hj1kWMilfOy8wDe+4q3++K9r7y+3dfkOT4Sgl83k+uRENXVbNUq6zn+8uVMy+cB6o1RX7tCp3L8qrTjp0ujV6XNZnN8fLztKwGdWMfU0tLy+uuv87u4zZs3s/37ZWk7YrGwgwfZ/PmMiPn5sfx8ps2r0iqNUaf7xZqoX/YFmp4jUV9fn5uby3eqpKSkL73ufrHkTp48+actW6wHjKtXs5oaWTbz4AEzGNiUKda2ewYD09ocCdXF6JUrVzB7UeU0N2N39FOq8NDgyeJtR6KjrQeMej2TqcKltVW7bfdUFKO8fhm1NJqglfox3u6LP/eNt/vCM1Pd0dfnobYjFRUsLc16/Lt2LdNI/ZgqYhSV3Rql8m4G5eXlS5cu5UfNWVlZF729slt2TU32tiMJCXK1HeHHv3Pm2NvuqT4NlI9Rxz5D69at87L6ZV/g2FsrISFBDb21GhsbbUNKTExUw5C8R3k5W7bM3nbkwgVZttLdzYqKrMe/4eGsuFjNbfeUjFF0vfQmKun02tvbazAYeGu7adOmGQwGjba2U7XRbUdkart35QrLy7NG9pIlTK2dXpWJUceTQR+pX/YFyj53gF8aioqKIi9qbadqvO1IYKDsbUfKylhKijVMs7OZ+ubteDpG8UQgr+f0FCzP3BavqqpauXIlPxbOzMz8xiefCKSMhoYRbUeOHJFlK4ODzGhkM2cyIhYUxAoLmZqeguXRGK2ursbzKX1EQ0PDxo0b+e9a1meytrW1+UhrO1VzajvS3CzLVu7cYb/4hbWKVE3Xuz0Uo3haum8ym82LFy+2lVHU19dLuHLHssWwsDAZyxZhMhzbjvADRplKbGtr2Y4dqppVKnuM8vpl3tqO1y9rtIgQ3MOLenkBOy/qlaSA3Ww2L1iwQPYmGuAqr2s7MhnyxqjZbF60aBHql4G3U+J1RJGRkSItZmpqalavXm1r6Xbq1ClphwoSOH9+RNsRz/+O+vpYSwurqWFXrzL5j9vkitH6+vqcnBy+rycnJx89elSmDYGGCCaghFkMsuNtRxYs8FzbkZ4e9vHHLD+f6XTWBLf9mTqV5eWxDz+UqZJV+hgdXb+MfR0cuXE+LtOVAZAdbzsSGmpvO9LfL/1WHj5kH3zAIiKc03P0n6lT2e9+x6TuFSdljA4ODtrql4OCglC/DE8z+u7QOF0QP//8c/nuU4EntLXJ2Hbk3j37BYRJ/omPZ42Nkg1AwhgtKytLSUmx7euoX4YJOc5VCggIKCgocJq/cfjwYT6ZX+5ZU+AJlZUsPd3edqSuToJ13r7NkpJGRGRqKvvjH1ldHevoYBYLu3uXXbzIjEa2cuWIl0VHM+ke1yhBjDY2Nubl5aF+GdxTVVXFu9IQ0fTp03lNcHNzs63Tgr+//549e9DazhvwtiNRUfa2IyKVZhYL27TJnoyhoezjj586McBiYf/9r3UCvy1wJSoUForR4eHhN998MygoiLe227Nnz+DgoCTDAp8yNDSk1+v5tGIimj17Nj9EJaLU1NRmmeZyg1Kkajvyj3+MuOhZUTHxW86dY+Hh9ncVFbmz3VFEj0a3bt2K+mWQRHt7u60bNL8/WVpaqvSgQDaNjfa2I4mJLlclPX7MFi60B+KHH072jfv22d8VFsbu3XN14KP5McZIQHt7e2dnZ3p6ushKAGyOHTtWWloaHR29c+dOpccC8isvpx076NIlIqLsbDIa6cktlgmUltKWLdbllBS6cIGenM1MbP16qqy0Ln/wAb31lisjHoNojAIACBkaor//nXbupJ4eCgqibdvo/fcpMnKCd/30p/TZZ9blP/+Zfv1rF7Z44AD97GfW5YwMqq11fdAjIEYBQAW6umjXLvroIxoeJp2O3nuPtm+ngICnvn7ePLp507rc2kpxcS5sa2CAwsPp0SMiosBA6u6madMEhk6TPgwGAJCPTkd//StdvEg5OdTVRTt20LJl9L//jf3i69ftGTpnjmsZSkQhIfSkxTg9fkzV1e4O2goxCgCqkZxMR4+S2UyLFlFDA+3dO/bLrl2zL9sC0SVpafblq1fdWYODQMH3AwBI7Cc/oexs2rOHtm4d+wVdXfblCa+ijsnxXd3d7qzBAWIUANQnNJR+//un/teeHvvyjBnurH/mTPuycIzipB4AtObxY/tyoFvHgsHB9uXBQcHhIEYBQGscjyX7+txZw/379mWdTnA4iFEA0JrwcPuyYyBOnuO7HNfmFsQoAGjNk142RO7eZ29uti9HRQkOB9PvAUBrBgZo5kwaGCAiCgig3l4KDXVtDY6z91taaP58keHgaBQAtCYkhDIyrMvDw3TmjGtvb2qyZ2hMjGCGEmIUADQpK8u+/O9/u/Zek8m+nJ0tPhac1AOABl2/TgkJNDxMRBQWRvX19OQBXxPo6qLkZOrosP71zBlyaM/oHhyNAoAGLVxImzdbl/v76Ve/okkeEb71lj1DV6wQz1BCjAKAVv3lL/bOTEeOkF4/Ylr+mHbupH/9y7ocGEh/+5skA0GMAoA2xcfT7t32v/7nP7RuHdXUjP3ipibKy6Ndu+w/+e1vaeVKSQaCa6MAoGUGA/3hDyN+kppKOTkUF0dz5tDdu3TrFpWXU3U1WSz21/zyl/TRR+P1M3UFYhQANO6f/6Tf/Ib6+yf14qAgeu+98fqeuA4xCgDa195OO3eSyURDQ099TUAAbdlC779PCQnSbhwxCgDeoreXTpyg8nJqaaHOTrp7l3Q6mj2b5s2jDRvolVfEu5CMCTEKACAEd+oBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIQgRgEAhCBGAQCEIEYBAIT8HxMWtkBj/gsrAAAAfnpUWHRyZGtpdFBLTCByZGtpdCAyMDIyLjA5LjQAAHice79v7T0GIOBlgABGIGYFYhYgbmBkY0gAiTFDaCYmVJqRmYNBA6SJiRuoj5GJgYkZqJGDSYRB3A1qFBiwPnRbZp+W9swOJjBr5sx9IBoovh8m/tBN7QBCXM0BxhYDAM6KGgqrMT5qAAAAwnpUWHRNT0wgcmRraXQgMjAyMi4wOS40AAB4nI1Qyw6DIBC88xXzA5KVRytHFdM0jZC0tv/Qu/+fLjGIHmrchWR2mZ3sIJDi6R/fGWsoLwRAB8c5h48mIjEiAXTD7R7QT22XO318h+kFC8MTnHtmO8Uxd2r0qKVyjnQDksayLg9IogVknmLe2q1IavpD1EysziiaRFTSuoaul0NJi4hTxCH4nbnFbheDL3ZTquKKC+iye83XlA3Tu92qb7VSnT+dsfgBQnVXqlgXtr4AAABQelRYdFNNSUxFUyByZGtpdCAyMDIyLjA5LjQAAHicc3Z2drb1V6jR0DXSM7W0MLDQ0TXQM9ax1jXUM7K0NDDRMdAzMdWxNoAKo4qiaNGsAQBlEg+iUMCE2gAAAABJRU5ErkJggg==\n"
          },
          "metadata": {},
          "execution_count": 51
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "Replace_mol = Chem.ReplaceSubstructs(Chem.MolFromSmiles('CCCCCO'), \n",
        "                                 Chem.MolFromSmiles('CO'), \n",
        "                                 Chem.MolFromSmarts('[#9]'),\n",
        "                                 replaceAll=True)\n",
        "Replace_mol[0]"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 167
        },
        "id": "y_ZMRcnoOnRS",
        "outputId": "9c0dcab5-0ab3-40e3-e5eb-54f5349d0c64"
      },
      "execution_count": 54,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "<rdkit.Chem.rdchem.Mol at 0x7fd0c7431430>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAPtklEQVR4nO3da0xU577H8QcERETl5g0v0FpAq9YqCtZLjxfSi9Jd4xab2rC7XzT2mPaAcbcHk8aQ2pzEbtsTqmkjiUkPb3qx3U2L0NqAtS3d9bZAQC6KWoZRVOSiKHecmfPiOXsOharAs5g1a+b7eVccZv1D4TdrPev5/5ePw+EQAIDh8jW6AAAwN2IUAJQQowCghBgFACXEKAAoIUYBQAkxCgBKiFEAUEKMAoASYhQAlBCjAKCEGAUAJcQoACghRgFACTEKAEqIUQBQQowCgBJiFACUEKMAoIQYBQAlxCgAKCFGAUAJMQoASohRAFBCjAKAEmIUAJQQowCghBgFACXEKAAoIUYBQAkxCgBKiFEAUEKMAoASYhQAlBCjAKCEGAUAJcQoACghRl0hOzs7Li5u9uzZX3/9tdG1wEP09PSkpqbGxMQsWbKkrq7O6HK8mo/D4TC6Bk92+vTpzZs3WywW51fmzZv31VdfxcTEGFcUTG/fvn07d+7s7OyU/+nr6/vOO+9kZGSMGjXK2MK8E2ejI6WlpeWpp55KTEyUGZqQkPD4448LISoqKuLi4p5//vmOjg6DS4QJnThxIjo6Oj09vbOz08/Pb8OGDWFhYXa7/a233po/f/73339vdIFeyYERkJubO3nyZPkTnjp1alFRkfx6Xl5eRESE/HpUVFR+fr6xdcJE2tvbMzMzneebq1evbm5ulv/02WefPfzww/LrycnJly5dMrZUb0OM6qykpOTJJ5+Uv9ARERH79u0b+Jrdu3c7wzQpKamystL1dcJE7Hb7oUOHoqKihBA+Pj7R0dEnTpzo95ru7u6srKxx48YJIQICAtLS0lpbWw2p1gsRo7ppampKS0uTJwvh4eFZWVl3796914t7enqysrImTJgghPD3909LS7t165Yrq4VZaJq2fPly+aEbHx/vvLL5Q/X19Vu3bvX19RVCREZGZmdn22w2l5XqtYhRHQw7E4eUvPA2V69edWbi1KlTB5+Jp0+fXrZsmUzexYsX//LLLyNdqpcjRlUVFBTMnTvXeYVeUVEx1HcoLi5euXKlfIeFCxf+/PPPI1EnTEReoY8fP975wTzUK3S5DjBz5ky5DpCSklJXVzdC1YIYHb7z58+vX79exl9sbOzhw4dV3i03Nzc6Otp5l6C2tlanMmEyubm5fe8XXbx4cdhvJe9KBQYGCiGCgoIyMzM7Ozt1LBUSMTocN2/ezMjICAgIEEKEhITs2bOnu7tb/W07Ojr27NkTHBwshBgzZkxGRsadO3fU3xZmUV1d/cwzz8gAnT179nfffafL21qt1tTUVPm2M2bMyMnJ0eVt4USMDo3NZsvJyZk0aZIQwtfXNzU1taGhQd9DXLlyJTU11cfHRwgxbdq0nJwcu92u7yHgbpqbm9PS0vz8/IQQYWFhWVlZvb29+h7i2LFjCxYskGG6atWq0tJSfd/fmxGjQ+DKX8STJ08uXbpUHishIeHXX38duWPBQL29vdnZ2XIDnJ+f39atW2/cuDFCx3LBSYB3IkYHxZDLIrvdnpOTM2XKFHmXIDU19dq1ay44LlymsLBw3rx58vdq7dq15eXlLjjoCC1JeTNi9AEMX6Rva2vLzMwcPXq0EGLs2LHcJfAMNTU1KSkpMkAfeeSRQ4cOubiAfjdI8/LyXFyAJyFG78mttoxcuHDB2L866GXg52JXV5dRxahv14ODGL0X99zAfPTo0fnz58uq1qxZ45prQOhFLk3KYQtyadIdVmloqFNHjPbn5u108o7ExIkTXXBHAjrqd8/w+PHjRlf0OzTUqSBG/5+Jhju0tLQ498eEhoaOxP4Y6MVEO9j6DtZZuHDhTz/9ZHRF5kCM/p9+rSOmGDVWXV397LPPOndrf/vtt0ZXhN9pb2939lMEBQWZpZ8iNzf3oYcecv4t/Pbbb0ZX5O6IUUdVVdXTTz8tf2nmzJlz5MgRoysamtzc3FmzZunSOwgdmbq7VzbUySsz2VB3+/Zto4tyX14do7J1RK4HydYRk64HybsEKpMsoKO+s2YWLVpk3lkzfZcjIiMj3Xk5wlheGqM9PT39WkcaGxuNLkqVnKsmPxXkXDWTfiqYV2Njo+fdqDl16tQTTzwhPxWWLFlCQ91A3hijBQUFztaRpKSks2fPGl2RnjRNW7FihfNU6P5TfqGXgRcEnrRtiIa6+/OuGO3bOhITE+Opm9j7PnNCLsxZLBaji/JkBQUFc+bMcX4wV1VVGV3RiJCNA7Kjj4a6vrwlRu/cueNsHQkODja2dcQ1ZBvrmDFjzHWb2FzOnTu3bt06GaBxcXHe8IzCvg11jN2TPD9GB7aOXL9+3eiiXOfy5cvOuwTTp0/nLoFe+m3d9bYBH327SFNefPFiR4fRFRnJw2P0xx9/lE+HF0IkJiYOfJ6ilzhx4kRiYiI/B130bSSTH8ze2Ujm/Dls+Pvfl2jartraFm/tAfHYGOUsrB8vPyvXy9GjRx977DHnWIOysjKjKzJYU3PzuxZLQnFxvKatLi39vKHhrvf9oXlgjLImeB9euEasl75rgrNmzfLU+5PDY+nsTLtwIV7T4jVtY0XFLx60S2EwPCpG+96hlqPtuEP9h7xkx4JeuEM9SCdbWzdVVMgw3VZTc8lrfkqeE6Oapi1fvlxGQ3x8PPslH6jv/tm1a9d62P5ZXbBfcqh67fZPGhqePHMmXtMSi4v3Wq13zN+A8ECeEKOye0eOtqN7Z0g8sptLL3TvDNvN3t69VuuS4uJ4TVtTWvpJQ4PNoxdMzR2j9JLrwmNmC+jFRKPt3Fl1e/sr587Ja/yXqqpKPPcWhYljlMlG+uo36Uqvh6SbC5ONdPfTrVvPlZfLMN1+4UK9J+6uNWWMMmdz5Jhx7qpemLM5Qjpttv+5dm1lSUm8pi0rKdl35Uq7Oz1RQp3JYlRefjL1fUSZ6CkAemHquws09PTsqq1drGnxmvZsWVleU5PHLJSYJkZ5BpGLufkzqfTCM4hcrKKt7a/V1fIa/y/V1eVtbUZXpANzxGhhYSFPxDSEez4hVRc8EdMoNocjr6npqbKyeE1brGm7amubenqMLkqJu8coz2c3nGxqmDlzprOpoa6uzuiiVBUUFDz66KPO0XY8n931Omy27Pr6J4qL4zVtRUlJdn19t2m3Q7hvjMrWEdm2KFtHaFs0kGyxlZ08QUFB5u3kOX/+/Pr162WAxsbG5uXlGV2RV7N2dWVcuiSv8TecPVvQ0mJ0RcPhjjFK64jbslqtqampfWdNmmhD5c2bNzMyMgICAoQQISEh3jbazp2dvH37hcpKGab/fv78BbON3XO7GD158uTSpUvlH2pCQsLx48eNrgj9HTt2bMGCBfL/0apVq0pLS42u6AHkdKtJkyY5p1s1NDQYXRR+567d/o/GxrWlpfGallBc/F8Wi4nG7rlRjNI6YiImCqYffvjBOdrOFKHvzVp7e/darc6xe5+YZOyeW8SobB0JDg52to4w2s4UBl4mu9X69cAlCKMrwqDUdnb+R02NvMb/c0XFP91+B4WPw+EQhjp8+HBaWprFYhFCJCcn79+/Pzo62tiSMCQ1NTU7duzIz88XQsTGxr7//vvJycnGltTe3r5379533323q6srKCjozTff3Llzp7w/BrP4ubX1vy9fvtLdLYRYOWHC32bMmD569GC+8eumpp6hxNoYX9/nwsOHWaUQQggjY7SkpGT79u1FRUVCiEWLFmVlZa1cudKoYqCosLBw+/btlZWVQoikpKSsrCzns3pcyeFwfPnll2+88YbVavXx8dm0adN7770nd2vBdO46HF80Nh64erXdZvPz8dk0ceK2yMixo0bd/7tWlZa22WyDP0q4v//3/1r2GR5flW8etubm5vT09ISEhKKiItk6curUKTLU1JKSks6cOSM3tBcWFi5cuDA9Pb21tdWVNWiatmLFis2bN1utVtks4NzxCjPy8/F5cdKkf8yduzEiwu5wfHbjxp8rK79qarIbfQ3dj6vPRnt7ez/66KPMzMzW1lZ/f/9t27bt3r1bdpLAMzQ3N+/evfvDDz+02Wzh4eG7du16/fXXRz3oDELR1atX33777YMHD9rt9sjIyMzMzFdeeUV2ssIzVHd0vHf5cllbmxBiTlDQ32bMeDw4+A9f6Twb9ffxmTmIlZxQP78DsbFKxblyIbZf60hlZaUrjw5Xctmwj4GDVBht56nsDkdBS0tyebnsIs24dOnqH+38/bczZ+TtqT+56oEOLorRc+fOrVu3Tv5RxcXF5efnu+a4MNZIj57r9/5eNdbPa3XabNn19ctKSuI1bXlJSXZ9fdfvh+Z4YIy2tLQ498SEhobSOuJtRmgQcr8h00eOHFF/T5hIfXf3f/6ri/S58vKiPpuiPCpG5Q5tOdrOzXdoY6TV19c7eysiIyNVeit45Amcim/ffrGqKl7T8pubnV/0nBjt2zqyevXqsrKyEToQTETxIXF9H8Dn7+/PA/jgcDhsdntBS0vfz2RPiNG6ujpaR3Av/ebOpKSkWK3WwXxjQUGBcyNqUlISj4PGvZg7RuVoO9krIkfbmXSWGkbakH5VampqnG1RMTExzJzF/Zk1Rod9igFv1vfCJSAgYOvWrf1eYLFYli5dKpdBg4ODmTmLwXB9jOqw/f706dPp6enHjx+XC14ffPCBc/0LeKBjx45t2bLl+vXrQojQ0NCDBw9u3Lixt7f31VdflXeihBApKSn79++fPHmy0cXCBJzb78eOGrU+LOwBLw4JSRg/XvGISjHa0NCwY8eOTz/91OFwTJ8+fc+ePVu2bJF3Y4HB6+rqevnll7/44gv52xgdHX39+vWuri4hREhIyIEDB1544QWja4RpDKmnPm3atL9MmaJ4RNVuufz8/MDAwIyMjKqqqpdeeokMxTAEBgZ+/vnnFy9eXLRokRDCYrF0dXX5+fm99tprLS0tZCjcnOpFfW5u7oIFC6KiovQqCF7um2+++fjjj/39/bOzs8MedEUGDOQ8Gw309V08btz9X/yniIg1ISGKRzR+3igA6MgZo9NGj/5m3jwXHJEROACghBgFACXEKAAoIUYBQAkxCgBKiFEAUEKMAoASYhQAlBCjAKCEGAUAJcQoACjxM7oAANBTuL//aF9fIUSYn4vyjdEkAKCEi3oAUEKMAoASYhQAlBCjAKCEGAUAJcQoACghRgFACTEKAEqIUQBQQowCgBJiFACUEKMAoIQYBQAlxCgAKCFGAUAJMQoASohRAFBCjAKAEmIUAJQQowCghBgFACXEKAAoIUYBQAkxCgBKiFEAUEKMAoASYhQAlBCjAKCEGAUAJcQoACghRgFACTEKAEqIUQBQQowCgBJiFACUEKMAoIQYBQAlxCgAKPlfcQV7L6UNxYEAAACQelRYdHJka2l0UEtMIHJka2l0IDIwMjIuMDkuNAAAeJx7v2/tPQYg4GWAAEYgZgViFiBuYGRjSACJMUNoJiZ0mlNAA6xJEqTdsSQ/F4Qzk/1KcxWlOKEmyjBocwONZWRiYGJmYGFmEGEQd4NaBAasD93UDsyaOXMfiPPQbdn+tLRndjBJJHF7mDhQvQNMXAwAKlcf4vkBZEcAAADCelRYdE1PTCByZGtpdCAyMDIyLjA5LjQAAHicjVBJDsIwDLznFfOBRs4GzbFtCkKoqQSFP3Dn/8IRpGkPVLUTyct45LFAslu4vt6YTQchANp43ns8DRGJASlA258vEd3UtLnSjY843eFgeYJ9jWymccgVhQ6Vls7XdDygIpmIiUck/YKM1AmppPaeTM196/4ADQPn6hajZeAeQocTdq3Yx7AS95XbjjEUucl10cQJTNlc8bdlP5VOuGRfcqU8H51j8QFBAFejpZ6oWAAAAE96VFh0U01JTEVTIHJka2l0IDIwMjIuMDkuNAAAeJxzdnZ2dlOo0dA10jO1tDCw0NE10DPWsdY11DOytDQw0THQMzHVsTaACqOKomjRrAEAUF4PXOCSrzgAAAAASUVORK5CYII=\n"
          },
          "metadata": {},
          "execution_count": 54
        }
      ]
    }
  ]
}