Hello. I am trying to cluster a huge fasta file containing using CD-HIT-EST with a threshold of 80%. According to the user's guide (http://www.bioinformatics.org/cd-hit/cd-hit-user-guide.pdf), I should use a word size (- n) of 5. However, it is taking forever. Could I change this parameter to -n 10 to speed up the process without changes in the final result, i. e., get the same result as -n 5?
This is my command:
cd-hit-est -i input -o output -d 0 -T 16 -g 0 -M 75000 -aL 0.97 -aS 0.97 -c 0.8 -n 5 -b 1
1 answer
In a word, no. The word size of 5 is already an upper limit for clustering at 80% identity. You will have to get a faster computer with more threads and memory, work with a smaller database, or just be patient. It may help to know that the process is not linear as the largest sequences are clustered first, so it will speed up as it goes along.
Log in to answer this question.
-
Why cd-hit-est not work when sequence identity threshold<0.95？
written by JZX •
-
How to use cd-hit-para.pl in SGE?
written by FelipeMSD •
-
duplicate gene IDs with different length after cd-hit-est dereplicate and cluster
written by phxu6780 •
-
Trying to use CD-HIT with a fasta file with 5 mers but not getting any output.
written by schlogl
-
how to speed up cdhit clustering?
written by bitpir
-
How can run cd-hit-est with a clstr threshold less than 0.8?
written by m.koohi.m
-
cd-hit threshold less than 0.7 is it possible??
written by casley_queiroz •
-
How stringent/loose should clustering be (CD-HIT-EST/USearch) for multi-kmer assemblies?
written by satshil.r •
-
Segmentation fault for cd-hit
written by Crystal •
-
Cd-Hit-Est-2D Segmentation Fault
written by Mariam Rizkallah •
