Hardness of cooked soybeans [Glycine max (L). Merr.] is an important attribute in food processing. We found one candidate gene, Glyma03g03360, to be associated with the hardness of cotyledons of cooked soybeans, based on a quantitative trait locus and fine-scale mapping analyses using a recombinant inbred line population developed from a cross between two Japanese cultivars, "Natto-shoryu" and "Hyoukei-kuro 3". Analysis of the DNA sequence of Glyma03g03360, a pectin methylesterase gene homologue, revealed three patterns of mutations, two of which result in truncated proteins and one of which results in an amino acid substitution. The truncated proteins are presumed to lack the enzymatic activity of Glyma03g03360. We classified 24 cultivars into four groups based on the sequence of Glyma03g03360. The texture analysis using the 22 cultivars grown in different locations indicated that protein truncation of Glyma03g03360 resulted in softer cotyledons of cooked soybeans, which was further confirmed by texture analysis performed using F2 populations of a cross between "Enrei" and "LD00-3309", and between "Satonohohoemi" and "Sakukei 98". A positive correlation between hardness and calcium content implies the possible effect of calcium binding to pectins on the hardness of cooked soybean cotyledons.