dev:=function(g,bb)
return Union(List(bb,x->Orbit(g,x,OnSets)));
end;

g:=Group((2,3,4)(5,6,7,8,9,10)(11,12,13,14,15,16),(1,5)(2,12)(3,15)(4,8)(6,14)(7,16)(9,10)(11,13));

bb2:=[[1, 2, 3, 4, 5, 6, 13], [1, 2, 3, 4, 5, 6, 14],  [1, 2, 3, 5, 6, 7, 11], [1, 2, 3, 5, 6, 8, 9],
[1, 2, 3, 5, 6, 9, 10],  [1, 2, 3, 5, 6, 9, 12], [1, 2, 3, 5, 6, 10, 15], [1, 2, 3, 5, 6, 14, 16],
[1, 2, 3, 5, 8, 11, 12], [1, 2, 5, 6, 7, 8, 16],  [1, 2, 5, 6, 7, 9, 14],  [1, 2, 5, 6, 7, 12, 13],
[1, 2, 5, 6, 7, 14, 15]];

bb5:=[[1, 2, 3, 4, 5, 6, 7], [1, 2, 3, 4, 5, 6, 14], [1, 2, 3, 4, 5, 6, 15], [1, 2, 3, 4, 5, 8, 12],
[1, 2, 3, 5, 6, 7, 10], [1, 2, 3, 5, 6, 7, 12], [1, 2, 3, 5, 6, 7, 15], [1, 2, 3, 5, 6, 8, 10],
[1, 2, 3, 5, 6, 8, 12], [1, 2, 3, 5, 6, 8, 13], [1, 2, 3, 5, 6, 9, 10], [1, 2, 3, 5, 6, 9, 12],
[1, 2, 3, 5, 6, 9, 13], [1, 2, 3, 5, 6, 10, 14], [1, 2, 3, 5, 6, 10, 16], [1, 2, 3, 5, 6, 12, 13],
[1, 2, 3, 5, 6, 13, 14], [1, 2, 3, 5, 6, 15, 16], [1, 2, 3, 5, 8, 11, 14], [1, 2, 3, 5, 8, 11, 16],
[1, 2, 3, 5, 9, 10, 11], [1, 2, 3, 5, 9, 12, 16], [1, 2, 3, 5, 9, 13, 15], [1, 2, 5, 6, 7, 8, 13],
[1, 2, 5, 6, 7, 8, 14], [1, 2, 5, 6, 7, 9, 10], [1, 2, 5, 6, 7, 9, 11], [1, 2, 5, 6, 7, 9, 15],
[1, 2, 5, 6, 7, 10, 15], [1, 2, 5, 6, 7, 10, 16], [1, 2, 5, 6, 7, 11, 13], [1, 2, 5, 6, 7, 12, 14],
[1, 2, 5, 6, 8, 13, 14], [1, 2, 5, 6, 10, 11, 15], [1, 2, 5, 6, 11, 14, 16]];

d2:=dev(g,bb2);

d5:=dev(g,bb5);

g1:=Group((2,3,4)(5,16,7,12,9,14)(6,11,8,13,10,15), (1,5,13)(2,8,6)(3,12,9)(4,15,16)(7,11,10));

bb3:=[[1, 2, 3, 4, 5, 6, 13], [1, 2, 3, 4, 5, 6, 15], [1, 2, 3, 4, 5, 7, 10], [1, 2, 3, 5, 6, 7, 14],
[1, 2, 3, 5, 6, 8, 10], [1, 2, 3, 5, 6, 8, 11], [1, 2, 3, 5, 6, 8, 15], [1, 2, 3, 5, 6, 9, 11],
[1, 2, 3, 5, 6, 9, 12], [1, 2, 3, 5, 6, 10, 15], [1, 2, 3, 5, 6, 10, 16], [1, 2, 3, 5, 6, 11, 12],
[1, 2, 3, 5, 6, 12, 13], [1, 2, 3, 5, 6, 14, 15], [1, 2, 3, 5, 6, 15, 16], [1, 2, 3, 5, 7, 8, 9],
[1, 2, 3, 5, 7, 8, 14], [1, 2, 3, 5, 7, 9, 13], [1, 2, 3, 5, 7, 10, 12], [1, 2, 3, 5, 7, 11, 13],
[1, 2, 3, 5, 8, 11, 13], [1, 2, 3, 5, 8, 11, 16], [1, 2, 3, 5, 8, 12, 16], [1, 2, 3, 6, 7, 8, 15],
[1, 2, 5, 6, 7, 8, 9], [1, 2, 5, 6, 7, 8, 11], [1, 2, 5, 6, 7, 10, 15], [1, 2, 5, 6, 7, 11, 12],
[1, 2, 5, 6, 7, 12, 13], [1, 2, 5, 6, 7, 13, 14], [1, 2, 5, 6, 7, 14, 15], [1, 2, 5, 6, 8, 9, 10],
[1, 2, 5, 6, 10, 11, 12], [1, 2, 5, 6, 11, 14, 16]];

bb4:=[[1, 2, 3, 4, 5, 6, 9], [1, 2, 3, 4, 5, 6, 12], [1, 2, 3, 4, 5, 6, 16], [1, 2, 3, 4, 5, 7, 9],
[1, 2, 3, 4, 5, 7, 15], [1, 2, 3, 5, 6, 7, 9], [1, 2, 3, 5, 6, 7, 13], [1, 2, 3, 5, 6, 7, 14],
[1, 2, 3, 5, 6, 8, 9], [1, 2, 3, 5, 6, 8, 11], [1, 2, 3, 5, 6, 8, 13], [1, 2, 3, 5, 6, 8, 15],
[1, 2, 3, 5, 6, 9, 11], [1, 2, 3, 5, 6, 9, 12], [1, 2, 3, 5, 6, 10, 12], [1, 2, 3, 5, 6, 10, 15],
[1, 2, 3, 5, 6, 10, 16], [1, 2, 3, 5, 6, 11, 14], [1, 2, 3, 5, 6, 11, 16], [1, 2, 3, 5, 6, 12, 14],
[1, 2, 3, 5, 6, 13, 16], [1, 2, 3, 5, 6, 14, 15], [1, 2, 3, 5, 7, 8, 12], [1, 2, 3, 5, 7, 9, 15],
[1, 2, 3, 5, 7, 11, 12], [1, 2, 3, 5, 7, 13, 15], [1, 2, 3, 5, 8, 11, 13], [1, 2, 3, 5, 8, 14, 16],
[1, 2, 3, 5, 9, 10, 15], [1, 2, 3, 5, 9, 12, 16], [1, 2, 3, 5, 9, 13, 15], [1, 2, 3, 6, 7, 8, 15],
[1, 2, 3, 6, 7, 13, 14], [1, 2, 5, 6, 7, 8, 10], [1, 2, 5, 6, 7, 8, 13], [1, 2, 5, 6, 7, 9, 10],
[1, 2, 5, 6, 7, 10, 16], [1, 2, 5, 6, 7, 11, 12], [1, 2, 5, 6, 7, 12, 13], [1, 2, 5, 6, 7, 14, 15],
[1, 2, 5, 6, 8, 11, 16], [1, 2, 5, 6, 9, 10, 12], [1, 2, 5, 6, 9, 14, 15], [1, 2, 5, 6, 10, 11, 12],
[1, 2, 5, 6, 10, 11, 15], [1, 2, 5, 6, 10, 13, 14], [1, 2, 5, 6, 11, 13, 14], [1, 2, 5, 6, 11, 14, 16]];

d3:=dev(g1,bb3);

d4:=dev(g1,bb4);