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

From ZhangLabWiki
Jump to navigation Jump to search
Line 15: Line 15:
#Incorporate an error model for sequencing;
#Incorporate an error model for sequencing;
#Should generate one raw read file and one methylHap file.
#Should generate one raw read file and one methylHap file.
  Perl Script: [[Media:mockReadsGenV1.txt]]

Revision as of 17:17, 24 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.
  Perl Script: Media:mockReadsGenV1.txt