Athurva Gore/LabNotes/2009-8-18: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Ajgore
>Ajgore
No edit summary
Line 17: Line 17:
** TDMR Set (Also check if out-of-memory happens again...maybe run on Miner?)
** TDMR Set (Also check if out-of-memory happens again...maybe run on Miner?)
*** Look into this tomorrow
*** Look into this tomorrow
* Check numbering of new probes
** Appears to be slightly inaccurate as well.
** For residue coordinate system:
*** End position should be increased by 1
*** Start position should be increased by 2
** For space coordinate system:
*** End position should be increased by 1
*** Start position should be increased by 1


=Exome SNPs - NA12878=
=Exome SNPs - NA12878=

Revision as of 18:29, 18 August 2009

Navigation

Probe Generation

  • DONE
    • KKESH72
    • RPLCACRD
    • CpG-SNP set (with 80 bp gap)
    • FlyDup936
    • A-to-I Probes for Erez and Billy
    • LeeCancer (on opposite strand to target cDNA)
    • LeeXGenes (on opposite strand to target cDNA)
    • ZhangSNP (on opposite strand to target cDNA)
  • CURRENTLY RUNNING:
    • FlyDup901
    • FlyDup5279
  • TO RERUN:
    • TDMR Set (Also check if out-of-memory happens again...maybe run on Miner?)
      • Look into this tomorrow
  • Check numbering of new probes
    • Appears to be slightly inaccurate as well.
    • For residue coordinate system:
      • End position should be increased by 1
      • Start position should be increased by 2
    • For space coordinate system:
      • End position should be increased by 1
      • Start position should be increased by 1

Exome SNPs - NA12878

  • Currently:
    • Running script to separate our SNP calls and 1KG SNP calls by in-range and out-of-range
    • NOTE: File Dr. Zhang provided uses a slightly inaccurate start position; needs to be incremented by 1 for the residue based coordinate system. End is accurate.
    • Check new probe design scripts; is this error still prevalent?
  • For our SNP calls:
    • FOUND THAT 3190/20762 were NONSPECIFIC
    • ~85% specific capture
  • Want to separate false positives from true positives in our SNPs.
  • NEXT STEP:
    • Look at IPS data; plot quality and coverage of SNPs.
    • Run bowtie on runs from NA12878 that have high false positives
    • See if there is a pattern; can we separate these false positives out somehow?

NA12878 High False Positive Reads

  • Locations
    • 07282009_HL002, Lanes 4-7
    • 08052009_HL003, Lanes 4-5
  • Storing data on genome-miner:/media/TmpStore1/AG_Scratch/Exome/NA12878
  • Using this sample for validation
    • First step, need to map these reads.
    • Start out with reads on a per experiment basis
      • Merge s_4 through s_7 in HL002, merge s_4 and s_5 in HL003
      • Run Bowtie-SAM script on these merged sets; check results
      • Currently running Bowtie-SAM on HL003, prepping HL002