-
Notifications
You must be signed in to change notification settings - Fork 589
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
Adding scripts to liftover gnomAD. Also bugfixes for Funcotator NIO. #5514
Conversation
* Added WDLs and some arguments to lift over gnomAD * Added chain file for b37->hg38 and arguments for liftover. * Limited to 1000 records in memory. * Added stack trace option to all wdls and sub tasks. * Fixed output to be consistent with local files for indexing. * Added timing information on wdls. * Added a wdl/json to create a TSV from gnomAD allele freq data. * Updated indexFeatureFile wdl, added params for run to index gnomAD. * Added json file for indexing a large gnomad file. * Fixed critical issues with NIO data sources. * Updates to the test script to save output and point to full cloud data. * Added some logging hooks to SeekableByteChannelPrefetcher
Codecov Report
@@ Coverage Diff @@
## master #5514 +/- ##
===============================================
- Coverage 87.076% 87.069% -0.006%
- Complexity 31320 31323 +3
===============================================
Files 1921 1921
Lines 144564 144577 +13
Branches 15951 15949 -2
===============================================
+ Hits 125880 125882 +2
- Misses 12895 12902 +7
- Partials 5789 5793 +4
|
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.
@jonn-smith If you answer my questions and make changes, feel free to merge.
@@ -0,0 +1,567 @@ | |||
#!/usr/bin/env bash |
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.
How was this verified to work? Did you give the results to the comms team to post somewhere?
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.
For now the resulting liftover file hasn't been validated. The liftover of gnomAD resulting from using this file is going to be QC'd, then it will have been verified.
I haven't yet given them the liftover file.
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.
I'm working with them now.
I haven't reviewed this since I made the changes to it to see what should stay, so it may need a fair bit of work.