Full-shell nanowires are hybrid nanostructures consisting of a semiconducting core encapsulated in an epitaxial superconducting shell. When subject to an external magnetic flux, they exhibit the Little-Parks (LP) phenomenon of flux-modulated superconductivity, an effect connected to the physics of Abrikosov vortex lines in type-II superconductors. We show theoretically that full-shell nanowires can host subgap states that are a variant of the Caroli-de Gennes-Matricon (CdGM) states in vortices. These CdGM analogs are shell-induced Van Hove singularities in propagating core subbands. We elucidate their structure, parameter dependence, and behavior in tunneling spectroscopy through a series of models of growing complexity. Using microscopic numerical simulations, we show that CdGM analogs exhibit a characteristic skewness towards higher flux values inside nonzero LP lobes resulting from the interplay of three ingredients. First, the orbital coupling to the field shifts the energy of the CdGM analogs proportionally to the flux and to their generalized angular momentum. Second, CdGM analogs coalesce into degeneracy points at flux values for which their corresponding radial wavefunctions are threaded by an integer multiple of the flux quantum. And third, the average radii of all CdGM-analog wavefunctions inside the core are approximately equal for realistic parameters and are controlled by the electrostatic band bending at the core/shell interface. As the average radius moves away from the interface, the degeneracy points shift towards larger fluxes from the center of the LP lobes, causing the skewness. This analysis provides a transparent interpretation of the nanowire spectrum that allows to extract microscopic information by measuring the number and skewness of CdGM analogs. Moreover, it allows to derive an efficient Hamiltonian of the full-shell nanowire in terms of a modified hollow-core model at the average radius.