So we have that
$$
t\left( {t + 1} \right) \cdots \left( {t + n} \right) = t^{\,\overline {\,n + 1\,} }
$$
and
$$
t^{\,\overline {\,n + 1\,} } = \sum\limits_{\,\left( {1\, \leqslant } \right)\,k\,\left( { \leqslant \,n + 1} \right)} {\left[ \begin{gathered}
n + 1 \\
k \\
\end{gathered} \right]t^{\,k} }
$$
Now
$$
P(k;\,q) = {1 \over {q!}}\left[ \matrix{
q \cr
k \cr} \right]
$$
is a discrete Probability distribution
(see for instance this post)
whose raw first and second moment are
$$
\eqalign{
& E\left[ k \right] = \sum\limits_{\left( {0\, \le } \right)\,g\,\left( { \le \,q} \right)} {{k \over {q!}}\left[ \matrix{
q \cr
k \cr} \right]} = H(q) = \gamma + \psi (q + 1) \cr
& E\left[ {k^{\,2} } \right] = \sum\limits_{\left( {0\, \le } \right)\,k} {{{k^{\,2} } \over {q!}}\left[ \matrix{
q \cr
k \cr} \right] = } \left( {H_{\,q} ^2 + H_{\,q} - {{\pi ^2 } \over 6} + \psi ^{\left( 1 \right)} \left( {q + 1} \right)} \right) \cr}
$$
Your question concerns the mode of $\left[ \matrix{ q \cr k \cr} \right]$ which coincides with that of $P(k;\,q)$,
but I could not find sufficient information on line about this distribution and its characteristics.
Plotting some examples of $P(k;\,q)$ it comes out that it has quite a peaked and symmetric shape concentrated in the lower values of $k$
and that the mode remains very close to the mean.
So ${\rm mode}\,P(k;q) \approx H(q)$.
If you need a more analytical result, then you can try and interpolate $P(k;\,q)$ by a continuous function,
which has a closed expression for the mode.