Daniel:Notebook/GenomeMiner/2013-8-13

From ZhangLabWiki
Revision as of 23:58, 13 August 2013 by >Djacobse
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

HL155

Back to Calendar

Continuing Matt's pipeline. Looking over the perl script.

I believe I fixed the perl script (uses my notation): File:Count mismatch.txt

Results

Error Rate: 2.08727984165367

  Total reads Filtered reads Perfect match 1 ins/del 1 sub 2 ins/del 2 sub 1 ins/del & 1 sub 3+
Number 20727275 15082383 4411040 3584377 2417322 1893421 754110 1647242 374871
Percent 100 72.77 21.28 17.29 11.66 9.13 3.64 7.95 1.81

Comparing this against Matt's results, it looks like the MYcroarray probes have more errors. Around 75% of Matt's reads were perfect, but mine tend to have at least 1 error (only around 20% perfect).

Also analyzed the results (partly). See here for analysis.