For each prime v, the characteristic polynomial of the matrix
rho(Frob_v) in GL_2(Z_2) or GL_2(Z/8Z) can be computed. This gives
plenty of conjugacy classes in GL_2(Z/8) that intersect your Galois
group. (At least when there are no roots or distinct roots, we can
identify the conjugacy class in this way.) On the other hand you know
an upper bound on the image from the isogeny over Q. Hopefully you can
prove that it is all of this. Sometimes one can also use Tate-curves
at multiplicative primes v to get the image of the decomposition group
D_v to generate more. But here c_5 is divisible by 2 so it is not
straight forward. Otherwise one gets elements like [1 1\\ 0 1] mod p
for free.
But as I said, I have not done the computations here at all and maybe
this is all irrelevant.
> (That's because Q(E[2])=Q
are we talking about the same curve ? For 20a1 it is Q(i), isn't it ?
Chris.