Human brain network is organized as interconnected communities for supporting cognition and behavior. Despite studies on the nonoverlapping communities of brain network, overlapping community structure and its relationship to brain function remain largely unknown. With this consideration, we employed the Bayesian nonnegative matrix factorization to decompose the functional brain networks constructed from resting-state fMRI data into overlapping communities with interdigitated mapping to functional subnetworks. By examining the heterogeneous nodal membership to communities, we classified nodes into three classes: Most nodes in somatomotor and limbic subnetworks were affiliated with one dominant community and classified as unimodule nodes; most nodes in attention and frontoparietal subnetworks were affiliated with more than two communities and classified as multimodule nodes; and the remaining nodes affiliated with two communities were classified as bimodule nodes. This three-class paradigm was highly reproducible across sessions and subjects. Furthermore, the more likely a node was classified as multimodule node, the more flexible it will be engaged in multiple tasks. Finally, the FC feature vector associated with multimodule nodes could serve as connectome "fingerprinting" to gain high subject discriminability. Together, our findings offer new insights on the flexible spatial overlapping communities that related to task-based functional flexibility and individual connectome "fingerprinting."