We present spectral and temporal properties of all the thermonuclear X-ray bursts observed from Aql X-1 by the Neutron Star Interior and Composition Explorer (NICER) between 2017 July and 2021 April. This is the first systematic investigation of a large sample of type I X-ray bursts from Aql X-1 with improved sensitivity at low energies. We detect 22 X-ray bursts including two short recurrence burst events in which the separation was only 451 s and 496 s. We perform time resolved spectroscopy of the bursts using the fixed and scaled background (f a method) approaches. We show that the use of a scaling factor to the pre-burst emission is the statistically preferred model in about 68% of all the spectra compared to the fixed background approach. Typically the f a values are clustered around 1-3, but can reach up to 11 in a burst where photospheric radius expansion is observed. Such f a values indicate a very significant increase in the pre-burst emission especially at around the peak flux moments of the bursts. We show that the use of the f a factor alters the best fit spectral parameters of the burst emission. Finally, we employed a reflection model instead of scaling the pre-burst emission. We show that reflection models also do fit the spectra and improve the goodness of the fits. In all cases we see that the disc is highly ionized by the burst emission and the fraction of the reprocessed emission to the incident burst flux is typically clustered around 20%.