Hello I have been trying to run Quake Error Corrector on my Illumina reads before assembly. I am following the manual at http://www.cbcb.umd.edu/software/quake/manual.html. I have a problem running the final step of error correction as Quake always gives an error

"Kmers are not of expected length 21"

As I understand it the Quake program has 3 sub-programs inside it called - 1)Count qmers 2)Coverage Cutoff 3)Correct reads I have tried searching for examples to run Quake but in all of the cases, nobody uses the second program which is Cov_model.py. The Quake.py program is directly run like

quake.py -f pairs.list -k 21

However, there are no commandline options for cut-off value or for providing a count-file.

My commands in brief were-

1) Using jellyfish to get kmer counts.

jellyfish count -m 21 -s 27G -t 48 -C -F 2 File_1.fastq File_2.fastq

This generated a file called mer_counts.jf

2)The second module in jellyfish computes a histogram of k-mer ocurrences.

jellyfish histo -f -t 48 -v -o reads.histo mer_counts.jf

This generated a file called reads.histo containing number of kmers and their multiplicity

3)Then I can plug it into Quake and use

cov_model.py --int reads.histo

This worked well and gave a Coverage Cutoff of 1.

4) I went ahead using reads.histo to the final Correction step

correct -f input_files.txt -k 21 -c 1 -m reads.histo -p 48 -z --log

This just gives me an error saying "Kmers are not of expected length 21"

Should I just go ahead with using Quake.py? It does not provide any options to provide a cutoff file or a count file. Please advise!