Hi everyone,
in order to get the SNPs of chromosome 1 ( homo sapiens,) I went to BioMart page on the site ensembl.org. I selected "Ensembl Variation 83"/"Homo Sapiens Short Variants(SNPs and indels excluding flanked variants)" as dataset and set the filter to chromosome 1. The data file I downloaded contains 1 488 472 records all with different rs numbers and different chromosome positions. As far as I know, that the chromosome 1 should have about 130K SNPs. So, my question is, did I do something wrong by obtaining the data? How should I reconcile this large discrepancy?
Thanks in advance
James
•
written
by
jamesxli2007 •
1 answer
I expect your 130k prediction is based on old data. There should be 11.6M short variants on human chr1. Perhaps your prediction is pre-1000 Genomes?
Note that BioMart is not capable of dealing with this volume of data, hence getting only 10% of the variants you ought to. I suggest downloading the VCF files and using tabix to get out only chr1.
Log in to answer this question.
-
How to use a variable to install an R package
written by nonaddldy •
-
Get gene names from rs SNP ids for mmusculus_gene_ensembl
written by julian.d.stamp •
-
Best practice for running GATK VQSR on X chromosome
written by samuelandjw •
-
Why do different users have different number of SNPs on the same 23andMe genotyping chip?
written by mprza •
-
How to tell a sequence found in NCBI ("Nucleotide") is from Human Reference Genome or not.
written by freeheart.steven •
-
dbSNP build in ENSEMBL biomaRt?
written by 09meh •
-
Query related SNP merger
written by madhura.b8 •
-
How to get a unique target gene for each SNP from different chromosomes/genes in a list at once?
written by cmcouto.silva •
-
Ensembl Variation REST API Phased Haplotype?
written by johntimcham •
-
How To Retrieve A Gff With Sequence Identifier Different Than Gb With Bioperl
written by fbrundu
