Kun:LabNotes/CpgSeq/2008-5-23: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
(New page: ==Informatics== ===Read mapping=== Should generate three files: #Haplotype file: xxx.methylHap.txt *One molecule per line; *target_id, offset1, methylotype1, offset2, methylotype2... #Meth...)
 
Line 3: Line 3:
Should generate three files:
Should generate three files:
#Haplotype file: xxx.methylHap.txt
#Haplotype file: xxx.methylHap.txt
*One molecule per line;
##One molecule per line;
*target_id, offset1, methylotype1, offset2, methylotype2...
##target_id, offset1, methylotype1, offset2, methylotype2...
#Methylation level file, average methylation level per site: xxx.methylFreq.txt
#Methylation level file, average methylation level per site: xxx.methylFreq.txt
#Methylation LD statistics file, all LD statistics: xxx.methylLD.txt
#Methylation LD statistics file, all LD statistics: xxx.methylLD.txt
 
===Mock read generation===
===Mock read generation===
#50% methylation at every CpG site;
#50% methylation at every CpG site;

Revision as of 21:59, 23 May 2008

Informatics

Read mapping

Should generate three files:

  1. Haplotype file: xxx.methylHap.txt
    1. One molecule per line;
    2. target_id, offset1, methylotype1, offset2, methylotype2...
  2. Methylation level file, average methylation level per site: xxx.methylFreq.txt
  3. Methylation LD statistics file, all LD statistics: xxx.methylLD.txt

Mock read generation

  1. 50% methylation at every CpG site;
  2. No correlation between adjacent CpG;
  3. Read position is evenly distributed;
  4. Can generate both single reads and pair-end reads;
  5. Incorporate an error model for sequencing;
  6. Should generate one raw read file and one methylHap file.