vcf and density script