I want to get allele frequencies of Illumina Omni2.5 in European population only and my file contains information on chr, pos, allele1, and allele2. I though there could be two ways of doing it:
- Annotate Reference and Alternate allele for each SNP (how??) and use Annovar to get 1000 genome population specific allele frequency
- Download allele frequency information from HapMap/1000 Genome directly.
Does any one has any good suggestions on how to go about it or recommend any software?
Thanks in advance!
1 answer
SeattleSNP can use a "custom" file format that will accept the FinalReport allele calls with a bit of reformatting (remove rows without allele calls):
http://snp.gs.washington.edu/SeattleSeqAnnotation138/
It will provide the reference sequence as well as specifically provide the European HapMap frequency.
That said, be careful how you export your data. I think the default in GenomeStudio is to use the "TOP" format. However, this often won't match the genomic reference. Once you know the built for which the genome coordinates are defined (e.g. hg18 or hg19), you should use the "Plus" format to get the alleles with respect to that reference. If you can only get hg18 coordinates, you'll need to use liftOver (in Galaxy, etc.) to convert them to hg19.
•
written
by
Charles Warden
Log in to answer this question.
-
How to convert Illumina SNP ids having mappings to multiple RSID in loci to rsid file ?
written by vasank1958 •
-
How do I annotate 1000 genomes VCF with superpopulation minor allele frequency?
written by Isobel •
-
Illumina final report file to PLINK
written by f.abbondanza94 •
-
How to Convert Illumina reports to Vcf files?
written by vikrantvicky.1987 •
-
Recalling SNPs count from Illumina SNP array raw data
written by nadne •
-
Ensembl Perl API - How to get genotype frequencies
written by entheologist33 •
-
What each column from Genotyping data is? (GWAS)
written by kkp5232 •
-
Handling Illumina Omni2.5 data from multiple array versions and multiple genome version (with a vie…
written by Nicholas Blackburn
-
Convert Illumina Final Report File To Plink
written by Kyle Oliveira •
-
Problem After Imputation From 1000 Genome
written by Raha •
I wouldn't use HapMap frequencies, since they didn't use the 2.5 chip in their genotyping (according to the genotyping platform filter in HapMart)
@Katie D'Aco Then what would you use?
I think 1000 genomes frequencies would be fine to use here.
I though the HapMap frequencies were deduced from the 1000genomes project. I hoped new projects, such as kaviar, would supersede the original 1000genomes.