This article investigates the source identification in the fractional diffusion equations, by performing a single measurement of the Cauchy data on the accessible boundary. The main results of this work consist in giving an identifiability result and establishing a local Lipschitz stability result. To solve the inverse problem of identifying fractional sources from such observations, a non-iterative algebraical method based on the Reciprocity Gap functional is proposed.