Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 9.1% intergenic (+73/‑32) 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,5910TG 9.1% 122.2 / 6.1 55intergenic (+73/‑32)yqhD/dkgANADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base G (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

ACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAA  >  NC_000913/3156450‑3156730
                                                                                                                                             |                                                                                                                                           
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:8204/143‑1 (MQ=255)
                                                                                                      tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                      <  1:258183/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)
                                                                                                                                             |                                                                                                                                           
ACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAA  >  NC_000913/3156450‑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: ‑