I am interested in accessing SNP annotation from dbSNP using a package like SNPlocs.Hsapiens.dbSNP-*. It seems however that these packages are not updated for every dbSNP build.
- Why are these pkgs not built more often?
- If I wanted to build a SNPlocs.Hsapiens.dbSNP-* package for the most recent build of the NCBI SNP database, what would I need to do? Would it suffice to run those tools indicated in the package's SNPlocs.Hsapiens.dbSNP***.GRCh38/inst/tools/README.TXT?
- If I do build a new SNPlocs.Hsapiens.dbSNP-* pkg, does it still work nicely with BSgenome, namely, could I still inject SNPs and have them landing at the correction locations?
Thank you very much in advance.
Thank you for your quick reply.
May I ask another thing? For each dbSNP build, NCBI provides this file RsMergeArch.bcp.gz that contains a translation table for SNPs that have been merged to new rsIDs. Is this type of resolution taken into account in SNPlocs.Hsapiens.dbSNP149.GRCh38, or do I need to translate my old SNP IDs before passing them to SNPlocs' functions?
Hi,
I didn't know about RsMergeArch.bcp.gz. Sounds like a valuable resource. Right now this information is not included in the SNPlocs packages and the
snpsById()
extractor doesn't perform any translation of the supplied ids.H.
Hi Hervé,
Thanks anyways. I am going to install R devel and your recent SNPlocs package.
RM