This file is indexed.

/usr/lib/ChimeraSlayer/ChimeraSlayer.pl is in chimeraslayer 20101212+dfsg-1.

This file is owned by root:root, with mode 0o755.

The actual contents of the file can be viewed below.

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
#!/usr/bin/env perl

use strict;
use warnings;

use Getopt::Long qw(:config no_ignore_case bundling);
#use FindBin;
#Hack needed because FindBin is braindead.
BEGIN{ $FindBin::Bin = "/usr/lib/ChimeraSlayer" };

use lib ("$FindBin::Bin/PerlLib");
use Fasta_reader;
use CdbTools;
use File::Basename;

my $db_NAST = "/usr/share/microbiomeutil-data/RESOURCES/rRNA16S.gold.NAST_ALIGNED.fasta";
my $db_FASTA = "/usr/share/microbiomeutil-data/RESOURCES/rRNA16S.gold.fasta";


my $usage = <<_EOUSAGE_;

##########################################################################################
#
#  Required:
#
#    --query_NAST      multi-fasta file containing query sequences in alignment format
#
#  Common opts:
#
#    --db_NAST        db in NAST format (default: $db_NAST)
#    --db_FASTA       db in fasta format (megablast formatted) (default: $db_FASTA)
#
#
#    -n       number of top matching database sequences to compare to (default 15)
#    -R       min divergence ratio default: 1.007
#    -P       min percent identity among matching sequences (default: 90)
#
#  ## parameters to tune ChimeraParentSelector:
#   
#  Scoring parameters:
#   -M match score   (default: +5)
#   -N mismatch penalty  (default: -4)
#   -Q min query coverage by matching database sequence (default: 70)
#   -T maximum traverses of the multiple alignment  (default: 1)

#
#  ## parameters to tune ChimeraPhyloChecker:
#
#
#    --windowSize                default 50
#    --windowStep                default 5
#    --minBS      minimum bootstrap support for calling chimera (default: 90)
#    --num_BS_replicates         default: 100
#    --low_range_finer_BS (default: 10)    If computed BS is between minBS and (minBS - low_range_finer_BS), then num_finer_BS_replicates computed.
#    --num_finer_BS_replicates (default: 1000)
#    -S           percent of SNPs to sample on each side of breakpoint for computing bootstrap support (default: 10)
#    --num_parents_test       number of potential parents to test for chimeras (default: 3)
#    --MAX_CHIMERA_PARENT_PER_ID    Chimera/Parent alignments with perID above this are considered non-chimeras (default 100; turned off) 
#    
#
#  ## misc opts
#
#   --printFinalAlignments          shows alignment between query sequence and pair of candidate chimera parents
#   --printCSalignments             print ChimeraSlayer alignments in ChimeraSlayer output
#   --exec_dir                      chdir to here before running
#
#########################################################################################

_EOUSAGE_
	
	;

## required options
my ($query_NAST);


# common opts
my $numSeqsCompare = 15;
my $minPerID = 90;
my $minDivR = 1.007;

# chimera parent selector opts:
my ($matchScore, $mismatchPenalty, $maxTraverse, $minQueryCoverage);

# chimera phylochecker opts:
my $minBS = 90;
my ($windowSize, $windowStep, $minPercentSampleSNPs, $num_parents_test, $num_BS_replicates,
	$low_range_finer_BS, $num_finer_BS_replicates,
);

# misc
my $printFinalAlignments = 0;
my $help;
my $exec_dir;
my $printCSalignments;
my $MAX_CHIMERA_PARENT_PER_ID;



&GetOptions ("query_NAST=s" => \$query_NAST,
			 "db_NAST=s" => \$db_NAST,
			 "db_FASTA=s" => \$db_FASTA,
			 
			 # common opts
			 "R=f" => \$minDivR,  # same as -R on chimeraMaligner
			 "n=i" => \$numSeqsCompare,
			 "P=f" => \$minPerID,
			 
			 ## ChimeraParentSelector parameters:
			 "M=i" => \$matchScore,
			 "N=i" => \$mismatchPenalty,
			 "Q=i" => \$minQueryCoverage,
			 "T=i" => \$maxTraverse,

			 ## ChimeraPhyloChecker parameters
			 			 
			 "windowSize=i" => \$windowSize,
			 "windowStep=i" => \$windowStep,
			 "minBS=i" => \$minBS,
			 'num_BS_replicates=i' => \$num_BS_replicates,
			 'low_range_finer_BS=i' => \$low_range_finer_BS,
			 'num_finer_BS_replicates=i' => \$num_finer_BS_replicates,
			 'S=i' => \$minPercentSampleSNPs,
			 "num_parents_test=i" => \$num_parents_test,
			 "MAX_CHIMERA_PARENT_PER_ID=f" => \$MAX_CHIMERA_PARENT_PER_ID,
			 
			 # misc

			 "printFinalAlignments" => \$printFinalAlignments,
			 "exec_dir=s" => \$exec_dir,
			 "printCSalignments" => \$printCSalignments,


			 "h" => \$help,
	);

if ($help) { die $usage; }

unless ($query_NAST) { die $usage; }



main: {
	
	if ($exec_dir) {
		chdir $exec_dir or die "Error, cannot chdir to $exec_dir";
	}
	
	
	## ensure that ncbi blast formatted database exists.
	unless (-e "$db_FASTA.nhr" && -e "$db_FASTA.nin" && -e "$db_FASTA.nsq") {
		my $cmd = "formatdb -i $db_FASTA -p F 2>/dev/null";
		&process_cmd($cmd);
	}
	

	## Run ChimeraParentSelector
	
	my $baseOut = basename($query_NAST);

	my $CPS_outfile = $baseOut . ".CPS";
	
	my $cmd = "$FindBin::Bin/ChimeraParentSelector/run_chimeraParentSelector.pl "
		. " --query_NAST $query_NAST "
		. " --db_NAST $db_NAST "
		. " --db_FASTA $db_FASTA ";

	$cmd .= " -n $numSeqsCompare " if defined($numSeqsCompare);
	$cmd .= " -P $minPerID " if defined($minPerID);
	$cmd .= " -R $minDivR " if defined($minDivR);
	
	$cmd .= " -M $matchScore " if defined($matchScore);
	$cmd .= " -N $mismatchPenalty " if defined($mismatchPenalty);
	$cmd .= " -T $maxTraverse " if defined($maxTraverse);
	$cmd .= " -Q $minQueryCoverage " if defined($minQueryCoverage);
		
	$cmd .= " > $CPS_outfile";
	
	
	&process_cmd($cmd); 
	

	#############################
	####  RE-NAST align based on CPS results

	my $query_RENAST_file = "$query_NAST.CPS_RENAST";
	
	$cmd = "$FindBin::Bin/ChimeraParentSelector/CPS_to_RENAST.pl "
		. "--CPS_output $CPS_outfile "
		. "--query_NAST $query_NAST "
		. "--db_NAST $db_NAST "
		. "> $query_RENAST_file";
	
	&process_cmd($cmd); 
	
	


	########################
	#### Run ChimeraPhyloChecker
	
	my $CPC_outfile = $CPS_outfile . ".CPC";
		
	## Run ChimeraSlayer
	$cmd = "$FindBin::Bin/ChimeraPhyloChecker/CPS_to_CPC.pl "
		. " --CPS_output $CPS_outfile "
		. " --query_NAST $query_RENAST_file "
		. " --db_NAST $db_NAST ";
	
	
	$cmd .= " --num_parents_test $num_parents_test " if defined($num_parents_test);
	$cmd .= " -P $minPerID " if defined($minPerID);
	$cmd .= " -R $minDivR " if defined($minDivR);
	$cmd .= " -S $minPercentSampleSNPs " if defined($minPercentSampleSNPs);
	$cmd .= " --minBS $minBS " if defined ($minBS);
	
	$cmd .= " --num_BS_replicates $num_BS_replicates " if defined($num_BS_replicates);

	$cmd .= " --low_range_finer_BS $low_range_finer_BS " if defined($low_range_finer_BS);
	$cmd .= " --num_finer_BS_replicates $num_finer_BS_replicates " if defined($num_finer_BS_replicates);
	
	
	$cmd .= " --windowSize $windowSize " if defined($windowSize);
	$cmd .= " --windowStep $windowStep " if defined($windowStep);

	$cmd .= " --printAlignments " if ($printCSalignments);
	
	$cmd .= " > $CPC_outfile ";
	
	&process_cmd($cmd);
	
	## Add species and chimera type information
	$cmd = "$FindBin::Bin/util/CS_add_taxonomy.pl < $CPC_outfile > $CPC_outfile.wTaxons";
	&process_cmd($cmd);
	

	if ($printFinalAlignments) {
		## Run ChimeraParentSelector again to print alignments and reestimate breakpoints
		my $outfile = "$CPC_outfile.align";
		
		$cmd = "$FindBin::Bin/ChimeraParentSelector/CPC_to_CPS.pl "
			. " --CPC_output $CPC_outfile "
			. " --query_NAST $query_RENAST_file "
			. " --db_NAST $db_NAST "
			. " > $outfile ";
		
		&process_cmd($cmd);
	}
	
	exit(0);
}


####
sub process_cmd {
	my ($cmd) = @_;
	
	print "CMD: $cmd\n";
	my $ret = system($cmd);

	if ($ret) {
		die "Error, cmd ($cmd) died with ret($ret)";
	}

	return;
}