Probability of finding single particle state in QM and QFT I was reading "Quantum Optics" by Walls and Milburn, and in chapter 10.1 describing about atom-radiation interaction, it says normal QM equation for single electron in an atom (here $\psi_j(x)$ is an energy eigenstates),
$$\psi(x)=\sum_j a_j \psi_j (x)  $$
could be quantized to a field operator,
$$\hat{\psi}^\dagger (x)=\sum _j \hat{a_j} ^\dagger \psi^*_j (x) $$
Here $\hat{a_j}^\dagger$ is a creation operator for Fock state. 
For normal QM equation, probability for energy state $\psi _j (x)$ being selected is $|a_j|^2$, by Born rule, and probability that this state measured to have specific position $x$ is $|\psi_j (x)|^2$. Therefore probability that $\psi(x)$ is being measured to have energy state $j$ and position $x$ should be $|a_j|^2|\psi_j(x)|^2$.
Now for the case of QFT, field operator $\hat{\psi}^\dagger(x)$ acts on vacuum state $\left|0\right>$, then 
$$\hat{\psi} ^\dagger (x) \left|0\right>= \sum _j  \psi^*_j (x) \hat{a_j} ^\dagger \left|0\right> = \sum _j  \psi^*_j (x) \left|1_j \right>$$
Now I believe $\left|1_j \right>$ is a state having energy same as QM state $\psi_j(x)$ and having position exactly $x$. Since we are talking about single electron, I believe $\left|1_j \right>$ is same as $\psi_j(x)$ with fixed position $x$. But instead of probability we calculated above, $\left|1_j \right>$ state has probability of only $|\psi_j(x)|^2$. 
What am I missing here?
EDIT: After some search, and thanks to the answer by Charles Francis, I have found the solution to my confusion. 
The whole confusion was from the notation. Here, $\psi(x)$ is a wavefunction of single electron state, and I thought $\hat{\psi} ^\dagger (x)$ is an operator generating this wavefunction $\psi(x)$ due to the resemblance of the notation, which is incorrect! Actually textbook has unrigorous approach to quantizing process.
From this lecture note or any other reference about second qunatization one could find basis transformation of creation and annihilation operator, given as 
$$\hat{a}^\dagger _{\tilde{\alpha}} = \sum_\alpha \left< \alpha | \tilde{\alpha}  \right> \hat{a}^\dagger _\alpha $$
Here $\{ \left|\alpha\right> \}$ and $\{ \left|\tilde{\alpha}\right> \}$ is some basis of Hilbert space. Putting basis $\{ \left|x\right> \}$ instead of $\{ \left|\tilde{\alpha}\right> \}$ gives
$$\hat{a}^\dagger _{x} = \sum_\alpha \left< \alpha | x  \right> \hat{a}^\dagger _\alpha= \sum_\alpha \psi_\alpha^* (x) \hat{a}^\dagger_\alpha  $$
One can now notice that field operator $\hat{\psi} ^\dagger (x)$ was actually an operator generating one particle at position $x$. Therefore $\hat{\psi} ^\dagger (x) \left|0\right>=\left|x\right>$, not $\left|\psi\right>$.
And as Charles Francis mentioned, creation operator for $\left|\psi\right>$ analogous to the equation $\psi(x)=\sum_j a_j \psi_j (x)  $, needs to be defined independent to the definition of $\hat{\psi} ^\dagger (x)$, and that will be
$$\hat{\psi} \equiv \sum_j a_j \hat{a_j}^\dagger$$
Of course, $a_j$ is just an probability coefficient, and $\hat{a_j}^\dagger$ is an creation operator for energy state $j$.
 A: I have to say I think the notation is horrible. I would not use it this way but it is clear that
$$ \left|1_j \right> = \hat{a_j} ^\dagger \left|0\right>  $$
and that this is not related to $\psi$. $\psi_j (x)$ are the coefficients of $\hat{\psi}  (x)$ in the energy basis.
One cannot think of $\psi_j (x)$ as referring to definite energy and position, because that is forbidden by the uncertainty principle. This can only be treated as a function of position at particular energy.
Part of the problem is that it is common in applied maths to be ambiguous about whether $\psi (x)$ means the function $\psi$ or whether it means the value of that function at the point $x$. Ket notation can be good for removing that ambiguity.
It may help to write your first equation
$$ \psi(x) = \langle x|\psi \rangle = \sum_j \langle x| 1_j \rangle \langle 1_j|\psi\rangle = \sum_j a_j(x)\psi_j. $$
I think this shows that the argument is incorrectly placed. The energy eigenstates have wave functions $a_j(x)$, and the $\psi_j$ are coefficients in the energy basis. It can be written
$$ \psi = |\psi \rangle = \sum_j | 1_j \rangle \langle 1_j|\psi\rangle = \sum_j a_j\psi_j. $$
when you quantise it, replacing $|1_j\rangle$ with the creation operator and ignoring $|0\rangle$  this becomes
$$ \hat\psi^\dagger  = \sum_j {\hat a_j}^\dagger\psi_j $$
and since $\psi$ is arbitrary, we can replace it with $\psi(x)$ so as to construct a field operator
$$ \hat\psi^\dagger(x)  = \sum_j {\hat a_j}^\dagger\psi_j(x) $$
