<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Matthieu<div class="">Thanks very much for the quick reply. I can see that there’s a bunch of parameters that can be passed to the script you mentioned but, as far as I can tell, these are not documented:</div><div class=""><div class=""> 'url=s' => \$url,</div><div class=""> 'mlss_id=s' => \$mlss_id,</div><div class=""> 'method=s' => \$method,</div><div class=""> 'ss_name=s' => \$ss_name,</div><div class=""> 'label=s' => \$label,</div><div class=""> 'stn_root_id=i' => \$stn_root_id,</div><div class=""> 'with_distances' => \$with_distances,</div><div class=""> 'ascii_scale=f' => \$ascii_scale,</div><div class=""> 'reg_conf=s' => \$reg_conf,</div><div class=""> 'compara_db=s' => \$compara_db</div><div class=""><br class=""></div><div class="">I couldn’t find a combination of url, reg_conf and compara_db that would work so I just hardcoded a DBA that seems to work. Then, by experimenting with combinations of species tree root and species set ID listed here: <a href="http://ensembl.org/info/docs/api/compara/compara_schema.html#species_tree_root" class="">http://ensembl.org/info/docs/api/compara/compara_schema.html#species_tree_root</a> I think I got the tree I wanted. I was also able to modify the ‘roll your own’ parameters without any problems.</div><div class=""><br class=""></div><div class="">Thanks again,</div><div class="">Greg</div><div class="">-- </div></div><div class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;">Dr Greg Slodkowicz<br class="">Structural Studies (Babu Group)<br class="">MRC Laboratory of Molecular Biology<br class="">Francis Crick Avenue,<br class="">Cambridge, CB2 0QH, UK</div></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On 5 Dec 2019, at 10:22, Matthieu Muffato <<a href="mailto:muffato@ebi.ac.uk" class="">muffato@ebi.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi Greg<br class=""><br class="">When we run our pipelines, e.g. TreeBest, we use internal node identifiers instead of species names. It's shorter and less ambiguous to parse. It's all part of the pipeline, but you can adapt the ensembl-compara/scripts/examples/species_getSpeciesTree.pl script to get it. Replace the string %{n} with %{o}%{-E"*"} .<br class=""><br class="">"o" tells the formatter to use node IDs, and the star character must be added for TreeBest not to penalise the species. Initially this was meant to cater for the "low-coverage" mammals assembly, but in our tests a few years ago it didn't seem useful any more, so we flag everything as "fully sequenced" by adding the star character.<br class=""><br class="">This script is not compatible with the e78 API / database schema, but I think you can use the latest species-tree and prune it. For protein-trees we use the NCBI taxonomy, and I don't think it's changed much (I remember the relationship between birds, turtles and other reptiles changed at some point, but can't remember when)<br class=""><br class="">Hope this helps,<br class="">Matthieu<br class=""><br class="">On 03/12/2019 11:05, Greg Slodkowicz wrote:<br class=""><blockquote type="cite" class="">Dear developers,<br class="">I was wondering if there is a way to access the species tree that is used for running TreeBEST? I have the ’normal’ Newick species tree from the Compara GitHub repository but it seems like TreeBEST is quite picky about the labelling of the tree nodes.<br class=""><br class="">What I would like to do is re-run gene-species tree reconciliation for a few gene trees of interest, get bootstrap replicates for that tree and then run some downstream analysis on them to get a measure of uncertainty introduced by the differences in the tree.I would ideally like the species tree from an archival (release 78) version of Ensembl (though I imagine it hasn’t changed that much).<br class=""><br class="">Many thanks,<br class="">Greg<br class=""></blockquote><br class=""></div></div></blockquote></div><br class=""></body></html>