Issue1506

Fixes #1506 

All `parse_fasta` calls have been removed from the script. All ways of filtering are now single pass over the sequence file. Minimal state in memory is now needed for filtering.
Can you improve the documentation? With a description of what `seqid_f` is?