calcscore<-function(seq,log2pwm){
# seq is a vector representing input DNA numerically
# log2pwm is a PWM (4xL) with elements as log base 2
score <- 0
for (j in 1:length(log2pwm[1,])){
score<-score+log2pwm[seq[j],j]}
return(score)
}