This paper investigates the numerical accuracy of implicit Large Eddy Simulations (iLES) in relation to compressible turbulent boundary layers (TBL). iLES are conducted in conjunction with Monotonic Upstream-Centred Scheme for Conservation Laws (MUSCL) and Weighted Essentially Non-Oscillatory (WENO), ranging from 2nd to 9th-order. The accuracy effects are presented from a physical perspective showing skewness, flatness and anisotropy calculations, among others. The order of the scheme directly affects the physical representation of the TBL, especially the degree of asymmetry and anisotropy in the sub-layers of the TBL. The study concludes that high-order iLES can provide an accurate and detailed description of TBL directly comparable to available DNS and experimental results.