#This script takes alignment files (BAMS) from pool-seq, subsamples to minimize effects of differential coverage and calculates a genome wide population statistic (Tajima's D) using Popoolation.
#This shell script was used to unzip a zip file of the alignment files and run popoolationStats_na_KL.py on a condor pool.
#This was the submit script created to submit the job and associated files (popoolationStats_na_KL.py, popoolationStats.sh, MtbNCBIH37Rv.fa, RNAseq_##.zip) to the condor pool.
#This is the Rscript used to create figures of the Tajima's D results