seqtk icon indicating copy to clipboard operation
seqtk copied to clipboard

`sample` command silently produces incorrect empty reads

Open rickymagner opened this issue 1 year ago • 1 comments

Hi, I have a fastq with some empty reads like so:

@name

+

@other
AAA
+
999

Running seqtk sample bad.fastq 0.9 on this produces:

>name

@other
AAA
+
999

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.

rickymagner avatar May 07 '24 17:05 rickymagner

Just saw this issue here but can confirm this still happens on the latest version (r130).

rickymagner avatar May 07 '24 17:05 rickymagner