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

Issue querying db's made with latest cyvcf2 #69

Open
matthdsm opened this issue Mar 17, 2022 · 0 comments
Open

Issue querying db's made with latest cyvcf2 #69

matthdsm opened this issue Mar 17, 2022 · 0 comments

Comments

@matthdsm
Copy link
Contributor

Hi @brentp,

I've come across a weird issue.
I've updated one of our bcbio installations, resulting in an environment with the latest vcf2db + cyvcf2=0.30.14.
When I try to query one of the db's generated with this setup using the gemini python API, I get a funky result when fetching genotypes.

when I try to print the gts field from a table row (gemini.GeminiQuery.GeminiRow) I get the following numpy array

["T" "" "" "T" "" "" "C" "" "" "/" "" "" "T" "" "" "T" "" "" "C" "" "" "T"
 "" "" "/" "" "" "T" "" "" "T" "" "" "C" "" "" "" "" "" "" "" "" "" "" "T"
 "" "" "/" "" "" "T" "" "" "T" "" "" "C" "" "" "" "" "" "" "" "" "" "" ""]

which should show

["TTC/TTC","T/TTC","T/TTC"]

e.g. the genotypes for three individuals.
This is the case for our older installs running cyvcf2=0.20.9.

I suppose this error may have something to do with brentp/cyvcf2#227.
When downgrading cyvcf2 to the older version and regenerating the db, everything seems to work again.

Any thoughts?
M

PS: Not sure if this is the right repo to post this, so feel free to move this issue to a more relevant place.

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