Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,581,618 G→A 7.5% V661V (GTG→GTA tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,581,6180GA 7.5% 174.7 / 4.5 53V661V (GTG→GTAtktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base A (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

ATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGG  >  NC_000913/2581482‑2581757
                                                                                                                                        |                                                                                                                                           
aTTACTGGTACAAGTATGTTGGTCTGAACGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGGTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAACGCGCATAAGGTGCTg                                                                                                                                          <  2:1084473/140‑1 (MQ=255)
aTTACTGCTACAACTATGTCGGTCTGCAAGGGGCAATTGTCGGCATGACGGGTCACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCATTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTg                                                                                                                                          <  2:816606/140‑1 (MQ=255)
 ttACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTgg                                                                                                                                         >  1:364579/1‑140 (MQ=255)
 ttACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGAGCTgg                                                                                                                                         >  1:364426/1‑140 (MQ=255)
       gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGa                                                                                                                                   >  1:689143/1‑140 (MQ=255)
       gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGa                                                                                                                                   >  1:1061016/1‑140 (MQ=255)
       gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGa                                                                                                                                   >  1:280048/1‑140 (MQ=255)
                 gTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTg                                                                                                                         >  1:14747/1‑140 (MQ=255)
                  ttGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGa                                                                                                                        >  2:86551/1‑140 (MQ=255)
                   tGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGAt                                                                                                                       <  1:607748/140‑1 (MQ=255)
                       cTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtg                                                                                                                   >  2:266384/1‑140 (MQ=255)
                        tGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtga                                                                                                                  <  1:1182033/140‑1 (MQ=255)
                        tGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtga                                                                                                                  <  1:322932/140‑1 (MQ=255)
                        tGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtga                                                                                                                  <  1:323083/140‑1 (MQ=255)
                            aGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgc                                                                                                              <  1:569889/140‑1 (MQ=255)
                               ggCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGATAATATTGTGGCAAAAGCGCATAAGGTACTGGGAGTGAAAGGTGCCTGATGGtgattgccgt                                                                                                           >  2:788426/1‑139 (MQ=255)
                               ggCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGATAATATTGTGGCAAAAGCGCATAAGGTACTGGGAGTGAAAGGTGCCTGATGGtgattgccgt                                                                                                           >  2:788796/1‑139 (MQ=255)
                               ggCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccgg                                                                                                           >  2:448091/1‑140 (MQ=255)
                                gCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccgga                                                                                                          >  1:288544/1‑140 (MQ=255)
                                    ttGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccggatgct                                                                                                      >  1:505119/1‑140 (MQ=255)
                                    ttGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccggatgct                                                                                                      >  1:505228/1‑140 (MQ=255)
                                                             aTCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCTGATGCGACGCTGAcgcg                                                                             >  2:537904/1‑140 (MQ=255)
                                                             aTCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGAcgcg                                                                             >  1:1060781/1‑140 (MQ=255)
                                                             aTCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGAcgcg                                                                             >  1:371728/1‑140 (MQ=255)
                                                                  cTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTa                                                                        >  2:49786/1‑140 (MQ=255)
                                                                  cTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTa                                                                        >  2:352499/1‑140 (MQ=255)
                                                                   tCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                       <  1:385747/140‑1 (MQ=255)
                                                                   tCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                       <  1:138534/140‑1 (MQ=255)
                                                                   tCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                       <  2:14747/140‑1 (MQ=255)
                                                                   tCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                       <  2:284864/140‑1 (MQ=255)
                                                                   accGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                       <  2:58942/139‑1 (MQ=255)
                                                                      ggcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCg                                                                    >  2:808965/1‑140 (MQ=255)
                                                                      ggcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCg                                                                    >  2:809013/1‑140 (MQ=255)
                                                                      ggcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCg                                                                    <  2:387511/140‑1 (MQ=255)
                                                                            tAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTAc                                                              >  1:959084/1‑140 (MQ=255)
                                                                                 tGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGt                                                         <  1:1096104/140‑1 (MQ=255)
                                                                                      ccGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTAt                                                    <  1:631121/140‑1 (MQ=255)
                                                                                      ccGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTAt                                                    <  1:195081/140‑1 (MQ=255)
                                                                                      ccGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTAt                                                    <  1:1002881/140‑1 (MQ=255)
                                                                                                gCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGTAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACa                                          <  2:1123604/140‑1 (MQ=255)
                                                                                                       cgccgAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCg                                   <  2:160504/140‑1 (MQ=255)
                                                                                                           gagaATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCAt                               <  1:1007977/140‑1 (MQ=255)
                                                                                                             gaATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCAttt                             <  1:1008521/140‑1 (MQ=255)
                                                                                                                     tgGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAAc                     <  2:112409/140‑1 (MQ=255)
                                                                                                                     tgGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAAc                     <  2:1109672/140‑1 (MQ=255)
                                                                                                                     tgGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAAc                     <  1:1125019/140‑1 (MQ=255)
                                                                                                                        cAAAAGCGCATAAGGTGCTGGGAGTGAATGGTGCCTGATGGGGATTGCCGGATGCTGATTGCCGGATGCGACGCTGTCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACggg                  <  2:339002/140‑1 (MQ=255)
                                                                                                                          aaaGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGTTGGTGATAGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTg                >  2:202876/1‑140 (MQ=255)
                                                                                                                               gcATAAGGTACTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGTATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCca           <  1:788796/140‑1 (MQ=255)
                                                                                                                               gcATAAGGTACTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGTATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCca           <  1:788426/140‑1 (MQ=255)
                                                                                                                                    aGGTGCTGGTAGTGAAAGGTGCCTGATGGTGATTGCCGGATTCTGATTGCCGTATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACg      <  1:1169498/140‑1 (MQ=255)
                                                                                                                                    aGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACg      <  1:1169415/140‑1 (MQ=255)
                                                                                                                                     ggTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGt     >  2:964999/1‑140 (MQ=255)
                                                                                                                                        gCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTggg  >  1:841539/1‑140 (MQ=255)
                                                                                                                                        gCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTAATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTggg  >  1:441767/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                           
ATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGG  >  NC_000913/2581482‑2581757

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