In this study, the statistics properties of the aftershock sequence and coulomb stress changes due to the mainshock of the 𝑀𝑤= 6.3, 9 April 2013 Shonbeh earthquake have been studied. The spatial distribution of aftershocks and focal mechanism event with 𝑀𝑤≥5.0 declared a NW–SE striking rupture with oblique thrust fault mechanism. The b-value is estimated as 0.77± 0.05 and the mainshock occurred in a highly stressed region and sequence comprised larger magnitude aftershocks due to the presence of large size asperities within the rupture zone. Also the high p-value (1.16 ± 0.05), implying fast decay rate of aftershocks, evidences high surface heat flux. A value of the spatial fractal dimension (𝐷𝑐) equal to 2.60 ± 0.02 indicates random spatial distribution and source in a two-dimensional plane that is being filled-up by fractures. Moreover, the evaluated coseismic Coulomb stress using the variable slip model for depth range 0–20 km reveals a ‘butterfly’ pattern and most of the aftershocks fall (95%) in the region of enhanced Coulomb stress.