Hi every one, I have a list of SNPs (one millions) for the below format:
SNPs name rs7552632 chr1:109457160 rs7552714 rs7552841 rs7553178
I need to create SNPs.map file with the below format:
SNPs name Chr position Rs678678 1 1561131345
I download from dbSNPS all the SNPS location.
I need to get one line from snps list and search the reference file and retrieve the snps chr and its position. Is there is any R package or script to do it?
Some snps name written in chr#: position, so I need to separate the chromosome number from its position.
I appreciate any help.