We study numerically the saddle point structure of two-dimensional (2D) lattice gauge theory, represented by the Gross-Witten-Wadia unitary matrix model. The saddle points are in general complex-valued, even though the original integration variables and action are real. We confirm the trans-series/instanton gas structure in the weak-coupling phase, and identify a new complex-saddle interpretation of non-perturbative effects in the strong-coupling phase. In both phases, eigenvalue tunneling refers to eigenvalues moving off the real interval, into the complex plane, and the weakto-strong coupling phase transition is driven by saddle condensation. Introduction: Path integral saddle points are physically important in quantum mechanics, matrix models, quantum field theory (QFT) and string theory, and are deeply related to the typical asymptotic nature of weak coupling perturbative expansions. Such relations are central to the concept of resurgence, whereby different saddles are intertwined by monodromy properties that connect them and account for Stokes phases. The theory of resurgence has recently provided new insights into matrix models and string theories [1][2][3][4][5][6], and has been applied to asymptotically free QFTs and sigma models [7][8][9][10], and localizable SUSY QFTs [11]. One motivation for such QFT studies is to find a practical numerical implementation of a semiclassical expansion that could provide a Picard-Lefschetz thimble decomposition of gauge theory, either in the continuum or on the lattice, especially for theories with a sign problem [12,13]. A unifying theme in these studies, and in related work [14][15][16][17], is the appreciation that complex saddles are important, in particular in the context of phase transitions, even though the original 'path integral' may be a sum over only real configurations.