Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,622,426 C→T 24.2% V740V (GTG→GTA maeB ← fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,622,4260CT24.2% 99.9 / 24.6 78V740V (GTG→GTAmaeBfused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base T (6/13);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TACAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGTCAC  >  minE/1622367‑1622494
                                                           |                                                                    
tACAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCg                                                             <  1:2238232/69‑1 (MQ=255)
tACAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCg                                                             <  1:634972/69‑1 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:1777530/1‑69 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:2373330/1‑69 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:34851/1‑69 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:135820/1‑69 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:1529833/1‑69 (MQ=255)
 aCAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGa                                                            >  1:1450823/1‑69 (MQ=255)
    gCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCg                                                         <  1:2953735/69‑1 (MQ=255)
              ttGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAacac                                               <  1:1107391/69‑1 (MQ=255)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:296159/1‑46 (MQ=39)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:1748883/1‑46 (MQ=39)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:3292668/1‑46 (MQ=39)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:3417998/1‑46 (MQ=39)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:2448068/1‑46 (MQ=39)
                gcgcTTCAACCACGGCCAGCGCCACCATGTTGACGATACGACGTAc                                                                    >  1:892585/1‑46 (MQ=39)
                     tCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCg                                                         <  1:2516679/52‑1 (MQ=255)
                                cAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCcac                             <  1:2742897/69‑1 (MQ=255)
                                  gcgcCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCacac                           <  1:297970/69‑1 (MQ=255)
                                  gcgcCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCacac                           <  1:3286472/69‑1 (MQ=255)
                                  gcgcCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCacac                           <  1:3227803/69‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCTATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:3392074/53‑1 (MQ=39)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACAGGTGAACCGGt                                     <  1:1054396/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:1676554/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:249947/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:2837491/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:2918916/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:382341/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:379491/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:1379258/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:2031975/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:3189831/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:2201188/53‑1 (MQ=255)
                                        ccaTGTTGACGATACGACGTACCGATGCGATCGGCGTTAACACGTGAACCGGt                                     <  1:2144389/53‑1 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTc                                  >  1:2057828/1‑51 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACAcc                          >  1:3026636/1‑59 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACGCATCAGCAcc                >  1:2267017/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:310570/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:559221/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:3123638/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:796289/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:2474289/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:10278/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:1331279/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:1309065/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:1136877/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:2114805/1‑69 (MQ=255)
                                             ttGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCAcc                >  1:2234089/1‑69 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGATCCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2106054/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:1559076/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:3307053/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:397670/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:1280040/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:202385/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2328693/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:734902/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2385991/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:854538/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2825470/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2464116/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:1111328/1‑67 (MQ=255)
                                                    tACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGcc           >  1:2541673/1‑67 (MQ=255)
                                                     acgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCGGAc        >  1:721709/1‑69 (MQ=255)
                                                      cgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACa       <  1:131038/69‑1 (MQ=255)
                                                      cgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACa       <  1:997050/69‑1 (MQ=255)
                                                      cgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACa       <  1:1811306/69‑1 (MQ=255)
                                                      cgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACa       <  1:2958535/69‑1 (MQ=255)
                                                      cgacgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACa       <  1:2389116/69‑1 (MQ=255)
                                                        acgCACCTATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:2624692/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:1384657/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:1700521/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:3278127/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:199850/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:121577/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:660323/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:2275939/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:2451/69‑1 (MQ=255)
                                                        acgCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGt     <  1:2486792/69‑1 (MQ=255)
                                                           cACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGTcac  <  1:3319379/69‑1 (MQ=255)
                                                           cACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGTcac  <  1:556850/69‑1 (MQ=255)
                                                           cACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGTcac  <  1:938154/69‑1 (MQ=255)
                                                           |                                                                    
TACAGCGGTTGGGTTTGCGCTTCTACCACGGCCAGCGCCACCATGTTGACGATACGACGCACCGATGCGATCGGCGTTAACACGTGAACCGGTTTCGCCACACCCATCAGCACCGGGCCGACAGTCAC  >  minE/1622367‑1622494

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: