Following recent experimental progress in the study of bottom baryons, we systematically calculate the mass spectra of Λb, Ξb, Σb, Ξb′, and Ωb baryons with a quark-diquark picture in the framework of a relativistic flux tube model with spin-dependent interactions in the j-j coupling scheme. Furthermore, we calculate the strong decay width of bottom baryons decaying into a bottom baryon and a light pseudoscalar meson. Good agreement is found between the calculated masses and the experimentally available masses of singly bottom baryons. By analyzing both mass spectra and strong decay widths, we interpret Σb(6097) as a 1P(3/2−) state and Ξb(6100) as a 1P(1/2−) state of the Ξb baryon. The Ξb(6227) is identified to be an orbital excitation 1P of the Ξb′ baryon with JP=3/2−. Further, we determine Ξb(6327) and Ξb(6333) as a 1P(3/2−) state and 1P(5/2−) state, respectively, of Ξb′ baryon. From the obtained mass spectra, we construct the Regge trajectories in the (J,M2) plane, which are found to be essentially linear, parallel, and equidistant. Our predictions for higher orbital and radial excited states can help experimentalists identify missing excited states of singly bottom baryons. Published by the American Physical Society 2024