In this study, variation in the characteristics of the pressure oscillations in a thermoacoustic engine is explored as the input heat flux is varied. A bifurcation diagram is plotted to study the variation in the qualitative behavior of the acoustic oscillations as the input heat flux changes. At a critical input heat flux (60 Watt), the engine begins to produce acoustic oscillations in its fundamental longitudinal mode. As the input heat flux is increased, incommensurate frequencies appear in the power spectrum. The simultaneous presence of incommensurate frequencies results in quasiperiodic oscillations. On further increase of heat flux, the fundamental mode disappears and second mode oscillations are observed. These bifurcations in the characteristics of the pressure oscillations are the result of nonlinear interaction between multiple modes present in the thermoacoustic engine. Hysteresis in the bifurcation diagram suggests that the bifurcation is subcritical. Further, the qualitative analysis of different dynamic regimes is performed using nonlinear time series analysis. The physical reason for the observed nonlinear behavior is discussed. Suggestions to avert the variations in qualitative behavior of the pressure oscillations in thermoacoustic engines are also provided.