We consider Benham's model for strand separation in negatively supercoiled circular DNA, and study denaturation as function of the superhelical density k < 0. We propose a statistical version of this model, based on bayesian segmentation methods of current use in bioinformatics; this leads to new algorithms with priors adapted to supercoiled DNA, taking into account the random nature of the free energies needed to denature hydrogen bonds.