def _plot_ids(ax, x, y, size, marker='.', color=None, **kwargs):
    """Plot points, return the used color"""
    # plot is faster than scatter, go figure...
    return ax.plot(x, y, color=color, markersize=size, marker=marker,
                   linestyle='None', **kwargs)[0].get_color()