Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,612,0670AG43.3% 88.5 / 14.6 60E161G (GAG→GGG) yjjJtype II toxin‑antitoxin system HipA family toxin YjjJ
Reads supporting (aligned to +/- strand):  ref base A (9/25);  new base G (26/0);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.04e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT  >  NZ_CP009273/4611989‑4612155
                                                                              |                                                                                         
gTCTCTGGCAGGAAGAACATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAAATATc                                                                                >  1:1020887/1‑90 (MQ=255)
 tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCa                                                                               <  2:345587/90‑1 (MQ=255)
  ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAg                                                                              >  1:178671/1‑90 (MQ=255)
  ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAg                                                                              >  1:349631/1‑90 (MQ=255)
   tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGaa                                                                                     >  1:411445/1‑82 (MQ=255)
   tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGaa                                                                                     <  2:411445/82‑1 (MQ=255)
   tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGATGGTTGGTCGGGGAGGGGAATTATCAGc                                                                             <  2:780459/90‑1 (MQ=255)
    ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCg                                                                            >  1:803622/1‑90 (MQ=255)
    ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGAATATTAGCg                                                                            >  2:441286/1‑90 (MQ=255)
     tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa                                                                           <  2:1047182/90‑1 (MQ=255)
        gaGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg                                                                        <  2:1490802/89‑1 (MQ=255)
        cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg                                                                        <  2:871724/90‑1 (MQ=255)
           gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGAtt                                                                     <  1:289359/90‑1 (MQ=255)
             agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                   <  1:677851/90‑1 (MQ=255)
             agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                   <  2:178671/90‑1 (MQ=255)
             agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                   <  2:909064/90‑1 (MQ=255)
              gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGGATTATTAGCGATGGATTTCt                                                                  >  2:984387/1‑90 (MQ=255)
                  aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGcac                                                              <  1:1061773/90‑1 (MQ=255)
                      gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATTAGCGATGGATTAATGCacaaaa                                                          >  2:1210192/1‑88 (MQ=255)
                      gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTTTCAGCGATGGGTTACTGCacaaca                                                          >  1:933736/1‑90 (MQ=255)
                             gCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCg                                                  <  1:275948/90‑1 (MQ=255)
                                cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGGTGGATTACTGCACAACACCCTGCggga                                               >  2:989637/1‑88 (MQ=255)
                                cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAAAACCCTGCggga                                               >  1:1267134/1‑88 (MQ=255)
                                cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGCggga                                               >  1:632870/1‑88 (MQ=255)
                                cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAACACCCCGCggga                                               >  2:339685/1‑88 (MQ=255)
                                cTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGGGGGATTACTGCACAAAACCCCGCGGaa                                               >  1:644563/1‑90 (MQ=255)
                                 tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag                                                                                          >  1:174957/1‑47 (MQ=255)
                                 tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCggggag                                                                                          <  2:174957/47‑1 (MQ=255)
                                 tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGgaa                                              >  2:1082582/1‑90 (MQ=255)
                                 tGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCTGCGGaaa                                              >  1:173716/1‑90 (MQ=255)
                                  gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAt                                             <  1:892408/90‑1 (MQ=255)
                                     cGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAAATTc                                           >  2:508698/1‑90 (MQ=255)
                                     cGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGAATTATTAGCGATGGGTTACTGCACAACACCCCGCGGGAAATTc                                           >  2:117672/1‑90 (MQ=255)
                                        aTTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTTTAAGGGAGGGATTATTGCCAAACACCCTCCGGAAATTCctct                                       >  2:113939/1‑90 (MQ=255)
                                         tttAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg                                      <  2:803622/90‑1 (MQ=255)
                                              cGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCa                                 <  1:1249960/90‑1 (MQ=255)
                                                gTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGGAATTTCTCTGGGTCaaa                               >  2:209493/1‑90 (MQ=255)
                                                gTGAATATACTGGCGGCTGGGTTGTCGGGGGGGGGAATTTTCAGCGATGGGTTAATGCAAAACACCCTCCGGGAAATACTCTGGGTTaaa                               >  2:1329633/1‑90 (MQ=255)
                                                 tGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCCCAACACCCCGCGGGAAATTCTCTGGGTCaaaa                              >  2:317513/1‑90 (MQ=255)
                                                  gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAATTCCTCTGGGTCaaaaa                             >  1:796609/1‑90 (MQ=255)
                                                  gAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAAATATCAGCGATGGATTACTGCACAACACCCCGCGGGAATTTCTCTGGGTCaaaaa                             >  2:148009/1‑90 (MQ=255)
                                                  gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                   <  2:408194/53‑1 (MQ=255)
                                                  gAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc                                                                   >  1:408194/1‑53 (MQ=255)
                                                  gAATATACTGGCGGTTGGGTTGTCGGGGGGGGGGATTTTCAGCGATGGGTTTCTGCACAAAACCCTGCGGGAATTTCTCTGGATCaaaaa                             >  1:1087142/1‑90 (MQ=255)
                                                  gAATATACTGGCGGTTGGGTGGTTGGGGGGGGGGATTTTTAGCGATGGGTTACTGCACAAAACCCCCCGGGAAATTCTCTGGGTTaaaaa                             >  2:1098614/1‑90 (MQ=255)
                                                     tataCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTc                          <  2:883829/90‑1 (MQ=255)
                                                       taCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAc                        <  2:678858/90‑1 (MQ=255)
                                                        aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc                       <  1:984387/90‑1 (MQ=255)
                                                          tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGGGATCGATTACTGGACAAAACCCCGCGGTAATTCCTCTGGATCAAAAACTCAcccc                     >  2:562984/1‑89 (MQ=255)
                                                          tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCa                     >  1:932975/1‑90 (MQ=255)
                                                          tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCa                     >  1:566996/1‑90 (MQ=255)
                                                          tGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCa                     <  1:151128/90‑1 (MQ=255)
                                                          tGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGGAAATCCTCTGGATCAAAAAATCACCCa                     >  2:658001/1‑90 (MQ=255)
                                                              ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaa                                                            <  1:249789/48‑1 (MQ=255)
                                                              ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaa                                                            >  2:249789/1‑48 (MQ=255)
                                                              ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAAc                            >  2:537028/1‑79 (MQ=255)
                                                              ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAAc                            <  1:537028/79‑1 (MQ=255)
                                                                  ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc             <  2:375822/90‑1 (MQ=255)
                                                                  ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc             <  2:1478295/90‑1 (MQ=255)
                                                                     tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGc          >  2:57720/1‑90 (MQ=255)
                                                                     tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGc          <  2:1273509/90‑1 (MQ=255)
                                                                      ggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCt         >  1:1373335/1‑90 (MQ=255)
                                                                          ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa     <  2:577206/90‑1 (MQ=255)
                                                                          ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa     <  2:1020887/90‑1 (MQ=255)
                                                                          ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa     <  2:463623/90‑1 (MQ=255)
                                                                          ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCa     <  1:1098841/90‑1 (MQ=255)
                                                                            ggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAg   >  1:1127217/1‑90 (MQ=255)
                                                                             gaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGt  >  1:49110/1‑90 (MQ=255)
                                                                              |                                                                                         
GTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGT  >  NZ_CP009273/4611989‑4612155

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