0. 简介
Weblogo图:即Sequence logo,生物信息学中对核酸或蛋白质序列保守性的可视化图像。
1.绘图数据
文件由长度相等的目标序列构成(不足长度的用"-"代替),前20行数据如下:
2. 绘图代码
library(ggplot2)
library(ggseqlogo)
data_source = read.table('C:\\Users\\19535\\Desktop\\sequence.csv')
data_source = as.vector(data_source$V1)
ggseqlogo(data_source,seqs_aa$AKT1,method='bits') +
theme_logo() +
theme(legend.position = 'none') +
theme(plot.title=element_text(hjust=0.5, size=15, face="bold")) +
ggtitle('sequence')