Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Move HCV resistance report into MiCall pipeline #412

Closed
15 tasks done
donkirkby opened this issue Aug 31, 2017 · 1 comment
Closed
15 tasks done

Move HCV resistance report into MiCall pipeline #412

donkirkby opened this issue Aug 31, 2017 · 1 comment

Comments

@donkirkby
Copy link
Member

donkirkby commented Aug 31, 2017

Convert the HCV scoring algorithm to use the Stanford ASI rules format with extensions.

The current algorithm is implemented here:

  • fdrive/RESIST/ruby/ckwlib/hcv_algorithm.rb

The driver script is fdrive/RESIST/MISEQ/miseq_gen_results.rb.

  • convert rules
  • change report title and small print from HIV to HCV
  • add drug names
  • display genotype
  • display multiple genotypes, one per page
  • decide how to display drugs with no rules for the given genotype (Not Indicated: genotype X)
  • add failure report (rules in comment below)
  • add unknown mutation rules
  • special rules for 6e: "resistance likely" and "resistance interpretation not available"
  • combine midi overlap
  • display "not indicated" even when region has bad coverage
  • compare results to Ruby version
  • update BaseSpace driver to use new HCV rules
  • display HCV rules in web page on MiCall web site
  • move key coverage positions into scoring algorithm (moved to Get key positions from resistance rules #431)
@donkirkby donkirkby added this to the near future milestone Aug 31, 2017
donkirkby added a commit that referenced this issue Jan 6, 2018
Also display Canadian drug names and updated HCV disclaimer.
@donkirkby
Copy link
Member Author

Summary of email discussion with Anita and others about failure scenarios:

  1. One genotype has good coverage on NS3, NS5a, NS5b, and other genotypes have bad coverage. Display the genotype with the good coverage, and ignore the genotypes with bad coverage.
  2. No genotypes have good coverage, and some genotypes have bad coverage. Display a one-page failure report with a message, "sequence does not meet quality-control standards".
  3. Nothing mapped. Same as sampleSheetParser fails to parse combined sample names #2
  4. One genotype has good coverage on some regions and bad coverage on others. For example, good coverage on NS3 and NS5a, bad coverage on NS5b. Display the resistance calls for drugs in the good regions, and display "sequence does not meet quality-control standards" for drugs in the bad regions.

There is a special rule for genotype 6e with Sofosbuvir (Epclusa): it should display the message, "resistance interpretation not available".

donkirkby added a commit that referenced this issue Jan 11, 2018
Also add new flags to HCV rules.
Generate failure report.
donkirkby added a commit that referenced this issue Jan 12, 2018
donkirkby added a commit that referenced this issue Jan 17, 2018
Change resistance level numbers to match Conan's.
Make sure reported mutations have wild type.
donkirkby added a commit that referenced this issue Jan 17, 2018
Add utility for rerunning resistance analysis.
donkirkby added a commit that referenced this issue Jan 18, 2018
Fix broken tests by upgrading pyvdrm to support bangs.
donkirkby added a commit that referenced this issue Jan 18, 2018
Use sample sheet to find matching sample names.
donkirkby added a commit that referenced this issue Feb 2, 2018
Rerun resistance interpretation on multiple run folders.
donkirkby added a commit that referenced this issue Feb 2, 2018
donkirkby added a commit that referenced this issue Feb 7, 2018
Add some microtest samples that will successfully combine main and midi.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant