  We study effects of a bounded and compactly supported perturbation on
multi-dimensional continuum random Schr\"odinger operators in the region of
complete localisation. Our main emphasis is on Anderson orthogonality for
random Schr\"odinger operators. Among others, we prove that Anderson
orthogonality does occur for Fermi energies in the region of complete
localisation with a non-zero probability. This partially confirms recent
non-rigorous findings [V. Khemani et al., Nature Phys. 11, 560-565 (2015)]. The
spectral shift function plays an important role in our analysis of Anderson
orthogonality. We identify it with the index of the corresponding pair of
spectral projections and explore the consequences thereof. All our results rely
on the main technical estimate of this paper which guarantees separate
exponential decay of the disorder-averaged Schatten $p$-norm of $\chi_{a}(f(H)
- f(H^{\tau})) \chi_{b}$ in $a$ and $b$. Here, $H^{\tau}$ is a perturbation of
the random Schr\"odinger operator $H$, $\chi_{a}$ is the multiplication
operator corresponding to the indicator function of a unit cube centred about
$a\in\mathbb{R}^{d}$, and $f$ is in a suitable class of functions of bounded
variation with distributional derivative supported in the region of complete
localisation for $H$.
