Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,558,446 T→C 40.0% G211G (GGT→GGC gpsA → glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,558,4460TC40.0% 24.6 / 44.1 65G211G (GGT→GGCgpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base C (14/12);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTG  >  minE/2558381‑2558504
                                                                 |                                                          
cTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCg                                                         <  1:662853/69‑1 (MQ=255)
  tGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCg                                                         <  1:34887/67‑1 (MQ=255)
                                 ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:745667/42‑1 (MQ=255)
                                 ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:2640090/42‑1 (MQ=255)
                                 ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:1408664/42‑1 (MQ=255)
                                 ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:2542503/42‑1 (MQ=255)
                                 ggTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:1674984/42‑1 (MQ=255)
                                 ggTGCGGGGATGGCCGACGGTATCGGTTTTGGTGCGAATgcg                                                   <  1:124379/42‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:582777/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:3039309/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:3253003/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:2488369/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:456660/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:465402/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:1587031/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:1281274/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:669609/68‑1 (MQ=255)
                                       gggATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGa                   <  1:1227560/68‑1 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:2286496/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:2318988/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:973152/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:3155020/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:3304575/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:386533/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:442720/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:583065/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1155232/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1723033/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1815391/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1188448/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1882745/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:144655/1‑69 (MQ=255)
                                             tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg            >  1:1424741/1‑69 (MQ=255)
                                                 aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc        >  1:1818865/1‑69 (MQ=255)
                                                 aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc        >  1:3070462/1‑69 (MQ=255)
                                                 aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc        >  1:1863774/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)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGgctggctg                    >  1:2590880/1‑51 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGt         >  1:1263267/1‑62 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:3243706/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:1330267/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:467302/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:1678912/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:70263/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:1697462/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:789875/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:855101/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:1764783/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:1167557/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:1773304/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:1801138/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:2233943/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:3045790/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:1849929/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:301407/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:3002511/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:2646968/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:2538764/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:2528285/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:2370741/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:2067837/69‑3 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  >  1:2315955/1‑67 (MQ=255)
                                                       tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg  <  1:2090468/69‑3 (MQ=255)
                                                                 |                                                          
CTTGGCGGCGCGGTGAAAAACGTTATTGCCATTGGTGCGGGGATGTCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTG  >  minE/2558381‑2558504

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: