We prove a theorem describing the limiting fine-scale statistics of orbits of a point in hyperbolic space under the action of a discrete subgroup. Similar results have been proved only in the lattice case with two recent infinite-volume exceptions by Zhang for Apollonian circle packings and certain Schottky groups. Our results hold for general Zariski dense, non-elementary, geometrically finite subgroups in any dimension. Unlike in the lattice case orbits of geometrically finite subgroups do not necessarily equidistribute on the whole boundary of hyperbolic space. But rather they may equidistribute on a fractal subset. Understanding the behavior of these orbits near the boundary is central to Patterson–Sullivan theory and much further work. Our theorem characterises the higher order spatial statistics and thus addresses a very natural question. As a motivating example our work applies to sphere packings (in any dimension) which are invariant under the action of such discrete subgroups. At the end of the paper we show how this statistical characterization can be used to prove convergence of moments and to write down the limiting formula for the two-point correlation function and nearest neighbor distribution. Moreover we establish a formula for the 2 dimensional limiting gap distribution (and cumulative gap distribution) which also applies in the lattice case.