Constructing the full network

Finally, putting the full network together with InPhyNet is simple:

species_net = inphynet(D, snaq_networks, namelist)

Our example has an outgroup taxa, so we can root the network at that outgroup.

rootatnode!(species_net, "OUTGROUP")

Now, we can use PhyloPlots.jl to plot the species network.

using PhyloPlots
plot(species_net);