The curated SNP alignment was then used to compute Neighbor Joining (NJ) and Maximum Parsimony (MP) trees with MEGA6 [54] and a Maximum Likelihood (ML) tree using PhyML 3.0 [55] with the GTR model used in previous Y. pestis work [7, 63], with 4 gamma categories and the best of ...