Motivation: Genome-scale metabolic models are widely constructed and studied for understanding various design principles underlying metabolism, predominantly redundancy. Metabolic networks are highly redundant and it is possible to minimize the metabolic networks into smaller networks that retain the functionality of the original network. Results: Here, we establish a new method, MINREACT that systematically removes reactions from a given network to identify minimal reactome(s). We show that our method identifies smaller minimal reactomes than existing methods and also scales well to larger metabolic networks. Notably, our method exploits known aspects of network structure and redundancy to identify multiple minimal metabolic networks. We illustrate the utility of MINREACT by identifying multiple minimal networks for 77 organisms from the BiGG database. We show that these multiple minimal reactomes arise due to the presence of compensatory reactions/pathways. We further employed MINREACT for a case study to identify the minimal reactomes of different organisms in both glucose and xylose minimal environments. Identification of minimal reactomes of these different organisms elucidate that they exhibit varying levels of redundancy. A comparison of the minimal reactomes on glucose and xylose illustrates that the differences in the reactions required to sustain growth on either medium. Overall, our algorithm provides a rapid and reliable way to identify minimal subsets of reactions that are essential for survival, in a systematic manner. Availability and implementation: Algorithm is available from https://github.com/RamanLab/MinReact. © The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com