A mathematical procedure is presented for the calculation of exact cumulative distribution statistics for a viscosity-free variant of Burgers nonlinear partial differential equation (PDE) in one space dimension and time subject to sinusoidal initial data with uncertain (random variable) amplitude or phase shift. Analytical solutions of nonlinear PDEs with uncertain initial and/or boundary data are invaluable benchmarks in assessing approximate uncertainty quantification techniques. The Burgers equation solution with uncertain initial data results in nonsmooth solution behavior in both physical and random variable dimensions which provides a severe test for approximate uncertainty quantification techniques. Mathematical proofs are provided to verify that exact cumulative distribution statistics can be systematically and robustly obtained for all forward time.