Kun:LabNotes/MONOD/2014-8-12

From ZhangLabWiki
Jump to navigation Jump to search

RRBS data analysis (continued)

1. Compile a list of RRBS targets

  • Take all the RRBS data that we generated from primary tumor samples, concatenate all the methylFreq files, and generate a single BED file.
 cat 6-T-1_1.methylFreq 6-T-1_2.methylFreq 6-T-2_1.methylFreq 6-T-2_2.methylFreq 6-T-3_1.methylFreq 6-T-3_2.methylFreq 6-T-4_1.methylFreq 6-T-4_2.methylFreq 6-T-5_1.methylFreq 6-T-5_2.methylFreq 7-T-1_1.methylFreq 7-T-1_2.methylFreq 7-T-2_1.methylFreq 7-T-2_2.methylFreq 7-T-3_1.methylFreq 7-T-3_2.methylFreq 7-T-4_1.methylFreq 7-T-4_2.methylFreq 7-T-5_1.methylFreq 7-T-5_2.methylFreq CTT-FFPE-100ng_1.methylFreq CTT-FFPE-100ng_2.methylFreq CTT-FFPE-5ng_1.methylFreq CTT-FFPE-5ng_2.methylFreq CTT-frozen-100ng_1.methylFreq CTT-frozen-100ng_2.methylFreq CTT-frozen-5ng_1.methylFreq CTT-frozen-5ng_2.methylFreq PC-T-1_1.methylFreq PC-T-1_2.methylFreq PC-T-2_1.methylFreq PC-T-2_2.methylFreq PC-T-4_1.methylFreq PC-T-4_2.methylFreq PC-T-6_1.methylFreq PC-T-6_2.methylFreq PC-T-7_1.methylFreq PC-T-7_2.methylFreq |  /home/dinh/scripts/methylFreq2BED.pl  20 > /home/kunzhang/CpgMIP/Data/MONOD/MONOD_primary_tumor_RRBS_combined.BED.txt
 cat MONOD_primary_tumor_RRBS_combined.BED.txt | sort -k 1 -k 2n > MONOD_primary_tumor_RRBS_combined.sorted.BED.txt  
  • Group the sites covered into clusters:
 ./bed2Clusters.pl MONOD_primary_tumor_RRBS_combined.sorted.BED.txt > MONOD_primary_tumor_RRBS_targets.BED.txt