cfe-lab / cfe-lab/MiCall

Improve criteria for when to stop remap iterations

Open
#405 2 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

enhancement
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

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. 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

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.