Skip to content
Snippets Groups Projects
Unverified Commit c4b13ec6 authored by Cats's avatar Cats Committed by GitHub
Browse files

Merge pull request #302 from biowdl/gridss_annotate

Add GT to gridss results in AnnotateSvTypes
parents 1c1cd589 c91a657c
No related branches found
No related tags found
No related merge requests found
......@@ -119,6 +119,8 @@ task AnnotateSvTypes {
gr <- breakpointRanges(vcf)
svtype <- simpleEventType(gr)
info(vcf[gr$sourceId])$SVTYPE <- svtype
# GRIDSS doesn't supply a GT, so we estimate GT based on AF (assuming CN of 2, might be inaccurate)
geno(vcf)$GT <- ifelse(geno(vcf)$AF > 0.75, "1/1", ifelse(geno(vcf)$AF < 0.25, "0/0", "0/1"))
writeVcf(vcf, out_path, index=~{index})
EOF
>>>
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment