I would like to know why Mod gives the correct result and PowerMod gives an incorrect result.
Mod[0^0,7] PowerMod[0,0,7]
I think that is indeed the real reason, i was trying to point out that 0 is one of those `fringe' cases where things might be different from what you expect...