Run VEP on the provided variants
runVEP(variants, plotDir, cpus = 1, genome = "hg19", vepCall = "vep",
forceRedoVEP = F, rareGermline = T)
variants: The variants to be VEPed.
character: The plot directory where outputSomatics have been run.
integer: The number of cpus to be used as most.
Logical: if VEP should be rerun even if saved data already exists.
This function calls VEP on the output from outputSomaticVariants. For this, VEP needs to be callable by system('vep').