We investigate the average inter-residue folding forces derived from mutational data of the 15 proteins: barstar, barnase, chymotrypsin inhibitor 2 (CI2), Src SH3 domain, spectrin R16 domain, Arc repressor, apo-azurin, cold shock protein B (cspB), C-terminal domain of ribosomal protein L9 (CTL9), FKBP12, alpha-lactalbumin, colicin E7 immunity protein 7 (IM7), colicin E9 immunity protein 9 (IM9), spectrin R17 domain, and ubiquitin. The residue-specific contributions to folding in most of the 15 protein molecules are highly non-uniformly distributed and are typically about 1piconewton (pN) per interaction. The strongest folding forces often occur in some of the helices and strands of folding nuclei which suggests that folding nucleation-condensation is partially directed by formation of some secondary structure interactions. The correlation of the energy changes of mutants with inter-residue contact maps of the protein molecules provides a higher resolution than assigning the mutant data to certain positions in the polypeptide strand alone. In contrast to previous Phi-value analysis, we now can partially resolve folding motions. Compaction of at least one alpha-helix along its axis mediated by internal hydrogen bonds and stabilized by diffuse tertiary structure interactions appears to be one important molecular event during early folding in barstar, CI2, spectrin R16 domain, Arc repressor, alpha-lactalbumin, IM7, IM9, and spectrin R17 domain. A lateral movement of at least two strands neighbored in sequence towards each other appears to be involved in early folding of the SH3 domain, cspB, CTL9, and FKBP12.