This research aims to comprehensively estimate seismicity parameters and create seismotectonic and isoseismal maps for the state of Meghalaya, located in the North-East (NE) region of India using an earthquake catalog collected between 1861 and 2022. The seismic influence zone considered for the study is a circular area with a radius of 350 km around Shillong city (latitude 25.57° N and longitude 91.88° E). Distinct magnitude classes, such as 3.50–3.99, 4.00–4.99, 5.00–5.99, 6.00–6.99, and ≥ 7.00, are considered for analysis, and corresponding completeness periods are obtained as 40, 60, 100, 160, and 130 years, respectively. The study area seismicity parameters a and b values of the Gutenberg–Richter (G–R) relationship are obtained with values ranging from 5.50 to 6.23 and 0.70 to 0.75, respectively. This study further estimates earthquake magnitude occurrence probabilities for 1, 50, and 100