mutations_caller_pipeline_aws 0.0.7 → 0.0.8
Sign up to get free protection for your applications and to get access to all the features.
@@ -4,7 +4,7 @@ class GatkCaller
|
|
4
4
|
def self.call(log_dir, gatk, index_fa, read_bam, read_vcf, job_prefix, account, debug)
|
5
5
|
cmd = "echo 'starting GATK for mutant at ' `date` >> #{log_dir}
|
6
6
|
qsub -cwd -b y -N #{job_prefix}_genotyper -l h_vmem=3G -hold_jid #{job_prefix}_recalibration #{account}\
|
7
|
-
|
7
|
+
#{gatk} -l INFO -R #{index_fa} -T UnifiedGenotyper \
|
8
8
|
-I #{read_bam} \
|
9
9
|
-o #{read_vcf} \
|
10
10
|
--genotype_likelihoods_model BOTH \
|
@@ -16,8 +16,8 @@ class GatkCaller
|
|
16
16
|
# Making recalibration table
|
17
17
|
def self.recalibrate_bam(log_dir ,gatk, index_fa, index_vcf, read_bam, recal_file, job_prefix, account, debug )
|
18
18
|
cmd = "echo 'starting recalibration table ' `date` >> #{log_dir}
|
19
|
-
qsub -cwd -b y -N #{job_prefix}_recalibration_table -l h_vmem=
|
20
|
-
|
19
|
+
qsub -cwd -b y -N #{job_prefix}_recalibration_table -l h_vmem=4G #{account} \
|
20
|
+
#{gatk} -knownSites #{index_vcf} -I #{read_bam} \
|
21
21
|
-R #{index_fa} -T CountCovariates \
|
22
22
|
-cov ReadGroupCovariate -cov QualityScoreCovariate -cov DinucCovariate \
|
23
23
|
-cov CycleCovariate \
|
@@ -29,8 +29,8 @@ class GatkCaller
|
|
29
29
|
# Using recalibration table
|
30
30
|
def self.table_calibration(log_dir, gatk, index_fa, read_bam, recal_bam, recal_file, job_prefix, account, debug)
|
31
31
|
cmd = "echo 'recalibrating bam_file at ' `date` >> #{log_dir}
|
32
|
-
qsub -cwd -b y -N #{job_prefix}_recalibration -l h_vmem=
|
33
|
-
|
32
|
+
qsub -cwd -b y -N #{job_prefix}_recalibration -l h_vmem=4G -hold_jid #{job_prefix}_recalibration_table #{account} \
|
33
|
+
#{gatk} \
|
34
34
|
-R #{index_fa} \
|
35
35
|
-I #{read_bam} \
|
36
36
|
-T TableRecalibration \
|
@@ -43,8 +43,8 @@ class GatkCaller
|
|
43
43
|
# Preparation realignement
|
44
44
|
def self.prepare_realigne(log_dir, gatk, read_bam, index_fa, target_intervals, job_prefix, account, debug)
|
45
45
|
cmd = "echo 'preparing realignement at ' `date` >> #{log_dir}
|
46
|
-
qsub -cwd -b y -N #{job_prefix}_prep_realignment -l h_vmem=
|
47
|
-
|
46
|
+
qsub -cwd -b y -N #{job_prefix}_prep_realignment -l h_vmem=4G -hold_jid #{job_prefix}_indexing #{account}\
|
47
|
+
#{gatk} \
|
48
48
|
-I #{read_bam} \
|
49
49
|
-R #{index_fa} \
|
50
50
|
-T RealignerTargetCreator \
|
@@ -56,8 +56,8 @@ class GatkCaller
|
|
56
56
|
# Realignment
|
57
57
|
def self.realigne(log_dir, gatk, read_bam, index_fa, target_intervals, realigned_bam, job_prefix, account, debug)
|
58
58
|
cmd = "echo 'preparing realignement at ' `date` >> #{log_dir}
|
59
|
-
qsub -cwd -b y -N #{job_prefix}_realignment -l h_vmem=
|
60
|
-
|
59
|
+
qsub -cwd -b y -N #{job_prefix}_realignment -l h_vmem=4G -hold_jid #{job_prefix}_prep_realignment #{account} \
|
60
|
+
#{gatk} \
|
61
61
|
-I #{read_bam} \
|
62
62
|
-R #{index_fa} \
|
63
63
|
-T IndelRealigner \
|