Skip to content

Instantly share code, notes, and snippets.

@greglandrum
Created December 15, 2021 05:19
Show Gist options
  • Save greglandrum/8fd229bc6bf6c734d1c21da7f2bebebb to your computer and use it in GitHub Desktop.
Save greglandrum/8fd229bc6bf6c734d1c21da7f2bebebb to your computer and use it in GitHub Desktop.
LewisMartin_FixingValences.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "from rdkit import Chem\nfrom rdkit.Chem import rdChemReactions\nfrom rdkit.Chem.Draw import IPythonConsole\nfrom rdkit.Chem import Draw\nimport rdkit\nprint(rdkit.__version__)",
"execution_count": 1,
"outputs": [
{
"output_type": "stream",
"text": "2021.09.3\n",
"name": "stdout"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "m1 = Chem.MolFromSmiles('COC(=O)C1=C(C=CC=C1)C1=CC=C(C[N+]#[N]=[N-])C=C1',sanitize=False)\nm1",
"execution_count": 19,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 19,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7f74e0093460>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3de1hU5b4H8O9wR1CUi0DewBQBBVK8BWnbSzsR8JKMFngpcZtuC3adik7nJHrO03M097O9VJapaFmWDCkXBUXTBBU3oJSK1xQ1uSiIXB0uM/OePxYthuGiAzOzFsPv85e8a2bWb3zk67ve913vkjDGQAghpLNMhC6AEEK6N4pRQgjpEopRQgjpEopRQgjpEopRQoyEUqlUKBRCV9ETUYwS0r0plcpTp05FR0cPHDjw448/fvfdd1UqldBF9SwSWvBESHekVCpPnDghk8kOHDhQWlrKNVpaWtbX14eFhe3Zs8fKykrYCnsOilFCuhOlUpmVlSWTyeLj40tKSrhGd3f30NBQqVRqamo6a9assrKyiRMnpqSkODo6ClttD0ExSkg30HF6vvDCC/wrL1++HBQUdPfuXS8vr7S0tCFDhghUcg9CMUqIeD19eqorLi6eOXPmr7/+6urqeujQodGjRxuw5J6IYpQQ0eHTc9++fffv3+can5ie6mpqasLCwo4cOWJrayuTyWbMmKHnkns0ilFCxKLN9Bw6dGhISEjH6alSqUxMNFfdNDQ0vPHGG3v37rWwsIiLi4uIiNBj6T0bxSghAut0etbV1R09elQmk508efLGjRsWFhYaL2CMrV27du3atRKJZPXq1WvWrNHft+jRGCFECAqFIjMzMyoqytnZmf99HDp0aFRUVGZmZgdvfPz4cUJCwoIFC2xsbLh3SSSSU6dOtff6LVu2cH3VZcuWNTY26uGr9HTUGyVEAA0NDUOGDOFnjby9vaVSqVQqHTlyZHtv4fueiYmJ1dXV6m9cuHDhsGHDOjhdYmJieHi4XC6fM2fO3r17ra2tdfhdCMUoIQIYPXr0vXv3bG1tZ82a9ZRX7p1IT3Vnz54NDQ0tKyubMGFCSkqKk5OTDr4GAUAxSojhlZSUuLq69unT59GjR62nhji6Sk91V65cCQoKunPnzrBhw9LS0jr3IaQ1M6ELIKTHycnJATB27NjWGSqXy48dO6bb9OR5eXllZWUFBwfn5eVNmjTp0KFDY8aM6coHEg7FKCGGlpubC2DcuHEa7UuWLElISHj8+DEAiUQyceJEqVQaFhY2ePBgXZ3a1dU1IyNDKpUePnz4xRdfjI+PDwoK0tWH91i0wxMhhsb3RjXaa2pqHj9+7O3tHRsbe+PGjaysrHfffVeHGcqxtbVNSUmJjIysqamZNWvW9u3bdfv5PRCNjRJiaM7Ozg8ePCgoKHBzc1Nvv3r1qo2NzaBBgwxQA6MlpbpDMUqIQd25c8fNzc3R0ZHf3U5AO3fuXLFihUKhWLp06bZt28zMaJSvM+iinhCD4q7oWw+MCiIyMjIhIcHa2jouLi4sLIwbliXaohglxKC4+aXWA6NCmT179okTJ5ycnJKSkqZMmSKGPnK3QzFKiEGJLUYBTJgwISMjw83NLTs7e9KkSUVFRUJX1M1QjBJiOIyxc+fOQWQxCsDT0zMrK8vT01OhUNy7d0/ocroZilFCDOfGjRsVFRUDBgx45plnhK5Fk4uLy8yZM2/evHnw4EGha+lmKEYJMRwRXtGr++233yDi8kSLYpQQwxFPjNbV1f344483b97kW0Q74CB+FKOEGI54Vjvl5eW99tprr7zyCt/y+++/V1RUuLq6inDAQeQoRgkxEKVSmZeXJ5FIxNDda90v5iJ+/PjxgtXUbVGMEmIgV65cqa2tdXd3d3BwELqWNmJUPAMO3Q7FKCEG0t6OJIJoPbzQ3r5T5InoFlpiDOrq6o4cOSKTyYYMGWJtbR0TE2Nubi50UZrE092rqam5du2ahYWFj48P16JSqX799VcAtANpJ1CMkm6MT8+UlJSqqioANjY2tbW1mZmZCQkJvXv3FrrAFsQzv3Tu3DmVSuXr62tpacm1XL58ubq62t3dnR4u0gkUo6T74R+wkZSUxKUn/twi3t/ff+XKlenp6YGBgampqQMHDhS2VF5DQ8OFCxckEsno0aOFroWu6HWMYpR0Gx2kZ3h4uIeHB9cyatSomTNnXrx4cdKkSWlpaZ6ensKV3OzixYv19fWenp52dnZC19LuNL0YBhy6I4pRInZPmZ48d3f306dPz549+9SpUwEBAUlJSZMmTTJ41ZrEc0UPmqbXOV0/+J4Q3ZDL5cnJyYsWLerTpw//z5V7wMa1a9ee+Pa6ujqpVArA0tJy3759Bii4Y5GRkQA2b94sdCHs4cOHEomkV69ejY2NXEtDQ4OVlZVEInn06JGwtXVT1Bsl4tJB3zMiImL48OFP+TmWlpY//PCDk5PT1q1bw8PDS0tLV61apbeqn0w83b3c3FzG2OjRo/m97i9cuFBXV+fp6dm3b19ha+umKEaJKOgqPYuLi52dnbkHF5uamn7xxRceHh7vvvvuW2+9df369Y0bN7b3XHi9ksvlly9fNjU19fPzM/zZNbSeTRJPxHdTFKNESHx6tn4su1bpySktLZ08ebKPj8/3339vbW3NNUZHR9vb20dGRm7ZsqWoqGjPnj1WVlY6/hpPkpeX19jY6OfnZ2NjY+BTt9Z6NolitKuEHlUgPddHH33Ex4pEIhk/fvyGDRsKCgo6/YFZWVncZenkyZPLy8vVDx07doybIp86dWpFRUVXS9fS5s2bAURGRmq0Z2dnv//++0ql0pDFcIvArl69yrdwfeTTp08bsgxjQjFKhLFx40buupJ/LLtOPjY/P597sLu3t/ft27fVD124cIFLkFGjRv3xxx86Od3TuH37NvdNNeaX6urquGcpz5s3Ty6XG6aY4uJiAHZ2dnx2y+Vyc3NzU1PTmpoaw9RgfChGiTACAwMB7NmzR+efXFRU9NxzzwFwdXU9f/68+qGCggJuGemAAQN+++03nZ9a3d27dzdt2hQYGCiRSAA4ODj4+vrev39f/TVnzpxxdHQEMHHixNLSUr3Ww0lOTua65Oo1APD19TXA2Y0VxSgRgEKhsLGxkUgkZWVl+vj86urql19+GYCtrW1aWpr6ofLycm4Zab9+/U6ePKnzUxcUFGzYsGH8+PFcenI1BAcHDxgwAMDw4cN///139dfz3WcvLy+N7rM+rF69GsAHH3zAt7Q34ECeHsUoEQD3sIpnn31Wf6eor68PDw8HYGFh8d1336kfqqurmz9/Pv5cFKWT0925c0e97wnA2to6JCTkm2++qa6uZoyVlJT4+/tz3VKNUciioiLuDlFXV9dz587ppJ72zJw5E0B8fDzfsmjRIgBffvmlXs9r3ChGtXbkCEtObv4xP5+lpgpXTfe0c+dOAAsWLNDrWVQqVWxsLDd/FRsbq3Hovffe4w59+umnnT7FE9NTXU1NTXBwMAAbG5uUlBT1Q9XV1TNmzGiz+6xb/fv3B6A+j+fl5QUgJydHfyc1ehSjWps7l5mYsLNnm37cupXNnCloQd3QypUrAfzzn//kWwoLC+vq6vRxri1btnBrRd966y2FQqF+aNOmTdyhqKgorabLtUpPdQqFYvny5QBMTU2/+uor9UP19fUREREAzMzMduzYodV3fEp3794F4OjoyLdUVVWZmJhYWFjo6S+/h6AY1drcuSwggPn7M+5XkovR6Gi2eDF74w2hi+smuCWKv/zyC98yY8YMc3Pz9PR0fZxu//793DLSOXPmPH78WP2QTCbjlpE+zXR5p9NTw7p167i3x8TEqFQqvr2D7rOuFBYWZmVl8T+eOHECwLhx4/Rxrp6DYlRrc+eyL75g48ezLVsYo96o9urr6y0tLU1MTCorK/lGbpvLO3fu6OmkWVlZ7c2JHz9+nFtSOmXKlDZvKr99+3Z76dnpRUJxcXHcxtKvv/56Q0OD+qHPPvuM6yNHRkbyt73ryYYNGwCsXLlSr2cxehSjWps7l23dys6eZX37ssJCilGtcXfReHt78y23bt0C4OTkpNfzXr58eciQIQCGDRumsUz14sWL3BLOkSNH3r17l2tsnZ69evXqYnqqO3LkCLex9EsvvVRVVaV+6MCBA1z3efbs2bW1tV0/V3sWLFgAYOfOnfo7RU9AMao1LkYZY8uXs/DwFjGalsYOHhSwtO5h69atABYvXsy37Nu3D0BwcLC+T83Pibu4uGjMiRcWFnI387i4uERHR3Oz6pw+ffpEREQcOHBA54vks7OzuTmfcePGlZSUqB/iu88TJkx48OCBbs/LGzp0KIALFy7o6fN7CIpRrfExWlbGnJzYokVNMXrzJrO1ZWZmbPt2YQsUu6VLlwLYwo2JMMYYe//99wHoaTRQg/qceGrLNRbl5eWTJ082Nzfnup+67Xu25+bNm9yWqe7u7uo3aDLGbty4MWzYMG5lmK7u8lL34MEDjR3zSOdQjGqNj1HG2I4dzMSkuTe6bh0DGMBiYpjazAFpgXuMmvpEx5QpUwBorAHSH35O3NLSUmPF++nTp7m57MTERIPdoFlWVhYQEADA3t7+1KlT6oeKi4u5Z8y5uLhwG9x1nVKpzMzMjIqKcnJyGjx48Lx583TysT0ZxajW1GNUpWIBAS3GRnftYubmDGBLlrCWMweEMcZqa2vNzMzMzMz4GXOVSsXtJ1JUVGSwMlQqVUxMzPr16zXav/zyS40BB8OoqakJCQnhusDJ6suSGausrJw+fTqA3r17HzlypNOnUCgUx44de/PNN9UfWjdhwoQu104oRrVXXc3U19g9fMg0fv3T01nv3gxg06cztblowhhjp06dAvDcc8/xLVeuXAEwaNAgAavitR5wMBiFQrFixQpuSanGPUWNjY3Lli0zNzfvxMp8hULB9T1dXFz49HR3d4+KisrMzNRd+T0axWiXKBTslVfYqFHsz9ndJtnZrH9/BrCxY1nLmQPd0GpntZYLJQW2ceNGAH/729/4lj179gCYO3eugFXxfH19AZw5c0aoAtatW8eNzEZFRWksKdXqRiNKT0OiGO2SkhLm5cUANngwy89vcejWLebhwQDm7s5azhx0Xl0dS05mixYxV1emMe1x/jwD2Lp1zS1OToz7vVu2TESdYm5Qctu2bXxLdHQ0gE8++UTAqjitBxwEsXv3bm5J6eLFixu0HBii9BQExWhXlZezyZMZwPr1YxobBpWVscBABjB7e9Zy5kA7cjlLSmILF7I+fZqmsACmMUp2/jyzsWGOjoy/W9rJiWVksNhYNmYM+/BDduhQ5wvQoREjRgBQX2zEza7o6f4lrbQecBBKeno6t6R0+vTplU/xfyCfns7OzpSehkcxqgN1dWzBAgYwS0umsWFQbS0LCWEAs7JiMpnWH8v1Pe3smtPT25vFxrbRvT1/nrm4sP/+bxYS0tTi5MSys1l5OVu8mN25I4pL+4qKChMTEysrq/r6eq5F3zvmaWXTpk0Ali1bJnQhjDGWk5PDZaKvr++9e/fafE2b6Tl06FBKTwOjGNUNlYq99x4DmETCNDYMUijYihUMYKamzVP8HdAqPXlcjNbUsEGD2E8/MaZ2Ub9zJzPU0p0n+PnnnzVmhw2wY97TW7hwIQCNHUMEdOvWLa7zrrGklNJTbChGdWnTJmZiwgAWFaU5C7RuHZNImg61t6Q0L49FRLS4cvf3Z+vWsZb7/GqqrGTp6U0xyhhLSGADB7KamuYYFY/169cDWLVqFd+yY8cOAK+++qqAVfFaDzgI7uHDh9xjAuzt7U+ePEnpKU4UozomkzErKwawefM0+4C7dzctKV28uO0lpSdOtOh7XrvW0Ylqa5s6rTY2zNSUHTvWFKOMsZdfZh9/LMYYlUqlAHbt2sW3cEt81HfME0plZaWJiYmlpSU/4CAStbW1oaGhANSfKjpy5MjY2NhLly4JXR1hjGJUH44fb7oenzKFaTyD8ujRps7mtGltzJ4rFOxf/2I3b3b04RUV7Ntv2axZzNKyKXNNTdlf/sISEppj9MYNZm/PrK1ZTo64bqZyc3MDoP7L33rHPKEcP34cwPjx44UupA0KhcLLy6tv376DBw+OjY3N11gUQoRGMaoXFy+yQYMYwEaO1FxSmpPDnJ1Znz7s4kUtPpDve/bq1ZSeJiYsMJBt2tS0+J+/qOesXs0A9ssvLDBQLLullJaWcnfp8Hsnt7ljnlBaDziICvckPl3dD0p0i2JUXwoLmZ8fA9gzz7Bff21x6OZNzaVR7Xn0iH3zDQsJae57aqQnr7SUqT9k8/Fjtm0bW7OGAczcnMXFdfn7dFlaWhqASZMm8S3Z2dlouWOegFoPOIgHt0e9paUl7VEvTmYg+vHMMzhxAnPmICMDf/kLEhPx4otNh+7cwb598PFBv34AcP8+/vd/8fnnze+tqEByMmQyHD2K+noAMDFBYCCkUsyfD1fXNk7n6IiFC5t/tLbG8uVgDIxh7VosXYpr1/DnnuvC4LYZ5Z7YzsnNzdVoERBXDDfIIDa5ubkqlcrPz8/S0lLoWkgbTIQuwJj164f0dLz6Kioq8PLL+PHHpvZLl7B9Oz78sOnHigrExTX94dtvERoKFxcsWYKDB9HYiMBAbNqEe/dw6hSio9vO0PZIJFizBjt3wswM69dj6VI0Nur0G2qjdU5xwSqG5CovL799+3avXr24a2exEXPEEwDUG9UvS0t8/z3698eWLYiIQP/+mDoVAKZPR1ISlixBQEDTK9etQ2wsGhoAwNQU06YhLAyvvIL+/dv98OJi7N+Pykp89FFHNSxdCkdHvPYadu1CYSESEtC7t66+nxZaZ8Hbb7/t5+f317/+VYBqWsrOzmaM+fv7m5mJ8TeCYlTkxPiPxsiYmGDzZgwdihMnmq/r7ezwP/+DlStx7lxTy7PPQqlsunJfsABqt0RrKitDaipkMhw+DIUCtrZ45x1YW3dUw6xZOHECoaFIT8e0aTh4sKN01ofi4uKioiI7OztuH2LO6NGjub3oBSfynGo9HkJEhWLUQKKjERWFPx/qAwDLlmH7dmzciFmzACA0FMXFUNsKUlNREfbvh0yGU6egUgGAlRVCQhAWBpOnGJsZPx5ZWZgxAzk5eP55pKXBw6OL30kLCQkJAPz9/SXqfwWiIeYYFfmAAwHFqCFpBIiJCb76CtOmgfvltbKClVUb79Loe3KvnD4dUinmzEGfPloUMHQosrIQGoqsLAQEIDm5eUhBT/Lz82UymUwmu3z5spub25UrV+7evTt48GD9nlV7Yu7ucfvjjRkzRpwDDgQUo8Ly90d4OFavbuNQm+k5Y0Zn0lOdgwPS0zF/PtLSsHr1D++80yc4OLjzX6Adly5dSkhI4NKTa7G3t5fL5ffv3w8ICEhLS+OeIyISJSUlrQccxEPMEU+aCL3iqifasoVJpU1/Li9n/fsza2vN13C7mXBbQ4WEsG++0eWeoY2N7KOPzgIwMzPbrrsn8F26dCk2Ntbb25v/19WvX79FixYlJyc3NDRUVVVxs0m2traHDx/W1Um7LikpCcC0adOELqRts2fPBvD9998LXQhpF8WoAH78kcXENP+4dy8LCNB8zcmTbO5ctncva/kAc11a9+c60piYGFUXbhptnZ729vZ8eqq/sr6+/rXXXgNgYWEhnlz4+OOPAXz44YdCF9K2AQMGALjW8Q4LRFAUoz3arl27uI3WlyxZou1G61x6enl5PTE91alUqtjYWAASicQwT1R+Iu55ywkJCUIX0obi4mIAdnZ2Sq2eG0MMi2K0p+M3Wn/ppZeqnqLr27n01LB582YTExMAb7/9tsECQqFQ3NXY4IAxxhj3pEyNJy2LhMgHHAiHYpSw7Ozs/v37Axg7dmxJO0/g00l6qtu/f7+1tTWAOXPm6PXZR/xj2V1dXf38/DSOFhQUAHByctJfAV3BDTjEqI8BEfGhGCWMMXbz5k0PDw+02mhd5+mp7syZM46OjgAmTpxYWlraxU/ToFAofv7555UrV6pvcjx8+PBHjx6pvyw+Ph7AzJkzdXt2XQkKCgIg0/b5M8SwKEZJk7KyMu7pcvb29nv27ImNjVVf763D9FSXn5/PLSP18vLSyWW1et+TL97NzY3bIr71TNoHH3wAYPXq1V0/tT5wVwniHHAgPIpR0qympoZbRmpra8sFkIODgz7SU11RURF3S6irq+v58+c79yF8eqo/WJhPzw7eOHXqVAApKSmdO69ecQMOjo6OQhdCnoCW35NmNjY2iYmJW7Zsqa6uLikpkUqlL774oqmpqV5P6urqmpGRERYWduTIkcmTJ8tkMm7q/GmoVKozZ87IZLL4+PiSkhKu0c3NbdasWVKp9IUXXmjvjXV1dUePHpXJZNnZ2WZmZv7+/jr4Jromqo0ESUeEznFCGGOsvr4+IiICgIWFxR71Dajbwj8aU73v+TSPZZfL5QcOHAgPD++ttsnVrxq7aotGTEwMRDzgQHjUGyWiwKXnsGHD1q5du3jx4t9//33NmjVtvrK8vNzT05N7JAmAESNGSKVSqVTq6+vb3ofzfc+kpKSqqiqu0dvbWyqVhoeHexhyjxZtiGc/VtIxCWNM6BoIafbZZ5/94x//UKlUy5Yt+/LLL9vcj2PMmDEVFRWhoaFPeeXevdKTwxizt7evqKgoKipy1WqzbmJwFKNEdBITE8PDw+Vy+ezZs3/44QfrVnuplpeX29vbt/f2bp2evOvXr48YMWLgwIF//PGH0LWQJ6CLeiI6c+bMOX78eGhoaFJS0pQpU1JSUpxa7sPaZoZ2kJ4RERHDhw83ROm6Q1f03QjFKBGjiRMnZmRkBAUF/fvf/37++ecPHz7c3i52RpaePDHvJE00UIwSkfLy8jp79mxwcPD58+cnTZp06NChMWPG8Ef59ExMTKyuruYajSA9lUrlyZMnBw0aRNuMdicCrxQgpEPV1dXcMlJbW9vU1FS5XJ6cnLxo0SL1FUve3t6xsbHXr18XutjO07j5Kjo62sbGRiKRlJWVCV0aeTLqjRJRs7W1TUpKioyM/O6770JDQy0sLORyOQCJRDJhwgSpVBoWFjZkyBChy+wkru8pk8n279//4MEDrtHDw8PS0rK2tnbo0KEODg7CVkieBsUoETsLC4tvv/322Wefzc/PT0hI4K7cFy5cKM5nfjyNJ958FRcXB6D7Dk30NLTgiXQbjLHCwsKBAwcKXUhnKZX45ZdfTp6c/9VX/O0DHh4e3O0Dfn5+/AvLy8tDQ0Pz8vLi4+NDQkIEKpc8LYpRQvRMqURWFmQyxMejpKTKx8fu4kV3d/cObh9QKpV///vfv/76a1NT0y+++OLNN980fNXk6VGMEqIfSiVOnIBMhgMH8GffEyNGsPnzL4WF+bR/66pcLj98+LBMJnNwcPj8888BxMTE8A/OIiJEMUqITrXsezY1ursjNBRSKbS5dfWdd97x8fFZvny5QqF4/fXXv/76a+7BWURsaIqJEN3JzkZwMMrKmn709IRUCqkUPj7tvUMul6elpclksoMHD9bU1HCN48aN4wZM3dzcBgwYEBYWtnv37sLCwp9++kl9pRcRCeqNEtIFmZlIT0dVFUaMQEQETEzg7AwXl070PdH+7QM5OTkhISEPHjwYN27cwYMHuS3xiXhQjBLSWTEx2L0bq1bBwQGHDyMvDxkZMDfHoEHtvkUuR2oqZDLZw4fzjx0DIJFIuL5nWFiYm5tbe++7detWUFDQ9evX3d3dDx8+3F02WOkhKEYJ6ZScHLzwAvLzwS9fXbgQ1dVISmrjxXV1OHoUMhkSE1FdDaDKx+evvXppdfvAgwcPQkJCcnJyggICUjduxPjxOvsupGsoRgnplA0bcPw40tKaW06fRlAQ/rxCB9pITwDw9oZUioULof3tA7W1tf/x5ptr09OdHz/Gjz+ClpSKA8UoIZ2yahUaG/H1180t9+5h0CCUlsLWVrfp2YJSiVWrsG0bTE3x+edYsaJLn0Z0gWbqCekUe3tcv96ipaICZmaws4NMhoiIpkZdpSfP1BRffQV3d/znf2LlSty+jf/7P0gkuvlw0ikUo4R0ysiR2LMHjY3g13JmZGDECJibIzQUkydj9myEhWHwYL2cPSYGzs5Yvhzr16OkBNu3g5aUCocu6gnpFLkcfn6YPh2ffgpbW2RmQirF+vVYssRwNRw9innzUF2N6dPx00/o08dwpyZqKEYJ6ay7d/H22zh6FGZmcHDAf/0Xli0zdA25uQgJwf378PVFaioGDDB0AYRilJCuYgxyOXr1EqyAggIEBeHaNbi7Iy0NI0YIVklPZSJ0AYR0cxKJkBkKwN0dmZmYMAEFBdi/v6mxvBznzuHqVVA/Sf+oN0qIUXj8GHFxTcuwVq2CTIaRI5v2Rvn2WwQGCl2fMaMYJcS4fPQRUlNx5AicnQHgX//CJ5/g+nXQ80j0hi7qCTEucXFYu7YpQwG88w5cXJCcLGhNRo5ilBAjUlWF+/fx3HPNLRIJ/Pw07xQgOkUxSogRUSoBwKTl77WZGRQKQcrpIShGCTEi/fqhXz/NvufVqzq7FZW0hWKUEOMybx7WrUN9fdOPqam4dIn2gtIrmqknxLg8eoSQEJSXY8oUlJTg55+xbRtefVXosowZxSghRkelwunTuHYNvXtj6lQ4OQldkJGji3pCjM6OHViyBHI5FiygDDUAilFCjE52NgoK6DZQg6EYJcTo5OYCwNixQtfRU9DYKCHGRS6HnR1UKlRWwsZG6Gp6BOqNEmJc8vLQ2IhRoyhDDYZilBDjQlf0BkcxSohx4WJ03Dih6+hBKEYJMS45OQD1Rg2KppgIMSLV1ejbF+bmqKyEpaXQ1fQU1BslxIjk5kKlgp8fZaghUYwSYkRofkkIFKOEGBGKUSFQjBJiPCbk56/08XlA0/SGRVNMhBiJ8vJyR0dHa2vryspKMzMzocvpQag3SoiRyMnJYYyNGTOGMtTAKEYJMRI5OTkAxtEVvcFRjBJiJHJzcwGMpfklg6MYJcRIUG9UKDTFRIgxKCkpcXV1tbOze/TokUQiEbqcnoV6o4QYg+zsbABjx46lDDU8ilFCjAE3MBQHprgAAABxSURBVEpX9IKgGCXEGHADozS/JAiKUUKMwblz50AxKhCKUUK6vdu3b5eWljo5OQ0ZMkToWnoiutuBkG7Pzc3t1q1bf/zxh9CF9FC04IkQQrqELuoJIaRLKEYJIaRLKEYJIaRLKEYJIaRLKEYJIaRL/h9SzOrrqiYC7QAAAYN6VFh0cmRraXRQS0wgcmRraXQgMjAyMS4wOS4zAAB4nHu/b+09BiDgZUAAESAWBeIGRjaGDJAAIzMHgwaIwcQGoVnQ+QgarIGZES8Dt1oUGSYmdgUtkPUMLOwKYBlWiMB/kMXomrmB6hiZGJiYOZiYWBhYWDmYWNkY2Ng5mNg5GDg4gTwuBi5uDiZuHgYeXg4mXj4GPn4GfgEOZgFBIE+IQUiYg4mThUGYi0GEiY2Fk4OdjZWNm4dXSJhL/BfYDbDAyWmMcPA+Wn4AxLl7Jtvh1C5JMLuEUcVh2//N+0Fsk6sT7W/N+wBmP/Uxd9gU0mgLYv9gum+v2d1uD2K/PcPkIHpB0QHEPv441mG5vAOYPeP6JIcdKQxg9gPDToe7GZxg9fyhi+24b4iC2TLNvvufJJ4CsxeW8x04eKgXzBY6rHiAxXPnPrD5N2sOnHgUCnbD6qBFB84rq9iB2BL3zh9YlrUArObM4c8HMrfpgfUGv9+8X/7EWrD6iBn9Nt8Dk8BsMQAGwGBW5Sp3aQAAAhF6VFh0TU9MIHJka2l0IDIwMjEuMDkuMwAAeJx9VNuOWyEMfM9X8ANBHhuDediHXFbbqt1EatP9h773/1WbKDlshcoJFpwztsEzzi7F+HH+9vtPeg4+73Yp0X9+vff0IUS0e0+xSMfXt6+XdLodjo83p+uvy+1nYkoM9/HnM/Zwu74/3iCdkmQxNo+7l2wNBS1RpjE2V05XB1Yr0hzIuRAT1wVQPCJnrazmQDiwcsECWDwiMqQLBtBat2ILoI6IVj3S+FwFSgtcdRxyK8zNfEF+ACkLXBvxyG+iFkdlE9YFzkZtSm3EzRdE6MILXHdcyWpmiHhEFWVVQtAAsjWRuIhKbboqIYIVyl4XogE0FV5VBuzAPWXzExr7ndUArM6IoMVpY66kgQQ1pWXyckcqNxVEdKkwWxUIwYyLpisZDaiJlr7MX9Ml7TVTY+2hBlTiLhMynb68vSB9HL6/lKdXC6+aS1er97NABDOpw0GfDhYOLVe0qpGmGpeGf9Ps73n46RYUDq1qc636SrQ576s+oMENoSn3cWVzXcoC+Xo5f2q1e/Mdr5fz1nzx8NZivkmyNZJvUtnaBb7VrSnYZ920D59tkzj7tE3J8Nk3wbJPzMIMb2ASIIcBT0rDMDIpKraul0k5GEYngWCYOulAwqBNHHMY2ETiPU6f+OFhZhoQ536m4lFKmi4R9Z+rHfvH36Svd38BOz8Gv7nxXHsAAAEUelRYdFNNSUxFUyByZGtpdCAyMDIxLjA5LjMAAHicJZBdisMwDISvsrAvDXWE/iUT8uT39AClN+nhV8pCcODzWDOj9VqP87UtOtdj8blW/9/X8/N7ne9r/2wFuI/66Of7EJDkOXaBDFIeh4CnSo6dQZFRxsFgzlYaKuKs4yAgmcJNMmZaa9IVYyCgC+nMFoUy5yBAy1lzEIQmYpbI0qSujh0hkSO5VJZEVIiB2dGKEIb9E+MwGaUWp2y/yjutRvUAMZ0VfDfA6JwI5MjtuDvotPT7JUkFKxbgFG4l82SNNuhiFtp9xMKjsyKFccbtkNX7XgSqhI/OkyK9K/WoFQkgUhtqtchqU0p0sgY1Q7Qri4f52L5/BuVVCOUv73IAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "tf1 = rdChemReactions.ReactionFromSmarts('[#6:1]-[N+:2]#[N:3]=[N-:4]>>[#6:1]-[N+0:2]=[N+1:3]=[N-:4]')\nwhile tf1.RunReactantInPlace(m1):\n pass\nChem.SanitizeMol(m1)\nm1",
"execution_count": 20,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 20,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7f74e0093460>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxUVf8H8O9lhhlWRXEnUeHnRrmSmCjKwJ0BFcw0zQ01E9wSUTNc0hK19LGeMkuzTJ60TNGeTHKJGQYHEKVQQHDpMVFyQQSVWEVm5vz+ODiMgIDAzL0zfN+v/ohzZ+790is/3nPuuecwhBBACCHUWBZcF4AQQqYNYxQhhJoEYxQhhJoEYxQhhJoEYxQhs1VSUsJ1CS0CxihCZkWj0Zw7d27Lli1SqbR9+/bLli1TKBRcF2XmGJzwhJAZ+Ouvv2JiYuRyuVKpLCwspI1CoVCtVltbWx85ckQmk3FboRnDGEXIVBUXF589ezY6Ojo6Ovr69eu6dhcXF5ZlWZaVSqVr16794osvxGLxoUOHAgMDOazWjGGMImRKNBpNWlqaQqFQKBQqlaqiooK2Ozo6+vj4sCzr5+fXrVs33ecJIUuWLNm+fbtIJDp8+DAmqSFgjCJkArKysmh0yuXygoIC2igUCgcMGMCybEBAgKenp4VF7Y86CCFhYWGff/65SCQ6dOjQuHHjjFh4i4AxihBP1dtn9/Pza9WqVUNORQhZunTptm3bRCJRVFTUq6++arCqWyKMUYR45Fl99nbt2kkkEpZl/f39nZ2d6z5JQUGBlZWVlZVVtfZVq1Zt3rxZJBIdPHhw/PjxBvkFWiSMUYR4Yc+ePb/88ktcXFxRURFtEYvFw4cPl0qlMpls4MCBz+qzU/r5Gx8f/6xbztWrV3/00UeYpM1LyHUBCCHYu3fv119/nZycDM/ZZ691npNIJMrKyqr18x9++CHDMB9++OHkyZMPHjz42muvNe8v0jLh3ShCHCOEODk55eTkbNq0adasWU5OTnV/vt4xU5lM1rp16zrO8N57723atMnS0vLAgQMTJkxonl+jBcMYRYhjmZmZ/fr169Kly61btxiGqfUzzzvPqV5r167duHGjQCDYu3fvtGnTmuHXaMGwU48Qx+jLmr6+vjUz9FnznNzd3eud51S3DRs2MAyzYcOGmTNnEkKmT5/exN+iJcMYRYhjsbGxAODr61utvX///hkZGbof3dzcZDKZVCodNWqUra1t068bERHBMExERMSsWbMIITNmzGj6OVsmjFGEuKRWq+Pj46G2GO3Ro0dOTk7D5zk1wvr16xmGWb9+/ezZswEAk7RxcGwUIS4lJSUNHz68T58+ly9frnbowYMHDg4OjeuzP5ctW7asXLlSIBBERkYGBQUZ+nLmB+9GEeKSbmC05qG2bdsap4bw8HAAWLly5ZtvvkkImTlzpnGuazYwRhHi0rMGRo0sPDycYZjw8PA5c+YQQmbNmsVtPaYFO/UIcaa0tLRt27ZqtTovL69NmzZclwNbt2599913LSwsvv32WzpaihoCV79HiDPx8fHl5eXu7u58yFAAWLFixdatW7Va7Zw5c7766iuuyzEZGKMIcYYnPXp977zzzscff0wIWbhw4c6dO7kuxzRgjCLEmTqeL3Fo+fLln3zyCSFk0aJFO3bs4LocE4AxihA37t+/f+HCBSsrK09PT65rqW7ZsmU0SRMSEriuxQRgjCLEjdjYWK1WO2LECGtra65rqcXSpUvbtGlz4MCBq1evcl0L3+GEJ2SGcnJy5HI5wzBDhgzp06cP1+XUjocDo/oyMjIePnzYpUuXnj17cl0L32GMIjNRVlZ2+vRpupDH+fPnCSGdO3e+d+/etm3bFi1axHV1teDVwOi+ffssLS1Hjx6tW2GPprxUKuW0LtOAMYpMGCEkPT2drlucmJj46NEj2m5nZ0fX7zh06NDixYsJIW+//Ta3pVaTnZ2dlZXl4OAwePBgrmsBAHj//fevX7+ekpLi7u5OW3h+s8wrGKPI9Ny7d0+lUikUimPHjt2+fZs2WlhY0OXjWJb18vISi8UA4OPjs2DBgtDQUELI4sWLOa36KXK5HAB8fHwEAgHXtcC1a9euX7/etm3bQYMG0Ra1Wk0fLkkkEk5LMw0Yo8g01Oyz0/ZOnTp5eXmxLBsYGNi5c+dq35o3bx7DMPPnz1+yZAkhJDQ01OiF145X93q0GB8fH90yKGfPni0sLOzbt+8LL7zAaWmmAWMU8VpWVlZ0dPSvv/6q32e3sbHx9PSkN56DBw9+1orxVEhICE3SsLAwQsiSJUuMUnhdCCFxcXEAwLIs17UA1JbptIUn5fEfxijindzc3Pj4+Ib02RsoODiYYZh58+aFhYWVlZWtXLnSMIU31IULF3Jzc52cnHr16sVtJaCX6TVjlCc3y/yHMYp4oe4+e0BAQEBAQFMWjps7dy7DMCEhIatWrSKErFq1qpkKbwxePQRPS0vLy8tzdnbWTWwqKSlJTk4WCAQjR47ktjZTgTGKuJeVleXm5lZeXk5/tLe3l0gkdH/2Rt+vZWdnz5o1KzIyskePHrTlrbfeYhgmODh49erVhJDVq1c3T/XPj1f3ejX77yqV6vHjx0OHDuXJgin8hzGKOHbz5s3+/fsLBAJdn33kyJEikaiJp12+fLlKpfL19Y2Li9NtmTlnzhwACA4OXrNmDSFkzZo1Ta3++fHtITgOjDYDghCnIiMjASAgIKB5T1tQUPDKK68AgLOz87Vr1/QP7d+/n04z2rBhQ/NetCHozksvvvii8S9dU3l5ua2tLcMwOTk5usYBAwYAgFKp5LAw04Lv1COOGWigsHXr1gqFYtSoUX///bdEIrl27Zru0NSpU7///nuhULh27doNGzY073Xrxase/ZkzZ0pKSl566aVOnTrRlvz8/IyMDCsrK/qXEGoIjFHEMcNN/bG1tT127Ji3t3fNJJ0yZcq+ffuEQuG6desiIiKa/dJ1eFaM5uXlabVaY1ZSazEKhUKr1Xp5efFzwRSe4vp2GLVomZmZANCpUyetVmugSxQXF3t7ewNA165dr169qn/owIEDQqEQAD744AMDXV2fRqNRqVRCoVAoFBYUFOgfunnzpqur65QpUyoqKoxQiQ5doy86OlrXMnfuXADYvHmzMcswdRijiEvbtm0DgBkzZuha8vLyCgsLm/cqxcXF9HlOp06dLl26pH/o4MGDNEnDw8Ob96I6OTk5UVFRISEh9CUrJycna2vrU6dO6X8mMTHRzs4OAGbOnKlWqw1USTWFhYWWlpZCofCff/7RNdKJDSkpKcapwTxgjCIujRs3DgAiIyN1LatWrRIKhZ9//nnzXqikpMTHx4cm6cWLF/UPRUVFNXuSFhUVRUdHL168uNoyfS4uLm5ubgBgb2+fkJCg/5XExER7e3sAeOONN4xzT3r06FEAGD58uK6Fjnu0adPGaFFuHjBGEWcqKirosmx///23rnHIkCEAcPLkyWa/nC5JO3bsWC1JDx06ZGlpCQDvvvtuo8+v0WhSUlI2b97Msqz+S1a2trYsy27evJne4qnVarp9sa2tbbWn4adPn27VqhUATJ482QhJSt+LXbduna5l165dADBx4kRDX9rMYIwiziQlJQFA7969dS0PHz4UCAQikai4uNgQVywpKaHPsjp27JiZmal/SJekK1aseK5z0j57UFCQ/ktWdBpseHi4XC4vLy+v9hWNRkOT1MbGJjY2Vv9QUlKS0ZL0pZdeAgCVSqVrmTx5MgDs2LHDoNc1PxijiDN0stHChQt1LT/99BMAjBo1ynAXLS0tfVaSHj58mCbpO++8U/dJiouL5XJ5eHi4bnVOXZ89JCQkKirqwYMHdZ9Bo9HQjeBrJmlKSgp9fWjSpEmGS9K7d+8yDGNjY/Po0SPaotVq27dvDwD/+9//DHRRc4UxijhDH6D/9NNPupaFCxcCQEREhEGvW1paSqepdujQISMjQ//Qr7/+Svvjy5cvr/YttVqt67Prv2RlZ2dH++zVBgrqpdVqFyxYQJNUoVDoH0pJSaH3tq+//vrjx48b92vW7fvvvweA0aNH61pSU1MBwNnZ2RCXM28Yo+YmOpocP171Y1oaefpPKF+UlJSIxWKBQHD//n1dY+/evQEgKSnJ0Fd/9OhRQEAATdILFy7oH9Il6bJlywghd+7c+e677yZNmvSsPntTYk6r1dK/OWxsbORyuf6hc+fO0StOnDjREEn65ptvAgDdkp7aunUrALz55pvNfi2zhzFqbiQSIhSStLTKH7dsIW+8wWlBz3Dy5EkAGDJkiK7l1q1b9BG2ge6/qikvLw8MDASA9u3bV0vSI0eO0PtN2snV6dmz56JFi44cOaI/Q6iJdElqbW0dExOjf+j8+fOOjo4AEBAQoOt6N5fFixe3a9cuNTVV1+Lv7w8A33//ffNeqCXAGDV5V6+SL78k48cTOl1HIiFeXmTECELns/M2RlesWAEAK1eu1LX85z//AYDAwECj1VBWVjZ69Ggal7m5ufqHvvvuO3t7e0tLS9pn/+yzz7KysgxUhlarpZvuicVi/ZnwRC9Jx44d2+xJqtFodG89lJeX29nZMQxz586d5r1KS4AxapKKiohcTkJDSY8eBKDyn549CSFEIiF795J+/cju3YTwOEbptj/6PdmgoCAA+Oyzz4xZRnl5+bhx49avX1+t/dtvvwUAX19f48yg1Gq1dNM9sVh89OhR/UOpqak0SceMGdPsSaqjUqkA4KWXXjLQ+c0bxqjJUKtJSgrZvJmwLLG0rEpPR0cyaRLZtYvcuEEIIRIJ+eEHolQSR0dy7x5PYzQ/P9/CwsLKyqq0tFTX6OTkBADVnp4bQa1BOW3aNABo9rcA6qDVaummeyKRqGaStmvXjj4RKisrM8TV161bBwB0nxX0vDBG+e7aNbJrF5k0iTg4VEWnUEjc3Ul4OElIIBrNU5+nMUoImT6dzJ3L0xiNioqi93q6losXL9JJSIZ7ub7htFotfXGz2pujRrgu3XRPJBL98ssv+ofS0tIMmqQ1X65HDYcxylOlpSQ4mHTrVhWdAMTNjYSFkWPHSB2T03UxeucOadOGTJvGxxidN28eAHz44Ye6ls8//xwApk+fzmFVOhkZGQDQpUsX419aq9XSl4tEItGRI0f0D+mS1N/fv3mT9ObNm3TBlGZ8dNaiYIzyV/fuBIC0a1fZZ8/Orv8rDx8Sb+/KGCWEbNtGLCz4GKOurq4AkJycrGt59dVXAWDPnj0cVqXz6aefAsDMmTO5KoBuuicSiX7++Wf99vT0dDp5YDcd+W4C/WmwlpaWrq6ukydPbuI5WyyMUc78/DNZsoToXhS8epXonlrv3EkIISdOkHPnqvfZa6o2ZuruXhWjajUZNIi88Qb58Ucyf379pzKOGzduAICDg4NuUFKtVjs4OADA9evXOS2tEp1S+t1333FYA910r2aSZmRkvP/++40+7dWrV3fs2DF+/Hi6mgElEonwPfqmwBjlAH3cuno1YRii69fGxxNnZ/L4MSksJBMmkMJCUvcj4v/9j3z5JXn1VdKqVVWvXyQiO3cS/Xe4S0tJVhaxtSUAZP58woOBR7J7924AeO2113QtZ86cobMyOaxKp6Kigr7VfvPmTW4roZvuWVpa/ve//23KeYqKiuRyeWhoqG53v2qvrlZb/BQ9L9zSzkjUakhPB4UCFAq4cwcuXgQAGDsWtm6FN94AF5fKj6Wmwi+/wOXLsHkzzJkDrq5PnaS4GM6ehehoiI6G69er2l1cgGWBZUEmA72bDAAAa2vo0QOOHYOxY+Grr6CiAr7+Giw43fSA53uo/f7774WFhX369HnhhRe4rWTTpk0Mw2zatGnSpEn79u2bOnVqw7+r0WjS0tLohtUqlaqiooK2Ozo6+vj4sCzr5+en2+kPNRHGqGFdvgxyOcTEgEoFxcWVjVZWkJsLANCjB/TrBwsXwsmTlYc8PMDDA/76CzZtqmzRaCAtrTJ/VSp48scB2rUDiQRYFvz8oN4/DqNGwfHjMHYsfPstEALffMNZkhJC6K4hNWOUJ9sTKRQK4E2mb9y4kWGYjRs3BgUFkSfTsOqQlZVFo1MulxcUFNBGoVBId10NCAjw9PS04PZvUbPE9e2wGcrPJ1FRJCSk+nN2FxcSEkKioghd3H31arJ4MSkqIl27kqioyk49petj7d1L7O2rziAWE4mEfPQRSUlpzEBnfDyxsyMAZM4czsZJ09PT4emH4GVlZVZWVhYWFnl5edzU9LSRI0cCQLURSW6tXbsWAAQCQa1vahYVFR09ejQkJKR79+76f7Sxz240eDfaPHR99uhoOHMGdFuTtW8P3t7AsjB6NHTtWssX7ezg008hLAy++aaqUdcx79oVioqq+ux+ftCqVeOL9PKC48dhzBjYswcIgd27ObgnrbkPaHx8/KNHj9zd3elsHm6VlpYmJycLBIJRo0ZxXUuViIgIhmEiIiJmzZpFCJkxY8az+uzt2rWTSCQsy/r7+zs7O3NbdsuBMdpU27fDyZOgUkFJSWWLlRV4eYFUCjIZ9O8PDFPPGSZOhN274d//ruXQ8OFw6xY4OTVbtV5ecOIEjBkDkZFQVgb79oHQuP8L1Oy/nzp1CnjTo4+Pjy8vL/fw8KArfvLH+vXrASAiImL27Nnbt2+/dOlS8ZNBIrFY7OPjI5VKZTLZwIEDsc/OAa5vh03eyy/X3mev25UrJDSULF5c+eNffxErK2K0lR4TEyvHCqZMIcbcibLWh+CPHz9OSEjgyVLB77zzDgCsXr2a60Jqt2rVqlatWjEMA3p9dpwzzzmM0fqdOUNYlpw7V/mjWk1YlujG8Q4dInv3kpyc+s+jGzOl8+r9/Ij+tr4ffURGjGju0p8tIaEySZcuvWu0/csSEhIAoG/fvsa5XCMMHDgQAKotR88fdMGUoUOH3r59m+taUBWM0fodPUosLcmQIZWPZSoqCACJjSU//kh+/JFcu1bXd8vLSVwcWb2aDBlCLCyqHhZ16EA2bjRO+c+UmEg8PR/Y2XU09LY/FRUV9IUZNzc3gUDg7+9vuGs1Ra0LpvAKfVK/fft2rgtBT8Gx0Qbp3x8EAti5ExYtqmwRiSqf9uhtJ1ElK6tyilJMDPzzT2WjUAjDhkFgILAsDBrE8eRNABg+HD7++Iq/f1lUVBQhZP/+/cJmHSi9fPlyTExMTEyMSqUqeTJybG1tLZfL9+/fX+/cHeOLjY3VarUjRoywtrbmupZaEEKUSiXwZhwZ6WCMNgjDwLZt4O8PEyYAXRC9R4/qT37y8yEuDhQK+O03yM6uatc9Z/f3B3t7o5Zdr2HDhimVSqlUeujQIZqkdE+3RsvPz4+Li1MoFL/99lu23n8FFxcXlmVZlr1y5cq6detmzpxZUVFBd8fkD15NX60pMzPz7t27Xbp06du3L9e1oKdgjDaUhwdMnAjLlsG+fdUPJSVBaCikplbNc+rYEaTSyn86dzZypc/H3d1dLpfLZLLDhw8DQCOSVK1Wp6enKxSK6OjoM2fOaJ/8V2jfvr23tzfLsmPGjNF/I8jKyurdd9+dM2cOAPAqSXk18b4mXr3rhZ7C8aCCKTh6lLz8MiGE3L9P2rcnCgUBILduVX3g4kUCQKytCcuSzZtJSgovXl1/Lo3YQO3atWu7du2aNGlSK725rNbW1nSbzJSUlDpWDqW7p1lYWERGRjbb79A0NRdM4Rs+LJiCaoUxWj9djBJCvvmGDBhQPUYJIbGxhK+PJRqqIdv+5OXlRUVFhYSEVJvarZt8U1RU1MDLffzxxwDAMMyOHTua75dovG+++QYAJkyYwHUhtePPgimoJozR+unHqEZDhg2rJUbNQ63b/lRUVCQkJISHh7u7u+tP7e7QocOkSZN27drV6D/Yn3zyCU3SL7/8svl+iUaaMmUKAPChklolJiYCQJ8+fbguBNUCx0br1749eHhU/ruFBezYAStWgFjMaU2GMXDgQIVCIZVKjx8/7ufnN2bMGKVSmZCQUFpaSj9gY2MzcuRIqVQqlUr79evXxMstW7aMYZhly5bR3dzoPsOcIE8WTOHtyCMOjPIa1zmOeIfek+pv0e7i4hIaGiqXyw2xC9Cnn37KMAzDMF988UWzn7warVabmpqaU+NlibS0NABwcnIydAGNxsMFU5AOxiiqRXp6ulKpnD179g8//FBtA3dD2LlzJ01SA00sz83NpUO6dPPRrVu3VvsAHaidPXu2Ia7edCUlJWKxWCAQPHjwgOtaUC2wU49q0b9/fwCQSCTGudz8+fMZhlmwYEFoaCghhO4z3ESlpaXx8fFyuTwmJiYzM1PX7uzsXHNGF89njKpUqvLy8qFDh/JtwRREYYwiXpg3bx7DMPPnz1+yZAkhhO4z3AhZWVnR0dG//vprYmLio0ePaKONjY2npyed/z948GDm6UW36NooYMS/Np4Xz1MeYYwivggJCaFJGhYWRgih+ww3RG5ubnx8vEKhOHbs2O3bt2mjhYUFXfKdZVkvLy/xs58JJicnFxcXv/jii07NuCJhs8IY5TmMUcQjwcHBDMPMmzcvLCysrKyM7jNcq7KystOnT9N1i8+fP08Ioe2dOnXy8vIKCAgYO3YsnbxVh7t378rlcjrvirchdf/+/QsXLlhZWQ0bNozrWlDtMEYRv8ydO5dhmJCQkFWrVpEn+wzXNGLEiPPnz9N/t7e3l0gkdN3iXr161X3+WvM3LCyMw+lWdVMoFFqt1svLi58LpiDAGEU89NZbbzEMExwcTJdPpvsMV8OyrKWlJY3OV155pe51AAgh6enpMTExcrlcf8zU1tbW29tbKpWOGzeu2ubD/IE9ev5jdL0hhHhlz549wcHBWq1248aNa9asacQZ7t27p1Kpao6ZDho0qCFjpjzh6uqalZX1xx9/vPzyy1zXgmqHMYr468cffwwKCtJoNBs2bHjvvfca8pW6x0xZlg0MDOzM80W39GRnZ3fv3t3BwSE/P18gEHBdDqodduoRf02dOpVhmKCgoLVr15In+wzXqhHznExCTEwMAPj6+mKG8hnGKOI1umJIUFDQunXrCCHr1q3THWriPCeTgAOjJgE79cgEHDx4cMaMGWq1es2aNd7e3nXMcwoICKALp5oBQkjnzp1zc3P//PPPemcgIA5hjCLTQMdJxWKxbrkpW1vbYcOG0RtPd3d3bsszhPT09IEDBzo5Od26dYvrWlBdsFOPTMPUqVMtLS0LCgp2794tk8mkUumwYcOadw8+vqGbmshkMq4LQfUw5/8LkZl5/fXXAWDu3LlcF2JYxcXFp06dcnV1xYFRU4ExihD36EKodMw3ISGhvLx8+fLlCQkJDMP4+PhwXR2qB8YoQpzJzs6mS/kplcr79+/TRqFQ6OnpKRAI6IIpJjTLtcXCGEXIuIqLVfHxh0+ckMvlf/75p67Z1dWVvtvq4+PTunXrDz74AABGjBjBWZ2owfBJPUKGp9VCaiooFKBQQEJC5JAhcxITQW+yQWBgoJubm/43Ll++HBAQIBAITp8+rb+hC+IhvBtFyGBu3AC5HGJiQKmEBw8qG4VC7+7dP2BZmUzm4eHxrNeT2rRpIxaLL1++7OvrGxsbi0nKZ3g3ilCzKimBM2cqbzzPnatqd3EBlq38p769QG7evBkTE3Pq1Km0tLTMzMw+ffoolUocJOUtjFGEmubxYyguBlvbyk23R46EhITKQ23agI8PSKUgk0F9C/GVlpYmJSVVe0FLpVK9/fbbGRkZmKR8hp16hBqrpASWLIGDB6FVKygogPHj4csvwc8PCAGZDKRSGDIE6lxSRKPRnDt3Ti6Xy+XypKSkiooK2t66dWuJRCKTydzc3GJjY319fTMyMiQSiVKp7NKli1F+N/Qc8G4UocaaOhXu3IGoKOjYER48gGnTAABOnqz/i9nZBfHxIdHRsbGxD56MmQqFQg8PD/qw3sPDQ/8Frby8PJZlL1y40Lt3b0xSHsIYRahR7t4FJye4fBl0i4bcugU9ekBmJvTuXcvnq42ZikSdLC1zS0pcXFzYJ+rYP/nhw4cymSwlJaVXr15xcXGYpLyCnXqEGuXKFXBwAP2Fl154AZyd4cqVqhjVaCAlpfJh/dmz8KTPDg4O4OOzf9y47l5eLi4uDblamzZtTp48ybJsWlqaj49PolLZDpOUNzBGEWqU4mKwsaneaGcHhYXw8CEcOgRyOcTGwsOHlYeEQhg+vHLM1MMDBILnfcfT0dExLi5OJpO5icVtfXwgNhb4uiN0S4OdeoQa5dIlGDAAiorAyqqyRaMBR0c4fhw6doT/+7/KRt08J6kUHByaftnChw/t/Pws/vgDevYEpRJeeKHp50RNhDGKUKNoNODqCkuWwNKllS2RkbByJWRng5UVLFgAgwaBVFrvPKfGKCgAPz/4/Xfo3h2USoNcAj0PjFGEGuvUKRg3DiZOhEGD4OJF2L8fDhyAsWONcel//gE/P0hOhm7dIC4Ok5RbGKMINcHt2xAVBXfvQocOMHEidO9uvEv/8w/4+8PZs+DsDHFx0LBHVcgQMEYRMlmYpPxgwXUBCKHGat0aFAoYNQr+/hskErh2jeuCWiiMUYRMma0tHDsG3t6YpBzCGEXIxNnawtGjMGIElJRAURHX1bREGKMImT57ezhxAlQqGDgQNBrYsgUGDoQOHWDwYPj3v0Gr5bo+M4ePmBAyLwsXQkIC7NwJffvChQswbx4EBsInn3BdljnDGEXIjNy5A926wfnz0K9fZUtyMnh5QU4OODpyWpk5w049QmYkMxM6dKjKUAAYOhRsbODSJe5qMn8YowiZkfz8WnYoadcO7t3jopqWAmMUITPStSvk5DzVotVCTg5068ZRQS0CxihCZmTAACgvhxMnqlp+/hnEYnjxRe5qMn+43ihCZqRVK/jXv2D6dPjgA+jXD1JTISICtm0Da2uuKzNn+KQeIbMTGws//AA3b0K3bjBzJowcyXVBZg7vRhEyL4TADz/AkCGwaxdYWnJdTYuAd6MImZfMTOjXD7p0gdu3uS6lpcBHTAiZF4UCAIBlua6jBcEYRci8xMYCAPj6cl1HC4KdeoTMiFoNjo5QWAg3b+Jud0aDd6MImeJihLIAAAFiSURBVJHkZCgshD59MEONCWMUITNCe/Q4MGpcGKMImREcGOUCjo0iZC5KS6FtW1CrIS+vlgVKkMHg3ShC5kKlgvJyePllzFAjwxhFyFxgj54jGKMImQuMUY5gjCJkDvLy8mY+fixnWfD05LqWFgdjFCFzoFQq9126tJVhwMqK61paHIxRhMxBbGwsAPhij54LGKMImQOFQgEYoxzBeaMImbwbN2706NHDwcEhPz9fIBBwXU6Lg3ejCJk8uVwOAL6+vpihnMAYRcjk4cAot7BTj5BpI4R07NgxLy/vzz//7NWrF9fltER4N4qQaUtPT8/Ly3NycsIM5QrGKEKmjfboZTIZ14W0XBijCJk2HBjlHI6NImTC1Gp127Zti4uLb9++3blzZ67LaaFwn3qETJhQKLx69ervv/+OGcohvBtFCKEmwbFRhBBqEoxRhBBqEoxRhBBqEoxRhBBqEoxRhBBqkv8Hqt1CDhDtBIYAAAGLelRYdHJka2l0UEtMIHJka2l0IDIwMjEuMDkuMwAAeJx7v2/tPQYg4GVAABEgFgXiBkY2hgQgzcjMwaABpJmZ2CA0C5zvAOEj6AwQzcyIl4FbLVQGbCsTE7sCmM/MrqAFcgUTC4TxH2Qvul5uBkagCgUmZg0mJhYFFtYMJla2BDb2DCZ2jgQOTiCPS4GLO4OJmyeBhzeDiZePgY+fgV9Ag0lAUIOJVyhBSDiDiZMlQZgrQYSJjYWTg52NlY2bh1dImEv8F8hueNCUrfB2OHyg7gCIY6Cc4qD3VhHMXv9H3qFmwZ79IPaLbW32J3jfg9mZ5RYOGSH79oDYFxa+sO9zaLEHsRed5nG481nOAcS+m5Pp0KRlAWavPznNQVX6JVhNUE+Xw/Pze+1A7PTavXaty5TB4gG3Dfcr8d4Es/klOA9cn78MzD45SfHAXCuJfSA2i0ndAduXXGA3OOsuPNDUxQxW86D2/IHY+q+2IHbdub8HQkN7wOon3ti139hxFli9H6/3ngvL/cFsMQBRAmahUlixhAAAAhF6VFh0TU9MIHJka2l0IDIwMjEuMDkuMwAAeJx9lEtuGzEMhvc+BS/ggfiQKC6y8CNIizY20Ka5Q/e9P/pzDHuURuiMRUgzHykN+dM7yuvH+dvvP/S45LzbEZX//CKC3rWUsnulnNDx+eXrhU5vh+P9yen66/L2k6SQMHxwf2QPb9fX+xOmE+nCXbg67XWJXipXKktZr81V6AqwNmvWaS9L1cKFJ6AioiwW4QGQF3NNl8+gISIvpYrl1rz05q3MwLpG7F2kIyLei4raBGwAEYfdtGdoVUESPnO+BmSkUiLP2rVN4/U1Oc1cWnLdvcaMC3C2OAtjOxxAmG2WGi4rqBING5dFvbj4DMyyrO9L5PumIW32JSwAkZIW5lWxd4vqrDMy64K6sUmJJFG/dJmQdiOzMN0yOltvbfpBWZpUjXNITbSaNJmijS60ryhKQxoBVNauMiM9yYayGMML23N0H0k6fXl5Yno/fH+yh1dPL19C2WrLo4hb6fGv2/7mJw+/LF7KtIZDBJgJOmDeAuWW7MKiSCFmnUHOFPF8OX9os1vjHa+X89Z4ecvWXliQbk2EBdnWKoxl3RrCMNome8PwTd2G0TcRG0ZsWjUMHjXJ8GYetGdpWAaRWRrWQUyWhm0QDa+mDtrg1bRBApKGfSi1pOE+lNHWODEUKJfI/lAHy3M/tpI8crYej/kfs53r+18k5ru/YloFxzD9tMkAAAERelRYdFNNSUxFUyByZGtpdCAyMDIxLjA5LjMAAHicHZA7bgMxDAWvknKNSIT4J7Fwtb19gMCVruHDh9xCKkYDPj5d7+t4vh8b9zE37b2P6/X8e/1+6pqfx95UZ2/8+R4MGIQ2JkPGUhong5qYjkmgvHCcBJIZYyKIcz2cCCWKNQlzu5UIohxzwSImZmsr0EVG2cxUgxewr5U+FhgnlVK6pbhySZbqNalSUWhlkQp3vUmlafRwlLDO620dk5qpkFEzrSDzHq/IwZU4DSQEqRBihjdySMZefgG5rMhi1UzTuQvRUsR7s4XEydxiIKJK18TM5FE/Emzcf2Xi5AXCXWuUgCNVYJUnRGnAlKXe5Vt9fP8BrldVeAP7oCgAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "m2 = Chem.MolFromSmiles('NC(N)=[NH]C1=NC(CSCC[NH]=CNS(=O)(=O)C2=CC=C(Br)C=C2)=CS1',sanitize=False)\nm2",
"execution_count": 15,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 15,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7f74e0093a00>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxUVf8H8O+wLy7sCCqbCwpuZK6EjybmU1Jp/sglsczCrSg1nawUdyfzSSy3waeSLDMsS1Ircc3HJR0QZBGUVRjZdwaY9fv7407DKsJsd4Tv+8XLF3Pn3HvOqHw4995zz+EgIhBCCFGXEdsNIISQJxvFKCGEaIRilBBCNEIxSgghGqEYJYQQjVCMEkKIRihGCSFEIxSjhBCiEYpRQgjRCMUoIYRohGKUEEI0QjFKCCEaoRglhBCNUIwSQohGKEYJIUQjFKOk6/vhBzh2rPFlSgqsW8dea0iXQzFKur5bt+DWrcaXDx/Cr7+y1xrS5Ziw3QBCuqmffwY7O5gyRfkyORn+9z9YulQvdRcXw5kzUFQEXl4wYwZYWeml1i6LeqOkW6iogPR05Vd+PtutAQCAU6fgf/9rfHn3Lnz7rV4qPn8eBg+GixdBLIavvwYfH8jM1EvFXRb1Rkm38OefkJ6u/L6qitWmsEsigYULYfduWLRIuWX5clixAv74g9VmPdkoRkm3MHcu/Oc/yu9jYyEsjNXW/EMigZoa5fcNDXqpMjUVampg4cLGLWFhMHIkyOVgbKyXFnRBFKOEsGbvXvjuO+X3dXUwYIDuq8zNhX79miWmuztIJCAUgpub7qvvmujaKCGsWbUKsrOVX3v36qVKOzuorGy2pbISOBywt9dL9V0T9UY7SiqVmpqast2Klq5fv15VVVVYWPjaa68ZYPMMxJQpwOE0vvTwgMWL2WsNu4YOhbIySE0FHx/llvPnYcAAsLZmtVlPNuqNPl5dXd28efPefffdiRMn/vbbb2w3RyklJeXVV1/19/dfsGDBokWLRo8efePGDbYbZaAQ4ccfQSRSviwpgbo6Vhukf6WlEBMDAODgAEuXwoIFcOsWVFXBmTOwdi2Eh7PdviccknbV19dPmzYNAMzNzZm/sQkTJvz+++8sNik9PX3u3LlGRkYA0KNHj5CQkIEDBwKAkZHRu+++W11dzWLbGMnJyXFxcWfPnmW7IUo7d6KxMa5Zo3z5ww84dSqrDUJERD4ff/ml8eXff+PGjTqrLCQEAXD7dkREmQx37sRRo7BfP/T3x+PHdVZrd0Ex2h6xWBwUFAQATk5OAoEgIiLCxcWFCdNRo0ZFR0crFAp9tufBgwehoaEmJiYAYGZmFhoaWlBQgIj19fXh4eFmZmYA4Orq+tNPP+mzVU3l5uaGhoYaGxt7eHgAwIwZM3Jzc9lqjMrOnTh7Njo54Z07iAYTo3/8gZs3o1yufHn9Oh45opuazp9HDgetrDAjQzcVdHcUo48kkUheeuklAHB0dExOTmY21tbWNg3TkSNH6idMS0pKuFyuhYUFAJiYmISEhGRlZbUoc+fOnfHjxzMNCwoKysvL03WrmhIKhcuWLWOuz5qZmQUGBvbu3RsAevbsuWfPHrkqLdiwcye+/TZ+8QX6+6NCYSgx+tFHyOHg/v3Kl/v34+zZOqhGJMIBAxAAd+5s3PjgAd69q4PKuimK0bbJZLI5c+YAgIODwx2mD9NEQ0MDn893dXVlMmvEiBG6C9Pq6moej9erVy8A4HA4wcHB6enpjyosl8v5fD5TuHfv3hEREXrIr7KyMi6Xa2lpyVxYCA4OzsjIQMSCgoKQkBDmr8jPz08gEOi6JY/CxKhUiiNH4qFDBhSjM2eioyMWFCDqLkY/+AABcMQIlEgaN86Ygebm+PPPOqivO6IYbYNMJps3bx4A2NjYtPPDz4Rp3759dRSmIpGIx+PZ2toyxw8MDIyPj+/IjkKhcPbs2cxe/v7+KSkp2mpSCzU1NTwez8bGhon4oKCgxMTEFmVOnTrl5ubGdKLDwsJqa2t11Jg2VVRgUpIyRhHx6lV0dsaDBw0lRteuxdWrcf58RB3FaGIimpqisTE2/W/8ww8IgL17o1Co7fq6KYrRlmQy2YIFC5je3M2bNx9bXiwWNw3T4cOHax6mzDFVlw78/f0vX77c2YNER0f36dMHAL76aurDh5sUCrEmTWqzhc7OzqqIv3Xr1qMKi0QiLpdrbGwMAF5eXvq59SQSYUQEOjmhtzfyeMoYRcTFi3HgQJw6FU+fxhdfxPx8PbSlbUyMVlWhqyvGxuogRmUyHD0aAXDVqsaNlZXo6ooA+N//arWybo1itBm5XL5w4UIA6NWr140bN1q8m5qaOnz48KioKJlM1uItJlb69evHxMqwYcOioqLUOKGWy+XR0dFeXl7MccaOHRsTE6P2x6moqNiw4b34eGuBAJKTfWpq/qf2oVQUCklx8YFZs55lWhgQEPDXX391ZMf4+PjRo0czewUHB5eUlGjemDaJxfjFF9inDwIgAD77LG7Y0BijpaXo4IDPPotDhiAA2tnh11+jfu8UYl0d5ucrYxQRjx7FIUNwzx4tx6ic/wUCoIcHNj0DWLwYAXDSJH1/5i6NYrQp+dat7zN3Ra5evdr67eXLlzMp4Ovre+zYsdYp2WaYts7cNikUipiYmBEjRjD7+vj4aOsSQU3NleTkoQIBCASc7OwQqbRM3SPJy8ujk5IGCQRw9Ohwpt/dqf2lUmlERIS1tTUA2NnZ8fl8dVvyiPbJMTpaeUMFAMeOReZ30OXLzS4DnjmDhw5hfj7OnKksGRCgpzsuEglGRaGnJ06b1hijiPjsszhiBM6ejadO4YYN2NCgaUVicXbCLZvqtZPxjyaD8y5dQg4Hzc0xNVXTCkgTFKMqitzc5bdvO7300rCLFy+2WUIikURFRTGDNJmke1TPNCoqasA/D0j7+vo+NkxjY2NVPTV3d3c+n9/B8O0ghUJSUMCLizMXCCAx0aW0NKqzB6ioOJGS4isQgEAAyclDy8uPI6oZ8RkZGYGBgcyHff7557Ozs9U7TrP2KTAmBkeMUMaijw9GR3eovxUdjU5OCIAWFhgejmJtXvloRi7HI0caI97PD1etaozR9HQ0N8dZs9DdHQFwyBDs/FWcZu7fDxIIICtrfuOmhgYcOhQBcOtWjQ5NWqEYZSgePHhHIID4eMuqqgvtF2XCdNCgQUwQeHl58fl8qVTaZjFVmD4qc69evTp58mSmTL9+/SIiIho074o8QkPDvfT0KUwU3r8fJBY/6MheNTVX0tICmL3u3HErKeErFC0/rBqio6MdHBwAwMrKisfjafJrIzZWeQ0QAN3dkc/HTh2sogJDQ5HDQQAcPhyvX1e7Ie21cNQoZQuHDMGoKJTL8ejRZgNF9+7FXbvw6lX09UUA5HAwJARLS9Wprqzse4EAEhLsJJJC1UZF+HoEwGHDdPi7oruiGEVEzM9fKxBAXJxZZeXpDu7SqTBtswN7586d4OBgZru9vT2Px6urq9PyB2uDorQ0KiHBTiCA27d7FRVFKBRM5Chqa6+Vl0fX1t5EVF6sqK29np4+lQnQhATHggKeXF6vxaYUFRWpRkSNGjWqIzf0Wrh06dIbb5xk4qlfP+Tzm43q6ZTLl5VXSzkcDA3Fqio1j9PClSsYEKAM0P79kc/HVv9HWpJIkMdDCwsEQGdnjOrkmYNUWpaQ4CQQQGnpN6qNdXVJqRdcxa9MxGvXOv0ZyONQjGJ+/rp/MvS3zu7L3BEaPHgwkwWenp6PCtNDhw55enoyxQYNGvTMM89wOBxmPMDmzZtramq09Gk6RCLJz8iYyeRjWpq/WJx3715gaqpfTk5oevqU8vIf6+qSMzODBQIO06kRCsNlMl09Y3rmzBnmkSdmRFQH/ypUv4SMjc1Hj677/HOs72TCy+UtE62uDtetQ1NTZeSdOaNR4ty4gVOnKgPU0RF5vM618P79xt1feAFzcjq6Y3b2GwIBpKf/q8lVF3la2kSBAB48eKcTLSAd1t1jVCj8RCCAuDjTioqTah+ECVNvb28mJT08PNrpmTIdWAcHBysrq7CwsKKiIs0+gfoqK2Pu3OmfkjKsuvpCQoKTQtHYkUtLe0YggPj4HkLhepmsUtctaToiytPT848//mincGpqanBwMPNLqEePHlwut7JSnRZ+8QWOHImte8B37uD48Thhwo+g7sNgCQkJK1ZcZxLQzg55PBSJ1GggKhQYGYk2NgiAPXvi2bP3/jl1eKTq6osCAScuzry+vvGWWVHRl8w1cZmsQp12kMfp1jEqFIYLBBAXZ1xefkzzo7UZppJWJ5lSqdTKygoAUg3gbqlMVllfn15TczUuzkQkahz4WVNzOS9vpVRarM/G3L59++mnn1aNiCoublm76oF95oqqJr+E5HLl/SgTE/zgg5YxJ5Ph/v1f9+zZEwBsbGwiIyM7OGQiOzubaaG1tbObm4zLxQqNg6uwEENC0NW1Mj6+b2qqn0j0yOdBFIqG5GRvgQAKCrarNkokwtu3ewsEUFHxy6N2JBrqvjFaWPgZk6FlZUe1eFiZTHbkyBFVmHp5ef33v/9tGqYPHz5keqNarFRjiry8NfHxFmlpAVr5jaI2ZkRUjx49AMDW1pbP5zP5VVxczOVymUm2TE1NQ0NDhRo/gVNXh1wumpggAHp6YusesFAofOWVV1RPQLT/MFhubu7ixYuZWWMsLCxWrlxZXFyuYQubys6Ov3OnP3PmlJ+/Ti5v4zK6TFaRlfVaSsqIpicWzNWbzExdPGdKlLppjBYWfv5PhupkUh2mZzp06FAAcHFxaXrv6NKlSwAwYcIEXdSrCZmsvLj4YEKCQ3HxXnZbkpWVNX36dFV+vfnmm8wD+8bGxq+//rpWBkipJCTgmDHKS5DBwdi6dxsTE8MMBDY1NeVyua3HUXRk1hitkMtF+fncuDhjgQCSkryqqtp+GEwub7y4XF5+nLmXKJGw97RWN9BdYrShISsvb01m5py8vDUiUVxSkpdAYNT54ZOdI5fLjx49euxYs/5dZGQkALz++us6rVptBQXb799/nu1WICJGR0c7Ojoy+fWoB/a1QirFzz5DKysEQFfX+u+++7HFKXxlZWVYWBgzweugQYMuXFAOievUrDHaIhLFp6aOZm4PZmYGS6XNHgZTKBpTXiarSkzsKxBAcfFBXbeqm+sWMSqTVSYm9hEKP66q+rO4eF9t7d9icR5bZ69r1qwBgK2GNARaLM6prb2GqJDLa+7d+3de3qrH76MXJSUlkZGRkZGRcXFxuq4rKwunT8fJkz8BgEmTJqWlpbUocOXKFeb0gsPhzJs3b8OGDWrMGqMVCoW0qCiCecY3IcGupISPKM/NXX77tm1iYp+kpEENDRmI2NCQkZr6VFraM6oRbERHukWMVldfTEiwZ7sVSi+//DIAdPYxSp2qrb2elhaQkOCYkOCQlTVPD7fmDdZ33x1lusCWlpY8Hq/FcIuGhob169ebmZlxOBxmqIB6s8ZoRUPD/fT0ZwUCSEiwLyk5lJw8lLkRLxbnqa6NKhRSqZS1oSDdBwcRoauTSB4kJw92cnrX2fl9U9O+7DbG19c3NTU1ISFh5MiR7LaEtKmiouLDDz88dOgQIo4YMSIyMnLcuHFNC6SkpJw6daqmpmbKlClTp05lq50AAIBlZYeNjKzE4uzKyl+HDLkOwHn8TkTbukWMAkBt7dXCwh3V1ed6957h4fGVsbENK81QKBTW1tZisbimpsaa1mI0YJcvX16yZEl6erqRkdFbb721a9cuZvyTYRKLM9LTA8zMPBwdQ21tXzUyov9aetVdYpQhkxVnZgZbWAx0d/+KlQZkZWUNGDCgf//+Dx48YKUBpOPq6+s//fTTHTt2WFhYpKamqqaUNUwyWXlFxbGyssNSaYm3919mZv3ZblE30r3WqTcxcbK1nV1VdZqtBty7dw8AVA+PEkNmaWm5cePGWbNmZWdnG3iGAoCJiZ2j43JHx2Xp6ZPLyr6xth4HwOnV6zl2W3X69OnMzMyampo333xTNQ1519MtYrS29mpl5c89ekxWKETFxRFOTqvYagnF6BNn5MiRhnoVG1VXQhUK0T8n8hwOxxjAKDf3LYkk39Y22M1tn4mJo/4bd+3atY8++ujy5cumpqZSqXTXrl2ffvrp22+/zdya62KM2G6APlhYDDIz86is/LWm5nzfvjwnpxUA7FzKKC9/ZvLkraNGvcRK7aTLQBSnpU0sLt4DoACAoqLPk5MHZ2XNSUsbJ5MVOzi87ez8gZFRj4qK48nJ3qWlkfr8D5+UlPTqq68yYxjs7e3XrFkzc+bMysrKJUuWBAQEpKam6q0l+sPuQAH9UyhkhYW7kpN9mj7soTfPPYcAeOqU/msmXQozpU5y8hDVeHuxOKe6+lxt7d+q6UsaGrLu3ZvODNRPT/9XfX3LkbBa186sMY99GOyJ1u1iFBHv3p0gEEBBwTb9V+3hgQB4757+ayZdR3393bg4c4HAqCOLa5WXRyckODJTkguF4dpd2VAlJyfnsbPGPOphsC6gO8ZoTc1lgQBu37bRYFUidTQ0oLExmpioP7UwIYjytDR/gQByc5erNslkVffuPfeoVJXJynNyQpnZY1NSRtTW/q3F1hQVFXVq1pgrV674+PgAAIfDCQkJKSvT68+gjnTHGEXEe/cCBQLIz/9In5UmJyMADh6szzpJV1NcvK/15KHMEjh377Y3301V1dmkJC9mRp4HD96Xy2vbKdwRpaWlXC6XmTXGyMgoODg4IyOjIztKJBIej8ckb58+faI6O7+/4emmMSoS3RIIOPHx1hJJgd4qPXECATAoSG8Vkq5GInl4+7aNQAAVFSdUG2tr/46LM46LMxGJbre/u1xeJxSGx8WZXbpkN2KE74kTJ9ov/yg1NTU8Hq93795Mp/Kxs8bk5+e3XiHm3r17U6ZMYe7QBAUFPXjQoZXBDFM3jVFEzMh4WSCABw/C9FYjj4cAuMpQ5v0gT56MjFnMcoSqLQqFNDV1lEAAQuHHHTyISHT7gw/mMvk1d+7cwsLCx+/zD2YJcWdnZ2b3wMBAgeCR00irzJw508jIKDQ0tLq62VI0CoUiKirKzs4OAKytrTVc2ZBF3TdG6+qSBAKjuDizhgadzA7Z2uLFCIAHDuinNtLVKE7/WsAf22Ly0IcPtwgEkJQ0qFOrDcrlcj6fr5rePyIiQi5/zCxQEomEz+erHkOYMGFCB+8RyWQyLpfLTGjt5uZ2qtU4lYKCAtXKhhMnTkxOTu74BzEQ3TdGETEr6zWBAHJyFuunOmaFyPPn9VMb6VoqK7FvXwSQnPxata2h4V58vIVAwKmuPqfGIfPz82fNmsXkV0BAwN27d9ssxsxBrlrddvjw4WrMT5aYmKia4SUoKCg/v+U00jExMf3791eNiKrv7AqFrOrWMdrQkBUXZxYXZ1xfr49lkfr0QQB8ki8BEfYsW4YAOG4cNnYbFczy1xr2A2JiYpg+poWFRXh4uLj5KvaxsbGqh7iGDBkSHR3dwWWpWntsF7iqqiosLIwZNTVgwIBz59T53cCKbh2jiJiTEyoQQGbmHF1XVF2NAGhpiY87eSKklRs30MgIzcywyQmv9OSR+795JSY6az5ur6KiIjQ0lBk2P2zYsGvXriFibGzsmDFjmAB1c3Nrc7FbNQiFwqZd4NYLO169etXX11c1Iqq0tFTzSnWtu8eoWJx37ZrTqlWTb99+zF3OjnjwAEtK2n7r1i0EwJEjNa+EdDNiMfr6IgBu2NC4sbAQbW3R3Lw+TmvrfV64cIFZ/dvIyMjT05NJOhcXl3379rXoompO1QVu86EmZkQUs7yVs7Oz4Y+I6u4xiojr1q0FgBkzZnR2x/JyvHIF+XwMC8PAQHR0RAD87LO2C1dU4M8/4y+0xi3prE2blOONm14unDMHAfB5La+aVV9fHx4ebmJi0rdvXxsbGx6PJ2qx9rT2VFRUqB5qUnWBm7p//75qVuwXXnghJydHRy3RHMUolpSUMKuStb8aRHFx8fnz57/88sslS5bMnbvDzk65nGTTL0dH5PGa7ZWbi7a2+MUXypcyGdraYrk2V94lXVp6OlpYIIfT7NbkmTMIgFZWqJslSCdMmAAAMTExujh4C3/99ZdqhavQ0NCqqqqm7yoUisjISGbNqx49ehhst5RiFBFx/fr1zJWaFtujoqJWrFgxZcoUBweHpvO52NoOAEAbGxw9GkNCkMfDmBh8+LCNI2dno7k59umDzJ1JmQwBsEs8/0b04p13EADffrtxi0iEnp4IgLt366jOsWPHAsCNGzd0dPwWmC4w81CTi4vL8ePHWxQoLCxkRkStWLGiU6Nc9YZiFBGxsrLS3t4eAM6ebbb298SJE1XRaWNjM3HixNDQ0D179pw7d66goEO3irKz0cYGN2zA//s/RIpR0llyOR44gBWNz33ie+8hAI4Zgzobqc7c4UlKStLR8duUlJTE9IKZEVF5eXktCnh7ewOAHtaIVQPFqBKPxwOAp59+uul4jm+//fY///nPn3/+2foftSPy85UxKhKhmxuePk0xSjRz86Zyehtdrufs4eEBAFm6uWLQDrlc/uWXXzIjop5++ukW7zKTnbde+NoQUIwq1dXVMbcO1X7QuMUdJycnBMC//0YbG0TEEyfQ0xNraylGScf8/jvOn4+BgRgSgpcuKTdOn44AyOXqtGZmiWnVTHcZGRmWlpZ+fn46rVQlLy/v5ZdfPt/qMRXmx1O9Do2uUYw22rt3LwD4+vp25MHe0tL6S5dw/35ctgz/9S+0t2/jjpODA8bEKGMUEV94ATdvphglHXDwIDo74zff4PXreOAA2tri0aOIiBUVuG4d1tXptHJmzdqaGuW85omJiczDSzqt9LGYG03lBnl/tnutDNo+qVQ6ZMiQrKysI0eOLFiwoOlbVVVVGRkZKSkpqampzJ/W1s8lJR1oWqZ3bxg4EHx8wNdX+aeXF+TkgJ8fVFQAAGRmwoQJUFICZWVgZ6fPT0aeKGIx9OkDP/wA//63csvRo7B2LeTlge4XMkJEExMTRJTJZMxopOvXr0+cOHH8+PHXr1/Xde3tMDc3l0gkDQ0NzM0og9ItlrTrIFNT0/Xr1y9atOjjjz/28PBIT09PTU1NTk5OSUkRCoUtCvft+9f48TBsGPj4wLBh4OsLrq6POf6AAbBiBWzcqKPmk64iKwvq62H69MYtM2fCa69BURH06aPryuvr6xUKhaWlJZOhACASiQCA6aKyRSaTSSQSY2NjA8xQ0DBG//wTFAp4/nnly+xsiImB995rLFBYCBER8NZbwExrUF8PmzbBtm1gbKxJtToUEhKyefPmgoKCgICAptvNzMwGDhzo6+vr4+PD/Dl06FCjDqwH6OwMv/zS+JLLBX9/6NlT2+0mXUlxMdjaNut4WlmBlRUUFuohRluHZl1dHbAdo4YQ5e3QKEZjY0Emaxaj+/c3i9HiYti5EwQCOHcOAKChAT79FLZsMdwYBQBEVCgU3t7eY8aM8f2Hp6enegvDWlrC5MmNLy0sIDBQWy0lXVT//lBaClIpmJoqt1RUQF0duLnpofLWgcVssbKy0kPtj9KVY7Qj+vWD6mr47jtofrHRQP3www85OTkDBgxITk5mZkgkRN/c3aFPHzh2DP6ZhROOHIFhw/RzQb1139MQIswQesTt0DQpZDKoq1N+Lxa3UYDDgYgImD0bZszQsCqdk8vl27ZtA4D169dThhLWGBvDwYMwfz5kZICvL8TFAZ8Pv/2mn8pb9z0NIUYNoUfcjg5c3mtXVBR4eyu/Xn+97TITJ8K//w3r1mlYlc79+OOPaWlpAwYMeO2119huC+neZsyA69dBJoNz58DcHOLjofnFet2ha6Nq0LTPtXgxREQov79wAZYta7vYrl0wdCi88oqGtemQQqHYsWMHAHz88cfUFSXs8/GBbdv0Xy1dG1WDpr3RDrK3h82b4cMP9VObOn788cfk5GR3d3fqipLu7FExStdG26H9GM3IgCVLYMkSyMxstj00FMzMtF6bdigUiu3btwPA+vXrzQy2lYToHhNYTfuehhBhhtAjbodGZ6+jRoFc3vjS2RleeAFu3YJNm0AkgrVr4fPPG+8sGRnB/v3A44FUCiYmengcoxOOHz/OdEVVKxQS0j3RSb0aNIrRFmOYfH1h927l95mZYGYG7u6wf39jgaeegtWrYfRoWL0a3npLk5q1CRGZG/SffPIJdUVJN2eYJ/WG0IZ26OTaqFQKa9e2fWs+OxvS0mD1anjwQBc1q+Onn35KSkpyc3NbuHAh220hhGWGGaOGcGGhHdqP0eJiCAmBlSthxIg23p07F2bPhupqWLwYDGFSFEQ4cGDCpEkr1q/fSF1RQujaqBq0H6O7dwMiHDkCn3zSdoEDB8DJCc6dg2++0XrlnXbiBFy82C87e+/ChYvYbgsh7KNro2rQ/gDJHTseU8DREXbvhtdeg5UrYdo06N9f603oKETYsgUA4KOPDHcUASH6RCf1atDTuNEW5s+HV16B6mp48002T+1/+QUSE6F/f1hEPVFCAMBQY9QQesTtYCdGAeDgQXB0hHPn4PBhdhqACNu3AwB8+CEY5ByGhLCgdb/PEHqChhDl7WAtRh0d4fPPAQBWroS8PBYacPIkxMWBiwu8+SYLtRNimFr3+y5cuHD+/PmerM6SSzH6SAsWwKxZUFUFS5eyUDvzvPLHH4OFBQu1E2KYWgeWn5/fs88+a8zqJMEUo+3Ztw/s7ODmTUV09A191hsTAwIBdUUJaYn1wDp16tSDVqPKWw/DMigsx6iLCxw4UITou3TpCw8fPtRbvVu3AgCsWweWlnqrk5AnABOjMplM/1Vfu3Zt8uTJL7744ubNm9tsFfVGH+nVV50DAoZUVFSEhobqp8ZTp+DWLejTx4AeSCXEQOzevdva2jowMPDbb7/VW6W3bt167rnn/P39L1++7OjoOHLkyBYFDDxGDWKd+ocPHzKLUH/77bd6qO6DDxAAIyL0UBUhT5js7OxJkyYx4RAUFJSbm6vT6lJTU4ODg5mFznr06MHlcquqqlqUqa+vZwK0tLRUp41Rm0HEKCJ+8803AGBjY2KnEvoAAAoZSURBVJOfn6/1g+fl4YULqFAoX5aU4OHDWFCAly41KxYbi7W1Wq+ckCeMQqGIioqyt7cHACsrKx6PJ5PJtF5LTk5OaGgoc+fKysoqLCysqKioRRmpVBoVFeXp6WliYmJhYTF8+PAbN25ovSWaM5QYRcSXX36Z+QWo9SPv24cAeOSI8mVMDI4bh1euoLt7s2I9emBqqtYrJ+SJVFhYqJo30s/P79atW9o6clFREZfLZVacNzU1DQ0NFQqFLcrI5fLvv/9+ILMyO4C3t3f//v0BwNjY+P3336+pqdFWY7TCgGJUKBQyp/bfffeddo+8bx+OGoUuLlhejkgxSkiHnT592t3dHQBMTEzCwsJqNTtfKy0t5XK5lpaWAGBkZBQcHJyRkdG6WGxsrJ+fHxOgnp6efD5fJpPV1dWFh4cz8we5urqeOHFCk5ZolwHFKCJ+9dVXAGBvb19QUKDFw+7bh/Pn49KluGQJIsUoIZ0hEom4XC5z9u3l5XX27Fk1DlJTU8Pj8Xr37g0AHA4nKCgoMTGxdbErV66orsz279+fz+dLpdKmBRITE8eNG6e6dKuLa4BqMKwYRcTnn38eAF588cX2i1VVVeXm5iYl1V2+jCdPYlQU7tmDmzbhypW4aBHOmoVTpqCfH3p64oQJyhgtKUEHB7x2rTFGzcxw1KjGLyMjilFC2hYfHz969Ggmv4KDg0tKSjq4o1gs5vP5zs7OzL6BgYECgaB1sRs3bgQGBjJlHBwceDxefX19mweUy+V8Pp95qsrGxiYiIkIul6v/wbTB4GJUdWr/xhtvbNmyZfXq1YsXL549e/bUqVNHjx7t5eVlZ2dnZKQcp+Xvnw2A7X/166eMUUSMjEQ/P/z1V2WM9uuHubmNX9bWFKOEPJJUKo2IiGBumtvZ2fH5/MfucvjwYVdXV+anddKkSVeuXGldJiUlRXWzvmfPnm3erG9NKBTOmjWLOXJAQMDdu3fV+UhaYnAxiogbN25kOv/t6NWrl5ub29y5CZMm4Usv4cKFGBaGGzbg55/j11/jiRN44QLGx2NWFlZUNMaoXI7jx+PLL9NJPSFqyszMnDZtGvNjOHny5PT09HYKM8vzDB8+PDo6uvW72dnZqpv11tbWXC63nLl90WFHjx5l+rkvThn7148Rcqmkcx9GSzhoCHPQNxccHPzTTz+NGjUqKCjIpjlbW1vmG1WHtCP274erV+H77wEAbt+GMWPg6adh1y5YsABychqL9ewJN2/C0KFa/jiEdD3Hjx9fvnx5aWmppaVleHj46tWrTUzamLxYJBKdPXt25syZnOZrWObn53/22WcHDx6USCRmZmZvvPHGxo0bXVxc1GhJeXn52rVrhhmlVxVkOXkMffGdnf28R6v5qdTGSni3Izk52cjIyMLCQosXj48cwfffb3zJ5eKsWXjzJo4f36yYtzfev6+tOgnp4srKylRPHo4cOfLmzZsd2aukpKTFzfrMzEzNG5ObfGPv0oCNQS4bX3T9be+aBlG15sfsOIOL0eDgYAB477332G4IIeTxLl68OGjQIPhnRFQ7IzqZm/W9evUCAA6HExwc3P4Fgc6Sihsufr9ry0y3jUEuuxaOSr16SosHb59hndSnpqYOHz7czMwsIyOjb9++Ld6VSCSVTVRVVTU09CksDKishNZfH30EK1aw8iEI6V7q6+s3bdq0a9cuuVzu6el54MCB6dOnNy1QV1d36NCh7du3FxcXA0BgYOCnn3761FNP6aIxxblpv325Jj89DgAGj5k2Yzmvl4M61wo6R2+B3RFz5swBgHfffbfpRpFI5Orq2uYcWU89xX3UDfoNG9j6EIR0RwkJCWPGjGF+MIODg5knOyUSCZ/PV1309Pf3v9TiEWwdUCgUgt+P7Hh18MYglx1zvG+cPKRQyEVVZQ/vJz68n1ian6H1O1EG1Bu9e/fusGHDTE1NMzIy+vXrp9qOiGZmZjKZzMTEpMUdJw+PZzmcZTY20PrL1pYmwSNEr2QyWURERHh4eF1dnb29/axZs2JjY3NzcwFg3LhxW7duVY0M1YPq0oIzB9el/30WANx8xg4eE3jtl4N9vHwl9XVVJfn/t/agm+84rVWm3VTWxNy5cwHgnXfeaf2WUCg0tKdoCSFtysrKYk7qmafgfXx8oqOjFaqZgfQr/e+zn7/x1M75vtd/4R/dtJDZeOX4l1EfB6vKyDTunGp/gWX13L9///jx4+bm5lwut/W7qhG8hBAD5+np+ccffxw9etTNzS03N3fevHmdGp6oXYPHTnPzHVeWn1GUc1e1US6TmJiaA0Dk+9PdfMclXTrh/38rJs5apnYthhKjGzdulMvlS5YsaXo6Twh5Qs2fPx8AnnnmGbYbAhbWvfp6P1WUczc/Pe7I+jniutrK4vxFvF8AQFxX08PG8YPvkloMa+0s9me/B4B798pjYn4zNzf/8MMP2W4LIaRrcug/KPD1j597c8PgMYEn96xEVACA+7AJGmYoGEhvdMsWO0vL+6tW/Y+5mEIIIVpnYdXLZeAIAHAZMHx78MDq0gJtHZn9GM3IgGPHgMNxXrRoNtttIYR0cYiK5CsnrXrZ9bRz1tYx2Y/RzZtBJoMlS8DDg+2mEEK6KGNTM+G921+8PUEmFds49Zv7yTdGxibmVj2NjI01PzjL40YzM2HIEOBwID0dPD1ZbAghhKiJ5VtMW7aATAZvvEEZSgh5UrHZG6WuKCGkC2CzN7ptG8hksHAhZSgh5AnGWm80NxcGDwa5HO7ehUGDWGkCIYRoAWu90S1bQCKBhQspQwkhTzZ2eqOqrmhqKgwerP/6CSFEa9gZNyqXw4wZ0Ls3ZSgh5ImnnZP68+dhz57Gl6WlsHhxswINDbBkCZw7p3zp5QXu7rBtm1YqJ4QQNmknRtPS4MKFxpc1NXDkSLMCUilERsJbb0FNjXLL4cNQW6uVygkhhE16vcU0bRps2KDPCgkhROe0dm20thbS0pTfC4Vtl9m8GYYPh4ULwc9PW9USQgjLtBajCQmN10PF4rbLODnBxo2wdClcv66tagkhhGVai9FnnoGTJ5XfZ2eDt3fbxZYtg2++gcOHtVUtIYSwTN8DnoyNYf9+mD0bpFI910wIITqhw1tMmZmwYAEsWAB//tls+7hx8MILIBLprmZCCNEf7fRG/fzA2rrxpY0NrF4NhYXw2WdgZwdTp8L588DlgmrJkx07wM4O7Oy0UjkhhLBJ5w+DVlTAnDlw9qxOKyGEENbodtyoTAZLl8LWrTqthBBC2KTDGC0shPnzISwMxo7VXSWEEMIyHZ7UHz6sHB/K4cDBgzqqhBBCWMbyknaEEPKkY3lJO0IIedJRjBJCiEYoRgkhRCMUo4QQohGKUUII0QjFKCGEaIRilBBCNEIxSgghGqEYJYQQjVCMEkKIRihGCSFEIxSjhBCiEYpRQgjRCMUoIYRohGKUEEI08v8+rYQO9IBcUQAAAeR6VFh0cmRraXRQS0wgcmRraXQgMjAyMS4wOS4zAAB4nHu/b+09BiDgZUAAKSCWBuIGRnaGDJAAIxMbgwaIwQIXYFewADNYEDJgmhnGZ4OoZGISgAgwwQXgDCQzwCLMMPuYGKGa2DigmjnghqCYzsyIxIDIKENoRmwqoAyYi7gZGIE+YWBk5mBiZmFgYeVgYmVjYGNnYOdg4OBk4ORi4OJm4ObhYOLhZeDlY+Dj52DiEwBiQQZBIQ4mIWEGYREOJhFRBhExBjFxDiY2CQ4mCUkGSRYGcUEGESZWFlY2CUk2IWERMXFBcSdGoJPg4XvTW+6gwYLgA2DnzBc+6Ov1ez+IrfBM4aDHZjMwWyns+4EnLu/B7Io/Dw8smssGZuuu2HTA1qt3H4j9atHqA7fnzLIHsc0lGg+87PgMZs9WMTxQvKAHzM6W27G/oecpmK1tZrBP1KMWzGZTmWIvMeE6mH1VVtrB1T8RzLY2LnDYMPkkmN0+f4FDSZgrmL0kcpdDQSKTA4id+bTNYc3Gv2A3FKiedChazr4XxF6ucNxhv+ARsDsnLHvvIFUrCvbjQgleR3W+S2DxadGKjn7VMmBxNmk+R7mzgmBzlDd8dlgSqwm2K2/BtQNtzB/A7GrzHwe+bPcFs8UAEIdyZOirCYwAAAJ/elRYdE1PTCByZGtpdCAyMDIxLjA5LjMAAHicfVXLbhsxDLz7K/QDFkiKz0MPcZwGRRsbaNL8Q69F/x8l13BXAdSuLGJXmOVSnBn50Or6fv7683f7e9H5cGgN/vOLiPY+AODw0uqmnZ6ev1za49vD6b7yeP1xeXttpI0s38nxEfvwdn25r2C7tGN0jyFo7Tj6AFWhBh22a3+X2mMhiXQQtiP2yBuQBXJUTsx1UfTRjtANh9pYQLmg1o2ClCupa1DAhGzvD98+8R0vVYR1UAKNyiwkzqvMWpmliw/Ip0SSSYgtkFY5pQ/NYQ07gQWuOuAF5A6D8+sJDCBFXQCjvbYjdVNziQRCuPIKiFApsfNwJk+kxRgCKyQWEjoagVdxkZQ5rpCUO8eOYoPr4+oMLv/sKI5MTJ1paFSJ7oDLTWFxNbrJQKo2ibHTslTJ7UsHDN44N8vvL4Hargl0GpwUUjbWMFYMpTKvjTuhCHB1gUM4lkUWR9rJpWpLJIyBGCtk3JBENrQ4yH1ZLIVfJFlnj+p88poUma5aT0WSd0eXbGvmVEZYao6onX618oizK1RWZkk9rbDjltVRwjcNJFnOKyRvlSqmK6o+zTtfOZQ2H2k3SG44yXQThGWhWlpOh4raxqbDEPcF8uly/nCy3M6a0/Vy3s+aGrQfKPnQxn5qUI7cwefnT7SfD5hTdvtTTt09jjltNzLm9N2umDN2U2LN2Xu4BZw8hlugyUtUoUwyl5XPuAWe3IFbkMkGVAF10vttxSZh47bik4CpAsYkVKxAsyBpCzgJDysQTQK7rYxJRrQFnuRC2TCSSRZYgXSiHyvskKo4X9wrLt5nluv5/m+U94c/IexJBpMecS4AAAFSelRYdFNNSUxFUyByZGtpdCAyMDIxLjA5LjMAAHicJZE7boZADISvkhKkxfLbXiGabJOKFH8Z5RKpc/jYpEDAx+x4xtxru/fr6/74XnTda1uvtfrtWvdruz73vhZfa11re//Z68b7tV709rsdE3KK8TgEBN1tnIWYXQoRzLpTIUJA80IIQeIhxQKCJ3vL0ifPB6EzRsuMLbVlBpaC0ozDpsXDxMV1EJR69gAFFK1RNRK5Bp0HQ3ikF8GZTy4ClfYkiCnSBIGCMXkgzKqR5URAFvWRwFMxx8mgLFGKTCTycQqECbWNhdbZ0wBpaksiMK1BsqgOrlBRAgUms26g03SWhwOn5bMPFCHKeBhzaO+j/KP2EaA5K0N1qbxR5xKS0qg1roRFerOpszWqVntqTVbOLlfJU9vHqfZZAb0eOuHhEFhhq0SGEf3/DvN4iqJY5th//wC7vW2W8MDwfwAAAABJRU5ErkJggg==\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "tf2 = rdChemReactions.ReactionFromSmarts('[#6:1]-[Nv4H1:2]=[C:3]>>[#6:1]-[NH0:2]=[C:3]')\nwhile tf2.RunReactantInPlace(m2):\n pass\nChem.SanitizeMol(m2)\nm2",
"execution_count": 22,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 22,
"data": {
"text/plain": "<rdkit.Chem.rdchem.Mol at 0x7f74e0093a00>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxU5f4H8O+wL6KIICACKpKikOKSJnm1wJ+lpJliXosys7F7M24/lzteTbHcRm8l92rmaP4KM60xK9HKosTMjTzsDILssovsDDAM53x/f5xpQDaBWZXv++WrF3PmmfN8MfnwnHOe8xwBIgIhhJC+MjF0AYQQ8mCjGCWEEI1QjBJCiEYoRgkhRCMUo4QQohGKUUII0QjFKCGEaIRilBBCNEIxSgghGqEYJYQQjVCMEkKIRihGCSFEIxSjhBCiEYpRQgjRCMUoIYRohGKUEEI0YmboAgjpp774AoYMgaefVr1MSICYGFi7Vi99374NUVFQWgojR8LixWBvr5deH1o0GiXEMH75BW7caH2ZlQVff62Xjs+dAz8/SE8HR0c4fx7GjoW0NL10/NCi0Sgh/YlCAa+9BhIJLFsGAPD227BuHaxZAxcuGLqyBxjFKCEGU1cHxcWqr6uq9NKlTAYKBSxd2rrljTdg/HhoaQEzSoM+or84Qgzm2DE4f171dU0NuLnpvsuCAhg2DEzanM1zdwelEoqLwcND990/nOjcKCEG8+abkJys+vP++3rp0tERKivv2VJRAQIBODrqpfuHE41Ge6qpqcnS0lIgEBi6kHv8/PPPcrm8sLBw1apV1tbWhi6HGD0fH6iuhqQkmDBBteWnn2DMGLCxMWhZDzYajd6fXC5fuHDh66+//uijjx47doxlWUNXBABw/fr1wMDAuXPnrlixIiwszNfXNzo62tBFGamoKHjxRaitVb28ehXee8+gBelfSQmcOAEA4OAA//u/sHw5xMRAYSF89RVs3Ajbtxu6vgcckm7V1dUFBAQAgJWVFf835ufnJ5VKWZY1VEkJCQnz58/nixkyZMiaNWv8/PwAQCAQrFix4u7du4YqTC02NvbKlSunT5/mOM7QtSAi7t2L5ub41luqlydPYmCgQQtCRMTISPz++9aXcXG4d6/OOluyBAHw3XcRETkODx3CWbNwzBh85hn88Ued9dpfUIx2Ry6XP/nkkwDg7u6enp4eGRnp5eXF59f48eMjIyNbWlr0WU96enpoaKiJiQkA2NraikSiqqoqRFQqlREREba2tgDg4OAgkUj0WVVbaWlpISEhAoFg2LBhADBz5sybN28aqhi1vXtx6VJ0ccEbNxCNJka//RbXrkWlUvXy4kU8dEg3PX3/PQKgnR3evq2bDvo7itEuNTQ0BAYGAsDw4cOzsrL4jc3NzZGRkaNHj+bDdNy4cfoJ09u3bwuFQjMzMwCwsLAQCoUlJSXt2mRnZ8+ZM4cvbPbs2bdu3dJ1VW1lZWW99NJL6oh/7rnnhg4dyo/id+zY0dzcrM9i2tm7F19/HSUSnDIFW1qMJUY3bUITE/zwQ9XLgwdx8WIddFNbi+7uCID//W/rxuxsjIvTQWf9FMVo5xQKBX/g7OzsnJaW1u7ddmE6atQoiUSiVI8rtKq8vFwkEvGnFMzNzUNDQ3NycrppL5VKHR0dAcDa2losFush4svKykQikaWlJV+hUCgsLi5GxKqqKqFQyF+U8/X1vXbtmq4r6QofoyyLjz2G+/cbUYyGhKCjIxYUIOouRtesQQB87DFU/0vgOAwMRDMz/PJLHfTXH1GMdkKhUDz77LMAMHToUJlM1lUzPky9vb11FKa1tbVisXjgwIEAYGJiEhIS0sMBZkVFhVAo5KuaMGHCH3/8oa2S2rl7965IJOJnCPAVqoftapcuXRo7dix/6lYoFNbU1OiomE6VluIff6hiFBEZBh0d8cABY4nRf/4TN23CRYsQdRSjsbFoaopmZpiQ0Lrx008RAB0d8c4dbffXT1GMttfc3Lxw4UIAcHJySklJafduxytLfJg+8sgjfGyNHDlS8zCtr68Xi8WDBw/m9xkUFJTQ9segZ3744YcRI0YAwPHjT92+HcaydZqU1E5dXZ1YLB40aBCfj8HBwUlJSV01bmxsDA8Pt7CwAABXV9evv/5ai5V0pbYWxWIcOBBHj8bdu1Uxioh//zt6emJgIH71Fc6bh3l5eqilc3yM1tejhweePauDGFUqceJEBMBNm1o3lpejkxMC4PHjWu2sX6MYvUdLS8uyZcsAwN7enmGYdu+mpKSMHDkyIiKiqamp3Vssy0ql0jFjxvDBN2LEiL6FqUKhkEgkLi4u/H4CAgIuXbrU52+nvr5+797N8fG2DAMpKSNran7q867UWLahtHTvc8/N4iucN29efHx8Tz6YnJw8ffp0/lPBwcEF/KGsDtTV4Y4daG+PACgQ4MKFGB7eGqM1NejqioGB6OODADhgAP7nP6jfK4VYU4PZ2aoYRcRvvsHRozEiQssxyh76AAHQ2xsbG1u3Ll+OADh3rjZ76vcoRltxXMt7770BAIMHD47r7AT8mjVr1EPOI0eOdLxs0mmY9vDqCv/ZUaNG8Z+dNm3aL7/8ooXvCrGhIfHmzakMAwwD2dkhSmVZ3/bDcc3l5ZKkpGEMA1Kpz6xZsy5fvtzLPXASiYQ/TTFo0KCIiAjtzhtTKFAiQVdXBEAADAhA/ndQbOw9s3ouXsSTJ7G0FENDVS39/bHDL02d4Ct0dsYnnmiNUUScNw/HjsXFi/HUKXz7bazT+MihsTEjgbGt2jYLYy60bj1/HgHQxgazszXtgLRBMarG5uS8FB9v/+KL/l2dTGRZNioqauLEiXzSeXh4RERENLb9Vf9nM6lUyp8QBABPT89OB7BqHMdJpVL1aYHx48dLpVLtzrjkOGVZWUR8/ACGgYSEweXlEsRe7Z+trJSmpIzms1gm86uslPa5mOLi4sWLF6uH292cfe5FfSxKpThqlCoWp03DHv4O+v579PREADQzw7AwrK/XvJbONTfjxx+jm5uqwpkzcd261hjNykIrK1y0SHVR3cMDz53TpDcuIyOQYSAvb2XrNrlc9Rf0/vua7Jp0RDHK4/LyhAwD8fG2dXX3OYjmw9Tf378nYerj49N9mEZHR0+aNKnt6FV319abmnJu3ZrLR2F6+l8aG9N78qmammiZbAL/qdRUn8pKaS8juHNRUVHDhw/nL+6LRKJufs10j+NQKsVHHlHF0/jxKJVir34HyeUoEqGpKQLgqFH48899K6RLfMSPHq2q0M8PpVJExG+/xdOnW5t99hkeOoRJSThtmqplcDAWFvalx/LyTxgGEhMdlcpy9UZu80YEwEmTUDdTSvozilFE5PLz/8YwEB9vU1t7sYef4cNUnYDu7u5dhWlXmXv58uVZs1RnGIcPH667KVPtVFZKExOdGAbi462LisI5ToGIHKesqTlfXn6kpiaa34KINTXR6rMBycke5eUSjtNmxFdXV4eFhfFTTb29vS9cuHD/z9zr7NmzL7zwJR863t544gT2+SRBfDxOnqzKr5AQLC+//0d6IjpadZkHAMeO7VHEsyxKJGhnhwBob48REb37ppTK8sRER4aBioqT6o1yeVzKJcfG0Fk0XVQXKEa5/Pw3/8zQXv8YcxzXNkydnZ3FYnFDQ0O7ZizLfvnll76+vurQVJ8ZGDp0aPeH/LqgVJbn5r785xH6BIUi7+bNx9PTnygoWJeVtaCq6nR9/bWMjKf4BomJTiUlYo7TVYW///77uHHj+Cv+oaGhFRUVPflUTEzMjBkzAMDExOyJJ2qPHOn1GEupxHZ/683NuGsXWlsjADo7c1KpRlfkoqNx6lRVgHp4oETSuwqLinDRotYzAB3mLncpJ2cZw8CtW0+rt3BcS1raZIaBgoL1vfkOSE/19xgtKNjAMBAXZ1ld/UOfd8KH6eTJk9XJ2GmYts1cJyenAQMGiEQiPc+jbKu29mJq6hiZbEJNza9JSS5tR5p8hiYmOpSUiFlWrutKmpubxWIxP4HfxcUlMjKym8Z//PFHcHAw/1ft6OjY6V91T+zdi6NH46+/tt+enY1z5uATT3wGAM8880xubm5v93z58uXXXrvEJ6CrKx44gApFHwpERDx5Ep2dEQCtrPCHH5I47j6XK6urf+DHBE1NrReRSkv/zTCQnOyp3UlvRK1fx2hh4UaGgbg4i+pqjc7n8/iUnDJlSvdhyrKsjY0NAHR/M5J+sGyDQpFXXx8bF2daWxuj3l5ff62oaGtLS7U+i7l16xa/ggE/I+p2hxvA1TfsA4CGv4RYVnUILxDg6tVYVXXPuxyHkZEnHRwcAMDW1vaDDz7o4QnrlJSUkJAQALC1HfrII8rwcKyt7VuBraqqMCwMR426Gx8/VCbzra+/2lVLlpWnpIxkGCgr26feqFDk8ZcWq6vPaloK6UL/jdHCwk1/ZmiUFnfLcdx3332nPsx3dXXdt29f22P2wsJCPmS12KnmiovfS0iwS0ubWF5+xIBlcBwXGRmpzi/1zax5eXlCodDU1BQAbGxswsLCysr6OG1LrbkZxWK0tEQAdHHBjiPgXt0Mlp6evnTpUj7i7ezstm7dWlOjcYK2UVycnJLizTDAMCa3b69paenk90dLS21+/hs3bz7W9sAiM3M+w0BOzktaLIa0009jtKhoC8NAXJx5VdV3OuoiOjr6scce4y8rKdoc1F24cAEAnnjiCR3122csK6+o+CIpybW01MATYoqKihYtWsTn19SpU5cuXWpubg4AFhYWb731Vsc1WTSRmYlPPdV6ZTw/v30D9c1gZmZmYWFhdR2mdPZk1RitYNnGoqLwuDgLhoGkJNfKylOdNmt7Frui4nOGgcTEIX2eLEx6or/EaGOjLD9/dWbmvPz81fX111JTfeLizCordX5X4rlz586cOdN2i0QiAYCVK1d29RHDKikRZ2Y+Y+gqEBGjoqLc3d0FAoGVlVVXN+xrBb/85qBB/JUl+ccff9LupgC5XC4Sifix8MiRI8+fP89v7+2qMVrR0JBy8+bj/NW/zMxghaL1ZjCOa1YqW1ebVSrvJiYOZRi4e/czXVfVz/WLGFUqKxITHUtK9tTXX62o+Fwuj1cqywx1qmj9+vUAsHv3boP03qmmpqyamp9YtlGpLMvIeLKwUGToilRqampOnDhx/PhxPSxaWlKCoaE4e/Z6AJg0aVLH29gSEhKmTp3Kj5EXLVr0zjvv9GHVGC3hysslCQkDGQYSEgaVlUVwnCI3NzQx0SElxSspya2pKRMRFYrc9PSZGRmBWpnnS7rRL2K0tvZCYqKDkfxjWrBgAQDoZ3mOHpLLb2Rmzk9KcktO9sjPX82yOruVx+h9990Z9U0BmzdvbjcRmF8ee8CAAQKBgD8NGhQU1MMlBbROoSjIylrAMJCU5HLnziGZzI+fU6FU3m1zbpRraanqZidEKwSICA87pbIoNXXMkCEvOTuvs7T0NmwxPj4+6enpycnJ/JM/iLGRy+Xbt29///33WZb18vI6dOhQUFBQ2wY5OTk//vhjXV3dU089xZ/+NqCqqq8FAoumpptVVdKxY2MFAnpIpQH0ixgFgIaGhLKyvdXVUQMGBIwc+YWZmZNByuBnO7W0tNTX19ODPI3ZtWvXhEJhamqqQCB49dVXP/jgA3t7e0MX1aXm5ryMjFmmpoMcHVc5OLxkZuZg6Ir6l/4SozyWrc3JWWJu7jJixDGDFJCVleXt7e3p6ZmXl2eQAkjPKZXKDz/8cNu2bba2tmlpafwzUYwWx9VXVX1bUfFpU1PGmDG/W1qOMnRF/Uj/OgQwNR04aND8mpofDVXArVu3AEC9mBMxZvyaKYsWLcrNzTXyDAUAE5MBQ4aEDhkSeutWUEXFMRsbfwDW3v55A5bEL12Wm5tbU1OzevXqkSNHGrAYneoXMVpXd7Gy8gs7u1ksKy8t/fewYdsMVQnF6APnkUceMc7/X4isQGDKf82yVaamgwBMABBRIRBY3L69Wqkss7df6O5+wMJiuP7LO3v27JYtW5KSkszMzFpaWg4cOPDuu+/+4x//4CfYPmRMDF2APtjYTBwwIEAuZ5qa0j09jzg6rkJsMUglFRVzZs2K8PdfYpDeyUOD4xpu3pxYWrqH/5dcXn44JWVEVlawTOaL2Dx06N/d3PaamTlUV5+RycaWlu5BZPVW25UrV2bPnr1gwYKkpKThw4fv3bt31apVcrl8/fr1U6ZMuXHjht4q0R/DThTQP45TFhWFJyeP0PMN47ygIATAH/q+CgohiIgFBesZBmSyR9WLlSiV5fX11xob09QT+5qbS3JzQ/mJ+jdvPt7Q0P7BYloXGxvb1aoxMTEx/MMfu7oZ7IHW72IUETMynmQYKCoK13/XHh4IgLq5GYf0Fw0NSXFx5nFxpnL5/Z98Ul19NjnZnb/1ubBQxLLtl8TVCplMpl41xs7OrtNVYxoaGjq9Gewh0B9jtK7uMsNAfPwAPd9o3NCAJiZoYUGrj5O+a7N46Fr1xpaW6vT0gJqazoOJZesLC0VxcaYMAykpXjU10Vqsp92qMSKRqLKyspv2iYmJ6pvBQkJCNF9ixhj0xxhFxMzMeQwDBQXr9NlpcrJqCXRC+qy09H3+YQRtFw/Ny1vFMJCRMaubD9bVXZHJxjMMMIwgL++15ubuwq4nCgsLw8LC+FVizc3NhUJhcXFxTz6ovhkMAAYPHiyRSLT75DH966cx2tCQxDAm8fFWCkX7RS115+uvEQAXLNBbh+Rh02bx0NbVHevqfmMYQVycZWPjfZbI57jmkhJxfLzVr786jxnj3f3y2N24e/euSCTi7x/pyaoxOTk5MTExHTfOnTuXH5b+5S9/SU/v0ZPBjFM/jVFEzM5eyjCQn79abz3u2oUAuJ6e40D6KjMzmGEgJ2e5egvHNaWm+jAMFBdv7+FOGhvTN258hc+v+fPn53dcHLBrdXV1YrF40KBBACAQCIKDg5OSku77qeDgYP4JMeUdHnEllUqdnJwAwNraOjw8XNHn5wQYVP+N0cbGjLg4s7g486YmPV3xefVVBECJRD+9kYcNd/qrO/sCEhMd2p7TLyrazDCQmjq2Vw/L4pfHHjJkCH9CU708djfkcnlERIT6NoSgoCCGuf8FLkRsaWnZvn07v5ygs7PzyZMn2zWorKwUCoX85alHH330+vXrPf9GjET/jVFEzMtbqc+FwQMCEAA7HNwQ0gMVFTh0KAIoz3yh3tbQkBIXZ8EwJnV1l/uwy9LS0tDQUD4W/f39b9y40Wmz5uZmiUQybNgwvuWMGTM6HqHfV1ZWlnqFl3nz5uXl5bVr8Ntvv40ZM4Y/SyAUCms1f/qKHvXrGFUo8uPiLBnGpKEhUQ/dDR2KAH188jjp71asQACcNavNA5rZ9PQAhoH8/Dc12fH333/v6empntFZX9+6TCLLslKpdPTo0Xz8+fn5SaXSPnfED4EdHR27GgI3NDSEh4dbWFgAwLBhw06fPt3370q/+nWMIuLt228xDGRlPafrjqqrEQBtbe//mHJC2rt4EQUCtLTENstXt3z5Sc7psUlJbp0+l6lX2i7vP2rUqJ9++gkRo6OjJ0yYwAeoj4+PVCrVyvX0tkPgiRMndhwCJyUlTZs2jW8QHBxc+CCMO/p7jDY3l/72m/vq1bNiYzU9I9PcjDIZdjXlIzYWAdDfX8NOSP/T1IRjxyIA7trVurGgAAcORDMzRYLWbomLjY3lc1MgELi6uvJBNmrUqGPHjrV7qorm2g2B293UxLKsRCKxs7MDAHt7+4iICK0XoF39PUYRcfPmTfwp895+sKgIo6JQLMbQUJw8Ga2tEQD37Om8cW0tRkfjzz9rWi3pd/71LwRAPz9sbvOQ+ueeQwBcvFi7XfEzOi0tLYcPH+7g4CAWi9s+1Fa7Oh0Ct9X2yYYzZ85MS7vPdC4DohjFqqqqwYMHA8Cvv/7aVRuWZbOyss6cObNr167ly5cvWBDOP5i37R8TE/Tywn377vlgbi4CoFisetnSggBYUaHL74c8TJKT0dwcTUzwypXWjfwM5IEDdXSiPSAgAAB+/PFHXey8nYSEhClTpqhvarpz5067BidPnnR2dgYAKyurw4cP66GkPqAYRUTcvn07ADz22GPtzv589NFHr7zyyuTJk21sbNqu5zJokCcAurpiUBCGhaFEgr//jvWdPcEoNxetrNDFBXNzESlGSW+9/TYC4FtvtW6pqUE3NwTAQ4d01Cd/s2ZsbKyO9t8OPwS2tbXt6qamqqqqsLAwgUCwatUqHT28WkMUo4iI9fX1/G+8s2fveVzozJkz1dHp7u4+d+7cdevWHT16NDY2tq6uR6fbc3PR3h537cLgYESKUdJbHIfHj2Pb2T+rVyMATp+OOjtdOH78eABITU3V0f47lZ2d/T//8z/8z9qsWbMyMjLaNfDx8QEAQz1AsHsUoyr79u3jp3S0PZktlUoPHz585cqV6uper6qnVGJ6uipGGxvRywtPn6YYJZq5dk21vI1MprtORowYAQA5OTn8S5Zl7zs5Xys4jvvss8/4mwKmTJnS7l1+nb2O8WoMKEZVmpqa3N3dAaDjXRY9VFSE0dEYEaG64mRjgwB47Rra2yMi/vADDh+ONTUUo6RnTp3ChQtx8mRcvBjPnVNtfPppBMDwcJ32zN+dqV57KTMzEwC8vb112qnanTt3XnrppcuX299Q4ObmBgAFBQX6KaNXKEZbHT58mP/norzfSnYcx+XlVZ09q7pMP2mS6jJ92z8CAY4ahWfPqmIUERctws2bKUZJD3z4Ibq746lTKJPh55+jkxMePYqIWF+PO3eizq6e8/grAep5+ImJifxtmjrt9L74J7N2vwqfofSvJ4N2j2VZX1/f9PT0Tz755LXXXmv7VlVVlUwmi4uLS0tLk8lkycnJnp7LUlMPt23j6grjx8O4car/TpwIAwZAXh74+0NVFQDA7dsweTLcvQsVFeBAT8AlXWlqgqFD4cwZePJJ1ZbTp2HNGiguBoFA150jIj8JqaWlxcTEBACuXr0aEBDw+OOPX716Vde9d8PCwkKpVCoUCv42J6PyED5eqs9MTU23bt26fPnybdu2DRs2LCMjIy0tLTU1NS0traampl3j5uaEOXPA1xfGjVP9d+DA++zfwwPWr4eNG3VVP3lIZGeDUgmzZ7dumTcPysqgpAT+vLFdd/jHftjY2PAZym8BgHaTVfRMqVQqlUozMzMjzFDQMEajooBl4c8ZspCZCadOwaZNrQ1KSyEiAlasgLFjAQAaG+Hdd2HnTjA11aRbHXrhhRe2bNlSWlo6b968ttsHDx48bty48ePH8//18/Pjr+zf17BhEBPT+nLtWggKun/gkn6togLs7e8ZeFpbg5UVlJfrJ0bh3tCUy+UAwE9IMhRjqKEbGsXopUvQ0tIaowUFEBl5T4zeuQN798KVK3DpEggE0NQEe/bA9u3GG6MsyyoUCoVC4e/vP23aND8/v3Hjxvn6+vLrKfSBhQVMnNj60twcJk/WTqnkoeXhAeXl0NwM6pHX3bvQ2AiennrovGNgGUOEGcOIuBs6P6h3dwcAOHoUVq3SdVdacOzYscLCQh8fH4Zh1Ac1hOiVhwd4eMCxY60/M//3f+DvD/b2eujcOGPUGGrohqYx2tAAZWWqr/kLKR1FRMAzz8DChWBm3GdiWZbds2cPAGzdupUylBiMiQkcOQLPPw8yGfj5wY0b8NVXcP68fjrvGFjGMBJ8yGP09Gm4fFn1dUMDWFp20mbyZFi8GDZsgH37NOxNt44dO5aZment7R0SEmLoWkj/FhgISUnw1VeQnAxeXpCaqoezojwajfaBpjEaGgoREaqvL1yAv/2t82Y7d4KPD7zwgoa96RDLsmKxGAC2bdtmarTnbkn/MWIEiET679Y4LzEZw4i4G3o6dHVwgD17YMMG/fTWF8ePH79165a3t/fSpUsNXQshBkMH9X2g/RiVyWDFCli5EtLS7tn+yivQ18vdOsey7O7duwFg69atZkZ+BpcQXaKD+j7QKDKmTQOOa305bBiEhEBGBuzbB3I5vPkmHDwIy5er3hUI4MAB2L8fGhrAzg6M6hLOiRMnMjIyRo8evWzZMkPXQoghUYz2gUYx2u5KzNixsGOH6us7d8DODtzcYPfu1ga+vvDXv8LEibB2Lbz1liY9axPLsjt37gSALVu20FCU9HP8IbyxxagxnFjohk7GhE1NsGEDbN7cyVt1dZCXBxs3QlaWLnruiy+//DIjI8PLy2u5euRMSH/Fh2bbwDKGCDOGKO+G9mO0sBBeegnCw8HHp5N3n30Wli+HhgZ4/XUwhkVRWBb27587e/bGrVu301CUEDqo7wPtx+inn8KAAXD4MLzzTucN/vtfcHGBixfh4EGtd95rX30FsbGOBQW7ly//q6FrIcTwjDNGjWFE3A3tj7+2bLlPgyFDQCKBhQtBJIKnnwYvL62X0FMcpzp1u3mzsd9hRYh+GOe8UWOooRuGuV6+YAH89a8glxv40F4qhdRU8PSEF180WA2EGBWaN9oHBpt2tH8/ODtDTAwcOmSYAjgOdu0CANiyBYxyDUNCDMA4D+qNoYZuGCxGhwyBw4cBADZsgOxsAxTw9deQkgIeHhAaaoDeCTFOHQPrxo0bDMMMNOgqucYwIu6GISfBL1gAL7wAcjkIhfo+tEeEnTsBAN55h4aihLTqGKPe3t6TJ0827JpnNBrtzv79MHQopKezx47F3L+19nz9NSQng4cHvPKKPrslxNgZdtzHcdzJkydzcnLabacY7Y6TExw5UlpfP/Kttxbm5+frp1NE1d1WmzbRUJSQe/CB1djYqP+uf/nll6lTpy5fvjw8PLzTqihGu7RggcvcudPr6upWrlypn8eUfvMNJCeDuzu8+qoeeiPkQbJ//34bG5unnnrqP//5D9d2yQxdiomJmTFjxpw5c+Lj4z08PJ566ql2DYz83KhRPKe+vLx86NChACCRSPTQ3ebNaGKCBw/qoStCHjAFBQXz58/nw2HWrFkZGRk67S42NjY4OJjvztHRUSwWNzY2tmtTW1vLB2hpaalOi+kzo4hRRPzyyy8BYODAgfn5+VrfeVYWniWBY94AAAntSURBVD6NHKd6WVSER49iU5PW+yHkIREVFeXm5gYA5ubmIpGoSQc/LTKZLCQkRCAQAICdnZ1IJKqpqWnXRqFQSCQSFxcXU1NTKysrLy+v6OhorVeiOWOJUURcsmQJAAQGBnLqwNOSjz5CADxyRPUyKgqnTdNuD4Q8bKqqqsLCwvgL9L6+vlevXtXWnvPy8oRCIf+MCRsbG5FIVFlZ2a5Nc3PzoUOHhg8fzg9U/f39vb29AUAgEKxcubJje8MyohhVH9ofUQeelnz0EU6dis7OeOcOIsUoIT126dIlHx8fPr+EQmHHAWOvFBQUhIWFWVpaAoCFhYVQKCwuLm7XhmVZqVTKhyaf4FKpFBGbm5vFYrGVlRUAODs7R0ZGalKJdhlRjCLiiRMndHFo/9FHuHw5/uMf+MoriBSjhPQGn1989rm6up46daoPO7l7965IJLK2tgYAExOTkJCQ7Ozsjs2io6MnTpzIB+jYsWMjIyNZlm3bICsrKygoiG8wb968vLy8Pn5XWmVcMYqIixcvBoCgoKBuDu3r6+sLCwtTU1OvX689dw6PH8cDB3DHDly/HletwiVLMCgIJ09GLy90dMRHH1XFaFUVurhgTAzFKCG9lpKS8vjjj/P5FRwcfPv27R5+sK6uTiwWDxo0iB/ShoSEpKend2z2+++/z5w5k9+/u7u7RCJRKpWd7pDjuMjISEdHR/6cgFgsbmlp6fs3pg1GF6N37txxcnICgOeff/6f//ynUCgMCQmZM2fO1KlTvb29nZyczM3N1dMMAgJyALD7P25uqhhFxMhI9PHBb76hGCWk1ziOk0gk/F2hgwYNioiIaDdU7Ojjjz8eMmQI/9M6f/78hISEjm2uX78eGBjIt3Fycur0Yn1HpaWloX/exz1x4sQbN2708bvSBqOLUUTcvn07/7urK7a2tm5ubuPHj1+xInH+fHzxRXzzTdy8Gf/9bzxyBE+dwuhoZBjMzMTyclQqW2OU4/Avf8Gnn6YYJaSPiouL+avBABAQECCTybpp/P777wPAjBkzLl682PHd1NRU9cV6BweH8PDw2traXhVz7tw5Dw8PAHh+zuPRn+1QKu6fv7ogQGNYg/5eCxYsOHv27PTp0xctWmRvbz948GD7e7UdkPbEwYNw5Qp88QUAgEwG/v4waRJcv66T4gnpD86ePfv3v/+9sLDQ3Nx87dq17777Ln/ytJ2mpqbLly+rz2aq5eXl7d69++jRoyzL2trarlmzZuPGjfb29n2opL6+fuvWLR61sTVleYNdRwT/fc+oiTP78i1pwiDh3Y34+HiBQGBra1tWVqatfX7+Ob79duvLrVtx0SJt7ZuQfqq6ulo9I2r06NG//vprTz51+/ZtoVDIP7CHv1hfUlKieTEl2amH3356W7DrtmBXqVgor76r+T57zuhidOHChQCwYcMGQxdCCLm/y5cvjxs3DgAEAkFoaOjdu13mV3l5uUgk4mcsmZmZhYaG5uTkaLEStkV5/cyRnUu8tgW7ipeNZX78XOsz0LtiXAf1iYmJkyZNsra2zs3N5eeQqnEcV11dXV1dXVVVVf0npXJEYWFgdTWo/9TUqL7Ytg3WrTPU90FIP6JUKj/88MPw8HCFQuHi4rJnz56XX365bYO6urqDBw/u2rWrtrZWIBAsWbJkx44djzzyiC6KqSrN//7gxuyE3wBghN+M4Df3DnEbpYuO7qGftO6h5557DgDWr1/fdqNcLu9qydhJk0RdXaDfssVQ3wQh/VFmZqZ6SZHg4GB+6nd9fb1YLB48eDC/PSgoKD4+Xg/FyH6P2vui77Zg1x3Pj4z54n1W2Vx7tyQ74VJ2wqWiW4mKxnrtdmdEo9HU1NQJEyZYWFjk5OS4urq2fcvCwoJlWf76UtsrTu7uASYmr9rbA/9n0CBQf23QtboJ6Y8Q8fDhw/zd8QMHDpw7d+6lS5fKysoAIDAwcOfOndOmTdNbMQ01Fec/CU+5+A0AuHr5jZsxP/b7Tz3GPdZUX1OWm7Zo3X4v/1la60y7qayJ559/HgDWrl3b8a3eToMghBhKSUkJP6Nz1KhRADBt2jQDrieSm3xl/+qA91+eeP3MJyfefZnfePWbjz/btETdRtFQ16Js1qQXY3mssEwm++6776ysrNZ1dkbTzs5O/yURQvrAxcXl2LFjISEhrq6uxcXFCxYsMGAxI/xmrP5PdGVxbtGtBPXGxvpqS2tbANgvnOHpOz0r/mLA4jenPftan3sxlhgNDw/nOO5vf/vbsGHDDF0LIURTzz77rKFLUDG3tHYeOa7oVkJ+Wuzh/31a2dTQWFe9cu8Z/l3nkeMXhH2oYReGX/0eAGSyygsXYqytrTds2GDoWgghD6dhXo+GiCRL/nnId9ZzX+99AzkWAIaNnqD5no0iRt9914Hjsjdt+qbdlSVCCNEWc0vrwS6eziPHBYZuLMlOqa0o1daeDX9Qn5YGp0+DhYX9q68+behaCCEPuZZmRfzPJwY4OA8c4qKtfRo+Rt97DzgOXn8d3NwMXQoh5CFlam5x++Yfe/7qY25pPdRz7PItkQITU0sbOxNTU813buB5ozdvgq8vmJtDVhb8+bwAQgh5kBj43Cg/FF21ijKUEPKgMuRoNDMTfHzA1BRu3QJPT0NVQQghGjHkaHTbNmBZWLWKMpQQ8gAz2Gg0MxPGjQOBADIzKUYJIQ8wg41G33sPWlrgtdcoQwkhDzbDjEazssDHBwQCuHULRozQf/+EEKI1hpk3amEBL74INjaUoYSQB552DurPnwexuPXlnTvwwgv3NGhshKVLISpK9dLDAwQC2LJFK50TQoghaSdGs7Ph2rXWl3I5fPvtPQ1aWuDUKXjzTaiuVm357juoq9NK54QQYkh6vcS0cCH861/67JAQQnROa+dGq6vhxg3V18XFnbcJD4fx4+GVV2D6dG11SwghBqa1GE1Pbz3X2djYeRsHB9i1C954A+LitNUtIYQYmNZidPp0OKNaTxpyc2HMmM6brVwJR4/C4cPa6pYQQgxM3xOeTEzg4EGYPx+am/XcMyGE6IQOLzGlpcGyZbBsGXz33T3b/f1hyRJoaNBdz4QQoj/aGY1OmwaOjq0vBw+G8HBoaIADB8DeHmbPhnnzQCwGkz9De/t2cHO75yOEEPKA0vnNoOXl8PLL8OOPOu2EEEIMRrfzRhUKWL36nhucCCHkIaPDGM3Ph+XL4Z13YIIWnmBKCCFGSocH9Z9+CtevAwAIBHDokI46IYQQAzPwI+0IIeRBZ+BH2hFCyIOOYpQQQjRCMUoIIRqhGCWEEI1QjBJCiEYoRgkhRCMUo4QQohGKUUII0QjFKCGEaIRilBBCNEIxSgghGqEYJYQQjVCMEkKIRihGCSFEIxSjhBCiEYpRQgjRyP8DFGqcgbj3ev0AAAHrelRYdHJka2l0UEtMIHJka2l0IDIwMjEuMDkuMwAAeJx7v2/tPQYg4GVAACkglgbiBkZ2hgwgzczIxMagAWKwwAXYFcACzGwOUAkHVD4bQwKQZmISYFAA0TAujIZrh5jHDLOJiRGigY0DYiETjIab6wDVgMSAyCiDNTJiUwBlCEBUMnEzMCowMikwMmswMbMosLBmMLGyJbCxM7BzMHBwMnByMXBxM3DzaDDx8Crw8jHw8XMw8QkAsSCDoFAGk5BwgrBIBpOIKIOIWIKYeAYTm0QGk4RkgiRLgrhggggTKwsrm4Qkm5CwiJi4oLgTI9BJ8HC96S130GBB8AEQh3m+8EFfr9/7QWyFZwoHPTabgdlKYd8PPHF5D2ZX/Hl4YNFcNjBbd8WmA7ZevftA7FeLVh+4PWeWPYhtLtF44GXHZzB7torhgeIFPWB2ttyO/Q09T8FsbTODfaIetWA2m8oUe4kJ18Hsq7LSDq7+iWC2tXGBw4bJJ8Hs9vkLHErCXMHsJZG7HAoSmRxA7MynbQ5rNv4Fu6FA9aRD0XL2vSD2coXjDvsFj4DdOWHZewepWlGwHxdK8Dqq810Ci0+LVnT0q5YBi7NJ8znKnRUEm6O84bPDklhNsF15C64daGP+AGZXm/848GW7L5gtBgA+D3nsg0ox3gAAAnZ6VFh0TU9MIHJka2l0IDIwMjEuMDkuMwAAeJx9VUtuWzEM3PsUuoAFkuJ30UVsp0FRxAaaNHfotuj9UfIZ7lMAoZJFSMI8itQM5UOr9uPy/def9q/R5XBoDf7zi4j2MQDg8Npq0k7PL9+u7fz+dHrsnG8/r+9vjbSR5TfZP2Of3m+vjx1s13aM7jEErR1HH6Aq1KDD1vZvqZ0LSaSDsB2xR05AFshRPjH3RdFHO0I3HGpjAeWCWjcKUi6nrkEBC6TU8dZBCTTKp5A4r3xq+ZQuPiBXiSSTEFsgrXxKH5rdGnYCC1zl7gXkDoPz9AQGkKIugNHe2pG6qblEAiFceQVEKJfYeTiTJ9JiDFkljlhI6GgEXsFFkuW4QlJmjh3FBtfh6gy+IghHuqTONDQqOHfAZTpY/IxuMpDqgsTYaRmkZOLSAYM3ns3y5CVQ2y2BToOTPMorNYwVN6nGW+NOKAJc+XMIxzLIYkc7uVRsiYQxEGOFjDuSyIbW7WdeFkuxFz3W2aPuPBlNckxXl05Fj3dHl7zW9KmMsFQbUTv9blUXzq5QXpkllbTCjrtXRwnf2E+ynFdI3iJVzHqo+DRnS9JpqyDtBskNJ5lugrAMVEvFWZWitrHpMMR9gXy+Xj69Jvf35XS7Xvb3pTrtj0gu2thfCsqeGXx9+UL7m4A5ZC/8XDTdq5tz2F7CmMP3QsUcsZcj1pirDjeDU3XhZmiqIipTRTKHlWvcDE/VgZuRqQyoDOqk9/uOTcLGbccnAXMZjEmoXIZmQXIZwkl4WIZoEhhvO2OSEZchnuTCeWEkkyy4DOlEP5fZIVYbPEVcvM8s1/rxD5Tzw19iuUY6Ta8byAAAAU16VFh0U01JTEVTIHJka2l0IDIwMjEuMDkuMwAAeJwlkTtuxDAMRK+S0ga0BP+kYGwT906xR1CfImn38CGdwrDxNBrO0Ne5XfvzWvS9tvN1ntfzvF7b82vvZ/Faa/v82dfiff3Sx3t7TMgpxuMhIOhu4yjE7FKIYNabChECmhdCCBIPKRYQPNlblj553gidMVpmbKktM7AUlGYcNi1uJi6ug6DUswcooGiNqpHINeh4MIRHehGceeciUGlPgpgiTRAoGJMHwqwaWU4EZFGHBJ6KOQ4GZYlSZCKRj0MgTKhtLLTuHgZIU1sSgWkNkkV1cIWKEigwmXUDnaazPBw4Le99oAhRxs2YQ3sf5R+1jwDNWRmqS+WNupeQlEatcSUs0ptNna1RtdpTa7JydrlKnto+TrXPCuj10QkfDoEVtkpkGNH/7zCPuyiKZY79/QfK52vF0xXhvAAAAABJRU5ErkJggg==\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"toc": {
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"base_numbering": 1,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
},
"language_info": {
"name": "python",
"version": "3.9.4",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "LewisMartin_FixingValences.ipynb",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment