Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,630,243 G→A 47.9% R517R (CGC→CGT aegA ← fused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,630,2430GA47.9% 4.7 / 63.0 75R517R (CGC→CGTaegAfused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (29/10);  new base A (17/19);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGTGA  >  minE/1630182‑1630310
                                                             |                                                                   
cTTGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCg                                                               <  1:339798/68‑1 (MQ=255)
  tGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttc                                                            >  1:264890/1‑69 (MQ=255)
  tGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttc                                                            >  1:2981068/1‑69 (MQ=255)
  tGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttc                                                            >  1:2139373/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:863070/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:2558502/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:2904188/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:2269477/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:401648/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:519256/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:1814744/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:670821/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:727909/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:1742186/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:850422/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:1657341/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:829070/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          >  1:1580160/1‑68 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCACCCTCTTCACGGGCGttcttc                                                         >  1:380285/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGACCCCTCTTCACGGGCGttcttc                                                         >  1:2387837/1‑69 (MQ=255)
     tcattcaGTTCAAGCGCCACCCGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttcttc                                                         >  1:1817996/1‑69 (MQ=255)
     tcattcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttctt                                                          <  1:400841/68‑1 (MQ=255)
      cattcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttcttc                                                         <  1:1662201/68‑1 (MQ=255)
      cattcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttcttc                                                         <  1:2098476/68‑1 (MQ=255)
      cattcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGttcttc                                                         <  1:3113134/68‑1 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:923832/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:2766331/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:849168/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:3099444/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:3127841/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:380642/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:801567/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:2455353/1‑69 (MQ=255)
         tcaGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCActt                                                     >  1:1688147/1‑69 (MQ=255)
                        ccGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTg                                      <  1:2689093/69‑1 (MQ=255)
                        ccGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTg                                      <  1:1982759/69‑1 (MQ=255)
                        ccGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTg                                      <  1:571234/69‑1 (MQ=255)
                        ccGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTg                                      <  1:3097179/69‑1 (MQ=255)
                                 cGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGtt                               <  1:2247366/67‑1 (MQ=255)
                                 cGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAg                             >  1:3040050/1‑69 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:734539/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:6030/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:1135685/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:1650084/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:855303/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:116121/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:1843713/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:868131/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2002909/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2028550/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:3433985/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2067932/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:3306472/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:3263848/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:3261043/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2944226/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2191219/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGGCGttctt                                                          <  1:2252180/38‑1 (MQ=255)
                                   ttAAATTCGAAGTTGGCCCCCTCTTCACGGCGttctt                                                          <  1:1909060/37‑1 (MQ=37)
                                          cGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCac                    >  1:2924537/1‑69 (MQ=255)
                                          cGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCac                    >  1:2800903/1‑69 (MQ=255)
                                          cGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCac                    >  1:1542198/1‑69 (MQ=255)
                                                 ggCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCa                                   >  1:3383477/1‑47 (MQ=255)
                                                 ggCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCa                                   >  1:23057/1‑47 (MQ=255)
                                                 ggCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCa                                   >  1:420803/1‑47 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAgg     >  1:821645/1‑68 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:2849108/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:2323554/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:2499524/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:1677870/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:1733013/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:2134350/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:3166818/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGt    >  1:3217749/1‑69 (MQ=255)
                                                          ttCGCGGGCGTTCTTCACTTCTTTCTTCGAGACTGGCATGTTAGCTtcatca                     >  1:3265835/1‑52 (MQ=255)
                                                            cgcgGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGTGa  >  1:1203943/1‑69 (MQ=255)
                                                             |                                                                   
CTTGTTCATTCAGCTCAAGCGCCACCGGCTGGACGTTAAATTCGAAGTTGGCCCCCTCTTCGCGGGCGTTCTTCACTTCTTTCTTCGAGCCTGGCATGTTAGCTTCATCACGACGATAAGCGCAGGTGA  >  minE/1630182‑1630310

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: