Comments (7)
Hi @sharmaa9,
Can you try pointing the --vcf path to the vcf and not the tabix index file?
Would be the following command instead:
python bin/PGxPOP.py --vcf /Users/sharmaa9/Desktop/PGx/PGXPOP/PGxPOP/test_vcf.vcf.gz -g CYP2D6 --phased --build grch38 -o test_new.txt
from pgxpop.
Thanks for your response.
So I try that and i get this error-
File "/Users/sharmaa9/Desktop/PGx/PGXPOP/PGxPOP/PGxPOP/bin/DawgToys.py", line 320, in fetch_genotypes
records = tb.query("%s" % chr, position-1, position+1)
tabix.TabixError: query failed
Thanks,
Anu
from pgxpop.
Hi Anu,
Sorry to hear that didn't work. Can you give me a little more information?
What do the top few lines of your VCF look like after the column header? I'd like to see how chromosome names are formatted, as this sometimes differs among VCFs.
Have you successfully tabix query'd this file outside of PGxPOP?
Are all the chromosomes present in the file?
from pgxpop.
Thanks Alan,
Really appreciate your help. Actually, I am trying this with your example data file "VcfReaderTest-phasing.vcf". I bgzipped it and then tabixed it and now running PGxPOP on it.
from pgxpop.
Hi Anu,
Sorry, that test file doesn't contain CYP2D6. Can you try running:
python bin/PGxPOP.py --vcf /Users/sharmaa9/Desktop/PGx/PGXPOP/PGxPOP/test_vcf.vcf.gz -g DPYD --phased --build grch38 -o test_new.txt
from pgxpop.
Hi Adam,
Got it! Its working! :) thanks so much for your help! Really appreciate your time.
And thanks for coding this tool. Also, I am just starting in this field so I might bug you again in near future with queries.
Best,
Anu
from pgxpop.
Glad to hear its working and welcome to the field, it's a good one :)
from pgxpop.
Related Issues (12)
- Create more informative error message when a gene is missing from provided VCF
- CYP3A4 HOT 2
- Can PGxPOP handle unphased VCFs? HOT 1
- different PGxPOP outputs with same vcf input HOT 2
- extra comma HOT 1
- Everything returns as normal
- ModuleNotFoundError: No module named 'tabix' HOT 3
- Diplotype frequency calculation
- UnicodeDecodeError
- tuple index out of range error HOT 1
- Recommendation for updated allele definitions? HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pgxpop.