[seb567]

Quick Question, does Ray support hybrid assembly?

Yes.

Ray: simultaneous assembly of reads from a mix of high-throughput sequencing technologies.
Sébastien Boisvert, François Laviolette, and Jacques Corbeil.
Journal of Computational Biology (Mary Ann Liebert, Inc. publishers).
November 2010, 17(11): 1519-1533.
doi:10.1089/cmb.2009.0238

See others here too: http://denovoassembler.sourceforge.n...lications.html

454 and illumina?

Yes. See paper above.

Can ray do a reference assembly?

No. Ray only does de novo stuff.

Is there GUI?

Not in the vanilla, but I know that Applied Maths Inc. is working on one.

See http://www.applied-maths.com/opensource/opensource.htm

I almost got to using RAY but had huge problems with openmpi.

Adrian

Open-MPI is very easy to install. If you have a decent GNU/Linux distribution, then you can install it right away as a package.

Otherwise:

wget http://www.open-mpi.org/software/omp...-1.4.3.tar.bz2
tar xjf openmpi-1.4.3.tar.bz2
cd openmpi-1.4.3
./configure --prefix=$(pwd)/build
make
make install

You now have Open-MPI in $(pwd)/build

In particular, $(pwd)/build/bin/mpic++ and $(pwd)/build/bin/mpirun.

Now add $(pwd)/build/bin to your path.

For example in your .bashrc:

export PATH=/home/boiseb01/software/openmpi-1.4.1/build/bin:$PATH


Anyway, you should use a decent GNU/Linux distribution -- it makes things easier.

[AdrianP]

Quick Question, does Ray support hybrid assembly? 454 and illumina? Can ray do a reference assembly? Is there GUI?

I almost got to using RAY but had huge problems with openmpi.

Adrian

[seb567]

Greetings to you all,

After making sure I know how to use velvet I will also try Ray and SOAP.

Hi !

I am the author of Ray so if you have any question, ask away.

Basically, with Ray, you will need to convert your two files to fasta or fastq format.

There is a script in maq for that.
http://maq.sourceforge.net/

maq-0.7.1/scripts/fq_all2std.pl export2std 100611_s_4_1_seq_GDR-7.txt > 100611_s_4_1_seq_GDR-7.txt.fastq
maq-0.7.1/scripts/fq_all2std.pl export2std 100611_s_4_2_seq_GDR-7.txt > 100611_s_4_2_seq_GDR-7.txt.fastq


Ray is available at http://denovoassembler.sf.net

Then, using Ray, you assemble these reads:

mpirun -np 8 Ray -k 31 -p 100611_s_4_1_seq_GDR-7.txt.fastq 100611_s_4_2_seq_GR-7.txt.fastq -o Ray-test-1.4.0

I encourage you to explore the files written by Ray:

ls Ray-test-1.4.0.*

see http://denovoassembler.sourceforge.n...00000000000000

[AdrianP]

Yeah actually my next genome to work with is a mitGenome and is about 70k, pretty cool.

I will start working with consed, not an easy program to work with but as I understand incredibly useful.

[tonybolger]

If Velvet assembles poor (small) contigs when other programs with same settings (coverage and insertsize) do much much better, what can be my conclusions?

Could be many things: appropriate vs inappropriate settings, more/less sensitive to data error vs coverage, or simply wrong/right tool for this particular job.

You also could have one tool making a decent size but completely incorrect assembly, with another making a cautious but correct assembly. N50/size isn't everything. Can you validate your results somehow, e.g. against another closely-related known genome?

Either way, getting the best out of a dataset may require months of trial and error, tweaking etc. Even getting a particular tool do run properly and produce decent output can take weeks and can make a massive difference - the DBG assemblers all seem to have glass jaws. Pre-filtering the data seems to make a massive difference to most though.

Incidentally, I don't have a lot of experience with velvet in particular, since it's simply too heavy for my project (>1GBase genome)

[Thorondor]

yes, it is the first hit when you google "definition: N50"

separate to "shortpaired" I assume. So you can use 2 different PE libraries with different insert sizes, but maybe I am wrong.

Total is their calculated base pairs total, but keep in mind that the bp length of eacht transcript is the "real" bp length minus kmer plus 1, as mentioned somewhere.^^

max might be the longest contig, but I don't remember it exactly, since you need to do more statistics anyway. ;-)

[AdrianP]

The Manual I have read it a few times but I guess what I was asking is what does that mean "separate insert size library" ? Separate from what?

As for the research paper, i had a look at it before, I can't find a defenition for n50, it jumps straightly to using that term, eve in the abstract. I would not have asked if I did not do the research myself.

I found the answer here after googling n50 http://seqanswers.com/forums/showthread.php?t=2332

[Thorondor]

come on, do a bit more research on your own. :P

read the velvet paper:

http://genome.cshlp.org/content/18/5/821.short

You don't want to use the graph and if you know what a graph is you should also know what nodes are. Nodes are the vertices of a graph.
shortpaired2 ist the same as shortpaired but for a separate insert size library (also stated in the manual).

[AdrianP]

I was wondering a bit more about Velvet's last line output.

What is n50? (I understand that it is a measurement of quality, the higher the better???)
What is max?
What is total?
What are nodes? (this isn't as important since i believe it is related to the graph that velvet builds, and I do not use the graph, just the final contigs.fa file. Should I use the graph?)

Also, what is the diffrence between shortpaired and shortpaired2 ? Something to do with inert libraries...

Thanks a lot!

[Thorondor]

i have no clue what genegenious is and on what algorithm it is based. So if it is a ovelap-based method, yes it is possible and depends on kmer, amount of reads you have, expected coverage, read length .....

[AdrianP]

One thing that is puzzeling me is that Genegenious takes a few days to assemble the data that velvet assembles in 10-15 mins. Is it normal that velvet runs so fast?

[Thorondor]

http://brianknaus.com/software/srtoolbox/fastq2fasta.pl

first hit in google. ;-) Also normally trimming "programs" takes fastq as input and output a fasta.

velvet needs a good coverage to do well because it's de brujin based. Since I don't know on what data you run velvet I and what you expect there is no help. Try smaller kmers, try different parameters, try multiple kmer....

[AdrianP]

If Velvet assembles poor (small) contigs when other programs with same settings (coverage and insertsize) do much much better, what can be my conclusions?

By the way most of those scripts are for fasta, and i got fastq, is there a script that converts?

Adrian

[Thorondor]

sure there are a lot of ways. ;-) Use a script, you even have the (length + kmer -1) in the id of the contig so it is really easy.

here are some perl scripts that might help:

http://wiki.bioinformatics.ucdavis.edu/index.php/Data_Analysis

[AdrianP]

Thank you for your previous reply.
Is there any way to see which contig is largest? Somehow sort contigs by their size?