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

TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAG  >  NZ_CP009273/687219‑687392
                                                                                    |                                                                                         
tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCt                                                                                      >  2:276450/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                     >  1:14449/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                     >  1:828943/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                     >  1:389943/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAGCGGCGGGGGCGTCTGcc                                                                                   >  2:829019/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGcc                                                                                   >  1:742365/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGcc                                                                                   >  1:150601/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                                   <  1:459405/90‑1 (MQ=255)
           tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCGACCAGACaaaa                                                                           >  1:600130/1‑89 (MQ=255)
             gtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCGAAAGc                                                                         >  2:404838/1‑90 (MQ=255)
              tcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGcc                                                                        <  1:656459/90‑1 (MQ=255)
                                 ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                                     <  1:264494/90‑1 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGTCAAAAGCCCGGCAAAATATTTTTTCTCGGCCAAGTAc                                           >  1:91224/1‑90 (MQ=255)
                                                 aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGAGGCTGCCCATCAAAAACCCCGGAAAATATTTTTTCTCCGCCATGGACGGAtgc                                     >  1:636300/1‑88 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                                 <  2:14449/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                                 <  2:459685/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                                 <  2:671909/90‑1 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                                 <  2:824907/90‑1 (MQ=255)
                                                      aaaCGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAATGCCTGGCAAATTATTCTTTTTCGGCAATGTACGGATCCTCATa                                >  2:506701/1‑90 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                            <  2:507276/90‑1 (MQ=255)
                                                                                    aCGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAg  <  1:556331/90‑1 (MQ=255)
                                                                                    |                                                                                         
TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAG  >  NZ_CP009273/687219‑687392

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