@prefix LoadListingRequirement: .
@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:position 3 ;
CommandLineBinding:prefix "-I"
] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:source .
cwl:source .
cwl:source .
cwl:inputBinding [ CommandLineBinding:position 2 ;
CommandLineBinding:prefix "-I"
] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.bam.basename)" ] ;
cwl:secondaryFiles [] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
a cwl:Workflow ;
rdfs:label "FASTQ to BQSR" ;
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:ScatterFeatureRequirement ] ;
cwl:requirements [ a cwl:MultipleInputFeatureRequirement ] .
a cwl:CommandLineTool ;
rdfs:label "align with bwa_mem and tag" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( [ CommandLineBinding:position 5 ;
cwl:valueFrom "$(runtime.cores)"
]
) ;
cwl:baseCommand ( "/bin/bash" "fastq_alignment_helper.sh" ) ;
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:coresMin 8 ;
ResourceRequirement:ramMin 20000
] ;
cwl:requirements [ a cwl:InitialWorkDirRequirement ;
cwl:listing [ cwl:entry "set -o pipefail\nset -o errexit\nset -o nounset\n\nFASTQ1=\"$1\"\nFASTQ2=\"$2\"\nREADGROUP=\"$3\"\nREFERENCE=\"$4\"\nNTHREADS=\"$5\"\n\n/usr/local/bin/bwa mem -K 100000000 -t \"$NTHREADS\" -Y -p -R \"$READGROUP\" \"$REFERENCE\" \"$FASTQ1\" \"$FASTQ2\" | /usr/local/bin/samblaster -a --addMateTags | /opt/samtools/bin/samtools view -b -S /dev/stdin\n" ;
cwl:entryname "fastq_alignment_helper.sh"
]
] ;
cwl:requirements [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/alignment_helper-cwl:1.0.0"
] ;
cwl:stdout "refAlign.bam" .
rdfs:comment "readgroup information for the BAM header" ;
cwl:inputBinding [ CommandLineBinding:position
3 ] ;
sld:type xsd:string .
cwl:source .
cwl:outputBinding [ CommandOutputBinding:glob
"bqsr.table" ] ;
sld:type cwl:File .
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-R"
] ;
sld:type xsd:string .
rdfs:comment "fastq of sequence data (first read if paired end)" ;
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type cwl:File .
cwl:source .
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:source .
cwl:inputBinding [ CommandLineBinding:position 2 ;
CommandLineBinding:prefix "-R"
] ;
sld:type xsd:string .
cwl:source .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.bam.nameroot).NameSorted.bam" ] ;
sld:type cwl:File .
sld:type [ sld:items xsd:string ;
sld:type sld:array
] .
sld:type sld:null ;
sld:type [ sld:items xsd:string ;
sld:type sld:array
] .
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type cwl:File .
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:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/gatk-cwl:3.6.0"
] ;
cwl:requirements [ a cwl:ResourceRequirement ;
ResourceRequirement:ramMin 18000
] .
cwl:source .
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 ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs ;
cwl:outputs ;
cwl:requirements [ a cwl:InitialWorkDirRequirement ;
cwl:listing "${ var f = inputs.bam; delete f.secondaryFiles; return f }"
] ;
cwl:requirements [ a cwl:ResourceRequirement ;
ResourceRequirement:ramMin 4000
] ;
cwl:requirements [ a cwl:ShellCommandRequirement ] ;
cwl:requirements [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/samtools-cwl:1.0.0"
] ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.name).merged.bam" ] ;
sld:type cwl:File .
sld:type cwl:File .
cwl:source .
cwl:in , , , ;
cwl:out ;
cwl:run .
cwl:outputSource ;
cwl:secondaryFiles [] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type cwl:File .
cwl:source .
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
cwl:source .
cwl:default "Final.bam" ;
sld:type xsd:string , sld:null .
cwl:default "chr9" , "chr21" , "chr2" , "chr13" , "chr6" , "chr14" , "chr17" , "chr7" , "chr19" , "chr18" , "chr8" , "chr5" , "chr22" , "chr20" , "chr11" , "chr10" , "chr12" , "chr15" , "chr16" , "chr1" , "chr3" , "chr4" ;
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type [ cwl:inputBinding [ CommandLineBinding:prefix
"-L" ] ;
sld:items xsd:string ;
sld:type sld:array
] .
cwl:inputBinding [ CommandLineBinding:position 3 ;
CommandLineBinding:prefix "-BQSR"
] ;
sld:type cwl:File .
cwl:source .
rdfs:comment "second read fastq for paired end data" ;
cwl:inputBinding [ CommandLineBinding:position
2 ] ;
sld:type cwl:File .
a cwl:CommandLineTool ;
rdfs:label "Samtools: merge" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( "$(inputs.name).merged.bam" ) ;
cwl:baseCommand ( "/opt/samtools/bin/samtools" "merge" ) ;
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:coresMin 4 ;
ResourceRequirement:ramMin 8000
] ;
cwl:requirements [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/samtools-cwl:1.0.0"
] .
cwl:outputBinding [ CommandOutputBinding:glob
"refAlign.bam" ] ;
sld:type cwl:File .
cwl:in ;
cwl:out ;
cwl:run .
cwl:in ;
cwl:out , ;
cwl:run .
cwl:in , , , ;
cwl:out ;
cwl:run .
cwl:outputSource ;
sld:type cwl:File .
cwl:source .
cwl:source , , .
a cwl:CommandLineTool ;
rdfs:label "sort BAM by name" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( "-t"
[ cwl:valueFrom "$(runtime.cores)" ]
"-m" "12G" "-n" "-o"
[ cwl:valueFrom "$(inputs.bam.nameroot).NameSorted.bam" ]
) ;
cwl:baseCommand ( "/usr/bin/sambamba" "sort" ) ;
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:coresMin 8 ;
ResourceRequirement:ramMin 16000
] ;
cwl:requirements [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/sambamba-cwl:0.6.4"
] .
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 ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs , , ;
cwl:outputs , ;
cwl:requirements [ a cwl:ResourceRequirement ;
ResourceRequirement:coresMin 8 ;
ResourceRequirement:ramMin 40000
] ;
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:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/mark_duplicates-cwl:1.0.1"
] .
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
a cwl:CommandLineTool ;
rdfs:label "apply BQSR" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( "-o"
[ cwl:valueFrom "$(runtime.outdir)/$(inputs.output_name)" ]
"-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:ResourceRequirement ;
ResourceRequirement:ramMin 18000
] ;
cwl:requirements [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "mgibio/gatk-cwl:3.6.0"
] .
cwl:source .
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:in , , , ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) ;
cwl:scatterMethod ScatterMethod:dotproduct .
cwl:in ;
cwl:out ;
cwl:run .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.bam.nameroot).mark_dups_metrics.txt" ] ;
sld:type cwl:File .
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:secondaryFiles [] ;
sld:type [ cwl:inputBinding [ CommandLineBinding:position 4 ;
CommandLineBinding:prefix "-knownSites"
] ;
sld:items cwl:File ;
sld:type sld:array
] .
cwl:source .
sld:type xsd:string .
rdfs:comment "bwa-indexed reference file" ;
cwl:inputBinding [ CommandLineBinding:position
4 ] ;
sld:type xsd:string .
sld:type xsd:string .
cwl:default "final.bam" ;
sld:type xsd:string , sld:null .
cwl:default "queryname" ;
cwl:inputBinding [ CommandLineBinding:position
5 ] ;
sld:type xsd:string .
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
cwl:source .
cwl:default "MarkedSorted.bam" ;
cwl:inputBinding [ CommandLineBinding:position
3 ] ;
sld:type xsd:string , sld:null .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.output_name)" ] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.output_name)" ] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:in , ;
cwl:out ;
cwl:run .