Skip to content
This repository has been archived by the owner on Jul 17, 2023. It is now read-only.

Filter bad-cigar reads in bam streamer #288

Open
wants to merge 1 commit into
base: master
Choose a base branch
from

Conversation

fbattke
Copy link

@fbattke fbattke commented Jul 5, 2022

These reads are ok according to SAM/BAM spec, but Manta can't work with them.

This fixes issue #137 "Manta crashes on input files realigned with ABRA #137"
as well as issue #184 "Unclear how to deal with CIGAR strings with InDels at the end of reads"

The fix is to remove offending alignments in the bam streamer instead of throwing an exception much later in the code. This has the same effect as pre-filtering the alignments using awk (as described in issue 184), but it is much faster.

These reads are ok according to SAM/BAM spec, but Manta can't work with them
@fbattke
Copy link
Author

fbattke commented Jul 6, 2022

We have tested this change on 100 cases and all output was identical to the original workaround using awk to filter alignments before calling manta.

@fbattke fbattke marked this pull request as ready for review July 6, 2022 08:07
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant