Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add PlantUML flow diagram #157

Merged
merged 7 commits into from
Apr 10, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension


Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion .github/workflows/render-puml.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,7 @@ jobs:

steps:
- name: Generate PUML diagrams
uses: uclahs-cds/[email protected].0
uses: uclahs-cds/[email protected].1
with:
github-token: ${{ secrets.GITHUB_TOKEN }}
ghcr-username: ${{ github.actor }}
Expand Down
3 changes: 3 additions & 0 deletions CHANGELOG.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,9 @@ This project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.htm
---

## [Unreleased]
### Changed
- Replace workflow diagram with PlantUML version
- Update PlantUML action to v1.0.1

---

Expand Down
2 changes: 1 addition & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -42,7 +42,7 @@ python submit_nextflow_pipeline.py \

## Flow Diagram

![call-gSNP flow diagram](call-gSNP-DSL2.png)
![call-gSNP flow diagram](docs/call-gsnp-flow.svg)

---

Expand Down
Binary file removed call-gSNP-DSL2.png
Binary file not shown.
87 changes: 87 additions & 0 deletions docs/call-gsnp-flow.puml
Original file line number Diff line number Diff line change
@@ -0,0 +1,87 @@
@startuml

skinparam ConditionEndStyle hline

start

if (Explicit intervals?) is (Yes) then
:==run_SplitIntervals_GATK
----
Split reference genome into up
to **scatter_count** interval lists,
without subdividing any of the
input intervals;
else (No)
:==run_SplitIntervals_GATK
----
Split reference genome into
**scatter_count** interval lists;
endif

split

:==run_HaplotypeCallerVCF_GATK
----
Generate VCFs for each split interval
using HaplotypeCaller;

:==run_MergeVcfs_Picard_VCF
----
Merge raw variants into whole VCF file;

#palegreen:Per-sample raw VCF + index files>

partition "Recalibrate Variants" {

:==run_VariantRecalibratorSNP_GATK
----
Generate Variant Quality Score Recalibration
(VQSR) table for SNPs;

:==run_ApplyVQSR_GATK_SNP
----
Filter SNPs based on VQSR table;

:==run_VariantRecalibratorINDEL_GATK
----
Generate VQSR table for INDELs;

:==run_ApplyVQSR_GATK_INDEL
----
Filter INDELs based on VQSR table;

}

#palegreen:SNP and INDEL recalibrated
variants + index files>

:==filter_gSNP_GATK
----
Filter ambiguous variants with
customized Perl script;

#palegreen:Filtered germline
variants + index files>

detach

split again

:==run_HaplotypeCallerGVCF_GATK
----
Generate GVCFs for each split interval
using HaplotypeCaller;

:==run_MergeVcfs_Picard_GVCF
----
Merge raw variants into whole GVCF file;

#palegreen:Per-sample GVCF + index files>

detach

endsplit


@enduml

96 changes: 96 additions & 0 deletions docs/call-gsnp-flow.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading