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

Enable use of indices for large genomes #151

Open
jd3234 opened this issue Apr 22, 2020 · 1 comment
Open

Enable use of indices for large genomes #151

jd3234 opened this issue Apr 22, 2020 · 1 comment

Comments

@jd3234
Copy link

jd3234 commented Apr 22, 2020

This might be a general issue for medaka but I have only used medaka_variant so far.

Is your feature request related to a problem? Please describe.
currently medaka only looks for bai etc indices. However, large genomes can't be stored in the default indices but only in csi

However, whatshap can also only handle small genomes as it expects vcf files indexed as tbi. They would also need to enable the use of vcf.gz.csi (I just ended up renaming the file and then restarted medaka to resume at the whatshap step again)

@cjw85
Copy link
Member

cjw85 commented Apr 23, 2020

Hi @jd3234,

As a piece of research software, I suspect this isn't something we will have time to implement. If you'd like to send a pull request making the changes necessary to output and use csi indexes I'm sure we can accomodate.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants