SNP integration: Feature and implementation details
Implemented due to Neha's suggestion, test feature and whether is functions flexibly enough.
Current implementation:
- Input: list of SNP positions and SNP genome assembly (rsid not enough) and GRN
- Overlap peaks with SNP positions and add (1) number of overlapping SNPs and (2) SNP annotation for each peak
- The returned modified GRN data frame then contains additional columns that can be used for filtering etc
Would it also be helpful to do the same with genes?