Connectivity analysis of resting state brain has provided a novel means of investigating brain networks in the study of neurodevelpmental disorders. The study of functional networks, often represented by high dimensional graphs, predicates on the ability of methods in succinctly extracting meaningful representative connectivity information at the subject and population level. This need motivates the development of techniques that can extract underlying network modules that characterize the connectivity in a population, while capturing variations of these modules at the individual level. In this paper, we propose a multi-layer raph clustering technique that fuses the information from a collection of connectivity networks of a population to extract the underlying common network modules that serve as network hubs for the population. These hubs form a functional network atlas. In addition, our technique provides subject-specific factors designed to characterize and quantify the degree of intra- and inter- connectivity between hubs, thereby providing a representation that is amenable to group level statistical analyses. We demonstrate the utility of the technique by creating a population network atlas of connectivity by examining MEG based functional connectivity in typically developing children, and using this to describe the individualized variation in those diagnosed with autism spectrum disorder.