Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8090TG32.6% 88.9 / 19.3 49T350P (ACC→CCC) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (9/24);  new base G (16/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATA  >  NC_000913/650677‑650946
                                                                                                                                     |                                                                                                                                         
ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTCACCTGTTTCTCTTGCGGGGGGCTTTTGa                                                                                                                                    >  2:381928/1‑140 (MQ=255)
ctgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACAACCCCAGCCTCAAGGGTGCCTGCAACGACGGACTGGTTTTTTTTTATTTTTCTGCTAATTCTTCAGGCGGGGTGCttttt                                                                                                                                     >  2:1170413/1‑139 (MQ=255)
   ctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTTACGa                                                                                                                                 >  2:694206/1‑140 (MQ=255)
         aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGcc                                                                                                                           <  2:795325/140‑1 (MQ=255)
             gCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGt                                                                                                                       <  1:621502/140‑1 (MQ=255)
              cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTc                                                                                                                      >  1:276483/1‑140 (MQ=255)
              cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTc                                                                                                                      >  1:1307513/1‑140 (MQ=255)
              cAGGTCGATATCCTGCGAATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCTGGCTCGTTTTTTTTCTCGTTTTACCCGGTTCTCGTGCGGGGGGTTTTTGATGGACCGCCGCGTc                                                                                                                      >  2:1119210/1‑140 (MQ=255)
                   cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGc                                                                                                                 <  1:1585893/140‑1 (MQ=255)
                   cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTCTTTTTCCACTGTTTCTCATACGGGGGGCTTTTGACGAGCTGCCGCGTTCTGGc                                                                                                                 >  2:1446379/1‑140 (MQ=255)
                           tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGc                                                                                                         >  1:1599639/1‑140 (MQ=255)
                           tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTTTCCTGCGGGGGGCTTTTTACGGGCTGCCGCGTCCTggggggggggc                                                                                                         >  2:336863/1‑135 (MQ=255)
                           tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggggggc                                                                                                         >  2:1585893/1‑135 (MQ=255)
                             cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTc                                                                                                       >  1:591137/1‑140 (MQ=255)
                             cGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGGTc                                                                                                       >  1:641112/1‑140 (MQ=255)
                               tATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCaa                                                                                                     >  1:1239106/1‑140 (MQ=255)
                                 tCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGc                                                                                                   <  1:381928/140‑1 (MQ=255)
                                  cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCa                                                                                                  >  1:1411843/1‑140 (MQ=255)
                                  cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAGCa                                                                                                  >  2:1257248/1‑140 (MQ=255)
                                   gAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAg                                                                                                 <  2:314790/140‑1 (MQ=255)
                                    aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTTCCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCTCCaggagg                                                                                                >  2:539523/1‑140 (MQ=255)
                                       ggggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAGCAgggtt                                                                                             >  1:1554393/1‑137 (MQ=255)
                                        gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                            <  1:539523/140‑1 (MQ=255)
                                        gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                            <  2:1599639/140‑1 (MQ=255)
                                        gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                            <  1:1257248/140‑1 (MQ=255)
                                         ggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGttttt                                                                                           >  2:658141/1‑137 (MQ=255)
                                         ggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTg                                                                                           >  2:834949/1‑140 (MQ=255)
                                                cGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAaataa                                                                                    <  2:1447312/140‑1 (MQ=255)
                                                  tATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATg                                                                                  <  2:519301/140‑1 (MQ=255)
                                                   atcatcACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGc                                                                                 <  1:795707/140‑1 (MQ=255)
                                                     catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAg                                                                               <  2:1693358/140‑1 (MQ=255)
                                                        cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTCTGCAAATAATGCAGCgg                                                                            >  1:885704/1‑139 (MQ=255)
                                                            tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTGCCAATAAAGCAGCggggct                                                                        >  2:1057245/1‑135 (MQ=255)
                                                            tCCCCAGACTCAAGGGGGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTCTGCAAAAAATGCAGCGTGGCt                                                                        >  2:1498647/1‑140 (MQ=255)
                                                                     tCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGAcc                                                               <  1:1057245/140‑1 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAg                                                             <  2:1411843/140‑1 (MQ=255)
                                                                       aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAACAGGGTCTGCAAAAAATGCAGCGTGACTGCAGGGGCCAg                                                             >  1:936601/1‑140 (MQ=255)
                                                                        aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                                            <  1:117847/140‑1 (MQ=255)
                                                                            tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcggcg                                                        <  1:1498647/140‑1 (MQ=255)
                                                                            tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTTTGCAAAAAATGCAGCGGGACTGCCGGGGCCAgcgggg                                                        >  2:670967/1‑138 (MQ=255)
                                                                             gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGgcgc                                                       >  1:257438/1‑140 (MQ=255)
                                                                              ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCAAATAAAGCAGCGTGACTGCAGGGACCAGCGGcgcg                                                      >  1:471328/1‑140 (MQ=255)
                                                                              ccTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGGAGGGTTTTCCAAAAAAGGAGCGGGAATGCCGGGGCCAGCGGggcg                                                      >  2:1107474/1‑140 (MQ=255)
                                                                               cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                                     <  1:1054306/140‑1 (MQ=255)
                                                                                tGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                                    <  1:1119210/140‑1 (MQ=255)
                                                                                  cAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCCAACAGGTTTTGCAAAAAATTCAGCGTGACTGCAGGGACCAGCGGCGCGATAg                                                  >  2:446875/1‑140 (MQ=255)
                                                                                      ggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGt                                              <  2:591137/140‑1 (MQ=255)
                                                                                                         aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagc                           >  2:1040628/1‑140 (MQ=255)
                                                                                                               cACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGa                     >  2:1095733/1‑140 (MQ=255)
                                                                                                                     tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc               >  1:661495/1‑140 (MQ=255)
                                                                                                                            aTGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGaa         <  2:782743/140‑1 (MQ=255)
                                                                                                                             tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg        <  1:670967/140‑1 (MQ=255)
                                                                                                                               cGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCg      >  2:831286/1‑140 (MQ=255)
                                                                                                                               cGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCg      <  2:936601/140‑1 (MQ=255)
                                                                                                                                gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg     >  1:1472405/1‑140 (MQ=255)
                                                                                                                                  gggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGat   <  1:834949/140‑1 (MQ=255)
                                                                                                                                   ggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGata  <  2:881174/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
CTGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATA  >  NC_000913/650677‑650946

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