Are the 'semi' trivial zeros of $\zeta(s) \pm \zeta(1-s)$ all on the critical line? The proof that $\Gamma(z)\pm \Gamma(1-z)$ only has zeros for $z \in \mathbb{R}$ or $z= \frac12 +i \mathbb{R}$ has been given here:
Are all zeros of $\Gamma(s) \pm \Gamma(1-s)$ on a line with real part = $\frac12$ ?
An obvious follow up question is whether $\zeta(s) \pm \zeta(1-s)$ also has zeros (other than its non-trivial ones that would induce 0+0 or 0-0). 
This is indeed the case and $\zeta(s)^2 - \zeta(1-s)^2$ has the following zeros:
$\frac12 \pm 0.819545329 i$
$\frac12 \pm 3.436218226 i$
$\frac12 \pm 9.666908056 i$
$\frac12 \pm 14.13472514 i$ (the first non trivial)
$\frac12 \pm 14.51791963 i$
$\frac12 \pm 17.84559954 i$
$\dots$
These 'semi' trivial zeros appear to all lie on the critical line. I wonder if anything is known or proven about their location (I guess not, since a proof that they must have real part of $\frac12$ would automatically imply RH, right?).
EDIT:
Two counterexamples found by Joro in the answers below. Both have real parts outside the critical strip, so I would like to rephrase my question as:
Are the 'semi' trivial zeros that are located within the critical strip all on the critical line?
 A: Here are two images to complement joro's image.
Zeroes of $ \zeta(s)^2 - \zeta(1-s)^2 $:

Zeroes of $\left ( 2^s \pi^{s−1} \sin \left ( \frac{\pi s}{2} \right ) \Gamma (1-s) \right ) ^2 - 1$:

In particular it seems to be the case that all zeroes inside the critical strip are on the critical line, and that all other zeroes don't stray too far from the origin, as hinted at by GH's answer.
A: Slightly off-topic : this article from Arxiv (in french) shows that zeros of functions f : s -> h(s) - h(1-s), where h is a meromorphic function satisfying appropriate growth conditions, are inclined to lie on the critical line, so this phenomenon is not specific to gamma functions.
I didn't check that the main result in the article apply to the function in GH's answer above.
A: If $\zeta(s)$ is nonzero, but $\zeta(s)\pm\zeta(1-s)=0$, then by the functional equation of the Riemann zeta function we have
$$ \pi^{-\frac{s}{2}}\Gamma\left(\frac{s}{2}\right)\pm \pi^{-\frac{1-s}{2}}\Gamma\left(\frac{1-s}{2}\right)=0.$$
That is, your question is just the Riemann Hypothesis plus a more elementary one similar to your earlier question about the zeros of $\Gamma(s)\pm\Gamma(1-s)$. I would expect that the exact same techniques work here, i.e. one can show by known estimates for the gamma function that all nonreal solutions of the displayed equation lie on $\Re s=1/2$.
EDIT 1. To keep up with new developments I now expect that within the critical strip all nonreal solutions of the displayed equation lie on $\Re s=1/2$. Moreover, it seems reasonable to believe that there are no nonreal solutions with $|\Re s|$ sufficiently large.
EDIT 2. It follows from a generalized Rouché's theorem and Stirling's approximation that there are no nonreal solutions with $|\Re s|$ sufficiently large. More precisely, consider the rectangular contour $C_n$ with vertices $2n\pm it$ and $2n+2\pm it$, where $n>0$ is a large integer and $t>0$ is sufficiently large in terms of $n$. It suffices to show that along $C_n$ we have 
$$ \left|\pi^{-\frac{1-s}{2}}\Gamma\left(\frac{1-s}{2}\right)\right|<\left|\pi^{-\frac{s}{2}}\Gamma\left(\frac{s}{2}\right)\right|,$$
because this implies that inside $C_n$ there is precisely one solution of the above displayed equation (which must be real by the reflection principle). One can show that the right hand side divided by the left hand side is
$ \gg n^{2n-\frac{1}{2}}(\pi e)^{-2n}$ on the vertical sides of $C_n$, while it is $\gg_n t^{2n-\frac{1}{2}}$ on the horizontal sides of $C_n$. The claim follows.
A: We now have strong results from GH for $\Gamma(s) \pm \Gamma(1-s)$ and $\zeta(s) \pm \zeta(1-s)$ (when both are nonzero), as well as a more generalized claim outlined in this article http://arxiv.org/abs/0712.1266, that zeros for $f(s) = h(s) - h(2a-s)$, with $h$ being a meromorphic function satisfying appropriate growth conditions, are all inclined to lie on a critical line $a$, with only a finite number of exceptions.
However, despite the encouraging outcomes that all zeros are on the critical line and only a finite number of exceptions residing outside the critical strip, there doesn't seem to be any way to extrapolate these results to the non trivial zeros ($\rho$). No new constraints seem to be imposed on their location and they really seem to originate "from deep down the unaccessible cave" of $\zeta(s)$ itself. They therefore could still reside anywhere on the critical strip.
Just to share some 'futile attempts' I made to make a link. Firstly I had hoped that somehow the results would force the following outcome for all $\Re s$ in the strip:
$\lim_{s\to\rho}\left|\frac{\zeta(s)}{\zeta(1-s)}\right|=1$
and thereby seriously limit the possible values $\Re(\rho)$ could assume. But they don't.
In an attempt to find another relation between $\zeta(s)2$ and $\zeta(1-s)$, I converted both of them to the alternating $\eta(s)$ function and then paired up the individual terms for each $n$ (this is allowed since $\eta(s)$ is valid for $\Re s>0$). This gives:
$\displaystyle  \sum _{n=1}^{\infty } \left({\frac { \left( -1 \right) ^{n-1}}{(1-{2}^{1-s}) {n}^{s}}} \pm{\frac { \left( -1 \right) ^{n-1}}{ (1-{2}^{s}) {n}^{1-s}}}\right)$
For each individual $n$, this yields a wave with a fixed frequency and amplitude, that only has zeros when $\Re s=1/2$. These waves nicely sum up to a curve that produces all the 'non & semi' trivial zeros from the OP. However, the $\rho$s obviously do arise from summing up the individual terms as well and I could not find a meaningful way to smartly swap out left and right terms, so that maybe new information about the non trivial zeros would be revealed.
As a last attempt, I also experimented with the PrimeZeta function $P(s)$. 
For $P(s)-P(1-s)$ I found:
$\displaystyle \sum _{k=1}^{\infty }  \frac{\mu \left( k \right)}{k} \ln  \left( {\frac {\zeta \left( ks \right) }{\zeta  \left( k \left( 1-s \right)  \right) }} \right)$
and for $P(s)+P(1-s)$:
$\displaystyle \sum _{k=1}^{\infty }  \frac{\mu \left( k \right)}{k} \ln  \left( {\zeta \left( ks \right) \zeta  \left( k \left( 1-s \right)  \right) } \right)$
For both functions, all zeros appear to lie on the critical line $\Re s=\frac12$, but for the first function the non-trivial zeros have now turned into poles. Interesting to note that the Wolfram site on the PrimeZeta states: "According to Fröberg (1968), very little is known about the roots", so maybe there is something new here to proof for $P(s) \pm P(1-s)$... :-) 
