Please get the binary latest version at https://github.com/armintoepfer/QuasiRecomb/releases
Description:
RNA viruses are present in a single host as a population of different but related strains. This population, shaped by the combination of genetic change and selection, is called quasispecies. Genetic change is due to both point mutations and recombination events. We present a jumping hidden Markov model that describes the generation of the viral quasispecies and a method to infer its parameters by analysing next generation sequencing data. We offer an implementation of the EM algorithm to find maximum a posteriori estimates of the model parameters and a method to estimate the distribution of viral strains in the quasispecies. The model is validated on simulated data, showing the advantage of explicitly taking the recombination process into account, and tested by applying to reads obtained from experimental HIV samples.
Features
- Reconstruction of haplotypes with their relative frequency
- SNV posterior computation
- Full BAM format support