Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,555,847 G→A 31.6% G123G (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,8470GA31.6% 73.5 / 39.1 82G123G (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 = 9.38e-01

TTTTTTCTGCGCCGCGTTCAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACCGC  >  minE/2555780‑2555912
                                                                   |                                                                 
ttttttCTGCGCCGCGTTCAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTAcgc                                                                  >  1:1574367/1‑69 (MQ=255)
 tttttCTGCGCCGCGTTCAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTAcgcc                                                                 <  1:1261689/69‑1 (MQ=255)
  ttttCTGCGCCGCGTTCAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTAcgcc                                                                 <  1:2525552/68‑1 (MQ=255)
                  cAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGa                                                 >  1:2887250/1‑68 (MQ=255)
                  cAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGa                                                 >  1:2125340/1‑68 (MQ=255)
                             tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCa                                             >  1:926177/1‑61 (MQ=255)
                             tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATgtgt                                     >  1:922227/1‑69 (MQ=255)
                             tCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATgtgt                                     >  1:2790440/1‑69 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGCCCCATAATGTGTa                                    >  1:478501/1‑69 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:540663/61‑1 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:525777/61‑1 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:315713/61‑1 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:136349/61‑1 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:3105184/61‑1 (MQ=255)
                              cTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            <  1:2566612/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:2474193/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:2269024/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:2557456/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:1601251/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:1251710/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:822425/1‑69 (MQ=255)
                                   aTCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCt                               >  1:912633/1‑69 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCGTGCGGAGAGCAGACCCAt                                            >  1:2819594/1‑53 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1432991/1‑53 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1728944/1‑53 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:1907522/1‑53 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:2755628/1‑53 (MQ=255)
                                      gCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCAt                                            >  1:956432/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:1033131/62‑1 (MQ=255)
                                            atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1488432/62‑1 (MQ=255)
                                            atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1479397/62‑1 (MQ=255)
                                            atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1718493/62‑1 (MQ=255)
                                            atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:1208095/62‑1 (MQ=255)
                                            atgTGATCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCttt                             <  1:2333540/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:1699349/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:2326685/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:814089/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2941956/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:1237611/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2893992/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:289063/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:835176/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:1918919/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:2789204/1‑69 (MQ=255)
                                                  tCTTCGTGGCTGTGTACACCACCTGCGGAGAGCAGACCCATAATGTGTACTGCTTTGCCTGCGTTTTTc                >  1:507142/1‑69 (MQ=255)
                                                    ttCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              >  1:1714229/1‑69 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:2490204/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:668923/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1317605/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1221931/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1211206/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:2005982/65‑1 (MQ=255)
                                                        tGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGc              <  1:1616458/65‑1 (MQ=255)
                                                                tACGTCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:454535/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGttt                   >  1:99796/1‑52 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGttt                   >  1:2128197/1‑52 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1350300/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:682528/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:3259122/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:748011/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:177931/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2904147/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1879360/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:444996/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2281709/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:1188843/1‑69 (MQ=255)
                                                                tACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACcgc  >  1:2467966/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)
                                                                   |                                                                 
TTTTTTCTGCGCCGCGTTCAGCTGCCAGTTCTACCATCGCCATGATGTGATCTTCGTGGCTGTGTACGCCGCCTGCGGAGAGCAGACCCATAATGTGTACCGCTTTGCCTGCGTTTTTCGCTTTATCTACCGC  >  minE/2555780‑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: