Plot box plot with anova model and different p value

I have created box plot which looks like this

Code -

df = pd.read_csv("file.csv")
import plotly.express as px
fig = px.box(df, x="Treatment", y="Fresh weight (g)", color='Treatment', facet_col='Crop', title='Fresh Weight', height=750)
fig.show()

I want to use anova model with different p value. Any idea how to do this ? It should be something like this

Data can be seen here