#' mapping quality
#'
#' @param bamfile path to bam file
#' @export
map_qua <- function(bamfile){
p_mapq <- ScanBamParam(what=scanBamWhat(),flag=scanBamFlag())
mapq = scanBam(bamfile, param=p_mapq)[[1]]$mapq
return(mapq)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.