Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 828,622 T→C 7.2% S226G (AGT→GGT)  narX ← sensory histidine kinase in two‑component regulatory system with NarL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE828,6220TC 7.2% 247.8 / 8.8 83S226G (AGT→GGT) narXsensory histidine kinase in two‑component regulatory system with NarL
Reads supporting (aligned to +/- strand):  ref base T (44/33);  new base C (2/4);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GAGAGGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATG  >  minE/828552‑828692
                                                                      |                                                                      
gagagGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACt                                                                        >  1:27178/1‑71 (MQ=255)
gagagGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACt                                                                        <  1:2817854/71‑1 (MQ=255)
gagagGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACt                                                                        >  1:3978622/1‑71 (MQ=255)
gagagGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACt                                                                        >  1:4507569/1‑71 (MQ=255)
 agagGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACtt                                                                       <  1:56889/71‑1 (MQ=255)
   agGATCTGATTTTTATGCTCCAGCCCGGCAGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACtt                                                                       <  1:3755905/69‑1 (MQ=255)
       tCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGc                                                                  >  1:4419915/1‑70 (MQ=255)
           aTTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGtt                                                             >  1:1270154/1‑71 (MQ=255)
           aTTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGtt                                                             >  1:4115679/1‑71 (MQ=255)
                    cTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACa                                                    >  1:1813071/1‑71 (MQ=255)
                     tCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACa                                                    >  1:824387/1‑70 (MQ=255)
                       cAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACtt                                                                       <  1:1938362/49‑1 (MQ=255)
                       cAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACatat                                                 >  1:2396701/1‑71 (MQ=255)
                          cccGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGt                                              <  1:3209979/71‑1 (MQ=255)
                           ccGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGtt                                             <  1:3009660/71‑1 (MQ=255)
                           ccGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGtt                                             <  1:3521455/71‑1 (MQ=255)
                               cggTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAg                                                       >  1:3262539/1‑57 (MQ=255)
                               cggTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAAc                                         <  1:3342373/71‑1 (MQ=255)
                                  ttttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCgg                                                                   <  1:235486/42‑1 (MQ=255)
                                  ttttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCgg                                                                   <  1:736516/42‑1 (MQ=255)
                                  ttttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGc                                                                  <  1:1692260/43‑1 (MQ=255)
                                  ttttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGc                                                                  <  1:416940/43‑1 (MQ=255)
                                  ttttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACg                                        >  1:15457/1‑69 (MQ=255)
                                   tttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCa                                      >  1:941332/1‑70 (MQ=255)
                                    ttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCa                                                                >  1:964960/1‑43 (MQ=255)
                                    ttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCa                                                                >  1:29379/1‑43 (MQ=255)
                                    ttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGt                                    >  1:1960156/1‑71 (MQ=255)
                                    ttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGt                                    >  1:123981/1‑71 (MQ=255)
                                    ttCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGt                                    >  1:2596554/1‑71 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAg                                                       <  1:1618684/47‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAg                                                       <  1:1472128/47‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAg                                                       <  1:1308715/47‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTc                                  <  1:1492957/68‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCaa                               <  1:4591372/71‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTATGCAg                                                       <  1:1142153/47‑1 (MQ=255)
                                         cTGAACCCGCTGCTCAAGTACGGAATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTc                                  <  1:4323864/68‑1 (MQ=255)
                                            aaCCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGc                             <  1:429943/70‑1 (MQ=255)
                                            aaCCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGc                             <  1:3801083/70‑1 (MQ=255)
                                            aaCCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGc                             <  1:2044496/70‑1 (MQ=255)
                                             aCCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCAt                           >  1:4512603/1‑71 (MQ=255)
                                             aCCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCAt                           >  1:2924335/1‑71 (MQ=255)
                                              cccGCTGCTCAAGTACTGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCTGTTc                                  <  1:3800756/63‑1 (MQ=255)
                                                  ctgctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTc                                            <  1:4736915/49‑1 (MQ=255)
                                                  ctgctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTc                                            <  1:479315/49‑1 (MQ=255)
                                                  ctgctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTc                                            <  1:149728/49‑1 (MQ=255)
                                                  ctgctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTc                                            <  1:3293539/49‑1 (MQ=255)
                                                    gctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTCaa                                          >  1:2129665/1‑49 (MQ=255)
                                                    gctCAAGTACGGCATAACCTTCGGCCAGTTCTGCAGACATATTGTTCaa                                          >  1:3723831/1‑49 (MQ=255)
                                                       cAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCAttt                   >  1:895466/1‑69 (MQ=255)
                                                        aaGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:4617135/69‑1 (MQ=255)
                                                        aaGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGt                <  1:3727376/71‑1 (MQ=255)
                                                           tACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:795715/1‑71 (MQ=255)
                                                           tACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:1773038/1‑71 (MQ=255)
                                                           tACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:4685525/1‑71 (MQ=255)
                                                           tACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:2343693/1‑71 (MQ=255)
                                                           tACGGCATAACTTTCGGCCAGTTCGGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:4245713/1‑71 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAg                                     >  1:3148000/1‑45 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAg                                     >  1:2806507/1‑45 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAg                                     >  1:3036477/1‑45 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAg                                     >  1:4312789/1‑45 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCAttt                   >  1:2391817/1‑63 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCAttt                   >  1:1911661/1‑63 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCAtta                   >  1:3683531/1‑62 (MQ=255)
                                                             cGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgcgc           >  1:779048/1‑71 (MQ=255)
                                                              ggCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTcc                                 >  1:3303711/1‑48 (MQ=255)
                                                              ggCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTcc                                 >  1:1637680/1‑48 (MQ=255)
                                                              ggCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTcc                                 >  1:1895110/1‑48 (MQ=255)
                                                              ggCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTcc                                 >  1:730773/1‑48 (MQ=255)
                                                              ggCATAACTTTCGGCCAGATCTGCAGACATATTGTTCAACGCAGTTcc                                 >  1:2082499/1‑48 (MQ=255)
                                                                cATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGc                                       <  1:1839719/40‑1 (MQ=255)
                                                                cATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGc                                       <  1:2953598/40‑1 (MQ=255)
                                                                cATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGc                                       <  1:2968357/40‑1 (MQ=255)
                                                                cATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGc                                       <  1:3530096/40‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCa                            <  1:2973315/49‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCa                            <  1:3509164/49‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCa                            <  1:2044824/49‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:553165/59‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:1126989/59‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:1823974/59‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:3766116/59‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTc                  <  1:2102164/59‑1 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCt      >  1:4497773/1‑71 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCt      >  1:50236/1‑71 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCt      >  1:1338051/1‑71 (MQ=255)
                                                                  tAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCt      >  1:2576984/1‑71 (MQ=255)
                                                                    aCTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGt                >  1:1247356/1‑59 (MQ=255)
                                                                    aCTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGa    >  1:2673159/1‑71 (MQ=255)
                                                                     cTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTgc             >  1:1439952/1‑61 (MQ=255)
                                                                     cTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGAt   >  1:4316319/1‑71 (MQ=255)
                                                                     cTTTCGGCAAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGACATTTCGTTGCGccc         >  1:512184/1‑65 (MQ=255)
                                                                      tttCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATg  >  1:4718846/1‑71 (MQ=255)
                                                                      tttCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATg  >  1:2013857/1‑71 (MQ=255)
                                                                      tttCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATg  >  1:3744977/1‑71 (MQ=255)
                                                                      tttCGGCCAGTTATGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATg  >  1:2425234/1‑71 (MQ=255)
                                                                      |                                                                      
GAGAGGATCTGATTTTTATGCTCCAGCCCGGCGGTTTTCTCCTGAACCCGCTGCTCAAGTACGGCATAACTTTCGGCCAGTTCTGCAGACATATTGTTCAACGCAGTTCCAAGCATCGCCATTTCGTTGCGCCCGCTGATG  >  minE/828552‑828692

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: