We study regularity of solutions of weakly singular Volterra integral equations of the first kind. We then study the numerical analysis of discontinuous piecewise polynomial collocation methods for solving such systems. The main purpose of this paper is the derivation of global convergent and super-convergent properties of introduced methods on the graded meshes. We apply relevant methods to a system of fractional differential equations and analyze them. The numerical experiments confirm the theoretical results.