lib/seqtrimnext/plugins/plugin_contaminants.rb in seqtrimnext-2.0.45 vs lib/seqtrimnext/plugins/plugin_contaminants.rb in seqtrimnext-2.0.46

- old
+ new

@@ -34,11 +34,11 @@ # TODO - Culling limit = 2 porque el blast falla con este comando cuando se le pasa cl=1 y dust=no # y una secuencia de baja complejidad como entrada blast = BatchBlast.new("-db #{@params.get_param('contaminants_db')}",'blastn'," -task blastn -evalue #{@params.get_param('blast_evalue_contaminants')} -perc_identity #{@params.get_param('blast_percent_contaminants')} -culling_limit 1") #get contaminants -max_target_seqs #{MAX_TARGETS_SEQS} - $LOG.info('BLAST:'+blast.get_blast_cmd(:xml)) + $LOG.debug('BLAST:'+blast.get_blast_cmd(:xml)) fastas=[] seqs.each do |seq| fastas.push ">"+seq.seq_name @@ -65,11 +65,11 @@ def exec_seq(seq,blast_query) if blast_query.query_id != seq.seq_name # raise "Blast and seq names does not match, blast:#{blast_query.query_id} sn:#{seq.seq_name}" end - $LOG.info "[#{self.class.to_s}, seq: #{seq.seq_name}]: looking for contaminants into the sequence" + $LOG.debug "[#{self.class.to_s}, seq: #{seq.seq_name}]: looking for contaminants into the sequence" #blast = BatchBlast.new('-db DB/formatted/contaminants.fasta','blastn',' -task blastn -evalue 1e-10 -perc_identity 95') #get contaminants # blast = BatchBlast.new("-db #{@params.get_param('contaminants_db')}",'blastn'," -task blastn-short -evalue #{@params.get_param('blast_evalue_contaminants')} -perc_identity #{@params.get_param('blast_percent_contaminants')} -culling_limit 1") #get contaminants -max_target_seqs #{MAX_TARGETS_SEQS} @@ -141,10 +141,10 @@ add_stats('contaminants_ids',c) end end else - $LOG.info('Contaminant ignored due to genus match: '+c.definition) + $LOG.debug('Contaminant ignored due to genus match: '+c.definition) end end reject=@params.get_param('contaminants_reject') # cond_if=false