Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 9.3% intergenic (+72/‑33) yqhD → / → dkgA NADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5900CA 9.3% 129.6 / 5.9 55intergenic (+72/‑33)yqhD/dkgANADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (3/2);  total (27/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 = 5.28e-01

TGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAA  >  NC_000913/3156448‑3156730
                                                                                                                                              |                                                                                                                                            
tGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                                                              <  1:114237/143‑1 (MQ=255)
  aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACATAAGTTCACAACACCTCCCCGGAGCCTg                                                                                                                                            >  2:583479/1‑143 (MQ=255)
    ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGcc                                                                                                                                          >  2:207957/1‑142 (MQ=255)
    ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATGTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCt                                                                                                                                          >  2:8204/1‑143 (MQ=255)
         tGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                     <  1:861940/143‑1 (MQ=255)
          gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                                                                    <  2:819289/143‑1 (MQ=255)
          gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTGCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                                                                    <  1:643178/143‑1 (MQ=255)
              gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                                                                <  2:246080/143‑1 (MQ=255)
              gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTCCGGGCCTTTCGTCCAGACTTCAAGTCACAACCCCTCACCGGAGACTGCACCGGTGAggc                                                                                                                                >  2:661223/1‑141 (MQ=255)
               aaaaTCATGACATTACGTTGGATCTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTc                                                                                                                               <  1:814731/143‑1 (MQ=255)
                  aTCATGACATTCCGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCTAGACTTAAGTTCACAACACCTCACCGTAGCCTGCTCCGGTGAGTTCata                                                                                                                            <  1:471472/143‑1 (MQ=255)
                   gcATGACACTACGCTGGATGTCAGCCGCCGTATATCCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAACCTTCCGTTTTCGGGCATTTCGCTCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCatat                                                                                                                           <  1:65245/142‑1 (MQ=255)
                       gACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                                                                       <  2:223745/143‑1 (MQ=255)
                          aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                    >  2:417998/1‑143 (MQ=255)
                                  ggATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAt                                                                                                            <  2:770935/143‑1 (MQ=255)
                                    aTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                          >  1:355410/1‑143 (MQ=255)
                                     gggcAGACGCCGGATATACGCCGCCGTCCGCGAAGCTGTGTACGGGTCTTTCGTGCGTTCGCCGGCCTTTCGTCCAGACTTAAGTTCACACCACCTCTCGTGAGCCTGCTCCGGTGAGTTCATATACGGGAGGAACGTATGGc                                                                                                         <  1:208784/140‑1 (MQ=255)
                                          gccgccGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCGTTTCGTCCAGACTTAAGGTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAACGAAGGAACGTCTGGCTAATc                                                                                                    >  1:366009/1‑143 (MQ=255)
                                              ccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                >  1:55466/1‑143 (MQ=255)
                                                 tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACACCACCTCACCGGAGCCTGCTCCGGTGAGTTCATATCAAGGAGGAACGTATAGCTCATCCAACCGt                                                                                             >  2:505658/1‑143 (MQ=255)
                                                       ccaaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                       <  2:79554/141‑1 (MQ=255)
                                                             cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGCGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGGACGGATGGCTAATCCCACCGCGATTAAGCtgac                                                                                 >  2:136110/1‑140 (MQ=255)
                                                                 cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCCCACTTAAGTTCACACCACCTCACCGGAACCTAAACAGAAAAGATAATATAAAGGAGGAACGTATGGCTAATCAAACCGTTATTAAGCTACAGGAt                                                                             >  1:386181/1‑143 (MQ=255)
                                                                 cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTCAGTTCACAACACCTCACCGGAACCTGCTCCGGAGAGTTCATATAGAGGAGGAAACAAAAGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                             >  1:382055/1‑143 (MQ=255)
                                                                 cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGTTGAGTTCATATAAAGGTGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                             <  1:64192/143‑1 (MQ=255)
                                                                 cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                             <  1:899753/143‑1 (MQ=255)
                                                                      cgCTTTTTCCGCCTCCAACTTTGGTCTTCGGGCATTTCGTCCGGACTTAAGTTCACAACACCTCACCGGAGCTTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                                                        <  2:551415/142‑1 (MQ=255)
                                                                                 ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                                                             >  1:157194/1‑143 (MQ=255)
                                                                                     aaaCTTTCGTTTTCGAGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGCGCCTGCTCCGGTGAGTTCATCTAAAGGAGGAACGTATGGCTAAACCAACCGTCATTAAGCAACAGAAAGGCAAAGACATGCCCCAGCt                                                         >  2:773010/1‑143 (MQ=255)
                                                                                       gcTTTCGTTCTCGGGCATTTCGGCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                       <  2:683836/142‑1 (MQ=255)
                                                                                         tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACCACACCTCCCCGGAGCCTGCTCCGGTGAGTTCAAATACAGGAGCACCGTATGGCTACTCCCACCCTTATTAACCTGCAGGATGGCAATGTCACGCCCCACCTgggc                                                     >  2:404111/1‑142 (MQ=255)
                                                                                         tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGTCAATGTCCTGCCCCAGCTGGGa                                                     >  2:487613/1‑143 (MQ=255)
                                                                                         tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                     >  1:858767/1‑143 (MQ=255)
                                                                                         tttCGTTTGCGCCCATTTCGCCCACGCTTAACTGCACACCACCTCTCCGCAACCTGCTCCCGTGAGTTCCTATAAAGGCGGACCGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                     <  2:495988/143‑1 (MQ=255)
                                                                                         gtgcGTGTTCGGGCATATCGTCCAGACTTAAGCTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                     <  2:490130/140‑1 (MQ=255)
                                                                                          tgccttttCGGGCATGGCGTCCAGCCTCACGTTCACAACACCTCCGCGGGGCCTGCTCCAGTGATCTCCTATAACGGAGTAACGTAGGGCTAATCCAACCGGTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                    <  2:106402/139‑1 (MQ=255)
                                                                                             gtgttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                 <  1:754451/140‑1 (MQ=255)
                                                                                                 tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGt                                             <  1:613013/143‑1 (MQ=255)
                                                                                                   gggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACCGGAAGGCAATGTCATGCCCCACCTGGGCCTGCGCGTCt                                           >  2:573292/1‑143 (MQ=255)
                                                                                                   gggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATAGCTAAACCAACCGATATTAAGCTACAGGATGGCAATGTCATGCCCCATCTGGGAAGGGGCGTcg                                           >  2:688963/1‑142 (MQ=255)
                                                                                                    ggCATTTCGTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGGGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                          <  2:395116/143‑1 (MQ=255)
                                                                                                     gCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTCCGCTACAGGATGGCAACGTCATGCCCCAGCTGGGACTGGGCGTCTgg                                         <  1:448877/143‑1 (MQ=255)
                                                                                                      cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgc                                        <  1:507849/143‑1 (MQ=255)
                                                                                                        tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                      <  1:258183/143‑1 (MQ=255)
                                                                                                        tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                      <  1:8204/143‑1 (MQ=255)
                                                                                                            gTCCAGCCATAGGTGCACAGCACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                  <  2:355410/143‑1 (MQ=255)
                                                                                                            gTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTGCTGCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                  <  1:769160/143‑1 (MQ=255)
                                                                                                                    ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGCGGAACGTATGGCTAATCCAACCGTACTTAAGCTACAGGATGGCAATGTCATGCCGCAGCTGGGACTGGGCGTCTGGCAAGCAAGTCATGAg                          >  2:872999/1‑143 (MQ=255)
                                                                                                                         ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGt                     >  1:260009/1‑143 (MQ=255)
                                                                                                                                   cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGcc           >  1:294573/1‑143 (MQ=255)
                                                                                                                                    tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa          >  1:518807/1‑143 (MQ=255)
                                                                                                                                    tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAGGTCATGCCCCAGCTGGGACTGGGAGCATGGCAAGCACGGAAGAAGTAAGCAATCACCGCCa          >  2:336554/1‑143 (MQ=255)
                                                                                                                                    tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAGCCAACCGTCATTAAGCTACAGGATGGCAATGTCATGCCCCANCTGGGACTGGGCGTCTGGCAAGCACGTGATAAGGACGTAATCACCGCCa          >  1:568615/1‑143 (MQ=255)
                                                                                                                                      aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGAAAGCGAGTAATGAGGAAGTAATCACCGCCAtg        >  1:463312/1‑142 (MQ=255)
                                                                                                                                      aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGGTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAGCGCCAtt        >  1:198042/1‑143 (MQ=255)
                                                                                                                                            gCCTGCTCCGGTGAGTTCATATACAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa  <  1:579231/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
TGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAA  >  NC_000913/3156448‑3156730

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