Connecting genetic variation to phenotypes and understanding the underlying biological mechanisms has been a fundamental goal of biological genetics. Here, we used the association analysis to identify a Vip3Aa resistance-associated genomic region in a strain of fall armyworm, JC-R, which exhibits >5000-fold resistance to the Bt toxin Vip3Aa. However, through various analytical approaches and fine-scale mapping across different populations, we demonstrated that this genomic region exhibits strong genetic linkage. The chromosome-level genome of JC-R and its parent strain JC-S were assembled, and extensive structural variations in the linkage regions were identified, which could be responsible for maintaining the linkage. To identify the causal variation within this linked region, a chromosome fragment stepwise knockout strategy based on CRISPR/Cas9 was developed. By crossing with the resistant strain and phenotyping segregating offspring on Vip3Aa-containing diet, we identified a chromosomal segment, KO8, containing the resistant gene. Subsequently, we conducted a comprehensive analysis of the variations in the KO8 region using multi-omics approaches, including genomic data, RNA-seq, proteomic, PacBio long read Iso-seq, and phosphoproteomic data. This analysis identified multiple variations in the chitin synthase geneCHS2, including amino acid substitution, alternative splicing, and changes in phosphorylation sites. After knocking out theCHS2, larvae exhibited over 6777-fold resistance to Vip3Aa. These results demonstrate that the chromosome fragment stepwise knockout strategy is a viable approach for studying complex genomic regions, and highlight the value of comprehensive analysis of genetic variations using multi-omics data. The identified candidate gene could potentially advance monitoring and management of pest resistance to Vip3Aa.