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
ADD COMMENT
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??

ADD REPLY
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.

ADD REPLY

Login before adding your answer.

Traffic: 2067 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6