scripts to go from fastq to vcf files for whole genome alignment