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

VCFs with chromosome names prefixed by "chr" throwing errors #13

Open
agilly opened this issue May 22, 2020 · 0 comments
Open

VCFs with chromosome names prefixed by "chr" throwing errors #13

agilly opened this issue May 22, 2020 · 0 comments

Comments

@agilly
Copy link

agilly commented May 22, 2020

While running the following command:

locuszoom --build hg38 --metal smallregion.assoc.txt --markercol rs --pvalcol p_score  --ld-vcf region.vcf.gz --refsnp chr8:1234 --chr chr8

An error is thrown, despite all chromosome names and rs being coherent:

Using user-specified VCF file to calculate LD with reference SNP chr8:1234..
Error: while calculating LD from VCF file: index SNP position 1234 does not exist in file..

First of all I did not understand, so I modified vcf_ld.py to be more explicit by printing the tabixed region, the error was:

Error: while calculating LD from VCF file: tabix returned null for region  8:1234-1234.

This fails because the chromosomes in the VCF are prefixed by 'chr', as in, chr8.

It would be useful to add smart handling of "chr" prefixes.

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

No branches or pull requests

1 participant