M := ModularForms(Gamma0(389),2);
M;
S := CuspidalSubspace(M);
S;
time N := Newforms(S);
fE := N[1][1];
// The newform attached to E.
qExpansion(fE,10);
fA := N[5][1];
// One of the 20 Galois-conjugate q-expansions of a newform attached to A.
qExpansion(fA,4);

// The q-expansions of the conjugates of f_A mod 5
time Reductions(fA,5);
// The last one is equal to the q-expansions of f_E mod 5
time Reductions(fE,5);