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

VS-849 - Use the annotation 'AS_MQ' for indels. #8261

Merged
merged 24 commits into from
Mar 29, 2023

Conversation

gbggrant
Copy link
Collaborator

@gbggrant gbggrant commented Mar 23, 2023

Monitoring script removed.

Add monitoring script usage throughout VQSR Classic and Lite.
@gbggrant gbggrant changed the title Use the annotation 'AS_MQ' for indels. VS-849 - Use the annotation 'AS_MQ' for indels. Mar 23, 2023
@codecov
Copy link

codecov bot commented Mar 23, 2023

Codecov Report

❗ No coverage uploaded for pull request base (ah_var_store@291bfd0). Click here to learn what that means.
The diff coverage is n/a.

Additional details and impacted files
@@               Coverage Diff                @@
##             ah_var_store     #8261   +/-   ##
================================================
  Coverage                ?   85.873%           
  Complexity              ?     35517           
================================================
  Files                   ?      2194           
  Lines                   ?    167039           
  Branches                ?     18005           
================================================
  Hits                    ?    143442           
  Misses                  ?     17216           
  Partials                ?      6381           

@gbggrant gbggrant marked this pull request as ready for review March 24, 2023 15:11
@gbggrant
Copy link
Collaborator Author

GvsCreateFilterSet Run on the 3K here.
Precision and Sensitivity run here and results stored here.

@gbggrant gbggrant requested a review from mcovarr March 24, 2023 21:07
@mcovarr
Copy link
Collaborator

mcovarr commented Mar 24, 2023

WDLs failing validation 😿

use_allele_specific_annotations = true,
monitoring_script = "gs://gvs_quickstart_storage/cromwell_monitoring_script.sh"
Copy link
Collaborator

Choose a reason for hiding this comment

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

is this going to work with the Cromwell Local backend test?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Do you mean our integration test? I modified it to be a parameter to JointVcfFiltering.wdl so that that wdl's cromwell test will work. But now I realize that this might also break our integration tests?

Comment on lines 57 to 60
-resource:hapmap,known=false,training=true,truth=true,prior=0 ~{hapmap_resource_vcf} \
-resource:omni,known=false,training=true,truth=true,prior=0 ~{omni_resource_vcf} \
-resource:1000G,known=false,training=true,truth=false,prior=0 ~{one_thousand_genomes_resource_vcf} \
-resource:dbsnp,known=true,training=false,truth=false,prior=0 ~{dbsnp_resource_vcf}
Copy link
Collaborator

Choose a reason for hiding this comment

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

naive question: what do these priors mean here?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

That's temporary - need to do this to use P&S best. I will make a comment - or ideally parameterize this change.

Comment on lines 181 to 183
if [ -s ~{monitoring_script} ]; then
bash ~{monitoring_script} > monitoring.log &
fi
Copy link
Collaborator

Choose a reason for hiding this comment

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

If this evaluates to false there won't be a monitoring.log but the output is not optional

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Yup - thanks.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Changed this to return a non-optional (zero length file) if no monitoring script is passed.

Comment on lines 231 to 233
if [ -s ~{monitoring_script} ]; then
bash ~{monitoring_script} > monitoring.log &
fi
Copy link
Collaborator

Choose a reason for hiding this comment

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

same here

@gbggrant gbggrant requested a review from mcovarr March 28, 2023 20:50
@gbggrant
Copy link
Collaborator Author

Cleaned up and back to just the AS_MQ changes.

@gbggrant gbggrant merged commit 529b078 into ah_var_store Mar 29, 2023
@gbggrant gbggrant deleted the gg_VS-849_Add_AS_MQ_IndelAnnot_To_VQSR_Lite branch March 29, 2023 17:36
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants