Shigellosis is the main cause of food and waterborne diarrhea and is an emerging threat to human health. The current study characterized the indigenous multidrug-resistant Shigella flexneri serotypes for their plasmid profiles and genetic diversity, to characterize the plasmid evolutionary patterns and distribution. In total, 199 identified S. flexneri isolates belonging to six different serotypes were analyzed for plasmid profiling, followed by an analysis of whole genome sequencing. All isolates of S. flexneri resistant to antibiotics harbored multiple copies of plasmids with sizes ranging from 1.25 kbp to 9.4kbp. These isolates were clustered into 22 distinct plasmid patterns, labeled as p1-p22. Among these, p1 (24%) and p10 (13%) were the predominant plasmid profiles. All S. flexneri strains were grouped into 12 clades with a 75% similarity level. Also, a significant association was observed among the plasmid patterns, p23 and p17 with the drug-resistant patterns AMC, SXT, C (19.5%) and OFX, AMC, NA, CIP (13.5%), respectively. Moreover, the most widespread plasmid patterns p4, p10, and p1 showed a significant association with the serotypes 1b (29.16%), 2b (36%), and 7a (100%), respectively. After plasmid sequence assembly and annotation analysis, a variety of small plasmids that vary in size from 973 to 6200bp were discovered. Many of these plasmids displayed high homology and coverage with plasmids from non-S. flexneri. Several novel plasmids of small size were discovered in multidrug-resistant S. flexneri. The data also showed that plasmid profile analysis is more consistent than antibiotic susceptibility pattern analysis for identifying epidemic strains of S. flexneri isolated in Pakistan.