Explore Workflows
View already parsed workflows here or click here to add your own
| Graph | Name | Retrieved From | View |
|---|---|---|---|
|
|
snaptools_create_snap_file.cwl
|
Path: steps/snaptools_create_snap_file.cwl Branch/Commit ID: bb023f95ca3330128bfef41cc719ffcb2ee6a190 |
|
|
|
Dockstore.cwl
|
Path: Dockstore.cwl Branch/Commit ID: b9eb2dcd745bf3cdd6e9c3396a507e11d60c87fe |
|
|
|
mutect panel-of-normals workflow
|
Path: definitions/pipelines/panel_of_normals.cwl Branch/Commit ID: ddb49a0951d9ad537269d7db3fe8f904495a8bf4 |
|
|
|
Variant calling germline paired-end
A workflow for the Broad Institute's best practices gatk4 germline variant calling pipeline. ## __Outputs__ #### Primary Output files: - bqsr2_indels.vcf, filtered and recalibrated indels (IGV browser) - bqsr2_snps.vcf, filtered and recalibrated snps (IGV browser) - bqsr2_snps.ann.vcf, filtered and recalibrated snps with effect annotations #### Secondary Output files: - sorted_dedup_reads.bam, sorted deduplicated alignments (IGV browser) - raw_indels.vcf, first pass indel calls - raw_snps.vcf, first pass snp calls #### Reports: - overview.md (input list, alignment metrics, variant counts) - insert_size_histogram.pdf - recalibration_plots.pdf - snpEff_summary.html ## __Inputs__ #### General Info - Sample short name/Alias: unique name for sample - Experimental condition: condition, variable, etc name (e.g. \"control\" or \"20C 60min\") - Cells: name of cells used for the sample - Catalog No.: vender catalog number if available - BWA index: BWA index sample that contains reference genome FASTA with associated indices. - SNPEFF database: Name of SNPEFF database to use for SNP effect annotation. - Read 1 file: First FASTQ file (generally contains \"R1\" in the filename) - Read 2 file: Paired FASTQ file (generally contains \"R2\" in the filename) #### Advanced - Ploidy: number of copies per chromosome (default should be 2) - SNP filters: see Step 6 Notes: https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/ - Indel filters: see Step 7 Notes: https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/ #### SNPEFF notes: Get snpeff databases using `docker run --rm -ti gatk4-dev /bin/bash` then running `java -jar $SNPEFF_JAR databases`. Then, use the first column as SNPEFF input (e.g. \"hg38\"). - hg38, Homo_sapiens (USCS), http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_hg38.zip - mm10, Mus_musculus, http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_mm10.zip - dm6.03, Drosophila_melanogaster, http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_dm6.03.zip - Rnor_6.0.86, Rattus_norvegicus, http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_Rnor_6.0.86.zip - R64-1-1.86, Saccharomyces_cerevisiae, http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_R64-1-1.86.zip ### __Data Analysis Steps__ 1. Trimming the adapters with TrimGalore. - This step is particularly important when the reads are long and the fragments are short - resulting in sequencing adapters at the ends of reads. If adapter is not removed the read will not map. TrimGalore can recognize standard adapters, such as Illumina or Nextera/Tn5 adapters. 2. Generate quality control statistics of trimmed, unmapped sequence data 3. Run germline variant calling pipeline, custom wrapper script implementing Steps 1 - 17 of the Broad Institute's best practices gatk4 germline variant calling pipeline (https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/) ### __References__ 1. https://gencore.bio.nyu.edu/variant-calling-pipeline-gatk4/ 2. https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels- 3. https://software.broadinstitute.org/software/igv/VCF |
Path: workflows/vc-germline-pe.cwl Branch/Commit ID: 12e5256de1b680c551c87fd5db6f3bc65428af67 |
|
|
|
Bacterial Annotation, pass 2, blastp-based functional annotation (first pass)
|
Path: bacterial_annot/wf_bacterial_annot_pass2.cwl Branch/Commit ID: 22ffe27d9d4a899def7592d75d5871c1856adbdb |
|
|
|
step_valuefrom5_wf_v1_1.cwl
|
Path: testdata/step_valuefrom5_wf_v1_1.cwl Branch/Commit ID: 77669d4dd1d1ebd2bdd9810f911608146d9b8e51 |
|
|
|
Trim Galore RNA-Seq pipeline single-read
The original [BioWardrobe's](https://biowardrobe.com) [PubMed ID:26248465](https://www.ncbi.nlm.nih.gov/pubmed/26248465) **RNA-Seq** basic analysis for a **single-end** experiment. A corresponded input [FASTQ](http://maq.sourceforge.net/fastq.shtml) file has to be provided. Current workflow should be used only with the single-end RNA-Seq data. It performs the following steps: 1. Trim adapters from input FASTQ file 2. Use STAR to align reads from input FASTQ file according to the predefined reference indices; generate unsorted BAM file and alignment statistics file 3. Use fastx_quality_stats to analyze input FASTQ file and generate quality statistics file 4. Use samtools sort to generate coordinate sorted BAM(+BAI) file pair from the unsorted BAM file obtained on the step 1 (after running STAR) 5. Generate BigWig file on the base of sorted BAM file 6. Map input FASTQ file to predefined rRNA reference indices using Bowtie to define the level of rRNA contamination; export resulted statistics to file 7. Calculate isoform expression level for the sorted BAM file and GTF/TAB annotation file using GEEP reads-counting utility; export results to file |
Path: workflows/trim-rnaseq-se.cwl Branch/Commit ID: 4f48ee6f8665a34cdf96e89c012ee807f80c7a3d |
|
|
|
trim-chipseq-pe.cwl
ChIP-Seq basic analysis workflow for a paired-end experiment with Trim Galore. |
Path: workflows/trim-chipseq-pe.cwl Branch/Commit ID: e284e3f6dff25037b209895c52f2abd37a1ce1bf |
|
|
|
step_valuefrom5_wf_with_id_v1_1.cwl
|
Path: testdata/step_valuefrom5_wf_with_id_v1_1.cwl Branch/Commit ID: 77669d4dd1d1ebd2bdd9810f911608146d9b8e51 |
|
|
|
revsort.cwl
|
Path: cwl/revsortlcase/revsort.cwl Branch/Commit ID: 19dee894600a1d265084f9187bfc1fe06b9f460b |
