mutations_caller_pipeline_aws 0.0.6 → 0.0.7

Sign up to get free protection for your applications and to get access to all the features.
@@ -44,7 +44,7 @@ options = { :bam_file_sorted => nil,
44
44
  :vcf => nil,
45
45
  :account => "",
46
46
  :project => "",
47
- :debug = 1,
47
+ :debug => 1,
48
48
  }
49
49
 
50
50
  optparse = OptionParser.new do |opts|
@@ -113,7 +113,7 @@ options[:account] = options[:project] if options[:account].empty?
113
113
  # Indexing
114
114
  SamtoolsIndexing.call(bam_file,
115
115
  job_prefix,
116
- account,
116
+ options[:account],
117
117
  options[:debug])
118
118
 
119
119
  # Realigne
@@ -10,7 +10,7 @@ class BwaCaller
10
10
  def self.call_paired_end(r1, r2, out_file, index, log_file, bwa, samtools, job_prefix,account, debug)
11
11
  cmd = "qsub -cwd -b y -N #{job_prefix}_bwa -l h_vmem=9G -pe make 4 #{account}\
12
12
  #{bwa} sampe -r '@RG\tID:foo\tSM:bar\tPL:Illumina' #{index} \
13
- <(#{bwa} aln #{index} #{r1} 2>>#{log_file} || exit 1) <(#{bwa} aln #{index} #{r2} 2>>#{log_file} ) \
13
+ <(#{bwa} aln #{index} #{r1} 2>>#{log_file}) <(#{bwa} aln #{index} #{r2} 2>>#{log_file} ) \
14
14
  #{r1} #{r2} 2>>#{log_file} | #{samtools} view -Su - 2>>#{log_file} | #{samtools} sort - #{out_file} 2>>#{log_file}"
15
15
  puts cmd
16
16
  system('bash','-c', cmd) if debug == 1
metadata CHANGED
@@ -1,7 +1,7 @@
1
1
  --- !ruby/object:Gem::Specification
2
2
  name: mutations_caller_pipeline_aws
3
3
  version: !ruby/object:Gem::Version
4
- version: 0.0.6
4
+ version: 0.0.7
5
5
  prerelease:
6
6
  platform: ruby
7
7
  authors: