Here's another pass at this:
Suppose for some $k$, we have a positive integer $m$ such that $2^k-1 = mk$. Or, $2^k \equiv 1 (m)$
Therefore, in the vector space $\mathbb{F}_2^k$, there are $mk$ non-zero vectors. Partition these into $m$ disjoint sets of $k$ vectors: $S_i$ for $i=1,2, \ldots m$. Choose vectors $v_i \in S_i$ to form $\{v_1, v_2, \ldots v_k\}$. Now for any other vector, $w$, not in this set, the collection $\{v_1, v_2, \ldots v_k\}\cup\{w\}$ is linearly dependent. Therefore, the set spans $\mathbb{F}_2^k$ and thus it is a basis. Therefore, the number of possible bases is $m^k$. But then, $o(GL_k(\mathbb{F}_2)) = (2^k-1)(2^k-2)\ldots(2^k-2^{k-1})$
Older attempt
Hope this is right. Suppose there is a $k$ such that $2^k - 1 \equiv 0 (k)$.
That is, $1+2+2^2+ \ldots +2^{k-1} \equiv 0 (k)$, or, adding the two up:
$2+2^2+\ldots +2^k \equiv 0 (k)$
But this can't be since $(2,k)=1$ as, by hypothesis, $2$ is invertible in $\mathbb{Z}/k\mathbb{Z}$