Assembly merge: Difference between revisions
Jump to navigation
Jump to search
| Line 33: | Line 33: | ||
MINMATCH = 20 # nucmer default; ideally would be read length/2 to align all reads that contain 1bp error/difference | MINMATCH = 20 # nucmer default; ideally would be read length/2 to align all reads that contain 1bp error/difference | ||
# for Pa, dropping the value to 17 made a big difference | |||
MINCLUSTER = 20 # contigs start breaking if this value is too high | MINCLUSTER = 20 # contigs start breaking if this value is too high | ||
MINOVL = 5 # noticed one case with Pa where an ~50bp insert was missed (adjacent contigs were merged) | MINOVL = 5 # noticed one case with Pa where an ~50bp insert was missed (adjacent contigs were merged) | ||
# AMOScmp has this value=10; for Pa we got about twice as many contigs using 10 instead of 5 | |||
MAXTRIM = 10 | MAXTRIM = 10 | ||
MAJORITY = 50 # reduced from 70 to avoid negative gaps | MAJORITY = 50 # reduced from 70 to avoid negative gaps | ||
ALIGNWIGGLE = 2 # duplications inserted if this value is increased | ALIGNWIGGLE = 2 # AMOScmp has this value=10; duplications inserted if this value is increased | ||
CONSERR = 0.06 | CONSERR = 0.06 | ||
Revision as of 13:33, 5 April 2008
Assemblers
Denovo
Minimus
* hash-overlap overlap:
40 bp default : too large for contig assemblies
20 bp minimum overlap; minimizer window length must be >=15bp; could these values be dropped lower?
very slow on large sequences (ex Ps.fasta, Ps.plasmid.fasta) even if USE_SIMPLE_OVERLAP=1 !!! WHY???
Velvet
* overlap:
18bp usually gives fewest contigs
24bp is also ok
15bp is too low => too many short contigs
* the coverage levels reported in the contig names look sort of small (by about 3 times smaller than they should be)
* some reads extend beyond the contig ends
* there are some reads shared by multiple contigs (status:D) ; at least one of the instances occurs at a contig end (goes beyond the contig end)
Edena
* contigs don't overlap
Comparative
AMOScmp
AMOScmp-shortReads
Defaults
MINMATCH = 20 # nucmer default; ideally would be read length/2 to align all reads that contain 1bp error/difference
# for Pa, dropping the value to 17 made a big difference
MINCLUSTER = 20 # contigs start breaking if this value is too high
MINOVL = 5 # noticed one case with Pa where an ~50bp insert was missed (adjacent contigs were merged)
# AMOScmp has this value=10; for Pa we got about twice as many contigs using 10 instead of 5
MAXTRIM = 10
MAJORITY = 50 # reduced from 70 to avoid negative gaps
ALIGNWIGGLE = 2 # AMOScmp has this value=10; duplications inserted if this value is increased
CONSERR = 0.06
Cases
No reference sequence
One data set, multiple denovo assemblers
Example:
* Solexa data * edena & velvet assemblers
Solution:
* merge 2 assembly contigs * run minimus on them
Multipls data sets, one(multiple) denovo assemblers
Example:
Solexa & 454 data velvet assemblers for each set
One reference sequence
Few indels, few rearrangements
Solution:
* AMOScmp * If there are many negative gaps try to further join contigs (fastaMerge.pl $PREFIX.fasta)
Many indels, few rearrangements
Few indels, many rearrangements
Multiple reference sequences
Examples
Pseudomonas_syringae
Reference:
Name Length %GC NC_004578.1 6397126 58.40 NC_004633.1 73661 55.15 NC_004632.1 67473 56.17
Repeats:
desc #repeats min max mean stdev sum 50bp+ 991 50 7362 393.73 792.41 390192 100bp+ 429 100 7362 815.36 1060.29 349793
Solexa reads
Type #reads min max mean Solexa 6340136 32 32 32 (~31x coverage)
Assemblies:
Assembler type input-data #reads #ctgs min max mean stdev ctgs-sum #singletons AMOScmp comparative Solaxa 6340136 187 20 577929 34863.06 91692.34 6519394 698638(11%) velvet denovo Solaxa 6340136 25161 45 5057 241.83 212.61 6084887 edena denovo Solaxa 6340136 14084 100 5075 210.92 145.68 2970720 4893301(77%)
Merged assemblies(contigs&singletons):
assemblers type input-data #reads #ctgs min max mean stdev ctgs-sum comments AMOScmp-merged ? AMOScmp(contigs) 187 166 20 804024 39272.2 121124 6519189 #merged 187-166=21 negative gaps out of a total of 32 minimus(ovl20) denovo velvet(contigs) 25161 19121 45 5057 311.3 297.27 5952381 #merged 25161-19121=6040 (25%) gaps minimus(ovl15) denovo velvet(contigs) 25161 16343 45 9903 361.32 359.78 5905143 #merged 25161-16343=8818 (35%) gaps minimus(ovl40) denovo edena+velvet(contigs) 39245 23644 45 6688 257.15 232.94 6080063 #very few 40bp overlaps are found minimus(ovl20) denovo edena+velvet(contigs) 39245 18603 45 6688 322.32 311.02 5996244
Simulated 32bp exact match reads
Type #reads min max mean Sim(ulated) 6538167 32 32 32 ( 32x coverage)
Single assemblies:
Assembler type input-data #reads #ctgs min max mean stdev ctgs-sum #singletons edena-sim denovo Sim 6538167 2068 100 47881 2994.03 4857.76 6191673 198699(3%) velvet-sim denovo Sim 6538167 2207 45 56810 2820.91 5348.36 6225757 123591(2%)
454 reads
Type #reads min max mean 454 77466 35 371 240
Pseudomonas aeruginosa b1 (PAb1)
References:
Name Length %GC PA14 6537648 66.29 PACS2 6492423 66.33 PAO1 6264404 66.56 ...
Solexa reads
Type #reads min max mean Solexa 8627900 33 33 33 (~43X coverage)
Assemblies:
All contigs:
Assembler type #ctgs min max mean stdev ctgs-sum #singletons comments AMOSCmp-PA14 comparative 2053 17 170485 3011.84 11917.53 6183320 1127399 AMOSCmp-PAO1 comparative 2797 17 75626 2161.19 5812.2 6044851 1592525 AMOScmp-PA2192 comparative 5816 17 133129 1072.8 3725.22 6239454 1601299 largest assembly velvet denovo 10684 45 16239 640.34 825.24 6841458 1241079 much better than Ps !!! edena denovo 11180 100 11300 552.36 610.52 6175460 3955865 (46%) much better than Ps !!! ssake denovo 185030 34 5490 77.21 141.23 14287079 3056893
200bp+ contigs:
Assembler type #ctgs min max mean stdev ctgs-sum AMOSCmp-PA14 comparative 428 203 170485 14262.09 22852.74 6104175 AMOSCmp-PAO1 comparative 865 200 75626 6893.96 8766.63 5963278 AMOSCmp-PA2192 comparative 1299 200 133129 4683.46 6735.52 6083817 velvet denovo 7382 200 16239 877.05 896.35 6474426 edena denovo 8316 200 11300 692.54 651.24 5759209 ssake denovo 12532 200 5490 486 329.93 6090567
Merged assemblies(contigs&singletons):
Assembler type input-data #reads #ctgs min max mean stdev ctgs-sum comments AMOSCmp-PA14-merge ? AMOSCmp-PA14(ctgs) 2053 1931 17 170485 3201.7 12981.57 6182486 2053-1931=122 gaps closed