R语言绘图-4-Weblogo图

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')

结果：