R语言中通过ggmsa包展示多序列比对结果

2022-07-01  本文已影响0人  队长的生物实验室

多序列比对表示不同序列中的氨基酸/核苷酸的位点同源性。 将序列比对用于进化分析时,处于相同位置的氨基酸/核苷酸位点则被认为在进化上是同源的,并且具有共同的祖先 。

其中,可视化多序列比对结果有多种方式,而ggmsa包简单方便,可视化效果令人眼前一新。官方网址及说明https://github.com/YuLab-SMU/ggmsa
需要准备多序列比对结果(fasta格式文件)。

作图code

devtools::install_github("YuLab-SMU/ggmsa") #安装ggmsa包
fai<-"d:/Documents/KUA1_T-KUA1.fas"#读入文件
pdf(file = "aligned_fasta3.pdf",width = 30,height = 2)
ggmsa(fai)+
  geom_seqlogo() + geom_msaBar()
dev.off()#出图,width与height可自定义改变
多序列比对效果
上一篇 下一篇

猜你喜欢

热点阅读