Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,279,2120GT36.8% 64.9 / 12.6 38A126S (GCC→TCC) agaDPTS galactosamine transporter subunit IID
Reads supporting (aligned to +/- strand):  ref base G (3/21);  new base T (14/0);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 7.03e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGC  >  NZ_CP009273/3279075‑3279347
                                                                                                                                         |                                                                                                                                       
cGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTg                                                                                                                                         <  1:850083/138‑1 (MQ=255)
   tGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccag                                                                                                                                     <  2:846969/139‑1 (MQ=255)
     gaagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTTGTTTTGTTTTTGGCCGATTATGGCGGGAATTTGCTCATCATTTgccagc                                                                                                                                    <  1:95406/138‑1 (MQ=255)
      aagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcca                                                                                                                                  <  2:56173/139‑1 (MQ=255)
        gaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAgg                                                                                                                                <  1:168897/139‑1 (MQ=255)
                agaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACctg                                                                                                                        <  2:814852/139‑1 (MQ=255)
                        gcgACACCATTAAAGGCCTCAAAGTGGCGCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGcc                                                                                                                <  1:171955/139‑1 (MQ=255)
                                    aaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                    <  1:551959/139‑1 (MQ=255)
                                    aaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                    <  1:244290/139‑1 (MQ=255)
                                      ggCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCg                                                                                                  <  2:296028/139‑1 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTTTCCTCATTTTCCCCCCAGGGGAACCTTCTTGGGGCGGTTCTTTTTTTTGCCg                                                                                               >  2:724746/1‑139 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTTTCATCATTTTCCACCCCGGGAAACCTGCTGGGGGCGCTTCTATTTTTTGCCg                                                                                               >  1:557831/1‑139 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTCTCCTCCTTTTCCCACCAGGGGAACCTGCTGGGGGCGGTTCTttttttttccc                                                                                               >  1:445801/1‑134 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCCTCATTTTTCAGCCAGGGGAAACTGCTGGGGGCGGTTTTAtttttttccg                                                                                               >  2:274239/1‑134 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAGTTTTTCATCATTTTCCAGCCAGGGGAAACTGCTGGGGGCGGGTCTATTTTTTGCCg                                                                                               >  2:751260/1‑139 (MQ=255)
                                         cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTCTCATCATTTTCCTGCCAGGGGAACCTGCTGGGGGCGGTTCTTTTTTTCTCCg                                                                                               >  2:307708/1‑139 (MQ=255)
                                             aaGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGATGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTa                                                                                           <  2:483609/139‑1 (MQ=255)
                                                        tttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTCTCTTTCTTTTTCAGGCCGGGGAACCCGGTGGGGGCGGGTCTTTTTTTTTCCCTTTTTCTTCTTAttt                                                                                >  2:510797/1‑137 (MQ=255)
                                                                aaTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGGTTTTGTCAACACTTTTCCGCCCGGGGAACCTGATGGGGGCGGTTCTTTTTTTTTCCGTTTTTCTGCTTTTCTTTTTTCtt                                                                        >  1:771002/1‑138 (MQ=255)
                                                                aaTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGTTCATCATTTTCCAACCAAGGAAACCCCCCGGGGCCGATTCTATTTTTTTCCGCTTTTCCGCTTTTCTTTTTTCTg                                                                        >  1:46077/1‑139 (MQ=255)
                                                                aaTCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg                                                                        <  1:751260/139‑1 (MQ=255)
                                                                  tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGCTCCTCATTTTCCAGCCCGGGGAAACCGCCGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTgcg                                                                      >  2:155853/1‑139 (MQ=255)
                                                                  tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCCTCCTTTTTCAGCCCGGGGAACCCTCTGGGGCCGGTTTTTTTTTTTTCCGTTTTCCTGCTTATTTTTTTTCTtcg                                                                      >  2:42934/1‑139 (MQ=255)
                                                                         gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGc                                                               <  1:95451/139‑1 (MQ=255)
                                                                              gcgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCCTCATTTTCCAGCCAGGGGAAACTGCTGGGGGCGGTTCTATTTTTTTCCGTTTTCCTGCTTATTTTTTTTCTTTGCGTCGGCTgggc                                                          >  2:544422/1‑137 (MQ=255)
                                                                              gcgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGa                                                           <  1:520026/138‑1 (MQ=255)
                                                                                  tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                      <  1:274239/139‑1 (MQ=255)
                                                                                      atTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCg                                                  <  2:502712/139‑1 (MQ=255)
                                                                                                 tACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGt                                        <  2:760/138‑1 (MQ=255)
                                                                                                 tACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGt                                        <  1:413049/138‑1 (MQ=255)
                                                                                                        ttgCCGATTATGGCGGGGAATTTCTCATCCTTTTCCCGCCCGGGGAAACTGCTGGGGGCGATTTTATTTTTTTCCGTTTACCCTCTTTTCTTTTTTCTTCGCGTCGGCTGGACCCACGTCGGGTTTTTAGTCGGCGTg                                 >  1:14658/1‑138 (MQ=255)
                                                                                                         tgCCGATTATGGCGGGAAATTGCTCATCATTTGCCAGCCAGGGGAAACTGCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTTCTTATTTTTTTTCTCCGCGTCGGCCGGAACCACGTCGGGTATTTAGTCGGCGTg                                 >  1:644552/1‑137 (MQ=255)
                                                                                                           ccGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAg                              <  1:724746/138‑1 (MQ=255)
                                                                                                            cGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTGTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGc                            <  2:450045/139‑1 (MQ=255)
                                                                                                                  tGGCGGGAATTTGCTCATCATTTTCCAGCCAAGGGAACCTGCTGGGGGCGGTTCTATTTTTTTCCGTTTTCCTGCTTATCTTTTTTCTTCGCGTCGGCTGGACCCCCCGCGGGTTTTCTGTCGGCGTGCAGGcgggtga                      >  2:389917/1‑134 (MQ=255)
                                                                                                                  tGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcga                      <  2:445801/139‑1 (MQ=255)
                                                                                                                    gCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATa                    <  2:163353/139‑1 (MQ=255)
                                                                                                                              gCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCga            >  1:288016/1‑137 (MQ=255)
                                                                                                                               cTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCga            >  1:313523/1‑136 (MQ=255)
                                                                                                                                  atcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgagaga        <  2:555803/137‑1 (MQ=255)
                                                                                                                                  atcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAAc      <  1:510797/139‑1 (MQ=255)
                                                                                                                                      tttGCCAGCCAGGGAAACCTGCTGGGGTCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCTCGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGGGCGAGAGAACTCGc  <  2:684100/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGC  >  NZ_CP009273/3279075‑3279347

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