Hi team,

this is possibly related to issue #550 which was recently closed, where it was noted that empty sequences where not handled correctly.

I have a file where some sequences are empty:

>seq1

>seq2
MPFLG...

If I have a pattern file listing all the headers that I want to extract (e.g. seq1, seq2, ...), the empty sequences are not be reported in the resulting fasta. Ideally they would appear in the output (and the user can filter them out based on length) - or at least the grep command should issue a warning.

Is it possible to look into this issue?

Thanks