The fluctuation of highly penetrated distributed generations (DGs) in distribution networks (DNs) increases branch overload risks, which makes the analysis uncertainty of branch power flow more complex. The probability analysis of branch power can grasp the power flow operation characteristics under uncertain conditions, which supports the power flow optimization management and ensures the safe and economic operation of DNs. Therefore, this paper proposes a data-model driven probability analysis method for branch power flow in DNs. An approximate branch power model is first derived to reveal the analytical relationship between branch power and node power. Then, based on the branch power approximation model and the forecasting error, the datasets of branch apparent power are constructed to capture the complex nonlinear characteristics of the power flow. The non-Gaussian distribution characteristics of branch apparent power and line loss are described by the optimal probability fitting method. Finally, the probability level of branch overload is evaluated and the probability distribution of line loss is analyzed. The case studies demonstrate that the branch power approximation model is highly accurate, and the probability distribution of the branch apparent power presents different characteristics. The proposed method can quickly and accurately calculate the branch overload probability level.