This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
archive:bioinformatic_tools:quake [2011/05/08 04:53] eyliaw [Running Quake] |
archive:bioinformatic_tools:quake [2011/05/09 18:02] eyliaw [Running Quake] |
||
---|---|---|---|
Line 15: | Line 15: | ||
Finally, correct the reads: | Finally, correct the reads: | ||
- | correct -f [fastq list file] -k [k-mer size] -c [cutoff] -m [counts file] -p [number of cores] | + | correct -f [fastq file list] -k [k-mer size] -c [cutoff] -m [counts file] -p [number of cores] -z (gzips the output) |
+ | |||
+ | In the file list, you should tab-separate paired end reads. Also, be sure that all .'s in the sequence are written as N's. | ||
[Kevin] I think that we want to select the k-mer size manually, rather than relying on quake. Their default cutoff is very conservative, and we'll probably do better over-correcting than under-correcting. | [Kevin] I think that we want to select the k-mer size manually, rather than relying on quake. Their default cutoff is very conservative, and we'll probably do better over-correcting than under-correcting. |