Dear,
I run BWA 0.7.5a
>bwa aln genome.fasta read1.fastq > read1.sai
>bwa aln genome.fasta read2.fastq > read2.sai
> bwa sampe -r '@RG\tID:s\tLB:ss\tPL:ILLUMINA\tPU:l1\tSM:' genome.fasta read1.sai read2.sai read1.fastq read2.fastq > alignment.sam
but during sampe step it returns me this error:
[bwa_sai2sam_pe_core] refine gapped alignments... bwa: bwase.c:182: bwa_refine_gapped_core: Assertion `re - rb == rlen' failed.
Aborted
anyone has any idea?
I run BWA 0.7.5a
>bwa aln genome.fasta read1.fastq > read1.sai
>bwa aln genome.fasta read2.fastq > read2.sai
> bwa sampe -r '@RG\tID:s\tLB:ss\tPL:ILLUMINA\tPU:l1\tSM:' genome.fasta read1.sai read2.sai read1.fastq read2.fastq > alignment.sam
but during sampe step it returns me this error:
[bwa_sai2sam_pe_core] refine gapped alignments... bwa: bwase.c:182: bwa_refine_gapped_core: Assertion `re - rb == rlen' failed.
Aborted
anyone has any idea?
Comment