The eternal inflation scenario predicts that our observable universe resides inside a single bubble embedded in a vast inflating multiverse. Collisions between bubble universes imprinted in the CMB sky provide a powerful observational test of this idea. I will describe a robust algorithm for non-Gaussian source detection in massive datasets, and present its application to the search for bubble collision signatures in CMB data from WMAP.