Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Annotation of merged Sniffles 2 VCF files does not work. #357

Closed
holtgrewe opened this issue Feb 21, 2024 · 0 comments · Fixed by #358
Closed

Annotation of merged Sniffles 2 VCF files does not work. #357

holtgrewe opened this issue Feb 21, 2024 · 0 comments · Fixed by #358
Labels
bug Something isn't working

Comments

@holtgrewe
Copy link
Contributor

Describe the bug
Annotating merged Sniffles 2 VCF files gives an error.

To Reproduce
See below for a transcripts

> mehari annotate seqvars --path-db /data/cephfs-1/work/groups/cubi/projects/biotools/mehari/0.21.0/db/ --path-input-ped test.ped --path-input-vcf test.multi.vcf --path-output-tsv test.tsv
2024-02-09T13:24:54.181726Z  INFO Mehari startup -- letting the dromedary off the leash...                                   
2024-02-09T13:24:54.181755Z  INFO config = Args {
    path_db: "/data/cephfs-1/work/groups/cubi/projects/biotools/mehari/0.21.0/db/",
    genome_release: None,
    path_input_ped: "test.ped",
    path_input_vcf: "test.multi.vcf",
    output: PathOutput {
        path_output_vcf: None,
        path_output_tsv: Some(
            "test.tsv",
        ),
    },
    transcript_source: Both,
    report_all_transcripts: true,
    transcript_picking: false,
    max_var_count: None,
}
2024-02-09T13:24:54.181794Z  INFO Loading HGNC map ...
2024-02-09T13:24:54.216160Z  INFO ... done loading HGNC map
2024-02-09T13:24:54.216501Z  INFO Loading pedigree...
2024-02-09T13:24:54.216736Z  INFO ... done loading pedigree
2024-02-09T13:24:54.216745Z  INFO Open VCF and read header
2024-02-09T13:24:54.220964Z  INFO Determined input assembly to be Grch38
2024-02-09T13:24:54.220973Z  INFO Opening frequency database
2024-02-09T13:24:54.249820Z  INFO Opening ClinVar database
2024-02-09T13:24:54.258946Z  INFO Opening transcript database
2024-02-09T13:24:55.340737Z  INFO Building transcript interval trees ...
2024-02-09T13:24:55.388625Z  INFO ... done building transcript interval trees
2024-02-09T13:24:55.388643Z  INFO Annotating VCF ...
Error: Parsing Error: Error { input: "<DEL>", code: TakeWhile1 }

Expected behavior
Annotation should succeed.

Additional context
N/A

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
Archived in project
Development

Successfully merging a pull request may close this issue.

1 participant