vcf <- read.vcf(system.file("extdata/clinvar_dbSNP138_example.vcf.gz", package = "bedr"));
vcf$header <- c(vcf$header, NOTE="vcf processed by bedr")
if (FALSE) {
write.vcf(vcf, filename = paste(tempdir(), "/bedr.example.vcf", sep = ""));
}
Run the code above in your browser using DataLab