The speed of sound of the matter within neutron stars may contain non-smooth structure related to first-or higher-order phase transitions. Here we investigate what are the observable consequences of structure in the speed of sound, such as bumps, spikes, step functions, plateaus, and kinks. One of the main consequences is the possibility of ultra-heavy neutron stars (with masses larger than 2.5 solar masses) and mass twins in heavy (with masses larger than 2 solar masses) and ultra-heavy neutron stars. These stars pass all observational and theoretical constraints, including those imposed by recent LIGO/Virgo gravitational-wave observations and NICER X-ray observations. We thoroughly investigate other consequences of this structure in the speed of sound to develop an understanding of how non-smooth features affect astrophysical observables, such as stellar radii, tidal deformability, moment of inertia, and Love number. Our results have important implications for future gravitational wave and X-ray observations of neutron stars and their impact in nuclear astrophysics.