This book presents research results concerning the distribution of prime numbers. The first major result discussed is the supremum for the maximal prime gaps. By an implementation of a binomial coefficient the maximal prime gaps supremum bound is proved, simultaneously establishing the infimum for primes in the short interval. Subsequently, a novel application of the theory of the primorial function establishes the tailored logarithmic integral, which is a superior adaptation of the classical Gauss' logarithmic integral. The tailored integral due to its radically improved accuracy over the Gauss' logarithmic integral, constitutes the supremum bound of estimation of the prime counting function. It presents the possibility to estimate the prime counting function with unprecedented accuracy.