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

Kraken2 seems to classify reads from read pairs separately #768

Open
SimonHegele opened this issue Oct 25, 2023 · 1 comment
Open

Kraken2 seems to classify reads from read pairs separately #768

SimonHegele opened this issue Oct 25, 2023 · 1 comment

Comments

@SimonHegele
Copy link

Tried to classify samples from paired-end sequencing using the command:

kraken2 --db {database} --paired {query} {mate} --report {report} --threads {threads} > {out} 2> {log}

My samples had 60-80 million read pairs, however Kraken2 reported 120-160 million fragments (two times as much).
Apparently each read is counted. At the same time same taxa have uneven numbers of associated fragments.
This indicates that the reads of a read pair could be mapped to different taxa.

@jenniferlu717
Copy link
Collaborator

I believe your command is incorrect. It should be

kraken2 --db {database} --report {report} --threads {threads} --paired {query} {mate} > {out} 2> {log}

Read files should be given last. and do not need to be right after the --paired argument

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