Improve criteria for when to stop remap iterations
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 21
- Forks
- 11
- Avg merge
- 12h 5m
- Merged PRs (30d)
- 24
Description
As an example, sample 88067A_HCV from the 29-May-2017.M01841 run has a drop in coverage at the start of the E2 region in HCV-2. That drop wasn't there when it was analysed with version 7.6; it appeared with version 7.7. When I investigated, I found that version 7.6 ran two remap iterations, and 7.7 only ran one. After the first iteration, version 7.6 still had the gap, so the second iteration was very useful.
Version 7.6 ran the second iteration, because one of the HIV seeds got dropped. It just happened to also help the HCV coverage. In 7.7, the HIV had already been mapped by the G2P step, and we already had over 95% mapped reads, so it didn't do a second remap iteration.
Ideas for better criteria:
- Always remap at least twice. (Too wasteful?)
- Don't consider reads with a lot of soft clipping as mapped.
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 by reproducing the 88067A_HCV result from the 29-May-2017.M01841 run and comparing remap behavior between versions 7.6 and 7.7. Investigate how mapped reads and dropped HIV seeds determine whether another iteration runs. Done means an agreed stopping criterion addresses the HCV coverage regression without unnecessary remapping.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- bioinformatics
- Issue type
- Feature
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100