这篇文章主要为大家展示了“如何使用CD-HIT合并pacbio转录本”，内容简而易懂，条理清晰，希望能够帮助大家解决疑惑，下面让小编带领大家一起研究并学习一下“如何使用CD-HIT合并pacbio转录本”这篇文章吧。
CD-HIT合并pacbio转录本
Pacbio三代全长转录组采用isoseq3 进行了转录本的筛选，聚类和校正之后得到高质量的转录本序列，可以对序列再进行一定的合并操作。
采用cDNA_Cupcaake 文档中提供的cd-hit方案，运行命令如下：
cd-hit-est -i <input> -o <output> -c 0.99 -T 6 -G 0 -aL 0.90 -AL 100 -aS 0.99 -AS 30
参数说明：
input: 是pacbio 的isoseq 3 输出的高质量转录本序列文件
output: 输出的合并序列文件
-c : 序列相似性，99%的相似性
-T ： 线程数，6个线程
-G： 采用的比对策略，默认是1，采用全局比对。这地方设置为0，将采用局部比对策略。
-aL: 长序列的覆盖度，设置为90%
-AL: 长序列覆盖区域的长度，设置为最小100bp
-aS: 短序列的覆盖度，设置为99%
-AS: 短序列覆盖区域的最小长度，设置为30bp
以上是“如何使用CD-HIT合并pacbio转录本”这篇文章的所有内容，感谢各位的阅读！相信大家都有了一定的了解，希望分享的内容对大家有所帮助，如果还想学习更多知识，欢迎关注亿速云行业资讯频道！
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
