Changing from "rs" to "chr_pos_ref_alt"
2
0
Entering edit mode
21 months ago
Karan ▴ 10

I'm attempting to change my SNP ids (in my bim file) from "rs_________" to "chr_pos_ref_alt" but am having some trouble.

Can someone provide some guidance? What awk scripts could I use?

EDIT: basically trying to switch from "rs11260745" to "1_24088128_G_A"

plink • 1.1k views
ADD COMMENT
2
Entering edit mode
21 months ago

plink 2.0's --set-all-var-ids flag provides one way to do this.

ADD COMMENT
0
Entering edit mode

awesome, thank you!

ADD REPLY
0
Entering edit mode

That's a very cool feature. Can I convert from CHR_POS to rsIDs by PLINK, vice versa?

ADD REPLY
1
Entering edit mode
21 months ago
Ram 43k

Use bcftools annotate -I. Your requirement is the exact example given in the manual.

ADD COMMENT
0
Entering edit mode

perfect, thanks for linking the manual!

ADD REPLY

Login before adding your answer.

Traffic: 1969 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