Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,280,7040AC52.4% ‑2.4 / 35.4 21T25P (ACC→CCC) radDputative DNA repair helicase RadD
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/11);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGG  >  NC_000913/2280642‑2280764
                                                              |                                                            
cacTACGCCCATATCAGAAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc                                                            <  1:1180829/65‑1 (MQ=255)
 acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc                                                           >  1:1396649/1‑65 (MQ=255)
         cagatCAACAAGACGCCGGGGATGCCACGCTCAACCATTTGCGTCGTCAAAAACCCCCTGCCGtt                                                   <  1:755076/62‑1 (MQ=255)
         aatagCAGCAAGAAGCTGTGGATGCCACGCTAAACCTTTTGCGTCGTCAAAAACCCCCTGCCGtt                                                   <  1:1012549/64‑1 (MQ=255)
              cAGAAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt                                              <  1:415135/65‑1 (MQ=255)
                 cAAGAAGCCGTGGAGGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgct                                           <  1:2668438/65‑1 (MQ=255)
                  aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg                                          >  1:954096/1‑65 (MQ=255)
                   agaagCCGTGGATGCCCCGCTCACCCATTTCCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                         <  1:3695775/65‑1 (MQ=255)
                   agaagCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                         <  1:3549967/65‑1 (MQ=255)
                   aaaagCCGTGGATCCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                         <  1:991188/63‑1 (MQ=255)
                     gagCCGTGAACGCGCCCCTCAGCCTTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGccc                                       <  1:3415312/64‑1 (MQ=255)
                       gCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCAc                                     >  1:3819189/1‑65 (MQ=255)
                         cGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCg                                   <  1:339178/65‑1 (MQ=255)
                           tGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGgc                                 <  1:3865031/65‑1 (MQ=255)
                           tGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgc                                 >  1:3298716/1‑65 (MQ=255)
                              aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                              >  1:2979729/1‑65 (MQ=255)
                              aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                              >  1:2846995/1‑65 (MQ=255)
                              aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                              >  1:1870868/1‑65 (MQ=255)
                              aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa                              >  1:3734852/1‑65 (MQ=255)
                                gCCACGCTCACCCATTTTCGTCTTAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg                            <  1:3063539/65‑1 (MQ=255)
                                 ccgcGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                           <  1:3769656/62‑1 (MQ=255)
                                 ccACCCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                           <  1:1997086/65‑1 (MQ=255)
                                      gtCAACCAGTTTCGCCGCCAAAAGCCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                      <  1:2818621/64‑1 (MQ=255)
                                      cTCATCCATGTGCCTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                      <  1:2910041/65‑1 (MQ=255)
                                      cTCACCCATTTTCCTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                      <  1:2448065/65‑1 (MQ=255)
                                      cTCAACTATCTTCGCCGTCATAACCCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                      <  1:2381947/65‑1 (MQ=255)
                                       tCACCCTTTTCCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGc                     <  1:2228880/65‑1 (MQ=255)
                                       ccAACTATTTCCGTCGTAAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGc                     <  1:1362781/64‑1 (MQ=255)
                                             aTTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTg               >  1:195256/1‑65 (MQ=255)
                                                  cgtcgtCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGc          <  1:3416488/65‑1 (MQ=255)
                                                     cgtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa       <  1:129493/65‑1 (MQ=255)
                                                          aaaacCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTgg  <  1:603456/64‑1 (MQ=255)
                                                              |                                                            
CACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGG  >  NC_000913/2280642‑2280764

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