Updating INDEL SNP positions
1
0
Entering edit mode
6.2 years ago

For microarray SNP data, I have a file containing the positions of all the SNPs. These SNPs originally came from Affy 6.0 chips but have been imputated since then. Thus, I have data on about 4.5 million SNPs and their positions... which is much more than the Affy 6.0 chips account for.

Looking at the position, it seems like the position stem from GRCh37 (hg19) instead of the more recent GRCh38 (hg38). By following Alex Reynolds advice here (how to find SNP positions (for non-bioinformaticians)) I was able to download the update SNP positions. Updating the SNP positions based on the rs# should be fairly straightforward using R.

The problem is, my data contains INDELs (e.g. chr1:730049:I at position 730049) and I am not really sure how to handle those positions. As these aren't coded with a rs#, I am a bit perplex about what to do with them.

Any advice would be appreciated!

SNP • 988 views
ADD COMMENT
1
Entering edit mode
6.2 years ago

Hello,

converting positions from one assembly to another is called liftOver. Depending on the data format you can use one of these online tools to convert the positions:

fin swimmer

ADD COMMENT

Login before adding your answer.

Traffic: 813 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6