ENH: output fastq or fasta