Class: CagnutBwa::Aln
- Inherits:
-
Object
- Object
- CagnutBwa::Aln
- Extended by:
- Forwardable
- Defined in:
- lib/cagnut_bwa/functions/aln.rb
Instance Method Summary collapse
- #aln_params_for_r1 ⇒ Object
- #aln_params_for_r2 ⇒ Object
- #fetch_filename ⇒ Object
- #generate_script ⇒ Object
-
#initialize(opts = {}) ⇒ Aln
constructor
A new instance of Aln.
- #job_params(script_name) ⇒ Object
- #queuing_options(previous_job_id = nil) ⇒ Object
- #run(previous_job_id = nil) ⇒ Object
Constructor Details
#initialize(opts = {}) ⇒ Aln
Returns a new instance of Aln.
8 9 10 11 12 13 14 15 16 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 8 def initialize opts = {} @order = sprintf '%02i', opts[:order] @input = opts[:input].nil? ? "#{seqs_path}" : opts[:input] @input2 = File. fetch_filename, File.dirname(@input) abort('Cant recognized sequence files') if @input2.nil? @output = "#{opts[:dirs][:output]}/#{File.basename(@input)}.sai" @output2 = "#{opts[:dirs][:output]}/#{fetch_filename}.sai" @job_name = "#{prefix_name}_#{sample_name}_Aln" end |
Instance Method Details
#aln_params_for_r1 ⇒ Object
45 46 47 48 49 50 51 52 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 45 def aln_params_for_r1 array = aln_params.dup array.insert 1, 'aln' array << "#{ref_fasta}" array << "#{@input}" array << "> #{@output}" array.uniq end |
#aln_params_for_r2 ⇒ Object
54 55 56 57 58 59 60 61 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 54 def aln_params_for_r2 array = aln_params.dup array.insert 1, 'aln' array << "#{ref_fasta}" array << "#{@input2} >" array << "#{@output2}" array.uniq end |
#fetch_filename ⇒ Object
18 19 20 21 22 23 24 25 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 18 def fetch_filename filename = File.basename(@input) if filename.match '_R1_' filename.gsub '_R1_', '_R2_' elsif filename.match '_1_' filename.gsub '_1_', '_2_' end end |
#generate_script ⇒ Object
63 64 65 66 67 68 69 70 71 72 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 63 def generate_script script_name = "#{@order}_bwa_aln" file = File.join jobs_dir, "#{script_name}.sh" template = Tilt.new(File. '../templates/aln.sh', __FILE__) File.open(file, 'w') do |f| f.puts template.render Object.new, job_params(script_name) end File.chmod(0700, file) script_name end |
#job_params(script_name) ⇒ Object
74 75 76 77 78 79 80 81 82 83 84 85 86 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 74 def job_params script_name { jobs_dir: jobs_dir, script_name: script_name, input: @input, input2: @input2, output: @output, output2: @output2, aln_params_for_r1: aln_params_for_r1, aln_params_for_r2: aln_params_for_r2, run_local: ::Cagnut::JobManage.run_local } end |
#queuing_options(previous_job_id = nil) ⇒ Object
34 35 36 37 38 39 40 41 42 43 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 34 def previous_job_id = nil threads = 2 { previous_job_id: previous_job_id, var_env: [ref_fasta], adjust_memory: ["h_vmem=adjustWorkingMem 5G #{threads}"], parallel_env: [threads], tools: ['bwa', 'aln'] } end |
#run(previous_job_id = nil) ⇒ Object
27 28 29 30 31 32 |
# File 'lib/cagnut_bwa/functions/aln.rb', line 27 def run previous_job_id = nil puts "Submitting bwaAln #{sample_name}" script_name = generate_script ::Cagnut::JobManage.submit script_name, @job_name, (previous_job_id) @job_name end |