We study the quench dynamics of an Anderson impurity model using the configuration interaction (CI) method. In particular, we focus on the relaxation behavior of the impurity occupation. The system is found to behave very differently in the weak-coupling and strong-coupling regimes. In the weak-coupling regime, the impurity occupation relaxes to a time-independent constant quickly after only a few oscillations. In the strong-coupling regime, the impurity occupation develops a fast oscillation, with a much slower relaxation. We show that it is the multi-peak structure in the manybody energy spectrum that separates these two regimes. The characteristic behavior, including the power-law decay and the period of oscillation, can also be related to certain features in the manybody energy spectrum. The respective advantages of several impurity solvers are discussed, and the convergence of different CI truncation schemes is provided.