Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,419:1 +T 42.9% intergenic (+23/+163) 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,4191.T42.9% 12.6 / 28.6 14intergenic (+23/+163)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base . (4/4);  new base T (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

CGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAA  >  NC_000913/3330337‑3330475
                                                                                   |                                                        
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaa                                  >  3:390993/1‑107 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaa                                  <  4:390993/107‑1 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCgggg                                                   <  1:393341/90‑1 (MQ=255)
cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCgggg                                                   >  2:393341/1‑90 (MQ=255)
         ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAg                                <  3:565336/100‑1 (MQ=255)
         ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAg                                >  4:565336/1‑100 (MQ=255)
                        aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                        <  1:973244/77‑1 (MQ=255)
                        aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                        >  2:973244/1‑77 (MQ=255)
                                     aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                        <  3:1038465/64‑1 (MQ=255)
                                     aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                        >  4:1038465/1‑64 (MQ=255)
                                                                     agccccGGCCATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGt                      <  3:367888/48‑1 (MQ=17)
                                                                     agccccGGCCATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGt                      >  4:367888/3‑50 (MQ=17)
                                                                      gccccGGCCATGTGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa  <  1:969424/68‑1 (MQ=17)
                                                                      gccccGGCCATGTGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa  >  2:969424/2‑69 (MQ=17)
                                                                                   |                                                        
CGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAA  >  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: