113 |
foreach my $i (7..12) |
foreach my $i (7..12) |
114 |
{ |
{ |
115 |
mkdir("$newD/Merged/$i",0777) || die "could not make $newD/Merged/$i: $!"; |
mkdir("$newD/Merged/$i",0777) || die "could not make $newD/Merged/$i: $!"; |
116 |
my $h = start(["FFB2_merge_oligos $newD/Kmers/$i/good.oligos $newD/PhyloSigs/$i/good.oligos"], "|", |
my $h = start(["FFB2_merge_oligos", "$newD/Kmers/$i/good.oligos", "$newD/PhyloSigs/$i/good.oligos"], "|", |
117 |
["FFB2_create_binary_kmers", "-l", $i, "-s", "4,2", "-", "$newD/Merged/$i/table.binary"]); |
["FFB2_create_binary_kmers", "-l", $i, "-s", "4,2", "-", "$newD/Merged/$i/table.binary"]); |
118 |
|
|
119 |
push(@merge_handles, [$h, $i]); |
push(@merge_handles, [$h, $i]); |
131 |
|
|
132 |
my $load1 = start(["FFB2_load_oligo_index", "$newD/function.index", "$newD/FRI.db"]); |
my $load1 = start(["FFB2_load_oligo_index", "$newD/function.index", "$newD/FRI.db"]); |
133 |
my $load2 = start(["FFB2_load_oligo_index", "$newD/setI", "$newD/setI.db"]); |
my $load2 = start(["FFB2_load_oligo_index", "$newD/setI", "$newD/setI.db"]); |
134 |
my $build = start(["build_prok_nonff_fasta $newD > $newD/extra_prok_seqs.fasta"]); |
my $build = start(["build_prok_nonff_fasta", "$newD"], ">", "$newD/extra_prok_seqs.fasta"); |
135 |
|
|
136 |
|
|
137 |
#foreach my $i (7..12) |
#foreach my $i (7..12) |