library("ape")
library("Biostrings")
library("ggplot2")
library("ggtree")

setwd("")

nwk <- ("tree.newick")
tree <- read.tree(nwk)

tipcategories = read.csv("tree.meta",
                         sep = "\t",
                         col.names = c("seq", "cat"),
                         header = FALSE,
                         stringsAsFactors = FALSE)

dd = as.data.frame(tipcategories)


p <- ggtree(tree, layout = "rectangular", right = FALSE, ladderize = TRUE) +
  xlim(0, 0.025) + # to allow more space for labels
  geom_treescale() # adds the scale

p %<+% dd +
  geom_tiplab(aes(fill=factor(cat)),
              color = "black", # color for labels
              geom = "label",  # labels not text
              label.padding = unit(0.15, "lines"),
              label.size = 0) +
  theme(legend.position=c(0.5,0.2),
        legend.title = element_blank(), # no title
        legend.key=element_blank()) # no keys