Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,558,488 A→G 32.0% E225E (GAA→GAG gpsA → glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,558,4880AG32.0% 63.1 / 39.4 78E225E (GAA→GAGgpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (38/15);  new base G (13/12);  total (51/27)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGG  >  minE/2558426‑2558547
                                                              |                                                           
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:2318988/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:973152/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:583065/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:442720/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:386533/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:3304575/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:3155020/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:2286496/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1882745/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1815391/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1723033/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:144655/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1155232/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1424741/1‑69 (MQ=255)
tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg                                                       >  1:1188448/1‑69 (MQ=255)
    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc                                                   >  1:3070462/1‑69 (MQ=255)
    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc                                                   >  1:1863774/1‑69 (MQ=255)
    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc                                                   >  1:1818865/1‑69 (MQ=255)
    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTAGCGTc                                                   >  1:358289/1‑69 (MQ=255)
    aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGATGAAATGTCGCGTc                                                   >  1:1887143/1‑69 (MQ=255)
         aTCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGt                                              >  1:2958844/1‑69 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGt                                                    >  1:1263267/1‑62 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:301407/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:2315955/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:2370741/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:2528285/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:2646968/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:3002511/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:855101/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:467302/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:789875/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:3045790/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:2538764/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:3243706/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:70263/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:1167557/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:2067837/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:1773304/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:1801138/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:1330267/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:1849929/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:1764783/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:1697462/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             <  1:2090468/69‑3 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:2233943/1‑67 (MQ=255)
          tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg                                             >  1:1678912/1‑67 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:2945863/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:1456102/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:1323265/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:1105895/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:304039/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:678659/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:1897895/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:2858010/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             <  1:1716311/69‑1 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:1699790/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:2578558/1‑69 (MQ=255)
                          tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa                             >  1:662648/1‑69 (MQ=255)
                           gcgcGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGAc                            <  1:3304314/69‑1 (MQ=255)
                           gcgcGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGAc                            <  1:1592503/69‑1 (MQ=255)
                                 aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc                      <  1:1420926/69‑1 (MQ=255)
                                 aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc                      <  1:2388427/69‑1 (MQ=255)
                                 aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc                      <  1:2478688/69‑1 (MQ=255)
                                    gcgcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCAcc                   >  1:2659863/1‑69 (MQ=255)
                                    gcgcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCAcc                   >  1:3400125/1‑69 (MQ=255)
                                      gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt                 <  1:2608601/69‑1 (MQ=255)
                                      gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt                 <  1:2339717/69‑1 (MQ=255)
                                      gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt                 <  1:881709/69‑1 (MQ=255)
                                       cTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt                 <  1:225322/68‑1 (MQ=255)
                                                cGTGGGCTGGCTTAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg       >  1:3370460/1‑69 (MQ=255)
                                                cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCAt        >  1:2648440/1‑68 (MQ=255)
                                                cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg       >  1:1214385/1‑69 (MQ=255)
                                                cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg       >  1:1427865/1‑69 (MQ=255)
                                                      ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg  <  1:1612056/68‑1 (MQ=255)
                                                      ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg  <  1:1587725/68‑1 (MQ=255)
                                                      ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg  <  1:2301751/68‑1 (MQ=255)
                                                       tggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg  <  1:1406831/67‑1 (MQ=255)
                                                           tgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg  <  1:754775/63‑1 (MQ=255)
                                                              |                                                           
TCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGG  >  minE/2558426‑2558547

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: