Code to calculate species trees from large genomic datasets