The COVID-19 pandemic worldwide has resulted in over 176 million cases and roughly 3.8 million deaths so far. We could analyze mutation dynamics across the genome from countries such as the USA, Italy, the UK, France, Brazil, and India considering the rapid mutations of the SARS-CoV-2 genome. The analysis would help us to understand the genome diversity, the implications of the mutations in protein stability, and viral transmission. Among the 11 genes, surface glycoprotein (S) was singled out because of its crucial function associated with the entry of virion into the human cell upon binding with the hACE2 receptor. 749 S protein sequences from India were retrieved from the NCBI database for our study. The S protein is an important antigenic component responsible for inducing host immune responses, neutralizing antibodies, and providing protective immunity against viral infection. During an epitope prediction from a mutation-prone S-protein region, it is necessary to ascertain how new mutations significantly change the S protein, such that our vaccine is effective against all the mutated strains as well. The S1 region of the S protein had been our prime focus for identifying immune epitopes against SARS-COV-2. Antigenic B- cell epitopes were YYPDKVF from NTD and LFRKSNLKP from RBD. Cytotoxic T-cell epitopes WTAGAAAYY (within NTD) and CVADYSVLY (within RBD) exhibited binding with a maximum number of MHC I alleles. The T-cell epitopes which showed a maximum affinity for MHC II alleles were FLPFFSNVT within NTD and YFPLQSYGF within RBD. Furthermore, the best epitopes were characterized in terms of their physicochemical properties to establish their potentiality.