On Thursday 07 April 2011 04:06:55 Rob Good wrote: > I've just used SOAPdenovo to generate an initial denovo assembly of a > Drosophila melanogaster paired end Illumina run. > This has generated 60,000 contigs from 32 to 40000 bp. Oh dear ... this is going to be such a mess. > Can I : > 1. Use MIRA to map these to a reference assembly MIRA has a size cap of 15 to 20k per "read", so you will need to split longer contigs. > 2. Use MIRA to joins these contains together with some 454 reads Same here: 15 to 20kb max as sequence size. Apart from that: I think it should work. B.