Brugia malayi: Difference between revisions
Jump to navigation
Jump to search
(→Data) |
|||
| Line 227: | Line 227: | ||
/fs/szattic-asmg4/brugia_malayi/Data/Frg/ | /fs/szattic-asmg4/brugia_malayi/Data/Frg/ | ||
=== Contaminant search === | === Contaminant & high copy repeat search === | ||
* nucmer -maxmatch -l 20 -c 65 | * nucmer -maxmatch "-l 20 -c 65" or "-l 12 -c 24" | ||
Sanger 454 | Sanger 454 | ||
jird(26,879 ctgs) 31,501 197,420 | jird(26,879 ctgs) 31,501 197,420 # we'd probably find more contaminated reads if we align all the reads to the whole mouse genome ?? | ||
Mt 1,507 2,634 | Mt 1,507 2,634 | ||
We 49,014 23,249 | We 49,014 23,249 | ||
UniVec ? 661,586 | |||
pUC19 134 562,107 # 99% avg identity, 99% of read length | |||
HhaI(~320bp) 16,336 69,400 | HhaI(~320bp) 16,336 69,400 | ||
mRNA(264bp) . 59,706 | mRNA(264bp) . 59,706 | ||
total ~910,000 | |||
== Assemblies == | == Assemblies == | ||
Revision as of 16:09, 25 February 2010
Articles
the copy number of the repeat in B. malayi was found to be about 30,000. The 320-base-pair Hha I repeated sequences are arranged in direct tandem arrays and comprise about 12% of the genome.
Brugia malayi Hha I-repeat family element >gi|156092|gb|M12691.1|BRPRSHA Brugia malayi Hha I-repeat family element GCGCATAAATTCATCAGCAAAATTAATAAAACTTTCAATTAATCATGATTTTAATTGAATGTAAGAATTT AAATTAAATTTAAATTCAAATTTAAATTTTTAATTTTTTAAAAATTTTAAAATTTGTTATAGTTTTCCTT CATTAGACAAGGATATTGGTTCTAATTTATCAATTTTAATTCTAATTAAGTGCCAAAACTACTAAAAAAA GCTTATTTTGAAATTAATTGACTACGTTAGCTGCATTGTACCAGTGCTGGTCGTGTATTGTGTTGTCATT TTATAGTTTAAATATTAAAATACGCTTTTGTAATTAAGTTTT
Genome Info
- 6 chromosomes: 1-5, XY ; diploit genome ~ 110M bp
- 30% GC,
- 32% coding, 15% repeats
Genome Project
Brugia malayi has a diploid genome of approximately 110 Mb, organized in 6 pairs of chromosomes (five pairs of autosomes and one pair of sex chromosomes). In addition to the nuclear genome, B. malayi has a mitochondrial genome of about 14kb, and the genome of the harbored bacterial endosymbiont Wolbachia sp (1-2Mb).
The B. malayi genome project has been completed by The Institute for Genomic Research. Whole Genome Shotgun sequencing was used to obtain more than eight-fold coverage of the genome. The complete genome was assembled into approximately 8200 scaffolds and deposited in GenBank. The accession for the WGS project is AAQA00000000 and consists of sequences AAQA01000001-AAQA01029808. File location:
/fs/szasmg3/dpuiu/Brugia_malayi/Data/Bm.fasta
ctgs min q1 q2 q3 max mean n50 sum
26,879 200 836 1005 1495 611,244 3241.17 18986 87,119,350
- TIGR Genome project (TRS strain)
Contamination
- mitochondrion finished: 13,657 bp; 24% GC
- Wolbachia endosymbiont strain TRS from Brugia malayi strain wMel 1.08 Mbp 34% GC
- complete: New England Biolabs
- progress: TIGR
- Rodent: some trace contamination; 44%GC
/fs/szasmg3/dpuiu/Brugia_malayi/Data/contam.fasta contaminants min q1 q2 q3 max mean n50 sum 2,929 200 527 675 820 8994 740.04 762 2,167,588
- pUC19c vector: 2686bp, 50.63% GC
- of some 454 libs
- what about Sanger reads???
Data
1.26M Sanger reads (original TA) : medLen=773bp; medGC=32.57% 1.26M Sanger reads (contamination free): medLen=771bp; medGC=32.36%
3.21M 454 reads (original sff) : medLen=274bp 3.29M 454 reads (linker free) : medLen=247bp; medGC=36.39%
Original Traces
- 1.26M Sanger reads & 15 Libraries:
- NCBI TA
- NCBI TA FTP
SEQ_LIB_ID INSERT_SIZE INSERT_STDEV TRACE_TYPE_CODE
1047113828118 1000 300 WGS 13500
1047113856575 1000 300 PRIMERWALK 325
1047111632737 1258 377 PRIMERWALK 3
1047111632737 1258 377 WGS 305,906
1047111540304 1415 424 WGS 51772
1047112577106 1415 424 WGS 337,789
1047111718946 3123 936 WGS 47597
1047113358719 3123 936 PRIMERWALK 173
1047113358719 3123 936 WGS 246,185
1047174912885 3123 936 TRANSPOSON 1437
1047113570927 6000 1800 WGS 3193
1047111814561 7158 2147 WGS 219,306
1047111480027 17168 5150 WGS 4087
1047111488095 17168 5150 WGS 3434
1047111495007 17168 5150 WGS 3716
1047111501919 17168 5150 WGS 3697
1047111480605 22419 6725 WGS 4638
1047111516154 22419 6725 WGS 4004
1047111523212 22419 6725 WGS 3766
1047111530126 22419 6725 WGS 5686
1047113855421 23000 6900 WGS 1
total 1,260,215
FRG file: (contaminant free)
- FRG.src : TI's
- FRG.acc: 2 ..
- DST.acc: 1260217, ... , 1260234
- Location
/fs/szasmg3/dpuiu/Brugia_malayi/Data/nucmer_seq/Bm-all.frg
DST 15
FRG 1178192
LKG 530930
seqs min q1 q2 q3 max mean n50 sum
len 1,178,192 65 645 771 850 1214 724 800 853,847,771 => 8X
gc% 1,178,192 0.00 29 32.36 35 100 32.41 33 .
Problems:
- All library insert sizes are underestimated ???
- The contaminant reads align at ~91-93% id to the contaminant ctgs while the Mt/We reads align at 99% id to Mt/We finished seq. What %id thold to use for contaminant?
BACS
8,000 BAC clones @Children's Hospital Oakland Research Institute. (!!! no NCBI TA submission)
PITT FTP data
- 3.21M 454 reads
- 3K insert flx libraries (estimated to 2K based on alignment to the existing assembly)
- 20K insert tit libraries (estimated to 28K ...)
CBCB Location:
/fs/szattic-asmg4/brugia_malayi/Data/ /fs/szattic-asmg4/brugia_malayi/Data/Sff/ # Sff files /fs/szattic-asmg4/brugia_malayi/Data/Frg/ # Frg files /fs/szattic-asmg4/brugia_malayi/Data/Seq/ # Seq files
FTP access:
lftp -u bma 136.142.191.201 pass: 6279 user: bma # empty as of --Dpuiu 12:04, 8 January 2010 (EST)
Elodie's table:
/scratch1/brugia_malayi/brugia-sequencing-summary.txt.csv
# elodie's date protocol platform type description run_name Reads Mates 1 01/17/2008 WGS Standard Full run (2/2) Mix of worms (calibration of the machine) R_2008_01_31_18_01_35_FLX10070260_adminrig_ghedintestsample 534822 0 2 07/01/2008 3Kb Standard Full single worm (pUC contamination) R_2008_08_06_13_52_29_FLX10070260_adminrig_080608_Ghedin-BrugiaLTPE1 492575 84341 3 09/11/2008 3kb Standard 4/8 wells single worm (pUC contamination) R_2008_09_19_14_17_55_FLX10070260_adminrig_091908_HATFULL-MIDrepeat_GHEDIN-LTPE1 263421 49258 4 10/01/2008 3Kb Standard Full Mix of worms (still pUC contamination) R_2008_10_14_15_06_50_FLX10070260_adminrig_101408_GHEDIN-Brugia-pool_LTPEtest 59711 5096 5 02/01/2009 WGS Standard 1/4 wells Mix of worms; regions 2 & 3 were myxoma R_2009_02_27_16_11_34_FLX10070260_adminrig_022709_GHEDIN 18025 0 6 04/06/2009 WGS Standard 1/4 wells Mix of worms; with comp. bio run R_2009_04_15_14_46_56_FLX10070260_adminrig_041509_GHEDIN_r1-WGS1_r2-LMW4_r3-pool2compbio_r4-pool3compbio 118490 0 7 05/01/2009 20Kb Titanium 7/8 wells Mix of worms R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1 631287 213524 8 10/28/2009 20Kb Titanium Full Mix of worms R_2009_10_22_15_30_12_FLX10070260_adminrig_102209_GHEDIN_Brugia20kb2 1095713 377547 9 Pending 3 Kb Titanium Full Mix of worms ? ? ? . Total 2542707 729766
- 22 Sff files:
run sffReads linker
1 R_2008_01_31_18_01_35_FLX10070260_adminrig_ghedintestsample/D_2008_01_31_18_01_35_FLX10070260_adminrig_FullAnalysis/sff/E4RA0X101.sff 272923 .
1 R_2008_01_31_18_01_35_FLX10070260_adminrig_ghedintestsample/D_2008_01_31_18_01_35_FLX10070260_adminrig_FullAnalysis/sff/E4RA0X102.sff 261899 .
2 R_2008_08_06_13_52_29_FLX10070260_adminrig_080608_Ghedin-BrugiaLTPE1/D_2009_02_12_22_12_04_j_SignalProcessing/sff/FEZH5RS01.sff 228204 flx
2 R_2008_08_06_13_52_29_FLX10070260_adminrig_080608_Ghedin-BrugiaLTPE1/D_2009_02_12_22_12_04_j_SignalProcessing/sff/FEZH5RS02.sff 264371 flx
3 R_2008_09_19_14_17_55_FLX10070260_adminrig_091908_HATFULL-MIDrepeat_GHEDIN-LTPE1/FHAVB5T02.sff 86862 flx
3 R_2008_09_19_14_17_55_FLX10070260_adminrig_091908_HATFULL-MIDrepeat_GHEDIN-LTPE1/FHAVB5T03.sff 87488 flx
3 R_2008_09_19_14_17_55_FLX10070260_adminrig_091908_HATFULL-MIDrepeat_GHEDIN-LTPE1/FHAVB5T04.sff 89071 flx
4 R_2008_10_14_15_06_50_FLX10070260_adminrig_101408_GHEDIN-Brugia-pool_LTPEtest/FIOXLOM01.sff 13695 flx
4 R_2008_10_14_15_06_50_FLX10070260_adminrig_101408_GHEDIN-Brugia-pool_LTPEtest/FIOXLOM02.sff 14197 flx
4 R_2008_10_14_15_06_50_FLX10070260_adminrig_101408_GHEDIN-Brugia-pool_LTPEtest/FIOXLOM03.sff 15515 flx
4 R_2008_10_14_15_06_50_FLX10070260_adminrig_101408_GHEDIN-Brugia-pool_LTPEtest/FIOXLOM04.sff 16304 flx
5 R_2009_02_27_16_11_34_FLX10070260_adminrig_022709_GHEDIN/FRLDXKV01.sff 18025 .
6 R_2009_04_15_14_46_56_FLX10070260_adminrig_041509_GHEDIN_r1-WGS1_r2-LMW4_r3-pool2compbio_r4-pool3compbio/D_2009_04_16_14_19_21_morty_fullProcessing/FT9KOI001.sff 118490 .
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY01.sff 73807 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY02.sff 91698 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY03.sff 93878 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY04.sff 90232 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY05.sff 97065 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY06.sff 94326 tit
7 R_2009_06_05_16_18_41_FLX10070260_adminrig_060509_GHEDIN_Brugia-gDNA-TI20kb1/D_2009_06_08_15_32_36_compute-0-2_fullProcessing/sff/FW1OXFY07.sff 90281 tit
8 R_2009_10_22_15_30_12_FLX10070260_adminrig_102209_GHEDIN_Brugia20kb2/F4H5CMB01.sff 551263 tit
8 R_2009_10_22_15_30_12_FLX10070260_adminrig_102209_GHEDIN_Brugia20kb2/F4H5CMB02.sff 544450 tit
total 3,214,044 .
- 22 Frg Libraries
. lib meanIns(orig) meanIns(est) #reads #mates linker medLen medGC
1 E4RA0X101 0 0 271066 0 . 250 37.04
1 E4RA0X102 0 0 260166 0 . 249 37.06
2 FEZH5RS01 3000 2000 181035 18676 flx 228 37.89
2 FEZH5RS02 3000 2000 211064 22270 flx 227 37.79
3 FHAVB5T02 3000 2000 68708 7850 flx 244 38.35
3 FHAVB5T03 3000 2000 69306 8227 flx 243 37.98
3 FHAVB5T04 3000 2000 70028 8353 flx 243 38.10
4 FIOXLOM01 3000 0 10921 0 flx 102 43.93 # no mates , shorter read length, highest GC !!!
4 FIOXLOM02 3000 0 11157 0 flx 103 43.75 # no mates , shorter read length, highest GC !!!
4 FIOXLOM03 3000 0 12197 0 flx 103 43.51 # no mates , shorter read length, highest GC !!!
4 FIOXLOM04 3000 0 12727 0 flx 103 43.56 # no mates , shorter read length, highest GC !!!
5 FRLDXKV01 0 0 17349 0 . 255 36.02
6 FT9KOI001 0 0 108826 0 . 256 36.10
7 FW1OXFY01 20000 28000 86127 15825 tit 276 35.40
7 FW1OXFY02 20000 28000 106911 19668 tit 275 35.37
7 FW1OXFY03 20000 28000 109874 20396 tit 271 35.23
7 FW1OXFY04 20000 28000 104797 18933 tit 269 35.23
7 FW1OXFY05 20000 28000 113716 20649 tit 265 35.09
7 FW1OXFY06 20000 28000 110693 20326 tit 270 35.05
7 FW1OXFY07 20000 28000 105931 19176 tit 271 35.18
8 F4H5CMB01 20000 28000 626046 109918 tit 241 36.08
8 F4H5CMB02 20000 28000 628432 118903 tit 256 35.87
. total . . 3,297,077 429,170 . . .
- Clr of the Sff seqs (good qual)
. seqs min q1 q2 q3 max mean n50 sum all 3,214,044 0 240 274 383 2042 294 326 947,254,956 => 9.4X
- Clr of the Frg seqs (good qual , no linker)
seqs min q1 q2 q3 max mean n50 sum all 3,297,077 3 156 248 301 2043 244 275 806,091,347 => 8X mated 858,340 64 107 156 223 612 171 201 147,070,298 unmated 2,438,737 2 207 261 335 2042 268 286 655,723,972
- Locations:
/fs/szattic-asmg4/brugia_malayi/Data/Sff/ /fs/szattic-asmg4/brugia_malayi/Data/Frg/
Contaminant & high copy repeat search
- nucmer -maxmatch "-l 20 -c 65" or "-l 12 -c 24"
Sanger 454 jird(26,879 ctgs) 31,501 197,420 # we'd probably find more contaminated reads if we align all the reads to the whole mouse genome ?? Mt 1,507 2,634 We 49,014 23,249 UniVec ? 661,586 pUC19 134 562,107 # 99% avg identity, 99% of read length HhaI(~320bp) 16,336 69,400 mRNA(264bp) . 59,706 total ~910,000
Assemblies
TIGR/NCBI
- 9X coverage, 856K Sanger traces => 8,200 scaff & 29,808 ctg (avg. scaff=~10K & avg ctg=~3K)
- "scaffolds totaling ~71 Mb of data with a further ~17.5 Mb of contigs not integrated into any scaffold (orphan contigs)" (Science 2007)
- NCBI AAQA00000000 AAQA01000001-AAQA01029808
* 26,879 good ctgs * 2,929 jird contaminants (Example: AAQA01001321 : mouse 99%id hits)
- Stats
. elem min q1 q2 q3 max mean n50 sum ctg.len(good) 26879 200 836 1005 1495 611244* 3241.17 18986 87,119,350 ctg.len(cont) 2929 200 527 675 820 8994 740.04 762 2,167,588
. elem min q1 q2 q3 max mean n50 sum
ctg.gc%(good) 26878 0.00 24.77 28.56 32.27 72.30 28.86 29 .
ctg.gc%(cont) 2929 18.09 39.16 43.59 48.35 75.96 44.10 44 .
- Location
/fs/szasmg3/dpuiu/Brugia_malayi/Assembly/TIGR/ <-> NCBI
PITT
- Date: 11/05/08
- Stats:
elem min q1 q2 q3 max mean n50 sum scf.len 3170 2000 2917 4483 14471 6534162* 22916 112914* 72,643,770 (66,051,795bp without gaps) scf.gc% 3170 15.70 25.53 28.17 30.95 66.60 28.46 28 .
- Location:
/fs/szasmg3/dpuiu/Brugia_malayi/Assembly/PITT/
CBCB CA 5.1 Sanger
- Assembler: wgs 5.1
- Date: 2008/08/26
- Input: filtered Sanger reads
- better assembly than the published one
- repeat Hha appears in a few dozen contigs but not in tandem
- Stats:
. elem min q1 q2 q3 max mean n50 sum scf 10317 935 1215 1538 3462 3890532 8018.85 41716 82,730,474 scf2K+ 3656 2001 3181 5733 18904 3890532 20189.57 50293 73,813,083 ctg 12753 273 1245 1632 3873 376744 6113.39 24748 77,964,006 deg 9661 65 858 949 1023 72494 1240.97 1008 11,988,997 singl 134119 (11.43%) reads 1178192(100%)
- Location:
/fs/szasmg3/dpuiu/Brugia_malayi/Assembly/CBCB/2008_0826_CA/
CBCB CA 6.0 454 (failed)
- Assembler: wgs 6.0-beta
- Input: 3,297,077 454 sffToCA processed reads
- Locations:
ginkgo:/scratch1/brugia_malayi/Assembly/454/CA.failed/ /scratch1/ -> umiacsfs01:/xraid03 ginkgo: 32 proc, 128G mem
genome6.umd.edu:/genome6/raid/dpuiu/Brugia_malayi/Assembly/CA.bog/ genome6: 32 proc, 256G mem
- Problem: high frequency contamination & repeats
obtMerThreshold, ovlMerThreshold set on auto (default) !!! runCA estimated them to: (see runCA.log) Reset OBT mer threshold from auto to 37235. Reset OVL mer threshold from auto to 43186.
=> olap-from-seeds very memory/cpu intensive!!! Example: 6 jobs: each is 2 thread, ~ 20G mem merOverlapperSeedConcurrency=6 => 6 jobs merOverlapperExtendBatchSize=20000
$ ps -C olap-from-seeds PID %MEM RSZ(KB) %CPU STIME TIME CMD 13158 0.0 1132 0.0 10:21 00:00:00 /bin/sh 1-overlapper/olap-from-seeds.sh 90 ... 13163 0.0 1136 0.0 10:21 00:00:00 /bin/sh 1-overlapper/olap-from-seeds.sh 95 13199 15.6 20675720 133 10:21 02:46:39 olap-from-seeds -a -b -t 2 -S 1-overlapper/asm.merStore -c 3-overlapcorrection/0092.frgcorr.WORKING -o 1-overlapper/olaps/0092.ovb.WORKING.gz asm.gkpStore 1820001 1840000 ... 13205 15.2 20139808 138 10:21 02:52:35 olap-from-seeds -a -b -t 2 -S 1-overlapper/asm.merStore -c 3-overlapcorrection/0091.frgcorr.WORKING -o 1-overlapper/olaps/0091.ovb.WORKING.gz asm.gkpStore 1800001 1820000
$ vmstat procs -----------memory---------- ---swap-- -----io---- --system-- -----cpu------ r b swpd free buff cache si so bi bo in cs us sy id wa st 9 0 9534576 121356 10124 2743568 5 23 7 23 36 69 16 0 83 1 0
$ free
total used free shared buffers cached
Mem: 132168632 132043920 124712 0 10064 2934076
-/+ buffers/cache: 129099780 3068852
Swap: 67108856 8842720 58266136
- 1-overlapper
#overlaps/read
0+ 3297077
0 1282138 !!! more than 1/3 of reads have no overlaps
1+ 2014939
100+ 561195
#seqs #maxOvl
100+ 561196 99376
UniVec 359340 99376
pUC19 355392 99376
contam 116857 97757
Hha1 50474 29186
mRNA 36296 12935
other 46214 4720
- 8-consensus failed on utg
- 9-terminator => asm.asm.FAILED
CBCB CA 6.0 454
- Input: 3,297,077 454 sffToCA processed reads
- obtMerThreshold=200, ovlMerThreshold=60
. elem min q1 q2 q3 max mean n50 sum scf.len 1406 585 1077 1178 1416 23313 1583.79 1354 2226813 ctg.len 1513 280 1074 1177 1404 4952 1309.21 1257 1980831 deg.len 192840 63 173 262 368 5688 290.91 339 56098824 utg.len 194394 63 174 263 373 5691 298.80 348 58084744 singl.len 1039561 64 162 249 277 1519 244.45 265 254119265 seq.len 2618840 64 140 239 287 1519 234.61 268 614416134 . elem min q1 q2 q3 max mean n50 sum scf.gc% 1406 17.45 29.91 32.75 35.47 69.26 33.21 33 . ctg.gc% 1513 17.45 30.01 32.54 34.93 69.26 32.94 33 . deg.gc% 192840 4.08 27.62 33.33 40.25 80.99 34.17 35 . utg.gc% 194394 4.08 27.65 33.33 40.20 80.99 34.16 35 . singl.gc% 1039561 0.00 34.62 43.11 49.35 86.76 41.58 45 . seq.gc% 2618840 0.00 29.44 37.01 45.48 86.76 37.36 40 .
CBCB CA 6.0 Sanger
- Input: 1,178,192 Sanger "clean" reads
- obtMerThreshold=200, ovlMerThreshold=60
. elem min q1 q2 q3 max mean n50 sum scf.len 10148 904 1235 1578 3564 2324103 8619.28 45146 87468465 scf.len2K+ 3701 2001 3123 6500 20644 2324103 21272.70 52830 78730278 ctg.len 12659 274 1265 1658 3811 565900 6478.56 27474 82012118 deg.len 8497 65 820 914 987 82776 1211.94 980 10297846 utg.len 40545 64 943 1247 1751 306650 2459.44 4231 99718105 singl.len 89048 64 545 702 812 1181 661.17 753 58875872 seq.len 1173341 64 716 827 915 1222 790.88 850 927974911 . elem min q1 q2 q3 max mean n50 sum scf.gc% 10148 13.21 25.02 28.34 32.07 65.75 29.11 29 . ctg.gc% 12659 13.21 24.86 28.22 31.70 65.75 28.78 28 . deg.gc% 8497 7.87 23.87 28.90 33.86 77.13 29.60 30 . utg.gc% 40545 1.39 25.26 29.02 32.72 77.13 29.29 30 . singl.gc% 89048 0.00 31.87 38.97 45.31 99.15 38.52 41 . seq.gc% 1173341 0.00 29.41 32.40 35.18 99.15 32.46 33 .
CBCB CA 6.0 Sanger+454 (Best so far)
- Input: 1,178,192 Sanger "clean" reads ; 3,297,077 454 sffToCA processed reads
- obtMerThreshold=400, ovlMerThreshold=120
. elem min q1 q2 q3 max mean n50 sum scf.len 10254 283 1149 1472 2552 1789338 9181.77 108218 94149853 scf.len2K+ 3193 2000 2728 4513 14385 1789338 26585* 150310* 84,886,293* (75,044,213bp without gaps) ctg.len 13607 66 1181 1578 3187 522243 6196.05 30726 84309668 deg.len 114780 63 195 279 478 25405 367.61 478 42193843 utg.len 157710 63 231 383 885 131968 901.05 1602 142104793 singl.len 904272 64 199 255 306 1690 282.94 278 255858016 seq.len 3849841 64 202 281 662 1690 410.15 699 1579010968 . elem min q1 q2 q3 max mean n50 sum scf.gc% 10254 13.80 24.59 28.34 32.95 69.26 29.44 29 . ctg.gc% 13607 13.80 24.54 28.17 32.02 69.26 28.95 28 . deg.gc% 114780 3.75 29.39 38.22 45.45 80.99 37.40 41 . utg.gc% 157710 0.00 27.63 34.59 43.08 80.99 35.28 38 . singl.gc% 904272 0.00 36.60 43.75 49.49 99.33 42.44 46 . seq.gc% 3849841 0.00 29.34 34.26 41.98 99.33 35.73 36 .
- Location:
ginkgo:/scratch1/brugia_malayi/Assembly/hybrid/CA/
- Rerun using bogBadMateDepth = 4 (default is 7) at Aleksey's advice; utgs are smaller; failed in cgw "scaffolder failed" message
CBCB newbler deNovo 454 (failed)
- Still running after 7 days (killed)
Detangling alignments... -> Level 2, Phase 8, Round 1...
PID %MEM RSZ %CPU STIME TIME CMD 4576 2.1 1427100 94.4 Feb16 7-03:26:12 /fs/szdevel/dpuiu/454/bin/runProject .
CBCB newbler refMapper 454
- Assembler: newbler 2.3
- Host: CBCB walnut server
- Input
# NCBI ref assembly
ctgs min q1 q2 q3 max mean n50 sum
All 26,879 200 836 1005 1495 611244 3241 18986 87,119,350
#Sff reads . seqs min q1 q2 q3 max mean n50 sum All 3,214,044 0 240 274 383 2042 294.72 326 947,254,956
- Output
#Ctg stats . ctgs min q1 q2 q3 max mean n50 sum All 101,286 100 236 323 530 7013 433.36 535 43,893,507
#Trimmed read stats . seqs min q1 q2 q3 max mean n50 sum All 3,898,373 1 45 163 264 1995 167.84 265 654319111 Full|Partial 1,085,167 20 119 216 285 706 214.13 271 232364804 Chimeric|Repeat|Unmapped|TooShort 2,015,920 20 111 221 276 1995 208.85 263 421032444 Deleted 797,286 1 1 1 1 19 1.16 1 921863
#Trimmed read counts
count %
All 3898373 100
Chimeric 25460 0.65
Deleted 797286 20.45 !!!
Full 1001745 25.7
Partial 83422 2.14
Repeat 406119 10.42 !!!
TooShort 14031 0.36
Unmapped 1570310 40.28 !!!
#Mate pair counts
count %
BothUnmapped 301390 42.86
OneUnmapped 110922 15.77
MultiplyMapped 108641 15.45
FalsePair 106249 15.11
TruePair 75992 10.81
Files
/fs/szattic/asmg1/adelcher/Genomes/Brugia : Art's files /fs/sztmpscratch/cole/tarchive_download/brugia_malay : Cole's files /fs/szasmg3/dpuiu/Brugia_malayi/ : Daniela's files /scratch1/brugia_malayi/Data/ : ftp PITT data /fs/szattic-asmg4/brugia_malayi : ftp PITT data (as well)