Ah! That's great, I didn't realise VCF's header provided that kind of detail...
Many thanks!
MrGraeme
Hi there,
I am working with a VCF that includes Copy Number Variants and hoped to confirm that I am reading the data correctly.
22 16050654 DUP_gs_CNV_22_16050654_16063474 A <CN0>,<CN2>,<CN3>,<CN4> 100 PASS AC=9,87,599,20;AF=0.00179712,0.0173722,0.119609,0.00399361;AN=5008;CS=DUP_gs;END=16063474;NS=2504;SVTYPE=CNV GT 3|0 0|0 0|0 0|0 0|0
From the line above extracted from my VCF I am assuming the following:
<CN3> copies of this block<CN3> means the block is copied 3 timesThe final assumption is my main concern, does <CNx> always mean there are x copies of the block?
If anyone can confirm that I am reading this correctly or otherwise set me straight I would be most grateful!
As an aside, as this is my first post, I'd love to thank the Biostars community for all the ongoing help and expertise! You guys are great!
Cheers
MrGraeme
You should find a definition of CN0,CN1, etc... in the VCF header in the ALT section.
Symbolic alternate alleles for imprecise structural variants:
##ALT=<ID=type,Description=description>
Ah! That's great, I didn't realise VCF's header provided that kind of detail...
Many thanks!
MrGraeme
Log in to answer this question.