I would like to know how to reverse the sizing in ggseqlogo() so that those loci with more than a nucleotide conflict are larger than those with identical nucleotides. Im working with some RNA editing sites and I'd like the A-I edited sites to stand out.
I have created an input for ggseqlogo() that uses four identical copies of unedited sequence (character string 'unalt1') and four strings identicle to those except for a G in place of an A. When plugging these into the program as follows...
mix1 <-c(unalt1,unalt1,unalt1,alt1 ,alt1 ,alt1 ,alt1 )
pmix1 = ggseqlogo(mix1, method = 'bits' )
... I get an image like the following:
Where the A/G site is much smaller than the surrounding nucleotides. Is there a simple way to reverse this? I know in the manual a description of how to generate a matrix with custom sizing is given, but I cant think how to incorporate a size matrix with these long nucleotide strings. Id appreciate any help!
