We calculate the quasi-normal mode complex frequencies of the Kerr-Newman black hole with arbitrary values of spin and charge, for the modes typically dominant during a binary black hole coalescence, ( , m, n) = {(2, 2, 0), (2, 2, 1), (3, 3, 0)}. Building analytical fits of the black hole spectrum, we construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant U (1) charge. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound q < 0.33 on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values q 0.5; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to q ∼ 0.3, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission.