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

Multiallelic VCF with annotation in INFO field doesn't query all available infromation #301

Open
KalinNonchev opened this issue Mar 12, 2024 · 3 comments

Comments

@KalinNonchev
Copy link

KalinNonchev commented Mar 12, 2024

Hello,

I noticed that in the case of annotated multiallelic VCF where annotation for each alternative is stored in the INFO field, variant.INFO.get("annotation") would only pick the first occurrence.
Example:

1 123 A T,TA ... annotation:0.5.. annotation:0.1

variant.INFO.get("annotation") would return only 0.5 instead of [0.5, 0.1]

Best

@brentp
Copy link
Owner

brentp commented Mar 12, 2024

Hi, can you show the actual line from the VCF?
the INFO field is key=value, so you can't have the same key twice.
Format fields are separated by : as you have shown above.

@KalinNonchev
Copy link
Author

That would be an example with multiple gnomad312_AF annotations... It looks like the annotations are separated by "ALLELE_END;"

chr1\t10340\t.\tT\tG,C\t;ANNOVAR_DATE=2020-06-08;gnomad312_AF=0.0014;ALLELE_END;ANNOVAR_DATE=2020-06-08;gnomad312_AF=0.0509;ALLELE_END\t

@brentp
Copy link
Owner

brentp commented Mar 12, 2024

That is not VCF format. It might be a text format output by annovar?

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

2 participants