ref.fa is the reference used for the mapping.
The output BAM file was processed with mpileup.
We filtered the output to keep the equivalent of 1000 reads.
