The computer code for the GARD model used in this paper is based on the work by Markovitch et al , and is available at github.com/ModelingOriginsofLife/GARD. The model is simulated with a Monte Carlo approach using a stochastic Gillespie algorithm (83), as described . The code was written and implemented in Matlab R2020b, whereby the original framework was modified to give specific stepwise data, including temporal compositions and fluxes, available at github.com/Amitmiti/GARD-model. Unless stated otherwise, the simulations run for 30 micellar generations with a repertoire of NG=100 lipid types and maximal size Nmax=100. All post-analyses were conducted in Python version 3.6, and the code is available at the latter GitHub link.