Banana Slug Genomics

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
===== ALLPATHS =====
Attached is the Allpaths3 version 1.0 documentation converted from .docx format to .pdf. You can download this file by clicking the following link: {{:bioinformatic_tools:allpathsv3_manual_r1.0_2.pdf|allpathsv3_manual}}.

====Potential Pitfalls====
Designed to work with 100+ bp paired end reads from a *minimum of one short and one long set of libraries. Also the program expects 40X coverage from each of those libraries! Additionally they say it requires a minimum of 32Gb of ram which I assume means shared memory, so it may not work on our cluster. Maybe it would be useful to run on small portions of our data though?

====High Level Overview====
ALLPATHS is the most recent (as of this writing) tool developed by the broad institute to assemble shotgun sequences[(cite:broad>http://www.broadinstitute.org/science/programs/genome-biology/computational-rd/computational-research-and-development)]. The broad institute claims that version 3 of the program can assemble up to mammalian sized genomes if the reads are at least 100+ base pairs[(cite:broad)]. Version 3 (currently 3.2) of the program may be downloaded from [[ftp://ftp.broad.mit.edu/pub/crd/ALLPATHS/Release-3-0/|here]] and that folder also contains some documentation on how to use the program. Also the program ships with test data that you can assemble. 

In class I mentioned that a fellow student, Amie, had a lot of trouble with assembling via Euler. Actually that program was ALLPATHS, but when she used it at the beginning of last quarter it was still in version 1.0. According to Broad version 1.0 was only expected to work on test data[(cite:broad)] (although Amie couldn't even get that working). Since then they have come out with two major revisions and one minor revision, and they claim that this version works on real data, so I think it is definitely worth a shot.  --- //[[jstjohn@soe.ucsc.edu|John St. John]] 2010/04/06 08:21//

===Expected Data Quality===
On the broad site they mention that they have tested the algorithm out on very short reads (30+bp) and short reads (100+bp). They mainly target sequencing strategies utilizing 100+ reads with the argument that in the future that will be the norm[(cite:broad)].
==Very Short Reads(30+bp paired reads)==
I am not sure exactly what this means so I'll quote it directly, it seems relevant.  --- //[[jstjohn@soe.ucsc.edu|John St. John]] 2010/04/05 17:50//
>We have developed and tested a method for assembling very short (~30 base) paired reads using the ALLPATHS algorithm. This method requires high coverage from two libraries, one from fragments of size 3-4 kb, and one from shorter fragments.[(cite:broad)]
==Short Reads (~100+bp reads)==
They are mainly developing the algorithm for slightly longer reads. According to the site they are targeting a combination approach consisting of 45x coverage of 100bp reads from 180bp fragments, 45x coverage of 100bp reads from 3000bp fragments, and "additional sequence from longer fragments for large genomes"[(cite:broad)].

Perhaps our data will not behave nicely with this algorithm? Also they didn't say whether or not they combine the reads from those three sources, or somehow analyze them separately and merge them at some point. Perhaps reading more into their documentation or publications will answer some of these questions. --- //[[jstjohn@soe.ucsc.edu|John St. John]] 2010/04/05 17:50//

===== Installation =====
Installing into ~/programs/allpaths

  Configure error, requires boost with at least the Boost.System binaries installed. Now I am installing that...
  
  Boost successfully installed. Now installing allpaths. Configure successful and currently building. --- //[[jstjohn@soe.ucsc.edu|John St. John]] 2010/04/07 00:53//


Build unsuccessful:

./ParallelVecUtilities.h: In function 'void ParallelSort(vec<T>&)':\\
./ParallelVecUtilities.h:27: error: '__gnu_parallel' has not been declared\\
./ParallelVecUtilities.h: In function 'void ParallelSort(vec<T>&, StrictWeakOrdering)':\\
./ParallelVecUtilities.h:35: error: '__gnu_parallel' has not been declared\\
./ParallelVecUtilities.h: In function 'void ParallelReverseSort(vec<T>&)':\\
./ParallelVecUtilities.h:42: error: '__gnu_parallel' has not been declared\\
./ParallelVecUtilities.h: In function 'void ParallelReverseSort(vec<T>&, StrictWeakOrdering)':\\
./ParallelVecUtilities.h:50: error: '__gnu_parallel' has not been declared\\
./ParallelVecUtilities.h: In function 'void ParallelWhatPermutation(const V&, vec<T3>&, C, bool)':\\
./ParallelVecUtilities.h:316: error: '__gnu_parallel' has not been declared

Reading deeper into the documentation (the PDF attached to this page), I see that it requires gcc-4.3+. Campusrocks currently has gcc-4.1 installed. Perhaps if we compile the latest gcc we can install this program?

===== References =====
<refnotes>notes-separator: none</refnotes>
~~REFNOTES cite~~
Banana Slug Genomics

User Tools

Site Tools

Page Tools
