1. calculated the wondows Fst of gene and chosen randomly Non-coding region that located at least 200 kb away from known gene
python Generate.Prir.Fst.py > fst.sh
sh fst.sh
sh Datadeal.sh
ls File/*.gene |while read line;do id=${line#*/};id=${id%%.*};Rscript T-test.R $id >> T-test.res;done
Rscript T-test.res.plot.R