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可自定义改变
