We present an extraction of the unpolarized transverse-momentum-dependent parton distribution and fragmentation functions that takes into account possible differences between quark flavors and final-state hadrons. The extraction is based on experimental measurements from Drell-Yan processes and semi-inclusive deep-inelastic scattering, whose combination is essential to distinguish flavor differences. The analysis is carried out at N3LL accuracy. The extracted flavor-dependent distributions give a very good description of the data (χ2/Ndat = 1.08). The resulting error bands take fully into account also the uncertainties in the determination of the corresponding collinear distributions.