Generate summary statistics/plots

[hyphaltip]

Need to generate summary statistics about number of reads, number of read mapped, and perhaps average depth of coverage across the genome, number of unmapped reads too.