Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,415 Δ1 bp 42.9% intergenic (+19/+167) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,330,4150A.42.9% 26.1 / 21.6 14intergenic (+19/+167)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base . (3/3);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAA  >  NC_000913/3330337‑3330475
                                                                              |                                                             
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACaa                                  >  3:390993/1‑107 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACaa                                  <  4:390993/107‑1 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                   <  1:393341/90‑1 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                   >  2:393341/1‑90 (MQ=255)
         ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAg                                <  3:565336/100‑1 (MQ=255)
         ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAg                                >  4:565336/1‑100 (MQ=255)
                        aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                        <  1:973244/77‑1 (MQ=255)
                        aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                        >  2:973244/1‑77 (MQ=255)
                                     aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                        <  3:1038465/64‑1 (MQ=255)
                                     aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                        >  4:1038465/1‑64 (MQ=255)
                                                                     agccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGt                      <  3:367888/48‑1 (MQ=17)
                                                                     agccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGt                      >  4:367888/3‑50 (MQ=17)
                                                                      gccccGGC‑CATGTGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa  <  1:969424/68‑1 (MQ=17)
                                                                      gccccGGC‑CATGTGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa  >  2:969424/2‑69 (MQ=17)
                                                                              |                                                             
CGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAA  >  NC_000913/3330337‑3330475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑