Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273906,3700AG56.2% 8.2 / 19.9 16S46P (TCC→CCC) poxBubiquinone‑dependent pyruvate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/7);  new base G (8/1);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGAC  >  NZ_CP009273/906284‑906436
                                                                                      |                                                                  
gCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCAt                                                                 >  2:591525/1‑90 (MQ=255)
      tCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTc                                                           <  2:141364/90‑1 (MQ=255)
               gACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATGCACTCGGGGGTggc                                                  >  2:727417/1‑88 (MQ=255)
               gACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCCCCATgggggc                                                  >  1:205552/1‑86 (MQ=255)
                   gCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATg                                              <  2:593000/90‑1 (MQ=255)
                             cGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGGATCCCCCCCATGGTGGCCCCCCGGTTAAGAc                                    >  1:94109/1‑90 (MQ=255)
                               cTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTg                                  <  2:368948/90‑1 (MQ=255)
                               cTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTg                                  <  2:408762/90‑1 (MQ=255)
                                        gCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCTCGATGGGGGCCATGCGGTTTAGACTGTCACTAAGa                         >  2:146134/1‑90 (MQ=255)
                                         cTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGAc                        <  1:246159/90‑1 (MQ=255)
                                             aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCCCCATGGGGGCCCCGCGATTAAAAACGGCACAAAGAACGtt                    >  1:306910/1‑90 (MQ=255)
                                             aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACTCCATGGTGCCCCTCCGGTTTAGAATGTCACTAAAACCGtt                    >  2:437734/1‑90 (MQ=255)
                                             aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCACCCCCCCCATGGGGCCCCCCCGCTTAAAACTGTCCCCAAAAcccct                    >  1:545725/1‑87 (MQ=255)
                                                 ccAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCaga                <  1:518909/90‑1 (MQ=255)
                                                  cAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCagag               <  2:451220/90‑1 (MQ=255)
                                                               cgccCCTTCCTCGGGGGGGGGGGGAATCCCCTTGGTGGGGGCCATGCCATTAAGGCCGTTACTAAGACCGTTCAGAGAGTCGCCTGTGAc  <  1:437734/90‑1 (MQ=255)
                                                                                      |                                                                  
GCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGAC  >  NZ_CP009273/906284‑906436

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