library(ape)
library(seqinr)
data<- read.dna('PrimateData.fasta.txt', format="fasta")
dist.dna(data, model="TN93", gamma=1)