SAM file flags show reads not aligned properly
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 21
- Forks
- 11
- Avg merge
- 12h 5m
- Merged PRs (30d)
- 24
Description
Sometimes a read shows up in prelim.csv with a flag without the 0x2 bit. This means that the read is not mapped in a proper pair. This flag is described in the manual (PDF) as, "each segment properly aligned according to the aligner".
The samtools mpileup command does not include reads in the pileup if this flag is not set, so we stopped including them in sam_to_conseq().
What exactly does this flag mean, and should we use the reads without this flag set? If we decide not to use them, filter them out of the remap.csv output, and check whether they need to get filtered anywhere else.
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start with sam_to_conseq() and the generation of prelim.csv and remap.csv. Read the linked SAM specification and Picard flag explanation, then trace where reads lacking 0x2 are handled and whether other outputs include them. Done means the flag behavior is understood, a filtering decision is recorded, and the affected outputs consistently follow it.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- bioinformatics
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 32/100