Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,238,425 A→G 43.5% T32T (ACT→ACC edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,238,4250AG43.5% 13.8 / 53.3 69T32T (ACT→ACCedd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (25/14);  new base G (19/11);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCGAAACCGTGTGCCAGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGTCTCGCGCG  >  minE/1238369‑1238474
                                                        |                                                 
gcagCGAAACCGTGTGCCAGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGt                                       <  1:1395415/69‑1 (MQ=255)
             gtgCCAGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCgg                              >  1:1498241/1‑65 (MQ=255)
             gtgCCAGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCgg                              >  1:3365327/1‑65 (MQ=255)
                  aGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGGGAGATa                      <  1:3459/68‑1 (MQ=255)
                  aGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATa                      <  1:871073/68‑1 (MQ=255)
                  aGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATa                      <  1:2877330/68‑1 (MQ=255)
                      tACCGGATGCCAACTGAGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:2760698/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1060189/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:2195764/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1837184/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:2262676/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1748492/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1665673/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:3041987/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:3156561/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:3374211/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1629716/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:3417863/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:342942/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:368389/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1284040/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1151847/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:701109/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1131239/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAGGTTTTCGCTTGTTCTATc                                 >  1:1008576/1‑53 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCag                 >  1:373941/1‑69 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCag                 >  1:956171/1‑69 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCag                 >  1:2159755/1‑69 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCag                 >  1:2064023/1‑69 (MQ=255)
                      tACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCag                 >  1:1632467/1‑69 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:812179/69‑1 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:2213628/69‑1 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:1818946/69‑1 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:2382978/69‑1 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:3341196/69‑1 (MQ=255)
                            aTGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:1661228/69‑1 (MQ=255)
                             tGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAg           <  1:3051970/68‑1 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:3257711/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1065135/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1145512/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:48280/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:467466/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:409852/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:3227927/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:2977130/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1988991/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1873665/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1849443/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:2810939/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:1661476/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtct        >  1:301222/1‑69 (MQ=255)
                               ccAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGtc         >  1:813580/1‑68 (MQ=255)
                                    tGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTc                                     <  1:1574055/35‑1 (MQ=255)
                                    tGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTc                                     <  1:1336601/35‑1 (MQ=255)
                                    tGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTc                                     <  1:1654862/35‑1 (MQ=255)
                                     gCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGTCTcgcgcg  >  1:2055635/1‑69 (MQ=255)
                                     gCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGTCTcgcgcg  >  1:880344/1‑69 (MQ=255)
                                     gCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGTCTcgcgcg  >  1:3094379/1‑69 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:2548825/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1114819/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:2209792/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:2213931/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1653116/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1171209/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:371877/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:348606/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1425416/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1515228/43‑1 (MQ=255)
                                                 cGGTCGAGGTTTTCGCTTGTTCTATCCGGGCGAGGTAAGCaga                <  1:1747400/43‑1 (MQ=255)
                                                        |                                                 
GCAGCGAAACCGTGTGCCAGGTTACCGCATGCCAACTGCGAACGATGAACGGTCGAAGTTTTCGCTTGTTCTATCCGGGCGAGATAAGCAGAGCGAGTCTCGCGCG  >  minE/1238369‑1238474

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: