3372 and 3376 had been identified as the highest number of core genes. There was only a small distinction within the estimated core between the two options. The default Roary pairwise id threshold is just too stringent for such a diverse dataset, and this is doubtless due to gene clusters being wrongly cut up into…