Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,743,478 T→C V249A (GTC→GCC)  yiaK → 2,3‑diketo‑L‑gulonate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,4780TC96.4% 182.6 / ‑3.0 56V249A (GTC→GCC) yiaK2,3‑diketo‑L‑gulonate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/27);  total (28/28)
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.99e-01

ATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTG  >  NC_000913/3743343‑3743612
                                                                                                                                       |                                                                                                                                      
aTTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCAcc                                                                                                                                    <  1:868396/140‑1 (MQ=255)
   tGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAg                                                                                                                                 <  1:307488/140‑1 (MQ=255)
          aGACCCTGGCGTTATCGAGAAGAATCGCTGCATTTTGCCGATGGGCTCCTGTAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTACCGACGGCGCATCCGTTGCCGAAGCCCCCCAGGacaaca                                                                                                                          <  2:250157/140‑1 (MQ=255)
               cTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGcgac                                                                                                                     >  2:358432/1‑140 (MQ=255)
                           agaagaATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCAtt                                                                                                         >  1:414648/1‑140 (MQ=255)
                                  tCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACaaa                                                                                                  >  2:112863/1‑140 (MQ=255)
                                     ccgcATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAAttt                                                                                               >  2:358414/1‑140 (MQ=255)
                                      cgcATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAAtttt                                                                                              >  2:625919/1‑140 (MQ=255)
                                         aTTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTAt                                                                                           <  1:388179/140‑1 (MQ=255)
                                          ttttGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCAACCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTAtt                                                                                          >  2:206505/1‑140 (MQ=255)
                                             tGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGcc                                                                                       >  2:2352/1‑140 (MQ=255)
                                              gCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCa                                                                                      <  1:206505/140‑1 (MQ=255)
                                                  aTGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGa                                                                                  >  1:999876/1‑140 (MQ=255)
                                                   tGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGaa                                                                                 <  1:1192098/140‑1 (MQ=255)
                                                       cTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTgg                                                                             >  1:612962/1‑140 (MQ=255)
                                                          cTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTACCATTGAAGTGGACa                                                                          <  1:358432/140‑1 (MQ=255)
                                                                      tGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACg                                                              <  2:1084609/140‑1 (MQ=255)
                                                                      tGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACg                                                              <  1:112863/140‑1 (MQ=255)
                                                                           tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTccc                                                         <  1:1042878/140‑1 (MQ=255)
                                                                           tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTccc                                                         >  1:779683/1‑140 (MQ=255)
                                                                            gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGt                                                            <  2:1215049/136‑1 (MQ=255)
                                                                            gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGt                                                            <  2:287154/136‑1 (MQ=255)
                                                                            gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGt                                                            >  1:1215049/1‑136 (MQ=255)
                                                                            gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGt                                                            >  1:287154/1‑136 (MQ=255)
                                                                            gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCa                                                        <  1:973699/140‑1 (MQ=255)
                                                                             tCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCAc                                                       >  1:969535/1‑140 (MQ=255)
                                                                                   gtgCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGa                                                 <  1:625919/140‑1 (MQ=255)
                                                                                        ggATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCa                                            <  1:1084003/140‑1 (MQ=255)
                                                                                        ggATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCa                                            <  2:393569/140‑1 (MQ=255)
                                                                                           tatGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGc                                         >  1:1118888/1‑140 (MQ=255)
                                                                                               aTCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCa                                     <  2:1118888/140‑1 (MQ=255)
                                                                                                tctcTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACAGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGGGATGCCAAGCTGCaa                                    <  2:168736/137‑1 (MQ=255)
                                                                                                tCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCaa                                    >  2:486897/1‑140 (MQ=255)
                                                                                                tCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCaa                                    <  1:2352/140‑1 (MQ=255)
                                                                                                tCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCaa                                    >  1:121536/1‑140 (MQ=255)
                                                                                                tCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCaa                                    <  2:999876/140‑1 (MQ=255)
                                                                                                   cTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGc                                 >  1:164309/1‑140 (MQ=255)
                                                                                                    tACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGGGGACAAGCTTATCAACGGTCCCACCCGCGATGCCAAGCTGCAACGca                                >  2:362901/1‑140 (MQ=255)
                                                                                                     aCTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGcat                               >  1:855210/1‑140 (MQ=255)
                                                                                                     aCTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGcat                               <  2:121536/140‑1 (MQ=255)
                                                                                                      ctctCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGcatc                              <  2:855210/140‑1 (MQ=255)
                                                                                                         tcCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATg                           >  1:85276/1‑140 (MQ=255)
                                                                                                         tcCTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATg                           <  2:923514/140‑1 (MQ=255)
                                                                                                           cTTTCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGa                         >  2:886419/1‑140 (MQ=255)
                                                                                                              tCCGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTa                      >  1:229558/1‑140 (MQ=255)
                                                                                                               ccGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTAc                     <  1:928350/140‑1 (MQ=255)
                                                                                                                cGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATgg                          <  1:764835/134‑1 (MQ=255)
                                                                                                                cGACGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATgg                          >  2:764835/1‑134 (MQ=255)
                                                                                                                   cGGCGCATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTa                 <  2:969535/140‑1 (MQ=255)
                                                                                                                     gcgcATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGcat                               <  2:1093836/124‑1 (MQ=255)
                                                                                                                     gcgcATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGcat                               >  1:1093836/1‑124 (MQ=255)
                                                                                                                     gcgcATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATgg                          <  2:444799/129‑1 (MQ=255)
                                                                                                                     gcgcATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATgg                          >  1:444799/1‑129 (MQ=255)
                                                                                                                        cATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGt            >  1:1196188/1‑140 (MQ=255)
                                                                                                                        cATCCGTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGt            >  1:1065450/1‑140 (MQ=255)
                                                                                                                             gTTGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGa       <  2:821966/140‑1 (MQ=255)
                                                                                                                              ttGCCGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAg      >  1:1212478/1‑140 (MQ=255)
                                                                                                                                  cGAAGCCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTg  <  1:886419/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTG  >  NC_000913/3743343‑3743612

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