$\DeclareMathOperator\prob{prob}$Alapan Das' clever argument may be rephrased on the probabilistic language. Write $[m]=\{1,2,\dotsc,m\}$. Choose a random non-empty subset $A\subset [n]$ (all $2^n-1$ possible outcomes having equal probabilities). Then choose a random element $\xi\in A$ uniformly. Denote $p=\prob (\xi=\max(A))$. On one hand, denoting $j=\lvert A\rvert$ we get $$ p=\sum_{j=1}^n \prob(\xi=\max(A)\mathrel||A|=j)\cdot \prob(|A|=j)=\sum_{j=1}^n\frac1j \cdot\frac{{n\choose j}}{2^n-1} $$
On the other hand, $$ p=\sum_{k=1}^n \prob(\xi=k \, \&\,A\subset [k])= \sum_{k=1}^n \prob(\xi=k|A\subset [k])\cdot \prob(A\subset [k])\\= \sum_{k=1}^n \frac1k\cdot \frac{2^k-1}{2^n-1}. $$$$ p=\sum_{k=1}^n \prob(\xi=k \, \&\,A\subset [k])= \sum_{k=1}^n \prob(\xi=k\mathrel|A\subset [k])\cdot \prob(A\subset [k])\\= \sum_{k=1}^n \frac1k\cdot \frac{2^k-1}{2^n-1}. $$