bioinformatics-centre / BayesTyper

A method for variant graph genotyping based on exact alignment of k-mers
86 stars 7 forks source link

genotyping error: intercluster_diploid_kmer_stats.at(sample_idx).at(bias _idx).front().getCount() == 0 #12

Open biobai opened 5 years ago

biobai commented 5 years ago

Hi When I run the BayesTyper v1.4, I got this error. Could anyone please help me to fix this? Thanks

[14/02/2019 06:13:18] You are using BayesTyper (v1.4)

[14/02/2019 06:13:18] Seeding pseudo-random number generator with 1550146398 ... [14/02/2019 06:13:18] Setting the kmer size to 55 ...

[14/02/2019 06:13:18] Parsed information for 1 sample(s)

[14/02/2019 06:13:18] Parsing reference genome ... [14/02/2019 06:14:49] Parsed 24 reference genome chromosomes(s) (3095677412 nucl eotides)

[14/02/2019 06:14:49] Parsing decoy sequence(s) ... [14/02/2019 06:15:03] Parsed 62 decoy sequence(s) (41777093 nucleotides)

[14/02/2019 06:15:14] Maximum resident set size: 3.37096 Gb

[14/02/2019 06:15:14] Parsing variant clusters ... [14/02/2019 06:16:04] Parsed 1829850 variant clusters (5052955 variants)

[14/02/2019 06:16:12] Parsing parameter kmers ... [14/02/2019 06:16:24] Parsed 1000000 kmers

[14/02/2019 06:16:24] Maximum resident set size: 25.0474 Gb

[14/02/2019 06:16:24] Counting kmers in variant cluster paths ... [14/02/2019 06:29:30] Counting kmers in inter-cluster regions and decoy sequence (s) ...

[14/02/2019 06:34:57] Parsing KMC table containing 8070384125 kmers for sample U DN765115 ...

[14/02/2019 07:48:11] Classifying kmers in variant cluster paths ... [14/02/2019 07:50:42] Out of 211574032 kmers:

    - 180492874 have a match to a single variant cluster
    - 17803661 have a match to single variant cluster group and multiple var

iant clusters

    - 1011777 have match to at least one variant cluster and has match to a

decoy sequence (not used for inference)

[14/02/2019 07:50:42] Maximum resident set size: 27.0064 Gb

[14/02/2019 07:50:42] Estimating genomic haploid kmer count distribution(s) from parameter kmers ...

bayesTyper: /isdata/kroghgrp/jasi/bayesTyper/code/releases/v1.4_static/BayesType r-1.4/src/bayesTyper/CountDistribution.cpp:92: void CountDistribution::setGenomi cCountDistributions(const std::vector<std::vector<std::vector > >&, c onst string&): Assertion `intercluster_diploid_kmer_stats.at(sample_idx).at(bias _idx).front().getCount() == 0' failed.

Best, Youhuang

jonassibbesen commented 5 years ago

Hi Youhuang,

Thanks for writing. This assert is the result of a bug introduced in v1.4. The latest version (1.4.1) fixed the bug.

Let me know if you are still having problems after switching to v1.4.1

Best,

Jonas