If you are interested in only solutions on subdomains, near the set $\{x = 0, y > 0\}$, then your mistake is that in the method of characteristics, the fact that $q = a$ **only holds along the characteristic**. (Your supposition that $q = a$ everywhere already contradicts the initial data.)

Everything you did after that is not part of the method of characteristics. 

Given 
$$ F = p^2 + 2q - x $$
The characteristic equations are
$$ \dot{p} = 1, \dot{q} = 0, \dot{u} = 2p^2 + 2 q, \dot{x} = 2p, \dot{y} = 2 $$
You obtain that 
$$ p(s) = s + p(0), q(s) = q(0), u(s) = u(0) + \frac23 s^3 + 2 s^2 p(0) + 2p(0)^2 s + 2 q(0)s + u(0)$$
and
$$ x(s) = s^2 + 2p(0) s, y(s) = 2 s + y(0) $$
The initial data are given in terms of $y(0)$:
$$ q(0) = - 2y(0), p(0) = \pm 2\sqrt{ y(0)}, u(0) = - y(0)^2 $$