Tag: bioinformatics

As a follow up to this previous post, I obtained the RefSeq gene models that have a CpG island within the 1,000 bp upstream region. There were 579 / 1009 GO terms. Note previously I identified ~780 RefSeq gene IDs, however curiously using the -wb and -f 1 options together with intersectBed, I got more…

Downloaded RefSeq gene models from UCSC Genome Browser. Total of 34,565 unique RefSeq ID Using this script refseq2go.pl join refseq to gene ontology (GO) terms 17,458 / 34,565 have GO terms, about half. In total there were 169,087 GO terms (10,316 unique) for the 17,458 RefSeq: 50,868 Component 53,144 Function 65,075 Process Top 10 Components…