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

Provide a bait file to Picard's CollectHsMetrics #3695

Open
wants to merge 2 commits into
base: master
Choose a base branch
from

Conversation

asee-imagia
Copy link

Hi!

I'm interested in providing proper bait files to Picard's CollectHsMetrics tool. From the code that I examined, bcbio currently forces both the bait file and the target file to be the same (eg: the file provided to variant_regions). This PR proposes to fix the small missing link that would allow providing a proper bait file.

  1. It creates a config:algorithm:bait_regions entry in the datadict, which allows for providing a path to the bait file in the yaml config.
  2. It gets the reference to this bait file properly passed over to the reports method of the PicardMetrics instance.

If anything is missing, I will gladly adapt this contribution as needed (updating the documentation, adding other missing elements that I oversaw, etc.).

Thank you!

Allow the ``bait_regions`` in the datadict, so that a path to a BED file can be specified. This will eventually be passed picard's CollectHsMetrics. (This part is not yet implemented).
Retrieve the bait_regions file from the datadict, and passed it on to picard.
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