Index reference genome

[Adamtaranto]

The current method loads a reference genome into memory and then uses string operations to extract sequences that match a hmm query. This is fine for most microbial genomes but can become clunky with larger sequences.

For the next major release - replace all sequence operations with pyfaidx which creates an index of the fast file and only loads minimal amounts of the seq into memory as required.