mutations_caller_pipeline_aws 0.0.7 → 0.0.8
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
|
@@ -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 \
|