A method for detecting and reporting, with a processor, a variant as a repeat of at least two nucleotide patterns in a genomic sequence of a patient sample, the method comprising: (a) identifying a reference repeat pattern P ref = N ∗ l as the repeat of l ( l >=2) genomic patterns N in a genomic region of a human genome reference sequence; (b) obtaining, with a next generation sequencer, n patient sets of next generation sequencing data reads covering the reference repeat pattern genomic region S= { S 1, S 2,..., Si,..., S n} from a pool of n enriched genomic patient samples, each set S i being associated with a patient sample, the number n of enriched genomic patient samples being at least 4; (c) for each patient sample i in the set S of patient samples, measuring the distribution P i of the length of the repeat pattern in the set of next generation sequencing reads S i; (d) for a possible pair of patient samples i and j, j > i, estimating a best-fit model V ij 1 | V ij 2 of the two allele variants for sample i relative to sample j, with a confidence level L ij; (e) for each possible triplet of patient samples i,j > i, k > j, comparing their respective best-fit models [ V ij 1 | V ij 2], V jk 1 | V jk 2, V ik 1 | V ik 2, grouping matching best-fit models into groups of best-fit variant models with an increased confidence level, and iterating the comparison until stable groups of best-fit variant models are formed; (f) identifying the most likely group carrying the wild type variant; (g) for each sample in the group carrying the wild type variant, reporting the sample variant as the wild type reference repeat pattern P ref = N ∗ l. (h) for each sample out of the group carrying the wild type variant, unbiasing the best-fit variant model of the group comprising this sample as a function of the best-fit variant model for the identified wild type group and reporting the sample variant as the unbiased best-fit model variant. (a) identifying a reference repeat pattern P ref = N ∗ l as the repeat of l ( l >=2) genomic patterns N in a genomic region of a human genome reference sequence; (b) obtaining, with a next generation sequencer, n patient sets of next generation sequencing data reads covering the reference repeat pattern genomic region S= { S 1, S 2,..., Si,..., S n} from a pool of n enriched genomic patient samples, each set S i being associated with a patient sample, the number n of enriched genomic patient samples being at least 4; (c) for each patient sample i in the set S of patient samples, measuring the distribution P i of the length of the repeat pattern in the set of next generation sequencing reads S i; (d) for a possible pair of patient samples i and j, j > i, estimating a best-fit model V ij 1 | V ij 2 of the two allele variants for sample i relative to sample j, with a confidence level L ij; (e) for each possible triplet of patient samples i,j > i, k > j, comparing their respective best-fit models [ V ij 1 | V ij 2], V jk 1 | V jk 2, V ik 1 | V ik 2, grouping matching best-fit models into groups of best-fit variant models with an increased confidence level, and iterating the comparison until stable groups of best-fit variant models are formed; (f) identifying the most likely group carrying the wild type variant; (g) for each sample in the group carrying the wild type variant, reporting the sample variant as the wild type reference repeat pattern P ref = N ∗ l. (h) for each sample out of the group carrying the wild type variant, unbiasing the best-fit variant model of the group comprising this sample as a function of the best-fit variant model for the identified wild type group and reporting the sample variant as the unbiased best-fit model variant.