Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,886,680 G→A 16.9% L714L (TTG→TTA barA → hybrid sensory histidine kinase, in two‑component regulatory system with UvrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,886,6800GA16.9% 145.2 / 16.0 88L714L (TTG→TTAbarAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base A (6/9);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

CAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATCCACCA  >  minE/1886612‑1886742
                                                                    |                                                              
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGGCGTTCGATTTg                                                                <  1:2555889/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:2314745/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:1431177/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:3285816/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:2658462/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:2522285/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:2405764/69‑1 (MQ=255)
cAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTg                                                                <  1:2403392/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGTCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:1864716/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:175151/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:2773389/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:2788341/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:2900377/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:30440/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3076917/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3112790/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3093201/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3118622/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:1396276/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:1291759/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3296869/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:3346130/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:539714/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:770576/69‑1 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:773791/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:930591/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:960248/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGACGCCGTTCGATTTGATCtt                                                           >  1:3150298/1‑69 (MQ=255)
     tgtgGAACTTTGCGATAGCGGCCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:1232419/1‑69 (MQ=255)
     tgtgGAACTTTGAGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           >  1:1351006/1‑69 (MQ=255)
       tgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTaa                                                         >  1:995386/1‑69 (MQ=255)
       tgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTaa                                                         >  1:2619468/1‑69 (MQ=255)
       tgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTaa                                                         >  1:152305/1‑69 (MQ=255)
       tgGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTaa                                                         >  1:162218/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:1936589/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:1993709/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:1411513/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:2244155/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:2692700/1‑69 (MQ=255)
         gAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTAATCTTAATg                                                       >  1:173884/1‑69 (MQ=255)
           aCTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGa                                                     <  1:1018544/69‑1 (MQ=255)
           aCTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGa                                                     <  1:2932802/69‑1 (MQ=255)
           aCTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGa                                                     <  1:2930853/69‑1 (MQ=255)
                  gATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCa                                               <  1:1270932/68‑1 (MQ=255)
                  gATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCa                                               <  1:2968004/68‑1 (MQ=255)
                  gATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCa                                               <  1:2357471/68‑1 (MQ=255)
                  gATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCa                                               <  1:1864517/68‑1 (MQ=255)
                      gCGGGCATCAGGCGGTTGAACGGGCGACACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGc                                          <  1:2810099/69‑1 (MQ=255)
                      gCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGc                                          <  1:2422382/69‑1 (MQ=255)
                              cAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatg                                  <  1:3292393/69‑1 (MQ=255)
                              cAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatg                                  <  1:1318013/69‑1 (MQ=255)
                                ggCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatgga                                >  1:617664/1‑69 (MQ=255)
                                 gCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCtt                                                           <  1:2818632/41‑1 (MQ=255)
                                   ggTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatggatgg                             >  1:1776777/1‑69 (MQ=255)
                                   ggTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatggatgg                             <  1:113732/69‑1 (MQ=255)
                                   ggTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatggatgg                             <  1:3103191/69‑1 (MQ=255)
                                   ggTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatggatgg                             <  1:1129391/69‑1 (MQ=255)
                                   ggTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatggatgg                             >  1:1114110/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatg                                  <  1:1917632/63‑1 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatg                                  <  1:222256/63‑1 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACatg                                  <  1:3323870/63‑1 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:1469787/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:325098/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:796494/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:3305044/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:1030962/1‑69 (MQ=255)
                                    gTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGc                            >  1:562363/1‑69 (MQ=255)
                                     ttGAACGGGCGAAACAGCTGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCa                           <  1:608912/69‑1 (MQ=255)
                                     ttGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCa                           <  1:2800866/69‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:3435861/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:956722/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:789571/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:1956068/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:1573462/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:1870783/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:2582412/36‑1 (MQ=255)
                                      tGAACGGGCGAAACAGATGCCGTTCGATGTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCag                          >  1:2466575/1‑68 (MQ=255)
                                      tGAAAGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:1694066/36‑1 (MQ=38)
                                       gAACGGGCGAAACAGATGCCGTTCGATTTAATCtt                                                           <  1:2834596/35‑1 (MQ=255)
                                            ggCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    >  1:974037/1‑69 (MQ=255)
                                            ggCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    >  1:2537667/1‑69 (MQ=255)
                                            ggCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    >  1:1162431/1‑69 (MQ=255)
                                            ggCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    >  1:1634407/1‑69 (MQ=255)
                                            ggCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGCACATGGATGGCATTCGGGc                    >  1:1155583/1‑69 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:3289303/63‑1 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:974055/63‑1 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:1730464/63‑1 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:777658/63‑1 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:381097/63‑1 (MQ=255)
                                                  aCAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGc                    <  1:157662/63‑1 (MQ=255)
                                                       tGCCGTTCGATTTGATCTTAATGGATATTCAAATGc                                          <  1:1814919/36‑1 (MQ=255)
                                                              cGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATccacca  <  1:1133244/69‑1 (MQ=255)
                                                              cGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATccacca  <  1:286365/69‑1 (MQ=255)
                                                              cGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATccacca  >  1:2972806/1‑69 (MQ=255)
                                                              cGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATccacca  <  1:3045816/69‑1 (MQ=255)
                                                                 tttGATCTTAATGGATATTCAAATGCCGGACATGGATGGCa                           <  1:3066213/41‑1 (MQ=255)
                                                                    |                                                              
CAACATGTGGAACTTTGCGATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACATGGATGGCATTCGGGCCTGCGAACTCATCCACCA  >  minE/1886612‑1886742

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: