User:Lindenb/Notebook/UMR915/20101203
From OpenWetWare
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