why (4*i+7) in Z[i] not a prime?norm(4*i+7) =65, 65mod4=1

edit

K.<a> = NumberField(x^2 +1);
(4*a+7).norm()
65
mod (65,4)
1

factor(4*a+7)
(-3*a - 2) * (a - 2)

(4*a+7).multiplicative_order();
+Infinity

K.elements_of_norm(65);factor(-8*a + 1);K.prime_factors(4*a -7);

-8*a - 1, 4*a - 7, 4*a + 7, -8*a + 1]
(a - 2) * (3*a - 2)
[Fractional ideal (3*a - 2), Fractional ideal (-a - 2)]
[Fractional ideal (a - 2), Fractional ideal (3*a - 2)]