使用cd-hit对核酸序列或氨基酸序列聚类


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必须和字长一起调整
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值