We compute the expectation values of BPS Wilson loops in the mass-deformed ABJM theory using the Fermi gas technique. We obtain explicit results in terms of Airy functions, effectively resumming the full 1/N expansion up to exponentially small terms. These expressions enable us to derive multi-point correlation functions for topological operators belonging to the stress tensor multiplet, in the presence of a 1/2-BPS Wilson line. From the one-point correlator, we recover the ABJM Bremsstrahlung function, confirming nicely previous results obtained through latitude Wilson loops. Likewise, higher point correlators can be used to extract iteratively new defect CFT data for higher dimensional topological operators. We present a detailed example of the dimension-two operator appearing in the OPE of two stress tensor multiplets.