Malaria transmission in Colombia is highly variable in space and time. Using a species distribution model, we mapped potential distribution of five vector species including Anopheles albimanus, Anopheles calderoni, Anopheles darlingi, Anopheles neivai, and Anopheles nuneztovari in five Departments of Colombia where malaria transmission remains problematic. We overlaid the range maps of the five species to reveal areas of sympatry and related per-pixel species richness to mean annual parasite index (API) for 2011-2014 mapped by municipality (n = 287). The relationship between mean number of vector species per municipality and API was evaluated using a Poisson regression, which revealed a highly significant relationship between species richness and API (p = 0 for Wald Chi-Square statistic). The results suggest that areas of relatively high transmission in Colombia typically contain higher number of vector species than areas with unstable transmission and that future elimination strategies should account for vector species richness.