What Causes The Bwa Error: "Assertion `Re - Rb == Rlen' Failed."
Neither bwa 6.0.2 nor bwa 7.0.5 worked for my alignment. The size of my genome is about 1.5M and the number reads are 25000000.
However, when I tried to reduce the number of reads to 15000000 it worked with bwa 6.0.2 (not sure about bwa 7.0.5).
Would anyone please help me figure out what the problem is?
ERROR MESSAGE is as follows:
[bwa_sai2sam_pe_core] refine gapped alignments... bwa: bwase.c:182: bwa_refine_gapped_core: Assertion `re - rb == rlen' failed.
1 answer
Log in to answer this question.
More posts like this
-
Bwa Mem -M Option
written by Varun Gupta
-
How To Get The Correct Alignment Statistics After Using Bwa Mem With -M Option
written by skm770 •
-
Bwa Mem: How To Specify The Fastq Phred Format?
written by Giovanni M Dall'Olio
-
Bwa Installation Failed After "Make" Command In Linux
written by Tonyzeng •
-
Bwa Alignments With Gaps
written by lob •
-
Do Different Versions Of Bwa Create Differently Named Index Files?
written by mathiasf52 •
-
Bwa: Should Sai/Sam Sequences Processed Match?
written by jomaco
-
Bwa Sampe Segmentation Fault
written by a11msp
-
Bwa Alignment: Gap Vs Soft-Clipped Sequences
written by michealsmith
-
Predict/Estimate/Find Bwa Parameters From Bam Or Sam File
written by Houkto •
if that's a bug, you'd better ask the samtools-dev mailing list.
At a guess it sounds like there is a malformed sequence in the input causing it to crash. Check if it is the number of sequences that are the problem or a specific one. Split the input set into many chunks, run each chunk and see if one is consistently responsible. Then subdivide it until you can track down the problem.