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

[feature request] Option to save (un)classified reads #30

Open
telatin opened this issue Dec 13, 2021 · 1 comment
Open

[feature request] Option to save (un)classified reads #30

telatin opened this issue Dec 13, 2021 · 1 comment

Comments

@telatin
Copy link
Contributor

telatin commented Dec 13, 2021

Kraken2 allows for --unclassified-out FILE (and --classified-out FILE) and this can be handy in contamination removal steps, and could benefit from the architecture of Sepia that could allow filtering multiple samples loading the database once.

@hcdenbakker
Copy link
Owner

Yes, that would be a nice feature! I have a read_filter subcommand that I am going to add to sepia that does something similar after a classification run. It is currently a stand alone Rust application living on my server. It takes the sepia outputs and the original read files and allows you to either include/exclude taxa (or unclassified reads) from a dataset and set a threshold for the kmer-similarity of the in-/excluded reads.

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

2 participants