The structures and band gaps of copper-zinc-tin selenosulfides (CZTSSe) are investigated for a range of anion compositions through experimental analysis and complementary first-principles simulations. The band gap was found to be extremely sensitive to the Sn-anion bond length, with an almost linear correlation with the average Sn-anion bond length in the mixed anion phase Cu2ZnSn(SxSe1-x)4. Therefore, an accurate prediction of band gaps using first-principles methods requires the accurate reproduction of the experimental bond lengths. This is challenging for many widely used approaches that are suitable for large supercells. The HSE06 functional was found to predict the structure and band gap in good agreement with the experiment but is computationally expensive for large supercells. It was shown that a geometry optimization with the MS2 meta-GGA functional followed by a single point calculation of electronic properties using HSE06 is a reasonable compromise for modeling larger supercells that are often unavoidable in the study of point and extended defects.