Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 41,783 G→A 29.8% A306A (GCG→GCA fixB → predicted electron transfer flavoprotein, NAD/FAD‑binding domain and ETFP adenine nucleotide‑binding domain‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE41,7830GA29.8% 86.5 / 38.4 87A306A (GCG→GCAfixBpredicted electron transfer flavoprotein, NAD/FAD‑binding domain and ETFP adenine nucleotide‑binding domain‑like
Reads supporting (aligned to +/- strand):  ref base G (44/17);  new base A (20/6);  total (64/23)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

AATGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCTGA  >  minE/41718‑41845
                                                                 |                                                              
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:3383659/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:68711/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:1781012/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:948441/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:1366517/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:1342692/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:3277789/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTg                                                             >  1:309959/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATACTTCCGGCGCTg                                                             >  1:1970532/1‑69 (MQ=255)
aaTGCGCCGATCTTCCAGTAAGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCGGGCGCTg                                                             >  1:3368097/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCGGCttt                                                  >  1:2065295/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGc                                                     >  1:1994262/1‑66 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCttt                                                  >  1:97628/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCttt                                                  >  1:2817911/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCttt                                                  >  1:3178346/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCttt                                                  >  1:2089086/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGAAGCttt                                                  >  1:2980801/1‑69 (MQ=255)
           cTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTACCGGCGCTGACCGCAGCttt                                                  >  1:559659/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:2541926/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:633605/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:2229156/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:617965/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:26335/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:2667937/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:2815889/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:1032300/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:1847353/1‑69 (MQ=255)
                    cgcgGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTg                                         >  1:155909/1‑69 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATc                                      <  1:2630492/69‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATc                                      <  1:3076643/69‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATc                                      <  1:284968/69‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATc                                      <  1:240733/69‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCTCAGCTTTAGCGCGTTGAt                                       <  1:1597436/68‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGAt                                       <  1:2640897/68‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGAt                                       <  1:2806198/68‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGAt                                       <  1:2818955/68‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGAt                                       <  1:3384603/68‑1 (MQ=255)
                       ggATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGAt                                       <  1:589713/68‑1 (MQ=255)
                          ttACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCAc                                   <  1:3369831/69‑1 (MQ=255)
                          ttACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCAc                                   <  1:1328052/69‑1 (MQ=255)
                          ttACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCAc                                   <  1:2180474/69‑1 (MQ=255)
                          ttACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCAc                                   <  1:2005427/69‑1 (MQ=255)
                            aCGGCATTGTTGGCGCCGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:790267/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGtt                                          >  1:886644/1‑60 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:2931326/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1227182/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1300206/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1201578/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1165684/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1089514/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:259421/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:2791091/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1368325/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:13883/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:2540740/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:410030/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:901882/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:876434/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1807144/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCActc                                 >  1:1976475/1‑69 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCAct                                  >  1:2569140/1‑68 (MQ=255)
                            aCGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCACTGACCGCAGCTTTAGCGCGTTGATCCAct                                  >  1:3072726/1‑68 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCTGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:980984/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:987577/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:685396/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:329421/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:3168147/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:3187465/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:1251705/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:1420349/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:2252810/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:2012346/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:1650989/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:1657647/1‑69 (MQ=255)
                                    gttgGCGACGCCGTGAAGATCATTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAggg                         >  1:2802266/1‑69 (MQ=255)
                                                        ccTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGcc       >  1:1420824/1‑67 (MQ=255)
                                                        ccTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCg      >  1:2557186/1‑68 (MQ=255)
                                                        ccTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCg      >  1:3413987/1‑68 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:1062446/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:1624125/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:1881242/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:2099761/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:2160103/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:239615/69‑1 (MQ=255)
                                                         cTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:2582700/69‑1 (MQ=255)
                                                          ttCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCt    <  1:2281487/68‑1 (MQ=255)
                                                           tCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCTGa  <  1:1948745/69‑1 (MQ=255)
                                                                 |                                                              
AATGCGCCGATCTTCCAGTACGCGGATTACGGCATTGTTGGCGACGCCGTGAAGATCCTTCCGGCGCTGACCGCAGCTTTAGCGCGTTGATCCACTCTGGCAGGGCTGCATTTTGGCCCTGCCGCTGA  >  minE/41718‑41845

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: