wgs_prep.pl from sra/
A good example for strain_bz.pl works with WGS
strain_bz.pl works with Ensembl