payload: {"action":"closed","issue":{"url":"https://api.github.com/repos/dellytools/delly/issues/154","repository_url":"https://api.github.com/repos/dellytools/delly","labels_url":"https://api.github.com/repos/dellytools/delly/issues/154/labels{/name}","comments_url":"https://api.github.com/repos/dellytools/delly/issues/154/comments","events_url":"https://api.github.com/repos/dellytools/delly/issues/154/events","html_url":"https://github.com/dellytools/delly/issues/154","id":435738023,"node_id":"MDU6SXNzdWU0MzU3MzgwMjM=","number":154,"title":"Delly finds identical deletions with different consensus sequences","user":{"login":"chapplec","id":16101661,"node_id":"MDQ6VXNlcjE2MTAxNjYx","avatar_url":"https://avatars3.githubusercontent.com/u/16101661?v=4","gravatar_id":"","url":"https://api.github.com/users/chapplec","html_url":"https://github.com/chapplec","followers_url":"https://api.github.com/users/chapplec/followers","following_url":"https://api.github.com/users/chapplec/following{/other_user}","gists_url":"https://api.github.com/users/chapplec/gists{/gist_id}","starred_url":"https://api.github.com/users/chapplec/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/chapplec/subscriptions","organizations_url":"https://api.github.com/users/chapplec/orgs","repos_url":"https://api.github.com/users/chapplec/repos","events_url":"https://api.github.com/users/chapplec/events{/privacy}","received_events_url":"https://api.github.com/users/chapplec/received_events","type":"User","site_admin":false},"labels":[{"id":65996285,"node_id":"MDU6TGFiZWw2NTk5NjI4NQ==","url":"https://api.github.com/repos/dellytools/delly/labels/bug","name":"bug","color":"fc2929","default":true}],"state":"closed","locked":false,"assignee":{"login":"tobiasrausch","id":5946458,"node_id":"MDQ6VXNlcjU5NDY0NTg=","avatar_url":"https://avatars1.githubusercontent.com/u/5946458?v=4","gravatar_id":"","url":"https://api.github.com/users/tobiasrausch","html_url":"https://github.com/tobiasrausch","followers_url":"https://api.github.com/users/tobiasrausch/followers","following_url":"https://api.github.com/users/tobiasrausch/following{/other_user}","gists_url":"https://api.github.com/users/tobiasrausch/gists{/gist_id}","starred_url":"https://api.github.com/users/tobiasrausch/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/tobiasrausch/subscriptions","organizations_url":"https://api.github.com/users/tobiasrausch/orgs","repos_url":"https://api.github.com/users/tobiasrausch/repos","events_url":"https://api.github.com/users/tobiasrausch/events{/privacy}","received_events_url":"https://api.github.com/users/tobiasrausch/received_events","type":"User","site_admin":false},"assignees":[{"login":"tobiasrausch","id":5946458,"node_id":"MDQ6VXNlcjU5NDY0NTg=","avatar_url":"https://avatars1.githubusercontent.com/u/5946458?v=4","gravatar_id":"","url":"https://api.github.com/users/tobiasrausch","html_url":"https://github.com/tobiasrausch","followers_url":"https://api.github.com/users/tobiasrausch/followers","following_url":"https://api.github.com/users/tobiasrausch/following{/other_user}","gists_url":"https://api.github.com/users/tobiasrausch/gists{/gist_id}","starred_url":"https://api.github.com/users/tobiasrausch/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/tobiasrausch/subscriptions","organizations_url":"https://api.github.com/users/tobiasrausch/orgs","repos_url":"https://api.github.com/users/tobiasrausch/repos","events_url":"https://api.github.com/users/tobiasrausch/events{/privacy}","received_events_url":"https://api.github.com/users/tobiasrausch/received_events","type":"User","site_admin":false}],"milestone":null,"comments":3,"created_at":"2019-04-22T14:03:42Z","updated_at":"2019-05-02T12:58:48Z","closed_at":"2019-05-02T12:58:47Z","author_association":"CONTRIBUTOR","body":"I am using delly to call structural variants on a human WGS sample. I noticed that it finds certain variants that are essentially duplicates. For example:\r\n\r\n    chr4    177594033       DEL00005339     G       <DEL>   .       PASS    PRECISE;SVTYPE=DEL;SVMETHOD=EMBL.DELLYv0.8.1;CHR2=chr4;END=177594095;PE=0;MAPQ=0;CT=3to5;CIPOS=-3,3;CIEND=-3,3;SRMAPQ=60;INSLEN=0;HOMLEN=3;SR=5;SRQ=0.993243;CONSENSUS=GACAAAAACATATAAACAAATACAAATTATTATTATATGAAGGGTGAAAAATTATGTTCAGCTTTAATATTATGTGGCCTGAAACAAATTTGTATAAGTCATATGATCCGTTTTGTATATTCATTAAAATAGTCATAAAAATATAAAC;CE=1.75466     GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV   0/1:-51.3797,0,-42.791:10000:PASS:1115:1821:1820:1:0:0:23:27\r\n    chr4    177594033       DEL00005340     G       <DEL>   .       PASS    PRECISE;SVTYPE=DEL;SVMETHOD=EMBL.DELLYv0.8.1;CHR2=chr4;END=177594095;PE=0;MAPQ=0;CT=3to5;CIPOS=-3,3;CIEND=-3,3;SRMAPQ=60;INSLEN=0;HOMLEN=3;SR=12;SRQ=0.994186;CONSENSUS=ATATAAACAAATACAAATTATTATTATATGAAGGGTGAAAAATTATGTTCAGCTTTAATATTATGTGGCCTGAAACAAATTTGTATAAGTCATATGATCCGTTTTGTATATTCATTAAAATAGTCATAAAAATATAAACAATGTAATAAGTATGTTTGAATAGTATTATGAT;CE=1.72637    GT:GL:GQ:FT:RCL:RC:RCR:CN:DR:DV:RR:RV   0/1:-51.3797,0,-42.791:10000:PASS:1115:1821:1820:1:0:0:23:27\r\n    \r\nThe two variants above have the same start position and the same end position (`END=177594095`). However, they differ in the `SR`, `SRQ` and `CONSENSUS` fields. These are defined in the VCF header as:\r\n\r\n    ##INFO=<ID=SR,Number=1,Type=Integer,Description=\"Split-read support\">\r\n    ##INFO=<ID=SRQ,Number=1,Type=Float,Description=\"Split-read consensus alignment quality\">\r\n    ##INFO=<ID=CONSENSUS,Number=1,Type=String,Description=\"Split-read consensus sequence\">\r\n    \r\nOf course, the consensus sequences are quite similar, but they aren't identical (I have simply added some spaces to the beginning of the 2nd sequence here to align them):\r\n\r\n    GACAAAAACATATAAACAAATACAAATTATTATTATATGAAGGGTGAAAAATTATGTTCAGCTTTAATATTATGTGGCCTGAAACAAATTTGTATAAGTCATATGATCCGTTTTGTATATTCATTAAAATAGTCATAAAAATATAAAC\r\n             ATATAAACAAATACAAATTATTATTATATGAAGGGTGAAAAATTATGTTCAGCTTTAATATTATGTGGCCTGAAACAAATTTGTATAAGTCATATGATCCGTTTTGTATATTCATTAAAATAGTCATAAAAATATAAACAATGTAATAAGTATGTTTGAATAGTATTATGAT\r\n\r\nSo what am I looking at here, exactly? If this is a deletion on chromosome 4 from position 177594033 to position 177594095, then how can the length of the consensus sequence be different? This is the \"split read consensus\", so presumably the sequence that appears to have been deleted based on how the reads covering that region are split. As I understand it, this means that these reads align to the genome with a gap and this gap indicates the deletion.\r\n\r\nTherefore, if the start and end position are the same, the consensus should also be identical. And, in any case, if the start and end positions are the same, then this should be the same structural variant, so why is it being reported twice? How can I interpret this result?\r\n\r\n  [1]: https://github.com/dellytools/delly/\r\n  [2]: https://www.researchgate.net/figure/Five-approaches-to-detect-CNVs-from-NGS-short-reads-A-Paired-end-mapping-PEM-strategy_fig1_260376412\r\n  [3]: https://i.stack.imgur.com/kodIq.png"}}
other: {"actor":{"display_login":"chapplec"}}