use the dbSNP bzip2 json file and a gtf as input, calculate the Tajima's D per gene over the exon region.
for num in $(seq 1 22) ;
do
./target/release/popgen_dbSNP --gtf ../chr${num}.gtf --dbsnp ../refsnp-chr${num}.json.bz2 --bin 0 --studyname 1000Genomes > chr${num}.popgen_1kG.tsv 2>chr${num}.popgen_1kG.err & ;
done