Fix a positive integer $k$. Let $b_n(q)$ be the polynomial defined by the recursive equation $$b_n(q)=\binom{n+k-1}k+(1+q^{n-1})b_{n-1}(q), \qquad n\geq1,$$ initializing with $b_0(q)=0$.

I run into the below through experiment.

QUESTION. Denote $(q)_m=(1-q)(1-q^2)\cdots(1-q^m)$. Is this true? $$\lim_{n\rightarrow\infty}\left(b_n(q)-\binom{n+k}{k+1}\right) =\sum_{m=1}^{\infty}\frac{q^{\binom{m+1}2}}{(q)_m(1-q^m)^{k+1}}.$$