BitMapperBS icon indicating copy to clipboard operation
BitMapperBS copied to clipboard

Segmentation fault in alignment

Open wbvguo opened this issue 2 years ago • 4 comments

Dear BitmapperBS developer,

I installed the latest bitmapperBS (v1.0.2.3) by git clone the source code and make (GCC version 9.4.0), the installation works well, I can see bitmapperBS and psascan in the folder. The genome indexing part is also good, with the following message:

.........
The length of Genome is 3217346917:
Open ../../test/test_idx/genome.index.bs sucessfully!
text_length=6434693834
rm: cannot remove 'tmp_ref.tmp.sa5': No such file or directory
the binary of psascan exists...
.........
Computation finished. Summary:
  elapsed time: 776.53s (0.1265s/MiB)
  speed: 7.90MiB/s
text_length =6434693834, file_length=6434693834
text_length = file_length...sucess...
SA has been generated!
bwt_length=201084183,occ_byte_length=50271047, bwt_length+occ_byte_length=251355230
Occ and BWT have been built!
write SA_length=2139726539, shapline=3845490917
SA_flag has been built
Occ has been writed!
BWT has been writed!
cc&sharp_line has been writed!
sparse_suffix_array_length=804336730
SA_flag_iterater=125677616
hash table has been writed!

But the alignment runs into trouble (Segmentation fault (core dumped)):

fileName[0]: ../../test/test_idx
fileName[1]: ../../test/test_idx.index
Read qualities are encoded by Phred+33...
GapOpenPenalty: 5, GapExtensionPenalty: 3, MistMatchPenaltyMax: 6
MistMatchPenaltyMin: 2, N_Penalty: 1
 Read files are in compressed FASTQ format...
Open ../../test/test_idx/genome.index sucessfully!
Start load hash table!
refGenLength = 3217346917
Open ../../test/test_idx/genome.index.bs.pac sucessfully!
shapline=3845490917
BWT has been loaded!
hash table has been loaded!
SA_length=2139726539
sparse_suffix_array_length=804336730
SA_flag_iterater=125677616
high_occ_table_length=196372
Segmentation fault (core dumped)

I also tried the install and run it with conda, the installation is good, but also run into the following error.

$bitmapperBS --help
Illegal instruction (core dumped)

I tested all above on a server running on ubuntu system (20.04) that supports sse4_2. I also tested on another cluster node running on centos 7, which supports AVX2, it has the same alignment issue.

May I ask if you have any idea to make it work?

Thanks!

wbvguo avatar Apr 29 '22 23:04 wbvguo