A rubber composition in which a relationship between: a maximum value Ξ c of a correlation length Ξ i of 10 nm to 100 μm obtained by performing curve-fitting with the following (Expression 1-2) to (Expression 1-7) on a scattering intensity curve I(q) obtained by performing X-ray scattering measurement or neutron scattering measurement in a region of q represented by the following (Expression 1); and a standard deviation σ c of the maximum value Ξ c is σ c ≤ Ξ c /2, and in which a filler distributivity index D obtained by the following (Expression 1-8) is not greater than 2, and wherein an amount of carbon black with respect to 100 parts by mass of the rubber component is not less than 1 part by mass and is not greater than 20 parts by mass, [Math. 1] q = 4 ⁢ π sin θ / 2 λ (θ: scattering angle, λ: a wavelength of X-rays or neutron rays) [Math. 2] I q = A 1 + q 2 ξ 2 g a r dr + ∑ i = 1 n ∫ B i 1 + q 2 Ξ i 2 2 g bi r dr ξ < Ξ i A = 8 πN a ρ 2 ξ 3 B i = 4 πN i ρ 2 Ξ i 2 g a r = 1 2 πσ a 2 exp − r − ξ 2 2 σ a 2 g bi r = 1 2 πσ i 2 exp − r − Ξ i 2 2 σ i 2 (A, B i , ξ, Ξ i : fitting parameters) (n: an integer) (q: the same as in (Expression 1)) (N a : the number of scatterers having a correlation length ξ per unit volume (number/cm 3 )) (N i : the number of scatterers having the correlation length Ξ i per unit volume (number/cm 3 )) (p: a scattering length density difference (cm -2 ) or an electron density difference (electron cm -3 ) between the scatterers and a deuterated solvent or air therearound) (g a(r) : a probability density function having an average correlation length ξ and a standard deviation σ a when a correlation length is assumed to be present in a Gaussian distribution) (g bi(r) : a probability density function having an average correlation length Ξ i and a standard deviation σ i when a correlation length is assumed to be present in a Gaussian distribution) [Math. 3] D = f A / f B f 0 , A / f 0 , B f 0,A : a volume fraction of a polymer A at time of addition for kneading f 0,B : a volume fraction of a polymer B at time of addition for kneading f A : a volume fraction of the polymer A obtained by gel fraction measurement f B : a volume fraction of the polymer B obtained by gel fraction measurement volume fraction of the polymer A at the time of addition for kneading > volume fraction of the polymer B at the time of addition for kneading. (A, B i , ξ, Ξ i : fitting parameters) (n: an integer) (q: the same as in (Expression 1)) (N a : the number of scatterers having a correlation length ξ per unit volume (number/cm 3 )) (N i : the number of scatterers having the correlation length Ξ i per unit volume (number/cm 3 )) (p: a scattering length density difference (cm -2 ) or an electron density difference (electron cm -3 ) between the scatterers and a deuterated solvent or air therearound) (g a(r) : a probability density function having an average correlation length ξ and a standard deviation σ a when a correlation length is assumed to be present in a Gaussian distribution) (g bi(r) : a probability density function having an average correlation length Ξ i and a standard deviation σ i when a correlation length is assumed to be present in a Gaussian distribution) f 0,A : a volume fraction of a polymer A at time of addition for kneading f 0,B : a volume fraction of a polymer B at time of addition for kneading f A : a volume fraction of the polymer A obtained by gel fraction measurement f B : a volume fraction of the polymer B obtained by gel fraction measurement <EOS>