This (in particular the 2nd equality) is incorrect. Let $\omega \in \Omega$ denote the sample realization. You state:
$$\int_{0}^{T} B_{t} (\omega) d t=T^{1+1/2} \int_{0}^{1} B_{s}(\omega) d s$$
For example assume $T>1$. This is claiming that in order to compute the integral you only need to consider realizations of $B_t$ between $[0,1]$ and then scale them up. However, the value of this integral also depends on the random realizations of $B_t$ between $[1,T]$, which cannot be deduced simply from its values between $[0,1]$. The LHS can take a value clearly different from the RHS depending on $\omega$. Intuitively if $W_t$ is e.g. a stock price, the integral represents a type of weighted average of the stock price over $[0,T]$. But you cannot calculate this average merely based on the stock price over $[0,1]$.
The probability that $\int_{1}^{T} B_{s}(\omega) d s$ takes any particular fixed value is zero (this integral is normally distributed). Hence the LHS and RHS are almost surely not equal. However, as @Kevin pointed out the LHS and RHS are equal in distribution that is these two random variables have the same moments.