Bam file to Hilber curve an then to R shiny

I would like to appear in the following graph to a framework made by R shiny. Where I must insert my piece of code? Thank you

library(chromstaR)
library(GenomicRanges)
library(HilbertCurve)
library(circlize)
bampath <-"C:/Users/Admin/Documents/R/win- 
library/3.6/GenomicRanges/extdata/wgEncodeUwRepliSeqBg02esG1bAlnRep1.bam"
bam <- bampath

reads <- readBamFileAsGRanges(bam, chromosomes='chr1', 
pairedEndReads=FALSE,
                    min.mapq=10, remove.duplicate.reads=TRUE)

hc = GenomicHilbertCurve(chr = "chr1", level = 5)
hc_points(hc, reads, gp = gpar(col = "red"))