Hi!
I have a set of CO1 sequences (n=18 and 657 sites) for 5 species from the same genus which I have aligned in MEGA X. I want to find out the nucleotide diversity (pi), number of segregating sites (S), parsimony informative sites (PIS) and to calculate Tajimas D. I have tried to use both MEGA and DNAsp to do this but they provide different results. To find the values in DNAsp I uploaded the alignment produced in MEGA as a fasta file. This produced copies of each sequence which i removed (I am also wondering why this is happening). I then used the Overwiev -> Polymorphism Data option. The values I get are presented below.
MEGA X:
- pi = 0.126
- S = 191
- D = 2.10
DNAsp
- pi = 0.138
- S = 182
- D = 0.763
- PIS = 168
Why is this happening? Do I need to calculate the values for each species separately? All help will be greatly appriciated!