Implementation of long‐lasting insecticide‐treated net (LLIN) programs to control human malaria transmission leads to substantial reductions in the abundance of Anopheles mosquitoes, but the impact on the population genetic structure of the malaria vectors is poorly known, nor has it been investigated in Papua New Guinea, where malaria is highly endemic and where several species of Anopheles have vector roles. Here, we applied Wright's F‐statistic, analysis of molecular variance, Bayesian structure analysis, and discriminant analysis of principle components to microsatellite genotype data to analyze the population genetic structure of Anopheles farauti between and within the northern and southern lowland plains and of Anopheles punctulatus within the northern plain of Papua New Guinea after such a program. Bottleneck effects in the two malaria vectors were analyzed using Luikart and Cornuet's tests of heterozygosity. A large, panmictic population of An. punctulatus pre‐LLIN program diverged into two subregional populations corresponding to Madang and East Sepik provinces post‐LLIN distribution and experienced a genetic bottleneck during this process. By contrast, the An. farauti population existed as two regional populations isolated by mountain ranges pre‐LLIN, a genetic structure that persisted after the distribution of LLINs with no further geographic differentiation nor evidence of a genetic bottleneck. These findings show the differential response of populations of different vector species to interventions, which has implications for program sustainability and gene flow.