powered by
write.vcf(x, filename = NULL, verbose = TRUE)
vcf <- read.vcf(system.file("extdata/clinvar_dbSNP138_example.vcf.gz", package = "bedr")); vcf$header <- c(vcf$header, NOTE="vcf processed by bedr") # write.vcf(vcf, filename = "bedr.example.vcf");
Run the code above in your browser using DataLab