sebhtml / ray

Ray -- Parallel genome assemblies for parallel DNA sequencing
http://denovoassembler.sf.net
Other
65 stars 12 forks source link

avoid deadends #233

Closed sebhtml closed 10 years ago

sebhtml commented 10 years ago

deadend

sebhtml commented 10 years ago

How is this happening anyway ?

sebhtml commented 10 years ago

http://maupiti.crchul.ulaval.ca:2280/client/?map=4&section=2&region=7&location=49073&depth=0&zoom=0.6072467142873316

sebhtml commented 10 years ago

dataset: AZPAE14881-1 Rank 01 1087 seeds (basic) 602 (linear detector) XXX (3-depth tree)

dataset: ANL NormalStool with raw code (-13): Rank 1 has 3583 seeds to register. (basic) with new patch: , Rank 1 has 2974 seeds (linear detector) XXX (3-depth tree)

sebhtml commented 10 years ago

before

Rank 1 Skipped paths because of dead end for head: 2259 Rank 1 Skipped paths because of dead end for tail: 2254 Rank 1 Skipped paths because of two dead ends: 3747 Rank 1 Skipped paths because of bubble weak component: 0 Rank 1 Skipped paths because of short length: 935989 Rank 1 Skipped paths because of bad ownership: 619 Rank 1 Skipped paths because of low coverage: 0 Rank 1 Eligible paths: 602

after

Rank 1 Skipped paths because of dead end for head: 2273 Rank 1 Skipped paths because of dead end for tail: 2276 Rank 1 Skipped paths because of two dead ends: 3758 Rank 1 Skipped paths because of bubble weak component: 0 Rank 1 Skipped paths because of short length: 935942 Rank 1 Skipped paths because of bad ownership: 619 Rank 1 Skipped paths because of low coverage: 0

after

sebhtml commented 10 years ago

09

before

Rank 9 Skipped paths because of dead end for head: 2346 Rank 9 Skipped paths because of dead end for tail: 2286 Rank 9 Skipped paths because of two dead ends: 3733 Rank 9 Skipped paths because of bubble weak component: 0 Rank 9 Skipped paths because of short length: 938686 Rank 9 Skipped paths because of bad ownership: 588 Rank 9 Skipped paths because of low coverage: 0 Rank 9 Eligible paths: 589

Rank 9 Skipped paths because of dead end for head: 2360 Rank 9 Skipped paths because of dead end for tail: 2301 Rank 9 Skipped paths because of two dead ends: 3750 Rank 9 Skipped paths because of bubble weak component: 0 Rank 9 Skipped paths because of short length: 938640 Rank 9 Skipped paths because of bad ownership: 588 Rank 9 Skipped paths because of low coverage: 0 Rank 9 Eligible paths: 589

after

sebhtml commented 10 years ago

wait for ANL NormalStool result, then close this.

sebhtml commented 10 years ago

[boisver1@ip03-mp2 ANL]$ pwd /mnt/scratch_mp2/corbeil/corbeil_group/nne-790-ab/projects/ANL [boisver1@ip03-mp2 ANL]$ qsub NormalStool-17.sh 22839.mp2.m [boisver1@ip03-mp2 ANL]$ showq | grep boisv 22839 boisver1 Running 11 1:23:56:13 Mon Feb 24 17:51:52

sebhtml commented 10 years ago

There is a reduction in the number of seeds that are deadends.

[boisver1@ip03-mp2 ANL]$ grep Eligible NormalStool-13.o25003|sort|head Rank 0 Eligible paths: 3641 Rank 100 Eligible paths: 3417 Rank 101 Eligible paths: 3599 Rank 102 Eligible paths: 3495 Rank 103 Eligible paths: 3518 Rank 104 Eligible paths: 3585 Rank 105 Eligible paths: 3428 Rank 106 Eligible paths: 3575 Rank 107 Eligible paths: 3548 Rank 108 Eligible paths: 3495 [boisver1@ip03-mp2 ANL]$ grep Eligible NormalStool-16.o27065|sort|head
Rank 0 Eligible paths: 2999 Rank 100 Eligible paths: 2823 Rank 101 Eligible paths: 2971 Rank 102 Eligible paths: 2865 Rank 103 Eligible paths: 2893 Rank 104 Eligible paths: 3005 Rank 105 Eligible paths: 2860 Rank 106 Eligible paths: 2934 Rank 107 Eligible paths: 2924 Rank 108 Eligible paths: 2903 [boisver1@ip03-mp2 ANL]$ grep Eligible NormalStool-17.o22839|sort|head
Rank 0 Eligible paths: 2999 Rank 100 Eligible paths: 2823 Rank 101 Eligible paths: 2966 Rank 102 Eligible paths: 2860 Rank 103 Eligible paths: 2888 Rank 104 Eligible paths: 3003 Rank 105 Eligible paths: 2859 Rank 106 Eligible paths: 2930 Rank 107 Eligible paths: 2920 Rank 108 Eligible paths: 2897