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.
data/bin/gatk_pipe_only_aws
CHANGED
@@ -44,7 +44,7 @@ options = { :bam_file_sorted => nil,
|
|
44
44
|
:vcf => nil,
|
45
45
|
:account => "",
|
46
46
|
:project => "",
|
47
|
-
:debug
|
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}
|
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
|