Experimentally Cr doping in the rutile phase of VO$_2$ is found to stabilize a charge ordered ferromagnetic insulating state in the doping range of 10% to 20%. In this work, we investigated its origin at 12.5% Cr doping using a combination of ab-initio electronic structure calculations as well as microscopic modeling. Our calculations are found to reproduce the ferromagnetic insulating state as well as a charge ordering at the V and Cr sites. The mapping of the ab-initio band structure onto a tight-binding Hamiltonian allows one to calculate the energy gain from different exchange pathways. This gain is quantified in this work for the first time and the role of charge ordering in stabilizing a ferromagnetic insulating state is understood.