I'm using nulcear genomic data from diploid organism. I want to build a phylogenetic tree by concatenating all sequences. One of the first steps is to estimate the model of molecular evolution for the concatenated sequences, for instance by using ModelTest-NG. However, I do not know how to deal wth heterozygotic sites.
Shall mark them as ambiguous sites or shall I use simply the consensus sequence?