  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. In the maximal supersymmetric case, 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.
