We carry out a comprehensive survey of electromagnetic form factors of all light, heavy and heavy-light ground-state pseudoscalar and scalar mesons. Our analysis is based upon a Schwinger-Dyson equations treatment of a vector $\times$ vector contact interaction. It incorporates confinement and ensures axial vector and vector Ward-Takahashi identities are satisfied along with the corresponding corollaries such as the Goldberger-Treiman relations. The algebraic simplicity of the model allows us to compute the form factors at arbitrarily large virtualities of the probing photon momentum squared with relative ease. Wherever possible and insightful, we compare our results for the electromagnetic form factors and the charge radii with those obtained earlier through Schwinger-Dyson equations, lattice and with experimental observations available. We also comment on the scope and shortcomings of the model.