To provide context for my question, I'll simplify it to consider three symmetric particles moving in one dimension. In my actual case, I'm dealing with $N$ generic particles in three dimensions.
Let $| \psi \rangle$ represent a state composed of three bosonic particles moving in one dimension, which can be expressed as follows: $$ | \psi \rangle = \int dx dy dz \ \phi(x,y,z) \ | x, y, z \rangle, $$ where $\phi(x,y,z)$ is the wave function of the state and $$ | x, y, z \rangle = \hat{a}^\dagger(x) \hat{a}^\dagger(y) \hat{a}^\dagger(z) \ |0 \rangle, $$ with the creation (or annihilation) operator given by the following commutation relations: $$ [\hat{a}(x), \hat{a}^\dagger(y)] = \delta (x-y), \\ [\hat{a}^\dagger(x), \hat{a}^\dagger(y)] = [\hat{a}(x), \hat{a}(y)] = 0. $$
My question is the following:
Using the commutation relations, the orthogonal relation of the basis $| x, y, z \rangle$, should be given by: $$ \langle x', y', z' | x, y, z \rangle = \delta (x-x') \delta (y-y') \delta (z-z') + \rm{permutations}. $$ However, I have seen in some papers the orthogonal relation of the Fock state given by $$ \langle x', y', z' | x, y, z \rangle = \delta (x-x') \delta (y-y') \delta (z-z'), $$ that is, without the permutations. I suspect that this simplification may be justified when the particles are distinguishable and independent from each other (in which case, the wave function can be factorized as $\phi(x,y,z) = \phi_1(x) \phi_2(y) \phi_3(z)$). However, I haven't been able to find any references or notes that specifically address this point.
Could someone please provide me with relevant bibliography or a better clarification on this matter?