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

KeyError: 'CIPOS' #103

Open
verne91 opened this issue Apr 29, 2019 · 1 comment
Open

KeyError: 'CIPOS' #103

verne91 opened this issue Apr 29, 2019 · 1 comment

Comments

@verne91
Copy link

verne91 commented Apr 29, 2019

Traceback (most recent call last):
  File "/home/cnsun/.local/bin/svtyper-sso", line 11, in <module>
    load_entry_point('svtyper==0.7.0', 'console_scripts', 'svtyper-sso')()
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/singlesample.py", line 848, in cli
    sys.exit(main())
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/singlesample.py", line 841, in main
    args.batch_size)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/singlesample.py", line 811, in sso_genotype
    genotype_parallel(src_vcf, vcf_out, sample, z, split_slop, min_aligned, sum_quals, split_weight, disc_weight, max_reads, max_ci_dist, debug, cores, batch_size, ref_fasta)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/singlesample.py", line 718, in genotype_parallel
    breakpoints = collect_breakpoints(src_vcf, max_ci_dist)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/singlesample.py", line 134, in collect_breakpoints
    brkpts = vcf.get_variant_breakpoints(variant, max_ci_dist)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/parsers.py", line 221, in get_variant_breakpoints
    breakpoints = self._default_get_breakpoints(variant, max_ci_dist)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/parsers.py", line 178, in _default_get_breakpoints
    ciA = confidence_interval(variant, 'CIPOS', 'CIPOS95', max_ci_dist)
  File "/home/cnsun/.local/lib/python2.7/site-packages/svtyper/parsers.py", line 12, in confidence_interval
    ci = map(int, var.info[tag].split(','))
KeyError: 'CIPOS'

I don't have 'CIPOS' field in my vcf file. How can I genotype my sample with this vcf? Thanks!

@dmckean
Copy link

dmckean commented Feb 4, 2021

#111 fixes the issue for me.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants