In the rapid shape change from spherical to deformed nuclei in the Z=40 Zr isotopic chain, recent work has identified shape coexistence in Zr96. Between Zr96 and the strongly deformed Zr100, Zr98 is expected to also exhibit coexistence of nuclear shapes. The degree of mixing between different configurations is mainly determined by the nucleon-nucleon interactions. For nuclear model predictions, experimental constraints are needed, but they are barely available for Zr98. To study low-lying transitions in Zr98, a Coulomb excitation experiment was conducted at the Argonne Tandem-Linac Accelerator System (ATLAS) facility using a Zr98 beam extracted from the Californium Rare Ion Breeder Upgrade (CARIBU) ion source and Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA) for γ-ray spectroscopy coupled to the compact heavy ion counter (CHICO2) for ion detection. This paper reports on the first decisive deduction of the B(E2;21+→01+) transition strength in Zr98 and on its interpretation