Skip to content

error with bayesTyper cluster #56

Open
@zhengluo-zz

Description

@zhengluo-zz

Hi,
when I use bayesTyper cluster, the following error occurs. Can you tell me how this problem is caused and how to resolve it?
command bayesTyper cluster -v ../../00_data/merge.c.vcf.gz -s samples.tsv -g ../../00_data/refgenome.fa -p 10
error

[09/10/2024 19:43:31] You are using BayesTyper (v1.5 1b04081ca4c64975270d267dbe9d80bd00bddaf9)

[09/10/2024 19:43:31] Seeding pseudo-random number generator with 1728474211 ...
[09/10/2024 19:43:31] Setting the kmer size to 55 ...

[09/10/2024 19:43:31] Parsed information for 1 sample(s)

[09/10/2024 19:43:31] Parsing reference genome ...
[09/10/2024 19:43:51] Parsed 10 reference genome chromosomes(s) (2131846805 nucleotides)

[09/10/2024 19:43:51] Parsing decoy sequence(s) ...
[09/10/2024 19:43:51] Parsed 0 decoy sequence(s) (0 nucleotides)

[09/10/2024 19:55:04] Setting the number of inference units to 1 across 3940713 variants ...

[09/10/2024 19:55:13] Maximum resident set size: 9.8092 Gb


[09/10/2024 19:55:13] Parsing variants in unit 1 ...
bayesTyper: /public/home/zhengluo/SOFTWARE/BayesTyper/src/bayesTyper/VariantFileParser.cpp:407: void VariantFileParser::parseVariants(ProducerConsumerQueue<std::vector<std::unordered_map<unsigned int, VariantCluster*>*>*>*, uint, const Chromosomes&): Assertion `count(alt_alleles.begin() + i + 1, alt_alleles.end(), alt_alleles.at(i)) == 0' failed.
/public/home/zhengluo/.lsbatch/1728474209.35253913: line 8: 660249 Aborted                 (core dumped) bayesTyper cluster -v ../../00_data/merge.c.vcf.gz -s samples.tsv -g ../../00_data/refgenome.fa -p 10

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions