svtyper icon indicating copy to clipboard operation
svtyper copied to clipboard

KeyError: 'CIPOS'

Open verne91 opened this issue 5 years ago • 1 comments

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!

verne91 avatar Apr 29 '19 14:04 verne91

#111 fixes the issue for me.

dmckean avatar Feb 04 '21 17:02 dmckean