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

Not every line was annotated? #49

Open
xiucz opened this issue Sep 13, 2022 · 2 comments
Open

Not every line was annotated? #49

xiucz opened this issue Sep 13, 2022 · 2 comments

Comments

@xiucz
Copy link

xiucz commented Sep 13, 2022

Hi,

I use duphold to calculate SV depth, but I find that not every vcf line was annotated with the tags, how to interpret this situation?

chr17	57986589	7_2	N	[chr17:41256075[N	0	.	SVTYPE=BND;STRANDS=--:8;IMPRECISE;CIPOS=-237,5;CIEND=-207,9;CIPOS95=-48,1;CIEND95=-29,3;MATEID=7_1;EVENT=7;SECONDARY;SU=8;PE=8;SR=0;GCF=0.437831	GT:SU:PE:SR:GQ:SQ:GL:DP:RO:AO:QR:QA:RS:AS:ASC:RP:AP:AB:DHFC:DHFFC:DHBFC	./.:8:8:0:.:.:.:0:0:0:0:0:0:0:0:0:0:.:0:0:0
chr17	59934382	9_2	N	]chr1:59536384]N	0	.	SVTYPE=BND;STRANDS=-+:4;CIPOS=-9,8;CIEND=-9,8;CIPOS95=0,0;CIEND95=0,0;MATEID=9_1;EVENT=9;SECONDARY;SU=4;PE=0;SR=4	GT:SU:PE:SR:GQ:SQ:GL:DP:RO:AO:QR:QA:RS:AS:ASC:RP:AP:AB	./.:4:0:4:.:.:.:0:0:0:0:0:0:0:0:0:0:.

Best,
xiucz

@brentp
Copy link
Owner

brentp commented Sep 13, 2022

Duphold doesn't annotate every BND, if the end is close enough it will. In this case, it looks like there were 0 reads with sufficient quality to evaluate the depth.

@xiucz
Copy link
Author

xiucz commented Sep 15, 2022

@brentp

Thanks for your quick reply, what does the the end is close enough mean?

And one more question, the three tags, DHFC, DHFFC, DHBFC, which one refers to the VAF value? Thank you.

best,
xiucz

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