calculate Tajima'D on haploid about vcf-kit

[danielecook]

Hi @AlexWanghaoming - I would recommend simply replacing your genotype calls from 0 or 1 to 0/0 and 1/1 in the VCF. Then you should be able to run VCF kit with a ploidy state of 1. I'm not sure whether this will affect your caculation - but I do not think it will.

[Rohit-Satyam]

Hi @danielecook, my VCF file contains genotype calls in 0/0 format. I got it from MalariaGen and now I was trying to calculate the Tajima's D (TD) for the two most polymorphic genes in Plasmodium (which is haploid) and yet I get negative TD values for them. I am not sure if the subcommand vk tajima is doing the right thing because highly polymorphic gene should have high TD values

Here are the subsetted VCF files for the same genes:
File1
File2

PF3D7_0930300 Pf3D7_09_v3:1201305-1207576
PF3D7_1133400 Pf3D7_11_v3:1292966-1296696