-
Notifications
You must be signed in to change notification settings - Fork 82
-
What depth of coverage do I need?
At a minimum, 20-30X coverage is required to get a usable result from a de novo assembly of short reads. I usually request 40X.
-
Can I assemble very short (<50bp) reads with SGA?
This is not recommended. For very short reads a de Bruijn graph assembler will give as good or better results with lower run time. SGA is currently designed for 100bp or greater reads.
-
The sga assemble/fm-merge steps are taking too much time/memory.
The cause of this problem is usually selecting an overlap length that is too short. If the overlap length is very short, many edges will be created in the graph between reads that have short repeats at their ends. The suggested parameters for 100bp reads are overlap length 65 for fm-merge and overlap length 75 for sga assemble. You should try this parameter combination first, then tweak the overlap parameter in sga assemble to try to improve the assembly.
-
I want to assemble a very large genome. How do I efficiently index the reads?
See this page which describes the construction algorithm for large data sets.
-
How do I build scaffolds using multiple libraries?
See here
-
What parameters should I tune to improve my assembly?
See here
-
Can I use bwa mem during scaffolding?
Yes, but you need to remove secondary alignments from the bam file. See this mailing list thread.
-
What is an ASQG file?
See here
-
Where can I get further help or advice?
There is a mailing list/google group for sga. You can sign up [here] (http://groups.google.com/group/sga-users)