The "Brownian bees" model describes an ensemble of N independent branching Brownian particles. When a particle branches into two particles, the particle farthest from the origin is eliminated so as to keep the number of particles constant. In the limit of N→∞, the spatial density of the particles is governed by the solution of a free boundary problem for a reaction-diffusion equation. At long times the particle density approaches a spherically symmetric steady-state solution with a compact support. Here, we study fluctuations of the "swarm of bees" due to the random character of the branching Brownian motion in the limit of large but finite N. We consider a one-dimensional setting and focus on two fluctuating quantities: the swarm center of mass X(t) and the swarm radius ℓ(t). Linearizing a pertinent Langevin equation around the deterministic steady-state solution, we calculate the two-time covariances of X(t) and ℓ(t). The variance of X(t) directly follows from the covariance of X(t), and it scales as 1/N as to be expected from the law of large numbers. The variance of ℓ(t) behaves differently: It exhibits an anomalous scaling (1/N)lnN. This anomaly appears because all spatial scales, including a narrow region near the edges of the swarm where only a few particles are present, give a significant contribution to the variance. We argue that the variance of ℓ(t) can be obtained from the covariance of ℓ(t) by introducing a cutoff at the microscopic time 1/N where the continuum Langevin description breaks down. Our theoretical predictions are in good agreement with Monte Carlo simulations of the microscopic model. Generalizations to higher dimensions are briefly discussed.