allGroupsColors<- c("pink", "orange", "green", "cyan", "blue", "purple", "seagreen", "steelblue", "red", "DarkTurquoise", "Sienna", "Chocolate", "BlueViolet", "Magenta", "brown", "gray", "darkred", "darkmagenta", "mediumvioletred")

p=plot_ordination(physeq.rarefy, physeq.rarefy_bray, shape ="Region", color = "Villages")
p = p + geom_point(size = 4, alpha = 0.75)
p = p + scale_color_manual(values = allGroupsColors)
p + ggtitle("PCoA root 16S data, Bray distance")