Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273762,1470CA56.0% ‑0.1 / 22.6 25L216L (CTC→CTAmngAPTS 2‑O‑a‑mannosyl‑D‑glycerate transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base A (0/14);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTC  >  NZ_CP009273/762059‑762217
                                                                                        |                                                                      
aaaaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCa                                                                       <  1:1168370/90‑1 (MQ=255)
  aaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAAt                                                                     <  2:1372537/90‑1 (MQ=255)
                     gCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTAttt                                                  >  2:1610629/1‑90 (MQ=255)
                         aaaTTTTTACGTAAAAAAACAAGGCGCTAACGGCATTGAAGGGCTTTTAACGGGGGATCAACTAAATAGTGCAACCGCCTGTATTTTTGc                                              <  2:542401/90‑1 (MQ=255)
                               ttACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGa                                        >  2:797257/1‑90 (MQ=255)
                                   gtaaaaaaaCAAGGCGCTAACGGCATTGAGGGGCTTTAAAGGGGGGACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTc                                    <  1:911797/86‑1 (MQ=255)
                                   gTTGAAAACCAGGGCCCAACGGGGATGGAGGGGGGTTAAAGGGCGAATAAACCAAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTc                                    <  2:1847482/90‑1 (MQ=255)
                                     taaaaaaCAGGGCGCAAACGGCTTTGAGGGGCTTTAAAGGGGGGATAAACAAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  1:716653/88‑1 (MQ=255)
                                     taaaaaaCAAGGCGCAAACGGCATTGAGGGGCTTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  2:618934/88‑1 (MQ=255)
                                     tGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  >  2:1512855/1‑90 (MQ=255)
                                     tGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAATAACTCAATAATGCAAACCCCTGTATTTTTTTGGCTGGAggggc                                  >  1:52832/1‑86 (MQ=255)
                                     aaaaaaaaAGGGCGCAACCGGCATTGAGGGGCTTTTAAGGGGGGATAACCTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  2:312860/88‑1 (MQ=255)
                                     aaaaaaaaAGGGCCCTACCGGAATGGAGGGGTTTTAAACGGGGAACAAACAAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  1:1610629/88‑1 (MQ=255)
                                     aaaaaaaCAGGGCGCAACCGGCATGGAGGGGCTTTAAGCGGGGGACCACCCAAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  1:1262775/88‑1 (MQ=255)
                                     aaaaaaaCAGGCCGCAAACGGAATGGAGGGGCTTTAAAGGGGGAACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc                                  <  1:797257/88‑1 (MQ=255)
                                      gAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTc                                                                        <  1:711208/51‑1 (MQ=255)
                                      gAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTc                                                                        >  2:711208/1‑51 (MQ=255)
                                         aaaCAGGGCGCAACGGGCATGGAGGGGTTTTAAGGGGGGGACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATc                              <  1:210782/90‑1 (MQ=255)
                                             aaGGCGCAACGGGAATGGAGGGGCTTTTAAGGGGGGATCAACTAAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAgg                          <  2:1072230/90‑1 (MQ=255)
                                                    aaagggCTTTGAGGGGTTTTAAACGGGGGATCAACTAAATAGTCCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGa                   <  2:1108087/86‑1 (MQ=255)
                                                      agggCTTTGAGGGGTTTTAAAGGGGGGATAAACTAAATAGTCAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGc                 <  2:1658011/88‑1 (MQ=255)
                                                      agggCATTGAGGGGCTTTAAACGGGGGATCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGc                 <  2:1414423/88‑1 (MQ=255)
                                                            ttGAAGGGCGTTTAACGGCGGATCAACTCAAAAATGCAAACCCCCGGATTTTTTCGGCTGGAGGCGCCAACCAGGGGAGGGGGCGgtttt           >  1:1370382/1‑89 (MQ=255)
                                                             tGAGGGGTGTTAAACGGCGAAAAAACTAAATGGTCCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTaa          <  1:1764244/90‑1 (MQ=255)
                                                                   gCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGGCGCCCCCAAGGGGAGGGGGGGGTTTTATGGCAt    >  1:1215568/1‑90 (MQ=255)
                                                                     gTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGGCGCCATCAAAGAGAGGGAGCGGTTTTATGGCAttt  >  1:458426/1‑89 (MQ=255)
                                                                     gTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTc  >  2:1309212/1‑90 (MQ=255)
                                                                                        |                                                                      
AAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTC  >  NZ_CP009273/762059‑762217

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