model{for (j in 1:min_FASTA_length){for (i in 1:2){matrix[i,1:length(all_of_letters), j] ~ dmulti(theta[i,1:length(all_of_letters), j], nS[j, i])theta[i,1:length(all_of_letters),j] ~ ddirch(vcounts[1:length(all_of_letters)] + 0.03)}}for (i in 1:length(all_of_letters)){ vcounts[i] ~ dgamma(0.1,0.1)}}
