In principle we have two problems, the first a hypothesis-testing problem, and the second a confidence interval problem.
For the hypothesis testing, the null hypothesis is that the probability of "success" (a $2$ or a $4$) is $\frac{1}{3}$. Under that hypothesis, the number $X$ nof successes in $9000$ independent trials has binomial distribution, $n=9000$, $p=\frac{1}{3}$.
This binomial has mean $np=3000$ and variance $np(1-p)=2000$. Since $n$ is very large, and $p$ is not close to either $0$ or $1$, probabilities for $X$ are well-approximated by probabilities for the normal distribution, mean $3000$, variance $2000$. So the standard deviation is about $44.72$. Our actual result of $3240$ is about $\frac{3240-3000}{44.72}$ stahdard deviation units up from the mean under the null hypothesis. Calculation shows that the result is about $5.36$ standard deviation units up from the mean. That is an enormous number of standard deviation units: A result of $3240$ or more successes for all practical purposes cannot happen with fair dice.
For the confidence interval for $p$, the probability of success, again we use the normal approximation.
The sample mean $\frac{X}{9000}$ can be assumed to be normally distributed, with standard deviation
$$s=\frac{a(1-a)}{\sqrt{9000}},\tag{1}$$
where $a=\frac{3240}{9000}$.
The number obtained using (1) is not necessarily the true standard deviation of $\frac{X}{9000}$, but it is good enough, particularly since $\sqrt{x(1-x)}$ is quite insensitive to smallish changes in $x$ when $x$ is not very far from $\frac{1}{2}$.
Unless there is good reason not to, one might as well go for a $95\%$ confidence interval. This will be the interval
$$\frac{3240}{9000}-1.96 s \le p \le \frac{3240}{9000}+1.96 s,$$
where $s$ is the number obtained in (1).