bug fixed: get_homologues-pl -X now removes previous DIAMOND results …

…of one vs others when new genomes are added
eead-csic-compbio · Aug 28, 2018 · 9b2ee98 · 9b2ee98
1 parent 17822b1
commit 9b2ee98
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 3 deletions.
diff --git a/CHANGES.txt b/CHANGES.txt
@@ -264,3 +264,4 @@
 06062018: fixed POCP computation; now it is [100(C1 + C2)/(T1 + T2)]
 06062018: fixed POCS computation; now it is [100(C1 + C2)/(T1 + T2)] and T1/T2 are #nr seqs
 23072018: updated Grid Engine instructions in manuals
+28082018: bug fixed: get_homologues-pl -X now removes previous DIAMOND results of one vs others when new genomes are added (thanks alexweisberg)
diff --git a/get_homologues.pl b/get_homologues.pl
@@ -1441,10 +1441,18 @@
         if(!-s $blast_file){ push(@tmp_blast_output_files,$blastout); }
         next;
       }
-      elsif($COMPRESSBLAST && -s $blastout.'.gz' && $current_files eq $previous_files) 
+      elsif($COMPRESSBLAST && -s $blastout.'.gz') 
       {
-        if(!-s $blast_file){ push(@tmp_blast_output_files,$blastout); }
-        next;
+        if($current_files eq $previous_files) # check previous runs
+        {
+          if(!-s $blast_file){ push(@tmp_blast_output_files,$blastout); }
+          next;
+        }
+        else
+        {
+          print "# delete previous output $blastout.gz\n";
+          unlink($blastout.'.gz');
+        }  
       }
 
       $command = format_DIAMONDblastp_command("$newDIR/$new_infile",