To clarify the neural systems deployed by individual subjects during working memory (WM), we collected functional neuroimaging data from healthy subjects, and constructed a model of 2-back WM using structural equation modeling (SEM). A group model was constructed, and models for each subject were validated against it. The group model consisted principally of regions in the prefrontal and parietal cortex, with considerable interindividual variance in the single-subject models. To explore this variance, subjects were split into two groups based on performance. Performance level and self-reported strategy scores were used in a correlation analysis against path weights between nodes of individual models. High performers utilized a left hemisphere sub-network involving inferior parietal lobule and Broca's area, whereas lower performers utilized a right hemisphere sub-network with interactions between inferior parietal lobule and dorsolateral prefrontal cortex. Further, we observed an interaction between the parahippocampal formation and the inferior parietal lobule that was related to the different strategies used by the individuals to perform the task. Strategy and performance level appear to be intricately related in this task, with neural systems supporting verbal processing producing better performance than those associated with spatial processing. These results demonstrate that individual behavioral characteristics are reflected in specific neurofunctional patterns at the system level and that these can be captured by analytical techniques such as SEM.