Yes. Since $W$ and $\xi$ are independent, we may write 

$$X_t = \int_{\mathbb R} \mathbb E[g(\xi, z)| \mathcal F^W_t] \, d\mu_Z,$$

as can be seen by, say, taking the [regular conditional probability](https://en.wikipedia.org/wiki/Regular_conditional_probability) with respect to $Z$.

We recognize that for every $z$, $Y^z_t: = \mathbb E[g(\xi, z)| \mathcal F_t^W]$ is a closable martingale with respect to the Brownian filtration. According to the results [here](https://mathoverflow.net/questions/402811/when-is-every-levy-martingale-of-a-process-a-continuous-martingale) and [here](https://math.stackexchange.com/questions/4011101/stopping-time-w-r-t-brownian-filtration-is-predictable), $Y^z_t$ is in fact continuous for every $z$. By the boundedness of $g$, $Y^z_t$ is also uniformly bounded.

Now the rest of the proof is analysis - we claim that $X_t$, being the average of continuous, uniformly bounded functions is also continuous almost surely. 

To see this, let

$$\phi(z, \delta, \omega) := \sup_{s, t; |s - t| < \delta} |Y_t^z (\omega) - Y_s^z (\omega)|$$

be a uniform modulus of continuity for $Y^z$, and let $M > 0$ be a uniform bound for $|g|$. 
By continuity of $Y^z_t$, we have for $\mu_z \times \mathbb P$-a.e. $(z, \omega)$ that $\lim_{\delta\to 0} \phi(z, \delta, \omega) = 0$. 

In other words, writing $E_{\varepsilon, \delta, \omega} := \{z \, | \, \phi(z, \delta, \omega) \leq \varepsilon\}$, we have for every $\varepsilon > 0$ that for $\mathbb P$-a.e. $\omega$ we have $\mu_Z(E_{\varepsilon, \delta, \omega}) \to 1$ as $\delta \to 0^+$.

If for each $\varepsilon > 0$, we write $N_{\varepsilon}$ for the $\mathbb P$-null set of exceptions to the above statement, we may set $N := \cup_{n \in \mathbb Z_+} N_{1/n}$ to find that for all $\omega$ in the $\mathbb P$-full measure set $\Omega \setminus N$, we have $\mu_Z(E_{\varepsilon, \delta, \omega}) \to 1$ as $\delta \to 0^+$, for every $\varepsilon > 0$.

Now let $\varepsilon > 0$ be arbitrary, and fix $\omega \in \Omega \setminus N$. Pick $\delta$ such that $\mu_z (E_{\varepsilon/2, \delta, \omega}) > 1 - \frac{\varepsilon}{2M}$. 

We then compute, for all $s, t$ with $|s - t| < \delta$,

$$|X_t (\omega) - X_s (\omega)|$$

$$= |\int_{\mathbb R} Y^z_t (\omega) - Y^z_s (\omega)\, d\mu_Z|$$

$$\leq \int_{\mathbb R}|Y^z_t (\omega) - Y^z_s (\omega)| \, d\mu_Z$$

$$ = \int_{E_{\varepsilon/2, \delta, \omega}} |Y^z_t (\omega) - Y^z_s (\omega)| \, d\mu_Z + \int_{E^c_{\varepsilon/2, \delta, \omega}} \mathbb |Y^z_t (\omega) - Y^z_s (\omega)| \, d\mu_Z$$

$$ < \frac{\varepsilon}{2} + M \frac{\varepsilon}{2M} = \varepsilon$$

and we conclude since $\varepsilon$ was arbitrary.