Improved understanding of the effects of variability in electrophysiological activity within the human heart is key to understanding and predicting cardiovascular response to disease and treatments. Previous studies have considered either regional variation in action potentials or inter-subject variability within one region. We hypothesize that the regional differences in morphology derive, not only from variation in dependence on individual conductance, but also on the relationship between multiple conductances. Using the Monte-Carlo sampling method and the Maleckar model, we created an in-silico population, varying each conductance +/-100%. The population was divided into regional groups based on the biomarkers: APD90, APD50, APD20, RMP, and APA. We observe the significance of the variation of 9 different conductances on the AP morphology for each atrial region. Results showed the significant dependence on various ion channel concentrations differed among atrial regions. GNaK concentration had a significant impact on APD50 for CT/BB, LA and LAA, but not the other atrial regions. CT/BB and LAA were the only regions where gKur lacked a significant impact on the APD20 biomarker. GITo, gKur, gKr, gKs and gK1 all had significant impacts on each atrial region, however the only significant cross-correlation between these ion channels, effecting APD90 was gKur and gk1, (p = 2.3x10^-80). Similar results were found for APD50 and RMP. The total population showed a significant relationship with the cross correlation of gINa and gNaCa for APA (p = 2.3x10^-88 ) and gINa and gK1 for APD50 (p = 4.11x10^-23 ), among others. This shows that not only are there relationships between individual ion channel concentrations and the AP morphology, but also between the combination of various conductances and AP morphology. The regionally based populations are to be applied to a three-dimensional model to determine if intra-regional variability significantly effects wave front propagation and re-entry patterns.