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