code for plotting the phylogeny created in sim 3325 and highlighting two clades (2486 and 3157) that are highlighted in Figure 3