messy.long <- gather(messy, taxa, count, -id, -trt)
head(messy.long)