Athurva Gore/LabNotes/2009-8-18: Difference between revisions
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
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
- TDMR Set (Also check if out-of-memory happens again...maybe run on Miner?)
- 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