Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,019,887 C→T Q169* (CAA→TAA)  dkgA → 2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,019,8870CT98.2% 182.4 / ‑3.9 57Q169* (CAA→TAA) dkgA2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/27);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GATTGATGAAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAACAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTA  >  minE/2019826‑2019949
                                                             |                                                              
gattgatGAAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGc                                                       >  1:355539/1‑71 (MQ=255)
gattgatGAAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGc                                                       >  1:2367492/1‑71 (MQ=255)
  ttgatGAAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCt                                                      >  1:2194192/1‑70 (MQ=255)
    gatgaAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTaca                                                   >  1:102033/1‑71 (MQ=255)
         aaCTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTacac                                                  <  1:1898530/67‑1 (MQ=255)
         aaCTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCt                                              <  1:649435/71‑1 (MQ=255)
          aCTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTg                                             >  1:2226675/1‑71 (MQ=255)
              gCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGaa                                          >  1:172307/1‑70 (MQ=255)
              gCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAAc                                         >  1:1970278/1‑71 (MQ=255)
                    cGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGaca                                   >  1:1142331/1‑71 (MQ=255)
                       cTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGacaca                                 >  1:122279/1‑70 (MQ=255)
                       cTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGacaca                                 >  1:1271647/1‑70 (MQ=255)
                       cTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGacaca                                 >  1:1370624/1‑70 (MQ=255)
                        tgtgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGcc                                                          <  1:913238/44‑1 (MQ=255)
                         gtgATAAACCAGATCGAACTTCATCCTCTGATGCAATAACGcc                                                          <  1:1250543/43‑1 (MQ=38)
                         gtgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGAcacaca                               >  1:676758/1‑70 (MQ=255)
                          tgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTgg                                            <  1:1092282/56‑1 (MQ=255)
                          tgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaa                              <  1:77061/70‑1 (MQ=255)
                          tgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaa                              >  1:736416/1‑70 (MQ=255)
                          tgATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaa                              >  1:2027288/1‑70 (MQ=255)
                           gATAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaaaa                            >  1:1141683/1‑71 (MQ=255)
                            aTAAACCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaaaa                            <  1:1230532/70‑1 (MQ=255)
                            aTAAACCAGATCGAACTTCATCCGCTGATGCAACAACGCCAGCTACACGCCTGGAACGCGACACACAAAAt                           <  1:1983309/71‑1 (MQ=255)
                              aaaCCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGac                                    >  1:1651449/1‑60 (MQ=255)
                              aaaCCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATc                          >  1:175439/1‑70 (MQ=255)
                               aaCCAGATCGAACTTCATCCGCTGATGCAATAACGc                                                           <  1:1487591/36‑1 (MQ=255)
                               aaCCAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATcc                         <  1:1179318/70‑1 (MQ=255)
                                 ccAGATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGa                      <  1:1407958/71‑1 (MQ=255)
                                    gATCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGAcc                    >  1:1704792/1‑70 (MQ=255)
                                     aTCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACg                                                 <  1:1531827/40‑1 (MQ=255)
                                     aTCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACg                                                 <  1:158187/40‑1 (MQ=255)
                                     aTCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAAt                           <  1:2415787/62‑1 (MQ=255)
                                     aTCGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGa                  <  1:971143/71‑1 (MQ=255)
                                       cGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGcc                                               >  1:1188602/1‑40 (MQ=255)
                                       cGAACTTCATCCGCTGATGCAATAACGCCAGCTACACGc                                                >  1:700109/1‑39 (MQ=255)
                                        gAACTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATc               <  1:164307/71‑1 (MQ=255)
                                         aaCTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAAt                           <  1:302719/58‑1 (MQ=255)
                                         aaCTTCATCCCCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAg                       <  1:2408739/62‑1 (MQ=255)
                                          aCTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATcc              <  1:1250385/70‑1 (MQ=255)
                                          aCTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATc               <  1:198655/69‑1 (MQ=255)
                                          aCTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCt             <  1:32890/71‑1 (MQ=255)
                                          aCTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCt             <  1:33555/71‑1 (MQ=255)
                                           cTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGAc                     >  1:1732122/1‑62 (MQ=255)
                                           cTTCATCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTg            >  1:1508194/1‑71 (MQ=255)
                                                tCCGCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGcc       <  1:2176485/71‑1 (MQ=255)
                                                   gCTGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCAt    <  1:2398254/71‑1 (MQ=255)
                                                    cTGATGCAATACCGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGAcc                    >  1:625579/1‑54 (MQ=39)
                                                     tGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaaaa                            >  1:1563849/1‑45 (MQ=255)
                                                     tGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAAt                           >  1:1826477/1‑46 (MQ=255)
                                                     tGATGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAAt                           >  1:2276198/1‑46 (MQ=255)
                                                       aTGCAATAACGCCAGCTACACGCCTGGAACGCGACACACaaaa                            <  1:2070978/43‑1 (MQ=255)
                                                       aTGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATc                          <  1:746485/45‑1 (MQ=255)
                                                       aTGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGa                      >  1:2274881/1‑49 (MQ=255)
                                                       aTGCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTa  <  1:2166746/69‑1 (MQ=255)
                                                         gCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCa                        >  1:1213486/1‑45 (MQ=255)
                                                         gCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCa                        >  1:1119045/1‑45 (MQ=255)
                                                         gCAATAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGAcc                    <  1:2118023/49‑1 (MQ=255)
                                                             |                                                              
GATTGATGAAACTGGCGTGACGCCTGTGATAAACCAGATCGAACTTCATCCGCTGATGCAACAACGCCAGCTACACGCCTGGAACGCGACACACAAAATCCAGACCGAATCCTGGAGCCCATTA  >  minE/2019826‑2019949

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