strobealign icon indicating copy to clipboard operation
strobealign copied to clipboard

Extension alignments longer than 2000 bp

Open marcelm opened this issue 2 years ago • 0 comments

Sometimes ssw_align() is asked to compute an alignment for a reference sequence longer than 2000 bp here: https://github.com/ksahlin/strobealign/blob/8d3ef7e113abe9159a38343d970e864fc5afaab4/src/aln.cpp#L35-L43

This gets triggered in the Drosophila test dataset for read @SRR6055476.50009. The chunk needs to start a couple of reads earlier. Putting reads @SRR6055476.50004 to @SRR6055476.50009 (paired end) in paired-end FASTQ files reproduces it reliably.

As discussed in #168

marcelm avatar Dec 05 '22 14:12 marcelm