This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Discrepancy Between Ncbi / Vcf And Ncbi / Html

in the VCF for dbsnp137, rs11412589 is said to be located at chr9:74300311

$ curl -s "ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/00-All.vcf.gz" |gunzip -c |  grep BUILD
##dbSNP_BUILD_ID=137

~$ curl -s "ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/00-All.vcf.gz" |gunzip -c | cut -d '      ' -f 1-5 |  awk -F '    ' '($3=="rs11412589")'
9    74300311    rs11412589    T    TA,TAA

but on its page on dbsnp the position is shifted by ~20 bases (I know it's an indel but...) : chr9:74300329:74300330 .

as well as on the UCSC/mysql:

+-------+------------+----------+
| chrom | chromStart | chromEnd |
+-------+------------+----------+
| chr9  |   74300329 | 74300329 |
+-------+------------+----------+
1 row in set (0.19 sec)

is it me or is it a bug ?

Pierre

ncbi vcf dbsnp error
• 3,500 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

It seems to me a bug in VCF. The base at chr9:74300311 is "A", not "T". Actually there is no "T" around this region.

• 11 views • link
written by lh3 lh3

1 answer

here is the answer from the NCBI;

Hi,

The reason for this difference is that the VCF format checks for repeats. For variations happening within a repeat region, it will left justify the position to the base that is NOT a repeat unit.

For this variation, the flanking clearly indicates it is within a poly A region:

rs11412589 [Homo sapiens] 
TACTCCCTAAAAAAAAAAAAAAAAAA[-/A/AA]GAAAAAGAAAAAAAATCAATTTTTA

The BLAST alignment shows this more directly:

SNP  1         TACTCCCTAAAAAAAAAAAAAAAAAAWGAAAAAGAAAAAAAATCAATTTTTA  52
               |||||||||||||||||||||||||| |||||||||||||||||||||||||
ch9  74300304  TACTCCCTAAAAAAAAAAAAAAAAAA-GAAAAAGAAAAAAAATCAATTTTTA  74300354

in which the T is at 74300311.

Regards,

• 2 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I misread 74300311 in your original post as 74300331. No wonder I have not found any T... Thanks.

• 0 views • link
written by lh3 lh3

Log in to answer this question.

More posts like this
  • dbSNP vcf file is corrupted?
    written by HB • HB
    Hi, I downloaded the file 00-All.vcf.gz from https://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh38p7/VCF/ by using the wget command: wget ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh38p7/VCF/00-All.vcf.gz When I try to gunzip the file: gunzip -k 00-All.vcf.gz …
  • Reference and dbSNP incompatibility issue (MuTect2)
    written by umn_bist • umn_bist
    When I try using MuTect2 (from GATK) I get this error Is there a link to an (old) dbSNP that is compatible with UCSC's hg19 …
  • Did dbSNP not get the memo about hg38 chromosome names?
    written by Jeremy Leipzig Jeremy Leipzig
    <p>I thought everyone was switching over to &#39;chr1&#39; style chromosome names for hg38. Why do I still see the Ensembl-style chromosome names here?</p> <pre> tabix …
  • Visualizing Snp Data From Dbsnp In Igv
    written by axl058 • axl058
    <p>I am working on importing <a href='http://www.ncbi.nlm.nih.gov/projects/SNP/snp_viewBatch.cgi?sbid=1050082'>this SNP data</a> into IGV and have been running into problems left and right. I have downloaded the file …
  • Does The Vcffilter Tool Work On String Type Info Tags?
    written by dwaggott dwaggott
    <p>Hey,</p> <p>vcflib looks great! To get started I need a bit of direction on filtering using vcffilter. Specifically, does it work on string type info …
  • Problems Extracting Non-Snps From A Vcf File
    written by GPR GPR
    <p>Hello,</p> <p>In an SNP analysis, I am trying to extract those editing sites no found in the dbSNPs vcf file I have downloaded a couple …
  • Download Dbsnp135
    written by Raony Guimaraes • Raony Guimaraes
    <p>Where I can find dbsnp135 ?</p> <p>This address still not updated: <a href="ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/v4.0/00-All.vcf.gz">ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/v4.0/00-All.vcf.gz</a></p>
  • How To Decompress 1000Genomes Bgzip-Compressed Files Using Java
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>(cross posted on <a href='http://stackoverflow.com/questions/5180001'>stack-overflow</a> )</p> <p>Ok. I'm puzzled. </p> <p>I'm want to use java to download <a href='ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.sites.vcf.gz'>ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.sites.vcf.gz</a> in order to annotate some VCFs …
  • Dbsnp132 1000Genomes Vcf File Info Field: Bitfield Structure
    written by Biomed Biomed
    <p>The new 1000Genomes+dbSNP132 resource (<a href='ftp://ftp.ncbi.nlm.nih.gov/snp/organisms/human_9606/VCF/v4.0/ByChromosomeNoGeno/00-All.vcf.gz'>ftp://ftp.ncbi.nlm.nih.gov/snp/organisms/human_9606/VCF/v4.0/ByChromosomeNoGeno/00-All.vcf.gz</a>) has dbSNP byte structured information about snps.</p> <p>One such example is VP=050100000a0105051a000100 </p> <p>Considering the information found here <a …
  • Top/Bottom Confusion For Illumina Snp Calls
    written by Perry • Perry
    <p>Here is the problem: Illumina calls their SNPs AA,AB,BB. The meaning of A and B depend on what they call "top" or "bottom" strand. One …