A. Neumaier said:
so that I can see the complete argument in the degenerate case. I don't ask questions if they are completely trivial!
OK, let me answer your question about the degenerate case, if I understood your question correctly.
A general measurement is described by a POVM. However, any POVM measurement can be viewed as a projective measurement in some higher-dimensional Hilbert space. So without loosing on generality, we can assume that we have a projective measurement defined by some projectors ##\pi_k## satisfying ##\sum_k \pi_k =1##. These projectors can be further decomposed as
$$\pi_k = \sum_l |k,l\rangle\langle k,l |$$
where the set of all ##|k,l\rangle## make a complete basis in the measured system. Such a measurement distinguishes different ##k## but not different ##l##, so ##l## can be said to label the degeneracy of the measured observable.
Before measurement, let the state of the measured system be
$$|\psi\rangle = \sum_{k ,l}c_{kl}|k,l\rangle$$
During the measurement, the measured system will become entangled with the apparatus and its environment, so the total state will be
$$|\Psi\rangle = \sum_{k ,l}c_{kl}|k,l\rangle |A_{k}\rangle |E_{k}\rangle$$
where ##|A_{k}\rangle## are macroscopically distinguishable states of the apparatus and ##|E_{k}\rangle## are macroscopically distinguishable states of the environment. Let me further introduce the notation ##|A_{k}\rangle |E_{k}\rangle \equiv |\Phi_{k}\rangle##, so that the equation above can be written as
$$|\Psi\rangle = \sum_{k ,l}c_{kl}|k,l\rangle |\Phi_{k}\rangle$$
This can also be written in terms of wave functions as
$$\Psi(x,y) = \sum_{k ,l}c_{kl}\psi_{kl}(x) \Phi_{k}(y)$$
where ##x## is the collective coordinate for particle positions of the measured system, while ##y## is the collective coordinate for particle positions of the apparatus+environment. For simplicity, I assume that there are no spin degrees.
Since ##\Phi_{k}(y)## are macroscopically distinguishable, we have an (approximate) equation
$$\Phi^*_{k'}(y) \Phi_{k}(y)=0 \;\; for \;\; k\neq k' \;\; (Eq. 1)$$
The probability density is
$$P(x,y) =|\Psi(x,y)|^2 \;\; (Eq. 2)$$
which is valid in both standard and Bohmian QM. Using (Eq. 1) this becomes
$$P(x,y) = \sum_{k}|\Phi_{k}(y)|^2 |\psi_{k}(x)|^2$$
where
$$\psi_{k}(x)=\sum_{l}c_{kl}\psi_{kl}(x)$$
The probability in the ##y##-space is given by average over ##x##
$$P(y) =\int dx P(x,y) = \sum_{k} |\Phi_{k}(y)|^2 \sum_{l}|c_{kl}|^2$$
where we have used orthogonality and standard normalization of the basis ##\psi_{kl}(x)##.
Now let ##\sigma_k## be the region in the ##y##-space in which ##\Phi_{k}(y)## is not negligible. Then the probability of finding the apparatus+environment inside ##\sigma_k## is
$$P_k=\int_{\sigma_k}dy P(y) = \sum_{l}|c_{kl}|^2$$
where we used standard normalization of ##\Phi_{k}(y)##. The obtained result ##P_k=\sum_{l}|c_{kl}|^2## is nothing but the Born rule in the degenerate case, so we have proved the Born rule in the ##k##-space with degeneracy from the Born rule (Eq. 2) in the position space. Q.E.D.
Note that the whole proof does not depend much on Bohmian mechanics, so in this sense there is nothing really new in it and can be considered "trivial" for an expert in standard QM, even he/she does not know much about Bohmian mechanics.