Daniel:Notebook/GenomeMiner/2013-9-9

From ZhangLabWiki
Revision as of 23:11, 9 September 2013 by >Djacobse (→‎Mock HL155)
Jump to navigation Jump to search

Mock HL155

At Dr. Zhang's recommendation Matt and I are working on mock sequencing data to investigate some of the discrepancies from sequencing, specifically why the reverse complement had such a large effect on the results from the MYcroarray data.

Generating the Mock Sequencing Data

The first run is going to be perfect copies of the sequencing data, since that is the easiest to make.

Sequences are taken from the ones Matt used for analysis on 9/3/2013. Those are taken from the original sequences from Dr. Zhang. They have been changed to the reverse complement and chopped down to 50bp from the 3' end.

Base qualities were extracted from the real sequencing data, specifically sam format (NAME OF SAMFILE). Being in sam format, they use the "Illumina 1.8+ Phred+33, raw reads typically (0, 41)" format [1].