Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,832 T→A 57.7% E316D (GAA→GAT eutG ← predicted alcohol dehydrogenase in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,614,8320TA57.7% 25.3 / 54.8 59E316D (GAA→GATeutGpredicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base A (13/21);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACA  >  minE/1614766‑1614878
                                                                  |                                              
cgcgTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGtt                                               <  1:3216040/68‑1 (MQ=255)
cgcgTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGtt                                               <  1:2402344/68‑1 (MQ=255)
cgcgTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGtt                                               <  1:1957400/68‑1 (MQ=255)
cgcgTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGtt                                               <  1:1573386/68‑1 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:2630033/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:996341/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:832287/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:788178/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:3284152/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:307836/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:3062612/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:2499816/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:2291892/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:2130547/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:1477823/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:1264599/1‑69 (MQ=255)
              cACGATCGTCGGATTTTTTAGTTCGCAGCGCCCGACCAATCTGACTAAACCGATCACGACAAACCATcc                                >  1:1549951/1‑69 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTa                           >  1:1070808/1‑46 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATGACCGTTGGCAGCaa    <  1:1569114/69‑1 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    >  1:928152/1‑69 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    >  1:2739105/1‑69 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    <  1:2623201/69‑1 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    >  1:2828117/1‑69 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    <  1:311930/69‑1 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    >  1:3278712/1‑69 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    <  1:3337290/69‑1 (MQ=255)
                                          tGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    >  1:1614215/1‑69 (MQ=255)
                                          tGCCCGACAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAAc   >  1:2997834/1‑69 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:1439596/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:305878/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:2840055/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:480666/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:1205614/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:2747598/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:1607465/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:2214978/65‑1 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACa  >  1:3073772/1‑69 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACa  >  1:555080/1‑69 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACa  >  1:397220/1‑69 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACa  >  1:1556207/1‑69 (MQ=255)
                                            cccGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACa  >  1:3306212/1‑69 (MQ=255)
                                             cccaCCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGgcagc      <  1:1649632/61‑1 (MQ=255)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:3286119/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1109031/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1586328/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:917660/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1647013/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1210828/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:524678/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:2006790/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:3299799/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:3292330/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:266981/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:2166080/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:128789/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:31785/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:2175269/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:3116640/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1339114/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:1517011/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:2441626/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTcc                     <  1:2493490/47‑1 (MQ=39)
                                               gACCAATCTGACTAAACCGATCACGACAAACCATCCGGTTAAATTca                     <  1:2182437/47‑2 (MQ=38)
                                                                  ttCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCaa    <  1:3333050/45‑1 (MQ=255)
                                                                  |                                              
CGCGTTAATAGCGTCACGATCGTCGGATTTTTTAGTTCGCAGTGCCCGACCAATCTGACTAAAGCGTTCACGACAAACCATCCGGTTAAATTCCATCACCGTTGGCAGCAACA  >  minE/1614766‑1614878

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: