The Nicolis—Puhl reaction, an isothermal, two-component, autocatalytic reaction exhibiting bistability, oscillatory solutons, and Takens-Bogdanov points in a perfectly micromixed CSTR has been shown to be sensitive to the degree of micromixing under non-premixed feed conditions. These properties make it an interesting theoretical model for exploring the effect of micromixing on the dynamics of nonlinear chemical reactions; in particular, for examining whether micromixing effects suffice to produce complex dynamics in an otherwise two-dimensional dynamical system. Here, the reaction is studied in conjunction with the interaction-by-exchange-with-the mean (IEM) model of micromixing. Numerical bifurcation and stability analysis techniques applicable to the specific class of integro-differential equations to which this model belongs are employed to carry out the analysis. For various values of the reaction rates, the micromixing rate and the inverse residence time have been used as bifurcation parameters to construct two-dimensional bifurcation sets containing saddle-node and Hopf bifurcation curves whose locations are shown to be strongly dependent on the degree of micromixing. As is observed experimentally by lowering the stirring rate for some real chemical reactions in flow reactors, regions of bistability and oscillatory behavior are shown to be accessible, by decreasing the micromixing rate, from regions exhibiting a single steady state. Manifolds of Takens-Bogdanov points for finite micromixing rates have also been computed; however, it is shown that the dynamical dimension of the IEM model remains the same for all values of the micromixing parameter.