file <- c(system.file("extdata/Strelka-ID-vcf",
"Strelka.ID.GRCh37.s1.vcf",
package = "ICAMS"))
if (requireNamespace("BSgenome.Hsapiens.1000genomes.hs37d5", quietly = TRUE)) {
catID <-
StrelkaIDVCFFilesToCatalogAndPlotToPdf(file, ref.genome = "hg19",
region = "genome",
output.file =
file.path(tempdir(), "StrelkaID"))}
Run the code above in your browser using DataLab