Pattern formation in biological aggregations is a topic of great interest, due to the complex spatial structure of various aggregations of cells/bacteria/animals that can be observed in nature. While many such aggregations look similar at the macroscopic level, they might differ in their microscopic spatial structure. However, the complexity of the non-linear and sometimes non-local interactions among individuals inside these aggregations makes it difficult to investigate these spatial structures. In this study, we investigate numerically the transitions between different spatial patterns of animal aggregations with various symmetries (even, odd or no symmetry) that characterise the microscopic distribution of individuals inside these aggregations. To this end, we construct a bifurcation diagram starting with perturbations of spatially homogeneous solutions with low, medium, and high amplitudes. For perturbations with low amplitudes, the bifurcating structures show transitions among even-symmetric, odd-symmetric, and non-symmetric solutions. For perturbations with large amplitudes, there are wide parameter regions with non-convergent solutions, characterised by oscillatory transitions between different relatively similar solutions. These numerical results emphasize: (i) the effect of nonlinear and non-local interactions on the microscopically different symmetric/non-symmetric structures of macroscopically similar ecological aggregations; (ii) the difficulty of developing continuation algorithms for this class of non-local models.