draw_svg_bit_ID

from rdkit.Chem.Draw import rdMolDraw2D
from IPython.display import SVG
from rdkit import Chem
from rdkit.Chem import Draw
mol = Chem.MolFromSmiles(Chem.MolToSmiles(l))
beez = mol_to_bits(mol)

def mol_to_bits(mol):
    """Convert the atoms in a molecule into bit IDs. 
    This uses the morgan fingerprint with a depth of zero -
    i.e. it only considers 1-atom fragments. Each fragment gets
    a unique ID following the Morgan algorithm. 
    """
    bi = {}
    fp = AllChem.GetMorganFingerprint(mol, 0, bitInfo=bi,useFeatures=True)
    bits = np.zeros(mol.GetNumAtoms(),dtype=int)
    for key, value in bi.items():
        for v in value:
            bits[v[0]] = key
    return bits

d = rdMolDraw2D.MolDraw2DSVG(500, 500) # or MolDraw2DSVG to get SVGs
for i in range(mol.GetNumAtoms()):
    mol.GetAtomWithIdx(i).SetProp('atomNote', str(int(beez[i])))


d.DrawMolecule(mol)
d.FinishDrawing()
SVG(d.GetDrawingText())