Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,280,7040AC62.1% 3.5 / 27.3 29T25P (ACC→CCC) radDputative DNA repair helicase RadD
Reads supporting (aligned to +/- strand):  ref base A (10/1);  new base C (0/18);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 5.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG  >  NC_000913/2280640‑2280768
                                                                |                                                                
tACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATaaaa                                                                  >  1:1360153/1‑65 (MQ=255)
  cccTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc                                                                <  1:3062641/63‑1 (MQ=255)
   acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc                                                               >  1:1138980/1‑65 (MQ=255)
        ccccATATCAGCAAGAGGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGcc                                                          <  1:2532382/64‑1 (MQ=255)
         cccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg                                                         >  1:1852610/1‑65 (MQ=255)
          ccAAATCAGAAAAAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGt                                                        <  1:2435765/65‑1 (MQ=255)
             tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTAt                                                     >  1:2620133/1‑65 (MQ=255)
                cAGAAAGAAGCGGGGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt                                                  <  1:2671036/65‑1 (MQ=255)
                cAGAAAAAAGCCGTGGATGCCCCGCTCACCCATTTCCGTCGTCATAAACCCCCTGCCGTTATCGt                                                  <  1:2635877/65‑1 (MQ=255)
                aaGCAAAAAGCCGTGGATGCCACGCCCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt                                                  <  1:3767906/64‑1 (MQ=255)
                   cAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgct                                               <  1:2726054/65‑1 (MQ=255)
                   aaagaggCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgct                                               <  1:2914152/64‑1 (MQ=255)
                    aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTCTCgtggtg                                              >  1:3063412/1‑65 (MQ=255)
                    aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg                                              >  1:592115/1‑65 (MQ=255)
                    aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg                                              >  1:3707277/1‑65 (MQ=255)
                     ggaagCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                             <  1:2650078/64‑1 (MQ=255)
                     agaagCCGTGGATCCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                             <  1:2379903/65‑1 (MQ=255)
                     aaaagCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                             <  1:3360789/63‑1 (MQ=255)
                             tGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTNATCGTGCTGCCCACCGgc                                     <  1:2785398/65‑1 (MQ=255)
                                aTGCCACGCTCAACCATTTTCGTCGTTATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                                  >  1:3331811/1‑65 (MQ=255)
                                aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                                  >  1:1136812/1‑65 (MQ=255)
                                aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                                  >  1:3405780/1‑65 (MQ=255)
                                  gCCACGCTCAACCATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGNAGCAgg                                <  1:2204918/65‑1 (MQ=255)
                                   ccACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                               <  1:2354708/65‑1 (MQ=255)
                                     aCGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTac                             >  1:3338281/1‑64 (MQ=255)
                                        cccAACCATTTTCGCCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:2458167/63‑1 (MQ=255)
                                        cTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:1893587/65‑1 (MQ=255)
                                        cTCAACCTTTTTCTGCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:3043073/65‑1 (MQ=255)
                                        cTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          >  1:1513104/1‑65 (MQ=255)
                                        cTCAACCATTTTCGGCGTAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:2736584/65‑1 (MQ=255)
                                        cTCAACCATTTGCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:445961/65‑1 (MQ=255)
                                        accAACCATTTTCTTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                          <  1:2681428/63‑1 (MQ=255)
                                             cctttttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTgg                     <  1:1172453/62‑1 (MQ=255)
                                              ctttttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGt                    <  1:2619437/63‑1 (MQ=255)
                                                       cgtcATAACCCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa           <  1:3291029/65‑1 (MQ=255)
                                                       cgtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa           <  1:2118573/65‑1 (MQ=255)
                                                            tAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTgg      <  1:2264507/65‑1 (MQ=255)
                                                            tAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGATGATAGCGGAACTgg      <  1:3198621/65‑1 (MQ=255)
                                                                aCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACg  >  1:1090439/1‑65 (MQ=255)
                                                                |                                                                
TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG  >  NC_000913/2280640‑2280768

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: