Skip to content

sample command silently produces incorrect empty reads #213

Open

Description

@rickymagner

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.

Activity

  1. rickymagner commented on May 7, 2024

    @rickymagner
    Author

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions