Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,476 G→A 6.9% G39G (GGC→GGT)  gntK ← D‑gluconate kinase, thermostable

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,4760GA 6.9% 146.9 / 3.2 58G39G (GGC→GGT) gntKD‑gluconate kinase, thermostable
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base A (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

ACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAAT  >  NC_000913/3577336‑3577617
                                                                                                                                            |                                                                                                                                             
actaTCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGcc                                                                                                                                             <  1:643587/140‑1 (MQ=255)
aCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGcc                                                                                                                                             <  1:129817/143‑1 (MQ=255)
aCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGGGGATGGAGGGAATCGcg                                                                                                                                             >  1:844631/1‑142 (MQ=255)
aCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTCTTTCGATATTGCGCCGTGGATGGAGGAAATCGcc                                                                                                                                             >  2:173381/1‑143 (MQ=255)
aCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGCTTGCGATCGTCGTCATTCAGTCGTTCGCCAGACGCCATCCTCTCGATACTGCGCCGTGGATGCCGGCAACCGcc                                                                                                                                             >  2:448097/1‑143 (MQ=255)
 cGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                                                                            >  1:609467/1‑143 (MQ=255)
 cGATCAGCGACACTTTATCAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                                                                            <  1:786236/143‑1 (MQ=255)
       gCGACACTTTAGTGGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGa                                                                                                                                      <  2:856627/143‑1 (MQ=255)
              tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCgg                                                                                                                               <  1:842268/143‑1 (MQ=255)
              tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCgg                                                                                                                               <  2:791958/143‑1 (MQ=255)
                   taGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGCAATCGCCATCAAGAAACGCGGCATGa                                                                                                                          >  1:572370/1‑143 (MQ=255)
                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGt                                                                                                                       <  2:43724/143‑1 (MQ=255)
                         gcTGCCTGGCAAACGCGGCGTCGTTGAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGa                                                                                                                    <  2:787123/143‑1 (MQ=255)
                           tGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGGGATTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGACGAAATGGCCATCAAGACACGCGGCATGAACTTAATg                                                                                                                  >  1:458460/1‑143 (MQ=255)
                              aTAGGAAACGCGGCGTCGTTCAGCGCCTGCAGCCGCGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATgcgc                                                                                                               <  2:509613/143‑1 (MQ=255)
                              aTAGCAAACGCGGCGTCGATAAGCGCCGACAGCCCCCGCTTCCGAGCGGCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATgcgc                                                                                                               <  1:390858/143‑1 (MQ=255)
                                aGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCa                                                                                                             <  2:104397/143‑1 (MQ=255)
                                  cAAACGCGGCGTCGTTCAGCGCCTGCAACCCCGGTTTGCGATCGTCGTCATTCCGTGGTTCGCCCGACGCCATTTTTTCGCTATTGCGCCGTGGATGGAGGAAATCGCCCGCCCGCAACGCGGCATGAAGTTGATGCGCCAcg                                                                                                           >  1:48480/1‑142 (MQ=255)
                                      cgcgGCGTCTTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCAc                                                                                                       <  1:565944/143‑1 (MQ=255)
                                         ggcgTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTgg                                                                                                    <  2:1033080/143‑1 (MQ=255)
                                               gtTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAGTCGCCATCAAGGAACGCGGCATGAAGTTGATGCGCCCCTTCGCTGGCGACCg                                                                                              >  1:1042226/1‑143 (MQ=255)
                                               gtTCAGCGCCTGCAACCACGGTTGGCGATCGTCGGCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCg                                                                                              <  1:825849/143‑1 (MQ=255)
                                               gtTCAGCGCCTGCAACCACGGCTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCg                                                                                              >  1:966084/1‑143 (MQ=255)
                                                  cAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAg                                                                                           <  2:609467/143‑1 (MQ=255)
                                                        cTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGc                                                                                     <  1:496443/143‑1 (MQ=255)
                                                         tGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGcc                                                                                    <  1:889344/143‑1 (MQ=255)
                                                               cACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGCCGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGcc                                                                              <  2:643721/143‑1 (MQ=255)
                                                                    cttGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAACCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATa                                                                         <  1:496723/142‑1 (MQ=255)
                                                                      tGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCCTGAAGTTGATGAGCCACTTCACTGGCGCCCGCAGATTTGCCGCTGCCCGATACg                                                                       >  2:665608/1‑143 (MQ=255)
                                                                         gATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATCCGCAACTTCACTGGCGCCGGCAGCTTTGCCGCGGCCCGCTACGgcc                                                                    >  2:731282/1‑143 (MQ=255)
                                                                                        gTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAAtg                                                     >  1:899900/1‑143 (MQ=255)
                                                                                         tGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGAAGTAAAtgt                                                    >  1:758695/1‑143 (MQ=255)
                                                                                            ttCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCAATCAAGAAACGCGGCATGGAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATgtggt                                                 >  1:489674/1‑143 (MQ=255)
                                                                                                    aCGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGAGACGCCCATCAAGACGAGAATGTGGGGagcgcgg                                         >  2:989195/1‑137 (MQ=255)
                                                                                                      gCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGtt                                       <  2:233992/143‑1 (MQ=255)
                                                                                                      gCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGCAACTGTGGTGCTCacgatt                                       >  2:869573/1‑138 (MQ=255)
                                                                                                         aTTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCAACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGCTACGCCCATCAAGACGTAACTGTGGTGATCATGGTTAGt                                    >  1:154029/1‑143 (MQ=255)
                                                                                                          ttttttCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTc                                   >  1:184840/1‑143 (MQ=255)
                                                                                                                cGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTc                             <  1:728486/143‑1 (MQ=255)
                                                                                                                   tatGGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTGGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCaaa                          <  2:293088/143‑1 (MQ=255)
                                                                                                                          ccGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATCTGGTGATCATGGTTAGTCGTACTCAAAGTTGTGc                   >  1:693050/1‑143 (MQ=255)
                                                                                                                           cGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGcc                  >  1:248244/1‑143 (MQ=255)
                                                                                                                               gatggaGGAAATCACCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa                                                               <  2:134167/94‑1 (MQ=255)
                                                                                                                               gatggaGGAAATCACCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa                                                               >  1:134167/1‑94 (MQ=255)
                                                                                                                               gatggaGGAAATCACCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa                                                               >  1:659983/1‑94 (MQ=255)
                                                                                                                               gatggaGGAAATCACCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa                                                               <  2:659983/94‑1 (MQ=255)
                                                                                                                                  ggaggaAAGCGCCAGCAGGAGGCGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACaa           <  2:512352/143‑1 (MQ=255)
                                                                                                                                   gaggaAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAAt          >  1:503447/1‑143 (MQ=255)
                                                                                                                                   gaggaAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGATGATCATGGTTAGTCGTGCTCACAGTGGCGcaccaccact          >  2:303304/1‑137 (MQ=255)
                                                                                                                                    aggaAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATa         >  1:429430/1‑143 (MQ=255)
                                                                                                                                       aaaTCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACaa      <  1:364860/143‑1 (MQ=255)
                                                                                                                                       aaaTCGCCAGCAAGACACGCGGCATGAAGTTGATGCGCCACGTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACaa      <  1:330944/143‑1 (MQ=255)
                                                                                                                                        aaTCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACacg     >  1:854379/1‑141 (MQ=255)
                                                                                                                                         aTCGCCATCAATGAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGa    <  2:537039/143‑1 (MQ=255)
                                                                                                                                         aTCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGa    >  1:703194/1‑143 (MQ=255)
                                                                                                                                          tCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGaa   <  1:688985/143‑1 (MQ=255)
                                                                                                                                          gcGCCATCAAGAACCGCGGCATGAAGCTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGaa   <  1:294098/142‑1 (MQ=255)
                                                                                                                                           cGCCACCAAGAAGCGCGGCATGAGGGTGATGCGCCACTTCATTGGCGACCGCAGATCTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAAt  <  2:630491/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAAT  >  NC_000913/3577336‑3577617

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