LineStacker is a new open access and open source tool for stacking of spectral lines in interferometric data. LineStacker is an ensemble of CASA tasks, and can stack both 3D cubes or already extracted spectra. The algorithm is tested on increasingly complex simulated data sets, mimicking Atacama Large Millimeter/submillimeter Array and Karl G. Jansky Very Large Array observations of [CII] and CO(3-2) emission lines, from $z\sim7$ and $z\sim4$ galaxies respectively. We find that the algorithm is very robust, successfully retrieving the input parameters of the stacked lines in all cases with an accuracy $\gtrsim90$\%. However, we distinguish some specific situations showcasing the intrinsic limitations of the method. Mainly that high uncertainties on the redshifts ($\Delta z > 0.01$) can lead to poor signal to noise ratio improvement, due to lines being stacked on shifted central frequencies. Additionally we give an extensive description of the embedded statistical tools included in LineStacker: mainly bootstrapping, rebinning and subsampling. Velocity rebinning {is applied on the data before stacking and} proves necessary when studying line profiles, in order to avoid artificial spectral features in the stack. Subsampling is useful to sort the stacked sources, allowing to find a subsample maximizing the searched parameters, while bootstrapping allows to detect inhomogeneities in the stacked sample. LineStacker is a useful tool for extracting the most from spectral observations of various types.

Outflows in quasars during the early epochs of galaxy evolution are an important part of the feedback mechanisms potentially affecting the evolution of the host galaxy. However, systematic observations of outflows are only now becoming possible with the advent of sensitive mm telescopes. In this study we use spectral stacking methods to search for faint high velocity outflow signal in a sample of [C II] detected, z ~ 6 quasars. We search for broad emission line signatures from high-velocity outflows for a sample of 26 z ~ 6 quasars observed with ALMA, with a detection of the [C II] line. The observed emission lines of the sources are dominated by the host galaxy, and outflow emission is not detected for the individual sources. We use a spectral line stacking analysis developed for interferometric data to search for outflow emission. We stack both extracted spectra and the full spectral cubes. We also investigate the possibility that only a sub-set of our sample contributes to the stacked outflow emission. We find only a tentative detection of a broad emission line component in the stacked spectra. When taking a region of about 2 arcsec around the source central position of the stacked cubes, the stacked line shows an excess emission due to a broad component of 1.1-1.5 sigma, but the significance drops to 0.4-0.7 sigma when stacking the extracted spectra from a smaller region. The broad component can be characterised by a line width of full width half max FWHM > 700 km/s. Furthermore, we find a sub-sample of 12 sources the stack of which maximises the broad component emission. The stack of this sub-sample shows an excess emission due to a broad component of 1.2-2.5 sigma. The stacked line of these sources has a broad component of FWHM > 775 km/s. Deeper ALMA observations are necessary to confirm the presence of a broad component in the individual spectra.