In this article we complete the classification of the umbilical submanifolds of a Riemannian product of two space forms, addressing the case of a conformally flat product Hk×Sn−k+1, which has not been covered in previous works on the subject. We show that there exists precisely a p-parameter family of congruence classes of umbilical submanifolds of Hk×Sn−k+1 with substantial codimension p, which we prove to be at most min{k+1,n−k+2}. We study more carefully the cases of codimensions one and two and exhibit, respectively, a one-parameter family and a two-parameter family (together with three extra one-parameter families) that contain precisely one representative of each congruence class of such submanifolds. In particular, this yields another proof of the classification of all (congruence classes of) umbilical submanifolds of Sn×R, and provides a similar classification for the case of Hn×R. We determine all possible topological types, actually, diffeomorphism types, of a complete umbilical submanifold of Hk×Sn−k+1. We also show that umbilical submanifolds of the product model of Hk×Sn−k+1 can be regarded as rotational submanifolds in a suitable sense, and explicitly describe their profile curves when k=n. As a consequence of our investigations, we prove that every conformal diffeomorphism of Hk×Sn−k+1 onto itself is an isometry.