- 论坛徽章:
- 0
|
详细如下:
#!/usr/bin/perl -w
use strict;
use File::Basename qw(basename dirname);
my $input_dir = "path/go";
my $max_len_dir="path/back";
my (@max_len,@len_out);
my $file=shift;
my $out5;
my $out6;
my ($out1,$out2,$next_file,$out3,$out4,$id,$wole,$avg_ins,$input_file_basename);
my ($out_single1,$out_single2,$out_pair1,$out_pair2);
open IN,$file || die "fail";
while(<IN>)
{
chomp;
split;
my @dir=split /\s+/,$_;
$id=$dir[0];
$avg_ins=$dir[1];
opendir DH,$input_dir;
foreach my $fiile(readdir DH){
if($fiile=~/^$id/){
$next_file=$fiile;
}
}
closedir DH;
my $whole="$input_dir/$next_file";
opendir DHI,$whole;
foreach $wole(readdir DHI){
if ($wole=~/single$/){
$out_single1="$input_dir/$next_file/$wole";
$out_single2="f=$out_single1\n";
my $whoe="$input_dir/$next_file/$wole";
`perl $max_len_dir/max_len.pl $whoe >$max_len_dir/lenn`;
open my $fh,"<","$max_len_dir/lenn" ||die "fail";
while( <$fh>){
my $kin=$_;
push @len_out, $kin;
}
close $fh; |
|