Why do "most" (basically all) statistics text books use $\hat{\sigma}^2=\hat{p} (1-\hat{p})$ as an estimator for the variance of a Bernoulli process which we know is biased. Should the first choice not be the bias corrected sample version $s^2 = \hat{p} (1-\hat{p}) \cdot n/(n-1)$, which is the default for non-binary data.
I am aware of the multiple related threads that seem similar to my question. However, I am not seeking a discussion on the merits of unbiased estimators or the origin/relevance of the $n/(n-1)$ correction.
My specific question is about teaching and communicating, namely why for binary data alone we seem to prefer the biased variance ($1/n$) when for all other cases, the default sample variance is chosen to be the unbiased version ($1/(n-1)$)