Fix a bug of clremovecontam and refine clmakeidentdb

astanabe · Apr 7, 2022 · 5fed601 · 5fed601
1 parent 271500b
commit 5fed601
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/clmakeidentdb.pl b/clmakeidentdb.pl
@@ -1,6 +1,5 @@
 use strict;
 use File::Spec;
-use Fcntl ':flock';
 use DBI;
 
 my $buildno = '0.9.x';

diff --git a/clremovecontam.pl b/clremovecontam.pl
@@ -980,7 +980,12 @@ sub estimateContaminationProbability {
 
 sub calculateTagJumpProbability {
 	print(STDERR "Calculating tag jump probability...\n");
-	my @samplenames = keys(%sample2blank);
+	my @samplenames;
+	foreach my $samplename (sort(keys(%samplenames))) {
+		if (exists($sample2blank{$samplename}) || exists($sample2sample{$samplename})) {
+			push(@samplenames, $samplename);
+		}
+	}
 	my @otunames;
 	foreach my $otuname (keys(%otunames)) {
 		if (!defined($ignoreotulist{$otuname})) {