#' cigar profile
#'
#' @param bamfile path to bam file
cigar_prof <- function(bamfile){
p_cig <- ScanBamParam(what=scanBamWhat(),flag=scanBamFlag(isUnmappedQuery = FALSE))
cig = scanBam(bamfile, param=p_cig)[[1]]$cigar
return(cig)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.