Entering edit mode
6.8 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