Burkholderia gladioli pv. cocovenenans is a serious safety issue in black fungus due to the deadly toxin, bongkrekic acid. This has triggered the demand for an efficient toxigenic phenotype recognition method. The objective of this study is to develop an efficient method for the recognition of toxin-producing B. gladioli strains. The potential of multilocus sequence typing and a back propagation neural network for the recognition of toxigenic B. cocovenenans was explored for the first time. The virulent strains were isolated from a black fungus cultivation environment in Qinba Mountain area, Shaanxi, China. A comprehensive evaluation of toxigenic capability of 26 isolates were conducted using Ultra Performance Liquid Chromatography for determination of bongkrekic acid and toxoflavin production in different culturing conditions and foods. The isolates produced bongkrekic acid in the range of 0.05–6.24 mg/L in black fungus and a highly toxin-producing strain generated 201.86 mg/L bongkrekic acid and 45.26 mg/L toxoflavin in co-cultivation with Rhizopus oryzae on PDA medium. Multilocus sequence typing phylogeny (MLST) analysis showed that housekeeping gene sequences have a certain relationship with a strain toxigenic phenotype. We developed a well-trained, back-propagation neutral network for prediction of toxigenic phenotype in B. gladioli based on MLST sequences with an accuracy of 100% in the training set and an accuracy of 86.7% in external test set strains. The BP neutral network offers a highly efficient approach to predict toxigenic phenotype of strains and contributes to hazard detection and safety surveillance.