I am wondering if anyone can help me figure out how to build a neighbor joining tree in R, where the pairwise distances (from sequence alignments) are estimated under the GTR model. I have tried the dist.dna() method from R's phanghorn package, but it does not have support for the GTR model. I am new to phylogenetics, and any help would be great.

This is great. I really appreciate it