Skip to content

Update wrapper script for paired WGS pipeline (battenberg_wgs.R)#174

Open
s-basrai wants to merge 5 commits into
Wedge-lab:masterfrom
s-basrai:update-wgs-pipeline-wrapper
Open

Update wrapper script for paired WGS pipeline (battenberg_wgs.R)#174
s-basrai wants to merge 5 commits into
Wedge-lab:masterfrom
s-basrai:update-wgs-pipeline-wrapper

Conversation

@s-basrai
Copy link
Copy Markdown

@s-basrai s-basrai commented May 7, 2026

This PR addresses a few small problems that can cause errors when trying to run battenberg_wgs.R.

  • Set default value of --max_allowed_state to 250 (default used in battenberg()). This prevents failure when the user does not set the parameter.
  • Add parameter to explicitly specify the chromosome naming convention. The previous approach of using scanBamHeader() prevents the user from running the script on CRAM files. Default is set to use "chr" prefix.
  • Add parameter for specifying path to Battenberg reference files. Removes the need for hard-coded paths.
  • Change probloci file name to probloci.txt.gz to match what is in the current hg38 Battenberg reference files (https://doi.org/10.48420/30406441).

Additionally, a deprecated JVM flag (UseParallelOldGC) was being used in impute.R that caused the pipeline to fail for newer Java versions (https://bugs.openjdk.org/browse/JDK-8235860). This was corrected.

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.

1 participant