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

TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTT  >  NZ_CP009273/687219‑687388
                                                                                    |                                                                                     
tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCt                                                                                  >  1:763359/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                 >  1:650411/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                 >  2:486747/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                 >  2:362883/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                 >  1:451860/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                 >  1:1133335/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                 >  2:173571/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTg                                                                                 >  2:1173914/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                               <  1:688242/90‑1 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                               <  1:310008/90‑1 (MQ=255)
      aaTTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGcccgg                                                                            >  1:482675/1‑87 (MQ=255)
         tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCCGCCAGTCaa                                                                         >  2:562434/1‑90 (MQ=255)
         tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTTCCCGCCaa                                                                         >  1:654362/1‑90 (MQ=255)
         tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCaa                                                                         <  2:1252682/90‑1 (MQ=255)
           tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTTCCCGTCaaaa                                                                       >  1:214243/1‑89 (MQ=255)
           tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGAGGGCTGCCAGTCAAAt                                                                       >  2:286252/1‑90 (MQ=255)
           tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGGTGCCCGTCcaaa                                                                       >  2:654965/1‑89 (MQ=255)
               cgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCt                                                                   <  2:1156451/90‑1 (MQ=255)
                 cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                                                 <  2:844916/90‑1 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGGCATTCAAAGCCCTGGcga                                                              >  1:364367/1‑88 (MQ=255)
                           cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCCGGCAAAAGCCCGGCAAagtattt                                                       >  1:170620/1‑89 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                                 <  1:603584/90‑1 (MQ=255)
                                    tcttATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTGCCAGTCAAATGCCTGGCAAACTATTCTTTCTCGGc                                              <  1:1115543/90‑1 (MQ=255)
                                    tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                              <  1:1297868/90‑1 (MQ=255)
                                      tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCaa                                            >  2:1252497/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTGCCAGTCAAAAAACTGGCAAATTATTTTTTTTCGGCAAAGTAc                                       >  1:115311/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCTGCAATGTAc                                       >  1:1508530/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTAc                                       >  2:1297868/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                       <  2:823162/90‑1 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                       <  2:1424023/90‑1 (MQ=255)
                                              gttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTCTTTCTCGGCGATGTACggg                                    >  1:1315137/1‑89 (MQ=255)
                                                taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                  <  1:963313/90‑1 (MQ=255)
                                                 aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAAATATTTTTTTTCGGCAAAGTACGGATcc                                 >  2:1457470/1‑90 (MQ=255)
                                                 aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCCGGCAAAATATTTTTTTTCGGCAAAGTACGGGTcc                                 >  2:1581332/1‑90 (MQ=255)
                                                 aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGACAAAAACCCGGCAAATTATTTTTTTTCGGCCATGTACTGAttc                                 >  1:1383955/1‑88 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                             <  1:1252497/90‑1 (MQ=255)
                                                       aaCGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCGAAGTACGGATCCCCCTAg                           >  1:1032546/1‑90 (MQ=255)
                                                         cGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGcc                         <  2:763359/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                        <  1:212668/90‑1 (MQ=255)
                                                           ccGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCa                       <  2:1383588/90‑1 (MQ=255)
                                                             gtcgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCaga                     <  2:1148522/90‑1 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc                   <  1:916231/90‑1 (MQ=255)
                                                                      cACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa            <  1:223726/90‑1 (MQ=255)
                                                                       aCCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCCGGCCAATTATTTTTTTTCCGCAATGTACGGGTCCTCATAGCCCAGAGCAAACATaa           >  1:1091586/1‑90 (MQ=255)
                                                                          aaCGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct        >  1:1424357/1‑90 (MQ=255)
                                                                                ggggGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTCTTTTTCGGCAAAGTACGGGTCCTCATAACCCCGAGCAAGCATAATCTCTGttt  >  1:1532825/1‑90 (MQ=255)
                                                                                    |                                                                                     
TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTT  >  NZ_CP009273/687219‑687388

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