Colors inverted with stats_summary in ggplot in R

I want the order of the colours in the graph to be the same as the order in the legend. In the legend, the green is on top, but it is at the bottom in the graph. How can I fix this? Below is my code and my graph. Thank you in advance!

asv_for_graph_taxon %>%
  mutate(., presence=case_when(total_counts>0~'Present',
                               total_counts==0~'Absent',
                               .default = 'WTF'))%>%
  filter(., sampling_point_simple!='make_up_water')%>%
  mutate(., taxon = factor(taxon, levels = c('*Nevskia* (ASV0000002)', 
                                             'Unclassified *Hyphomicrobiales* (ASV0000011)',
                                             '*Methyloversatilis* (ASV0000013)',
                                             '*Methylobacterium* (ASV0000016)',
                                             'Unclassified *Bacteria* (ASV0000021)',
                                             'Unclassified *Hyphomicrobiaceae* (ASV0000036)',
                                             '*Acidovorax* (ASV0000046)',
                                             'Unclassified *Betaproteobacteria* (ASV0000072)',
                                             'Unclassified *Pseudomonadota* (ASV0000088)',
                                             '*Sphingobium* (ASV0000106)',
                                             '*Rubribacterium* (ASV0000107)',
                                             'Unclassified *Betaproteobacteria* (ASV0000127)',
                                             'Unclassified *Bacteria* (ASV0000171)',
                                             '*Novosphingobium* (ASV0000179)',
                                             '*Mesorhizobium* (ASV0000192)',
                                             '*Novosphingobium* (ASV0000193)',
                                             'Unclassified *Bacteria* (ASV0000214)',
                                             'Unclassified *Bacteria* (ASV0000277)',
                                             'Unclassified *Alphaproteobacteria* (ASV0000313)')))%>%
  filter(., asv%in%asvs_maaslin3_lp_logistic)%>%
  ggplot(., aes(x = lp_qpcr_2, y = presence, color = taxon))+
  stat_summary(fun.data = median_hilow, geom = 'pointrange', 
               fun.args = list(conf.int = 0.5),
               position = position_dodge(width = 0.8))+
  scale_color_manual(name = 'ASV',
                     values = as.vector(p21))+
  labs(x = 'Log(GU/L)',y='ASV Presence')+
  scale_x_continuous(breaks = scales::pretty_breaks(n = 10))+
  theme_classic()+
  theme(text = element_text(size = 12, family = 'Times New Roman'),
        legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.key.height = unit(0.4, 'cm'), #change legend key height
        legend.key.width = unit(0.1, 'cm'), #change legend key width
        legend.spacing.x = unit(0.1, "cm"),
        legend.spacing.y = unit(0.1, "cm"),
        legend.position = 'right',
        legend.title = element_text(size = 9),
        legend.text = element_markdown(size = 9))

stat_summary graph where colors do not match between legend and graph

I'm not sure if I'm understanding this correctly. Can't reproduce it as no data was provided. Not sure if guide_legend(reverse = TRUE) is an answer you want.

library(tidyverse)

data <- tibble(x = c(1,2,3), y =c(4,5,6), class = c('a', 'b', 'c'))

ggplot(data = data, aes(x = x, y = y, color = class))+
  geom_point()

library(tidyverse)

data <- tibble(x = c(1,2,3), y =c(4,5,6), class = c('a', 'b', 'c'))

ggplot(data = data, aes(x = x, y = y, color = class))+
  geom_point()+
  guides(color = guide_legend(reverse = TRUE))