Background: We aim to implement an immune cell score model in routine clinical practice for resected non-small-cell lung cancer (NSCLC) patients (NCT03299478). Molecular and genomic features associated with immune phenotypes in NSCLC have not been explored in detail. Patients and methods: We developed a machine learning (ML)-based model to classify tumors into one of three categories: inflamed, altered, and desert, based on the spatial distribution of CD8+ T cells in two prospective (n = 453; TNM-I trial) and retrospective (n = 481) stage I-IIIA NSCLC surgical cohorts. NanoString assays and targeted gene panel sequencing were used to evaluate the association of gene expression and mutations with immune phenotypes. Results: Among the total of 934 patients, 24.4% of tumors were classified as inflamed, 51.3% as altered, and 24.3% as desert. There were significant associations between ML-derived immune phenotypes and adaptive immunity gene expression signatures. We identified a strong association of the nuclear factor-κB pathway and CD8+ T-cell exclusion through a positive enrichment in the desert phenotype. KEAP1 [odds ratio (OR) 0.27, Q = 0.02] and STK11 (OR 0.39, Q = 0.04) were significantly co-mutated in non-inflamed lung adenocarcinoma (LUAD) compared to the inflamed phenotype. In the retrospective cohort, the inflamed phenotype was an independent prognostic factor for prolonged disease-specific survival and time to recurrence (hazard ratio 0.61, P = 0.01 and 0.65, P = 0.02, respectively). Conclusions: ML-based immune phenotyping by spatial distribution of T cells in resected NSCLC is able to identify patients at greater risk of disease recurrence after surgical resection. LUADs with concurrent KEAP1 and STK11 mutations are enriched for altered and desert immune phenotypes.