Adding static visualizations on top of a slider

Hello,

I am giving a tutorial on MLE in R and would love to be able to show my normal distributions moving over 3 points on the X axis. However, I am only able to get the points to appear in the first active frame and am not sure how to adjust for this.

x <- seq(0, 10, length.out = 1000)
df_temp <- data.frame(x = c(-1,0,1), y = c(0,0,0))
aval <- list()
for (step in 1:6) {
  aval[[step]] <- list(
    visible = FALSE,
    name = paste0('v = ', step),
    x = x,
    y = dnorm(x, step+1)
  )
  
}

aval[3][[1]]$visible = TRUE

steps <- list()
fig <- plot_ly()

for (i in 1:6) {
  fig <-
    add_lines(
      fig,
      x = aval[i][[1]]$x,
      y = aval[i][[1]]$y,
      visible = aval[i][[1]]$visible,
      
      name = aval[i][[1]]$name,
      type = 'scatter',
      mode = 'lines',
      hoverinfo = 'name',
      
      line = list(color = '00CED1'),
      showlegend = FALSE
    )

step <- list(args = list('visible', rep(FALSE, length(aval))),method = 'restyle')
  step$args[[2]][i] = TRUE
  steps[[i]] = step
  
}


fig <- fig %>% add_markers(x = c(4.5,5,5.5), y = c(0,0,0))
# add slider control to plot

fig <- fig %>%
  layout(sliders = list(list(
    active = 0,
    currentvalue = list(prefix = "Frequency: "),
    steps = steps
  )))


fig
If you run this code you get the three points I would like to graph on the first "slide plot" however they disappear after and I am not sure how to keep them for every plot

I would appreciate any help.