Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8090TG34.0% 92.6 / 14.4 51T350P (ACC→CCC) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base G (17/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.54e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGA  >  NC_000913/650678‑650944
                                                                                                                                   |                                                                                                                                       
tgctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGAc                                                                                                                                 <  1:1079006/140‑1 (MQ=255)
 gctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACg                                                                                                                                >  2:495225/1‑140 (MQ=255)
        aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGcc                                                                                                                         <  1:357016/140‑1 (MQ=255)
        aTTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGcc                                                                                                                         >  1:1195328/1‑140 (MQ=255)
          ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTTATTTTTCACCTGTTTCTCCTGCGGGGGGCTTTTGACGAGCTGCcgc                                                                                                                       >  2:217561/1‑140 (MQ=255)
             cAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTTCTTTTTCACCTGTTTCTCATGCGGGGGGGTTTTGACGAGCTGCCGCGTc                                                                                                                    >  2:870655/1‑140 (MQ=255)
                 tCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTgg                                                                                                                <  1:688604/140‑1 (MQ=255)
                    atatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCgg                                                                                                             <  2:570505/140‑1 (MQ=255)
                    atatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCgg                                                                                                             <  2:1124581/140‑1 (MQ=255)
                     tatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCggg                                                                                                            <  1:495225/140‑1 (MQ=255)
                      atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                           <  1:671747/140‑1 (MQ=255)
                      atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                           <  2:1147823/140‑1 (MQ=255)
                      atCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCAGTTTTCACCTGTTTCGCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGa                                                                                                           <  1:1256310/140‑1 (MQ=255)
                       tCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAg                                                                                                          <  2:237637/140‑1 (MQ=255)
                        ccTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGt                                                                                                         <  1:466411/140‑1 (MQ=255)
                          tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggggggc                                                                                                       >  1:460029/1‑135 (MQ=255)
                          tGCGTATCGAGTGGGGCGATGCGTATCATCACACCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGc                                                                                                       >  1:983264/1‑140 (MQ=255)
                           gCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCt                                                                                                      <  2:624705/140‑1 (MQ=255)
                                 cGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCa                                                                                                >  2:902189/1‑140 (MQ=255)
                                       gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt                                                                                          <  2:897340/140‑1 (MQ=255)
                                          cGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTCTGCa                                                                                       >  1:146196/1‑140 (MQ=255)
                                                      tcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTTGCGGGGGGGCTCAAACCGGGTCTCCAAATAATGCAGCg                                                                           >  2:338185/1‑140 (MQ=255)
                                                       cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTCTGCAAATAAAGCAGCgg                                                                          >  2:1230592/1‑139 (MQ=255)
                                                       cacaTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGGTCTGCAAATAATGCAGCgg                                                                          >  1:502970/1‑139 (MQ=255)
                                                         caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                                        <  1:1024644/140‑1 (MQ=255)
                                                         caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTGCAAAAAAAGCAGCGGGa                                                                        >  2:450055/1‑140 (MQ=255)
                                                           tCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGGTCAAGCAGGGTTTGCAAATAAAGCAGCGTGGCt                                                                      >  2:128652/1‑140 (MQ=255)
                                                            ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTTTTCAAAAAAAGCAGCGGGACTg                                                                     >  2:129497/1‑140 (MQ=255)
                                                              ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAGCAGGTTTTGCAAATAATGCAGCGTGACTTCa                                                                   >  2:618152/1‑140 (MQ=255)
                                                                 gACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAggg                                                                >  1:212144/1‑140 (MQ=255)
                                                                    tCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGAcc                                                             <  1:491346/140‑1 (MQ=255)
                                                                      aaGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGTGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAg                                                           >  1:242279/1‑140 (MQ=255)
                                                                           tGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCCGGGGGCTTTTGACAAGCTGCCGCGTCCTGGGGGGGGCGCCCAACCACGGTTTGCACAAAATGCGGCGAGACTGCAGGGACCCggggcg                                                      >  2:491346/1‑140 (MQ=255)
                                                                            gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACAGGGTCTTCAAAAAAAGCCGCGGGGCTGCAGGGGCCCGCGgggc                                                     >  2:477993/1‑140 (MQ=255)
                                                                             ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACAGGGTCTGCCAATAAAGGAGCGTGGCTGCAGGGGCCCGCGGcgcg                                                    >  2:1058817/1‑140 (MQ=255)
                                                                             ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTGCAAAAAATGCAGCGGGGCTGCAGGGGCCAGCGGcgcg                                                    >  2:1100125/1‑140 (MQ=255)
                                                                             ccTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCCGCCGCGTCCTGGCGGGGGGGCTCAAACAGGTTTTGCAAAAAATTCCACGTGGCTGCAGGGGCCAgcgggggg                                                    >  2:276071/1‑136 (MQ=255)
                                                                              cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCCCCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGGCCTGGCGGGGGGGCTTAAACAGGGTTTGCAAATAAAGCAGCGTGACTGCCGGGACCAGCGGcgcgg                                                   >  2:187393/1‑139 (MQ=255)
                                                                              cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGGTCAAGCAGGGTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa                                                   >  1:372980/1‑140 (MQ=255)
                                                                               tGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                                  <  2:730121/140‑1 (MQ=255)
                                                                                 cAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAg                                                <  1:902189/140‑1 (MQ=255)
                                                                                    cggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCg                                             <  2:992705/140‑1 (MQ=255)
                                                                                                     ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc                            <  1:276071/140‑1 (MQ=255)
                                                                                                             tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg                    <  1:187393/140‑1 (MQ=255)
                                                                                                                ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg                 <  1:477993/140‑1 (MQ=255)
                                                                                                                  tGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTa               <  2:983264/140‑1 (MQ=255)
                                                                                                                   gTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAc              >  2:506186/1‑140 (MQ=255)
                                                                                                                    tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc             >  2:844675/1‑140 (MQ=255)
                                                                                                                        tcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCggg         <  1:1058817/140‑1 (MQ=255)
                                                                                                                          aTGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGaa       >  2:717104/1‑140 (MQ=255)
                                                                                                                           tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg      >  2:902651/1‑140 (MQ=255)
                                                                                                                              gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCgg   >  2:356198/1‑140 (MQ=255)
                                                                                                                              gggggCCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATAGTTTTTCGCCAGCAGTTGACGTACCCGGGGAGCgg   >  1:1250452/1‑140 (MQ=255)
                                                                                                                               ggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGa  >  2:938180/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TGCTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGA  >  NC_000913/650678‑650944

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