Skip to content

Commit

Permalink
Fixing human gatk4 var calling
Browse files Browse the repository at this point in the history
  • Loading branch information
mwhamgenomics committed Jul 1, 2019
1 parent a33c522 commit 6920bc1
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 3 deletions.
4 changes: 2 additions & 2 deletions analysis_driver/pipelines/human_variant_calling_gatk4.py
Original file line number Diff line number Diff line change
Expand Up @@ -198,11 +198,11 @@ def stage(cls, **params):
output_vcf_file=hard_filter_indels.hard_filtered_vcf,
previous_stages=[hard_filter_snps, hard_filter_indels])

gender_val = stage(qc.GenderValidation, vcf_file=merge_hf.hard_filtered_vcf, previous_stages=[merge_hf])
sex_val = stage(qc.SexValidation, vcf_file=merge_hf.hard_filtered_vcf, previous_stages=[merge_hf])

vcfstats = stage(qc.VCFStats, vcf_file=merge_hf.hard_filtered_vcf, previous_stages=[merge_hf])

final_stages = [contam, blast, geno_val, gender_val, vcfstats, verify_bam_id, samtools_depth, samtools_stat,
final_stages = [contam, blast, geno_val, sex_val, vcfstats, verify_bam_id, samtools_depth, samtools_stat,
gather_gcvf]

output = stage(common.SampleDataOutput, previous_stages=final_stages, output_fileset='gatk4_human_var_calling')
Expand Down
2 changes: 1 addition & 1 deletion etc/output_files.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -291,7 +291,7 @@ gatk4_human_var_calling:
basename: '{sample_id}-chr22-vbi.selfSM'
new_name: '{user_sample_id}-chr22-vbi.selfSM'

gender_call:
sex_validation:
location: ['gatk4']
basename: '{user_sample_id}_hard_filter.sex'
new_name: '{user_sample_id}.sex'
Expand Down

0 comments on commit 6920bc1

Please sign in to comment.