Statistical Inference of Phylogenetic Networks from Unlinked Bi-allelic Markers