I downloaded a VCF file from dbSNP, and I'm curious if the indels in the file are left-aligned for GRCH37 genome. The documentation doesn't say anything.
How can we tell if a VCF file has left or right aligned indels?
Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.Sign up to join this community
One possibility is to run bcftools norm on the file. At the end it will print out a statistic about how many variants were realigned. I did this for dbSNPv151 (GRCh38) only for chromosome 22 like this:
$ bcftools view -i 'TYPE="indel"' dbSNP_GRCh38.p7_b151.vcf.gz 22|bcftools norm -Ou -f GRCh38.fa > /dev/null Lines total/split/realigned/skipped: 901503/0/490/0
By definition in the VCF format specification the POS field is left aligned
POS int32 t 0-based leftmost coordinate
To confirm that you VCF is following the standard you would have to look at the reference and the match up the entries in the VCF. However if it is from a reliable source like dbSNP you can probably assume they follow the spec