Let
$$ I = \int_0^1W_tdt, $$
where $W_t$ is a Brownian motion. From Integral of Brownian motion w.r.t. time we have that
$$ \mathbb{E}[I]=0, $$ by Fubini's theorem. And that $$ \mathbb{V}\text{ar}[I] = \mathbb{E}[I^2] = \mathbb{E}\left[\int_0^1\int_0^1W_sW_tdsdt\right] = \int_0^1\int_0^1\mathbb{E}[W_sW_t]dsdt = \frac{1}{3}. $$ Here I do not understand how we are able to bring the expectation into the integral, Ito isometry says that $$ \mathbb{E}\left[\left(\int_0^tX_sdW_s\right)^2\right] = \int_0^t\mathbb{E}[X_s^2]ds $$ so I am unsure how it is possible to use this here?
Furthermore, I am required to find $$ \mathbb{C}\text{ov}[I, W_1], $$
I know that $$ \mathbb{C}\text{ov}[I, W_1] = \mathbb{E}[IW_1], $$ since $\mathbb{E}[I] = \mathbb{E}[W_1] = 0$, my best guess here is then that $$ \mathbb{E}[IW_1] = \min\left(\frac{1}{3},1\right) = \frac{1}{3}, $$ however, this is probably wrong.
Thank you in advance!