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

AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTTACCTT  >  ACB143/6067222‑6067487
                                                                                                                                    |                                                                                                                                        
accGCCTTGACCCCACATGGCTGCAGCCCACGGGCCAGGCGCCGCGGGCCCGGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCt                                                                                                                                       <  5:216983/134‑1 (MQ=255)
   gCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGGCGCCAccaaga                                                                                                                                    >  3:136067/1‑132 (MQ=255)
     cTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCATGccc                                                                                                                                  >  2:25954/1‑136 (MQ=255)
         aTCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGACCGcgt                                                                                                                              >  7:258751/1‑134 (MQ=255)
          tCACACATGGCTGCATGCCACTGGCCAGGCGCAGCGGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTc                                                                                                                            <  5:295192/136‑1 (MQ=255)
            acacaTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCgg                                                                                                                          <  6:276931/136‑1 (MQ=255)
            acacaTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCgg                                                                                                                          <  7:196539/136‑1 (MQ=255)
            acacaTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCGCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCgg                                                                                                                          <  3:363495/136‑1 (MQ=255)
                aTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCtt                                                                                                                      <  6:326284/136‑1 (MQ=255)
                 tGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCCCCCCCGTTGCCGC‑GGCGGGCttt                                                                                                                     >  8:293381/1‑135 (MQ=255)
                   gCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCCCCCCCATGACCCC‑GTCGGGCGTGcc                                                                                                                   >  8:150619/1‑136 (MQ=255)
                      gCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCCCCGTGCCCGC‑GTCGGGGGTGccccc                                                                                                                >  4:216446/1‑134 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGGGCCGCCCCCATGACCGC‑TTCGGGGTTGCGCGCAgac                                                                                                            >  3:369506/1‑134 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGC‑GTCGGGCGTGCCCGCAGt                                                                                                             >  6:331971/1‑135 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCCCCCTGAACCC‑GTCCGGCGTGCCCGccgg                                                                                                             >  3:61166/1‑132 (MQ=255)
                          gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTGGGCGCCGCCCCCGGGACCGC‑GTCGGGCGTGCCCGCCGTc                                                                                                            >  7:33852/1‑136 (MQ=255)
                           ccACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCa                                                                                                           <  4:93527/136‑1 (MQ=255)
                           ccACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCCCCACCATTCCCCC‑GTCGGGGTTGccccccct                                                                                                             >  4:2459/1‑129 (MQ=255)
                           ccACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGGCGGGGCCCCCCCCATGGCCGC‑GTCCGGCGTGCCCGcccg                                                                                                             >  2:357786/1‑131 (MQ=255)
                              cTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTGAgcg                                                                                                        <  8:199335/136‑1 (MQ=255)
                               tGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGGGCCGCCCCCCTGACCGC‑GTCGGGCTTTGCCGCCGTGAgcgc                                                                                                       >  4:53476/1‑136 (MQ=255)
                                    aGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCACCTTGCCCCC‑TTCGGGCGTGCCCGCCGGCAGCGCGGTc                                                                                                   >  2:232159/1‑135 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGCCCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCCCCGTGTCCGC‑GGGGGGCTTGCCCGCGCGGCGCGCGGTCAACGAg                                                                                             >  3:39142/1‑136 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCCCCTTGACCGC‑GTCGGGGTTGCCCGCAGTCCGCGCGGTCCACGAg                                                                                             >  8:55390/1‑136 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGC‑TTCGGGGGTGCCCGCAGTCAGCGCGTTCAACgcc                                                                                             >  4:196253/1‑134 (MQ=255)
                                         cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCGTGACCGC‑TGCGGGGTTGCGCGCAGGGCTCGCGTTCAACgcg                                                                                             >  6:69669/1‑134 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCCGCCCCATGCCCGC‑GTCGGGCGTGCCCGGAGTGAGCGCGTTGAACGAGTCG‑tt                                                                                       >  2:230163/1‑136 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCCCCAAGGCCGC‑GTGGGGCGTGCGCGCCGTCAGCGGGGTCAACGcgtcg‑tg                                                                                       >  4:130406/1‑135 (MQ=255)
                                              ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTGGGCGCCGCCCCCGTGCCCGG‑TTCGGGCGTGACCGCCGTCAGCGCGGGCAACGAgggggt                                                                                        >  5:290190/1‑131 (MQ=255)
                                               gccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTg                                                                                      <  2:64223/136‑1 (MQ=255)
                                                       cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCGGGTGCCGGTGTCGGCGCCGCCCCCCTGACCCC‑GTCGGGGGTGCCCGCCGTCAGCGCGTGCAAAGAGTCG‑TGGGGAAA‑GGc                                                                             >  1:305642/1‑136 (MQ=255)
                                                        aCGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGcc                                                                            <  5:108086/136‑1 (MQ=255)
                                                             cGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGC‑TTAGGGGTTGCCCGCAGTCAGCGCGTTGAACGCGGCG‑GTGGGAAA‑GGCCggggt                                                                       >  4:308836/1‑134 (MQ=255)
                                                             cGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGtt                                                                       <  5:278192/136‑1 (MQ=255)
                                                                  ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTGGGGGCCGCCCCCATTACCGC‑GTCGGGGGTGGCCCCCGGCCGCGCGGTGGACGAGGGG‑GTGGGAAAGGGCCGGGTggggg                                                                   >  8:162387/1‑134 (MQ=255)
                                                                  ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGGGCCGCCCCCATGGCCGCGTTCGGG‑TTTCCCGCCGTCAGCGCGGTGAACGAGTCGATGGGGAAG‑GGCCGGGTggggc                                                                   >  6:68458/1‑134 (MQ=255)
                                                                                 ggTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATgg                                                   <  6:21884/136‑1 (MQ=255)
                                                                                    gCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCa                                                <  4:74191/136‑1 (MQ=255)
                                                                                    gCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑GTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCa                                                <  4:86672/136‑1 (MQ=255)
                                                                                        gcgcCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGC‑TTGGGGCGTGGCCGCCGGCCGCGCGTTGAAAGAGTCG‑GTGGGAAA‑AGGCGGGGTGGGTCGCGCCGGGGCGAGGGcccggcc                                            >  6:283335/1‑131 (MQ=255)
                                                                                                   aaaGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCt                                 <  6:108748/136‑1 (MQ=255)
                                                                                                             gcggtgccggtgTCGGCGCCGCGACCTTGACCGC‑TTCGGGGGTGCGCGCAGTCAGCGCGTTCTACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCAGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCtt                       <  3:338363/135‑1 (MQ=255)
                                                                                                             ccggtgccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCtt                       <  2:217033/136‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa                  <  1:131086/135‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa                  <  8:33852/135‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa                  <  3:198623/135‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCa                  <  6:126226/135‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg                 <  7:162387/136‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg                 <  3:266851/136‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg                 <  8:35976/136‑1 (MQ=255)
                                                                                                                   ccggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg                 <  8:41157/136‑1 (MQ=255)
                                                                                                                    cggtgTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAgg                <  1:34513/136‑1 (MQ=255)
                                                                                                                                  ccACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTTACCtt  <  4:136067/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                        
AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGC‑TTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCG‑TTGGCAAA‑GGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTTACCTT  >  ACB143/6067222‑6067487

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: ‑