I want to add PARID into the INFO field of an SV vcf file with breakend notation
> StructuralVariantAnnotation:::breakpointGRangesToVCF(
+ breakpointRanges(
+ VariantAnnotation::readVcf('SV.delly.somatic.vcf'),
+ inferMissingBreakends = TRUE
+ )
+ )
Error: Elements 1, 2, 3, 4, 5, ... of width(gr) == 1 are not true