Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,4890GA25.5% 87.5 / 29.9 51G825S (GGT→AGT) aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (6/7);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGC  >  NC_000913/125341‑125628
                                                                                                                                                    |                                                                                                                                           
gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                                                                                                                             <  3:131486/149‑1 (MQ=255)
gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                                                                                                                             >  6:110574/1‑149 (MQ=255)
                cTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCgg                                                                                                                             <  5:108839/149‑1 (MQ=255)
                           gcgcgTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa                                                                                                                    >  6:82860/1‑147 (MQ=255)
                             gcgTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGAc                                                                                                                   <  1:166538/146‑1 (MQ=255)
                             gcgTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGAc                                                                                                                   >  2:166538/1‑146 (MQ=255)
                               gTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTAc                                                                                                                                         <  1:217715/122‑1 (MQ=255)
                               gTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTAc                                                                                                                                         >  2:217715/1‑122 (MQ=255)
                                 tCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTg                                                                                                            <  6:23289/149‑1 (MQ=255)
                                      atatCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTg                                                                                                            <  4:74883/144‑1 (MQ=255)
                                      atatCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTg                                                                                                            >  3:74883/1‑144 (MQ=255)
                                      atatCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGt                                                                                                                            <  7:38140/128‑1 (MQ=255)
                                      atatCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGt                                                                                                                            >  8:38140/1‑128 (MQ=255)
                                               aGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTcac                                                                                              <  5:21247/149‑1 (MQ=255)
                                                      gaACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTc                                                                                        >  8:130487/1‑148 (MQ=255)
                                                        acgacgCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGaa                                                                                     <  5:110574/149‑1 (MQ=255)
                                                               tCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGa                                                                                <  3:64666/147‑1 (MQ=255)
                                                               tCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGa                                                                                >  4:64666/1‑147 (MQ=255)
                                                                  ggCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt                                                                                                                        <  2:48821/104‑1 (MQ=255)
                                                                   gCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACa                                                                                                                  <  2:53954/109‑1 (MQ=255)
                                                                   gCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACa                                                                                                                  >  1:53954/1‑109 (MQ=255)
                                                                       tGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTAt                                                                      <  5:82860/149‑1 (MQ=255)
                                                                       tGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTAt                                                                      <  5:173238/149‑1 (MQ=255)
                                                                         gCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGcttct                                                                         <  6:119212/144‑1 (MQ=255)
                                                                         gCATCTACAGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCTGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGcttct                                                                         >  5:119212/1‑144 (MQ=255)
                                                                               aCCGACAATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAg                                                                                    <  6:174234/127‑1 (MQ=255)
                                                                                ccGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGc                                                                                                   >  1:60142/1‑111 (MQ=255)
                                                                                ccGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGc                                                                                                   <  2:60142/111‑1 (MQ=255)
                                                                                ccGACTATATGAAACTGTTCGCTGAGCAGGGCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTgg                                                               <  7:130487/147‑1 (MQ=255)
                                                                                 cGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTcc                                                                                                                      >  4:14473/1‑91 (MQ=255)
                                                                                 cGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTcc                                                                                                                      <  3:14473/91‑1 (MQ=255)
                                                                                   aCTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg                                                                                                                              >  1:147300/1‑81 (MQ=255)
                                                                                   aCTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg                                                                                                                              <  2:147300/81‑1 (MQ=255)
                                                                                   aCTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTgcg                                                          >  2:219532/1‑149 (MQ=255)
                                                                                       tatGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTg                                                                >  1:49919/1‑139 (MQ=255)
                                                                                       tatGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTg                                                                <  2:49919/139‑1 (MQ=255)
                                                                                         tGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTc                                                                                                                       >  4:146695/1‑82 (MQ=255)
                                                                                         tGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTc                                                                                                                       <  3:146695/82‑1 (MQ=255)
                                                                                          gAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                                                                                                                             >  7:56458/1‑59 (MQ=255)
                                                                                          gAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                                                                                                                             <  8:56458/59‑1 (MQ=255)
                                                                                          gAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTc                                                                   >  2:174550/1‑133 (MQ=255)
                                                                                          gAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTc                                                                   <  1:174550/133‑1 (MQ=255)
                                                                                            aaCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGc                                                 >  3:41796/1‑149 (MQ=255)
                                                                                                        aGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTaaa                                     <  1:219532/149‑1 (MQ=255)
                                                                                                             gTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTAAACGTgg                                >  2:5394/1‑149 (MQ=255)
                                                                                                              tCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTAAACGTGGc                               <  1:5394/149‑1 (MQ=255)
                                                                                                                     ttACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGc                                         >  6:107295/1‑132 (MQ=255)
                                                                                                                     ttACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGc                                         <  5:107295/132‑1 (MQ=255)
                                                                                                                             cGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGTTgg                                                   >  6:12584/1‑114 (MQ=255)
                                                                                                                             cGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGTTgg                                                   <  5:12584/114‑1 (MQ=255)
                                                                                                                              ggCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTgg                                          >  5:194468/1‑122 (MQ=255)
                                                                                                                              ggCTGACGACTACCGCGTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTgg                                          <  6:194468/122‑1 (MQ=255)
                                                                                                                               gCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGcttctt                                                                        >  6:169011/1‑91 (MQ=255)
                                                                                                                               gCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGcttctt                                                                        <  5:169011/91‑1 (MQ=255)
                                                                                                                                            gcgTACTGAGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGc  >  8:165701/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAACCTGCGTCACCACTTCGAAGTTGATGCTTCTTATGTCGTGGTTGCGGCGCTGGGCGAACTGGCTAAACGTGGCGAAATCGATAAGAAAGTGGTTGCTGACGC  >  NC_000913/125341‑125628

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