J2:=PrimitiveGroup(100,33);
IsSimple(J2);
orb:=[ o : o in Orbits(Stabilizer(J2,1)) | #o eq 36 ];
Gamma:=OrbitalGraph(J2,1,Set(orb[1]));
IntersectionArray(Gamma) eq [ 36, 21, 1, 12 ];
CompositionFactors(AutomorphismGroup(Gamma));