Having been really pleased with results from the Ray assembler I was a little disappointed to see that there are no coverage statistics in the fasta headers (as there would be in Velvet).
What would be the best way to assess the coverage of each contig assembled?
and for bonus points; what would be the easiest way to incorporate the coverage into the fasta headers?
currently looking into mapping reads against assemblies using bwa/samtools etc.
[Working with HiSeq 101bp PE data, 2 lanes. Not using Velvet because I literally don't have the computing resources to run it on all the samples]