Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

sample command silently produces incorrect empty reads #213

Open
rickymagner opened this issue May 7, 2024 · 1 comment
Open

sample command silently produces incorrect empty reads #213

rickymagner opened this issue May 7, 2024 · 1 comment

Comments

@rickymagner
Copy link

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
Copy link
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
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant