Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,224,379 G→A 13.7% E81K (GAA→AAA)  yebU → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,224,3790GA13.7% 203.6 / 20.6 73E81K (GAA→AAA) yebUpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base A (4/6);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

GCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTT  >  minE/1224310‑1224448
                                                                     |                                                                     
gCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATg                                                                       >  1:1902130/1‑70 (MQ=255)
gCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGa                                                                      <  1:3851518/71‑1 (MQ=255)
gCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGa                                                                      <  1:2391948/71‑1 (MQ=255)
  tCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGa                                                                      <  1:3445765/69‑1 (MQ=255)
    cTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGAt                                                                  >  1:3865326/1‑71 (MQ=255)
     ttATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATg                                                                 >  1:932772/1‑71 (MQ=255)
       aTGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATg                                                                 <  1:3223874/69‑1 (MQ=255)
       aTGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCa                                                               >  1:1012465/1‑71 (MQ=255)
       aTGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCa                                                               >  1:1482530/1‑71 (MQ=255)
        tGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCAt                                                              >  1:1979609/1‑71 (MQ=255)
        tGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCAt                                                              >  1:806245/1‑71 (MQ=255)
        tGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCAt                                                              <  1:1091330/71‑1 (MQ=255)
              gACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCa                                                         <  1:1075616/70‑1 (MQ=255)
              gACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCAt                                                        <  1:2717784/71‑1 (MQ=255)
              gACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCAt                                                        <  1:2026093/71‑1 (MQ=255)
                cGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCAtt                                                       >  1:3116829/1‑70 (MQ=255)
                cGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTg                                                      >  1:566602/1‑71 (MQ=255)
                 gCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTgg                                                     >  1:1213104/1‑71 (MQ=255)
                    tACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTggg                                                    >  1:1428429/1‑69 (MQ=255)
                     aCGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 >  1:3171790/1‑71 (MQ=255)
                      cGCCAATTCCGTGGTGTGAGGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 <  1:3527718/70‑1 (MQ=255)
                          aaTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGc                                                                >  1:3052230/1‑51 (MQ=255)
                          aaTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGc                                                                >  1:208253/1‑51 (MQ=255)
                          aaTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGc                                                                >  1:2709742/1‑51 (MQ=255)
                          aaTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGc                                                                >  1:156468/1‑51 (MQ=255)
                          aaTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTAccg                                            <  1:26125/71‑1 (MQ=255)
                           aTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGaa                                                                     <  1:1767240/45‑1 (MQ=255)
                           aTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGaa                                                                     <  1:1753012/45‑1 (MQ=255)
                           aTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGaa                                                                     <  1:565314/45‑1 (MQ=255)
                           aTTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAccgc                                           >  1:3063471/1‑71 (MQ=255)
                            ttCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAccgcc                                          >  1:2788995/1‑71 (MQ=255)
                             tCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGa                                                                   >  1:4029756/1‑45 (MQ=255)
                             tCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGa                                                                   >  1:214965/1‑45 (MQ=255)
                                gtggtgTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGc                                      <  1:3418053/71‑1 (MQ=255)
                                gtggtgTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTACCGCCGAGc                                      <  1:2161829/71‑1 (MQ=255)
                                gtggtgTGAAGAAGGTTTCTGGATTGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTACCGCCGAGc                                      <  1:1292011/71‑1 (MQ=255)
                                   gtgtGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCa                                                               >  1:1924687/1‑43 (MQ=255)
                                    tgtgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAcc                                             <  1:983904/60‑1 (MQ=255)
                                     gtgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCa                                     >  1:310547/1‑67 (MQ=255)
                                     gtgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCa                                     >  1:372187/1‑67 (MQ=255)
                                      tgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTg                                                      <  1:989500/49‑1 (MQ=255)
                                      tgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTg                                                      <  1:283490/49‑1 (MQ=255)
                                      tgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTg                                                      <  1:1945323/49‑1 (MQ=255)
                                      tgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGTATTTa                                 <  1:2066205/70‑1 (MQ=255)
                                      tgAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCAttt                                  <  1:774051/69‑1 (MQ=255)
                                       gaagaaGGTTTCTGGATTGAACGCGACAATGAAGATGCAtt                                                             <  1:3393685/41‑1 (MQ=255)
                                        aagaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 >  1:1870051/1‑52 (MQ=255)
                                        aagaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 >  1:2593084/1‑52 (MQ=255)
                                        aagaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 >  1:2510167/1‑52 (MQ=255)
                                        aagaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtag                                                 >  1:420646/1‑52 (MQ=255)
                                        aagaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGt                              >  1:620943/1‑71 (MQ=255)
                                         agaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGc                                                           >  1:4032159/1‑41 (MQ=255)
                                         agaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGc                                                           >  1:3164932/1‑41 (MQ=255)
                                         agaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGt                              <  1:3501187/70‑1 (MQ=255)
                                          gaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCAttt                                  <  1:343946/65‑1 (MQ=255)
                                          gaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCAttt                                  <  1:436669/65‑1 (MQ=255)
                                          gaagGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCAttt                                  <  1:3073388/65‑1 (MQ=255)
                                           aagGTTTCTGGATTGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTAccgcc                                          <  1:1527379/56‑1 (MQ=255)
                                               tttCTGGATTGAACGCGACAATGAAGATGCATTGCCAtt                                                       <  1:4103906/39‑1 (MQ=255)
                                                 tCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGttt                     >  1:3806406/1‑71 (MQ=255)
                                                  cTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGttt                     <  1:1990605/70‑1 (MQ=255)
                                                   tGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTa                                 <  1:3984056/57‑1 (MQ=255)
                                                     gATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTta                   <  1:1188475/69‑1 (MQ=255)
                                                     gATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTta                   <  1:2464165/69‑1 (MQ=255)
                                                     gATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTta                   <  1:782196/69‑1 (MQ=255)
                                                        tGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATAtt               >  1:2111407/1‑70 (MQ=255)
                                                        tGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATAtt               >  1:1296230/1‑70 (MQ=255)
                                                        tGAACGCGACAATAAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATAtt               >  1:2752825/1‑70 (MQ=255)
                                                         gAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAccgccg                                         >  1:2930351/1‑43 (MQ=255)
                                                             gcgACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAgg           >  1:2934827/1‑69 (MQ=255)
                                                             gcgACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAgg           >  1:1033251/1‑69 (MQ=255)
                                                              cgACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTatat                <  1:3372845/63‑1 (MQ=255)
                                                              cgACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAg        >  1:839940/1‑71 (MQ=255)
                                                                    tGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGtt  <  1:2865985/71‑1 (MQ=255)
                                                                     |                                                                     
GCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTT  >  minE/1224310‑1224448

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: