def _find_rct_position(rxn, mol):
    """Find the position of a reactant in a reaction"""
    react_pos = -1
    for pos, rct in enumerate(rxn.GetReactants()):
        if mol.HasSubstructMatch(rct):
            react_pos = pos
    return react_pos