-
Notifications
You must be signed in to change notification settings - Fork 31
Add create_vds function that only supports creating from gvcfs #716
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
Conversation
gnomad/utils/file_utils.py
Outdated
""" | ||
Combine gVCFs into a single VDS. | ||
:param gvcfs: Path to file containing gVCF paths with no header. |
There was a problem hiding this comment.
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
There was a problem hiding this comment.
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
gnomad/utils/file_utils.py
Outdated
def create_vds( | ||
gvcfs: str, | ||
output_path: str, | ||
temp_path: Optional[str] = None, |
There was a problem hiding this comment.
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?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ah good catch
gnomad/utils/file_utils.py
Outdated
save_path = temp_path + "combiner_plan.json" | ||
|
||
gvcfs = read_list_data(gvcfs) | ||
intervals = read_list_data(intervals) if intervals else None |
There was a problem hiding this comment.
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
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thank you!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
lgtm
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.