59 |
mkdir("$dir/$ssN/$roleN",0777) || die "could not make $dir/$ssN/$roleN"; |
mkdir("$dir/$ssN/$roleN",0777) || die "could not make $dir/$ssN/$roleN"; |
60 |
&FIG::run("split_and_trim_sequences $dir/$ssN/$roleN/split_info < $tmpF"); |
&FIG::run("split_and_trim_sequences $dir/$ssN/$roleN/split_info < $tmpF"); |
61 |
|
|
62 |
if (-s "$dir/$ssN/$roleN/split_info") |
if (-s "$dir/$ssN/$roleN/split_info/set.sizes") |
63 |
{ |
{ |
64 |
open(SZ,"<$dir/$ssN/$roleN/split_info/set.sizes") || die " could not open $dir/$ssN/$roleN/split_info"; |
open(SZ,"<$dir/$ssN/$roleN/split_info/set.sizes") || die " could not open $dir/$ssN/$roleN/split_info/set.sizes"; |
65 |
while (defined($_ = <SZ>)) |
while (defined($_ = <SZ>)) |
66 |
{ |
{ |
67 |
if (($_ =~ /^(\d+)\t(\d+)/) && ($2 > 3)) |
if (($_ =~ /^(\d+)\t(\d+)/) && ($2 > 3)) |