DNA supercoiling, the under or overwinding of DNA, is a key physical mechanism both participating to compaction of bacterial genomes and making genomic sequences adopt various structural forms. DNA supercoiling may lead to the formation of braided superstructures (plectonemes), or it may locally destabilize canonical B-DNA to generate denaturation bubbles, left-handed Z-DNA and other functional alternative forms. Prediction of the relative fraction of these structures has been limited because of a lack of predictive polymer models that can capture the multiscale properties of long DNA molecules. In this work, we address this issue by extending the self-avoiding rod-like chain model of DNA so that every site of the chain is allocated with an additional structural degree of freedom reflecting variations of DNA forms. Efficient simulations of the model reveal its relevancy to capture multiscale properties of long chains (here up to 21 kb) as reported in magnetic tweezers experiments. Well-controlled approximations further lead to accurate analytical estimations of thermodynamic properties in the high force regime, providing, in combination with experiments, a simple, yet powerful framework to infer physical parameters describing alternative forms. In this regard, using simulated data, we find that extension curves at forces above 2 pN may lead, alone, to erroneous parameter estimations as a consequence of an underdetermination problem. We thus revisit published data in light of these findings and discuss the relevancy of previously proposed sets of parameters for both denatured and left-handed DNA forms. Altogether, our work paves the way for a scalable quantitative model of bacterial DNA.