Visualize the similarity between two motifs, Very ineficent,
1 | motifs2View(m1, m2, reg, Sequences, nearHeights = FALSE)
|
1 2 3 4 5 6 7 | FastaFile<-"~/Dropbox/UTX-Alex/jan/combined.fasta"
Sequences <- readDNAStringSet(FastaFile, "fasta")
data<-loadHeightFile("~/Dropbox/UTX-Alex/jan/combined_heights.bed")$data
reg<-as.matrix(ascore(data,c(1,1,1),c("top","bottom","middle"),c(1,1,1)))
m<-motifs2View("CANNTG","TGACCT",reg[,3],Sequences)
k<-motifs2View("CANNTG","GATAAG",reg[,1],Sequences)
t0<-unlist(lapply(seq(min(k),max(k)),function(i)length(which(k==i))))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.