Aims. We introduce a framework to measure the asphericity of Sun-like stars using a₁, a₂, and a₄ coefficients and constrain their latitudes of magnetic activity. Methods. We evaluated systematic errors on the inferred coefficients in function of key physical and seismic parameters (inclination of rotation axis, average rotation, height-to-noise ratio of peaks in power spectrum). The measured a-coefficients account for rotational oblateness and the effect of surface magnetic activity. We used a simple model that assumes a single latitudinal band of activity. Results. Using solar SOHO, VIRGO, and SPM data, we demonstrate the capability of the method to detect the mean active latitude and its intensity changes between 1999 and 2002 (maximum of activity) and 2006 and 2009 (minimum of activity). We further applied the method to study the solar-analogue stars 16 Cyg A and B using Kepler observations. In 16 Cyg A, we detected an equatorial band of activity exhibiting an intensity that could be comparable to that of the Sun. However, 16 Cyg B exhibits a bimodality in a₄ that is challenging to explain. We suggest that this could be a manifestation of the transition between a quiet and an active phase of activity. Validating or invalidating this hypothesis may require new observations.