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!