We developed a dedicated statistical test for a massive detection of spot- and facula-crossing anomalies in multiple exoplanetary transit lightcurves, based on the frequentist $p$-value thresholding. This test was used to augment our algorithmic pipeline for transit lightcurves analysis. It was applied to $1598$ amateur and professional transit observations of $26$ targets being monitored in the EXPANSION project. We detected $109$ statistically significant candidate events revealing a roughly $2:1$ asymmetry in favor of spots-crossings over faculae-crossings. Although some candidate anomalies likely appear non-physical and originate from systematic errors, such asymmetry between negative and positive events should indicate a physical difference between the frequency of star spots and faculae. Detected spot-crossing events also reveal positive correlation between their amplitude and width, possibly owed to spot size correlation. However, the frequency of all detectable crossing events appears just about a few per cent, so they cannot explain excessive transit timing noise observed for several targets.