--- a/mutation_analysis.r	Mon Aug 18 04:42:59 2014 -0400
+++ b/mutation_analysis.r	Wed Sep 10 10:33:29 2014 -0400
@@ -294,35 +294,15 @@
 cat(length(dat$Sequence.ID), file="n.txt")
 
 
-
-
-
-
-
-
-
+weblogo = dat[,c("Sequence.ID", "VGene")]
+weblogo$VGene = gsub("\\*.*", "", weblogo$VGene)
 
-
-
-
-
-
-
-
+rs12 = read.table("HS12RSS.txt", sep="\t", header=TRUE)
+rs23 = read.table("HS23RSS.txt", sep="\t", header=TRUE)
 
-
-
-
-
-
-
-
+result12 = merge(weblogo, rs12, by.x="VGene", by.y="Gene")
+result23 = merge(weblogo, rs23, by.x="VGene", by.y="Gene")
 
-
-
-
+write.table(x=result12$Sequence, file="weblogo_in_rs12.txt", sep=",",quote=F,row.names=F,col.names=F)
+write.table(x=result23$Sequence, file="weblogo_in_rs23.txt", sep=",",quote=F,row.names=F,col.names=F)
 
-
-
-
-