Subjet distributions in deep inelastic scattering at HERA