Ontology is not yet a thing in this SE—perhaps the philosophy SE could be more forgiving to it. In any case, it is assumed that you are in complete and unequivocal peace with the uncertainty principle (UP), the way it has been settled in the last 95 years or so, and you merely wish to understand how it presents itself in deformation quantization, without expecting new handles or insights into it, somehow unavailable to the operators-in-Hilbert-space traditional formulation.
The knowables in QM are mathematical (=logical) facts, and the expectation values of measurements, and, significantly, nothing else. There is no answerable "why" there. QM is a gimmick producing reliable expectation values.
Can we know momentum and position simultaneously according to uncertainty principle? if not, how we talk about c-number functions on phase space? and why do we need to restrict ourselves to expectation values from theoretical point of view ?
In that sense, the answer to your first question here is: of course not, whenever "know" dares go beyond expectations of measurable quantities. The third was already summarily dispatched above: it is a cornerstone of the theory; those not at peace with it have been walking in hapless circles for a century, mumbling. So, I will try to explain the second one.
Deformation quantization works in phase space, where you may "think" about points (x,p) in it, without the naive, fallacious association with feasible measurements. One may think about them the way thinks about wavefunctions, that is formal ingredients leading to the answers.
The UP is baked into the theory from its foundations, perhaps counterintuitively, which lots of physicists (including Dirac) struggled to understand in the late forties, and Moyal triumphantly succeeded in. (His free online biography, by his wife details how). For starters,
The Wigner function f(x,p) characterizing the state cannot be a δ function, since $|f(x,p)|\leq 2/h$, a trivial consequence of its definition. In fact, since it is normalized to 1 in phase space, you may see that its boxy-skyscraper height L=1/A where A is its footprint area; so the particle is condemned to spread over a phase-space area larger than h/2. We can talk about arbitrary xs and ps, as long as we don't picture a particle as something localized in an area smaller than h/2. Moreover,
The Wigner function f is not an arbitrary function in phase space. For a pure state, it satisfies an almost "mystical" property (in the sense it used to mystify novices in the area). For an arbitrary function g,
$$
\langle g^* \star g\rangle = \int\!\! dx dp ~ f (g^* \star g) \geq 0.
$$
This is the intrinsic constraint serving to prove the UP in this formulation, tactically equivalent to the positivity of the norm in Hilbert space. Finally,
It is inappropriate to think of the Wigner function as a probability versus mutually exclusive contingencies. Two nearby points in phase space necessarily participate in describing the state (particle) and cannot represent independent locations. Dramatically, the Wigner function can go negative in such small areas, $O(h)$ (only). The uncertainty principle prevents you from visualizing a "particle galvanometer" in such small areas, to "locate" ("know the phase-space location") of a particle. The UP is the invisible cloak that prevents direct observation of negative f regions and localizing particles.
The upshot is that, contrary to the intuition you developed from classical statistical mechanics, the peculiar rules of the nonlocal star product endow phase-space variables, observables, and distributions with completely different properties than one's naive imagination would suggest.
You saw how different. Any phase-space function, such as your oscillator hamiltonian, is poised to be star-multiplied, so, then, perform the function of an operator. In that sense, it means something quite different than in classical mechanics. (Even though, even there, it is but an input to the Poisson-Bracket evolution operator...)