linux cd-hit下载安装,CD-HIT学习

CD-hit

参数解读

-i 设置输入文件

-o 设置输出文件，可以将每次分析的ID阈值放到名称中，方便以后使用，如clean90，就是被清洗后，使用-c 0.90的分析结果

-c 设置ID阈值

-n 在ID各个范围内，作者给了一些设置mer值得建议

Choose of word size:

-n 5 for thresholds 0.7 ~ 1.0

-n 4 for thresholds 0.6 ~ 0.7

-n 3 for thresholds 0.5 ~ 0.6

-n 2 for thresholds 0.4 ~ 0.5

-d clstr文件中的描述字符长度，默认为20，设置为0的时候只去除第一个空格前面的字符

-M 设置使用内存大小，单位为MB

-T 设置核心数

-G 设置全局比对还是局部比对，如果这个参数设置为0，那么最好是配合-aS使用，避免过段的高质量匹配产生的无意义信息扰乱下游分析

-aS 比对序列的长度至少需要占较短序列(即冗余序列，因为较长的序列是代表序列)的比例

-AS 设置alignment coverage of short sequence，即短序列中，未在alignment中的序列必须短于这个值

-g 精确模式，但是会更慢，但是无论是否开启，代表representative序列不会变

-sc 对输出的cluster进行排序，大的cluster先输出

层级聚类

先对各个层级聚类生成的序列文件再次聚类，生成很多clstr文件

使用clstr_rev.pl聚合多个聚类文件

如果需要设置很低的ID值，那么使用psi-cd-hit可以达到这个目标，常用的代码：

cd-hit -i db -o db_90 -c 0.9 -n 5 -g 1 -G 0 -aS 0.8 -d 0 -p 1 -T 16 -M 0 > db_90.log

cd-hit -i db_90 -o db_60 -c 0.6 -n 4 -g 1 -G 0 -aS 0.8 -d 0 -p 1 -T 16 -M 0 > db_60.log

./psi-cd-hit.pl -i db_60 -o db_30 -c 0.3 -ce 1e-6 -aS 0.8 -G 0 -g 1 -exec local -core 16

clstr_rev.pl db_90.clstr db_60.clstr > db90-60.clstr

clstr_rev.pl db90-60.clstr db_30.clstr > db90-60-30.clstr

Here,

-ce 1e-6 and -c 0.3 means cutoff at either 30% identity or 1e-6 e.value

-G 0 means global identity

-aS 0.8 means that alignment must cover 80% of shorter (redundant) sequence

-g 1, slow but accurate mode, allowing sequences to be grouped to its most similar cluster

-core 16, use 16 threads for blast search

-clstr_rev.pl will combined all the cd-hit runs, see Hierarchically clustering