We present first principles calculations of the magnetic and orbital properties of Ba2NaOsO6 (BNOO), a 5d 1 Mott insulator with strong spin orbit coupling (SOC) in its low temperature emergent quantum phases. Our computational method takes into direct consideration recent NMR results that established that BNOO develops a local octahedral distortion preceding the formation of long range magnetic order. We found that the two-sublattice canted ferromagnetic ground state identified in Lu et al., Nature Comm. 8, 14407 (2017) is accompanied by a two-sublattice staggered orbital ordering pattern in which the t2g orbitals are selectively occupied as a result of strong spin orbit coupling. The staggered orbital order found here using first principles calculations asserts the previous proposal of Chen et al., Phys. Rev. B 82, 174440 (2010) and Lu et al., Nature Comm. 8, 14407 (2017), that two-sublattice magnetic structure is the very manifestation of staggered quadrupolar order. Therefore, our results affirm the essential role of multipolar spin interactions in the microscopic description of magnetism in systems with locally entangled spin and orbital degrees of freedom.