Kimimila: A new model to classify ngs short reads by their allele origin

Next generation sequencing (NGS) technologies, often referred to as massively parallel sequencing, are having a huge impact on genomics and clinical applications. These technologies generate billions of short sequences (reads) that are consequently mapped to their corresponding reference genome to find out known and/or novel genomic variants potentially correlated to patients phenotype. DNA fragment library is usually derived from a diploid genome: we refer to genotyping on NGS data as the analytical process to assign the zygosity of identified variants. Current algorithms typically rely on data of the single genomic locus where variants have been called and are based on the condition of independence between variant locus and reads. These strong assumptions might bring to possible inaccuracies throughout the genotyping process. We have therefore developed an efficient assumption-free algorithm based on a kinetic model approach and distance geometry (Kimimila) that delivers the belonging allele for each read using the inference provided by the measure of differences (i.e. variants) among overlapping reads.