  Galactic ultra compact binaries are expected to be the dominant source of
gravitational waves in the milli-Hertz frequency band. Of the tens of millions
of galactic binaries with periods shorter than an hour, it is estimated that a
few tens of thousand will be resolved by the future Laser Interferometer Space
Antenna (LISA). The unresolved remainder will be the main source of ``noise''
between 1-3 milli-Hertz. Typical galactic binaries are millions of years from
merger, and consequently their signals will persist for the the duration of the
LISA mission. Extracting tens of thousands of overlapping galactic signals and
characterizing the unresolved component is a central challenge in LISA data
analysis, and a key contribution to arriving at a global solution that
simultaneously fits for all signals in the band. Here we present an end-to-end
analysis pipeline for galactic binaries that uses trans-dimensional Bayesian
inference to develop a time-evolving catalog of sources as data arrive from the
LISA constellation.
