An example for $n=7$ is given in J. W. S. Cassels, Factorization of polynomials in several variables, Proc. Fifteenth Scandinavian Congress (Oslo, 1968), vol. 118, Lecture Notes in Mathematics, Springer, Berlin, pp. 1-17. This reference is from https://oeis.org/A112090 and is quite technical (using topology of Riemann surfaces).
here I type the example from the paper in a sort of computer-readable format
l=(1+sqrt(-7))/2
m=(1-sqrt(-7))/2
# t a nonzero(?) parameter
f(x)-g(y)=(x^3+l*x^3*y-m*x*y^2-y^3-(3*l+2)*t*x+(3*m+2)*t*y+t)*
(x^4-l*x^3*y-x^2*y^2-m*x*y^3+y^4+2*(m-l)*t*x^2-
7*t*x*y+2*(l-m)*t*y^2+(3-l)*t*x-(3-m)*t*y-7*t^2)
Sage script (with input f and (-)g taken from the paper
z=QQ['z'].0
K.<s>=NumberField(z^2+7)
l=<l>=NumberField(1+sz^2-z+2)/2
m=m=l.conjugate(1-s)/2
R.<x,y,t>=K[]
# t a nonzero(?) parameter
f=x^7-7*l*t*x^5+(4-l)*t*x^4+(14*l-35)*t^2*x^3-(8*l+10)*t^2*x^2+(3-l+7*(3*l+2)*t)*t^2*x
g=-y^7+7*m*t*y^5+(4-m)*t*y^4-(14*m-35)*t^2*y^3-(8*m+10)*t^2*y^2-(3-m+7*(3*m+2)*t)*t^2*y-7*t^3
(f+g).factor()
outputs
(x^3+x^3 + (1/2*s+1/2l)*x^2*y+*x^2*y + (1/2*sl - 1/2)*x*y^2 -y^3+ y^3 + (-3/2*s3*l -7/ 2)*x*t + (-3/2*s3*l + 7/25)*y*t + t) *
(x^4 + (-1/2*s - 1/2l)*x^3*y - x^2*y^2 + (1/2*sl - 1/2)*x*y^3 + y^4 +
(-2*s4*l+2)*x^2*t -
7*x*y*t + (2*s4*l-2)*y^2*t + (-1/2*sl + 5/23)*x*t + (-1/2*s l- 5/2)*y*t - 7*t^2)
which is the factorisation that should be like the one above. You are welcome to spotBut it is not - there is a typo, if any above (and in the paper: in the 1st factor the monomial lx^3y should be lx^2y. After this change everything checks out.