Does it necessary to use both forward and reverse fastq input files in kmergnie?
Entering edit mode
3 months ago
nakandalau • 0

For kmergenie, I used the following command

kmergenie hiseq_reads_ecoli_R1.fastq -l 21 -k 121 -s 2 -o out_file.txt -t 16

But If I want to use both forward and reverse input files, how can I do that? I have tried importing the two files giving space in between them, but it doesn't work. Can you please tell me is it okay to use only one input file to choose the best kmer?

kmer • 230 views
Entering edit mode
3 months ago
liorglic ▴ 870

In k-mer analyses (including KmerGenie), paired-end information is not used, so you can just use all your files - order doesn't matter. You have two options with KmerGenie:
1.Concatenate your fastq files:

cat hiseq_reads_ecoli_R1.fastq hiseq_reads_ecoli_R2.fastq > hiseq_reads_ecoli_R12.fastq
kmergenie hiseq_reads_ecoli_R12.fastq -l 21 -k 121 -s 2 -o out_file.txt -t 16

2.Use a reads list:

echo -e "hiseq_reads_ecoli_R1.fastq\nhiseq_reads_ecoli_R2.fastq" > reads.list
kmergenie reads.list -l 21 -k 121 -s 2 -o out_file.txt -t 16

Option 2 has the advantages that you don't need to duplicate your data and you can include as many files as you need.


Login before adding your answer.

Traffic: 852 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6