We first compute m trees, each on n−2 leaves using the double anchored method (Corollary 1) and then combine these m trees using a supertree method. Using a compatibility supertree (i.e., one that given a set of compatible input trees, outputs a tree that refines all input trees) would make the approach statistically consistent (Theorem S2, Additional file 1).