TCGA_ATAC-seq 比对文件整理 （Adjusting alignment file of ATAC-seq in TCGA）

#!/usr/bin/perl -w
use strict;
use warnings;
use Cwd;
use File::Copy;

my $usage=<<USAGE;
Usage:
    perl $0 filename
    将该目录下的所有文件放在一个叫filename的文件夹下
USAGE
if(@ARGV==0){die $usage};

my $file1=$ARGV[0];
my $dir=getcwd;
my @filename;

opendir(DIR,$dir) || die $!;
my @temp=readdir(DIR);
for my $i (0..$#temp){
    if ((-d "$dir/$temp[$i]") & ($temp[$i] ne $file1)){
        opendir(DIR2,"$dir/$temp[$i]") || die $!;
        my @temp1=readdir(DIR2);
        for my $j (0..$#temp1){
            if((-f "$dir/$temp[$i]/$temp1[$j]") & (($temp1[$j])=~/.bam$/) | ($temp1[$j]=~/.bai$/) ){
                push(@filename,"$dir/$temp[$i]/$temp1[$j]");
            }
        }
    }
}
close(DIR2);
close(DIR);

#print @filename;

my $lujing="$dir/$file1";

if(!(-e $lujing)){
    mkdir ($lujing) || die "can not create $file1 doctory";
}

for my $i (0..$#filename){
    move("$filename[$i]","$lujing");
}