  In this note, we prove a general identity between a $q$-multisum $B_N(q)$ and
a sum of $N^2$ products of quotients of theta functions. The $q$-multisum
$B_N(q)$ recently arose in the computation of a probability involving modules
over finite chain rings.
