Question: How to convert Illumina top/bot SNPs to forward/reverse?
1
gravatar for Heuer
4.6 years ago by
Heuer10
Heuer10 wrote:

My goal is to convert Illumina data to forward/reverse strand data (AffyMetrix). I have 10K, 60K, 80K and 660K data from both sources.

I have the following types of files:

***.ref
snp CHROMOSOME POSITION COUNT_ALLELE OTHER_ALLELE
snp1 17 190170 G A
snp2 17 469495 A G
--- One line for each SNP


***.dat
pig_id snp1 snp2 snp3 snp4 --- snp58448
11641206 1 0 1 0 0 
11324561 1 0 2 0 0 
14561322 2 1 1 0 1 
13513507 0 0 2 0 0 
--One line per animal, one row per snp

I also have an Illumina HD annotation file which looks as follows:

Illumina     Inc.                                                                   
[Heading]                                                                       
Descriptor File Name    GGP HD Porcine.csv                                                                  
Assay Format    Infinium HD Ultra                                                                   
Date Manufactured   5/15/2013                                                                   
Loci Count  68516                                                                   
[Assay]                                                                     
IlmnID  Name    IlmnStrand  SNP AddressA_ID AlleleA_ProbeSeq    AddressB_ID AlleleB_ProbeSeq    GenomeBuild Chr MapInfo Ploidy  Species Source  SourceVersion   SourceStrand    SourceSeq   TopGenomicSeq   BeadSetID
LD-Porcine80K_ALGA0000022-0_T_F_2164561890-0_T_F_2165597341 ALGA0000022 TOP [A/G]   19808437 *seq*     10.2 1   865364  diploid Sus scrofa  rs80958395  0   TOP *seq*    [A/G]877

When I compare SNPs from the illumina annotation file with the affymetrix annotation file, I find inconsistency in the SNP callings.
For example after manually checking some values I find the following:
sometimes illumina TOP A/G is called as T/C in affymetrix, it sometimes is also called as A/G in affymetrix.
sometimes illumina TOP A/C is (sometimes?)called as T/G in affymetrix,
sometimes illumina BOT T/G is (sometimes?)called as A/C in affymetrix

I hope somebody will be able to help me out here

affymetrix convert snp illumina • 4.8k views
ADD COMMENTlink modified 3.9 years ago by Floris Brenk970 • written 4.6 years ago by Heuer10
0
gravatar for Floris Brenk
3.9 years ago by
Floris Brenk970
USA
Floris Brenk970 wrote:

Check this out:

http://gengen.openbioinformatics.org/en/latest/tutorial/coding/

ADD COMMENTlink written 3.9 years ago by Floris Brenk970
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1642 users visited in the last hour