Ns126:Heyn2016

From ZhangLabWiki
Revision as of 00:33, 5 March 2016 by >Shicheng (→‎Bam to Hapinfo)
Jump to navigation Jump to search

Background

  • Epigenomic analysis detects aberrant super-enhancer DNA methylation in human cancer
  • WGBS by 101 bp pair-end sequencing.
  • Phred Score: 33
  • 11 primary tumor tissue, 2 metastasis tissue and 9 normal tissues

Method and Procedure

Sample

SRA Download

SRA to Fastq

Fastq to Bam

Bam to Hapinfo

  • Bam File Directory
TSCC-1: /home/ddiep/dinh_working/Bellvitge_BRI_WGBS/methylfiles/BAMfiles
TSCC-2: /oasis/tscc/scratch/ddiep/BAMfiles
Genome-miner:  /media/LTS_60T/Dinh/WGBS_LTS33/Hg19/Estellar_Bellvitge/BAMfiles
  • Bam to hapinfo
cd /oasis/tscc/scratch/ddiep/BAMfiles
perl ~/bin/SaminfoPre4hapinfo.pl > ~/oasis/Estellar2016/SaminfoPre4hapinfo.txt
cd /home/shg047/oasis/Estellar2016/hapinfo
perl ~/bin/bam2hapInfo2PBS.pl ../SaminfoPre4hapinfo.txt
  • Merge hapinfo in different chrosome to one file by sample ID

Haploinfo to Methylation Haplotype Block (MHB)

Hapinfo to Methylation Haplotype Load

trim_glore

Alignment