how to change indel to smallest length in vcf
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.
• 1,294 views
•
link
0 answers
No answers yet.
Log in to answer this question.
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