The equation \$\mathrm{LSB} = \frac{\mathrm{FS}}{2^N - 1}\$ is commonly repeated, but is not in fact correct. For an ADC with \$N\$ bit resolution, each bit must therefore represent a span of \$\frac{1}{2^N}\$ of the input range. As such the equation \$\mathrm{LSB} = \frac{\mathrm{FS}}{2^N}\$ is correct.
So the first input will span from \$\left[0 \space \frac{1}{2^N}\right)\$, the second \$\left[\frac{1}{2^N}\space\frac{2}{2^N}\right)\$, and so fourth.
What this means is that the voltage you are measuring could be anywhere within that range, so you must decide how you want to represent it - for example you could floor and pick the bottom of the range (1), ceil and pick the top of the range (2), or pick the mid value (3). That would give calculations of:
$$\begin{align}\\
V_{in} &= \frac{V_{ref}\times \mathrm{Code}}{2^N}\tag1\\\\
V_{in} &= \frac{V_{ref}\times (\mathrm{Code}+1)}{2^N}\tag2\\\\
V_{in} &= \frac{V_{ref}\times (\mathrm{Code}+0.5)}{2^N}\tag3\\
\end{align}$$
To add a real example to the mix, take the ADC of the ATMega1284. If you refer to page 253 of the datasheet, the equation relating voltage to ADC code is given. Rearranging that into your form, it is:
$$V_{in} = \frac{V_{ref}\times \mathrm{Code}}{2^N}$$