India had witnessed unprecedented surge in SARS-CoV-2 infections and its dire consequences during the second wave of COVID-19, but the detailed report of the epidemiological based spatiotemporal incidences of the disease is missing. Here in, we have applied various statistical methods like correlation, hierarchical clustering to know the pattern of pathogenesis of the circulating VoCs. B.1.617.1 (Kappa) was the predominant VoC during the early phase of second wave. Delta (B.1.617.2) or Delta-like (AY.x) VoC constitutes majority (>90.17) of the cases during the peak of second wave. The correlation plot showed Delta/Delta-like lineage is inversely correlated with other lineages including B.1.617.1 (kappa), B.1.1.7, B.1, B.1.36.29 and B.1.36. Delta/Delta-like surge coincided with second wave whereas all other lineages (B.1.617.1, B.1.36.29, etc.) occurred during the prior phase of the second wave. The spatiotemporal analysis showed that most of the Indian states were affected during the peak of the second wave due to delta surge and fall under the same cluster. The second cluster populated mostly by north-eastern states and islands of India were minimally affected. The presence of signature mutations (T478K, D950N, E156G) along with L452K, D614G and P681R within the spike protein of Delta or Delta-like might cause elevation in host cell attachment, increased transmission and altered antigenicity which in due course of time has replaced the other circulating variants.. The timely assessment of new VoCs will provide a rationale for updating the diagnostic, vaccine development by medical industries and decision making by various agencies including government, educational institutions, and corporate industries.