Hi, I have a fastq with some empty reads like so:
Running seqtk sample bad.fastq 0.9 on this produces:
which is no longer valid fastq, as the @name read has changed format. This breaks lots of downstream tools. Is it possible to keep it in the original format? I'd like to preserve these if sampled because there might be a read pair which is nonempty.
Hi, I have a fastq with some empty reads like so:
Running
seqtk sample bad.fastq 0.9on this produces:which is no longer valid fastq, as the
@nameread has changed format. This breaks lots of downstream tools. Is it possible to keep it in the original format? I'd like to preserve these if sampled because there might be a read pair which is nonempty.