A new QCD analysis of Deep Inelastic Scattering (DIS) data is presented. All available neutrino and anti-neutrino cross sections are reanalysed and included in the fit, along with charged-lepton DIS and Drell-Yan data. A massive factorisation scheme is used to describe the charm component of the structure functions. Next-to-leading order parton distribution functions are provided. In particular, the strange sea density is determined with a higher accuracy with respect to other global fits.