For every lattice \(\Lambda\), we construct functions which are arbitrarily close to the Gaussian, do not agree up to global phase but have Gabor transform magnitudes agreeing on \(\Lambda\). Additionally, we prove that the Gaussian can be uniquely recovered (up to global phase) in \(L^2(\mathbb{R})\) from Gabor magnitude measurements on a sufficiently fine lattice. These two facts give evidence for the existence of functions which break uniqueness from samples without affecting stability. We prove that a uniform bound on the local Lipschitz constant of the signals is not sufficient to restore uniqueness in sampled Gabor phase retrieval and more restrictive a priori knowledge of the functions is necessary. With this, we show that there is no direct connection between uniqueness from samples and stability in Gabor phase retrieval. Finally, we provide an intuitive argument about the connection between directions of instability in phase retrieval and Laplacian eigenfunctions associated to small eigenvalues.