Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,798,765 T→A 13.5% R262R (CGT→CGA proV → glycine betaine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,798,7650TA13.5% 167.5 / 13.5 80R262R (CGT→CGAproVglycine betaine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base A (9/2);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCG  >  minE/1798698‑1798831
                                                                   |                                                                  
aTGGTGAAGTGGTACAGGTCGGCACCCCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                                                   >  1:1490224/1‑69 (MQ=255)
aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                                                   >  1:1108219/1‑69 (MQ=255)
aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                                                   >  1:2098361/1‑69 (MQ=255)
aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                                                   >  1:200454/1‑69 (MQ=255)
aTGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTa                                                                   >  1:1585516/1‑69 (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:2943286/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:2308075/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:2413363/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:1554601/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:2581030/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:1440763/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:2746942/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:1114292/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:2930697/1‑69 (MQ=255)
        gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCttcttc                                                           >  1:3051825/1‑69 (MQ=255)
          ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                                                         >  1:476301/1‑69 (MQ=255)
          ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                                                         >  1:165504/1‑69 (MQ=255)
          ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                                                         >  1:2123881/1‑69 (MQ=255)
          ggTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCg                                                         >  1:431225/1‑69 (MQ=255)
                 gTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGtt                                                  <  1:1772256/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1281284/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1427428/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1617855/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1619030/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1676144/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1762318/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1922830/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1983349/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:1116361/69‑1 (MQ=255)
                   cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                <  1:2073258/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:2914025/1‑69 (MQ=255)
                            ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa                                       >  1:3006026/1‑69 (MQ=255)
                            ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCa                                       >  1:2736172/1‑69 (MQ=255)
                              aTGAAATTCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGtc                                                  <  1:2961118/56‑2 (MQ=39)
                               tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                                                    <  1:1160785/53‑1 (MQ=255)
                               tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                                                    <  1:2161050/53‑1 (MQ=255)
                               tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                                                    <  1:2476254/53‑1 (MQ=255)
                               tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                                                    <  1:1363036/53‑1 (MQ=255)
                                gAAATTCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGtc                                                  <  1:3314869/54‑2 (MQ=38)
                                gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa                                   <  1:2265536/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:3059046/69‑1 (MQ=255)
                                gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa                                   <  1:2336749/69‑1 (MQ=255)
                                gAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa                                   <  1:2269534/69‑1 (MQ=255)
                                 aaaTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa                                   <  1:2565709/68‑1 (MQ=255)
                                 aaaTTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTa                                   <  1:3179727/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:2862477/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:1009403/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:249475/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:1078269/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:1391489/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:389718/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:3347082/1‑69 (MQ=255)
                                    ttCTTAATAATCCGGCGAATGATTATGTCCGAACCTTCTTCCGTGGCGTCGATATTAGCCAGGTGTTCa                               >  1:3350898/1‑69 (MQ=255)
                                     tctcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAg                              >  1:2840694/1‑69 (MQ=255)
                                     tctcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAg                              >  1:818237/1‑69 (MQ=255)
                                               ccGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGa                     <  1:1286228/68‑1 (MQ=255)
                                               ccGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGa                     <  1:1516646/68‑1 (MQ=255)
                                               ccGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGa                     <  1:2838222/68‑1 (MQ=255)
                                                              gTCCGTACCTTCTTCCGTGTCGTTGATATTAGTCAGGTa                                   <  1:534572/39‑1 (MQ=255)
                                                               tCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGccc             >  1:107628/1‑60 (MQ=255)
                                                               tCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGaca     >  1:1348116/1‑68 (MQ=255)
                                                               tCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGaca     >  1:801957/1‑68 (MQ=255)
                                                               tCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGaca     >  1:446964/1‑68 (MQ=255)
                                                               tCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGac      >  1:93251/1‑67 (MQ=255)
                                                               tCCGTACCTTCTTCAGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGaca     >  1:1249854/1‑68 (MQ=255)
                                                               tCCGTAACTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCg                          >  1:441872/1‑47 (MQ=255)
                                                                 cGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac    >  1:572458/1‑67 (MQ=255)
                                                                 cGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCg  >  1:1031028/1‑69 (MQ=255)
                                                                   |                                                                  
ATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCG  >  minE/1798698‑1798831

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: