R script to read .VCF files