PAUP can perform maximum-likelihood estimations on phylogenetic trees. You can save all output trees to a file using:

savetree file=output.phy format=phylip

where format=phylip is optional.

If you want to save a subset of the trees, add the arguments from and to; for example, from=1 to=1 will save the best tree only