We study new supergravity solutions related to large-N c N=1 supersymmetric gauge field theories with a large number N f of massive flavors. We use a recently proposed framework based on configurations with N c color D5 branes and a distribution of N f flavor D5 branes, governed by a function N f S(r). Although the system admits many solutions, under plausible physical assumptions the relevant solution is uniquely determined for each value of x ≡ N f /N c . In the IR region, the solution smoothly approaches the deformed Maldacena-Núñez solution. In the UV region it approaches a linear dilaton solution. For x < 2 the gauge coupling β g function computed holographically is negative definite, in the UV approaching the NSVZ β function with anomalous dimension γ 0 = −1/2 (approaching − 3/(32π 2)(2N c  − N f )g 3)), and with β g  → −∞ in the IR. For x = 2, β g has a UV fixed point at strong coupling, suggesting the existence of an IR fixed point at a lower value of the coupling. We argue that the solutions with x > 2 describe a"Seiberg dual" picture where N f  − 2N c flips sign.