Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,453 G→C 6.9% intergenic (‑78/+139) 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,4530GC 6.9% 183.1 / 3.2 58intergenic (‑78/+139)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAA  >  minE/1614384‑1614523
                                                                     |                                                                      
cTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGt                                                                       <  1:590887/71‑1 (MQ=255)
   ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTa                                                                      <  1:1234003/69‑1 (MQ=255)
   ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                    >  1:585577/1‑71 (MQ=255)
   ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                    >  1:1964364/1‑71 (MQ=255)
         aGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGa                                                               >  1:3048216/1‑70 (MQ=255)
          gttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                                                              <  1:2624477/70‑1 (MQ=255)
          gttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                                                              <  1:49593/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)
               tatttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGc                                                                   >  1:3003371/1‑60 (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:1054854/71‑1 (MQ=255)
                                         gCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGc                              <  1:2316897/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:3246438/69‑1 (MQ=255)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                 <  1:1599352/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:1820989/1‑52 (MQ=255)
                                                         aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                 >  1:3735406/1‑52 (MQ=255)
                                                         aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                 >  1:1215143/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:1277451/71‑1 (MQ=255)
                                                                   ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg    <  1:3237139/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)
                                                                     |                                                                      
CTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAA  >  minE/1614384‑1614523

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: