gdtools_subtract

gdtools_subtract creates a refined Genome Diff file by excluding mutations present in one or more subtracted Genome Diff files while preserving evidence and validation entries for downstream genomic analysis.

Key Features:

• Genome Diff refinement: Produces an output Genome Diff that retains only entries present in the input file after excluding mutations found in subtracted Genome Diff files.
• Preservation of evidence: All evidence and validation entries from the input are maintained in the output file.
• Integration with breseq: Operates as part of the breseq computational framework for analyzing DNA resequencing data from haploid microbial genomes.
• Structural variation detection: Identifies structural variations using split-read alignments to a reference genome, including events involving repeat sequences.
• Statistical validation: Applies a statistical model to evaluate read coverage evenness to support mutation predictions.
• Biological description of mutations: Combines predictions of new sequence junctions and deleted chromosomal regions to provide biologically relevant mutation descriptions.
• Read-depth applicability: Capable of predicting structural variations with modest read-depth coverage (over 40-fold).

Scientific Applications:

• Microbial epidemiology: Identifies structural variations that contribute to microbial evolution and pathogenic traits.
• Experimental evolution studies: Detects spontaneous mutations such as transposon insertions and chromosomal changes in evolution experiments.
• Synthetic biology and genetics: Supports detection of mutations responsible for significant genomic changes relevant to genome engineering and gene-function studies.

Methodology:

Creates a refined Genome Diff by excluding mutations found in subtracted Genome Diff files while retaining evidence entries; leverages the breseq pipeline to analyze DNA resequencing data; evaluates support for new sequence junctions via split-read alignments to a reference genome (including repeat sequences); employs a statistical model to validate read coverage evenness; predicts structural variations with modest read-depth coverage (over 40-fold).