plot.bammshifts function

Plot distinct rate shift configurations on a phylogeny