How to facet correctly in ggline functon in r?


My data frame is like this data frame

I used ggline to graph the plot with add = c("mean_se", "jitter"). My plan was to separate the graph using the facet.by function inside the gland to separate my graph in the "Protein" column that corresponds to the "Peptide" column. However, the facet.by function that displays the x-axis label with all peptide labels, including the top panel, is labeled "Protein" for the peptide. I will show my codes here.

p_221_anova_plot = function(anova_plot_df) {
  anova_scatter = ggline(anova_plot_df, x = "Peptide", y = "Norm_Counts", color = "Group", shape = "Group", add = c("mean_se", "jitter"),
                         add.params = list(size = 0.7), font.label = list(color = "Protein"), repel = TRUE, point.size = 2, plot_type = "p",
                         error.plot = "errorbar", facet.by = "Protein") +
    theme(axis.text.x = element_text(angle = 90, size = 8, hjust = 1), plot.title = element_text(hjust = 0.5)) + 
    scale_color_manual(values = c("#440154FF", "#21908CFF")) +
    ylab("Counts normalized to beads only controls") + 
    xlab("Peptide") + 
    ggtitle("Mean comparison between VSG221 infected mice (3 infected mice) and mock mice (2 mock mice)") +
    stat_compare_means(aes(group = Group), label = "p.signif", method = "anova",
                       symnum.args = list(
                         cutpoints = c(0, 0.000001, 0.00001, 0.0001, 0.001, 1), symbols = c("****", "***", "**", "*", "ns")))

  return(anova_scatter)
}

p221_D0_plot = p_221_anova_plot(p221_D0_anova)
p221_D0_plot

I want to make such an argument. Desired plot
After using the facet.by function, the plot looks like this. Incorrect plot
Please help! Thank you!