#nexus begin paup; cd *; [allows specification of pathnames relative to this file] execute mammals447.nex; include 1-37280/only; [first 10 genes] rm mammal.tre; lset nst=6 rmatrix=estimate rates=gamma shape=estimate basefreq=estimate; showmodel; beginbysubset charpartition=loci; ["loci" was defined in the sequence-data file] raxml; savetrees file=mammal.tre format=phylip append; endbysubset; end;