Ns126:Calendar/NOTES/2015-8-31: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Shicheng
No edit summary
>Shicheng
No edit summary
Line 15: Line 15:
*[[Samples List]]
*[[Samples List]]


[[File:GSI.22.tissues.png|400px]]
*[[File:GSI.22.tissues.png|1600px]]
[[File:Code.20190903.TSI.22.tissues.R.txt]]
[[File:Code.20190903.TSI.22.tissues.R.txt]]



Revision as of 04:21, 23 September 2015


  • liftover epigenetic annotation from hg18 to hg19
./liftOver Hic.topological.domain.hESC.hg18.bed hg18ToHg19.over.chain Hic.topological.domain.hESC.hg19.bed tmp
./liftOver Hic.topological.domain.IMR90.hg18.bed  hg18ToHg19.over.chain Hic.topological.domain.IMR90.hg19.bed tmp
./liftOver Hic.boundary.IMR90.hg18.bed hg18ToHg19.over.chain Hic.boundary.IMR90.hg19.bed tmp
./liftOver Hic.boundary.hESC.hg18.bed hg18ToHg19.over.chain Hic.boundary.hESC.hg19.bed tmp
./liftOver Hic.common.boundary.hESC.IMR90.hg18.bed hg18ToHg19.over.chain Hic.common.boundary.hESC.IMR90.hg19.bed tmp

File:Code.20190903.TSI.22.tissues.R.txt



compare the cluster analysis with raw methylation signal

achieve raw methylFreq files

I found the overlapped CpG sites between methylation haplotype and raw methylFreq (depth>5) was only 49. The reason when we calculate the methylation haplotype we did not get rid of low coverage reads. Therefore, I need collect the raw methylation methylFreq file and merage them again without low coverage reads discarding. After discuss with Dinh, the methylFreq files of 106 samples were collected (10 N37, 36 Salk, 57 MONOD and 3 Heyn2013Age)