The output can contain the original values when a properly scaled linear phase interpolation filter is used that has an odd number of coefficients. This is because the group delay for a linear phase filter is $(N-1)/2$ where $N$ is the number of coefficients in the filter. For the same reason, an even length filter will have a half sample offset (at the output rate) in every output sample, so therefore cannot include any of the input samples directly.
To realize this condition, the filter should be designed as a $1/M$-band filter where $M$ represents the interpolation factor. (So a half-band filter for interpolate by 2, quarter band filter for interpolate by 4, etc.) When designed in this fashion as an odd length filter, the center coefficient will be maximum, and then spaced off of the center coefficient, every $M$th coefficient will be zero. This means for the condition when the input sample is aligned with the center coefficient of the filter, all other coefficients will be zero'd due the zero-stuffing prior to the filter, and all other samples will be aligned with the filter zero coefficients. Thus only the center tap, as the input sample, will be passed to the output. (This also assumes the filter coefficients are scaled such that the center tap is 1).
This is illustrated with the determined coefficients for a quarter-band filter (for an interpolate by 4 application) shown below. Such filter coefficients are easily determined by centering the frequency targets for passband and stopband on $1/M$ (in this case 0.25) using Least-Squares design tools such as
firls in Matlab, Octave and Python scipy.signal (or Parks-McClellan, which is not recommended for high order interpolation and decimation filters due to the constant stop-band with consideration to aliasing).
coeff = 4 * sig.firls(31, [0, .24, .26, 1], [1,1,0,0])
Often, multi-band filters are used for interpolation (and decimation) filter designs which concentrate rejection only where needed (at the signal images for the case of interpolation). These can also be designed as an odd-length filter such that the center coefficient is $1$ and every $M$th coefficient is zero, resulting in the input sample being passed to the output without modification.