#!/usr/bin/env cwl-runner cwlVersion: v1.0 class: Workflow label: "exome alignment and variant detection" requirements: - class: SubworkflowFeatureRequirement inputs: reference: string bams: type: File[] readgroups: type: string[] mills: type: File secondaryFiles: [.tbi] known_indels: type: File secondaryFiles: [.tbi] dbsnp_vcf: type: File secondaryFiles: [.tbi] bqsr_intervals: type: string[]? bait_intervals: type: File target_intervals: type: File per_target_intervals: type: File per_target_bait_intervals: type: File per_base_intervals: type: File per_base_bait_intervals: type: File omni_vcf: type: File secondaryFiles: [.tbi] picard_metric_accumulation_level: type: string varscan_strand_filter: type: int? default: 0 varscan_min_coverage: type: int? default: 8 varscan_min_var_freq: type: float? default: 0.1 varscan_p_value: type: float? default: 0.99 varscan_min_reads: type: int? default: 2 maximum_population_allele_frequency: type: float? default: 0.001 vep_cache_dir: type: string? synonyms_file: type: File? annotate_coding_only: type: boolean? default: true hgvs_annotation: type: boolean? default: true variants_to_table_fields: type: string[]? default: [CHROM,POS,REF,ALT,set] variants_to_table_genotype_fields: type: string[]? default: [GT,AD,AF,DP] vep_to_table_fields: type: string[]? default: [Consequence,SYMBOL,Feature_type,Feature,HGVSc,HGVSp,cDNA_position,CDS_position,Protein_position,Amino_acids,Codons,HGNC_ID,Existing_variation,gnomADe_AF,CLIN_SIG,SOMATIC,PHENO] sample_name: type: string docm_vcf: type: File secondaryFiles: [.tbi] custom_gnomad_vcf: type: File? secondaryFiles: [.tbi] qc_minimum_mapping_quality: type: int? qc_minimum_base_quality: type: int? readcount_minimum_mapping_quality: type: int? readcount_minimum_base_quality: type: int? outputs: cram: type: File outputSource: alignment_and_qc/cram mark_duplicates_metrics: type: File outputSource: alignment_and_qc/mark_duplicates_metrics insert_size_metrics: type: File outputSource: alignment_and_qc/insert_size_metrics insert_size_histogram: type: File outputSource: alignment_and_qc/insert_size_histogram alignment_summary_metrics: type: File outputSource: alignment_and_qc/alignment_summary_metrics hs_metrics: type: File outputSource: alignment_and_qc/hs_metrics per_target_coverage_metrics: type: File? outputSource: alignment_and_qc/per_target_coverage_metrics per_target_hs_metrics: type: File? outputSource: alignment_and_qc/per_target_hs_metrics per_base_coverage_metrics: type: File? outputSource: alignment_and_qc/per_base_coverage_metrics per_base_hs_metrics: type: File? outputSource: alignment_and_qc/per_base_hs_metrics flagstats: type: File outputSource: alignment_and_qc/flagstats verify_bam_id_metrics: type: File outputSource: alignment_and_qc/verify_bam_id_metrics verify_bam_id_depth: type: File outputSource: alignment_and_qc/verify_bam_id_depth varscan_vcf: type: File outputSource: detect_variants/varscan_vcf secondaryFiles: [.tbi] docm_gatk_vcf: type: File outputSource: detect_variants/docm_gatk_vcf annotated_vcf: type: File outputSource: detect_variants/annotated_vcf secondaryFiles: [.tbi] final_vcf: type: File outputSource: detect_variants/final_vcf secondaryFiles: [.tbi] final_tsv: type: File outputSource: detect_variants/final_tsv vep_summary: type: File outputSource: detect_variants/vep_summary tumor_bam_readcount_tsv: type: File outputSource: detect_variants/tumor_bam_readcount_tsv steps: alignment_and_qc: run: exome_alignment.cwl in: reference: reference bams: bams readgroups: readgroups mills: mills known_indels: known_indels dbsnp_vcf: dbsnp_vcf bqsr_intervals: bqsr_intervals bait_intervals: bait_intervals target_intervals: target_intervals per_target_intervals: per_target_intervals per_target_bait_intervals: per_target_bait_intervals per_base_intervals: per_base_intervals per_base_bait_intervals: per_base_bait_intervals omni_vcf: omni_vcf picard_metric_accumulation_level: picard_metric_accumulation_level qc_minimum_mapping_quality: qc_minimum_mapping_quality qc_minimum_base_quality: qc_minimum_base_quality out: [cram, mark_duplicates_metrics, insert_size_metrics, insert_size_histogram, alignment_summary_metrics, hs_metrics, per_target_coverage_metrics, per_target_hs_metrics, per_base_coverage_metrics, per_base_hs_metrics, flagstats, verify_bam_id_metrics, verify_bam_id_depth] detect_variants: run: detect_variants/tumor_only_detect_variants.cwl in: reference: reference cram: alignment_and_qc/cram interval_list: target_intervals varscan_strand_filter: varscan_strand_filter varscan_min_coverage: varscan_min_coverage varscan_min_var_freq: varscan_min_var_freq varscan_p_value: varscan_p_value varscan_min_reads: varscan_min_reads maximum_population_allele_frequency: maximum_population_allele_frequency vep_cache_dir: vep_cache_dir synonyms_file: synonyms_file variants_to_table_fields: variants_to_table_fields variants_to_table_genotype_fields: variants_to_table_genotype_fields vep_to_table_fields: vep_to_table_fields sample_name: sample_name docm_vcf: docm_vcf hgvs_annotation: hgvs_annotation custom_gnomad_vcf: custom_gnomad_vcf readcount_minimum_mapping_quality: readcount_minimum_mapping_quality readcount_minimum_base_quality: readcount_minimum_base_quality out: [varscan_vcf, docm_gatk_vcf, annotated_vcf, final_vcf, final_tsv, vep_summary, tumor_bam_readcount_tsv]