User Tools

Site Tools


lecture_notes:04-20-2015

This is an old revision of the document!


A PCRE internal error occured. This might be caused by a faulty plugin

======Team 1 report: assembly with Meraculous====== =====Basic features===== Published by the Joint Genome Institute, part of the US Department of Energy. Meraculous was initially designed for haploid assembly, but currently supports diploid assembly as well. The advantages of this assembler include multi-threaded and parallelized computation, absence of error-correction for faster processing, paired-end short reads compatibility (e.g., Illumina), efficient and conservative traversal of subgraphs of the de Bruijn graph, selection of kmer set, production of a set of maximal linear sub-paths of the de Bruijn graph, and alignment of reads to assembly in order to identify useful read-pair information and closure of gaps. Meraculous has been used to assemble the //Pichia stipitis// genome, a 15.4 Mb genome, using 75 bp paired reads with 425x coverage. The resulting assembly covered 95% of the genome and had an N50 of 101 kb. =====Meraculous algorithm===== * Counts occurrences of each kmer in the data set. * Removes kmers whose frequency are below a threshold provided by the user. * For each kmer, counts the number of high-quality single-base extensions * Classifies the 5' and 3' ends of each kmer as U, F, or X, corresponding to having zero, one, or multiple high-quality single-base extensions * Stores the extensions of kmers with a classification in a hash * Removes non-reciprocal U-U extensions between kmers (i.e. an extension where the end of one mer is marked as U but the other is marked F). * Stores the linear subgraph of U-U extensions * Selects kmers at random and extend outwards to produce contigs * Aligns all reads to contigs via BLAST * Assembles contigs into scaffolds using paired-end data * Searches unaligned reads as potential gap-closers using mate-pair data =====Meraculous limitations===== * The assembler relies on data with high quality in order to avoid error correction * Initial release did not support polyploid genome assembly due to allowing for linear subgraphs of the de Bruijn graph only * Low memory footprint =====User experience===== * Requires an array of other scripts in other languages * Most of high level scripts are written in perl * Tested the program in small dataset and obtained contigs =====Installation===== * Main issue was get all dependencies together * There was one non-standard perl mode needed * Some scripts contain error but they aren't hard to fix. =====Running Meraculous===== * Execute run_meraculous.sh scripts along with the configuration file * Configuration file contains info on where where data is and what format it comes in * It creates a timestamped folder that includes directories containing results of each step and executables to modify the run * Then you can check the errors that made a run fail and resume the run * Logs are informative =====Overall impression===== * Straightforward to figure out what went wrong just requiring a basic understanding of Perl * Handles all directory creation for you * Logs are very useful =====Error correction===== * Meraculous requires error correction and adapter removal. Trimming is unnecessary. * High error rates stop the assembler. Need to be removed. * Kmer size chosen directly affects assembly quality =====KamerGenie===== Meraculous requires an optimal kmer size for runs. KmerGenie is a program used to give optimal assembly kmer size by generating abundance histograms for many abundance histograms for many values of k. Here is a link that helped me understand KmerGenie: http://kmergenie.bx.psu.edu/. =====Musket===== =====Previous analysis===== =====Future directions=====

You could leave a comment if you were logged in.
lecture_notes/04-20-2015.1429930733.txt.gz · Last modified: 2015/04/25 02:58 by calef