Using IQTREE I received the most optimum models for my data to analyse phylogenetic tree-reconstruction. IQTREE is able to use, aside from the Gamma distribution for heterogeneity across sites also the FreeRate model:
"FreeRate model (Yang, 1995; Soubrier et al., 2012) that generalizes the +G model by relaxing the assumption of Gamma-distributed rates. The number of categories can be specified with e.g. +R6 (default 4 categories if not specified). The FreeRate model typically fits data better than the +G model and is recommended for analysis of large data sets."
Does anyone know if it is possible to use this FreeRate model in MrBayes (or in RevBayes), and if yes, how to I put that in the script?