We analyze a new numerical relativity data set of spinning but non-precessing binary black holes on eccentric orbits, with eccentricities from approximately 0.1 to 0.5. Spinning black holes with dimensionless spins of up to 0.75 are included at mass ratios q = m 1 /m 2 = (1, 2), and further non-spinning binaries at mass ratios q = (1.5, 3, 4). Comparison of the final mass and spin of these simulations with non-eccentric data extends previous results in the literature on circularisation of eccentric binaries to the spinning case. For the (l, m) = (2, 2) spherical harmonic mode we construct eccentric hybrid waveforms that connect the numerical relativity data to a post-Newtonian description for the inspiral, and we discuss the limitations in the current knowledge of post-Newtonian theory which complicate the generation of eccentric hybrid waveforms. We also perform a Bayesian parameter estimation study, quantifying the parameter biases introduced when using quasicircular waveform models to estimate the parameters of highly eccentric binary systems. We find that aligned-spin quasicircular models with higher order modes produce lower biases in certain parameters, like the mass ratio and the luminosity distance, than aligned-spin quasicircular models with only the (l, m) = (2, ±2) modes and precessing models.