download and index genome for alignment and variant calling
Side note that even though BWA could index .fa.gz, it will generate extra files that will end up taking up similar amount of disk space as .fa.
To format gff for bcftools csq, see notes here: samtools/bcftools#530 (comment)