User Tools

Site Tools


lecture_notes:04-16-2010

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:04-16-2010 [2010/04/19 18:41]
jlong Mapping code, .rdb code for Makefile, 454NewblerMetrics
lecture_notes:04-16-2010 [2010/04/20 03:09] (current)
karplus fixed punctuation that was messing up final list
Line 13: Line 13:
   * Currently, tools are installed under /​campusdata/​BME235/​bin/​old_Newbler/​.   * Currently, tools are installed under /​campusdata/​BME235/​bin/​old_Newbler/​.
   * Tools with prefix "​gs"​ are not supposed to be run directly.   * Tools with prefix "​gs"​ are not supposed to be run directly.
-  * Kevin has written several scripts in Python (version 2.6) which aid in building and analyzing genomes. Currently, these scripts do not work on Campusrocks,​ as the version of Python installed is 2.4 and it is under the process of being updated to version 2.6.+  * Kevin has written several scripts in Python (version 2.6) which aid in building and analyzing genomes. Currently, these scripts do not work on Campusrocks,​ as the version of Python installed is 2.4 and it is under the process of being updated to version 2.6. (Python2.6.5 has now been installed in /​campusdata/​BME235/​bin/ ​ --- //​[[karplus@soe.ucsc.edu|Kevin Karplus]] 2010/04/19 20:03//)
   * Newbler assembly tools take .sff (color space and quality data) files as input and converts them into .fna (fasta file with nucleotide information) files.   * Newbler assembly tools take .sff (color space and quality data) files as input and converts them into .fna (fasta file with nucleotide information) files.
   * Good only with 454 data, and is not good on reads with length < 50.    * Good only with 454 data, and is not good on reads with length < 50. 
Line 36: Line 36:
   * Where -e 25 is specific to the Vibrio sequence coverage, and -noace prevents the building of an ace file (large file) which is used with CONSED.   * Where -e 25 is specific to the Vibrio sequence coverage, and -noace prevents the building of an ace file (large file) which is used with CONSED.
   * Slug is AT-rich, so Illumina data may be better than 454.   * Slug is AT-rich, so Illumina data may be better than 454.
-  * rdb files were described as useful simple to create relational databases. ​ An example of rdb file generation with a makefile is given below as implemented by Kevin'​s in /​pluck/​rachel/​combined_cleaning1/​Makefile .+  * rdb files were described as useful simple to create relational databases. ​ An example of rdb file generation with a makefile is given below as implemented by Kevin'​s in /​pluck/​rachel/​combined_cleaning1/​Makefile ​. Note that this example was **not** given in class, and is intended for pulling out a subset of the contigs, not making an rdb file for all contigs.
 <​code>​ <​code>​
 %.stats: %.ids  %.stats: %.ids
Line 53: Line 53:
     ​     ​
 ====== Things to remember while running assembly tools ====== ====== Things to remember while running assembly tools ======
-  ​ + 
-    * All the assemblies should be listed under /​campusdata/​BME235/​assemblies.  +  * All the assemblies should be listed under /​campusdata/​BME235/​assemblies.  
-    * Include .cshrc file in your path.  +  * Include .cshrc file in your path.  
-    * Its better to run the tool in the current working directory.  +  * Its better to run the tool in the current working directory.  
-    * Create a README file in each new directory and it should contain all the necessary stuff required to run the assembly tool.  +  * Create a README file in each new directory and it should contain all the necessary stuff required to run the assembly tool.  
-    * Create Makefile for each assembly tool. (Makefile for newbler_assembly tool is in /​campusdata/​BME235/​assemblies/​Pog//​newbler-assembly1/​). You can use it as a template and modify the data source and the expected coverage as required. Makefile should be considered as "a book for lab protocols"​.  +  * Create Makefile for each assembly tool. (Makefile for newbler_assembly tool is in /​campusdata/​BME235/​assemblies/​Pog/​newbler-assembly1/​ ). You can use it as a template and modify the data source and the expected coverage as required. Makefile should be considered as "a book for lab protocols"​. 
-    Its always better to say append to make.log in Makefile.  +  It is always better to say append to make.log in Makefile.  
-    * Wiki page for assembly tools should contain a summary of how to run the tool and other things that might be useful to look at.+  * Wiki page for assembly tools should contain a summary of how to run the tool and other things that might be useful to look at.
  
  
lecture_notes/04-16-2010.1271702484.txt.gz · Last modified: 2010/04/19 18:41 by jlong