@prefix LoadListingRequirement: . @prefix SchemaDefRequirement: . @prefix CommandLineBinding: . @prefix xsd: . @prefix DockerRequirement: . @prefix Workflow: . @prefix sld: . @prefix rdfs: . @prefix CommandOutputBinding: . @prefix ns1: . @prefix NetworkAccess: . @prefix rdf: . @prefix ScatterMethod: . @prefix cwl: . @prefix ResourceRequirement: . cwl:inputBinding [ CommandLineBinding:prefix "-1" ] ; sld:type cwl:File , sld:null . cwl:source . sld:type [ sld:items ; sld:type sld:array ] . cwl:source . a cwl:Workflow ; rdfs:label "wgs alignment with qc" ; ns1:original_cwlVersion "v1.0" ; Workflow:steps , ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs , , , , , , , , , , , , , , ; cwl:outputs , , , , , , , , , , , , , , , , ; cwl:requirements [ a cwl:SubworkflowFeatureRequirement ] ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types , ] . cwl:source . cwl:inputBinding [ CommandLineBinding:position 2 ; CommandLineBinding:prefix "-I" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:in , , ; cwl:out , ; cwl:run . cwl:inputBinding [ CommandLineBinding:position 1 ] ; sld:type cwl:File . cwl:source . cwl:in , , , , , , , , , ; cwl:out ; cwl:run ; cwl:scatter ( ) ; cwl:scatterMethod ScatterMethod:dotproduct . cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . sld:type xsd:int , sld:null . sld:type [ sld:items ; sld:type sld:array ] . sld:type [ sld:items ; sld:type sld:array ] . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; sld:type xsd:string . cwl:source . a cwl:CommandLineTool ; rdfs:label "create BQSR table" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "-o" [ cwl:valueFrom "$(runtime.outdir)/bqsr.table" ] "--preserve_qscores_less_than" "6" "--disable_auto_index_creation_and_locking_when_reading_rods" "--disable_bam_indexing" "-dfrac" ".1" "-nct" "4" ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/opt/GenomeAnalysisTK.jar" "-T" "BaseRecalibrator" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , ; cwl:outputs ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/gatk-cwl:3.6.0" ] . cwl:inputBinding [ CommandLineBinding:prefix "R=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:source ; cwl:valueFrom "$(self.sequence.hasOwnProperty('fastq2')? self.sequence.fastq2 : null)" . cwl:source . cwl:outputSource ; sld:type cwl:File . cwl:in , ; cwl:out ; cwl:run . cwl:source . cwl:in , , , , ; cwl:out ; cwl:run . cwl:inputBinding [ CommandLineBinding:prefix "--bam" ] ; sld:type cwl:File . cwl:source . cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:outputSource ; sld:type cwl:File . cwl:default "out" ; sld:type xsd:string , sld:null . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.basename).flagstat" ] ; sld:type cwl:File . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . sld:type [ sld:items ; sld:type sld:array ] . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).VerifyBamId.depthSM" ] ; sld:type cwl:File . cwl:source , , . cwl:inputBinding [ CommandLineBinding:prefix "R=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "R=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "R=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type xsd:string , cwl:File . cwl:source . cwl:in , , , ; cwl:out ; cwl:run . cwl:source . rdfs:comment "the unaligned sequence data with readgroup information" ; sld:type . cwl:source ; cwl:valueFrom "$(self.file)" . cwl:outputSource ; sld:type cwl:File . cwl:source . sld:type [ sld:items ; sld:type sld:array ] . cwl:in ; cwl:out , ; cwl:run . cwl:source . cwl:source ; cwl:valueFrom "$(self.file)" . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).mark_dups_metrics.txt" ] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "PER_TARGET_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerTargetCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; cwl:inputBinding [ CommandLineBinding:prefix "PER_TARGET_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerTargetCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; cwl:inputBinding [ CommandLineBinding:prefix "PER_TARGET_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerTargetCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; sld:type xsd:boolean , sld:null . cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; sld:type xsd:string . cwl:outputSource ; sld:type cwl:File . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.output_name).bam" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:source . a cwl:CommandLineTool ; rdfs:label "collect HS metrics" ; ns1:original_cwlVersion "v1.2" ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ) ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ) ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx48g" "-jar" "/usr/picard/picard.jar" "CollectHsMetrics" ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx48g" "-jar" "/usr/picard/picard.jar" "CollectHsMetrics" ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx48g" "-jar" "/usr/picard/picard.jar" "CollectHsMetrics" ) ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , , , , , , , ; cwl:outputs , , ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] ; cwl:requirements [ a cwl:InlineJavascriptRequirement ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] ; cwl:requirements [ a cwl:InlineJavascriptRequirement ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 52000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 52000 ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 52000 ] ; cwl:requirements [ a cwl:InlineJavascriptRequirement ] . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).VerifyBamId.selfSM" ] ; sld:type cwl:File . cwl:source ; cwl:valueFrom "base-$(self.label)" . cwl:inputBinding [ CommandLineBinding:prefix "INPUT=" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:secondaryFiles [] ; sld:type cwl:File . sld:type cwl:File . sld:type xsd:int , sld:null . a cwl:Workflow ; rdfs:label "WGS QC workflow" ; ns1:original_cwlVersion "v1.0" ; Workflow:steps , , , , , , ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , , , , , , , , ; cwl:outputs , , , , , , , , , , , , , , ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types ] ; cwl:requirements [ a cwl:SubworkflowFeatureRequirement ] . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:source . sld:type xsd:string . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "--vcf" ] ; sld:type cwl:File . cwl:default "final" ; sld:type xsd:string , sld:null . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).AlignmentSummaryMetrics.txt" ] ; sld:type cwl:File . cwl:source ; cwl:valueFrom "$(self.file)" . cwl:source . sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:position 3 ; CommandLineBinding:prefix "-I" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.sample_name).GcBiasMetricsChart.pdf" ] ; sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . cwl:secondaryFiles [] ; sld:type cwl:File . a cwl:Workflow ; rdfs:comment "Some workflow engines won't stage files in our nested structure, so parse it out here" ; rdfs:label "adapter for sequence_align_and_tag" ; ns1:original_cwlVersion "v1.0" ; Workflow:steps ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs , ; cwl:outputs ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types ] ; cwl:requirements [ a cwl:InlineJavascriptRequirement ] ; cwl:requirements [ a cwl:StepInputExpressionRequirement ] ; cwl:requirements [ a cwl:SubworkflowFeatureRequirement ] . cwl:source ; cwl:valueFrom "$(self.file)" . sld:type xsd:int , sld:null . cwl:source . a cwl:CommandLineTool ; rdfs:label "collect gc bias metrics" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.sample_name).GcBiasMetrics.txt" ] "CHART=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.sample_name).GcBiasMetricsChart.pdf" ] "S=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.sample_name).GcBiasMetricsSummary.txt" ] ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/usr/picard/picard.jar" "CollectGcBiasMetrics" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , ; cwl:outputs , , ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/cle:v1.3.1" ] . cwl:source . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:secondaryFiles [] ; sld:type cwl:File . rdfs:comment "bwa-indexed reference file" ; cwl:inputBinding [ CommandLineBinding:position 4 ; CommandLineBinding:prefix "-r" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . sld:type cwl:File . cwl:source ; cwl:valueFrom "$(self.readgroup)" . cwl:default "final" ; sld:type xsd:string , sld:null . cwl:default "queryname" ; cwl:inputBinding [ CommandLineBinding:position 5 ] ; sld:type xsd:string . sld:type , . cwl:source . sld:type xsd:int , sld:null . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:outputBinding [ CommandOutputBinding:glob "bqsr.table" ] ; sld:type cwl:File . a cwl:CommandLineTool ; rdfs:label "collect alignment summary metrics" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "OUTPUT=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).AlignmentSummaryMetrics.txt" ] ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/usr/picard/picard.jar" "CollectAlignmentSummaryMetrics" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , ; cwl:outputs ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:source . sld:type cwl:File . cwl:in , , ; cwl:out ; cwl:run . cwl:outputSource ; sld:type cwl:File . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:in , ; cwl:out , ; cwl:run . a cwl:CommandLineTool ; rdfs:label "Samtools: merge" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "$(inputs.name).merged.bam" [ CommandLineBinding:prefix "--threads" ; cwl:valueFrom "$(runtime.cores)" ] ) ; cwl:baseCommand ( "/opt/samtools/bin/samtools" "merge" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , ; cwl:outputs ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:coresMin 4 ; ResourceRequirement:ramMin 8000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/samtools-cwl:1.0.0" ] . a cwl:CommandLineTool ; rdfs:label "Mark duplicates and Sort" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( [ CommandLineBinding:position 2 ; cwl:valueFrom "$(runtime.cores)" ] [ CommandLineBinding:position 4 ; cwl:valueFrom "$(inputs.bam.nameroot).mark_dups_metrics.txt" ] ) ; cwl:baseCommand ( "/bin/bash" "markduplicates_helper.sh" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , ; cwl:outputs , ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/mark_duplicates-cwl:1.0.1" ] ; cwl:requirements [ a cwl:InitialWorkDirRequirement ; cwl:listing [ cwl:entry "set -o pipefail\nset -o errexit\n\ndeclare MD_BARCODE_TAG\nif [ ! -z \"$6\" ]; then\n MD_BARCODE_TAG=\"BARCODE_TAG=$6\"\n/usr/bin/java -Xmx16g -jar /opt/picard/picard.jar MarkDuplicates I=$1 O=/dev/stdout ASSUME_SORT_ORDER=$5 METRICS_FILE=$4 QUIET=true COMPRESSION_LEVEL=0 VALIDATION_STRINGENCY=LENIENT \"$MD_BARCODE_TAG\" | /usr/bin/sambamba sort -t $2 -m 18G -o $3 /dev/stdin\nelse\n /usr/bin/java -Xmx16g -jar /opt/picard/picard.jar MarkDuplicates I=$1 O=/dev/stdout ASSUME_SORT_ORDER=$5 METRICS_FILE=$4 QUIET=true COMPRESSION_LEVEL=0 VALIDATION_STRINGENCY=LENIENT | /usr/bin/sambamba sort -t $2 -m 18G -o $3 /dev/stdin\nfi\n" ; cwl:entryname "markduplicates_helper.sh" ] ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:coresMin 8 ; ResourceRequirement:ramMin 40000 ] . cwl:source . cwl:source . sld:type xsd:int , sld:null . cwl:outputSource ; sld:type cwl:File . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; sld:type xsd:string . cwl:in , , , , , , ; cwl:out , ; cwl:run . a cwl:CommandLineTool ; rdfs:label "collect WGS metrics" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.sample_name).WgsMetrics.txt" ] ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/usr/picard/picard.jar" "CollectWgsMetrics" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , ; cwl:outputs ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] . a cwl:CommandLineTool ; rdfs:label "apply BQSR" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "-o" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.output_name).bam" ] "-preserveQ" "6" "-SQQ" "10" "-SQQ" "20" "-SQQ" "30" "-nct" "8" "--disable_indel_quals" ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/opt/GenomeAnalysisTK.jar" "-T" "PrintReads" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , ; cwl:outputs ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/gatk-cwl:3.6.0" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] . cwl:source . cwl:default false . cwl:source . cwl:default false . cwl:source . cwl:source . cwl:source . cwl:source . cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:default "final" ; sld:type xsd:string , sld:null . cwl:default "final" ; sld:type xsd:string , sld:null . cwl:outputSource ; sld:type cwl:File . cwl:source . cwl:secondaryFiles [] ; sld:type [ cwl:inputBinding [ CommandLineBinding:position 4 ; CommandLineBinding:prefix "-knownSites" ] ; sld:items cwl:File ; sld:type sld:array ] . rdfs:comment "bwa-indexed reference file" ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:source ; cwl:valueFrom "$(self.sequence.hasOwnProperty('bam')? self.sequence.bam : null)" . cwl:inputBinding [ CommandLineBinding:position 1 ] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:source . cwl:outputSource ; sld:type cwl:File . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerBaseCoverage.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerBaseCoverage.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerBaseCoverage.txt" ] ; sld:type cwl:File , sld:null . cwl:source ; cwl:valueFrom "$(self.file)" . cwl:outputSource ; sld:type cwl:File . a cwl:CommandLineTool ; rdfs:label "samtools index" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "/opt/samtools/bin/samtools" "index" "$(runtime.outdir)/$(inputs.bam.basename)" "$(runtime.outdir)/$(inputs.bam.basename).bai" [ CommandLineBinding:shellQuote false ; cwl:valueFrom " && " ] "cp" "$(inputs.bam.basename).bai" "$(runtime.outdir)/$(inputs.bam.nameroot).bai" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs ; cwl:outputs ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/samtools-cwl:1.0.0" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 4000 ] ; cwl:requirements [ a cwl:ShellCommandRequirement ] ; cwl:requirements [ a cwl:InlineJavascriptRequirement ] ; cwl:requirements [ a cwl:InitialWorkDirRequirement ; cwl:listing "${ var f = inputs.bam; delete f.secondaryFiles; return f }" ] . cwl:source ; cwl:valueFrom "$(self.sequence.hasOwnProperty('fastq1')? self.sequence.fastq1 : null)" . sld:type cwl:File . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "I=" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . sld:fields , ; sld:type sld:record . cwl:default true . sld:type xsd:string . cwl:source ; cwl:valueFrom "$(self.file)" . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputSource ; sld:type cwl:File . cwl:valueFrom "ALL_READS" . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputSource ; sld:type cwl:File . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "TARGET_INTERVALS=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "TARGET_INTERVALS=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "TARGET_INTERVALS=" ] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "-2" ] ; sld:type cwl:File , sld:null . sld:type [ sld:items ; sld:type sld:array ] . cwl:source . sld:type [ sld:items ; sld:type sld:array ] . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.sample_name).WgsMetrics.txt" ] ; sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).InsertSizeHistogram.pdf" ] ; sld:type cwl:File . sld:type cwl:File . sld:type xsd:string , sld:null . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-HsMetrics.txt" ] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "R=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:in , , , , , , ; cwl:out , , , , ; cwl:run . cwl:outputSource ; sld:type cwl:File . cwl:source . rdfs:label "a file with a label" ; sld:fields , ; sld:type sld:record . cwl:default "chr21" , "chr16" , "chr13" , "chr15" , "chr4" , "chr6" , "chr3" , "chr1" , "chr18" , "chr10" , "chr17" , "chr8" , "chr20" , "chr22" , "chr5" , "chr12" , "chr7" , "chr19" , "chr11" , "chr2" , "chr14" , "chr9" ; cwl:inputBinding [ CommandLineBinding:position 1 ] ; sld:type [ cwl:inputBinding [ CommandLineBinding:prefix "-L" ] ; sld:items xsd:string ; sld:type sld:array ] . cwl:outputSource ; sld:type cwl:File . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputSource ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.basename)" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:in ; cwl:out ; cwl:run . cwl:source . cwl:in ; cwl:out ; cwl:run . cwl:inputBinding [ CommandLineBinding:prefix "REFERENCE_SEQUENCE=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:inputBinding [ CommandLineBinding:prefix "-g" ] ; sld:type xsd:string . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . a cwl:Workflow ; rdfs:label "Raw sequence data to BQSR" ; ns1:original_cwlVersion "v1.0" ; Workflow:steps , , , , , , ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , , , , ; cwl:outputs , ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types ] ; cwl:requirements [ a cwl:ScatterFeatureRequirement ] ; cwl:requirements [ a cwl:MultipleInputFeatureRequirement ] ; cwl:requirements [ a cwl:SubworkflowFeatureRequirement ] . cwl:default 1 ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_BASE_QUALITY=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_BASE_QUALITY=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_BASE_QUALITY=" ] ; sld:type xsd:int , sld:null . cwl:inputBinding [ CommandLineBinding:position 1 ] ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:inputBinding [ CommandLineBinding:prefix "METRIC_ACCUMULATION_LEVEL=" ] ; sld:type xsd:string . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.sample_name).GcBiasMetrics.txt" ] ; sld:type cwl:File . cwl:source . cwl:valueFrom "ALL_READS" . cwl:source . a cwl:CommandLineTool ; rdfs:label "samtools flagstat" ; ns1:original_cwlVersion "v1.0" ; cwl:baseCommand ( "/opt/samtools/bin/samtools" "flagstat" ) ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs ; cwl:outputs ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 4000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/samtools-cwl:1.0.0" ] ; cwl:stdout "$(inputs.bam.basename).flagstat" . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.sample_name).GcBiasMetricsSummary.txt" ] ; sld:type cwl:File . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:source . cwl:source . sld:fields ; sld:type sld:record . sld:type xsd:string . cwl:inputBinding [ CommandLineBinding:position 1 ] ; cwl:secondaryFiles [] ; sld:type cwl:File . sld:type sld:null ; sld:type [ sld:items xsd:string ; sld:type sld:array ] . a cwl:CommandLineTool ; rdfs:label "collect insert size metrics" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "O=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).InsertSizeMetrics.txt" ] "H=" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).InsertSizeHistogram.pdf" ] ) ; cwl:baseCommand ( "/usr/bin/java" "-Xmx16g" "-jar" "/usr/picard/picard.jar" "CollectInsertSizeMetrics" ) ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs , , ; cwl:outputs , ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 18000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/picard-cwl:2.18.1" ] . sld:type [ sld:items ; sld:type sld:array ] . sld:type xsd:int , sld:null . cwl:inputBinding [ CommandLineBinding:prefix "-b" ] ; sld:type cwl:File , sld:null . cwl:inputBinding [ CommandLineBinding:prefix "BAIT_INTERVALS=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "BAIT_INTERVALS=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "BAIT_INTERVALS=" ] ; sld:type cwl:File . cwl:in , ; cwl:out ; cwl:run ; cwl:scatter ( ) ; cwl:scatterMethod ScatterMethod:dotproduct . cwl:default false . cwl:source . cwl:source . cwl:default "MarkedSorted.bam" ; cwl:inputBinding [ CommandLineBinding:position 3 ] ; sld:type xsd:string , sld:null . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "refAlign.bam" ] ; sld:type cwl:File . cwl:source . cwl:inputBinding [ CommandLineBinding:position 1 ; CommandLineBinding:prefix "-R" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:source . cwl:default false . cwl:source . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:inputBinding [ CommandLineBinding:prefix "REFERENCE_SEQUENCE=" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:source . cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:outputSource ; sld:type cwl:File . cwl:outputSource ; sld:type cwl:File . a cwl:CommandLineTool ; rdfs:label "sort BAM by name" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "-t" [ cwl:valueFrom "$(runtime.cores)" ] "-m" "22G" "-n" "-o" [ cwl:valueFrom "$(inputs.bam.nameroot).NameSorted.bam" ] ) ; cwl:baseCommand ( "/usr/bin/sambamba" "sort" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs ; cwl:outputs ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/sambamba-cwl:0.6.4" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:coresMin 8 ; ResourceRequirement:ramMin 26000 ] . cwl:inputBinding [ CommandLineBinding:position 2 ; CommandLineBinding:prefix "-R" ] ; cwl:secondaryFiles [] ; cwl:secondaryFiles [] ; sld:type cwl:File , xsd:string . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.output_name)" ] ; cwl:secondaryFiles [] ; sld:type cwl:File . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).InsertSizeMetrics.txt" ] ; sld:type cwl:File . cwl:in , , , ; cwl:out ; cwl:run . cwl:source ; cwl:valueFrom "target-$(self.label)" . cwl:source . cwl:secondaryFiles [] ; sld:type cwl:File . cwl:source ; cwl:valueFrom "summary-$(self.label)" . sld:type sld:null ; sld:type [ sld:items xsd:string ; sld:type sld:array ] . a cwl:CommandLineTool ; rdfs:comment "Due to workflow runner limitations, use sequence_align_and_tag_adapter.cwl subworkflow to call this" ; rdfs:label "align with bwa_mem and tag" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( [ CommandLineBinding:position 5 ; CommandLineBinding:prefix "-n" ; cwl:valueFrom "$(runtime.cores)" ] ) ; cwl:baseCommand ( "/bin/bash" "sequence_alignment_helper.sh" ) ; cwl:cwlVersion ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:inputs , , , , ; cwl:outputs ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:coresMin 8 ; ResourceRequirement:ramMin 20000 ] ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/alignment_helper-cwl:1.0.0" ] ; cwl:requirements [ a cwl:InitialWorkDirRequirement ; cwl:listing [ cwl:entry "set -o pipefail\nset -o errexit\nset -o nounset\n\nwhile getopts \"b:?1:?2:?g:r:n:\" opt; do\n case \"$opt\" in\n b)\n MODE=bam\n BAM=\"$OPTARG\"\n ;;\n 1)\n MODE=fastq\n FASTQ1=\"$OPTARG\"\n ;;\n 2) \n MODE=fastq\n FASTQ2=\"$OPTARG\"\n ;;\n g)\n READGROUP=\"$OPTARG\"\n ;;\n r)\n REFERENCE=\"$OPTARG\"\n ;;\n n)\n NTHREADS=\"$OPTARG\"\n ;;\n esac\ndone\n\nif [[ \"$MODE\" == 'fastq' ]]; then\n /usr/local/bin/bwa mem -K 100000000 -t \"$NTHREADS\" -Y -R \"$READGROUP\" \"$REFERENCE\" \"$FASTQ1\" \"$FASTQ2\" | /usr/local/bin/samblaster -a --addMateTags | /opt/samtools/bin/samtools view -b -S /dev/stdin\nfi\nif [[ \"$MODE\" == 'bam' ]]; then\n /usr/bin/java -Xmx4g -jar /opt/picard/picard.jar SamToFastq I=\"$BAM\" INTERLEAVE=true INCLUDE_NON_PF_READS=true FASTQ=/dev/stdout | /usr/local/bin/bwa mem -K 100000000 -t \"$NTHREADS\" -Y -p -R \"$READGROUP\" \"$REFERENCE\" /dev/stdin | /usr/local/bin/samblaster -a --addMateTags | /opt/samtools/bin/samtools view -b -S /dev/stdin\nfi\n" ; cwl:entryname "sequence_alignment_helper.sh" ] ] ; cwl:stdout "refAlign.bam" . cwl:source . cwl:default "final" ; sld:type xsd:string . sld:type [ sld:items ; sld:type sld:array ] . cwl:source . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).NameSorted.bam" ] ; sld:type cwl:File . cwl:source . cwl:source . cwl:source . cwl:source . cwl:source . cwl:valueFrom "ALL_READS" . cwl:inputBinding [ CommandLineBinding:prefix "PER_BASE_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerBaseCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; cwl:inputBinding [ CommandLineBinding:prefix "PER_BASE_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerBaseCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; cwl:inputBinding [ CommandLineBinding:prefix "PER_BASE_COVERAGE=" ; cwl:valueFrom "${\n if(self) {\n return inputs.bam.nameroot + \".\" + inputs.output_prefix + \"-PerBaseCoverage.txt\"\n } else {\n return false;\n }\n}\n" ] ; sld:type xsd:boolean , sld:null . cwl:source . sld:type xsd:string . cwl:source . cwl:inputBinding [ CommandLineBinding:prefix "INTERVALS=" ] ; sld:type cwl:File , sld:null . cwl:source . cwl:in , , , , , , , , , , ; cwl:out , , , , , , , , , , , , , , ; cwl:run . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.name).merged.bam" ] ; sld:type cwl:File . a cwl:CommandLineTool ; rdfs:label "verify BAM ID" ; ns1:original_cwlVersion "v1.0" ; cwl:arguments ( "--out" [ cwl:valueFrom "$(runtime.outdir)/$(inputs.bam.nameroot).VerifyBamId" ] ) ; cwl:baseCommand ( "/usr/local/bin/verifyBamID" ) ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs , ; cwl:outputs , ; cwl:requirements [ a cwl:DockerRequirement ; DockerRequirement:dockerPull "mgibio/verify_bam_id-cwl:1.1.3" ] ; cwl:requirements [ a cwl:ResourceRequirement ; ResourceRequirement:ramMin 4000 ] . sld:type [ sld:items ; sld:type sld:array ] . cwl:inputBinding [ CommandLineBinding:position 3 ; CommandLineBinding:prefix "-BQSR" ] ; sld:type cwl:File . sld:type cwl:File . cwl:source . cwl:outputSource ; sld:type cwl:File . cwl:in , , , ; cwl:out ; cwl:run . cwl:default 1 ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_MAPPING_QUALITY=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_MAPPING_QUALITY=" ] ; cwl:inputBinding [ CommandLineBinding:prefix "MINIMUM_MAPPING_QUALITY=" ] ; sld:type xsd:int , sld:null . cwl:outputSource ; sld:type [ sld:items cwl:File ; sld:type sld:array ] . cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerTargetCoverage.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerTargetCoverage.txt" ] ; cwl:outputBinding [ CommandOutputBinding:glob "$(inputs.bam.nameroot).$(inputs.output_prefix)-PerTargetCoverage.txt" ] ; sld:type cwl:File , sld:null . cwl:source . cwl:outputSource ; sld:type cwl:File . cwl:default true . cwl:outputSource ; sld:type cwl:File . cwl:source . cwl:in , , , , , , , , , ; cwl:out , ; cwl:run ; cwl:scatter ( ) ; cwl:scatterMethod ScatterMethod:dotproduct . cwl:secondaryFiles [] ; sld:type cwl:File . sld:type [ sld:items ; sld:type sld:array ] . cwl:in , , , , , , , , , ; cwl:out , ; cwl:run ; cwl:scatter ( ) ; cwl:scatterMethod ScatterMethod:dotproduct . cwl:outputSource ; sld:type cwl:File . cwl:default "ALL_READS" ; sld:type xsd:string . cwl:source . cwl:in , , , ; cwl:out , , ; cwl:run . rdfs:label "sequence data with readgroup information" ; sld:fields , ; sld:type sld:record . a cwl:Workflow ; rdfs:label "HS Metrics workflow" ; ns1:original_cwlVersion "v1.0" ; Workflow:steps , , ; cwl:cwlVersion ; cwl:hints [ a ; LoadListingRequirement:loadListing "deep_listing" ] ; cwl:hints [ a ; NetworkAccess:networkAccess true ] ; cwl:inputs , , , , , , ; cwl:outputs , , , , ; cwl:requirements [ a cwl:SchemaDefRequirement ; SchemaDefRequirement:types ] ; cwl:requirements [ a cwl:ScatterFeatureRequirement ] ; cwl:requirements [ a cwl:StepInputExpressionRequirement ] ; cwl:requirements [ a cwl:SubworkflowFeatureRequirement ] . cwl:source . cwl:in ; cwl:out ; cwl:run .