Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,280,7040AC41.4% 38.1 / 13.1 29T25P (ACC→CCC) radDputative DNA repair helicase RadD
Reads supporting (aligned to +/- strand):  ref base A (12/5);  new base C (0/12);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.23e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAAC  >  NC_000913/2280638‑2280761
                                                                  |                                                         
tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAc                                                          <  1:2959642/68‑1 (MQ=255)
 ttACCCTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc                                                         <  1:1815780/68‑1 (MQ=255)
 ttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc                                                         <  1:3027321/68‑1 (MQ=255)
 ttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc                                                         <  1:2803133/68‑1 (MQ=255)
 cagcacTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCACAAAAcc                                                         <  1:2055389/65‑1 (MQ=255)
  tACACTTCGCCCATATCAGCAAGAAGCCGCGGATGCCACGCTCAACCATTTTCGTCGGCATAGAAccc                                                        >  1:1357302/1‑68 (MQ=255)
     acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTg                                                     >  1:1414765/1‑68 (MQ=255)
     acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTg                                                     >  1:3161953/1‑68 (MQ=255)
           cccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTa                                               >  1:2098436/1‑68 (MQ=255)
               taccAGAAAGAAGCGGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt                                           <  1:2701094/65‑1 (MQ=255)
               atacAGAAAGAAGCCGGGGATGCCACGCTCAACCATTTTCGTCGTAAAAAACCCCCTGCCGTTATCGt                                           <  1:3359663/65‑1 (MQ=255)
                 acAGAAAAAAGCCGTGTATGCCACGCTCAACCATTTGCGTCGTCAAAAACCCCCTGCCGTTATCGtgc                                         <  1:1864973/67‑1 (MQ=255)
                  cAGAAAGAAGCCGTGGATGCCACTCTCAACCATTTGCGTCGTCAAAAACCCCCTGCCGTTATCGtgct                                        <  1:1405626/68‑1 (MQ=255)
                    gaaagaggCCGTGGATGCCGCGCTCACCCTTTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                      <  1:2963894/66‑1 (MQ=255)
                    gCAAGAAGCCGTGGATGCCCCGCTCACCCATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGtgctgc                                      <  1:2585748/68‑1 (MQ=255)
                    gCAAAAAGCCGTGGATGCCCCGCTCACCCATTTTCGTCGGCAAAAACCCCCTGCCGTTATCGtgctgc                                      <  1:986695/68‑1 (MQ=255)
                    acAAGAGGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc                                      <  1:3456190/67‑1 (MQ=255)
                      aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGccc                                    >  1:779009/1‑68 (MQ=255)
                      aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGccc                                    >  1:501543/1‑68 (MQ=255)
                           gCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCgg                               >  1:973215/1‑68 (MQ=255)
                           gCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCgg                               >  1:482107/1‑68 (MQ=255)
                                 gctGCCCCGCTCAACCATTTTGGTCTTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg                         <  1:987232/66‑1 (MQ=255)
                                 gATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg                         <  1:1218389/68‑1 (MQ=255)
                                 aaTGCCCCCCTCACCCATTTTCGTCGTAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg                         <  1:2005434/67‑1 (MQ=255)
                                  ttGCCAGGCGAAACCATTCTGGTCTTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                        <  1:380427/67‑1 (MQ=255)
                                  atcccACGCTCACCCTTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                        <  1:2934105/65‑1 (MQ=255)
                                  aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                        >  1:2070425/1‑68 (MQ=255)
                                  aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                        >  1:61058/1‑68 (MQ=255)
                                  aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt                        >  1:1709304/1‑68 (MQ=255)
                                       gggCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                   <  1:520098/66‑1 (MQ=255)
                                       gctcTCACCCATTTTCGTCGTCATAAAGCCCCGGCCGTTATCGTGTTGCCCACCGGCGCAGGTAAAAg                   <  1:2309757/65‑1 (MQ=255)
                                       ccccTCACCCATTTTCTTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                   <  1:3096422/65‑1 (MQ=255)
                                       ccGCTCACCCTTTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg                   <  1:2812919/67‑1 (MQ=255)
                                       aCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaagg                   <  1:1664141/68‑3 (MQ=255)
                                         gCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGcc                 <  1:3105417/68‑1 (MQ=255)
                                              aCCTTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCGCCGGCGCAGGTAAAAGCCTGGTg            <  1:782021/68‑1 (MQ=255)
                                              aCCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTg            >  1:2613451/1‑68 (MQ=255)
                                                   tttCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGc       >  1:853037/1‑68 (MQ=255)
                                                   tttCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGc       >  1:2942699/1‑68 (MQ=255)
                                                      cgtcgtCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa    <  1:276132/68‑1 (MQ=255)
                                                        tcgtcATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAAc  >  1:78376/1‑68 (MQ=255)
                                                                  |                                                         
TTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAAC  >  NC_000913/2280638‑2280761

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