This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Next revision Both sides next revision | ||
archive:bioinformatic_tools:quake [2011/05/09 17:53] eyliaw [Running Quake] |
archive:bioinformatic_tools:quake [2011/05/09 18:02] eyliaw [Running Quake] |
||
---|---|---|---|
Line 17: | Line 17: | ||
correct -f [fastq file list] -k [k-mer size] -c [cutoff] -m [counts file] -p [number of cores] -z (gzips the output) | 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. | + | 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. |