Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,761,918 G→A 62.3% Q244Q (CAG→CAA aroE → dehydroshikimate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,761,9180GA62.3% 24.7 / 32.4 53Q244Q (CAG→CAAaroEdehydroshikimate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (5/15);  new base A (14/19);  total (19/34)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAGCAGCGAGGCTCAAAGCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCAATTGC  >  minE/2761861‑2761985
                                                         |                                                                   
gtgtGAGCAGCGAGGCTCAAAGCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGc                                                          <  1:1671792/69‑1 (MQ=255)
gtgtGAGCAGCGAGGCTCAAAGCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGc                                                          <  1:3398343/69‑1 (MQ=255)
gtgtGAGCAGCGAGGCTCAAAGCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGc                                                          <  1:1899841/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:3024171/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1843772/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1886714/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1909281/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1991476/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:2202499/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:2971554/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:2400011/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:3067986/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:3072039/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:3230991/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:3427529/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:37932/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:882152/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1753441/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1785472/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1585827/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1129685/69‑1 (MQ=255)
  gtgAGCAGCGAGGCTCAAAACGTACTGCTGATGGTTTAGGAATGCTGGTGGCACAAGCGGCTCATGCCt                                                        <  1:1038919/69‑1 (MQ=255)
                     gCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:1511426/69‑1 (MQ=255)
                     gCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:872131/69‑1 (MQ=255)
                     gCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:3340929/69‑1 (MQ=255)
                     gCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:1062070/69‑1 (MQ=255)
                      cGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:1289029/68‑1 (MQ=255)
                      cGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGgtgt                                     <  1:2820781/68‑1 (MQ=255)
                          aTGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTttct                                                    <  1:2947887/49‑1 (MQ=255)
                          aTGCTGATGGTTTAGGAACGCTGGTGGCACAGGCGGCTCATGCCTttct                                                    <  1:243358/49‑1 (MQ=255)
                                              cTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTta              >  1:1144376/1‑67 (MQ=255)
                                              cTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTta              >  1:3126813/1‑67 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:399911/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:427636/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:608718/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:67670/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:822295/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:1015705/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:1177823/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:3062305/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:2989751/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:1011023/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:1292368/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:2178295/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:2082810/1‑68 (MQ=255)
                                                 gtggCACAAGCGGCTCATGCCTTCCTTCTTTGGCATGGCGTTCTGCCTGATGTAGAACCGGTCATcaa          >  1:1632540/1‑68 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       >  1:1134705/1‑69 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       <  1:1283957/69‑1 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       <  1:271974/69‑1 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       >  1:578492/1‑69 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       >  1:2191845/1‑69 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       <  1:1429927/69‑1 (MQ=255)
                                                   ggCACAGGCGGCTCATGCCTTTCTTCTATGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCa       <  1:2806834/69‑1 (MQ=255)
                                                         ggcggcTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCAATTGc  >  1:360531/1‑68 (MQ=255)
                                                         |                                                                   
GTGTGAGCAGCGAGGCTCAAAGCGTAATGCTGATGGTTTAGGAATGCTGGTGGCACAGGCGGCTCATGCCTTTCTTCTCTGGCACGGTGTTCTGCCTGACGTAGAACCAGTTATAAAGCAATTGC  >  minE/2761861‑2761985

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: