We developed a dedicated statistical test for a massive detection of spot- A nd facula-crossing anomalies in multiple exoplanetary transit light curves, based on the frequentist p-value thresholding. This test was used to augment our algorithmic pipeline for transit light curves 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 due 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.