We establish global-in-time existence and non-uniqueness of probabilistically strong solutions to the three dimensional Navier-Stokes system driven by space-time white noise. In this setting, solutions are expected to have space regularity at most −1/2 − κ for any κ > 0. Consequently, the convective term is ill-defined analytically and probabilistic renormalization is required. Up to now, only local well-posedness has been known. With the help of paracontrolled calculus we decompose the system in a way which makes it amenable to convex integration. By a careful analysis of the regularity of each term, we develop an iterative procedure which yields global non-unique probabilistically strong paracontrolled solutions. Our result applies to any divergence free initial condition inand implies also non-uniqueness in law. Contents 1. Introduction 2 1.1. Singular SPDEs 3 1.2. Convex integration 5 1.3. Decomposition 6 1.4. Final remarks 8 2. Preliminaries 9 2.1. Function spaces 9 2.2. Paraproducts, commutators and localizers 10 3. More regular stochastic perturbations 12 3.1. The case of γ ∈ (1/2, 3/2] 13 3.2. The case of γ ∈ (1/6, 1/2] 13 3.3. The case of γ ∈ (0, 1/6] 14 4. Paracontrolled solutions 15 4.1. Stochastic objects 15 Date: December 30, 2021. 2010 Mathematics Subject Classification. 60H15; 35R60; 35Q30. Key words and phrases. stochastic Navier-Stokes equations, probabilistically strong solutions, paracontrolled calculus, non-uniqueness in law, convex integration.This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 949981). The financial support by the DFG through the CRC 1283 "Taming uncertainty and profiting from randomness and low regularity in analysis, stochastics and their applications" is greatly acknowledged. R.Z. is grateful to the financial supports of the NSFC (No. 11922103).