Skip to content

Allele frequency sum above 1 #241

Closed
Closed
@mbdabrowska1

Description

In my merge_output.vcf.gz file most of the time I'm getting just a single variant per position. However I also get one position where there are two variants, but their allele frequency does no add up to 1, nor do the allelic depths add up to the total read depth for the position. Am i misunderstanding the values or is there something wrong happening with the calculations?

A snippet of my merge_output.vcf.gz looks like this:

#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	SAMPLE
K03455.1	2229	.	C	A	22.46	PASS	F	GT:GQ:DP:AD:AF	1/1:22:19666:255,18257:0.9284
K03455.1	2259	.	G	A	36.43	PASS	F	GT:GQ:DP:AD:AF	1/1:36:19721:182,18951:0.9610
K03455.1	3161	.	G	A,T	18.93	PASS	F	GT:GQ:DP:AD:AF	1/2:18:20313:190,10274,17655:0.5058,0.8691

Where the last row is the one with multiple variants.

Any help would be greatly appreciated

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions