Working wrapper script

Computes the species summary and gene abundance, along with functional gene annotations. The SNP quantification is not currently implemented (mostly due to the extremely large amount of output it produces)