Tools for extracting reads by coverage

See https://github.com/ctb/khmer-recipes/blob/master/001-extract-reads-by-coverage/index.rst

Why only look at the coverages of the last 1000 seqs? Instead of continuously `sum`ing you could use a counter to store the total coverage and divide by the number of reads so far. More accurate, faster, and less memory used.