Skip to content

[Fix] correct CIGAR string in SAM format as well as correct header in…#157

Merged
h-2 merged 2 commits intoseqan:lambda3from
sarahet:fix_sam_output
Jan 30, 2020
Merged

[Fix] correct CIGAR string in SAM format as well as correct header in…#157
h-2 merged 2 commits intoseqan:lambda3from
sarahet:fix_sam_output

Conversation

@sarahet
Copy link
Copy Markdown
Member

@sarahet sarahet commented Jan 30, 2020

… bam format

Already known problem with macth/record qLength entries.
Wrong number of sequences used for the BAM header

@sarahet sarahet requested a review from h-2 January 30, 2020 05:02
Copy link
Copy Markdown
Member

@h-2 h-2 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We should not switch to m.qLength, that this the deprecated member. Using r.qLength is the right thing to do, we just need to make sure it has the right value. Or we need to do std::ranges::size(CORRECTseqs[i]). Whatever CORRECT means in the context.

seqan::resize(subjSeqLengths, globalHolder.indexFile.seqs.size());
SEQAN_OMP_PRAGMA(parallel for simd)
for (size_t i = 0; i < globalHolder.redSbjSeqs.size(); ++i)
for (size_t i = 0; i < globalHolder.indexFile.seqs.size(); ++i)
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

wow, good catch, this looks really weird!

@sarahet sarahet requested a review from h-2 January 30, 2020 16:01
@h-2 h-2 merged commit b6d3403 into seqan:lambda3 Jan 30, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants