Skip to content


A Phylogenetics program design to estimate the joint posterior distribution of gene trees and species.

BEST requires an environment module

In order to use BEST, you must first load the appropriate environment module:

module load intel openmpi

BEST (Bayesian Estimation of Species Trees) is a bioinformatics program designed to analyze phylogenetic trees built out of multilocus molecular sequence data and estimate their joint posterior distributions in a Bayesian framework. This software is designed for use with the MrBayes phylogenetic tree building program. The program can account for deep coalescence but is not designed to address horizontal gene transfer or gene duplication issues among others.

Using BEST on RCC Resources#

For Interactive Jobs with BEST#

BEST can be run in a couple of different modes. In order to start an interactive session in BEST, all you need to do is load the necessary modules and call BEST from the command line:

$ module load intel openmpi
$ best

For Non-Interactive Jobs with BEST#

If you want to run BEST non-interactively with a genome file, you can run it on the HPC login nodes by running the following commands, where EXT is a bioinformatics file format like Nexus (.nex) and FILE is the name of the data file. Example data files are available on the official website.

$ module load intel openmpi
$ best FILE.EXT

Alternatively, you can run BEST in parallel using the Slurm scheduler on multiple nodes. Below is an example Slurm script to run BEST. The script must be saved as a file with the .sh suffix. Again, EXT is a bioinformatics file format like Nexus (.nex) and FILE is the name of the data file.

#SBATCH -J BestJob   # Can rename the job whatever you want
#SBATCH -n 16        # Can use more or less processes if you want
#SBATCH -p genacc_q  # Can use any queues you have access to
#SBATCH -t 1:00:00   # Can request more time in HH:MM:SS format if needed
module load intel openmpi
srun best FILE.EXT

Then submit your script using the following command, replacing YOURSCRIPT with the name of your script file:

$ sbatch