Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,682,068 A→G 38.1% V969A (GTG→GCG)  yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,682,0680AG38.1% 55.2 / 69.5 105V969A (GTG→GCG) yfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/40);  new base G (24/16);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-02

Only 100 of 105 total aligned reads displayed.

GCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCAC  >  minE/1682004‑1682122
                                                                |                                                      
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:2286254/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:1981859/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:3175825/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:587748/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:734382/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:186881/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:2113659/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:2122726/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:785354/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:1425819/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:3325467/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:331349/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:1191081/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:447262/69‑1 (MQ=255)
gCCAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:1222846/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:3418734/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:1162639/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:1868801/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:1856670/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:1350059/69‑1 (MQ=255)
gCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAc                                                    <  1:766335/69‑1 (MQ=255)
 ccAATTGGCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGAc                                                    <  1:217824/68‑1 (MQ=255)
  cAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACaa                                                  <  1:2395666/69‑1 (MQ=255)
  cAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACaa                                                  <  1:2988268/69‑1 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:2272582/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:31148/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:272667/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:56315/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:2455394/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:2006159/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:2851411/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:1999259/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:1799897/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:3276041/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAata                                                >  1:1155144/3‑68 (MQ=255)
     ttggCAGTGTAACCGAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAat                                                 >  1:1905426/3‑67 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:312276/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:962471/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:3316970/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:1360707/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:3300227/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:3111897/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:3246139/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:1627865/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:740418/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:717767/69‑1 (MQ=255)
           gtgtAACCGAGCCTTCGACCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:1906723/69‑1 (MQ=255)
            tgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGGCTGCTGCACGACAATATTGACGt                                         <  1:2213315/68‑1 (MQ=255)
            tgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:529428/68‑1 (MQ=255)
            tgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:769540/68‑1 (MQ=255)
            tgtAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:2484589/68‑1 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAatat                                               >  1:421539/1‑60 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:3168738/1‑67 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:2107193/1‑67 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:1075079/1‑67 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:953588/1‑67 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:1797827/1‑67 (MQ=255)
              tAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTg                                        >  1:2809154/1‑67 (MQ=255)
                aCCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTg                                    <  1:492243/69‑1 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACaa                                                  >  1:1586086/1‑52 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:667579/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:1344907/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:879937/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:1456291/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:609203/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:581672/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:3059795/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:2469486/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:2476055/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:2858981/1‑64 (MQ=255)
                   gAGCCTTCGCCCTGTTCATTGAGCGTTACCGGCAGCGCCTGCTGCGCGACAATATTGACGTGAt                                      >  1:2555819/1‑64 (MQ=255)
                      ccTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:1925855/58‑1 (MQ=255)
                             ccTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:981432/69‑1 (MQ=255)
                             ccTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:3421364/69‑1 (MQ=255)
                             ccTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:889435/69‑1 (MQ=255)
                             ccTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:1986259/69‑1 (MQ=255)
                             ccTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:1959678/69‑1 (MQ=255)
                              cTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTac                       <  1:1474594/68‑1 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:3244854/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:1174895/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:1267509/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:1372858/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:2574950/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:2511478/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:2407671/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccacc                   >  1:1519550/1‑69 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccac                    >  1:3319446/1‑68 (MQ=255)
                                 ttCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAAAATATTGACGt                                         >  1:1028295/1‑47 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:1941281/1‑45 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:1889510/1‑45 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:606360/1‑45 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:2907923/1‑45 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:2964691/1‑45 (MQ=255)
                                   cGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         >  1:3216376/1‑45 (MQ=255)
                                       gAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGt                                         <  1:1215284/41‑1 (MQ=255)
                                         gCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACTGGCGGTTTACCACCACGTTTc            <  1:3317094/68‑1 (MQ=255)
                                            ttACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCt        <  1:2047640/69‑1 (MQ=255)
                                                  gCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCAc  <  1:1811356/69‑1 (MQ=255)
                                                  gCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCAc  <  1:1509339/69‑1 (MQ=255)
                                                  gCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCAc  <  1:942947/69‑1 (MQ=255)
                                                                |                                                      
GCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCAC  >  minE/1682004‑1682122

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