Given the general composition of plastomes (coding regions with triplett structure, tRNA and rRNA genes, non-coding regions, I wonder if a distance-based tree (how were the distances calculated?) provides branch lengths that allow proper comparison between the main clades.

Maybe you can use the matrix you used to infer pairwise distances for a quick, partitioned ML analysis (e.g. with ExaML: ht...