Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,555,850 G→A 31.6% G122G (GGC→GGT gpmI ← phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,555,8500GA31.6% 72.7 / 39.0 82G122G (GGC→GGTgpmIphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base G (42/14);  new base A (16/10);  total (58/24)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

CAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACCGC  >  minE/2555798‑2555912
                                                    |                                                              
cAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGa                                                 >  1:2887250/1‑68 (MQ=255)
cAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGa                                                 >  1:2125340/1‑68 (MQ=255)
           tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCa                                             >  1:926177/1‑61 (MQ=255)
           tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATgtgt                                     >  1:2790440/1‑69 (MQ=255)
           tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATgtgt                                     >  1:922227/1‑69 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGCCCCATAATGTGTa                                    >  1:478501/1‑69 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:540663/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:136349/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:2566612/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:525777/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:3105184/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:315713/61‑1 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAAt                                         >  1:2504812/1‑64 (MQ=255)
            cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTa                                    >  1:1486650/1‑69 (MQ=255)
             tACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:454740/60‑1 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:2269024/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:2474193/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:912633/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:822425/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:1601251/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:2557456/1‑69 (MQ=255)
                 aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:1251710/1‑69 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCGTGCGGAGAGCAGACCCAt                                            >  1:2819594/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1728944/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1907522/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:956432/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:2755628/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1432991/1‑53 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTg                            >  1:3425106/1‑69 (MQ=255)
                    gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTg                            >  1:2132892/1‑69 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCg                      >  1:984221/1‑69 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCg                      >  1:1435065/1‑69 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCg                      >  1:75220/1‑69 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCATACCCATAATGTGTACTGCttt                             <  1:3203661/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:2560981/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:37967/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1057093/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:2333540/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1718493/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1488432/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1479397/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1208095/62‑1 (MQ=255)
                          atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1033131/62‑1 (MQ=255)
                          atgTGATCATCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCg                      >  1:1836677/1‑69 (MQ=255)
                           tgtgATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGt                     >  1:1987068/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTa                                    >  1:741422/1‑49 (MQ=39)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:835176/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:1237611/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:814089/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:507142/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:3395778/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:3223933/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:3176310/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:1699349/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2941956/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2326685/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:1918919/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2893992/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2789204/1‑69 (MQ=255)
                                tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:289063/1‑69 (MQ=255)
                                  ttCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              >  1:1714229/1‑69 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1616458/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:668923/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1317605/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:2005982/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1221931/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1211206/65‑1 (MQ=255)
                                      tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:2490204/65‑1 (MQ=255)
                                              tACGTCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:454535/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGttt                   >  1:2128197/1‑52 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGttt                   >  1:99796/1‑52 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1350300/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:682528/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1879360/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:748011/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:444996/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2281709/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2467966/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:3259122/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2904147/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1188843/1‑69 (MQ=255)
                                              tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:177931/1‑69 (MQ=255)
                                                cgccgcCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTAt         >  1:1569668/1‑60 (MQ=255)
                                                cgccgcCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:3147040/1‑67 (MQ=255)
                                                cgccgcCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcg   >  1:438403/1‑66 (MQ=255)
                                                    |                                                              
CAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACCGC  >  minE/2555798‑2555912

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: