This paper proposes, validates and analyzes the dynamics of the susceptible exposed infectious recovered (SEIR) model for the propagation of COVID-19 in Saudi Arabia, which recorded the largest number of cases in the Arab world. The model incorporates a saturated incidence rate, a constant vaccination rate and a nonlinear treatment function. The rate of treatment is assumed to be proportional to the number of infected persons when this number is low and reaches a fixed value for large number of infected individuals. The expression of the basic reproduction number is derived, and the model basic stability properties are studied. We show that when the basic reproduction number is less than one the model can predict both a Hopf and backward bifurcations. Simulations are also provided to fit the model to COVID-19 data in Saudi Arabia and to study the effects of the parameters of the treatment function and vaccination rate on disease control.