pyani是一个很方便可以计算平均核酸一致性的软件,用于从基因组水平比较多个基因组序列之间的相似性。
推荐使用conda安装:conda install pyani
安装完成后,将所需要进行分析的基因组文件放到同一个文件夹中即可。
命令行例子:
average_nucleotide_identity.py -i 基因组所在文件夹地址 -o 想要输出的文件夹地址(自动新建) -m 所选算法 -g
具体例子:
![](https://img.haomeiwen.com/i19462250/8e4f0cfd5b878cc4.jpg)
算法有4种:ANIm (使用MUMmer去比对输入序列)、ANIb (使用BLASTN+去比对输入序列的1020nt片段)、ANIblastall (使用旧版BLASTN比对输入序列的1020nt片段)、TETRA ( 计算每个输入序列的四核苷酸频率)
-g 为输出图像文件
本次使用了6株从NCBI上下载的大肠杆菌菌株基因组进行分析
软件运行完后,可得到类似的结果:
![](https://img.haomeiwen.com/i19462250/c0679fdb59a78eb1.png)
除了图片,还有每个图所对应的比对数据结果,可以具体查询到相应的相似度/覆盖度等。
![](https://img.haomeiwen.com/i19462250/b71593070873a915.jpg)
网友评论