ASTRAL-Pro on ACCESS1.16Species tree reconstructionChao Zhang and Siavash MirabChao Zhang and Siavash Mirarab (2022) ASTRAL-Pro 2: Ultrafast species tree reconstruction from multi-copy gene family trees. Bioinformatics 38(21) 4949–4950, https://doi.org/10.1093/bioinformatics/btac620Phylogeny / Alignmentastral_pro_accessastralpro_invoke1perl$num_gtrees < 10000 && $num_taxa < 300perl"astral-pro_1.16_expanse -t 16"0astralpro_invoke2perl$num_gtrees < 10000 && $num_taxa > 299 && $num_taxa < 2000perl"astral-pro_1.16_expanse -t 48 "0astralpro_invoke3perl$num_gtrees < 10000 && $num_taxa > 1999perl"astral-pro_1.16_expanse -t 4"0astralpro_invoke4perl$num_gtrees > 9999perl"astral-pro_1.16_expanse -t 8"0number_cores12scheduler.confperl$num_gtrees < 10000 && $num_taxa < 300 perl
"ChargeFactor=1.0\\n" .
"nodes=1\\n" .
"mem=32G\\n" .
"node_exclusive=0\\n" .
"cpus-per-task=16\\n" .
"threads_per_process=16\\n"
number_cores22scheduler.confperl$num_gtrees < 10000 && $num_taxa > 299 && $num_taxa < 2000 perl
"ChargeFactor=1.0\\n" .
"nodes=1\\n" .
"mem=96G\\n" .
"node_exclusive=0\\n" .
"cpus-per-task=48\\n" .
"threads_per_process=48\\n"
number_cores32scheduler.confperl$num_gtrees < 10000 && $num_taxa > 1999 perl
"ChargeFactor=1.0\\n" .
"nodes=1\\n" .
"mem=8G\\n" .
"node_exclusive=0\\n" .
"cpus-per-task=4\\n" .
"threads_per_process=4\\n"
number_cores22scheduler.confperl$num_gtrees > 9999 perl
"ChargeFactor=1.0\\n" .
"nodes=1\\n" .
"mem=16G\\n" .
"node_exclusive=0\\n" .
"cpus-per-task=8\\n" .
"threads_per_process=8\\n"
infileInputgenetrees.tre1infile_invokeInputperl"-i genetrees.tre"99runtime1scheduler.confMaximum Hours to Run (up to 168 hours)0.5The maximum hours to run must be less than 168perl$runtime > 168.0 The maximum hours to run must be greater than 0.05perl$runtime < 0.05perl"runhours=$value\\n"The job will run on 16 processors as configured. If it runs for the entire configured time, it will consume 16 x $runtime cpu hoursperl$num_gtrees < 10000 && $num_taxa < 300 The job will run on 48 processors as configured. If it runs for the entire configured time, it will consume 48 x $runtime cpu hoursperl$num_gtrees < 10000 && $num_taxa > 299 && $num_taxa < 2000 The job will run on 4 processors as configured. If it runs for the entire configured time, it will consume 4 x $runtime cpu hoursperl$num_gtrees < 10000 && $num_taxa > 1999 The job will run on 8 processors as configured. If it runs for the entire configured time, it will consume 8 x $runtime cpu hoursperl$num_gtrees > 9999 Estimate the maximum time your job will need to run. We recommend testing initially with a < 0.5hr test run because Jobs set for 0.5 h or less depedendably run immediately in the "debug" queue.
Once you are sure the configuration is correct, you then increase the time. The reason is that jobs > 0.5 h are submitted to the "normal" queue, where jobs configured for 1 or a few hours times may
run sooner than jobs configured for the full 168 hours.
num_gtreesNumber of Gene trees in your input filePlease specify the number of gene trees in your sampleperl!defined $num_gtreesnum_taxaNumber of taxa in your input filePlease specify the number of unique taxa in your samplesperl!defined $num_taxaspecify_searchroundsSet the number of search rounds (-r)perl(defined $value) ? "-r $value":"" 48specify_subsamplesSet the number of subsamplings (-s)perl(defined $value) ? "-s $value":"" 48mapping_fileSelect a mapping file (-a)mapping_file.txtperl(defined $value) ? "-a multitree_genename.map":"" 8outfile_nameOutput file name perl(defined $value) ? "-o $value":"" multitree.nw.outputPlease specify a name for your outputfileperl!defined $outfile_name8Specify a filename for storing the output species tree.resolved_speciestreePlace taxa on this fully resolved species tree (-c)perl(defined $resolved_speciestree) ? "-c resolved_speciestree.tre":"" resolved_speciestree.tre18resolved_speciestree2Score and annotate the fully resolved species tree (-C)perl( $value) ? "-C ":"" 8This option scores the full species tree file after `-c` without exploring other topologies (`-r 1 -s 0`)specify_seedSpecify a seed value (-seed)perl(defined $value) ? "--seed $value":"" 12Seed for pseudorandomnessspecify_rootRoot at this species (-root)perl(defined $value) ? "--root $value":"" 12Root at the given speciesspecify_supportoptOutput support option (-u)0123perl(defined $value) ? "-u $value":"" 110Select No output support if you want to compute species tree topology only; Select Detailed if you support and local-PP for all three resolutions of each branch..specify_downweightrepeatNumber of trees sampled for each locus. (-w)perl(defined $value) ? "-w $value":"" 112specify_lambdaSet the lambda parameter for the Yule prior. (-l)perl(defined $value) ? "-l $value":"" 0.512Set the lambda parameter for the Yule prior used in the calculations of
branch lengths and posterior probabilities. Set to zero to get ML branch
lengths instead of MAP. Higher values tend to shorten estimated branch
lengths and very high values can give inaccurate results (or even result
in underflow). (default: 0.5)select_constraintfileSelect a constraint file30constraint.treperl(defined $value) ? "-c constraint.tre":"" Newick file containing a binary species tree to place missing species onselect_guidetreefileSelect a guidetree file30guidetree.treperl(defined $value) ? "-g guidetree.tre":"" Newick file containing binary trees as guide trees results*