[mira_talk] miraSearchESTSNPs aborts

Hi Bastien,

i've tried miraSearchESTSNPs with version 2.9.45 and i get a segmentation 
fault early in step 1 :

Analysing hstat files:
 [0%] ....|.... [10%] ....|.... [20%] ....|.... [30%] ....|.... [40%] 
....|.... [50%] ....|.... [60%] ....|.... [70%] ....|.... [80%] ....|.... 
[90%] ....|.... [100%] 
Hash statistics:
=========================================================
Measured avg. coverage: 12207165

Deduced thresholds:
-------------------
Min normal cov: 6103582
Max normal cov: 18310748
Repeat cov: 23193614
Crazy cov: 97657320
Mask cov: 0

Repeat ratio histogram:
-----------------------
1       3
=========================================================

Assigning values:
Localtime: Mon May 18 14:40:34 2009
 [0%] Segmentation fault


So i've tried with previous version 2.9.44, step 1 runs without problem 
but the process aborts at step 2 with following error :

Building new contig 6162
Localtime: Mon May 18 14:30:24 2009
Unused reads: 9731
miraSearchESTSNPs: contig.C:5686: void Contig::updateCountVectors(int32, 
int32, __gnu_cxx::__normal_iterator<const char*, std::vector<char, 
std::allocator<char> > >, uint32, bool): Assertion `(from) >= (0) && 
(from) < (CON_counts.size())' failed.

Started with non-multicopy.


Any idea where the problem comes from ? Can you fix it ?

Here are my OS details (RedHat) :

Linux sg5.biogemma.fr 2.6.18-128.1.6.el5 #1 SMP Tue Mar 24 12:05:57 EDT 
2009 x86_64 x86_64 x86_64 GNU/Linux

Cheers
jd.

Other related posts: