Calculating GC content of all CpG islands
0
1
Entering edit mode
23 months ago

Hi Biostars, I have two files containing the genomic ranges of all CpG islands and promoter sequences of Human genome hg19. I would like to do an association test between CpG islands and promoters. So I need to calculate the GC contents of the two files. However, I am not conversant working with genomic coordinates. Could somebody please help with a command that I can use to calculate the GC contents of the two files. Thanks.

Genomic coordinates alignment gene GC content • 593 views
0
Entering edit mode

Please use the search function. One approach is to get the sequences of the coordinates using bedtools getfasta and then use e.g. How To Calculate A Region‘ Gc Content In Reference Genome？？

0
Entering edit mode

Hi

If you are using R; can try seqinr package. On that you could find function "GC" and can do whatever you are looking for.