Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,798,738 C→T 9.9% L253L (CTC→CTT proV → glycine betaine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,798,7380CT 9.9% 223.2 / 7.4 109L253L (CTC→CTTproVglycine betaine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (47/51);  new base T (9/2);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 5.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 113 total aligned reads displayed.

TTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAG  >  minE/1798671‑1798803
                                                                   |                                                                 
ttGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATtctc                                                                   >  1:2137611/1‑68 (MQ=255)
ttGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATtctc                                                                   >  1:2334919/1‑68 (MQ=255)
  ggCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaat                                                                >  1:2006415/1‑69 (MQ=255)
  ggCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaat                                                                >  1:1446099/1‑69 (MQ=255)
  ggCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaat                                                                >  1:1797937/1‑69 (MQ=255)
  ggCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaat                                                                >  1:1196595/1‑69 (MQ=255)
     gACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaataat                                                             <  1:2372273/69‑1 (MQ=255)
     gACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaataat                                                             <  1:2805670/69‑1 (MQ=255)
     gACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaataat                                                             <  1:3199384/69‑1 (MQ=255)
     gACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaataat                                                             <  1:1849851/69‑1 (MQ=255)
           aTTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCg                                                       <  1:1637434/69‑1 (MQ=255)
           aTTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCg                                                       <  1:1248709/69‑1 (MQ=255)
           aTTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCg                                                       <  1:736476/69‑1 (MQ=255)
                 tttATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3062687/66‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:2577422/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:2470118/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:1504006/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:2474053/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:1060579/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:684965/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3290055/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3066600/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:132832/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:2884445/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3185744/65‑1 (MQ=255)
                   tATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:2946114/65‑1 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaataat                                                             >  1:2412195/1‑54 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGaa                                                     >  1:2043969/1‑62 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTa                                               >  1:1673049/1‑68 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTa                                               >  1:2386065/1‑68 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:379313/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:415760/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:2981870/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:2084272/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:2472929/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:2786421/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:62744/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:744454/1‑69 (MQ=255)
                    aTGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                              >  1:1712694/1‑69 (MQ=255)
                           aTGGTGAAGTGGTACAGGTCGGCACCCCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                       >  1:1490224/1‑69 (MQ=255)
                           aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                       >  1:200454/1‑69 (MQ=255)
                           aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                       >  1:1585516/1‑69 (MQ=255)
                           aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                       >  1:2098361/1‑69 (MQ=255)
                           aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                       >  1:1108219/1‑69 (MQ=255)
                              gTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3055171/54‑1 (MQ=255)
                               tGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:1500294/53‑1 (MQ=255)
                                   gTGGTACAGGTTGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:2208047/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGt                                        >  1:338126/1‑60 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACcttctt                                >  1:1379168/1‑68 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:3051825/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:2746942/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:2308075/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:2930697/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:2581030/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:1114292/1‑69 (MQ=255)
                                   gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                               >  1:1554601/1‑69 (MQ=255)
                                     ggTACAGGTCGGCACACCGGATGAAATTCTCaataaa                                                             >  1:1819377/1‑36 (MQ=255)
                                     ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                             >  1:2123881/1‑69 (MQ=255)
                                     ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                             >  1:476301/1‑69 (MQ=255)
                                     ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                             >  1:165504/1‑69 (MQ=255)
                                       tACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATg                                                   <  1:3355812/45‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1281284/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1116361/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1983349/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1427428/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1619030/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1922830/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1676144/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:2073258/69‑1 (MQ=255)
                                              cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                    <  1:1762318/69‑1 (MQ=255)
                                                       ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                             >  1:3124045/1‑51 (MQ=255)
                                                       ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa           >  1:3127500/1‑69 (MQ=255)
                                                       ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa           >  1:2103837/1‑69 (MQ=255)
                                                       ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa           >  1:3006026/1‑69 (MQ=255)
                                                       ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa           >  1:2914025/1‑69 (MQ=255)
                                                         aTGAAATTCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGtc                      <  1:2961118/56‑2 (MQ=39)
                                                          tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                        <  1:2161050/53‑1 (MQ=255)
                                                          tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                        <  1:1160785/53‑1 (MQ=255)
                                                          tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                        <  1:1363036/53‑1 (MQ=255)
                                                          tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                        <  1:2476254/53‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:2336749/69‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:3056952/69‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:137500/69‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:1453018/69‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:1527549/69‑1 (MQ=255)
                                                           gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:2269534/69‑1 (MQ=255)
                                                            aaaTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:3179727/68‑1 (MQ=255)
                                                            aaaTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa       <  1:2565709/68‑1 (MQ=255)
                                                              aTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTc    <  1:779985/69‑1 (MQ=255)
                                                              aTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTc    <  1:596212/69‑1 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTTGCGTCGATATTAGCCAGGTGTTCa   >  1:405852/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:389718/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:1009403/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:249475/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:3350898/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:3347082/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:1391489/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:1078269/1‑69 (MQ=255)
                                                               ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa   >  1:2862477/1‑69 (MQ=255)
                                                                tctcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAg  >  1:2840694/1‑69 (MQ=255)
                                                                   |                                                                 
TTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAG  >  minE/1798671‑1798803

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