[Mapass-help] sh: line 1: 6472 Segment default
Status: Beta
Brought to you by:
lh3lh3
|
From: Liu T. <ltn...@16...> - 2007-06-15 08:57:40
|
Hi, mapass-admin:
When I run mapass2, using the command:
mapass2.pl easyrun -d assembly_rst reference.fa sequence_rst.fastq
Because of my misunderstanding of the introduction of fastq's format. I got the following errors:
===================================================================================================
[easyrun] Converting the reference sequence 'reference.fa' to binary fasta format (.bfa)...
89 sequences have been converted.
[easyrun] Converting fastq file 'sequence_rst.fastq' to binary fastq format (.bfq)...
2463084 sequences were loaded.
[easyrun] Aligning reads-1.bfq to the reference...
[ma_load_reads] loading reads...
[ma_load_reads] set length of the first read as 33.
[ma_load_reads] 2463084*2 reads loaded.
[ma_longread2read] encoding reads... 4926168 sequences processed.
[match_core] round 1/3...
[match_core] making index...
sh: line 1: 6472 Segmentation fault /usr/local/bin/mapass2 match reads-1.map ref.bfa reads-1.bfq
Died at /usr/local/bin/mapass2.pl line 109.
=====================================================================================================
My fault is : the quality of the sequence was converted by the code
you pubilished on the home site of mapass:
======================================================================================================
my$std_qseq ='';
while($sol_qseq =~ /(\S)/g){
$std_qseq .=chr(int(33 +10*log(1+pow(10,(ord($1)-64)/10.0))/log(10)+.499));
}
======================================================================================================
It may be not a bug of program though, supply it here to let you know that
it will report the imformation above, and the reason is ambiguous.
And now I have corrected the error, then go through the program and get some assemblies.
The result seems good.
I'd glad to help to improve this project. Thank you.
Best Regards!
Liu Tao |