generic script to use bwa and samtools to pileup reads to a reference.
parses a blastn output file (in outfmt 3: query-aligned flat) mode.