Overview
Transitions to other unoccupied states are possible but extremely unlikely, more likely that the photon will not be absorbed.
Introduction
The Pauli exclusion principle prevents a third electron occupying the $2s$ state. Even if there was space in the $2s$ state a $1s\to 2s$ transition is unlikely due to selection rules and a $1s\to2p$ transition is significantly more likely if there is space in the $2p$ orbital.
Other answers here have stated that the transition to other energy levels is forbidden. Now while the probability of a transition is extremely small, it is non-zero.
A quick note on notation: I will be using a bold typeface for vectors as opposed to an over arrow so that vector operators are clearer.
Quantisation of the Electromagnetic Field
Minimal coupling of the electron to the electromagnetic field using the Coulomb potential adds a perturbation of:
$$\hat H_1=\frac{e}{m_e}\hat{\boldsymbol p}\cdot\hat{\boldsymbol A}\left(\boldsymbol r,t\right)$$
to the Hamiltonian. Where $e$ and $m_e$ are the charge and mass of the electron, $\hat{\boldsymbol p}$ is the momentum operator acting on the electron and the vector potential operator has the form:
$$\hat{\boldsymbol A}\left(\boldsymbol r,t\right)=\sum_{\lambda,\boldsymbol k}\sqrt{\frac{\hbar}{2v\epsilon_0\omega\left(\boldsymbol k\right)}}\left(\hat a_\lambda\left(\boldsymbol k\right)\boldsymbol s_\lambda\left(\boldsymbol k\right)e^{i\left(\boldsymbol {k}\cdot\boldsymbol r-\omega t\right)}+\text{h.c.}\right)$$
where $\text{h.c.}$ is the Hermitian conjugate of the preceding terms, $v$ is the volume of the cavity in which the experiment is taking place; $\omega\left(\boldsymbol k\right)$ is the angular frequency of the photon mode as a function of the wavevector $\boldsymbol k$; $\lambda$ labels the two polarisations; $\boldsymbol s_\lambda\left(\boldsymbol k\right)$ is the polarisation vector of the mode; $\hat a_\lambda\left(\boldsymbol k\right)$ is the annihilation operator for the mode; and $\boldsymbol r$ is the position of the atom (assuming the wavelength is larger than the atom the uncertainty in the electron's position can be ignored).
If we have a single wavelength and polarisation then:
$$\hat H_1=\frac{e}{m_e}\sqrt{\frac{\hbar}{2v\epsilon_0\omega}}\hat{\boldsymbol {p}}\cdot\boldsymbol s\hat a e^{i\left(\boldsymbol k\cdot\boldsymbol r-\omega t\right)}+\text{h.c.}$$
Thus, let:
$$\begin{align}\hat V&=\frac{e}{m_e}\sqrt{\frac{\hbar}{2v\epsilon_0\omega}}\hat{\boldsymbol {p}}\cdot\boldsymbol s\hat a e^{i\boldsymbol k\cdot\boldsymbol r}\\\implies\hat H_1&=Ve^{-i\omega t}+\hat V^\dagger e^{i\omega t}\end{align}$$
Then using first-order time-dependent perturbation theory which holds in the limit $\frac{t}{\hbar}\left|\langle f|\hat V|i\rangle\right|\ll1$ for all $n\ge2$. We find the probability of a transition having occured if the atom is measured after a time $t$ since the electromagentic field was applied is:
$$\begin{align}P\left(t\right)=\frac{t^2}{\hbar^2}\Bigg|&\overbrace{e^{i\left(\Delta\omega-\omega\right)t/2}\operatorname{sinc}\left(\frac{1}{2}t\left(\Delta\omega-\omega\right)\right)\langle f|\hat V|i\rangle}^\text{absorption}\\+&\underbrace{e^{i\left(\Delta\omega+\omega\right)t/2}\operatorname{sinc}\left(\frac{1}{2}t\left(\Delta\omega+\omega\right)\right)\langle f|\hat V^\dagger|i\rangle}_\text{emission}\Bigg|^2\end{align}\tag{1}$$
where $\Delta E=\hbar \Delta \omega$ be the difference in the energy levels of the initial $|i\rangle$ and final $|f\rangle$ states. This is in general non-zero even when $\Delta \omega\ne\omega$. However, we can make one more approximation to aid in understanding: if the final state $|f\rangle$ has absorbed a photon then in the limit $t\Delta\omega\gg2\pi$ the $\operatorname{sinc}$ functions do not overlap and so we need only retain the absorption term:
$$P\left(t\right)=\left(\frac{\left|\langle f|\hat V|i\rangle\right|}{\hbar}\right)^2t^2\operatorname{sinc}^2\left(\frac{1}{2}t\left(\Delta\omega-\omega\right)\right)\tag{2}$$
Further approximations from here will give you Fermi's Golden rule, one of these approximations is taking the limit such that $t\operatorname{sinc}^2$ tends to a delta function and so removes the possibility for a transition when the energy of the photon is not exactly equal to the energy gap: and so this is an inappropriate approximation to make in this case.
Energy Conservation in Quantum Mechanics
While the expectation value of the energy is conserved in the evolution of a system as described by Schrödinger equation, there may be a discontinuous jump in the energy of the system when a measurement is performed. Consider a system in a superposition of energy eigenstates, when you measure the energy the state will collapse into an energy eigenstate which in general will not have the same energy as the expectation value for the energy - the energy of the system has increased or decreased!
The energy may be transferred to or from the measurement device or surroundings to compensate.
In previous edits this section also contained a discussion of the many-words type interpretation which in my naivety I included. I apologise for anyone I have mislead and for more details you can see this question:
"Conservation of energy, or lack thereof," in quantum mechanics
@Jagerber48's answer is the most relevant to this question giving additional details that will likely be of interest to any reader of this question.
@benrg's answer gives a good explanation of why energy is conserved.
@NiharKarve's comment includes a blog post which explains why the paper may be misleading.
Putting this all Together
Equation (1) shows, in general, the when an atom is illuminated by light of a single specific wavelength and polarisation, a transition is possible even if the energy of the photons is not equal to the energy gap, which would violate energy conservation (but this is allowed); however, the probability is extremely small.
Equation (2) makes a further approximation which we can now use to find an expression for the probability:
$$P\left(t\right)=\frac{e^2}{2v\epsilon_0m_e^2\hbar\omega}\left|\langle f|\hat{\boldsymbol {p}}\cdot\boldsymbol s\hat a |i\rangle\right|^2t^2\operatorname{sinc}^2\left(\frac{1}{2}t\left(\Delta\omega-\omega\right)\right)$$
As $|i\rangle\equiv|i\rangle_e|f\rangle_{EM}$ and $|f\rangle_e|f\rangle_{EM}$ where subscript $e$ is the electrons states and subscript $EM$ are the states of the electromagnetic field. Without detail $_e\langle f|\hat{\boldsymbol {p}}\cdot\boldsymbol s|i\rangle_e\equiv\boldsymbol d_{fi}\cdot\boldsymbol s$ where $\left\{\boldsymbol d_{fi}\right\}$ are the dipole matrix elements and are zero for transitions between certain orbitals independent of the energy supplied (for more details see selection rules). Finally, $_{EM}\langle f|\hat a|i\rangle_{EM}=\sqrt{N}$ if the state $|i\rangle_{EM}$ is the state for $N$ photons of the given wavelength and polarisation - but other states such as coherent states are also posible.
$$\implies P\left(t\right)=\frac{e^2N}{2v\epsilon_0m_e^2\hbar\omega}\left|\boldsymbol d_{fi}\cdot\boldsymbol s\right|^2t^2\operatorname{sinc}^2\left(\frac{1}{2}t\left(\Delta\omega-\omega\right)\right)\tag{3}$$
which holds in the limit:
$$t\ll\frac{\hbar}{\left|\,_e\langle f|\left(\hat{\boldsymbol {p}}\cdot\boldsymbol s\right)|i\rangle_e\right|}\sim 10^{-25}\text{s}$$
As the limit $t\Delta\omega\gg2\pi$ is not needed when the state $|i\rangle_{EM}$ is the state for $N$ photons of the given wavelength and polarisation because the creation operator causes the emission term to vanish anyway. However, the time is of the order of $10^{-25}\text{s}$ give or take a few orders of magnitude for Neon (obtained using the only data I could find for reduced matrix elements for dipole transitions), which is not a practical time scale to measure on.
Finally, considering your given case, given selection rules, the most likely case if the $1s$ electron did absorb a photon is a transition to the $3p$ state (as $2p$ is occupied and $3s$ is forbidden to first order by selection rules). Substituting values into equation (3) gives an order of magnitude estimate for the probability of transitioning from $1s$ to $3p$ in Neon of $10^{-12}\%\text{ per }\left(\text{photon }m^{-3}\right)$ for $t=10^{-25}\text{s}$ which is the point the approximation of first order perturbation breaks down.