Hi everyone!

I'm running Kraken2 on different metagenomes (Illumina paired end). I'm currently using it this way:

kraken2 --output Sample_1 --paired Sample_1_R1.fastq  Sample_1_R2.fastq

With this I get everything running ok and get a single file with both classified and unclassified sequences. I would like to keep just the classified ones so I've tried:

kraken2 --classified-out Sample_1_classified --paired cseqs#.fq Sample_1_R1.fastq  Sample_1_R2.fastq

but it Kraken tells me

--paired requires positive and even number filenames

I'm not really sure what this means since I'm giving it two files.

On an unrelated matter but taking advantage of the question. I would like to run this analysis on several samples which should be easy doing a bash script but I haven't been able to find a way to merge the different outputs into one table to the use for example bracken or pavian on a single file. Does anyone now a way to do this?

Thanks a lot!!!

EDIT: I've come up with a not so elegant workaround:

awk -F "\t" '$1=="C"' Sample_1> Sample_1_Classified

I would still like to know if there's a better solution to this. Again, thanks a lot!