Skip to content

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