Good morning!

I would like to build neighbor-joining tree, but can't decide what nucleotide substitution model should be used.

In MEGA I choose Models > Find best DNA/Protein models (ML) (yes, ML, I would like to use NJ, but futher there is NJ option? ->), then Tree to use - Automatic (NJ tree) and obtain table with GTR, GTR+I, GTR+G, GTR+G+I, K2, K2+G, K2+I, K2+G+I, JC, JC+I, JC+G, JC+G+I, T92, T92+G, T92+I, T92+G+I, HKY, TN93+I, TN93+G, TN93, HKY+I, HKY+G, TN93+G+I, HKY+G+I, whereas in building phylogeny window (NJ) there are such options: No. of differences, P-distance, Jukes-Cantor model (JC), Kimura 2-parameter model (K2), Tajama-Nei model (where is estimation?), Tamura 3-parameter model (T92), Tamura-Nei model (where is estimation?), Maximum Composite Likelihood (where is estimation?), LogDet (Tamura-Kumar) (where is estimation?).

I have also tried modeltest (phangorn), but can't find estimation for Tajama-Nei model, Tamura 3-parameter model, Maximum Composite Likelihood, LogDet (Tamura-Kumar), and can't determine Tamura-Nei = TrN or TrNe in it's output (separate post).

Could you please help me to choose the nucleotide substitution model for NJ?

Thank you in advance!

Best regards, Poecile