Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,459 T→C 10.7% intergenic (‑84/+133) eutH ← / ← eutG predicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,614,4590TC10.7% 141.6 / 9.0 58intergenic (‑84/+133)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base C (4/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614393‑1614528
                                                                  |                                                                     
aGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGa                                                                    >  1:3048216/1‑70 (MQ=255)
 gttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                                                                   <  1:49593/70‑1 (MQ=255)
 gttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                                                                   <  1:2624477/70‑1 (MQ=255)
  ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                                                                 <  1:1615146/71‑1 (MQ=255)
  ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCGGGATaa                                                                 <  1:2091896/71‑1 (MQ=255)
    gttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAg                                                                <  1:3029307/70‑1 (MQ=255)
    gttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGa                                                               >  1:3297469/1‑71 (MQ=255)
     ttatttatTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTg                                                                      <  1:3441438/63‑1 (MQ=255)
      tatttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTg                                                                      <  1:2279511/62‑1 (MQ=255)
          tattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCgg                                                          >  1:1788477/1‑70 (MQ=255)
              ggCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAg                                                      >  1:4088015/1‑70 (MQ=40)
               gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                    <  1:3648030/71‑1 (MQ=39)
                     gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCAt                                              >  1:2160004/1‑71 (MQ=14)
                     gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCAt                                              >  1:2548088/1‑71 (MQ=39)
                         cgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCgg                                          >  1:1725374/1‑71 (MQ=14)
                          gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                         <  1:3225980/71‑1 (MQ=39)
                          gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                         <  1:149226/71‑1 (MQ=39)
                          gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                         <  1:1481487/71‑1 (MQ=39)
                           tAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                         <  1:3780611/70‑1 (MQ=39)
                           tAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCa                                        <  1:2721790/71‑1 (MQ=38)
                             aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                       >  1:2087881/1‑70 (MQ=37)
                             aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                      >  1:3785453/1‑71 (MQ=38)
                             aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                      >  1:2293539/1‑71 (MQ=38)
                                gCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGc                                   <  1:2316897/71‑1 (MQ=255)
                                gCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGc                                   <  1:1054854/71‑1 (MQ=255)
                                  cTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAc                                 >  1:2086139/1‑71 (MQ=255)
                                    tATCCGCCCTACATGTGCAATCCCGTAGGCAGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACcg                               >  1:1002052/1‑71 (MQ=255)
                                     aTCCGCCCTACATGTGCAATCCCGTAGGCTGGATAATACGCGGCAAGCGTCGCATCCGGCAATTGCACcgc                              >  1:15835/1‑71 (MQ=255)
                                     aTCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAtt                                     >  1:3390398/1‑64 (MQ=40)
                                           ccTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCAc                         <  1:441987/70‑1 (MQ=17)
                                             tACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCa                                        <  1:935385/53‑1 (MQ=2)
                                             tACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                         <  1:194964/52‑1 (MQ=35)
                                              aCATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGt                                                   >  1:1057612/1‑41 (MQ=38)
                                              aCATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGt                                                   >  1:1608772/1‑41 (MQ=38)
                                              aCATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                  >  1:3041866/1‑58 (MQ=255)
                                               cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                      <  1:1599352/69‑1 (MQ=255)
                                               cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                      <  1:3246438/69‑1 (MQ=255)
                                               cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                    >  1:3317216/1‑71 (MQ=255)
                                               cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                    <  1:509452/71‑1 (MQ=255)
                                               cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                    >  1:1241331/1‑71 (MQ=255)
                                                aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGc                                                     <  1:920064/37‑1 (MQ=40)
                                                aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                      >  1:3735406/1‑52 (MQ=255)
                                                aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                      >  1:1215143/1‑52 (MQ=255)
                                                aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                      >  1:1820989/1‑52 (MQ=255)
                                                  gtgCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCa                                        <  1:1707396/48‑1 (MQ=38)
                                                  gtgCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                         <  1:3527178/63‑1 (MQ=255)
                                                   tgCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTc                                                  >  1:659091/1‑37 (MQ=40)
                                                       aTCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggc            >  1:746497/1‑71 (MQ=255)
                                                        tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg           >  1:2300886/1‑71 (MQ=255)
                                                         cccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGa                  >  1:622255/1‑63 (MQ=255)
                                                          ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg         <  1:3237139/71‑1 (MQ=255)
                                                          ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg         <  1:1277451/71‑1 (MQ=255)
                                                          ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg         <  1:1959502/71‑1 (MQ=255)
                                                            gTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgga       >  1:1908908/1‑69 (MQ=11)
                                                            gTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGaa       >  1:107921/1‑71 (MQ=255)
                                                               ggCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc    <  1:987206/71‑1 (MQ=255)
                                                               ggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc    >  1:1985257/1‑71 (MQ=14)
                                                                 cTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt  <  1:1344270/71‑1 (MQ=255)
                                                                  tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                    <  1:3679221/52‑1 (MQ=255)
                                                                  |                                                                     
AGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614393‑1614528

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: