Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1436,067,3540AC40.0% 79.8 / 15.8 56V162G (GTG→GGG) cheB_5Protein‑glutamatemethylesterase/protein‑glutamin e glutaminase
Reads supporting (aligned to +/- strand):  ref base A (3/30);  new base C (22/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.77e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTTAC  >  ACB143/6067222‑6067484
                                                                                                                                    |                                                                                                                                  
aaCGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCtt                                                                                                                                <  2:125852/137‑1 (MQ=255)
     cTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCg                                                                                                                           <  8:244387/137‑1 (MQ=255)
         aTCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCGTc                                                                                                                       >  8:3809/1‑137 (MQ=255)
           gacacaTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCgg                                                                                                                     <  5:178213/136‑1 (MQ=255)
           cacccaTGGCGGCAGGCCACTGGCCAGGCGCGGGGGGCCGTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCgg                                                                                                                     <  1:172677/137‑1 (MQ=255)
           cacacaTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCgg                                                                                                                     <  3:327934/137‑1 (MQ=255)
           cacacaTGGCTGCATGCCACTGGCCAGGCGCAGCGGGCCCGGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCgg                                                                                                                     <  6:68315/137‑1 (MQ=255)
                 tGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGGGCCGCCCCCTTGACAGCATCGGGCGTGc                                                                                                               >  5:160598/1‑137 (MQ=255)
                  ggCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGcc                                                                                                              <  4:78913/137‑1 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGAACCCTTAGGGCTTGCCCGCAGTcc                                                                                                      >  6:178213/1‑136 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGCCCGCTTCGGGCTTGCCCGCAGtaa                                                                                                      >  2:280306/1‑135 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTAGGCGCCGGCCCCTTTACCGCTGGGGGCTTGCGCGCAGTcc                                                                                                      >  1:7087/1‑136 (MQ=255)
                                    aGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCCCCTGGACCGCGTCGGGCGTGCCCGCCGTCAGcgcggtgca                                                                                            >  7:314924/1‑132 (MQ=255)
                                    aGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGGTCaa                                                                                            >  8:222170/1‑137 (MQ=255)
                                        gcAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAg                                                                                        <  4:213398/137‑1 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCGTCGGGCGTGCCCGCCGTCCGCGCGGTGAACGCGt                                                                                       >  2:69837/1‑137 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCGGCGGGGTTGCCCGCAGTGAGCGCGTTGAACGCGt                                                                                       >  6:377210/1‑137 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCTTTCCCGCGGAGGGCTGGGCCGCCGGCCGCGCGTGCAAcacgg                                                                                       >  2:278423/1‑133 (MQ=255)
                                            ctggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCGTCGGGCGTGCCCGCCGTCCGCGCGGTCAACGcgtcgt                                                                                    >  8:78710/1‑137 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCTTGCCCGCGTCGGGCTTGCCCGCCGTGCGCGCGTGCAACGCGTCGTTg                                                                                  >  2:345744/1‑137 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCTTGACCGGGTCGGGCTTGCCCGCCGGCAGCGCGTTGAACGAGTCGTTg                                                                                  >  2:350546/1‑137 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTg                                                                                  <  3:82310/137‑1 (MQ=255)
                                                   tggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCCCCCTGACCCCGTCGGGGGTGGGCGCCGGCAGCGCGGTGAACGAGGCGTTgggcaa                                                                             >  1:29928/1‑133 (MQ=255)
                                                   tggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCTTGACCGCGTCGGGCGTGCCCGCCCGGAACGCGTGAAACGAGTCGGTTGCaaa                                                                             >  4:212774/1‑137 (MQ=255)
                                                   tggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCCCTTCGGGCGTGCCCGCAGTCAGCGCGTTCCACGGGTCGTTGGCaaa                                                                             >  4:220506/1‑137 (MQ=255)
                                                       cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGcc                                                                         <  4:85622/137‑1 (MQ=255)
                                                       cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTGGGCGCCGCCCCCGTGACCGCGTCGGGGGTGGCCGCCGGCAGCGCGGGGAACCAGGCGTTGGGACAgggc                                                                         >  3:92266/1‑135 (MQ=255)
                                                       cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCTTGACCGCGTCGGGGGTGCCCGCAGTGCGCGCGTTCCACGCGTCGGTGGCCAAGGcc                                                                         >  2:236026/1‑137 (MQ=255)
                                                       cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTAGGCGCCGCCGCCATGACCGCAGCGGGCATGCCCGCAGTCAGCGCGGTCAACCAGGCGGTGGCAAAGGcc                                                                         >  7:240176/1‑137 (MQ=255)
                                                            ccGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGGCCGCGGCGGGCGTGCCCGCAGTCAGCGCGTTGACCCAGGCGTGTGGAAAGGCCGGGtt                                                                    >  5:356247/1‑137 (MQ=255)
                                                            ccGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGACCCCTTCGGGCTTTCCCGCCGGCAGCGCGTTGAACGAGTCGTTGGCAAAGGCCcgggt                                                                    >  6:123191/1‑135 (MQ=255)
                                                            ccGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGtt                                                                    <  6:276100/137‑1 (MQ=255)
                                                            ccGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGtt                                                                    <  3:326803/137‑1 (MQ=255)
                                                            ccGCGCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGtt                                                                    <  1:261183/137‑1 (MQ=255)
                                                                  ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCACCTTGACCGCGTCGGGGTTGCCCGCACGCAGGGCGTTCAACGAGTAGTTGGCAAAGGCCGGGTTGGCTcg                                                              >  8:58461/1‑137 (MQ=255)
                                                                          ggcggcAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGc                                                      <  7:237162/137‑1 (MQ=255)
                                                                           gcggcAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGCCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGcc                                                     <  2:30510/137‑1 (MQ=255)
                                                                                tggTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATgg                                                <  6:138122/136‑1 (MQ=255)
                                                                                   tGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCa                                             <  2:9132/137‑1 (MQ=255)
                                                                                   tGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCa                                             <  3:6076/137‑1 (MQ=255)
                                                                                         cgccgcGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGc                                       <  6:382908/137‑1 (MQ=255)
                                                                                         cgccgcGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGc                                       <  8:314924/137‑1 (MQ=255)
                                                                                          gccgcGACCAAAGGCTTGGCCGGTGCCGGGGTGGGCGCCGCCCCCGTGGCCCCGTGGGGCGTGCCCGCCGGCAGCGCGGGGAACGCGGGGGTGGGCAGGGCCCGGGTGGCGCGCGCCGGGCCGCTGGGCCGCAGGcg                                      >  2:8993/1‑136 (MQ=255)
                                                                                               gACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCATTACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTGAACGCGGCGGTGGGAAAGGCCGGGTTGGCGCGCGCCCGGCCGGTGGGCCGAAGGCTGTCgc                                 >  2:77340/1‑136 (MQ=255)
                                                                                                aCCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCTTGACCGCTTCGGGCTTTCCCCCCGGCAGCGCGTTCAAAGTGTGGGTGGGAAAGGCCGGGTTGGCTCGCGCCGGGCCGGTGGGCATACGGCTggggtg                                >  4:29025/1‑132 (MQ=255)
                                                                                                  cAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCt                              <  1:69837/137‑1 (MQ=255)
                                                                                                       gcgtGGCCGGTGCCGGTGTAGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCg                         <  8:77670/134‑1 (MQ=255)
                                                                                                            gccggtgccggtGTCGGCGCCGCCACCTTGACCTCTTCTGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCATGCTGTCGAGCTTACCGATCtt                    <  1:157177/137‑1 (MQ=255)
                                                                                                                 tgccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa               <  2:182265/137‑1 (MQ=255)
                                                                                                                 tgccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa               <  7:170670/137‑1 (MQ=255)
                                                                                                                 ggccggtgGCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa               <  3:7713/136‑1 (MQ=255)
                                                                                                                  gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg              <  7:105624/137‑1 (MQ=255)
                                                                                                                  gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg              <  6:141965/137‑1 (MQ=255)
                                                                                                                            gcgcCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGTCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGctttcttt    <  8:240176/137‑1 (MQ=255)
                                                                                                                            gcgcCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGctttcttt    <  2:29928/137‑1 (MQ=255)
                                                                                                                              gccgccACCTTTACCGCGTCGGGCTTTCCCGCAGTCAGCGCGTTCAAACAGTCGTTGGCAAAGGCCGGGTTTGCTCGCGCCGGGCCGCTGGCCATCAGGCTGTGGAGCTTAACGATCTTTGTGAGGGCTTTCTTTaa  >  6:127608/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                  
AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTTAC  >  ACB143/6067222‑6067484

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