chrX 1 1000000
chrX 1000001 2000001
chrX 2000001 3000001
...etc.
####
open (INTERVAL, "/Users/logancurtis-whitchurch/Dropbox/thesis_folder/galaxy_chrX_data/chrX_1Mbwindow_nonoverlapping.interval") or die "can't open file\n";
while (my $interval = ){
chomp($interval);
my @find_interval = split(/\t/, $interval);
my $start = $find_interval[1];
my $end = $find_interval[2];
####
my $switch = 1;
while ($switch == 1) {
open (CG, "/Users/logancurtis-whitchurch/Dropbox/thesis_folder/CompleteGenomics/28_males_inAll/CGS.inall.28.chr.23.txt") or die "can't open CG file\n";
my $output_file = "/Users/logancurtis-whitchurch/Desktop/temp_$count.txt";
open(OUT, ">$output_file");
####
my @SNPs = ;
my @get_SNP = split(/\t/, $SNPs[$placeholder]);
my $position = $get_SNP[3];
####
my $switch = 1;
while ($switch == 1) {
if (($position < $end) && ($position >= $start)) {
print OUT "$SNPs[$placeholder]\n";
$placeholder++;
}
else {
$switch = 0;
$count++;
}
}
}