Flight Paths in Leaflet

Question

I'd like to plot headquarter locations as circles and lines connecting the headquarter locations to their branches in leaflet. The data is structured exactly in the same way as the flight route example for plotly (see example:

https://plot.ly/r/lines-on-maps/ ).

I need an interactive map that shows the geography (rivers, streets, towns, etc) of Kansas, parts of Nebraska, parts of Colorado, and Parts of Missouri.I'd like for popups to appear describing the locations. So I think Leaflet might be better than plotly for this task.

This what I currentltly have for plotly:

This is my code for the plotly map:

library(plotly)
library(dplyr)


# HQ locations
HQ <- read.csv('HQ.csv', stringsAsFactors=FALSE, fileEncoding="UTF-8-BOM")
j <- c(17:18)
HQ[ , j] <- apply(HQ[ , j], 2,            # Specify own function within apply
                       function(x) as.numeric(as.character(x)))

# Storage Locations from HQ
Storage <- read.csv('Storage.csv', stringsAsFactors=FALSE, fileEncoding="UTF-8-BOM")

#COLUMN NUMBER

Storage$id <- seq_len(nrow(Storage))
i <- c(16:18)
Storage[ , i] <- apply(Storage[ , i], 2,            # Specify own function within apply
                       function(x) as.numeric(as.character(x)))

colnames(HQ)

# map projection
geo <- list(
  scope = 'usa',
  projection = list(
    type = 'albers usa'),
  showland = TRUE,
  landcolor = toRGB("lightgray"),
  showlaKes = T,
  lakecolor = toRGB('white'),
  lonaxis = list(range = c(-108, -95)),
  lataxis = list(range = c(37, 42)))


p <- plot_geo(locationmode = 'USA-states', color = I("purple")) %>%
  add_markers(
    data = HQ, x = ~Lng_HQ, y = ~Lat_HQ, text = ~paste('</br><b> Co-op Name:</b> ',

                                                       Cooperative,
                                                        '</br><b> Co-op HQ Location:</b> ', Headquarter,
                                                        '</br><b> Grain Capacity (bu):</b> ', Grain.Capacity..in.bushels.,
                                                        '</br><b> Grain Storage Locations:</b> ', CNGL
                                                        ),
    size = ~CNGL, hoverinfo = "text", alpha = 0.5
  ) %>%
  add_segments(
    data = group_by(Storage, id),
    x = ~Lng_HQ, xend = ~Lng,
    y = ~Lat_HQ, yend = ~Lat,
    alpha = 0.3, size = I(1), hoverinfo = "text", text =~paste('</br><b> Co-op Name:</b> ',

                                                               Cooperative,
                                                               '</br><b> Co-op Branch Location:</b>', Storage$Branch,
                                                               '</br><b> Grain Capacity (bu):</b> ', Storage$Grain.Capacity..in.bushels.)) %>%
  layout(
    title = ' <b>Map of All Cooperatives with Grain Locations</b><br>(Hover for Cooperative & Grain Storage Details)',
    geo = geo, showlegend = F, autosize = F, height=1000, width=1000, margin = m 
  )

p