Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,2660AC30.8% 41.4 / 12.4 26H491P (CAC→CCC) recNDNA repair protein RecN
Reads supporting (aligned to +/- strand):  ref base A (15/3);  new base C (0/8);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGCAGTTGTCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGATG  >  NC_000913/2753201‑2753333
                                                                 |                                                                   
aGCGGCAGTTGTCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCAcc                                                                   >  1:560852/1‑68 (MQ=255)
aGCGGCAGTTGTCGGCAAACTGCTGCGTCAACTCGGCGAATCAACGCAGGTGATGTGTGTTACCCAcc                                                                   <  1:1940889/68‑1 (MQ=255)
aGCGGCAGTTGTCGGCAAACTGCTGCGTCAACGCGGCGAATAAACTCAGGTGATGTGTGTTACCCAcc                                                                   <  1:608657/68‑1 (MQ=255)
         tGTCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAg                                                          >  1:1466931/1‑68 (MQ=255)
            cGGAAAACGGCGGCGTCAACCCGGGGAATCAACTAGGGTGATGTGTTTTCCCCCCCTGCCACAAGTcg                                                       <  1:2412303/68‑1 (MQ=255)
                cagctgctcCGTCAACTCGGCAAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACGAGTCGCggg                                                   <  1:627927/65‑1 (MQ=255)
                    tgctgcTTCAACCCGGCGAATCAACTCAGGTGATGTGTGTTCCCCCCCTGCCACAAGTCGCNGGAtgt                                               <  1:2044346/68‑1 (MQ=255)
                    tgcggcGTCAACCCGGCGAACCACCTCAGGTGATGTGTTTTCCCCCCCTGCCACAAGTCGCGGGAtgt                                               <  1:760393/68‑1 (MQ=255)
                       tgcGTCAACGCGGCGAATCAACTCAGGTGATGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGt                                            <  1:2818090/68‑1 (MQ=255)
                         ctccAACTCAGCAAAACACCGCGGGTGATTTGTGTGCCCCCCCTGCCACAAGTCGCGGGCTGTGGtca                                          <  1:1702684/65‑1 (MQ=255)
                         cGTCACCTGGGCGAATCACCTCAGGTGATGTGTGTTCCCCCCCTGCCACAAGTCGTGGGATGTGGtca                                          <  1:2149475/68‑1 (MQ=255)
                          gTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGtcat                                         >  1:2589282/1‑68 (MQ=255)
                              acccGGCGAATCAACTCAGGTGGTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAac                                     <  1:2191420/65‑1 (MQ=255)
                              aCTCGGCGAATCAACTCAGGTGATGTGTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAac                                     <  1:1709958/68‑1 (MQ=255)
                                  ggCAAATCATCTCAGTTGTTGTGTGTGCCCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTa                                 <  1:2816266/68‑1 (MQ=255)
                                    cGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGCCATCAACACTAtt                               >  1:1745154/1‑68 (MQ=255)
                                     gAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTAttt                              <  1:2945115/68‑1 (MQ=255)
                                        tCAACTCAGGTGATGCGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGt                           >  1:1818809/1‑68 (MQ=255)
                                             tCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCa                      >  1:1785643/1‑68 (MQ=255)
                                               aGGTGATGTGTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCaaa                    <  1:109779/68‑1 (MQ=255)
                                                 gTGATGTGTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGa                  <  1:1767803/68‑1 (MQ=255)
                                                     tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAcc              >  1:2522153/1‑68 (MQ=255)
                                                     tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAcc              >  1:450177/1‑68 (MQ=255)
                                                       tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa            >  1:2817003/1‑68 (MQ=255)
                                                       tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa            >  1:3101847/1‑68 (MQ=255)
                                                       tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa            >  1:1593090/1‑68 (MQ=255)
                                                       tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa            >  1:871765/1‑68 (MQ=255)
                                                       tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGCCAGCAAAGAAACCGa            >  1:2145991/1‑68 (MQ=255)
                                                          gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATgg         >  1:34194/1‑68 (MQ=255)
                                                            ccccccccTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg       <  1:943877/66‑1 (MQ=255)
                                                               ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGa    >  1:1406667/1‑68 (MQ=255)
                                                               ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGa    >  1:673423/1‑68 (MQ=255)
                                                               ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGa    >  1:1042820/1‑68 (MQ=255)
                                                                 aCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGATg  >  1:1332611/1‑68 (MQ=255)
                                                                 |                                                                   
AGCGGCAGTTGTCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTGCGATG  >  NC_000913/2753201‑2753333

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: ‑