Tautomer canonicalization bug

[driesvr]

Hi,

I have noticed some erratic behaviour in the tautomer canonicalization procedure. In this specific example, double bonds jump from one ring to the next, erasing stereocenters as they go.

from rdkit import Chem
from rdkit.Chem.Draw import IPythonConsole, MolsToGridImage
import molvs


#MWE failing molecule
mol = Chem.MolFromSmiles("NC(=O)c1[nH]nc(c1O)[C@@H]2O[C@H](CO)[C@@H](O)[C@H]2O")

smi_std = molvs.Standardizer()
standardized_tautomer = smi_std.tautomer_parent(mol)
mols = [mol, standardized_tautomer]

MolsToGridImage(mols)

Output attached for reference.

Cheers,

Dries