I have a peptide variant sequences, I want to use seqlogo to visualize the amino acids compositions of each position. But I also want to exclude the WT sequence, I do not know how to achieve that with seqlogo in R. I will be very appreciate if anyone can help me. The WT sequence is GNNRPVYIPQPRPPHPRL. Thank you! Here is my R script in MacOS.
Hi,
I have a peptide variant sequences, I want to use seqlogo to visualize the amino acids compositions of each position. But I also want to exclude the WT sequence, I do not know how to achieve that with seqlogo in R. I will be very appreciate if anyone can help me. The WT sequence is GNNRPVYIPQPRPPHPRL. Thank you! Here is my R script in MacOS.
library(ggplot2) library(ggseqlogo) library(stringr) library(ggsci) library(tidyverse) header=T ggseqlogo(as.character(seqlog$seq), method = 'prob') + theme_bw() + scale_y_continuous(labels = scales::percent)