/// <summary>
/// Generates the Scaffold.
/// </summary>
public void GenerateScaffold()
{
Output.WriteLine(OutputLevel.Information, Resources.ScaffoldStarting);
if (this.FileNames.Length != 2)
{
Output.WriteLine(OutputLevel.Error, "\nError: A reference file and 1 query file are required.");
Output.WriteLine(OutputLevel.Required, Resources.ScaffoldHelp);
return;
}
if (!string.IsNullOrEmpty(this.CloneLibraryName))
{
CloneLibrary.Instance.AddLibrary(this.CloneLibraryName, (float)this.MeanLengthOfInsert, (float)this.StandardDeviationOfInsert);
}
TimeSpan algorithmSpan = new TimeSpan();
Stopwatch runAlgorithm = new Stopwatch();
FileInfo refFileinfo = null;
using (GraphScaffoldBuilder scaffoldBuilder = new GraphScaffoldBuilder())
{
refFileinfo = new FileInfo(this.FileNames[0]);
long refFileLength = refFileinfo.Length;
runAlgorithm.Restart();
IEnumerable <ISequence> contigs = AssembleArguments.ParseFile(this.FileNames[0]);
runAlgorithm.Stop();
algorithmSpan = algorithmSpan.Add(runAlgorithm.Elapsed);
if (this.Verbose)
{
Output.WriteLine(OutputLevel.Verbose);
Output.WriteLine(OutputLevel.Verbose, "Processed contigs file : {0}", Path.GetFullPath(this.FileNames[0]));
Output.WriteLine(OutputLevel.Verbose, " Read/Processing time: {0}", runAlgorithm.Elapsed);
Output.WriteLine(OutputLevel.Verbose, " File Size : {0}", refFileLength);
Output.WriteLine(OutputLevel.Verbose, " k-mer Length : {0}", this.KmerLength);
}
refFileinfo = new FileInfo(this.FileNames[1]);
refFileLength = refFileinfo.Length;
runAlgorithm.Restart();
IEnumerable <ISequence> reads = AssembleArguments.ParseFile(this.FileNames[1]);
runAlgorithm.Stop();
algorithmSpan = algorithmSpan.Add(runAlgorithm.Elapsed);
if (this.Verbose)
{
Output.WriteLine(OutputLevel.Verbose);
Output.WriteLine(OutputLevel.Verbose, "Processed reads file : {0}", Path.GetFullPath(this.FileNames[1]));
Output.WriteLine(OutputLevel.Verbose, " Read/Processing time: {0}", runAlgorithm.Elapsed);
Output.WriteLine(OutputLevel.Verbose, " File Size : {0}", refFileLength);
}
runAlgorithm.Restart();
ValidateAmbiguousReads(reads);
runAlgorithm.Stop();
if (this.Verbose)
{
Output.WriteLine(OutputLevel.Verbose);
Output.WriteLine(OutputLevel.Verbose, "Time taken for Validating reads: {0}", runAlgorithm.Elapsed);
}
runAlgorithm.Restart();
IList <ISequence> scaffolds = scaffoldBuilder.BuildScaffold(reads, contigs.ToList(), this.KmerLength, this.Depth, this.Redundancy);
runAlgorithm.Stop();
algorithmSpan = algorithmSpan.Add(runAlgorithm.Elapsed);
if (this.Verbose)
{
Output.WriteLine(OutputLevel.Verbose);
Output.WriteLine(OutputLevel.Verbose, "Compute time: {0}", runAlgorithm.Elapsed);
}
runAlgorithm.Restart();
WriteContigs(scaffolds);
runAlgorithm.Stop();
algorithmSpan = algorithmSpan.Add(runAlgorithm.Elapsed);
if (this.Verbose)
{
Output.WriteLine(OutputLevel.Verbose);
Output.WriteLine(OutputLevel.Verbose, "Write contigs time: {0}", runAlgorithm.Elapsed);
Output.WriteLine(OutputLevel.Verbose, "Total runtime: {0}", algorithmSpan);
}
}
}
