cd-hit简要原理
cd-hit是一款用于蛋白质序列或者核酸序列聚类的工具,由 Dr. Weizhong Li开发。其主要目的是在一堆较相似的序列找出几条代表序列,这几条代表序列就能分别代表这堆序列中的一个小类。
cd-hit算法的核心在于首先找出较长的序列作为代表序列,之后将剩下的序列与这些代表序列比对。某条序列与某代表序列比对达到自行设置的相似性阈值后则直接归为此代表序列一类,默认情况下某序列第一次找到符合条件的代表序列后就不再进行后续与其它代表序列的比对。
从其设计原理来看此工具也具有一些缺点,因为其首先找到长的代表序列,剩下的序列只与这些代表序列比对,而剩下的序列之间不用互相比对,某一类序列之间也没有互相比对,但剩下的序列间相似性差距可能也很大,这样就可能导致结果不够可靠。但是也正是因为这样比对的方式,不用所有序列间进行两两比对,使得这个工具运行速度很快。因此实际应用时要根据需要选择使不使用此工具。
cd-hit的参数
cd-hit一般在命令行下操作,只需要一行代码即可,十分方便,其有很多参数可以帮助改变聚类方式及聚类的相似性阈值等。我仅列了几个个人认为比较重要的参数,当然有需要的时候可以自己敲出cd-hit -h查看帮助文档哈!
-i input filename in fasta format, required
# 输入文件,必需。为fasta格式。基本格式为一行以'>'开头的注释行和单独的一列序列行反复重复。
-o output filename, required
# 输出文件名,必需。
-c sequence identity threshold, default 0.9
# 相似性阈值,默认为0.9,c必须和字长一起调整