Skip to content

convert --gvcf2vcf adds one to end of block #242

Description

@mwonge

Current --gvcf2vcf adds an extra position to the end of block.

Example below: End block is position 12782 but after conversion, end block is 12783 which then clashes with the real 12783 and you get duplicated locations.

gVCF (truncated to end block) :
chr1 12768 . C . . PASS END=12782;BLOCKAVG_min30p3a
chr1 12783 . G T 63 PASS SNVSB=0;SNVHPOL=2
chr1 12784 . T . . PASS END=12800;BLOCKAVG_min30p3a

After conversion:
chr1 12768 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12769 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12770 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12771 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12772 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12773 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12774 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12775 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12776 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12777 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12778 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12779 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12780 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12781 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12782 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF 0/0:33:12:1
chr1 12783 . C . . PASS BLOCKAVG_min30p3a GT:GQX:DP:DPF
chr1 12783 . G T 63 PASS SNVSB=0;SNVHPOL=2 GT:GQ:GQX:DP:DPF:AD

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