User:Lindenb/Notebook/UMR915/20101203

cleaning day.

belgium
get 1000G data: curl -s --proxy cache.XXXXXXXXXXX:1111 "ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/2010_11/ALL.2of4intersection.20100804.sites.vcf.gz" | gunzip  -c |\ awk -F '   ' '{if(($1=="X1" && int($2)>=1 && int($2) <=2) || ($1=="X2" && int($2) >= 1  && int($2)<=2)) print $0;}'  > /home/lindenb/belgium.1000G.vcf

wc /home/lindenb/belgium.1000G.vcf 1810743 14485944 112524319 /home/lindenb/belgium.1000G.vcf