internal static IWittyerVariant Create([NotNull] IVcfVariant baseVariant, string sample, double percentageDistance, uint basepairDistance, IReadOnlyList <uint> bins, WittyerVariantType svType) { if (svType == WittyerVariantType.Invalid) { throw new InvalidDataException( $"Invalid {VcfConstants.SvTypeKey} in variant: \n{baseVariant}\nNot sure why you got here though. Check with a witty.er developer!"); } var end = baseVariant.Position + baseVariant.GetSvLength(); var baseInterval = BedInterval.Create(baseVariant.Position, end); var borderInterval = baseVariant.Position.CalculateBorderInterval(baseInterval, baseVariant.ParseCi(WittyerConstants.Cipos), percentageDistance, basepairDistance); // wittyerVariant should all have end border, it's a matter of how to find it, // either END key in INFO field, sort out through SVLEN or other ways, details can be defined in FindEndBorder() later var endInterval = (end - 1).CalculateBorderInterval(baseInterval, baseVariant.ParseCi(WittyerConstants.Ciend), percentageDistance, basepairDistance); var posContigAndInterval = ContigAndInterval.Create(baseVariant.Contig, borderInterval.Start, borderInterval.Stop + 1); var endContigAndInterval = ContigAndInterval.Create(baseVariant.Contig, endInterval.Start, endInterval.Stop + 1); return(Create(baseVariant, baseInterval, svType, posContigAndInterval, Winner.Create(svType, baseInterval, bins), new List <OverlapAnnotation>(), WittyerSample.CreateOverall(baseVariant, sample, svType == WittyerVariantType.CopyNumberReference), endContigAndInterval)); }