Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,404,986 T→G 25.2% T825T (ACA→ACC cyaA ← adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,404,9860TG25.2% 137.8 / 54.8 103T825T (ACA→ACCcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (34/43);  new base G (13/13);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGA  >  minE/2404930‑2405053
                                                        |                                                                   
tGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          <  1:1668857/68‑1 (MQ=255)
tGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          <  1:2706517/68‑1 (MQ=255)
      aGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCAc                                                    <  1:3068113/68‑1 (MQ=255)
      aGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCAc                                                    <  1:1946164/68‑1 (MQ=255)
              aTCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACg                                           <  1:1198368/69‑1 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcgga                                                              >  1:2608606/1‑49 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:476195/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:238339/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:2624485/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:1571035/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:32044/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:1505585/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:651183/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:931097/1‑64 (MQ=255)
               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCa                                             >  1:647273/1‑66 (MQ=255)
                 gTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGAcc                                        >  1:2751557/1‑69 (MQ=255)
                    aTCCTGATTGGCAGTCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATc                                     <  1:2532139/69‑1 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:2131739/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:3132327/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:3176585/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:2790597/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1547821/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1524797/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:604976/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1221195/1‑48 (MQ=255)
                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1163715/1‑48 (MQ=255)
                    aTCCTGATTGGAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:567355/1‑47 (MQ=255)
                      gcTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:243939/67‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2727753/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2482668/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2038116/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1868656/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1856178/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1707244/68‑1 (MQ=255)
                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:321248/68‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:2123791/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:303901/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3043317/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:2745012/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3189284/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3223547/57‑1 (MQ=255)
                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:849830/57‑1 (MQ=255)
                           ttGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCtt                               >  1:103609/1‑68 (MQ=255)
                           ttGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCtt                               >  1:2707007/1‑68 (MQ=255)
                            tGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTc                              >  1:3269990/1‑68 (MQ=255)
                            tGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTc                              >  1:3149616/1‑68 (MQ=255)
                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:787507/68‑1 (MQ=255)
                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:531540/68‑1 (MQ=255)
                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:2194728/68‑1 (MQ=255)
                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:3151045/69‑1 (MQ=255)
                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:1513856/69‑1 (MQ=255)
                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:2140016/69‑1 (MQ=255)
                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:859666/69‑1 (MQ=255)
                                   cggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:1438289/68‑1 (MQ=255)
                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:832882/1‑67 (MQ=255)
                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:336464/1‑67 (MQ=255)
                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:2373654/1‑67 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3084391/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2409409/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3389026/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:42716/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3127075/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3080485/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:1311517/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2983003/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2925019/48‑1 (MQ=255)
                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTGCACGAc                                         <  1:8072/48‑1 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCTAAACGGAATCACCTGTTCACGAc                                         <  1:2711608/45‑1 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGt                                                >  1:1211266/1‑38 (MQ=38)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:270163/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:640615/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:374941/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:70361/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2504626/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1551112/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2202101/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2877556/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1948967/1‑68 (MQ=255)
                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1768686/1‑68 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:1400339/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:3124877/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:34046/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:210785/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:1253834/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:650246/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:80633/66‑1 (MQ=255)
                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:785718/66‑1 (MQ=255)
                                              cGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGc            <  1:1310753/68‑1 (MQ=255)
                                              cGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGc            <  1:736466/68‑1 (MQ=255)
                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTg                      >  1:1147242/1‑55 (MQ=255)
                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1061324/1‑59 (MQ=255)
                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:85542/1‑59 (MQ=255)
                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1629389/1‑59 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  <  1:800696/69‑1 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:1939650/1‑69 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  <  1:2183908/69‑1 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2270563/1‑69 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2287194/1‑69 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2453555/1‑69 (MQ=255)
                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2471830/1‑69 (MQ=255)
                                                        |                                                                   
TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGA  >  minE/2404930‑2405053

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: