# the file of VCF
vcf.fn <- seqExampleFileName("vcf")
# or vcf.fn <- "C:/YourFolder/Your_VCF_File.vcf"
# convert
seqVCF2GDS(vcf.fn, "tmp.gds")
# display
(f <- seqOpen("tmp.gds", FALSE))
seqDelete(f, info.varname=c("HM2", "AA"), format.varname="DP")
f
# close the GDS file
seqClose(f)
# clean up the fragments, reduce the file size
cleanup.gds("tmp.gds")
# remove the temporary file
unlink("tmp.gds", force=TRUE)
Run the code above in your browser using DataLab