bioforensics / MicroHapulator

Tools for empirical microhaplotype calling, forensic interpretation, and simulation.
https://microhapulator.readthedocs.io/
Other
6 stars 1 forks source link

Dynamic cutoff calculation for genotyping #51

Closed standage closed 5 years ago

standage commented 5 years ago

The previously implemented automatic cutoff calculation works well for genotyping loci with high effective coverage (few discarded reads). But it has failed in a few cases where loci have a low effective coverage (microhaps with a large extent, resulting in many reads that do not span the entire microhap and thus must be discarded). In these cases, a fixed threshold performs better.

Rather than choosing the automatic approach or the fixed threshold approach, this update introduces a new dynamic approach that chooses the automatic or fixed threshold at each locus based on effective coverage.


codecov-io commented 5 years ago

Codecov Report

Merging #51 into master will increase coverage by <.01%. The diff coverage is 100%.

Impacted file tree graph

@@            Coverage Diff             @@
##           master      #51      +/-   ##
==========================================
+ Coverage   97.18%   97.18%   +<.01%     
==========================================
  Files          26       26              
  Lines         817      818       +1     
  Branches      151      151              
==========================================
+ Hits          794      795       +1     
  Misses         11       11              
  Partials       12       12
Impacted Files Coverage Δ
microhapulator/cli/type.py 100% <ø> (ø) :arrow_up:
microhapulator/type.py 97.1% <100%> (ø) :arrow_up:
microhapulator/genotype.py 96.42% <100%> (+0.01%) :arrow_up:

Continue to review full report at Codecov.

Legend - Click here to learn more Δ = absolute <relative> (impact), ø = not affected, ? = missing data Powered by Codecov. Last update 12d7958...e13e0f5. Read the comment docs.