Entering edit mode
5.9 years ago
kumajis
•
0
I want to annotate VCF with HGMD database information by sql, but there is some problem with indel records, such as :
In VCF the variant is
5 13868103 . TA T
, but in HGMD, it becomes5 13868104 A .
So for annotation correct, should I change variants like 5 13868103 . TA T
to 5 13868104 . A .
in VCF (if so, what tools could do this), or there are some other ways to solve it.
Thanks.
Hello 1097049297 ,
if possible you should try to make a vcf out of the sql database and normalize the result with bcftools norm.
fin swimmer