match CpG positions to genes

Hello, I want to map the CpGs I have to their correspondant genes. I have the following information:

chr     CpG_position     sample1        sample2    ...
1            1058335       .....          ....
1            1058345       .....          ....
3              ....        .....          ....

How can I map those cpgs to a gene name if I only know that hg37 has been used? They are TAB seq and BS seq results (EPIC)

meth-seq cpg epigenomics rna-seq

• 2,041 views • link

updated 8.4 years ago by Devon Ryan

• written 8.4 years ago by Pin.Bioinf

1 answer

bedtools intersect or bedtools annotate would come to mind. Download or make a BED file with the bounds of genes to use with that.

• 1 views • link

written 8.4 years ago by Devon Ryan

More posts like this

Time-series RNA-seq analysis using DESeq2
written 3.0 years ago by Debebe •

Hello All, I am working on time-series RNA-seq data. I have four individuals and collected RNA-seq samples for each individual at 4 different time points …
t-SNE plot dots that are close together assigned to different clusters.
written 3.0 years ago by leranwangcs •

Hi, I have been having trouble to understand why some dots are so close together in t-SNE plot but they are assigned to different clusters …
R removes 1st column (gene-id) from featureCounts count.txt table
written 3.4 years ago by Pegasus •

Hi all, I generated a count.txt for sorted.bam files using featureCounts on Linux following the RNA-SEQ data analysis steps. 1- Using txt.editor, I checked the …
Metagenome binning of samples with individual assembly.
written 3.5 years ago by jylee •

Hello, I have been working with human fecal metagenome samples. The samples were individually assembled using metaSPAdes, resulting in the following files: `sample1.fastq` to `sample1.contigs.fasta` …
creating a txt file with absolute path to each directory
written 3.9 years ago by arshad1292 •

I have 3200 files in a directory. I want to generate a single .txt file containing file name one column and absolute path to each …
Interpreting ATAC-SEQ motif enrichment output
written 5.4 years ago by emalekos •

After aligning and peak calling ATAC SEQ data, I used the homer getDifferentialPeaks script to do pairwise comparisons between samples. This produces a file that …
Too low recall for insertion
written 5.7 years ago by 275237370 •

Dear developer. I am working on measuring the sv calling performance of vg software using nanopore data (80× coverage) and Illumina Platinum Genome data (50× …
Differential methylation with 1 sample only
written 8.4 years ago by Pin.Bioinf

Hello, I have a smoothed, filtered methylation matrix that has : - 8 samples for TAB sequencing data - each sample is a different time …
Coverage of methylation txt files
written 8.5 years ago by Pin.Bioinf

Hello, I have been provided with methylation txt files for target_regions that look like this: chr position meth unmeth total degree.of.meth chr1 1058346 1 1 …
Lux oxBS-seq and BS-seq quantificator of 5mc and 5hmC. Which are the input files?
written 8.5 years ago by Pin.Bioinf

Ive been reading the Lux documentation, which is a quantifier of 5mc and 5hmc and they say the input files for this program are the …