Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG32.6% 63.0 / 28.0 43T255P (ACC→CCC) citGtriphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (2/27);  new base G (14/0);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 1.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGC  >  NC_000913/646608‑646879
                                                                                                                                           |                                                                                                                                     
aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGt                                                                                                                                       <  2:158352/140‑1 (MQ=255)
             gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCt                                                                                                                         >  1:1345279/1‑140 (MQ=255)
              cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGcccccct                                                                                                                        >  1:1384168/1‑138 (MQ=255)
                 aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATTCCCCCctttt                                                                                                                     >  1:605338/1‑140 (MQ=255)
                      atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCcttttttcat                                                                                                                >  2:1202792/1‑135 (MQ=255)
                        tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCCTTTTGCaataa                                                                                                              >  1:1270329/1‑140 (MQ=255)
                             cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGttt                                                                                                         <  1:418396/140‑1 (MQ=255)
                               gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCTATAAAATTTgt                                                                                                       >  2:382005/1‑139 (MQ=255)
                                tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCCATAAAGTTTtcg                                                                                                      >  2:439253/1‑140 (MQ=255)
                                tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTAATAAAGTTTtcg                                                                                                      >  1:1035763/1‑140 (MQ=255)
                                    aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                  <  1:559340/140‑1 (MQ=255)
                                       aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                               <  1:921151/140‑1 (MQ=255)
                                                    tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTTATAATTTTTTCGCCCCCCGCTGTAGCCCgcgc                                                                                  >  2:905544/1‑140 (MQ=255)
                                                      aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTTAATAATGTTTTCGCCCCGCGCCGTGGCCAgcgccg                                                                                >  2:540686/1‑138 (MQ=255)
                                                      aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATATTTTCGCCCCGCGCCGTAGCCAgcgccg                                                                                >  1:1057523/1‑138 (MQ=255)
                                                      aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGGCCCCCTTTTGTAATAAAATTTTCGCCCCCCGCTGTAGCCGgcgccg                                                                                >  2:1081515/1‑138 (MQ=255)
                                                         aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                             <  1:439253/140‑1 (MQ=255)
                                                         aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                             <  1:546666/140‑1 (MQ=255)
                                                         aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                             <  2:1345279/140‑1 (MQ=255)
                                                            tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGAGGCCCCCCTTTTTTAATAAATTTTTTGCCCCCCCCTCCAGCCCGCGCCCGccccc                                                                          >  2:1110483/1‑140 (MQ=255)
                                                             cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTAGCCAGCGCAGGcccccc                                                                         >  1:1144304/1‑140 (MQ=255)
                                                               gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAATATTTTTGCCTCCCGCTGTGGCCCGCGCAGGcccccccc                                                                       >  2:1150670/1‑138 (MQ=255)
                                                               gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTCTAGCCCACGCAGGcccccccc                                                                       >  2:612041/1‑138 (MQ=255)
                                                                cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATTCCCCCCTTTTTCTATAATGTTTTTGCCCCCCCGTGTAGCCCGGGCAGGccccccccc                                                                      >  2:1393415/1‑137 (MQ=255)
                                                                 aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGc                                                                     <  1:905544/140‑1 (MQ=255)
                                                                  cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAAAAATATTTTTGCCTCCCCCCGTAGCCCGCGCCGGccccccccccc                                                                    >  2:1234255/1‑135 (MQ=255)
                                                                        ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                              <  1:100176/140‑1 (MQ=255)
                                                                        ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                              <  2:1270329/140‑1 (MQ=255)
                                                                            gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGAt                                                          <  1:1393415/140‑1 (MQ=255)
                                                                             aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                         <  1:382697/140‑1 (MQ=255)
                                                                             aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                         <  1:266401/140‑1 (MQ=255)
                                                                              cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                        <  1:1284980/140‑1 (MQ=255)
                                                                              cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                        <  1:819695/140‑1 (MQ=255)
                                                                               tGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCa                                                       <  2:448058/140‑1 (MQ=255)
                                                                                      tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCCCCCCGTAGCCCCCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGCTgg                                                >  2:428758/1‑140 (MQ=255)
                                                                                         ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATATTTTCGCCTCCCGCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGATGCAACGGTGGGAt                                             >  2:1230975/1‑140 (MQ=255)
                                                                                           cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                           <  1:413041/140‑1 (MQ=255)
                                                                                                  cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                    <  1:1150670/140‑1 (MQ=255)
                                                                                                  cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                    <  1:1268145/140‑1 (MQ=255)
                                                                                                    cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCTAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                  <  1:333292/140‑1 (MQ=255)
                                                                                                     ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                 <  2:469999/140‑1 (MQ=255)
                                                                                                      cTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAAAGTTTTCGCCCCCCCCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGCTGCCACGGTGGGATTGCCGGTGGTCGc                                >  2:719354/1‑140 (MQ=255)
                                                                                                       tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc                               <  1:981097/140‑1 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                              <  2:1057523/140‑1 (MQ=255)
                                                                                                            aaCTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGGCCCCCTTTTTTAATAAATTTTTCGCCCCCCGCCGTTGCCAGCGCCGGGCCCCCCCCCCCCCCCCCGCTGCAACGTTGGTATTGCCGGTTTTCGCCCTCac                          >  2:298123/1‑139 (MQ=255)
                                                                                                             aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                         <  1:719354/140‑1 (MQ=255)
                                                                                                               tGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGGCCCCCTTTTTCCATAAATTTTTCGTCTCGCCCTGGTGACAAGGCACGTCCCCCTCGCCACCGCGCCCTGCCACATTGGGAACGCCCTTTATCCCCCTCaccag                       >  2:562746/1‑136 (MQ=255)
                                                                                                                   ggAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAAAGTTTGTATCGCCCTTGGTCGCCCTCAGgaggagg                   >  2:1191522/1‑139 (MQ=255)
                                                                                                                    gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                  <  1:540686/140‑1 (MQ=255)
                                                                                                                      gaTAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAAAGTTTTTGCCCCGCCCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACCGTGGGATCCCCCCTGATCGCCCTCCGGAGGAGGAAg                >  2:825159/1‑140 (MQ=255)
                                                                                                                          aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt            <  2:1384168/140‑1 (MQ=255)
                                                                                                                               gaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc       <  2:493567/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc  <  2:1144304/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
AAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGC  >  NC_000913/646608‑646879

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