Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927316,7340AC54.5% 19.7 / 28.8 22intergenic (+177/+17)BW25113_RS00080/mokCIS4‑like element IS421 family transposase/type I toxin‑antitoxin system toxin MokC
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/12);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCT  >  NZ_CP009273/16651‑16823
                                                                                   |                                                                                         
gaaCGGGGGGCATCGTCTTAACGAGGCACGGGGGCGTCGCTTTCTTCAGATGTTTCACCCCTTAAGTTAGCGTTTATGGGATcccccccc                                                                                     <  2:983156/90‑1 (MQ=11)
  aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc                                                                                   >  1:867961/1‑90 (MQ=39)
  aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc                                                                                   >  1:627791/1‑90 (MQ=39)
      gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGTTTATGGGATCACTCCCCGCCGtt                                                                               <  1:439659/90‑1 (MQ=255)
            tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCtt                                                                         >  1:969161/1‑90 (MQ=255)
             cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGTTTATGGGATCCCTCCCCGCCGTTGCTCTTa                                                                        <  1:868205/90‑1 (MQ=255)
                        ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGt                                                             >  1:34852/1‑90 (MQ=255)
                         gCACGGAGGCGTCGCATTCTTCAGTTGGTTCACCCTTTAAGTTAGCGTTTATGGGACCCCTCCCCGCCGTTGCTCTTACTCGGATTCGTa                                                            <  1:866360/90‑1 (MQ=21)
                         gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTa                                                            >  1:633318/1‑90 (MQ=255)
                                    tCGCATTCTTCAGGGGGTTCCACCCTTAAGTTGGCCTTTTGGGGACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGaaaa                                                 <  2:136210/90‑1 (MQ=255)
                                                   ggTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCt                                  >  2:534928/1‑90 (MQ=255)
                                                       caccccTTAAGTTGGCGTTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  2:222304/87‑1 (MQ=38)
                                                       caccccTTAAGTTAGCGTTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  2:572389/87‑1 (MQ=38)
                                                       caccccTTAAGTTAGCGTTTAGGGGACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  2:627791/87‑1 (MQ=38)
                                                       caccccTTAAGTTAGCGCTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  1:145855/87‑1 (MQ=39)
                                                       caccccTAAATTTAGCGCTTAGGGGACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  2:34852/87‑1 (MQ=37)
                                                       aaaCCCTTAAGTAAGCCCTTGTGAAAACCCCCCCCGCCGTTGCTCTTCCTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                              <  2:637645/89‑1 (MQ=25)
                                                            cTTATGTTGGCGCTTGTGGACCCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTc                         <  2:798601/90‑1 (MQ=25)
                                                                   tAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGAtgtg                  >  2:922514/1‑90 (MQ=255)
                                                                              ggATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACaga       <  1:909880/90‑1 (MQ=255)
                                                                              ggATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACaga       <  1:760780/90‑1 (MQ=255)
                                                                               gATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACagag      >  2:481974/1‑90 (MQ=255)
                                                                                  cACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTc   >  1:139812/1‑90 (MQ=255)
                                                                                   aCTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCt  >  1:1011577/1‑90 (MQ=255)
                                                                                   |                                                                                         
GAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCT  >  NZ_CP009273/16651‑16823

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