diff --git a/gridss.wdl b/gridss.wdl index 9a09bddee130edf4b1c9545a6b30e8f18d848535..8e1474c1024e114e2c81be1f14d66149606d1b9e 100644 --- a/gridss.wdl +++ b/gridss.wdl @@ -120,7 +120,7 @@ task AnnotateSvTypes { svtype <- simpleEventType(gr) info(vcf[gr$sourceId])$SVTYPE <- svtype # GRIDSS doesn't supply a GT, simply set it to 0/1 - geno(vcf)$GT <- "0/1" + geno(vcf)$GT <- as.matrix(sapply(row.names(vcf), function(x) {"0/1"})) # Select only one breakend per event (also removes single breakends): # sourceId ends with o or h for paired breakends, the first in the pair # end with o the second with h. Single breakend end with b, these will