WebLogo

WebLogo: https://weblogo.berkeley.edu/logo.cgi

输入：

输出： 

序列分析图（sequence logo）

R包：https://cloud.tencent.com/developer/article/1172910

install.packages("ggseqlogo")

library(ggplot2)

library(ggseqlogo)

setwd("D:/biodata/01-浙大数据/")

getwd()

read <- read.table("a.txt")

read <- data.frame(t(read))

colnames(read) <- c("1","2","3","4","5","6","7","8")

rownames(read) <- c("A","C","G","T")

matrix_input <- read.table(all, header=T, row.names=1)

matrix_input <- as.matrix(matrix_input)

ggseqlogo(as.matrix(read),method="bits",seq_type="DNA")

ggseqlogo(as.matrix(read),method="prob")