Extracting SNPs Based on Position Information

264 views
Skip to first unread message

Sue

unread,
Apr 2, 2024, 6:37:28 PM4/2/24
to plink2-users

Hi Christopher,

I'm reaching out to seek guidance regarding the extraction of SNPs from VCF files. Currently, I'm dealing with VCF files that contain SNP information tagged with axiom IDs, rather than the standard rs# identifiers. Additionally, I possess a list of target SNPs for genetic risk score analysis, identified by their rs#.

Given this scenario, I'm facing the challenge of not having a reference to match axiom IDs with their corresponding rs# identifiers. Consequently, I'm interested in extracting SNPs based on the positional information of the target SNPs.

To elaborate, my objective is to extract SNPs surrounding the target SNPs using the positional information available in the VCF files. Could you please advise if there is a method to achieve this using either PLINK 1.9 or PLINK 2?

I greatly appreciate your assistance with this matter.


Christopher Chang

unread,
Apr 3, 2024, 12:47:42 PM4/3/24
to plink2-users
With PLINK 2, "--extract bed0 <filename>" extracts SNPs within the positional intervals defined in the provided UCSC BED file.
Reply all
Reply to author
Forward
0 new messages