Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG41.0% 89.2 / 28.7 79intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (9/37);  new base G (32/0);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 5.70e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGT  >  NZ_CP009273/687219‑687386
                                                                                    |                                                                                   
tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCt                                                                                >  2:1426592/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:424145/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:2087942/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                               >  1:2468950/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                               >  1:1736486/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                               >  1:890671/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                               >  2:523658/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                               >  2:1522157/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTg                                                                               >  1:197907/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCGCATGGCGTCGTATGGGGTCTCATGTTAGGTAATCAAAACGCCGTCGGTAATCACCAACGGGGGTGTCTTCTg                                                                               >  2:1700823/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGcc                                                                             >  1:1178496/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGcc                                                                             >  1:1617296/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                             >  1:945941/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCACGGCGGGGGCGTCTGCCa                                                                            >  1:2224609/1‑90 (MQ=255)
    tgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCa                                                                            >  1:1675558/1‑90 (MQ=255)
            cgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCAGGCaaaag                                                                    >  2:240993/1‑88 (MQ=255)
                 cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                                               <  1:721993/90‑1 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTGCCCGTCAAAAGCCTGGcta                                                            >  1:2461553/1‑88 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAATGCCTGGCaa                                                            >  1:1666884/1‑90 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGTTGCCCGTGAAAAACCCGGCaa                                                            >  1:1320840/1‑90 (MQ=255)
                      ggCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAat                                                          <  2:661638/90‑1 (MQ=255)
                          tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGGAAattatt                                                      >  2:1229578/1‑90 (MQ=255)
                          tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAATGGCTGGCAAattttt                                                      >  2:793455/1‑90 (MQ=255)
                          tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGGCCGGGAAaatatt                                                      >  2:679309/1‑90 (MQ=255)
                          tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGGCAAAAGCCTGGCAAattatt                                                      >  2:1591546/1‑90 (MQ=255)
                          tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGAAAattatt                                                      <  2:955594/90‑1 (MQ=255)
                          tcgtATGGGCTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGAAAattatt                                                      <  2:371911/90‑1 (MQ=255)
                           cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAAATAttc                                                     >  2:721993/1‑90 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTTCCCGTCAAAAGCCTGGGAAATTATTTTTtttc                                               >  1:351006/1‑90 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                               <  1:296736/90‑1 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                               <  2:1617296/90‑1 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATAACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                               <  1:894474/90‑1 (MQ=255)
                                  gTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTCTCg                                              >  2:1261113/1‑90 (MQ=255)
                                   ttCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCAAAGCAAAAGCCTGGCAAAGTATTTTTTTTCgg                                             >  1:1815061/1‑90 (MQ=255)
                                   ttCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAACCCGGCAAAATATTTTTTTTCgg                                             >  1:1134958/1‑90 (MQ=255)
                                    tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCCGTCAAAAGCCCGGCCAAATATTTTTTTTCGGc                                            >  1:1567261/1‑90 (MQ=255)
                                    tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                            <  2:890671/90‑1 (MQ=255)
                                    tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                            <  1:1261113/90‑1 (MQ=255)
                                    tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                            <  1:976087/90‑1 (MQ=255)
                                      tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCaa                                          >  2:141246/1‑90 (MQ=255)
                                      tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTTGGAAAATATTTTTTTTCCGCaa                                          >  2:2177163/1‑90 (MQ=255)
                                      tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCaa                                          <  1:1700823/90‑1 (MQ=255)
                                       cATGTTAGTTAATCAAAACGCCGTCGTTAATCCCCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAAt                                         <  1:852276/90‑1 (MQ=255)
                                          gttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGGAAATTATTTTTTTTCGGGAATgga                                      >  1:1217636/1‑88 (MQ=255)
                                          gttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATgga                                      >  2:372041/1‑88 (MQ=255)
                                            tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTTCCCGTCAAAAACCCGGCAAATTATTTTTTTTCGGCAATGTACg                                    >  1:1453899/1‑90 (MQ=255)
                                            tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACg                                    >  2:1333493/1‑90 (MQ=255)
                                              gttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGa                                  >  1:434311/1‑90 (MQ=255)
                                               ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGAt                                 <  2:945941/90‑1 (MQ=255)
                                               ttaATCAAAACGCCGTCGTTAATCACCACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  2:2224609/90‑1 (MQ=255)
                                                taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  2:459893/90‑1 (MQ=255)
                                                taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  2:2239890/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                           <  2:2385378/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                           <  1:1333493/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                           <  1:785742/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCCACGGCGGGGAAGTCTGCCCGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                      <  1:2507517/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                      <  1:528497/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                      <  2:1147140/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                      <  1:118003/90‑1 (MQ=255)
                                                           ccGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGGTCCTCCTAAcccc                     >  2:528497/1‑89 (MQ=255)
                                                           ccGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCa                     <  2:157504/90‑1 (MQ=255)
                                                             gtcgtTAATCACCAACGGCGGGGGCGGCTGCCAGGCAAAAGCCTGGCAAATTATTTTTTTTCCGGCATGGACGGATCCTCATAACCCcga                   >  1:157504/1‑90 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGGAAATTATTCTTTCTAGGCAATGTACGGATCCTCATAGCCCagaga                 <  1:171410/90‑2 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc                 <  1:1522157/90‑1 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc                 <  2:434311/90‑1 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc                 <  1:1201045/90‑1 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGAAATGTACGGATCCTCATAGCCCAGagc                 <  2:1178496/90‑1 (MQ=255)
                                                                gtTAATCACCAACGGCGGGGGCGTCTGCCAGTCCAAAGCCTGGCAAAATATTCTTTTTCGGGAAAGTACGGGACCCCATAACCCCGagca                >  1:1995376/1‑90 (MQ=255)
                                                                 ttAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagcaa               <  2:1675558/90‑1 (MQ=255)
                                                                 ttAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagcaa               <  2:653917/90‑1 (MQ=255)
                                                                     tCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCAt           <  1:424145/90‑1 (MQ=255)
                                                                      cACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa          <  1:1591546/90‑1 (MQ=255)
                                                                       aCCAACGGCGGGGGCGGCTGCCCGTCAAAAGGCTGGCAAAATATTTTTTTTCGGCGAAGTACGGGTCCCCATAACCCCAAGCAAGCATaa         >  1:2043959/1‑90 (MQ=255)
                                                                       aCCAACGGCGGGGACGTCTTCCCGTCAAAAGCCTGGCAAATTATTTTTTCTCGGCAATGTACGGATCCTCATAACCCAGAGCAAGCATaa         >  2:1501305/1‑90 (MQ=255)
                                                                       aCCAACGGCGGGGACGGCTGCCAGTCAAAATCCTGGCAAATTATTCTTTTTCGGCAAAGTACGGATCCCCATACCCCAGAGCAAGCATaa         >  2:2290789/1‑90 (MQ=255)
                                                                          aaCGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct      <  1:2059500/90‑1 (MQ=255)
                                                                           caggcggGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtctc     <  2:1821953/88‑1 (MQ=255)
                                                                           aCGGCGGGGACGTCTGCCAGTCAAATGCCTAGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtcta     <  1:1410384/90‑2 (MQ=255)
                                                                              gcggGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtctcctg  <  1:2034069/90‑4 (MQ=255)
                                                                                    |                                                                                   
TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGT  >  NZ_CP009273/687219‑687386

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