本篇内容介绍了“kallisto怎么使用”的有关知识，在实际案例的操作过程中，不少人都会遇到这样的困境，接下来就让小编带领大家学习一下如何处理这些情况吧！希望大家仔细阅读，能够学有所成！kallisto 是2016年发布的一款无须比对的转录本定量工具，采用了名为pseudo-alignment的算法。传统的定量算法是根据reads的比对位置来确认其属于哪个转录本或者基因，而pseudo-alignment 算法不关系reads具体的比对位置，而是通过reads的kmer特征来判断其属于哪一条转录本，示意图如下首先将每个转录本序列划分为kmer, 利用所有转录本的kmer序列构建de Bgujin Graph, 简称T-DBG，在这个图中，每个节点是一个kmer, 每条路径代表一个转录本， 由于转录本序列的冗余，实际上每个kmer对应多条路径，也就是对应多个转录本; 然后将测序的reads也划分为kmer, 并将其映射到T-DBG中。最终定量时，将该reads的所有kmer对应的转录本取交集，就能够分析出reads可能属于哪些转录本序列。官网有编译好的可执行文件，下载解压即可。代码如下解压之后，在文件夹下可以看到名为kallisto的可执行文件。从算法也可以看到，软件的运行需要两步，第一步对转录本的序列划分kmer, 构建T-DBG, 也称之为建索引；第二步对reads 定量。kallisto支持读取gzip压缩的转录本序列，用法如下只需要提供转录本的fasta格式的序列即可。-k参数指定kmer的长度，-i参数指定输出的索引的名字，注意kallisto建立的索引为一个文件。kallisto 支持单端和双端数据的定量，双端数据用法如下-i参数指定转录本的索引文件，-o参数指定输出结果的目录，-t参数指定线程数，kallisto支持gzip压缩的序列文件。
单端数据用法如下对于单端数据而言，必须指定fragment长度的均值和方差，分别对应-l和-s参数。
在输出目录，会生成以下3个文件run_info.json 文件为JSON格式，保存了运行的命令和参数。前缀为abunda 香港云主机nce 的文件，保存了转录本的定量信息。其中h6为HDF5格式的文件，当转录本数量较多时，相比纯文本，这种格式的文件大小会小很多；tsv为纯文本的文件，内容如下对于HDF5的文件，可以采用如下命令转换为tsv格式的文件-o参数指定输出结果的目录，最终生成的文件名称为abundance.tsv。“kallisto怎么使用”的内容就介绍到这里了，感谢大家的阅读。如果想了解更多行业相关的知识可以关注开发云网站，小编将为大家输出更多高质量的实用文章！
这篇文章给大家介绍如何解决DataFlow部署K8s应用的问题，内容非常详细，感兴趣的小伙伴们可以参考借鉴，希望对大家能有所帮助。基于各种原因，团队的Kubernetes被加了限制，必须在特定的Node才可以部署。而之前没有指定，所以Spring Cloud …
免责声明：本站发布的图片视频文字，以转载和分享为主，文章观点不代表本站立场，本站不承担相关法律责任；如果涉及侵权请联系邮箱：360163164@qq.com举报，并提供相关证据，经查实将立刻删除涉嫌侵权内容。
微信扫一扫 
