Custom Code

filter.py selects the best hit for each marker protein in each genome based on the highest scoring alignment region. SubsampleAlignmentRadomly.py randomly deletes 50% of alignment columns. phylobayes_convergence_statistics.r calculates Rhat and ESS values for molecular clock outputs.