Kun:LabNotes/ExonomeSeq/2009-11-15

From ZhangLabWiki
Revision as of 21:42, 15 November 2009 by KunZhang (talk | contribs) (New page: ==Analysis of iPS exome sequencing data== *Data used: **Libraries were constructed with the 1-adaptor protocol. **DF6/Foreskin: HL022, two lanes per line, captured with probe set #1-7; **C...)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Analysis of iPS exome sequencing data

  • Data used:
    • Libraries were constructed with the 1-adaptor protocol.
    • DF6/Foreskin: HL022, two lanes per line, captured with probe set #1-7;
    • CV-iPS/CV-fibroblast:
      • HL020, two lanes per line, captured with probe set #8-9;
      • HL022, three lanes for CV-iPS, captured with probe set #1-7;
  • Read mapping and variant calling.
    • Reads that contain potential H1/H2 capturing arms are identified and removed by mapping to a database that contains H1/H2 with 20bp flanking sequences.
    • The remaining reads were mapped, down-sampled and fed to Samtools.
    • This is the current mapping script: variantCallerBowtieSam.pl).
    • Potential mutations were identified by comparing two pileup files ( pileup2variantsPair.pl);
    • To be considered as a candidate mutation, I required that no single read contains the second allele in the reference sample.
  • Results: