Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG52.5% 45.3 / 74.9 59intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (8/20);  new base G (31/0);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-08
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGT  >  NZ_CP009273/687217‑687386
                                                                                      |                                                                                   
actGTGTGAATTGTCGTCGCTCATGGCGTCGGATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGt                                                                                  >  2:1282803/1‑90 (MQ=255)
  tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCt                                                                                >  2:1051679/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:145768/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:388995/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:29497/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                               >  2:159135/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                               >  2:863773/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                               >  2:833957/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                               >  1:1659946/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCtt                                                                               >  2:901027/1‑89 (MQ=255)
     gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCCCCAACGGCGGGGACGTCTGcc                                                                             <  1:355019/90‑1 (MQ=255)
     gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGcc                                                                             >  2:1110363/1‑90 (MQ=255)
     gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                             <  1:1145812/90‑1 (MQ=255)
        aaTTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGt                                                                          >  2:1404982/1‑90 (MQ=255)
          ttGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCa                                                                        >  2:1485635/1‑90 (MQ=255)
             tcgtcaCTCATGGCGTCGTATGGGTTCTCATGTTAGGTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGCCGGCCAGTCaaag                                                                     >  1:1258804/1‑89 (MQ=255)
               gtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTTAAAAGc                                                                   >  1:895305/1‑90 (MQ=255)
                   cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                                               <  2:934315/90‑1 (MQ=255)
                   cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                                               <  2:1051780/90‑1 (MQ=255)
                    tCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGGCAAATGCCTGGc                                                              >  2:1822663/1‑90 (MQ=255)
                    tCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAAtggctggc                                                              >  2:352818/1‑90 (MQ=255)
                            tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGcc                                                                             >  2:984875/1‑67 (MQ=255)
                            tcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                             <  1:984875/67‑1 (MQ=255)
                            cagtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAattatt                                                      <  2:1705358/88‑1 (MQ=255)
                             cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAACCTGGCAAattattt                                                     >  1:778887/1‑89 (MQ=255)
                                   ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAACCTGGGAAATTATTTTTtttc                                               >  1:1842473/1‑90 (MQ=255)
                                     ttCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCgg                                             >  2:675017/1‑90 (MQ=255)
                                      tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTTCCCGTCAAAAGCCTGGCCAAATATTTTTTTTCCGc                                            >  1:952092/1‑90 (MQ=255)
                                      cagcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                            <  2:1315449/87‑1 (MQ=255)
                                       ctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTGCCAGTCAAATGCCCGGCAAATTATTTTTTTTCGGCa                                           >  2:1210909/1‑90 (MQ=255)
                                        tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCAGTCAAAAACCGGGCCAAATATTTTTTTTCGGCaa                                          >  1:57675/1‑90 (MQ=255)
                                        tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTTTCGGCaa                                          >  2:289947/1‑90 (MQ=255)
                                         cATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCTGGCAAATTATTTTTTTTCGGCAAt                                         >  2:124290/1‑90 (MQ=255)
                                             ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTTTTTTTTTTCGGCAATGTAc                                     >  1:261716/1‑90 (MQ=255)
                                             ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                     <  1:1051679/90‑1 (MQ=255)
                                              tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCCGGCAAAAGCCCCGCAAAATATTCTTTCCCGGCAATGTACg                                    >  1:1515987/1‑90 (MQ=255)
                                               agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACgg                                   <  2:1513220/90‑1 (MQ=255)
                                                gttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTAAGGa                                  >  2:1285311/1‑90 (MQ=255)
                                                 ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGGTGCCAGTCCAATAACTTGCAAAATATTTCTTCTTCGCCAAGGACGGAt                                 >  1:1226819/1‑90 (MQ=255)
                                                 ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGAt                                 <  1:1548500/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  1:836073/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  1:933157/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  2:284302/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                <  1:291592/90‑1 (MQ=255)
                                                   aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGGCAAATGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGGTcc                               >  2:1006584/1‑90 (MQ=255)
                                                   aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGGTcc                               >  1:1731805/1‑90 (MQ=255)
                                                   aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAAATATTTTTTTTCGGCAATGTACgggctg                               >  1:1315449/1‑86 (MQ=255)
                                                   aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAAATATTTTTTTTCGGCAAAGTACGGGTcc                               >  2:1772259/1‑90 (MQ=255)
                                                       aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                           <  2:154514/90‑1 (MQ=255)
                                                            gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                      <  1:1285311/90‑1 (MQ=255)
                                                               gtcgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCaga                   <  1:1190448/90‑1 (MQ=255)
                                                               gtcgtTAATCACCAACGGCGGGGACGGCTGCCCGGCAAAAGCCTGGCAAAATATTCTTTTTCGGCAAAGGACGGGTCCCCATAGCCCcga                   >  1:308881/1‑90 (MQ=255)
                                                                    tAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTTTCGGCAATGTACGGATCCCCATAACCCAGagcaaa              >  2:1803494/1‑89 (MQ=255)
                                                                        caccacCGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa          <  1:1515333/90‑1 (MQ=255)
                                                                        cACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa          <  1:1772607/90‑1 (MQ=255)
                                                                          ccAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAt        >  2:695738/1‑90 (MQ=255)
                                                                            aaCGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGATCCTCATAGCCCCGAGCAAACATAAtct      >  1:592642/1‑90 (MQ=255)
                                                                            aaCGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct      <  2:653616/90‑1 (MQ=255)
                                                                                gcggGGACGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTTGGCAAAGTACGGATCCTCATAGCCCACAGCAAGCATAATCTCTGt  >  1:1664409/1‑90 (MQ=255)
                                                                                      |                                                                                   
ACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGT  >  NZ_CP009273/687217‑687386

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