#!usr/bin/perl
open(IN, "ALL.chr16.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf");
open (OUT, ">HMOX2_1000genomes1");
@a = ();
while (<IN>){
        chomp $_;
	$_ =~m/^\d+\t(\d+)\t.*/;
	if ($1 >= 4524719 && $1 <= 4560348){
                print OUT $_,"\n";
        }
}