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

Update output names for ReblockGVCF #1443

Open
wants to merge 14 commits into
base: develop
Choose a base branch
from
Open
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
6 changes: 3 additions & 3 deletions pipeline_versions.txt
Original file line number Diff line number Diff line change
Expand Up @@ -2,19 +2,19 @@ Pipeline Name Version Date of Last Commit
CheckFingerprint 1.0.22 2024-10-28
RNAWithUMIsPipeline 1.0.18 2024-11-04
AnnotationFiltration 1.2.7 2024-11-04
UltimaGenomicsWholeGenomeGermline 1.1.2 2024-11-04
UltimaGenomicsWholeGenomeGermline 1.1.3 2024-12-05
WholeGenomeGermlineSingleSample 3.3.3 2024-11-04
ExomeGermlineSingleSample 3.2.3 2024-11-04
JointGenotypingByChromosomePartTwo 1.5.2 2024-11-04
JointGenotypingByChromosomePartOne 1.5.2 2024-11-04
ReblockGVCF 2.3.2 2024-11-04
ReblockGVCF 2.4.0 2024-12-05
JointGenotyping 1.7.2 2024-11-04
UltimaGenomicsJointGenotyping 1.2.2 2024-11-04
VariantCalling 2.2.4 2024-11-04
UltimaGenomicsWholeGenomeCramOnly 1.0.23 2024-11-04
GDCWholeGenomeSomaticSingleSample 1.3.4 2024-11-04
BroadInternalRNAWithUMIs 1.0.36 2024-11-04
BroadInternalUltimaGenomics 1.1.2 2024-11-04
BroadInternalUltimaGenomics 1.1.3 2024-12-05
BroadInternalArrays 1.1.14 2024-11-04
BroadInternalImputation 1.1.14 2024-11-04
Arrays 2.6.30 2024-11-04
Expand Down
Original file line number Diff line number Diff line change
@@ -1,3 +1,8 @@
# 2.4.0
2024-12-05 (Date of Last Commit)

* Updated output names for ReblockGVCF workflow from output_vcf and output_vcf_index to reblocked_gvcf and reblocked_gvcf_index respectively

# 2.3.2
2024-11-04 (Date of Last Commit)

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ import "../../../../../../tasks/broad/Utilities.wdl" as utils

workflow ReblockGVCF {

String pipeline_version = "2.3.2"
String pipeline_version = "2.4.0"


input {
Expand Down Expand Up @@ -70,8 +70,8 @@ workflow ReblockGVCF {
}

output {
File output_vcf = Reblock.output_vcf
File output_vcf_index = Reblock.output_vcf_index
File reblocked_gvcf = Reblock.output_vcf
File reblocked_gvcf_index = Reblock.output_vcf_index
}
meta {
allowNestedInputs: true
Expand Down
Original file line number Diff line number Diff line change
@@ -1,3 +1,8 @@
# 1.1.3
2024-12-05 (Date of Last Commit)

* Updated the name of the output for ReblockGVCFs; this does not affect this pipeline

# 1.1.2
2024-11-04 (Date of Last Commit)

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -50,7 +50,7 @@ workflow UltimaGenomicsWholeGenomeGermline {
filtering_model_no_gt_name: "String describing the optional filtering model; default set to rf_model_ignore_gt_incl_hpol_runs"
}

String pipeline_version = "1.1.2"
String pipeline_version = "1.1.3"


References references = alignment_references.references
Expand Down Expand Up @@ -202,8 +202,8 @@ workflow UltimaGenomicsWholeGenomeGermline {

# Outputs that will be retained when execution is complete
output {
File output_gvcf = ReblockGVCF.output_vcf
File output_gvcf_index = ReblockGVCF.output_vcf_index
File output_gvcf = ReblockGVCF.reblocked_gvcf

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do we need to change the outputs here to? I don't think we (GVS) ever direct consume the output of this WDL, so it's likely that maintaining its output naming scheme was deliberate. But I just wanted to verify that here

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This change is necessary since the output of the reblock workflow is being referenced. I left the output name for this workflow the same, but needed to update the reference to the output of ReblockGVCF

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

yup, totally makes sense why this change was needed. As I don't know the users of this particular wdl, I mostly wanted to ensure that it wasn't something any end users would be running before sending the data into GVS (in which case we'd want to update its outputs as well). It sounds as though this is not expected behavior though, and this one has only been touched because it consumes ReblockGVCF for other reasons. Thanks!

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It's a good question. I'm not sure whether we expect the outputs of this workflow to be used in joint calling.
@rickymagner would you be able to comment on whether we expect the outputs of UltimaGenomicsWholeGenomeGermline to be used for joint calling in GVS?

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hm I don't know if GVS has been used on Ultima samples before. Does anyone from the Variants team know if that's an existing/planned application for GVS?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks @rickymagner!
Sounds like this isn't currently a known issue @koncheto-broad, let me know if this becomes an issue in the future and we can change it here too.

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for clarifying, y'all! I think we pretty exclusively use Illumina samples at this point, so the Ultima pipeline is unlikely to affect us. And if we discover that it does, we can always just reach out then!

File output_gvcf_index = ReblockGVCF.reblocked_gvcf_index
File output_vcf = ConvertGVCFtoVCF.output_vcf
File output_vcf_index = ConvertGVCFtoVCF.output_vcf_index

Expand Down
Original file line number Diff line number Diff line change
@@ -1,3 +1,8 @@
# 1.1.3
2024-12-05 (Date of Last Commit)

* Updated the name of the output for ReblockGVCFs; this does not affect this pipeline

# 1.1.2
2024-11-04 (Date of Last Commit)

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ import "../../../../../../../pipelines/broad/qc/CheckFingerprint.wdl" as FP

workflow BroadInternalUltimaGenomics {

String pipeline_version = "1.1.2"
String pipeline_version = "1.1.3"

input {

Expand Down
8 changes: 4 additions & 4 deletions verification/test-wdls/TestReblockGVCF.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -54,8 +54,8 @@ workflow TestReblockGVCF {
# Collect all of the pipeline outputs into single Array[String]
Array[String] pipeline_outputs = flatten([
[ # File outputs
ReblockGVCF.output_vcf_index,
ReblockGVCF.output_vcf,
ReblockGVCF.reblocked_gvcf_index,
ReblockGVCF.reblocked_gvcf,
],

])
Expand Down Expand Up @@ -86,13 +86,13 @@ workflow TestReblockGVCF {
if (!update_truth){
call Utilities.GetValidationInputs as GetGvcf {
input:
input_file = ReblockGVCF.output_vcf,
input_file = ReblockGVCF.reblocked_gvcf,
results_path = results_path,
truth_path = truth_path
}
call Utilities.GetValidationInputs as GetGvcfIndex {
input:
input_file = ReblockGVCF.output_vcf_index,
input_file = ReblockGVCF.reblocked_gvcf_index,
results_path = results_path,
truth_path = truth_path
}
Expand Down
Loading