Calculating GC content of all CpG islands
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.
• 2,107 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Please use the search function. One approach is to get the sequences of the coordinates using
bedtools getfastaand then use e.g. How To Calculate A Region‘ Gc Content In Reference Genome??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.