Since you mention DNA seqs, you might want to take a look at [nim-hts](https://github.com/brentp/hts-nim/). It's easy on unix systems to get htslib, and not too difficult on windows using msys2
You could probably write an iterator around his BGZ type.
