resonance1

{
 "metadata": {
  "name": "",
  "signature": "sha256:6f1c56bd83c5b48e3eba9511139db05a46dddc7f87aa2a56b3d175470c3839be"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import sys\n",
      "import rdkit\n",
      "from rdkit import Chem\n",
      "from rdkit.Chem.Draw import IPythonConsole"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 38
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol = Chem.MolFromSmiles('[NH2+]=C1C=CC(C=C1)=CN=CN')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 39
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAAbAklEQVR4nO3de3RU5b3G8W+ClIAg\nSJCboCJqFZBLogVCQshkZgIhEdCAWCCCPaKUKtSzLF6wtj21RY72LNt6AY+XoiJV4YiMQCaTnQs3\nqQ0ErXIpxRAIdxJuCeTCzPljmigaJCR7MhP281kryyUzs/cPl3uevHu/7+8N8/l8PkRERCwqPNgF\niIiIBJOCUERELE1BKCIilqYgFBERS1MQioiIpSkIRUTE0hSEIiJiaQpCERGxNAWhiIhYmoJQREQs\nTUEoIiKWpiAUERFLUxCKiIilKQhFRMTSFIQiImJpCkIREbE0BaGIiFiaglBERCxNQSgiIpamIBQR\nEUtTEIqIiKUpCEVExNIUhCIiYmkKQhERsTQFoYiIWJqCUERELE1BKCIilqYgFBERS1MQioiIpSkI\nRUTE0hSEIiJiaQpCERGxNAWhiIhYmoJQREQsTUEoIiKWpiAUERFLuyzYBYg0hfLycvLy8vj8889x\nOp0MGDAg2CWJSIgI8/l8vmAXIRIIn332GW63m4yMDNauXQtA//792bJlC0uWLGHs2LHBLVBEQoJu\njcolo6ioiIULFzJhwgQ6derEoEGD+PDDD7Hb7axdu5ZTp06xceNGXnrpJSZMmMCiRYuCXbKIhADd\nGpVm6+TJk6xevRqPx4PH42HXrl306NGD5ORkXn75ZRISEujUqdN3PnfffffRunVrpk6ditfrZerU\nqU1fvIiEDAWhNBs+n49NmzbVBt+6devw+XzY7XbmzJmD3W7n+uuvr9ex7rnnHsLCwkhPT+f06dPM\nmDEjwNWLSKhSEEpIKykpweVy4XK5yMvL49ChQ0RFRdWG37Bhw2jdunWDjj1x4kTCw8OZPHkyXq+X\nmTNnmly9iDQHCkIJKZWVleTl5dWO+jZv3kz79u2x2+3Mnz8fp9NJ165dTTvfhAkTCA8PZ9KkSXi9\nXh566CHTji0izYOCUIKusLCQ5cuX43K52LBhA5WVlcTHxzN+/HgWLFjAwIEDadGiRcDOn5aWRkRE\nBOPHj8fr9TJr1qyAnUtEQo+WT0iTO378OG63+5xJLtdeey1JSUmkpKSQkJBA27ZtTTuf1+slPz+f\nvn370qZNm/O+b+XKldx111089dRTPPHEE6adX0RCm4JQAs7r9bJ582ZWrFiBy+Viy5YttG7dmpEj\nR2K32y9qkkt9FRcXk5GRURu4paWleDweEhISvvdzq1evZty4cTz55JPMnTvX1JpEJDQpCCVgKioq\neOihh3C73ezevZsbb7wRp9NJUlKS6aO+Y8eOkZmZec4os1evXqSmppKamsrQoUO5/PLL63WsjIwM\nxo0bx2OPPcYvf/lL02oUkdCkIJSAmT59Ol9++SWTJ0/G6XSaOuqrrq4mJyenNvgKCgpo27YtTqfT\nlFFmbm4uKSkpzJw5k3nz5plWt4iEHgWhBITP56NLly68+uqrjBkzxpRjFhUV1S6gNwyD0tJShg4d\nSmpqKna73fRJNWvWrCE5OVlhKHKJUxBKQHz22WdERUVx9OhR2rdv36Bj1NU5pmfPnowaNQq73X7e\nzjFmWrt2LcnJyTz44IPMnz8/oOcSkeDQ8gkJiKysLKKioi4qBL/dOaamUXZDOseYJTY2ltWrVzNq\n1Ci8Xi/PPfdck55fRAJPQSgBYRgGNpvtgu/bt28fLpcLj8dDbm4uhw8fNq1zjFliYmIwDAOn00lZ\nWRkvvfQSYWFhQa1JRMyjW6NiuqqqKiIjI/nggw9wOp3feb2oqIi//OUvuN1uPvnkEyIiIrDZbDid\nTpxOJzfeeGMQqr6wTZs24XQ6SUtL4+WXX1YYilwiNCIU0+Xn51NZWUlsbGydrxcXF7N8+XKcTie/\n/e1viYmJoWXLlk1c5cWLiooiMzMTh8OB1+vllVdeITw89HcyO3nyJNnZ2WRkZGAYBvHx8Tz33HOm\nLl8Rac40IhTTPfPMM7jdbnJzc4NdSkAUFBTgcDgYM2YMCxcuDLkwrGlgUNNQYP369bRq1YoRI0YQ\nExPD4sWL6dChAytXrqRdu3bBLlck6BSEYrrExESGDx/O008/HexSAmbr1q0kJiaSkJDAokWLAtoL\ntT6++uorPvroo9p+rVVVVQwfPrx2TeWgQYNqA/v48eMkJSVRWVlJZmYmkZGRQa1dJNgUhGKq06dP\n07FjR9xuN3FxccEuJ6C2bdtWG/pvvfUWl13WdE8a6urX2qdPn9o1lRfqpHP8+HFGjhzJmTNnyMzM\nDPgyFJFQpiAUUxmGwZgxYzh69Cg/+MEPgl1OwG3fvh2bzUZsbCzvvPNOwMKwrn6tbdu2xeFwNLiT\nTllZGaNHj+bYsWN4PB6FoViWglBMNXfuXPLz81m1alWwS2kyO3bsIDExkSFDhrB48WLTJv4UFxfz\n8ccf4/F4yMnJobS0lBEjRtQGnxmddMrKykhNTeXAgQNkZWXRrVs3U2oXaU4UhGKqoUOHcuedd/Lo\no48Gu5QmVVhYiM1mo2/fvnzwwQe0atXqoo9RUVHBmjVrakd939yeym63Y7PZAvI8r7y8nNTUVPbt\n24dhGApDsRwFoZjmxIkTREZG8sknnxAdHR3scprc7t27sdls3HzzzSxdupSIiIgLfiY/P7/2Od/6\n9eu57LLLaoOvKTvplJeXc8cdd7B3714Mw6B79+5Ncl6RUKAgFNO4XC7uvfdeDh8+HHJLCppKUVER\nNpuNm266iWXLln0nDGv2RfR4PGRmZlJYWFjbSSclJYUhQ4aY+pyxurqajRs3kpGRwdixY4mKijrv\ne8vLyxkzZgxFRUUYhsHVV19tWh0ioUxBKKZ55JFHKCwsZNmyZcEuJaj27NmDzWajd+/evP/++2zc\nuPGc7aK6detGcnIydrudESNGcNVVV5l6/q+++gq32127gL5mKcWTTz553iYHNSoqKkhLS2P79u0Y\nhkGPHj1MrU0kFCkIxTQDBgzg/vvv52c/+1mwSwm6wsJC4uLiCAsLY//+/fzoRz+q3ZT49ttvN3Xd\nYUlJCS6XC5fLRV5eHkePHm3UpJqaZ4YHDhzG48mnW7fQ7/oj0hgKQjHFwYMH6datG//4xz/o06dP\nsMsJCRMnTqSiooI33niDDh06mHbcyspK8vLyakeZmzdvpnv37rWjTDMm1ZSXlzNrVhF5eTdjGKC7\npHIpU69RMUVubi5du3bllltuCXYpIcHn85GXl8fzzz9vSghu3bqVjz76CI/Hw4YNG6iuriYuLo4p\nU6bw3nvvmT6ppk2bNrz44s3cfTfExIBhQO/epp5CJGRoRCimePDBBzl16hRvv/12sEsJCVu3bqVv\n374cOHCAzp07X/Tn6+ocEx0dXXu7MyYmhjZt2gSg8nNVVcHEifDpp/4wvOGGgJ9SpMlpRCimyMrK\n4vHHHw92GSHDMAz69OlT7xD0er2sX7++dm/GgoIC2rVrh8PhYM6cOTgcDnr16hXgqr+rZUt47z1I\nT4eEBMjOVhjKpUcjQmm03bt3c91117F7926uueaaYJcTEu666y569OjBCy+8cN737Nmzh1WrVuHx\neMjOzqakpIShQ4fW9gs1o3OMWc6ehalTweOBrCzQY2C5lGhEKI2WnZ3N9ddfrxD8N6/XS05ODq+/\n/nqdr+/cuZNJkyaRn5/PFVdcgd1uZ968eSQlJYXscoUWLeDNN2HaNLDZ/GHYt2+wq7p427Zto2fP\nnt/bkFysR0EojWYYBjabLdhlhIyCggKOHz/OiBEj6ny9ZobnCy+8YPpSikBq0QLeeAN+8pOvw7Bf\nv2BX9f1KSkrIysrC7XbjdrspLi6mS5cuLFy4kNGjRwe7PAkR1mz/IabKyspSEH5DVlYWUVFRtG/f\nvs7X27Rpw9NPP82QIUOaTQjWaNECXn8dkpMhMRE+/zzYFZ3rzJkzrFixggceeIDevXsTGRnJs88+\nS2RkJK+99honTpzgV7/6FePGjWPx4sXBLldChEaE0ijbtm1j//79JCYmBruUkJGdnX1J//cID/eH\n4cyZ/jD0eKB//+DU4vP52LRpU+3s2nXr1hEREYHdbj/vJKP777+f8PBwpk6ditfrZfLkycEpXkKG\nglAa5WJnR17qaha7//znPw92KQEVFgYvvugPxREjICMDbr+9ac596tQpsrOza9vI7dy5k/79++N0\nOpkzZw6xsbEXbHj+k5/8hDZt2jBt2jR8Ph9TpkxpmuIlJCkIpVEMw7ikRz8X629/+xvV1dUMGzYs\n2KUEXFgY/PGPUF4OK1YELggrKyEvzz/y9HjgBz/wsG/fbJKSknj++edJSEigbdu2F33ce+65h7Cw\nMKZOnYrP5yM9PT0A1UtzoCCUBvN6vWRnZ/Paa68Fu5SQYRgGQ4YMaZLF7qEgPBz+93/9oQhQtQOS\nhkDlo7DyMbji339euRUcY2DqKphWjw41+/eD2+3/ycyE0lJ/h5s774RRo1IZOHAMYTUnbYSJEycS\nHh5Oeno6Pp+Pe++9t9HHlOZHQSgNdqHZkVZkxRm039lx6yys+x38MRGe/BHUJ66OH/eHXs2o76uv\nICoK7HZYvBiGDYPWrWvebe4EowkTJhAeHs7kyZPxer1MmzbN1ONL6FMQSoMZhkFUVJSpDaWbs7Ky\nMjZs2MAzzzwT7FKCqxOM6wV/mgX21TCkjsmzXi+sXw8ulz/4CgqgXTtwOODpp8HphK5dzS1r3759\n591wOC0tjdatW5OWlsbp06f56U9/au7JJaQpCKXBrDj6+T7r16+nVatWDB48ONilBNdlkPxr6HY/\nPPw8ZPwavr18PSYGNm70d6gZORJ+9zsYPhwuMMflopw+fZq8vDwyMjJYvXo1W7duZe/evefdcHj0\n6NEsXbqUtLQ0fD4fM2fONK8YCWkKQmmQmtmRs2fPDnYpIcMwDGJjY03dYb65atEVnv4fsE2B+XZ4\n+lt7D8+fD716Qc+e5p3T6/WyefNmVqxYgcvloqCgoLZzz+zZs3E4HOcNwRrJycksW7aMO++8E5/P\np701LUJXrFzQiRMnyMjIqF2rNWvWLKKioqiurr7gjudWYhgG48ePD3YZIaOzA15Ihx/PhhF/PPe1\n4cPNOce+fftqG5Xn5ORw7Ngx4uPjGT9+PAsWLGhQv9aRI0fy17/+lbvvvpsrr7yNSZOGmFOshCwF\nodRp+/btuN1uVq9eTW5uLlVVVcTExHD//fczcuRIlixZwuDBgy0zO/JCSktLyc/P55VXXgl2KaEj\nHBIeh6lZMPv3EHa28Yc8exby8/0Ta3bt+j/eemsCkZGROBwO/vCHP+BwOOjSpUujz5OamspHH/2d\nsWP7cOoUPPBA42uX0KUgFOC7v1kfOXKEqKgo7HY7jz/+OEOGDDnnlp+eD54rLy+PK6+8kgEDBgS7\nlJASHglzXgBPEmxq4PO/PXu+XkqRlQVlZRAbC6mp8Tz88KcMGDDAlKUU32a392HVKn87uSNH4Mkn\nTT+FhAgFoUXVPOPzeDysWLGCL7/8km7dupGamsqLL75IfHz8ebvF1MyO/M1vftPEVYeu7Oxs4uPj\nCf/OWgLpGAd/+gU4/lS/95eW+meSrljh/+fBgzB0KKSmwi9+AQMGgP93so7//gmcuDhqw9Drhaee\nCujpJEi0H6GF7Nvnb4X1xRdv8+absyktLWXQoEE4nU6SkpKIiYmhZcuWFzyO1+tl2LBh9OvXjwUL\nFujLH7j11luZMWOGpt03woIF/k2A163zL6Ww2/3LKBwOCPbuVOvWwahR8Mgj8KtfBbcWMZ+C8BJ2\n5gysXesPv4wM/04BXbvCpEmfMWjQ5zgcjgb3CN2+fTs2m424uDjefvttS8+UPHDgAN27d2fr1q38\n8Ic/DHY5zdZ//zdUVfnDLyqqjoX6QbZ+vT8MZ8yAefOCXY2YSUF4CfF6YfNm/y0llwu2bIG2bf2/\nUdvt/p/rrzfvfDt27MBmsxETE8M777xTr9HkpWjJkiU8+uij7NmzJ9ilSIDl5/uvp//4D/8SELk0\nKAibueJi+Phj/7OUnBw4dgzi478OvoED/XvImeXMmTPk5eXhdrux2+3cfPPNJCQkEB0dzbvvvmvJ\nMJw+fTpnzpxh0aJFwS5FmkB+vn/UOm0aPPdcsKsRM1j3flYztns3/OUv/tudGzdCp07+31Kff97/\nT7NbU3355ZdkZGTgdrvJzc0lIiKCxMREIiIiuO6668jJySEhIYFx48axdOlSWrVqZW4BIc4wDObO\nnRvsMqSJREf7G4E7nf4ZrC+99HXTcWmeNCJshjZu9M9eczj8F2P//uZeiHv37mXlypV4PB6ys7M5\nefIkcXFxpKSkkJqayvV13F8tKirCZrNx0003sWzZsgvuB3epKCwspFevXhQVFdHTzDYpEvI2b/Zf\ng2lp/jAMtWeaUn8KQhMczoT7lsMbf4ZOwS6mASoqKlizZk1ta6pdu3YRHR2N3W7HbrcTExNTr4Xz\ne/bswWazccMNN7Bs2TJaf71dwCXr9ddfZ968eezYsSPYpUgQFBT4w3DcOHjlFYVhc6Vboxa1Z49/\nEfjy5X/G4/Fw5swZ4uLimDlzJklJSfTt2/eij9mzZ0/WrFlDYmIiY8eO5cMPP7zkwzA7O1uNBSxs\n4MCv1xk+8ID/tqk0P/r9pRGOr4NbW0JnJ7hehKvCoO9/wqGtEH8jvPqtQULFFxB7Pby09fzHrNoB\nto4Q+3s48Y2xeuVWiL8J3vhXw2otK/NPqpk1C265Ba65BpYtq6BHjx68++67HD16lIyMDB555JEG\nhWCNrl27YhgGxcXFjBw5klOnTjX4WKHO5/ORlZWlILS4226Df/3LH4L1uX4rG/H9IIGhIGyE9sPg\n8yo45IaUmXDYB188Dx0ae+CajU0/hYbet66s9M8kfewx/4Xavj08+KB/beFvfgOHD8PSpf7+jElJ\nSaaO3Lp06YJhGJSUlJCcnMzJkydNO3Yo2bZtGwcPHiQxMTHYpUiQtWv3jX8x4fqVpqUgDEWdYNxg\n/8amG4/X/2P5+fDss/5nFh07QkqK/8+mTIEdO/y3QxcsgPHj/TNNA6lz585kZWVx7NgxRo0axYkT\nJwJ7QhMdOnSIxYsXc++995KQkHDe9xmGwa233kpkZGQTVichr4HXrwSPgtAEVzlghZkTZf69sWna\ncf/GpqXn+bWyrAxWroSf/xz69vWP/F58Ea67Dt58E/bv90/znjXL3IX09dW5c2dyc3OpqKgI6TCs\nrKwkJyeHJ554gujoaLp168bDDz9MRUUFU6ZM4XzzydR4XOpUz+tXQoeCMFAOwvQf+pc11PxE9IN1\nR+r38ZqNTctfgflrwFvHe3JzYdIk/6L6WbNg1y4oKoJXX/VP6b7ySlP/Rg1y5ZVXkpmZSVVVFTab\njZKSkmCXhM/nIz8/n2effZbY2Fguv/xyUlNT2blzJw888AD//Oc/OXLkCEuWLOG+++6rc2cDr9dL\nTk6OglDqdMHrt5HfD2IuBWGgdIGF28Hn+/rnzD9g2EUMG2s2Nn19NuQe/e7rDod/1PfeezB9un/H\n71DUoUMH3G43l112GXa7naNH6/jLBFhxcTELFy5kwoQJdOnShdtvv53333+f2NhYVq1axaFDh3jv\nvfeYPn16neskv62goICTJ08SHx/fBNVLc/S9168J3w9iHi2fCGUX2Ni0ZUv/T3PQoUMHMjIyGDly\nJHa7nczMTDoF8EFlzdrI820zNXz48EZt4JqVlUV0dDTtzpklIfINAdiYWAJDI8JgqYQVc+CaNhAW\nAfH/Cbsrvvu2mo1N2xiwNfh3FRulffv2eDweOnTogN1u5/Dhw6Yev+Z2p8PhoGPHjowePZr8/HzS\n09P5+9//zt69e1mwYAHjx49v9C7m2dnZmi0qF9Tg67ee3w9iDgVhkJwsgD//HRZsg2OfQx8D/ueT\nut9bs7Fpm0ugn+Hll1+Oy+WiY8eODB8+nP379zf4WKWlpbz//vukp6fTvXt3brvtNhYtWkR0dDTL\nly+npKSEzMxM5syZQ3R0dKP3TayurmbNmjXMnTuXnTt36vmg1EtDrt+L+X6QxlOLtVBQBa9Ngv2/\ngLm3BbuYplFeXs4dd9zB3r17MQyD7t27X/AzVVVV5Obm4vF48Hg8FBQU0K5dOxwOR207uPo837sY\nu3btwu12k5GRgWEYlJWVMXjwYGbMmMHkyZNNPZdInSz4/dDUFITB5oVNr8D8Q/DHp6CziVsmhbry\n8nLGjBlDUVERhmFw9dVXf+/709PTeeutt+jXrx9Op5OkpCTi4uJMbQZQUlKCy+XC5XKRl5fHwYMH\nz+m7OmzYsEu+bZyEEAt/PzQlBWEQ+cph2X9BRnf4/U8h0oL/k1dUVJCWlsb27dsxDIMePXqc973b\ntm2jXbt2FwzMi1FZWUleXl7tKHPz5s20b98eu91OSkoKTqeTrmbvayVSD/p+aDoKwmApg1dnwb/u\ngF/fAdbawe9clZWVjB8/ni1btpCdnU2vAK8DKSwsZPny5bhcLjZs2EBlZSXx8fG1o76BAwfSwszd\njEUulr4fmpSCMEh2vwuDfwwHv/FnP3bBO6ODVlJQVVZWcvfdd7Np0yYMw6B3796mHfvEiRNkZGTU\njvp27drFtddeS1JSEikpKSQkJNC2bVvTzifSWPp+aFoKQgkZlZWV3HPPPXz66acYhsENN9zQoOP4\nfD42bdpUu7/ili1baN26de0axkBMqhGR5ktBKCHl7NmzpKenYxgGWVlZ9OnTp16fKy4u5uOPP8bj\n8ZCTk8ORI0eIioqqDb7Y2FgiIiICXL2INEcKQgk5Z8+eZerUqWRmZpKVlVXn/ojf1znGbrcTHx9P\n586dg1C9iDQ3CkIJSWfPnuW+++6rfbbXr18/8vPza5/zrV+/nurqaoYPH1476hs0aFCjF82LiPUo\nCCVkVVdXM3nyZNatW8cVV1xRO+pLSkrC6XTicDgC2q9URKxBQSgh7ezZs7jdbr744gucTif9+/cP\ndkkicolREIqIiKXpgYqIiFiaglBERCxNQSgiIpamIBQREUtTEIqIiKUpCEVExNIUhCIiYmkKQhER\nsTQFoYiIWJqCUERELE1BKCIilqYgFBERS1MQioiIpSkIRUTE0hSEIiJiaQpCERGxNAWhiIhYmoJQ\nREQsTUEoIiKWpiAUERFLUxCKiIilKQhFRMTSFIQiImJpCkIREbE0BaGIiFiaglBERCxNQSgiIpam\nIBQREUtTEIqIiKUpCEVExNIUhCIiYmkKQhERsTQFoYiIWJqCUERELE1BKCIilqYgFBERS/t/tsCH\nE4//KYMAAAAASUVORK5CYII=\n",
       "prompt_number": 40,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x259ff30>"
       ]
      }
     ],
     "prompt_number": 40
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# shouldn't this be considered aromatic?\n",
      "for b in mol.GetBonds():\n",
      "    sys.stdout.write('{0:d},{1:d}'.format(b.GetBeginAtomIdx(), b.GetEndAtomIdx()) \\\n",
      "                     + ' arom = ' + str(b.GetIsAromatic()) \\\n",
      "                     + ', conj = ' + str(b.GetIsConjugated()) + '\\n')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0,1 arom = False, conj = True\n",
        "1,2 arom = False, conj = True\n",
        "2,3 arom = False, conj = True\n",
        "3,4 arom = False, conj = True\n",
        "4,5 arom = False, conj = True\n",
        "5,6 arom = False, conj = True\n",
        "4,7 arom = False, conj = True\n",
        "7,8 arom = False, conj = True\n",
        "8,9 arom = False, conj = True\n",
        "9,10 arom = False, conj = True\n",
        "6,1 arom = False, conj = True\n"
       ]
      }
     ],
     "prompt_number": 41
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol=Chem.MolFromSmiles('NC1=CC=C(C=C1)C=NC=[NH2+]')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 42
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAAc7ElEQVR4nO3de1RU5d4H8O8AioKK\ngAHel6KnvGVIZioyDg6ghqIlKErk8YKVlbXKLI/Z8X1Xmb6ndY6apXY0o+sRPRYiylz2ACqWidc0\nTI8BpZgXEBCFAeZ5/5jkZDIj4txgfz9rsTL3s/f+sZby9Zn97N+jEEIIEBERyZSbswsgIiJyJgYh\nERHJGoOQiIhkjUFIRESyxiAkIiJZYxASEZGsMQiJiEjWGIRERCRrDEIiIpI1BiEREckag5CIiGSN\nQUhERLLGICQiIlljEBIRkawxCImISNYYhEREJGsMQiIikjUGIRERyRqDkIiIZI1BSEREssYgJCIi\nWWMQEhGRrDEIiYhI1hiEREQkawxCIiKSNQYhERHJGoOQiIhkjUFIRESyxiAkIiJZYxASEZGsMQiJ\niEjWGIRERCRrDEIiIpI1BiEREckag5CIiGSNQUhERLLm4ewCiBzh+vXryMnJwfHjxxEVFYXBgwc7\nuyQichEKIYRwdhFE9nDs2DFoNBpkZmZi7969AIAHH3wQR48exZdffolJkyY5t0Aicgn8aJRajKKi\nImzYsAHx8fHo1KkTQkJC8NVXX0GtVmPv3r24du0avv32W7z//vuIj49HSkqKs0smIhfAj0ap2aqo\nqMDu3buh0+mg0+lw9uxZdOvWDePHj8cHH3wAlUqFTp063XberFmz0LZtW8ycORMmkwkzZ850fPFE\n5DIYhNRsCCFw6NCh+uDbt28fhBBQq9VYtGgR1Go1evfu3ahrJSQkQKFQICkpCTdu3MAzzzxj5+qJ\nyFUxCMmllZSUID09Henp6cjJycHFixcxZMiQ+vAbOXIk2rZt26RrT5s2DW5ubkhMTITJZML8+fNt\nXD0RNQcMQnIpRqMROTk59bO+w4cPw8fHB2q1GitXrkRUVBSCgoJsdr/4+Hi4ublhxowZMJlMeP75\n5212bSJqHhiE5HQFBQX4+uuvkZ6ejv3798NoNEKpVCIuLg7r16/HQw89BHd3d7vdf8qUKWjTpg3i\n4uJgMpmwYMECu92LiFwPX58ghysrK4NGo7llkUvPnj0RHR2NmJgYqFQqtGvXzqb3zM/PR69eveDp\n6WlxTEZGBp544gm88cYbWLx4sU3vT0Sui0FIdmcymXD48GHs2LED6enpOHr0KNq2bYuxY8dCrVbf\n1SKXphBCYPjw4fD19cX27dvRpk0bi2N3796NyZMn4y9/+QuWLFlit5qIyHUwCMluqqur8fzzz0Oj\n0aCwsBB9+/ZFVFQUoqOj7TLrs+bixYtQq9Xw8fFBRkYG2rdvb3FsZmYmJk+ejNdeew1Lly51WI1E\n5BwMQrKb5ORknDx5EomJiYiKirLrrK8xLl26BLVajfbt2yMjIwMdOnSwODY7OxsxMTGYP38+3nnn\nHQdWSUSOxiAkuxBCIDAwEB9++CFiY2OdXU690tJSREVFAQA0Gg18fX0tjt2zZw/Gjx/PMCRq4dhi\njezi+PHjKCkpwejRo51dyi18fX2h1Wrh5uaGyMhIlJSUWBw7atQo7Nq1C++//z5effVVB1ZJRI7E\nICS70Ov1GDJkCHx8fJxdym06duyIzMxMeHh4QK1W48qVKxbHhoWFYffu3Vi/fj1eeeUVB1ZJRI7C\nICS7kCQJERERzi7Dopth6OnpCbVajcuXL1scO2LECEiShI8++gjPPPMM+DSBqGVhEJLN1dTUIDs7\n22IQFhUV4dtvv3VwVbfz8fGBTqdDx44doVarcenSJYtjQ0NDodVqkZqayjAkamEYhGRzeXl5MBqN\nCAsLa/D4J598gpdfftnBVTXM29sb6enp8PPzg1KpRHFxscWxQ4YMgVarxdatWzFv3jyYTCYHVtp0\nFRUVSEtLw/z589GvXz88/fTTuHbtmrPLInIZDEKyOb1ej2HDhsHLy6vB45IkYcyYMQ6uyrKbYdi5\nc2eoVCqcP3/e4tiQkBDodDps374dycnJLhmGJpMJeXl5ePvttzF69Gj4+/tjxowZKCoqQlJSEvbt\n24dx48ahoqLC2aUSuQS+PkE2N2bMGISHh+PNN9+87diNGzfg5+eH3bt3Q6lUOqE6y65fv47Y2FgU\nFRVBkiR07drV4tgffvgBY8aMgUqlQkpKil17oTbGTz/9hLS0tPp+rTU1NQgPD6/v3BMSEgI3N/O/\ne8vKyhAdHQ2j0QitVgt/f3+n1k7kbAxCsqmbQafRaDBq1Kjbjuv1esTGxqKkpAStW7d2QoXW3bhx\nA7GxsSgoKIAkSejWrZvFsfn5+fWh/8knn8DDw3E97Bvq19q/f39MmDABarUaw4cPh7e3t9Xzx44d\ni6qqKmi12gY3MCaSCwYh2ZQkSYiNjcWVK1caDLolS5bgwIED0Gg0TqiucaqrqxEXF4fjx49DkiT0\n6tXL4thTp04hIiICYWFh+Oyzz+wWhg31a23Xrh0iIyOb3K+1srISjz32GK5evQqdTscwJNniM0Ky\nKUmSEBYWZnG2p9frXfq1CgDw9PTE1q1bMXjwYKhUKpw9e9bi2Pvvvx8GgwG5ublISEhATU2Nzeo4\nd+4cNmzYgPj4eAQFBeHRRx/Fvn37EBcXh2+++QaXL1/Gli1bkJyc3KT2dd7e3ti5cyf8/PwQHh5u\ndaEQUUvGGSHZ1PDhw/H4449j4cKFtx0rLy+Hv78/cnNzMXToUCdUd3eMRiOmTZuGvLw8SJKE4OBg\ni2MLCgoQERGBAQMGYOvWrVa3e7Kkuroae/bsqZ/1/X57KrVajYiICLs8z7t+/TomTJiA8+fPQ5Ik\ndO7c2eb3IHJlDEKymZtB98033yA0NPS24+np6UhKSsKlS5ecvriksWpqajBt2jR89913kCQJffr0\nsTi2sLAQEREReOCBB7Bt2zar2z3dlJeXV/+cLzc3Fx4eHvXBZ8vtqc6ePWv1WtevX8fEiRPxyy+/\nQJIkdOnSxSb3JWoO+NEo2UxOTg46dOiAkJCQBo9LkgSlUtlsQhAAWrVqhS1btiA8PByjRo3CyZMn\nLY7t2bMnDAYDTp06hccffxxVVVW3jSktLUVqairmzZuH3r17Y+jQoUhNTUVoaCgyMzNx5cqVe/q4\nsyGSJGHQoEGQJMniGC8vL6SlpaF79+5QqVQ4d+6cTe5N1CwIIht56aWXxOTJky0eHzx4sFi9erUD\nK7Kd2tpa8eSTT4rAwEDx/fffWx1bVFQk+vTpI6Kjo0V5ebnQarVi0aJFIjQ0VLi7u4tu3bqJ5ORk\nsWXLFnHx4kWH1P+3v/1NeHp6iq+++srquKqqKhETEyP69u0rfv75Z4fURuRs/GiUbGbw4MGYO3cu\nnnvuuduOXbp0CYGBgTh+/DgGDBjghOruXV1dHWbPno1du3ZBr9dj4MCBFscWFBRg1KhRUCgUKC4u\nxiOPPFK/KfHQoUOdMitet24dXnzxRXzxxReYPHmyxXE3nxleuHAJOl0eOndu5cAqiZzA2UlMLcOF\nCxeEQqEQJ06caPD4v/71L9G5c2dhMpkcXJlt1dXViZkzZ4qAgABx7Ngxq2OnTp0qJk2aJEpLSx1T\nXCOsW7dOeHp6iq1bt1odV1lZKebM+UH86U9C/PKLg4ojchI+IySbyM7ORlBQEPr169fgcYPBAJVK\nBYVC4eDKbMvNzQ2bNm3ClClToFQqcfDgwQbHCSGQk5OD+Ph4dOzY0bFFWjFv3jysWbMGCQkJ+PTT\nTy2O8/Lywtq1D6B/f2DECOA//3FgkUQO5rhWGNSi3dx2yVLQ6fV6LFq0yMFV2YdCocDq1auhUCgs\n9hrNz8/HhQsXXKqn6k1z586Fm5sbZs2aBZPJhKSkpAbHtW4NbNkCTJsGqFSAJAFWFs0SNVucEZJN\nWHtRvrCwEKdPn3b5F+nvhru7O9577z088sgjDR6XJAn9+/dHQECAgytrnNmzZ2Pz5s2YO3cuPv74\nY4vjWrUyh+GoUeYwPHPGgUXa2bFjx7B06VJuqUWcEdK9KywsxJkzZ6BWqxs8npWVhV69elltVdbS\nuNoOGw2ZPn063Nzc8NRTT6Gurg6zZs1qcJy7O5CSAsycaQ5EvR7o39+xtdpDVVUV1qxZg6tXr2LV\nqlXN/mN7ajrOCOmeGQwG9O7dGz169GjweHMIBVsymUzIyspqFjPgadOmYd26dfj00xB88onlce7u\nwObNQGQkEBEBnDjhsBJtKj8/H7W1tQCARx55BLm5uUhNTUViYiLq6uqcXB05C4OQ7tnN54OW6HS6\nZhEKtnLkyBGUlZVh9OjRzi6lUf785z/j2WdDMGcOsHat5XHu7sBHHwFjx5rD8PvvHVdjU5WUlCA1\nNRVz585Fz549MXDgQBw7dqz+eL9+/WAwGJCVlYXExMT6kCR5YRDSPbP2fPDHH39EcXGxrGaEer0e\nQ4YMgY+Pj7NLabQpU4Bt24BXXgHee8/yOHd3YNMmYPx4YMwY4Phxx9XYGFVVVdixYwfmzZuH4OBg\n+Pv7Y8WKFfD398fGjRtRXl6OIUOG3HLOAw88AEmSkJOTgxkzZjAMZYhBSPckPz/fatDp9XqXXjRi\nDwaDoVkGf0yMOQwXLgTWrLE8zs3NHIZPPGEOw99NsBxOCIG8vDysWLECkZGR8PPzw1NPPYXS0lIs\nWrQIZ8+excGDBzFnzhz069cPXl5eDV7HnruIUDPg5PcYqZlbu3atGDBggMXjcXFx4oUXXnBgRc5V\nXV0tvL29hUajcXYpTZaRIUTbtkKsWmV9nMkkxPz5Qvj6CnHggGNqE0KIiooKkZaWJp577jnRt29f\noVAoxODBg8XChQuFVqsVN27cuO2cWbNmiT59+oiioiKr1/7pp59Er169RExMjKiqqrLXt0Auhi3W\n6J5MmTIFXbt2xapVq247ZjKZEBAQgI0bNyI2NtYJ1Tne3r17oVarUVJSYnH20RxkZgKTJwPLlwML\nFlgeZzIBc+YA3boB//M/9qnFaARycgCdzvzVuvVXOH/+RURHRyMmJgYqlQrt2rWzeo2amhokJCTg\nwIEDMBgMVrfUasouItTMOTuJqfmqq6sTfn5+Yvv27Q0eP3z4sHB3dxdXr151bGFOtGzZMqFUKp1d\nhk1kZQnRrp0Qb79tfVxdnXl2KIQQxlNCqHyFGPm2EGW/66ZXfVKI8L5CbDrTuHufPy/E5s1CTJ8u\nxH33CeHhIUR4uBBvvSXEoUO1TWrVV1tbKxITE0X37t3F6dOnrY4tLCwUwcHBYty4cQ3OMKll4TNC\narI7rY6UJKnZLRq5V3daQducKJXAzp3A22+bvyxxcwNueQWvDtj3NrD6O6CxHzeVlQGpqcC8eUBw\nMNC1q/k5ZffuwOefA+XlQHY2sHgxEBLi3qR3/tzd3bF582aMHj36jltq9ejRAwaDAadPn8akSZNw\n48aNu74fNR8MQmqym0FnqZem3N4frKysxP79+1vU9xweDvz73+YgbPQK0U7A5GHAmgXAt2UNDzGZ\ngL17gddeAx5+GPD3B5KTgdJS4M03gfPngYMHgXfeAdRqoG1b23w/7u7u+OijjxAdHY2IiAirYdi9\ne3fs2bMHP//8M8aNG4dr167ZpghyOQxCajJrs5+amhrk5OS0mNlRY+Tm5sLT0xPDhg1zdik2FRlp\nbq02aFAjT/AAxi8DppQBL7wLlDYwLRwxwtylZscO88wzIwMoLja3c0tKAoKCbPot3MLd3R2bNm3C\n+PHjERERgRNWugMEBQVh9+7daNNmNJKS2qCBvZapBWAQUpMYjUarQXfw4EEYjUaMHDnSwZU5jyRJ\nCAsLg4dHy+tceLfB5B4EvPl34Po6YOUe4I+tyVeuBIqKzB1q3n0XiIoCHLkmxc3NDRs3bsSUKVMQ\nERGB7610B+jevTtSUv6KM2c8EBUFVFQ4rk5yDAYhNcmBAwdQW1uLsLCwBo9LkoRHH320Wa+cvFst\n6fmgLQREAquSgE0vAtlXbj0WHm5+/udMCoUCa9aswdSpUzFmzBgct/LZb0CAucdqWZm5mQDDsGVh\nEFKTSJKEYcOGWQw6uYVCaWkp8vLyWtTzwXvmBqheB2YK4MXlwCUXbOWpUCiwatUqJCQkYMyYMbe0\nX/uj++4DsrKA6mpg3DjzAh5qGRiE1CR3CroFCxYgISHBgRU5V05ODnx9fTF48GBnl+JS3PyBRasA\nLwn4ocTZ1TRMoVDgH//4B+bMmYPRo0db3GwZAHx9AY0GqKlhGLYkDEK6a0eOHEFubi5UKpXFMRMn\nTkTfvn0dWJVzGQwGKJVKuLnxr9Qf+Y0C1rwKeLn4LkdvvfUWEhMTMX/+e1Z31+jY0RyGdXXm5uMl\nLhrw1HjsLEN3dPXqVej1emg0Gmg0GhQVFaFHjx5Qq9VYv349f/gDGDRoEJ555hk8++yzzi6F7oEQ\nAosXG7Fxoyd0OuDBBy2PLSszzwqrqwGtFvDzc1ydZFv8CUa3qa6uhk6nw4IFCxAcHAxfX18sXboU\nvr6++PDDD1FeXo7du3cjIyMD06dPl323/gsXLuDEiRN8PtgCKBQKLF/uiWefNS/o+eYby2N9fMwB\n2KGD+RWTK1csjyXXxhkhQQiBQ4cOQafTQafTITc3F61bt0ZkZCTUajXUajV69+5923k//vgjIiIi\nMGLECHz22Wdo1aqVE6p3vi+//BILFy7Ezz//7OxSyIZWrDD3Wt21Cxg+3PK4ykpg4kRzMwCt1twc\ngJoXBqFMVVZWwmAw1H/ceerUKQwaNAhRUVGIiorCqFGj0LYR7TwKCgqgUqkQGhqKL774QpZhmJyc\njKqqKqSkpDi7FLKx//s/YNky84v/Vh6J4/p1IDYWuHzZHIadOjmuRrp3DEKZ+GMHfw+PnSgufhZj\nx46FWq1GREQE/Jv4T9nCwkKoVCr0798f27Ztg6enp42rd219+vTBkiVLMHPmTGeXQnbw7rvmtm9p\naebFMZZcuwY89hhw44Z59w5fX8fVSPfIWd2+W7rqk0KE9xFiw6lbf7/qeyFG9hJi7UnL59qqg39x\nsRAffyzEjBlCBAQI0aqVEEqleTeBgwdrRV1d3V1/X5bItVv/Tz/9JADccZ87at4++EAIb28htFrr\n465dE+L1183/peaDi2VclQ06+HfpAqxebd4rbts288c3WVnA668DoaHuNl3t+ftu/ZMnT5ZNt35J\nktC3b190d3abFLKrp582zwwnTTK/OmGJt7e5Qbm3t8NKIxtgELqqJnbwf+kl87F33gEuXPhvB/+w\nMMDeLTBvdusvKiqSzdY1BoNBVh105GzePODvfzc/C0xLc3Y1ZEsMQld1hw7+JSVAz57A6NHmfdoe\ne8z834ICYP16IC7O3B/R0YKCgiBJEs6dO4exY8e26K1rhBDQ6/UMQhmZOxd47z1g6lTgq6+cXQ3Z\nCoPQnn4Fku83b1p686vNQGDf5cadbq2Dv58f8MEHwKVLwP795pVtI0faf9bXGIGBgZAkCSUlJRg/\nfjwqWmiH4vz8fPz66698f1BmZs8GNm4Epk8Htm93Xh2XtMCE54BG/jghKxiE9hQIbDgFCPHfr6rv\ngZF3sbTaWgf/mBjXXZkWEBAAvV6Pq1evYty4cShvRk0ZL168iM8//xxPPfWU1TZykiRh0KBBTV5t\nS83X9Onm5++rVpn/XlPzxiB0dc2gg78lAQEByM7ORnV1tUuHodFoRFZWFhYvXozQ0FB07twZL7zw\nAqqrq/Hkk09CWPhJJ7cdNuhWc+YAkmT+pOf3an4EIvyAsOVA+e/+6Bh/AJR/Aj76z2+/7gt8+OOt\n51afAMJ6A+//YPm+ZfuAQa2AgCggfS1wnwIY8DJw1WbfmfwwCJ3JCOxYBPTwAhRtAOXLQGH17cOa\nQwd/S3x9faHValFTU4OIiAiUuECHYiEE8vLysGLFCoSFhcHb2xsTJkzAmTNnMG/ePJw+fRqXL1/G\nl19+iVmzZkHxx590AEwmE7KyshiEMmdx4XUTVn03ls9I4HgNcFEDxMwHLgngxLtARxvfR04YhE5U\ncQR47yCwPh+4ehzoLwF/t9DbsLl08G9Ix44dodFo4OHhAbVajStOaMp47tw5bNiwAfHx8QgMDMTQ\noUORmpqKsLAw7Nq1CxcvXsSWLVuQnJzcYDu5Pzpy5AgqKiqgVCodUD01O41Y9U2uwwWWVrRMrfsB\n2adv/33PAcDes//9/0z9b7+oAR7uCxR7A63+BEin/nCiAhixDKhcZq+K7atjx47IzMys72Sj1WrR\nyY59qKqrq7Fnzx7odDrs2LEDJ0+eROfOnTFhwgSsXbsW4eHhCAwMbPL19Xo9QkND0b59extWTS3G\nb6u+O881r/rOXAbY+tXC+yKBHZE2vqhMMQhdgQk49CGg7Q+sDnF2Mfbj4+MDnU6HmJiY+jC87777\nbHb9vLy8WxqH19bWIjw8HElJSVCr1QgJCbFZEwGDwcDVomTVzVXfEU8CK9XAm3/8o/7bqvLkP57Y\nHpjuoBrJjB+NOpm4Dmz7C7CuDlj7BhDg7uyK7Mvb2xvp6enw8/NDeHg4iouLm3yt0tJSpKamIikp\nCV26dMHDDz+MlJQUhIaG4uuvv0ZJSQm0Wi0WLVqE0NDQew7B2tpa7NmzB0uWLMGZM2f4fJDuyNqq\nb1usKifb4IzQmSqBfy4A/jMRWDMRkEur6pthOHHiRKhUKkiShC5dutzxvJqaGmRnZ9fP+o4cOYL2\n7dsjMjISf/3rXy1uF3Uvzp49C41Gg8zMTEiShMrKSgwbNgxLly5lENKd3Vz1rTev+lY0o1XfcsIZ\noRMVpgFvbARWxAJtfnvhfsZOZ1flGF5eXkhLS0P37t2hUqlw7ty5O54ze/ZsREZGYufOnVAqlcjI\nyMD58+fvapHLnZSUlCAlJQXx8fEICgpCcHAw/vnPf+L+++/Htm3bUFFRgX379iExMfGe70Xy0ORV\n341cVU73jjNCJ+qZAFxIcHYVzuPl5YX09HRMmTKlfmbYrVs3i+MXL16M5cuXo2vXrjarwWg0Iicn\np36WefjwYfj4+ECtVmPlypWIiopCUFCQze5H8nRz1Xfkmsaf8/tV5SOqgdfizavK/8GFyjbH/QjJ\n6YxGI+Li4nD06FEYDAb06tXLrvcrKCjA119/jfT0dOzfvx9GoxFKpRJqtRpqtRoPPfQQ3N1b+MNa\nal5qgI0zgOJXgSUPO7uYlodBSC7BaDRi6tSpOHToECRJQnBwsM2uXV5ejszMzPpZ39mzZ9GzZ09E\nR0cjJiYGKpUK7dq1s9n9iGzKBBxaB6y8CKyWwYI6Z2AQksswGo1ISEjAd999B0mS0KdPnyZdRwiB\nQ4cOYceOHUhPT8fRo0fRtm3b+ncY7bGohsgexHXg3/8LZHYBlj8L+DME7YJBSC6lrq4OSUlJkCQJ\ner0e/fv3b9R5586dw86dO6HT6ZCVlYXLly9jyJAh9cEXFhaGNm3a2Ll6IhuqBD78bVX5MhmtKncG\nBiG5nLq6OsycORNarRZ6vR4DBgy4bYy1zjFqtRpKpRIBztiQkchGCr8Ahk0Hfv3d701PBz57zGkl\ntVgMQnJJdXV1mDVrVv2zvYEDB1rsHHNz1mfLzjFEJB8MQnJZtbW1SExMxL59+9ChQ4f6WV90dDSi\noqIQGRlp136lRCQPDEJyaXV1ddBoNDhx4gSioqLw4IMPOrskImphGIRERCRrfKBCRESyxiAkIiJZ\nYxASEZGsMQiJiEjWGIRERCRrDEIiIpI1BiEREckag5CIiGSNQUhERLLGICQiIlljEBIRkawxCImI\nSNYYhEREJGsMQiIikjUGIRERyRqDkIiIZI1BSEREssYgJCIiWWMQEhGRrDEIiYhI1hiEREQkawxC\nIiKSNQYhERHJGoOQiIhkjUFIRESyxiAkIiJZYxASEZGsMQiJiEjWGIRERCRrDEIiIpI1BiEREcka\ng5CIiGSNQUhERLLGICQiIlljEBIRkawxCImISNb+H1+PCVDLWFWnAAAAAElFTkSuQmCC\n",
       "prompt_number": 43,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x27c6360>"
       ]
      }
     ],
     "prompt_number": 43
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# if this is aromatic, shouldn't the resonance form above be considered aromatic too?\n",
      "# I think both this and the above should match a 'C=[NH2+]' query twice, both on (0,1) and on (9,10), but not on (7,8)\n",
      "for b in mol.GetBonds():\n",
      "    sys.stdout.write('{0:d},{1:d}'.format(b.GetBeginAtomIdx(), b.GetEndAtomIdx()) \\\n",
      "                     + ' arom = ' + str(b.GetIsAromatic()) \\\n",
      "                     + ', conj = ' + str(b.GetIsConjugated()) + '\\n')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0,1 arom = False, conj = True\n",
        "1,2 arom = True, conj = True\n",
        "2,3 arom = True, conj = True\n",
        "3,4 arom = True, conj = True\n",
        "4,5 arom = True, conj = True\n",
        "5,6 arom = True, conj = True\n",
        "4,7 arom = False, conj = True\n",
        "7,8 arom = False, conj = True\n",
        "8,9 arom = False, conj = True\n",
        "9,10 arom = False, conj = True\n",
        "6,1 arom = True, conj = True\n"
       ]
      }
     ],
     "prompt_number": 44
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol = Chem.MolFromSmiles('C=C1C=CC(C=C1)=C')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 45
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAAL/klEQVR4nO3dbUyV5QPH8Z9lG22u\nRYq0GCy3RkCiKRhUPoyKwhEvgLAXhrlMe8QQUB46yUqZoyGuzMy0EsK2iKRNUIyZlJLBMAqJp2xZ\nxgYrJmPGRLLzf6Fs/vUYZniuc7y+n5fCi98rv9zn3Pd1j3M6nU4BAGCp60wPAADAJEIIALAaIQQA\nWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAAqxFC\nAIDVCCEAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUIIALAa\nIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAA\nqxFCAIDVCCEAwGrjTQ8AvE1NTY2Ki4tNz3ApIyNDcXFxpmcAXoUQAv/Ct99+q4SEBD3//PPy9/c3\nPef/9PT0KCEhQY2NjZoxY4bpOYDXGOd0Op2mRwDewOl06r777lN4eLjeffdd03NcWrp0qX744QfV\n19dr3LhxpucAXoHvCIHLtGPHDnV0dGjt2rWmp1xSQUGB2tra9NFHH5meAngNQghchpMnTyonJ0cO\nh0OTJ082PeeSJk+eLIfDoezsbP3555+m5wBegRACl6GwsFATJkzQ8uXLTU8Z1fLlyzVhwgQVFhaa\nngJ4Bb4jBEZx7NgxhYaGqqKiQvHx8abnXJaqqiotWLBAbW1tuv32203PATwaIQRGsWDBAp08eVK7\nd+92+fMTJ06ovb3dzavOCg0Nla+vr8ufzZ8/XzfddJM+/vhjN68CvAshBP5BXV2dHn74YbW0tCgk\nJMTl71RXVysxMdHNy86qrKy85FVqe3u7pk+frtraWs2bN8/NywDvQQiBSzhz5owiIiIUExOjDRs2\nmJ5zRdLT0/Xll1+qqalJ119/vek5gEfiZhngEt577z11d3crPz/f9JQrlp+fr99++03vv/++6SmA\nxyKEgAv9/f1yOBx67bXXdPPNN5uec8V8fX316quvyuFwqL+/3/QcwCPx0SjgQmZmpvbt26fDhw97\n/UeKZ86c0cyZMxUbG6uioiLTcwCPQwiBC3R2dmratGmqqalRTEyM6TljYv/+/YqLi1NLS4vuvPNO\n03MAj0IIgQs8+uij8vHxUUVFhekpYyo5OVmnT5/Wrl27TE8BPAohBM6zZ88eJSUlqa2tTVOmTDE9\nZ0z9/PPPCgsL086dOzV//nzTcwCPQQiBc4aHhzVt2jQlJSWpoKDA9JyrIi8vT5WVlWppadENN9xg\neg7gEbhrFDhn06ZNGhgYUG5urukpV01eXp4GBgb09ttvm54CeAyuCAFJv//+u4KDg/XGG29o0aJF\npudcVSUlJVqxYoW6uro0adIk03MA4wghIOm5555Tc3OzDh06dM2/0NbpdCo6OloRERFcGQIihIC+\n//57RUZG6sCBA4qOjjY9xy2++eYbzZkzR01NTZo+fbrpOYBRhBBWczqdmjt3rqZMmaLS0lLTc9wq\nNTVV3d3d+uKLL0xPAYwihLDap59+qsWLF6ujo0MBAQGm57hVd3e3QkJCVFJSoqSkJNNzAGO4axTW\nOnXqlFauXKmcnBzrIihJAQEBys7OVlZWlk6dOmV6DmAMIYS11q9fL+nsuaK2ysrKktPpVHFxsekp\ngDF8NAor/frrrwoNDVVJSYkee+wx03OM+uSTT/TUU0+ps7NTt912m+k5gNsRQlhp0aJFOn78uPbv\n3296ikeIiYlRUFCQSkpKTE8B3I4QwjoNDQ2aPXs2jw6cZ+QRkoMHDyoqKsr0HMCtCCGs8vfff+ue\ne+5RZGSk3nnnHdNzPMozzzyjlpYWff3119f8oQLA+QghxtSRI0d05MgR0zMu4u/vrwcffFClpaV6\n6aWX1NXVJT8/P9OzPMrIMXNvvvmmUlNTtW/fPvX29pqedZHw8HCFh4ebnoFryHjTA3Btqaur08aN\nG03PuEhUVJRmzZqlVatWafXq1UTQBT8/P73yyivKzc1VYmKitm/froaGBtOzLpKWlkYIMaa4IoQ1\nHA6HPvvsM3333XcaP56/AV0ZHh5WeHi4UlJStGbNGtNzALcghLDC0aNHNXXqVFVWVvJS2lHs3r1b\nycnJ1+TLiQFXCCGssGLFCjU3N6uurs70FK8wb948RURE8KA9rMDJMgAAq3FFCCuMHDC9fft2JScn\nm57j0SoqKrRkyRJ1dXXJ39/f9BzgquOKEFYICAhQTk6OVq5cyQHT/2BwcFAZGRnKzs4mgrAGV4QY\nUxs3bvTYxye2bt2qsLAwLVmyRC+//LLpSR5p7dq1KikpUWtrq55++mmPfXwiLS3N9AxcQwghxpSn\nP1BfUVGhJ598Uu3t7QoKCjI9y6OMHEReWlqq5ORkHqiHNQghrBMTE6PAwEDr3kg/mtTUVPX29urz\nzz83PQVwK0II63DA9MUOHjyoBx54QM3NzbrrrrtMzwHcihDCSs8++6yamprU2Nio666z+56xkYPI\nZ82apc2bN5ueA7id3f8DwFpr1qzRTz/9pLKyMtNTjPvwww917NgxFRQUmJ4CGEEIYSU/Pz+tXr1a\nq1at0sDAgOk5xgwMDCg7O1sOh0O33HKL6TmAEYQQ1kpLS9OkSZP0+uuvm55iTGFhofz8/PTiiy+a\nngIYw3eEsNqePXuUmJio1tZW3XHHHabnuNWPP/6oqVOnqrq6Wg899JDpOYAxhBDWi4+P14033qiK\nigrTU9wqOTlZQ0NDqqqqMj0FMIoQwnojr2jatWuXYmNjTc9xi9raWiUkJFh5JQxciBACkjIzM7V3\n714rXtr7119/6e6771ZcXJyKiopMzwGM42YZQFJ+fr76+vq0bds201Ouuq1bt6qvr0/5+fmmpwAe\ngStC4JwtW7YoLy9PXV1dmjhxouk5V0VfX5+Cg4O1bt06LVu2zPQcwCMQQuCckRNW5syZow0bNpie\nc1Wkp6ervr5eDQ0N1p+oA4wghMB56uvrFRMTc02eudna2qoZM2aorq5O999/v+k5gMcghMAFUlJS\n1N/fr9raWtNTxlRsbKx8fX1VXl5uegrgUQghcIHjx48rJCRE5eXlio+PNz1nTFRVVenxxx9XR0eH\nAgMDTc8BPApfEgAXCAwMVEZGhtLT0zU0NGR6zn82NDSk9PR0ZWZmEkHABUIIuJCbm6vTp09r06ZN\npqf8Z2+99ZaGh4eVm5tregrgkfhoFLiEsrIyvfDCC+rs7NStt95qes4V6enpUXBwsDZv3qyFCxea\nngN4JEIIXILT6dTcuXMVFhamLVu2mJ5zRZYtW6b29nZ99dVXGjdunOk5gEcihMA/OHz4sKKjo3Xo\n0CFFRka6/J3q6molJia6edlZlZWVl7yhp6mpSffee68aGho0c+ZMNy8DvAchBEaxePFiHT16VAcO\nHHB5VXXixAm1t7cbWCaFhobK19f3on93Op2aPXu2goOD9cEHHxhYBngPQgiMore3V8HBwdq2bZtS\nUlJMz7ks5eXlWrp0qbq6uuTv7296DuDRuGsUGIW/v79ycnKUmZmpwcFB03NGNTg4qKysLOXm5hJB\n4DIQQuAyZGZmysfHR+vXrzc9ZVRFRUXy8fFRRkaG6SmAV+CjUeAy7dy5U0888YQ6OjoUFBRkeo5L\nv/zyi0JDQ7Vjxw5jN/AA3oYQAv/CI488Ij8/P5WVlZme4tLChQv1xx9/aO/evaanAF7j2n4VNzDG\n1q1bp6ioKE2cONHjvn/r6elReXm5GhsbTU8BvApXhMC/VFNTo+LiYtMzXMrIyFBcXJzpGYBXIYQA\nAKtx1ygAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUIIALAa\nIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAA\nqxFCAIDVCCEAwGqEEABgNUIIALAaIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDVCCEAwGqEEABgNUII\nALAaIQQAWI0QAgCsRggBAFYjhAAAqxFCAIDV/gcQ1ehrDRBotgAAAABJRU5ErkJggg==\n",
       "prompt_number": 35,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x27c6050>"
       ]
      }
     ],
     "prompt_number": 35
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# I think this shouldn't be considered aromatic, because the ring bonds cannot be delocalized\n",
      "for b in mol.GetBonds():\n",
      "    sys.stdout.write('{0:d},{1:d}'.format(b.GetBeginAtomIdx(), b.GetEndAtomIdx()) \\\n",
      "                     + ' arom = ' + str(b.GetIsAromatic()) \\\n",
      "                     + ', conj = ' + str(b.GetIsConjugated()) + '\\n')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0,1 arom = False, conj = True\n",
        "1,2 arom = True, conj = True\n",
        "2,3 arom = True, conj = True\n",
        "3,4 arom = True, conj = True\n",
        "4,5 arom = True, conj = True\n",
        "5,6 arom = True, conj = True\n",
        "4,7 arom = False, conj = True\n",
        "6,1 arom = True, conj = True\n"
       ]
      }
     ],
     "prompt_number": 47
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol = Chem.MolFromSmiles('O=C1C=CC(C=C1)=O')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 48
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAANWElEQVR4nO3dfUzVdcPH8Q+V942b\nl0aKtBwsNyMgxVQUKrVhWRhxb2DodhvmXWrPRoDy0FGuUudNIa4HMtMeIG2TUNsAxZxpKRkMUrlI\nHqRpGRtcxSUxIxHt3H+ou80Lk0DP9xy/79dfzqPHz+Zvvj3n/H6/4+V0Op0CAMBS15keAACASYQQ\nAGA1QggAsBohBABYjRACAKxGCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRACAKxG\nCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRACAKxGCAEAViOEAACrEUIAgNUIIQDA\naoQQAGA1QggAsBohBABYjRACAKxGCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRAC\nAKxGCAEAViOEAACrEUIAgNVuMD0A8DSlpaXKyckxPaNbSUlJioqKMj0D8CiEEPgLvvnmG8XExOiZ\nZ56Rn5+f6Tl/0NzcrJiYGFVUVGjMmDGm5wAew8vpdDpNjwA8gdPp1N13361Ro0bp3XffNT2nW/Pm\nzdO3336rsrIyeXl5mZ4DeAQ+IwR6aMOGDaqrq9OyZctMT7mk5cuX69ChQ/r4449NTwE8BiEEeuDE\niRNKS0uTw+HQ0KFDTc+5pKFDh8rhcCg1NVW//vqr6TmARyCEQA9kZWVpwIABWrBggekpl7VgwQIN\nGDBAWVlZpqcAHoHPCIHLOHr0qIKDg1VYWKjo6GjTc3qkuLhYM2bM0KFDh3TrrbeangO4NUIIXMaM\nGTN04sQJbd26tdvHjx8/rtraWhevOis4OFg+Pj7dPjZt2jQNHDhQGzdudPEqwLMQQuBP7N69Ww88\n8ICqq6sVFBTU7a8pKSlRbGysi5edtWXLlku+Sq2trdXo0aO1Y8cO3XvvvS5eBngOQghcwpkzZzRu\n3DhFRkZq1apVpuf0SmJior744gtVVlbq+uuvNz0HcEucLANcwnvvvaempiZlZmaantJrmZmZ+vHH\nH/X++++bngK4LUIIdKOtrU0Oh0OvvPKKbrzxRtNzes3Hx0cvv/yyHA6H2traTM8B3BJvjQLdSE5O\n1s6dO1VVVeXxbymeOXNGY8eO1dSpU5WdnW16DuB2CCFwkfr6eoWGhqq0tFSRkZGm51wRu3btUlRU\nlKqrq3X77bebngO4FUIIXOThhx+Wt7e3CgsLTU+5oqZPn65Tp06pqKjI9BTArRBC4ALbtm1TXFyc\nDh06pOHDh5uec0UdOXJEISEh2rx5s6ZNm2Z6DuA2CCFwTldXl0JDQxUXF6fly5ebnnNVZGRkaMuW\nLaqurla/fv1MzwHcAmeNAufk5uaqvb1d6enppqdcNRkZGWpvb9fbb79tegrgNnhFCEj66aefFBgY\nqNdff12zZ882PeeqysvL04svvqiGhgYNGTLE9BzAOEIISHr66ae1f/9+7du375r/Qlun06mIiAiN\nGzeOV4aACCGggwcPKiwsTHv27FFERITpOS7x9ddfa9KkSaqsrNTo0aNNzwGMIoSwmtPp1OTJkzV8\n+HDl5+ebnuNSCQkJampq0ueff256CmAUIYTVNm3apDlz5qiurk7Dhg0zPcelmpqaFBQUpLy8PMXF\nxZmeAxjDWaOw1smTJ7Vw4UKlpaVZF0FJGjZsmFJTU5WSkqKTJ0+angMYQwhhrZUrV0o6e19RW6Wk\npMjpdConJ8f0FMAY3hqFlX744QcFBwcrLy9PjzzyiOk5Rn3yySd6/PHHVV9fr1tuucX0HMDlCCGs\nNHv2bB07dky7du0yPcUtREZGKiAgQHl5eaanAC5HCGGd8vJyTZw4kUsHLnD+EpK9e/cqPDzc9BzA\npQghrPL7779rwoQJCgsL0zvvvGN6jlt58sknVV1dra+++uqav6kAcCFCiL47/S8pf4n0vx9Lh3+R\nhtwmzUyWMudIvudu7FxRITU29vqP+KF/f+397bde/34/Pz/dd999ys/P1wsvvKCGhgb5+vr2+vmu\nRedvM/fGG28oISFBO3fuVEtLS6+fL6x/fwX24e9MI0ZIEyac/XFPjjGgl24wPQAezvmL9Gqc9P6N\n0ptfSZMDpB+/lF56XopulLaukIZcJxUVSRs39vqPOR4aqr9XV/f694eHh2v8+PFatGiRlixZQgS7\n4evrq8WLFys9PV2xsbH68MMPVV5e3uvnWxIaqsA+/J1p5syzIezpMQb0Eq8I0TeH35OiVksf7ZTu\nHvT/P9/+tTRlpjR7q7TgDnP7LuBwOPTpp5/qwIEDuuEG/g/Yna6uLo0aNUrx8fFaunSp6TlnedAx\nBs/Ef6PQB07pH1ulAVOkOwb98aGBI6X/8pW2VklnzKy7UGNjo7Kzs/Xaa68RwT/Rr18/5eTkKDs7\nW0eOHDE9R550jMFzEUL0QZf0z1bpb7dI/3nxY97SsEHSLy3SKRPb/ig3N1cRERF8M3sPPPTQQ5ow\nYYLefPNN01PkSccYPBchRB/0k24aJP36z27+ITopNf0iDRwq/YeJbbg2cIzh6iOE6AMv6c4YqW2b\ndPD4Hx9qr5GKfpIeGiddb2bdhVJSUlRVVaVNmzaZnuL2CgsLdeDAAaWmppqeIk86xuC5CCH65rZH\npLmDpDn/I5XWS7+dlL7bKc2dLTlnSP8dYnqhpLM3mE5LS9PChQu5wfSf6OjoUFJSklJTU+Xn52d6\nzlkecozBc3HWKPru9L+kDUulrA1S/c+Sz23SjBelvz8hDT13jdfixX26fOJgaKji+3j5xNq1axUS\nEqInnnhCL730Uq+f61q2bNky5eXlqaamRnPnzu3z5ROP9vXyifNnrvbkGAN6iRDCNdzkgvrCwkI9\n9thjqq2tVUBAQK+f71p0/kbk+fn5mj59untdUA9cRYQQ1omMjJS/v79130h/OQkJCWppadFnn31m\negrgUoQQ1uEG0/9u7969mjJlivbv36877uDidNiFEMJKTz31lCorK1VRUaHrrrP7nLHzNyIfP368\nVq9ebXoO4HJ2/wsAay1dulTfffed1q9fb3qKcR999JGOHj2q5cuXm54CGEEIYSVfX18tWbJEixYt\nUnt7u+k5xrS3tys1NVUOh0M33XST6TmAEYQQ1nr++ec1ZMgQvfrqq6anGJOVlSVfX18999xzpqcA\nxvAZIay2bds2xcbGqqamRiNGjDA9x6UOHz6skSNHqqSkRPfff7/pOYAxhBDWi46OVv/+/VVYWGh6\niktNnz5dnZ2dKi4uNj0FMIoQwnqNjY0aOXKkioqKNHXqVNNzXGLHjh2KiYmx8pUwcDFCCEhKTk7W\n9u3brfjS3tOnT+vOO+9UVFSUsrOzTc8BjONkGUBSZmamWltbtW7dOtNTrrq1a9eqtbVVmZmZpqcA\nboFXhMA5a9asUUZGhhoaGjR48GDTc66K1tZWBQYGasWKFZo/f77pOYBbIITAOefvsDJp0iStWrXK\n9JyrIjExUWVlZSovL7f+jjrAeYQQuEBZWZkiIyOvyXtu1tTUaMyYMdq9e7fuuece03MAt0EIgYvE\nx8erra1NO3bsMD3lipo6dap8fHxUUFBgegrgVgghcJFjx44pKChIBQUFio6ONj3niiguLtbMmTNV\nV1cnf39/03MAt8KHBMBF/P39lZSUpMTERHV2dpqe02ednZ1KTExUcnIyEQS6QQiBbqSnp+vUqVPK\nzc01PaXP3nrrLXV1dSk9Pd30FMAt8dYocAnr16/Xs88+q/r6et18882m5/RKc3OzAgMDtXr1as2a\nNcv0HMAtEULgEpxOpyZPnqyQkBCtWbPG9JxemT9/vmpra/Xll1/Ky8vL9BzALRFC4E9UVVUpIiJC\n+/btU1hYWLe/pqSkRLGxsS5edtaWLVsueUJPZWWl7rrrLpWXl2vs2LEuXgZ4DkIIXMacOXPU2Nio\nPXv2dPuq6vjx46qtrTWwTAoODpaPj8+//bzT6dTEiRMVGBioDz74wMAywHMQQuAyWlpaFBgYqHXr\n1ik+Pt70nB4pKCjQvHnz1NDQID8/P9NzALfGWaPAZfj5+SktLU3Jycnq6OgwPeeyOjo6lJKSovT0\ndCII9AAhBHogOTlZ3t7eWrlypekpl5WdnS1vb28lJSWZngJ4BN4aBXpo8+bNevTRR1VXV6eAgADT\nc7r1/fffKzg4WBs2bDB2Ag/gaQgh8Bc8+OCD8vX11fr1601P6dasWbP0888/a/v27aanAB7j2v4q\nbuAKW7FihcLDwzV48GC3+/ytublZBQUFqqioMD0F8Ci8IgT+otLSUuXk5Jie0a2kpCRFRUWZngF4\nFEIIALAaZ40CAKxGCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRACAKxGCAEAViOE\nAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRACAKxGCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1\nQggAsBohBABYjRACAKxGCAEAViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABYjRACAKxGCAEA\nViOEAACrEUIAgNUIIQDAaoQQAGA1QggAsBohBABY7f8A5190vCczv4EAAAAASUVORK5CYII=\n",
       "prompt_number": 50,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x27c6600>"
       ]
      }
     ],
     "prompt_number": 50
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# This is not considered aromatic (correct)\n",
      "for b in mol.GetBonds():\n",
      "    sys.stdout.write('{0:d},{1:d}'.format(b.GetBeginAtomIdx(), b.GetEndAtomIdx()) \\\n",
      "                     + ' arom = ' + str(b.GetIsAromatic()) \\\n",
      "                     + ', conj = ' + str(b.GetIsConjugated()) + '\\n')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0,1 arom = False, conj = True\n",
        "1,2 arom = False, conj = True\n",
        "2,3 arom = False, conj = True\n",
        "3,4 arom = False, conj = True\n",
        "4,5 arom = False, conj = True\n",
        "5,6 arom = False, conj = True\n",
        "4,7 arom = False, conj = True\n",
        "6,1 arom = False, conj = True\n"
       ]
      }
     ],
     "prompt_number": 49
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol = Chem.MolFromSmiles('NC=CC=O')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 51
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAAN60lEQVR4nO3df2zUdZ7H8RdMW0r5\n2R/IUWV72F8UcNNqK1QCVNKLt2vx9rwtf5At2VtqMe7udN0zW045p7cb3GaD2Xa5Mw45SYrcqeWy\niSBsNwRDuRTbjeciFGxxy1qIUNBKtbL0B533/dHg8VNnoDND/TwfSRPS+X5m3uWP7zPf73fmO+PM\nzAQAgKPGR3sAAACiiRACAJxGCAEATiOEAACnEUIAgNMIIQDAaYQQAOA0QggAcBohBAA4jRACAJxG\nCAEATiOEAACnEUIAgNMIIQDAaYQQAOA0QggAcBohBAA4jRACAJxGCAEATiOEAACnEUIAgNMIIQDA\naYQQAOA0QggAcBohBAA4jRACAJxGCAEATiOEAACnEUIAgNMIIQDAaYQQAOA0QggAcBohBAA4jRAC\nAJxGCAEATiOEAACnEUIAgNMIIQDAaYQQAOA0QggAcBohBAA4jRACAJxGCAEATiOEAACnEUIAgNMI\nIQDAaYQQAOA0QggAcBohBAA4jRACAJxGCAEATiOEAACnEUIAgNMIIQDAaYQQAOA0QggAcBohBAA4\njRACAJxGCAEATiOEAACnEUIAX/jwww+jPQIQcYQQgCTp1KlTyszMVHl5ubq7u6M9DhAxhBCAJCk1\nNVXvvvuuPvvsM6Wnp2vdunXq6+uL9lhA2BFCOGt4eFjDw8PRHiPqLv8/yMzMVENDg15//XXt2rVL\nOTk52rx5swKBQBQnBMKLEMI5/f392rBhg7Kzs5WZmant27dHe6SoaWhoUFZWltra2q74fXFxsd55\n5x09++yzeuaZZ7Ro0SI1NzdHaUogvAghnBEIBLR582alp6dry5Yt2rhxoyorK1VeXq78/HyndvTN\nzc3Kz8/XY489Jq/Xq+zs7Gu2iY2NVUVFhTo6OlRYWKiioiKtXLlSXV1dUZgYCCMDHNDc3Gz5+fk2\ndepUq62ttYGBgS8eO3XqlFVUVFhMTIyVlpZaV1dXFCcNrw8++MBKS0stJibGKioq7NSpU0GvPXr0\nqD300EOWkJBgVVVV1tfXF8ZJgcjhiBBfaydOnNDKlSu1bNky3XvvvWpvb1dlZaXi4uK+2GbWrFny\n+/1qaWnR6dOnlZOTo+rqavX390dx8tHV39+v6upqzZs3T6dPn1Zra6v8fr9eeOEFNTY2BvUcOTk5\namxs1KuvvqqGhgbl5ORo69atMrMwTw+EWbRLDITDhQsXzOfzWUJCgi1evNjefvvtoNYFAgFraGiw\ntLQ0mz17ttXX11sgEAjztOETCASsvr7e7rrrLktLS7MdO3Zc8dhTTz1lcXFx9vDDD1tHR0fQz9vX\n12dPP/20xcfH25OlpWaHD4djfCAiCCG+durr62327NmWlpZmDQ0NNxWy8+fPm8/ns4kTJ1pRUZEd\nPHgwDJOG18GDB23ZsmU2efJkq6mpsfPnz193u5MnT1pZWZl5PB4rKyuzM2fOBP0ax48ft84f/9gs\nJsbsiSfMPv54tMYHIoYQ4mvj4MGDVlRUZBMnTjSfz3fDHX8oLkVi/PjxIUciWs6cOXPFzCdPngxq\nXWtrqxUWFlpiYqLV1NRccR31K/3hD2YPPGA2aZKZz2d24cLNDQ9EASHEmHfu3Gf25JNPWmxsrJWU\nlIR0ii9Y+/bts9zcXJs+fXrokYiQ/v5+q6mpsenTp9vChQutpaUl5Oe4dCp11qxZlp2dbbt27Qpl\nsVlDg9ns2WYZGSP/BsYAQogx6+JFsxdfNEtL67PFix+0xsbGkNbv2LHD9u7dG/T2w8PDVl9fbzNn\nzrSsrCx74403Qh05bHbu3GlZWVk2c+ZMq6+vt+Hh4Vt6vs8//9x8Pp/Fx8dbcXGxHTlyJJTFI0eF\n8fFmy5ebHTp0S7MA4UYIMSbt3GmWlWU2c6ZZfb3Zzez3q6urLSYmxh599FE7fvx40Ov6+vrM5/PZ\nhAkTrLi42Nra2kJ/8VHS1tZmxcXFNmHCBPP5fKP+kYYTJ05YWVmZxcXFmdfrtd7e3lAWm5WVjVw/\nrKgwO3t2VGcDRgshxJjS1mZWXGw2YcLIQcet7vfff/99Ky0ttdjYWPN6vXbu3Lmg1x47dsxKSkq+\nWBtSJG5Rb2+veb1ei42NtdLSUjt27FhYX6+lpcUWLVpkycnJVltbaxcvXgx+8Ztvmn3zm2aJiWa1\ntWZDQ+EbFLgJhBBjQm+vmddrFhtrVlpqNtr7/b1799o999xjSUlJIe/o9+zZY/Pnz7+5SITo4sWL\nVltba8nJybZgwQLbs2dP2F7rakNDQ7Zp0yYrzc214aIisyA/kmJmI4fs9fVmM2aYzZ1rtnt3+AYF\nQkQIcVu7eHHkICI52WzBArNw7veHhobM7/fbjBkzLDc315qamoJeOzg4aH6/31JSUiwvLy+ktcFq\namqyvLw8S0lJMb/fb4ODg6P+GkHp6TH70Y/M4uLMfvADs9Ong1/b3W1WXm7m8Zhd9plGIJoIIW5b\nTU1mublmKSlmfn/kzqh98sknVlVVZXFxcVZSUhLS9cOenh7zer3m8XhCXnsjx48ft5KSEvN4POb1\neq2np+eWn3NUvP/+yOF5QoJZVVVo56n/+Eez2/Cdt3ATt1jDbefPf5ZWrJCWL5eWLpU6OqSKCikm\nJjKvn5iYqJqaGh06dEiBQEDz588P+rv5kpKSVFdXp8OHD2toaCiktVfr6+vTunXrNH/+fA0NDenw\n4cOqq6tTUlLSzfxZoy8jQ2pokHbulHbvlubOlbZulYK55VpurnTZbe6AqIp2iTF2DBw1W5phtvmq\nj+n1t5ktnmP270dvvHaww+zBRLPFz5l9etmNXgaOmi3NNNvyJ7PBQbMNG8wmTzYrLDRrbQ3P3xGq\nPXv22Lx58+zOO+8M+ZZrO3bssDlz5oS09tJn+VJTU+3uu+++4rZot62hoZHD9hkzzO6/3+ytt4Jc\n12P20g/NMhPNNN4sJdvsh5vNzkbptC+cRAgRtFsO4VQzTTb7RavZpRxcHsLhYbOSErOXXx75bPbt\nZHBw0Gpra23atGlWUFBgbwW7ozezgYEBq62ttSlTplhBQYEdOHDghtseOHDACgoKbMqUKdd8S8aY\n0NPz/+9qKisbuSZ4I4Fesw3LzNL/zmz3e2afnzdr/53ZP2SYFfzM7KNb+ywkECxOjSJyUqS/Xyht\nqpRaP7324fHjR86yfe970rhxkR/vy8TGxqqyslKdnZ0qLCzU0qVLtXr1anV3d3/l2ri4OFVWVqq9\nvV15eXlasmTJNWu7u7u1evVqLVmyRHl5edf9lowxISlJqquTDh2SPvpo5PRpdbU0MHDttn/6b+ml\nz6Wt9dK35kqTEqTsv5W2vCwFXpX+672Ijw83EUJEToz07X+Vvvup5H1eOjcGv70nOTlZdXV1am1t\nVVdXlzIyMlRdXa2B6+3or5KamvrF1z11dnYqIyND69ev1/r165WRkaHOzk61tLTI7/crNTU1An9N\nGM2dO3Ld0O+XXnpJKiiQ9u27bAOTDu+WJi+X5k+7cu3UBdIjM6Td/ysNR3JouIoQIjRnpIrskSO2\nSz/xC6Tmj4Nb7vkryfdr6S8vSr/6HykQ3mnDJi8vT01NTXrllVe0detWLViwQNu3bw9qbX5+vvbv\n36+NGzdq27Zt2rZtm55//nnt379f+fn5YZ48gsaNk1atktrbpe98R+rsvOzBIelsjzQlVZpw9cJ4\n6c5p0qdnpMHIjQt3EUKEZqa0uWPkjYGXfvrbpMUpwT/FHX8j1a2WtvxEauoJ26QRsWLFCh05ckTl\n5eVas2aNiouL1dbW9pXrPB6PHn/8cXV2dqqzs1Nr166Vx+OJwMRRMGmS9POfS2vWXPbLWClpmnT+\n7HVi1y99+Kk09Q5pjJ0ZxthECDG6BqWdVdI3EqRx8dKyf5K6rj5rOF568J+l75v0k19KH43x018T\nJ05UVVWV2tvblZ6ervvuu09r167Vxx9/9WGyx+P5+gbwS42TcldIvb+T3j135UOftUk7P5K+fZ/k\n4n8NIo4QYlT1HZT+7W3J3y71HpbmvSn9uuXa7cYnS1V1UsKb0nufRHzMsLh0DbC5uVltbW3Kzs5W\nXV2dhofHeOnDJfO7Uvk06fv/KDV2SBf6pc69UvlqyVZKq+ZFe0I4ghBiVE25X/r9Xulb35Cm/bWU\nnymlTLr+tklLpE0/kxJus3eI3qpL1wCfe+45bdiwQQ888IB6e3ujPdbtZ9x0qeq30rNzpJ8ukSYn\nSAufkFJ+Ku1+TprB7gmRMc4smNtAACEKSO+8KP3qrPSbf5HucPQUV29vr1577TWtXbs22qMAuAFC\niFFnf5F++wvp96nSL5+Qkh2NIICxIUJ3b4Qzzkv/USl1PiJteuQ674wHgNsMR4QYVV2vSAtXSWcu\n+92qN6T/fDhqIwHAlyKEAACn8bYsAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAA\npxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEII\nAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4j\nhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADg\nNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggB\nAE4jhAAAp/0fNQEfvedJNR0AAAAASUVORK5CYII=\n",
       "prompt_number": 52,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x27c60c0>"
       ]
      }
     ],
     "prompt_number": 52
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol = Chem.MolFromSmiles('[NH2+]=CC=C[O-]')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 53
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "mol"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAYAAABNcIgQAAAN2UlEQVR4nO3de0yVd57H8Y9yKxYU\nAWtL28ELoCJtZMRYO146U2edKttON8v+0Q2kbfTUzDYHd+yUbrQFs7GQpruDOzNrz9RtB+umiSYz\nCbq6s+pEu6LYXVsbqSKKkeiKVPFO53A5fPePE6nWSzlyOcDv/UpIKjy/c772D988v+c55wwzMxMA\nAI4aHu4BAAAIJ0IIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQA\nOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFC\nAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAa\nIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAA\npxFCAIDTCCEAwGmEEADgNEIIAHAaIQQAOI0QAgCcRggBAE4jhAAApxFCAIDTCCEAwGmEEADgNEII\nAHAaIQQAOI0QAgCcRgjhnDNnzoR7BAADCCGEM5qamrRkyRKlpaURQwBdCCGGvGvXrumNN97QhAkT\ndPnyZX3xxRdKSUkJ91gABojIcA8A9BUz00cffaQVK1YoISFBlZWVevrpp286JhAIKCIiIkwTDgyB\nQECSnP//AHdxRoghqbq6WrNmzdLy5cv15ptv6vPPP78lgocOHVJGRoY2bdoUpinDb+PGjUpPT9ek\nSZNUWloqv98f7pGAfkcIMaQ0NTWpoKBAc+fO1cyZM3X06FF5PB5FRt66+TF58mR5vV4tXrxYOTk5\nqqqqCsPE4VFVVaWcnBwtWbJEhYWFevfdd7Vu3TpNnDhR69evl5mFe0Sg/xgwBPj9fisuLrb777/f\nnnnmGTty5Ei31545c8Y8Ho9FRkZaXl6eNTQ09OGk4XXy5EnLy8uzyMhI83g8dubMma6ftba2Wnl5\nuY0cOdJmzJhh+/btC+OkQP/hjBCD3qZNm5SVlaUNGzbo448/1tatWzV58uRurd22bZvWrl0rn8+n\n6upqNTY2asqUKSopKRlS24R+v18lJSXKzMxUY2Oj9u/fL5/Pp4ceeqjrmOjoaBUWFqq2tlbZ2dma\nM2eOCgoK1NTUFMbJgX4Q7hID9+r48XP23HPPWWxsrL311lvW0tLS7bW1tbW2cOFCi46Ottdee806\nOzvNzKyzs9M2btxoqamp9uijj1pFRUXXzwajzs5Oq6iosEceecRSU1OtsrKy22sPHDhgs2fPtri4\nOCsuLja/39+HkwLhQwgx6Fy5Yvb662ZxcW3m8SwPaSuzqanJ8vPzbfjw4Zafn2+nTp267XEtLS1W\nXFxssbGx9tRTT9nBgwd7a/x+c/DgQZs3b57FxcVZWVlZSL8o3KiystLGjRtn6enptnnz5l6eEgg/\nQohBo73drLzcbPRos9mzzQ4c6P5av99vZWVllpCQYNOmTbNdu3Z1a92pU6duCmdTU9M9Tt9/uhv7\nULS0tFhZWZnFxcXZ/PnzraamphcmBQYGQohBYetWs8mTzcaNMwthd8/MzDZv3mwZGRk2duxYq6io\nsEAgEPLz79q1y6ZNm2YJCQlWVlZmra2tIT9GX7sx9jNnzrTq6upef47Tp09bfn6+RUdHm9frtUuX\nLvX6cwD9jRBiQDt+3Cw31ywuzqyszCyU3b2amhqbP3++xcTEWHFxsV29erVHswQCAauoqLCxY8da\nRkaGbdmypUeP15t6EvudO3eGvOW5Y8cOe+yxx+yDRYvMNmwwG8TXUQFCiAHp8mUzr9csOtosP9/s\n9Onur7148aJ5vV6LioqyvLw8q6ur69XZrl69asXFxRYTExP2bcKexL6+vt6ef/55i4yMtFWrVoX8\n3O3t7fZnn88sKcls1iyz/ftDfgxgICCEGFA6OoLXAZOTzZ54IrR/Wzs6Oqy8vNwSExMtKyvLduzY\n0XeDmlldXZ3l5uZaVFRUv28TXrp06Z5jf+7cOfN4PF1rjx071rNhrl0zKy42i4kJnr6fONGzxwP6\nGSHEgPHJJ2bf/77Zgw+aVVSEtttWXW2Wn7/J4uPjrbS0tF9v9d++fbtNnTrVkpKSrLy83Do6Ovrs\nua7HPikpybKysmz79u0hr01MTLTHH3/cdu7c2bvDHTtmlpdnFhtrVlQUvL0XGAQIIcLu4kWzl18O\nboO++qpZc3P3154+Hdw6jYgwe/nlDmtsPNtnc95NW1ub+Xw+S05OtuzsbNu9e3evP8fu3bstOzvb\nkpOTzefzWVtbW7fXXr/ZZ8yYMebz+ay9vb3X5+uyY4dZVpZZSoqZz2d2DzcnAf2JECLs/P5gCEO5\n1HblyjfXEH/yE7MQ3lGtTzU3N5vX67WIiAjLzc21E72wTXjixAnLzc21iIgI83q91hzCbwrX10ZH\nR1tRUZFduHChx/N0S1tbcI87IcEsJ8dsz57+eV7gHhBCDCqBQPAkIyXFbOLE0F9K0V8OHz5sCxYs\nsNjYWCsqKrIr97BNeOXKFSsqKrLY2FhbsGCBHT58+J7WLlq0yGpra0N+/l7R3Bz8jSUyMrhtevJk\neOYA7oIQYtDYsyd4cjFyZPBkI4SdwbCprKy08ePH28MPP9ztt2u7/rZoKSkpNmHChJDeFu3GtZmZ\nmSFdQ+xTn31mNneu2YgRwRtrvv463BMBXQgh7uir/zLL/Tuzc2Geo6EheDIRGWnm8Zjd8IEJg8L1\nT3WIj4+3GTNm2N69e+947N69e23GjBkWHx9v5eXlIb1w//raUaNGWXl5eUjXEPtFIGD24YfBu6Ey\nM4N74sAAwKdPoNvajkjz0qX3627+fuuX0uwJ0r8eufPa9jrpR4nS7FLpyg0fddd2RJqXIX1Yf+sa\nv18qKZGmTJEaG6Xqasnnk274wIRB4U6f6nD27NmuY86ePauCggLNmTNH2dnZqq2tVWFhoaKjo7/z\n8a+vnTdvnmbNmqX6+noVFhYqKiqqL/9aoRs+XHrxRamuTlq9WoqJCX6/44L0watSRqI0LEIaM1l6\n9X3pXHtYx4U7CCFucblKeixKeuAvpC2/kcYMk6Yuly719IEDUtXb0r/8j3S3j301k9avl9LTpd/9\nLvj1ySfS9Ok9HSC8UlJSuj7uqb6+XmlpaVq5cqVWrlyptLQ01dfXq7q6Wj6fTykpKd/5eK2trSop\nKVFaWpoaGhq0f/9+rVmzRklJSf3wt+mB+Hjppz8N/rddlt75K+nt09KavdK1q9KecunsO9KildL5\nznBOCleE+5QUA9e3t0ZbD5vNTTP77dGbj/PXmP1gvNlv7nIvR9tRsx9OMHv+abMHnjDbd+mGx0w3\n++D4N8f+4Q/Bu0F/8YvgO8wMRR0dHbZ27VpLTU211NRUe++990J6/eHGjRstLS3Nxo8fH9I1xAGn\nbp3ZhOlmVd96M4LL+8ymf89szQGzf37STDJTqtm2gf+m5xh8IsMdYjgkUlq4SnpoieT9J+mPq6T7\nb3NYbq50+LA0cWK/T9hvIiIitHTpUi1ZsqTrz91RU1OjZcuW6dNPP9WKFSvk9XoVGxvbl6P2IZMO\nbZXifiRNHXXzj0ZmSc+OkbbWSP9RJf19eCaEG9gaxR2N+bG0+ddS8o3fbJI8k6Rhw775ui9Lqjrf\nvceMeFAq/qX09XvSO/8t3W7jKzJyaEfwRhEREd2K4Pnz5/XKK69o+vTpmjhxompra1VUVDSIIyhJ\n7dJXzVJ8ihTz7Z/dJz08SrrcJLWFYza4hBAiNGOl3x4NXse7/uWvkX6Q/N1Lr3vgx9KaAumDZdLu\n5j6bdEgIBAJas2aNJk2apC+//FJVVVXdvoY48EVJiaOklq9uEzu/9H+XpZEPSN99vxDQI4QQ/W+4\n9MN/kF40aVmpdC4Q7oEGposXL+rJJ5/U6tWrVVpaqt27dysnJyfcY/WiYdK0v5QubZO+uHjzj67U\nSJvPSQunS93bNQbuGSFE72qTNhdJ3xshDbtPmrdcami99bDhSVLRGmnEn6QjF/p/zMFg9OjReuml\nl1RXVyePx9Pt64iDSvpfS4tHSS++JP3nUenPfql+p7S4QLK/kV7IDPeEcAAhRK+6elD69f9Kvlrp\n0iEp80/SL6tvf2ziHOlXr0sjhvXriIPK0qVLlZCQEO4x+s6wBKno99Jb46Wfz5HiRkgzfyYl/1za\n+rY0hn+i0PeGmdndXtIF3Lt26d/+Vmp8XVo5lHb0AAwp/LqFvtEpffa+tD1T8mSHexgAuDPOCNHr\n7Gvp9/8o/TFFKv2ZlDQEL20BGDp4QT16V4u0rlCqf1b61bO3eXkYAAwwnBGiVzV8LM18QWq64Xsv\nbJH+fVHYRgKAuyKEAACncbMMAMBphBAA4DRCCABwGiEEADiNEAIAnEYIAQBOI4QAAKcRQgCA0wgh\nAMBphBAA4DRCCABwGiEEADiNEAIAnEYIAQBOI4QAAKcRQgCA0wghAMBphBAA4DRCCABwGiEEADiN\nEAIAnEYIAQBOI4QAAKcRQgCA0wghAMBphBAA4DRCCABwGiEEADiNEAIAnEYIAQBOI4QAAKcRQgCA\n0wghAMBphBAA4DRCCABwGiEEADiNEAIAnEYIAQBOI4QAAKcRQgCA0wghAMBphBAA4DRCCABwGiEE\nADiNEAIAnEYIAQBOI4QAAKcRQgCA0wghAMBphBAA4DRCCABwGiEEADiNEAIAnPb/PV3CgIZIDzcA\nAAAASUVORK5CYII=\n",
       "prompt_number": 54,
       "text": [
        "<rdkit.Chem.rdchem.Mol at 0x27c6210>"
       ]
      }
     ],
     "prompt_number": 54
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# In principle both of the above resonance structure should match a 'C=[NH2+]' query on (0,1)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}