User Tools

Site Tools


lecture_notes:05-27-2011

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
lecture_notes:05-27-2011 [2011/05/27 20:25]
karplus mapping 454_run3 not needed.
lecture_notes:05-27-2011 [2015/09/14 11:40] (current)
68.180.230.228 ↷ Links adapted because of a move operation
Line 1: Line 1:
-====== To Do ======+====== To Do (lecture notes) ​======
 We went over some things that need updating in the wiki and plans for this weekend. We went over some things that need updating in the wiki and plans for this weekend.
  
-  * Document scripts added in [[computer_resources:​bin|bin]] folder.+  * Document scripts added in [[archive:computer_resources:​bin|bin]] folder.
  
  
Line 12: Line 12:
      * Only on Illumina data.      * Only on Illumina data.
   * try 2:   * try 2:
-     * Rank 1: All Illumina data for contig building+     * Rank 1: All Illumina data as rank 1
-     * Rank 2: 454 data for scaffolding.+     * Rank 2: 454 data as rank 2 (both for contig building and scaffolding).
   * try 3:   * try 3:
      * Illumina + 454 data as rank 1.      * Illumina + 454 data as rank 1.
-  * If insert length is negative, don't treat them as PE reads. +  * If insert length is negative, don't treat them as PE reads.   
-  * Newbler assemblies.+    * This was an error in assumptions about what SOAPdenovo wants. ​ The number it wants is the total fragment length, which is what we are already estimating. ​ We just need to look at the average length for the pairs (based on the histogram) after SeqPrep. 
 +    * If Quake changes the distribution of reads (by trimming and discarding uncorrectable reads) it may be important to remap the new set of pairs to the 454 reads to get an improved estimate of fragment length
 +  * Newbler assemblies—no need for a new one, as there is no new 454 data.
   * Take another look at Barcode of Life mapping: [[http://​www.ncbi.nlm.nih.gov/​Taxonomy/​Utils/​wprintgc.cgi#​SG5|Invertebrate mitochondrial translation table]]   * Take another look at Barcode of Life mapping: [[http://​www.ncbi.nlm.nih.gov/​Taxonomy/​Utils/​wprintgc.cgi#​SG5|Invertebrate mitochondrial translation table]]
 +    * It turns out that blastall does not seem to have any documented way to tell tblastx to use a different genetic code, though the NCBI web server has the option.
 +    * The work Kevin did so far is now in assemblies/​slug/​barcode-of-life.
 +    * Next step is to use BWA to find all the SeqPrep+Quake treated Illumina reads that map to what has been found so far of the barcode.
 +
  
   * Ed presenting paper on Wednesday. ​ Get paper over the weekend.   * Ed presenting paper on Wednesday. ​ Get paper over the weekend.
lecture_notes/05-27-2011.1306553117.txt.gz · Last modified: 2011/05/27 20:25 by karplus