Skip to content

Conversation

mike-w-wilson
Copy link
Contributor

Adds a function to create a VDS from GVCFs only. This is mainly a wrapper around hails function with limited utility. We can expand this function's utility to ingest VDS when the time comes and that is needed.

"""
Combine gVCFs into a single VDS.
:param gvcfs: Path to file containing gVCF paths with no header.
Copy link
Contributor

Choose a reason for hiding this comment

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

can drop the types after params since they are already defined above

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Ah copilot you silly AI

def create_vds(
gvcfs: str,
output_path: str,
temp_path: Optional[str] = None,
Copy link
Contributor

Choose a reason for hiding this comment

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

looks like temp path is required?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Ah good catch

@mike-w-wilson mike-w-wilson requested a review from klaricch July 17, 2024 18:41
save_path = temp_path + "combiner_plan.json"

gvcfs = read_list_data(gvcfs)
intervals = read_list_data(intervals) if intervals else None
Copy link
Contributor

Choose a reason for hiding this comment

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

can replace to use hl.import_locus_intervals and then collect the intervals here

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Thank you!

@mike-w-wilson mike-w-wilson requested a review from klaricch July 18, 2024 19:01
Copy link
Contributor

@klaricch klaricch left a comment

Choose a reason for hiding this comment

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

lgtm

@mike-w-wilson mike-w-wilson merged commit 0a91d10 into main Jul 23, 2024
@mike-w-wilson mike-w-wilson deleted the mw/create_vds branch July 23, 2024 13:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants