Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,2660AC24.4% 106.1 / 15.5 41H491P (CAC→CCC) recNDNA repair protein RecN
Reads supporting (aligned to +/- strand):  ref base A (30/1);  new base C (0/10);  total (30/11)
Fisher's exact test for biased strand distribution p-value = 9.81e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTG  >  NC_000913/2753212‑2753328
                                                      |                                                              
tCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACaa                                                      >  1:2449654/1‑65 (MQ=255)
   gCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTc                                                   >  1:326266/1‑65 (MQ=255)
   gCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGAGTGTTACCCACCTGCCACAAGTc                                                   >  1:1113501/1‑65 (MQ=255)
        ctgctgCTTCAACGCGGCGAATCAACTCAGGGGTTGTTTTTTCCCCCCCTGCCACAAGTCGCggg                                              <  1:231897/65‑1 (MQ=255)
        ctgctgCCTCAACTCGGCGAATCAACTCAGGTGGTGTGTTTTACCCCCCTGCCACAAGTCGCggg                                              <  1:1961173/65‑1 (MQ=255)
          gctgctTCAACTCGGCGAATCACCTCAGGTGATTTGTTTTCCCCCCCTGCCACAAGTCGCGGGAt                                            <  1:1220538/65‑1 (MQ=255)
               gTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCCCAAGTCGCGGGATGTGGt                                       >  1:1352049/1‑65 (MQ=255)
               gTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGt                                       >  1:2698130/1‑65 (MQ=255)
               gTCAACTCGGAGAATCAACTCAGGTGATGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGt                                       <  1:2744682/65‑1 (MQ=255)
               gTCAACCCGGCAAAGCAACTCAGGTCATGTCTGTTACCCACCTGCCACAAGTCGCGGGATGTGGt                                       <  1:1169755/65‑1 (MQ=255)
                 cagcTCGGCAAAAAAACTCAGGTGGTGTTTGTTACCCGCCTGCCACAAGTCGCGGGATGTGGtca                                     <  1:3011056/62‑1 (MQ=255)
                 cAACGCGGCGAATCAACTCAGGTGTTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtca                                     <  1:1459130/65‑1 (MQ=255)
                    cgcGGCGAATCAACTCGGGTGGTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtcatca                                  <  1:2500900/63‑1 (MQ=255)
                    cgcGGCGAATCAACGCAGGTGAGTTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtcatca                                  <  1:968250/63‑1 (MQ=255)
                      cGGCGAATCAACTCAGGTGATGTGTTTTCCCCACCTGCCACAAGTCGCGGGATGTGGTCATCAac                                <  1:546423/65‑1 (MQ=255)
                        gcaaaTCAGCTCAGGTGTTTTTTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAacac                              <  1:2201408/62‑1 (MQ=255)
                             tCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTAttt                         >  1:2870846/1‑65 (MQ=255)
                                  tCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCa                    >  1:1329235/1‑65 (MQ=255)
                                       tGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCaaa               >  1:2504699/1‑65 (MQ=255)
                                        gATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAg              >  1:2355803/1‑65 (MQ=255)
                                         aTGTGTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGa             <  1:437857/65‑1 (MQ=255)
                                         aTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGa             <  1:1900293/65‑1 (MQ=255)
                                          tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa            >  1:411899/1‑65 (MQ=255)
                                          tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa            >  1:1157999/1‑65 (MQ=255)
                                          tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa            >  1:294053/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:1521880/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:44470/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:3750494/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:1256590/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:1390969/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:702743/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:2749502/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:1902208/1‑65 (MQ=255)
                                            tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc          >  1:2134268/1‑65 (MQ=255)
                                             gtgtTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAcc         <  1:249018/65‑1 (MQ=255)
                                              tgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCg        >  1:1929937/1‑65 (MQ=255)
                                               gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa       >  1:1564072/1‑65 (MQ=255)
                                               gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa       >  1:1762364/1‑65 (MQ=255)
                                               gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa       >  1:644519/1‑65 (MQ=255)
                                               gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa       >  1:2128741/1‑65 (MQ=255)
                                                  aCCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATgg    >  1:849636/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  >  1:1091704/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  >  1:3851169/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  >  1:3600288/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  >  1:3005351/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  >  1:1684948/1‑65 (MQ=255)
                                                    ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg  <  1:1010764/65‑1 (MQ=255)
                                                      |                                                              
TCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTG  >  NC_000913/2753212‑2753328

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