public static void WriteMultiSampleSegments(string outVcfPath, ISampleMap <List <CanvasSegment> > segments, List <double> diploidCoverage,
string wholeGenomeFastaDirectory, List <string> sampleNames, List <string> extraHeaders, List <PloidyInfo> ploidies, int qualityThreshold, int?denovoQualityThreshold, int?sizeThreshold, bool isPedigreeInfoSupplied = true)
{
using (BgzipOrStreamWriter writer = new BgzipOrStreamWriter(outVcfPath))
{
var genome = WriteVcfHeader(segments.Values.First(), diploidCoverage.Average(), wholeGenomeFastaDirectory, sampleNames,
extraHeaders, writer, qualityThreshold, denovoQualityThreshold, sizeThreshold);
WriteVariants(segments.Zip(), ploidies, genome, writer, denovoQualityThreshold);
}
}
private static IEnumerable <ISampleMap <OverlappingSegmentsRegion> > GetOverlappingSegmentsRegionSampleLists(ISampleMap <List <OverlappingSegmentsRegion> > sampleRegions)
{
return(sampleRegions.Zip());
}