Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,417,802 A→G 6.0% intergenic (+83/‑250) proL → / → narP tRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,417,8020AG 6.0% 181.5 / 4.0 67intergenic (+83/‑250)proL/narPtRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base A (41/22);  new base G (2/2);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

AAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTC  >  minE/1417732‑1417872
                                                                      |                                                                      
aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                                                        <  1:1860026/71‑1 (MQ=255)
aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                                                        <  1:2053566/71‑1 (MQ=255)
    accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                    >  1:1799754/1‑71 (MQ=255)
    accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                    >  1:3546563/1‑71 (MQ=255)
        accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                                                 >  1:3396603/1‑70 (MQ=255)
         cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                                               <  1:2069672/71‑1 (MQ=255)
          ccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                                              <  1:1948386/71‑1 (MQ=255)
            ttACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                            >  1:1231866/1‑71 (MQ=255)
              aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                    >  1:2684234/1‑61 (MQ=255)
              aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                                           <  1:630439/70‑1 (MQ=255)
              aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                                           <  1:2298954/70‑1 (MQ=255)
              aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                            <  1:3051711/69‑1 (MQ=255)
              aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                          >  1:3109944/1‑71 (MQ=255)
                 tttggttTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                                             <  1:2331348/64‑1 (MQ=255)
                      tttttttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACt                                                  <  1:3233242/71‑1 (MQ=255)
                       ttttttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                            <  1:3981545/60‑1 (MQ=255)
                        tttttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGt                                                >  1:3485906/1‑71 (MQ=255)
                        tttttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGt                                                >  1:87407/1‑71 (MQ=255)
                        tttttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGt                                                >  1:361242/1‑71 (MQ=255)
                          tttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGtt                                               >  1:1281458/1‑70 (MQ=255)
                           ttATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGttt                                              <  1:1886374/70‑1 (MQ=255)
                             atatCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAAGTGTTTa                                             >  1:1198623/1‑69 (MQ=255)
                                  tGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                          >  1:3653086/1‑51 (MQ=255)
                                   gCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTaa                                                    >  1:1655322/1‑56 (MQ=255)
                                   gCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTaa                                                    >  1:2037346/1‑56 (MQ=255)
                                    cAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATg                                     >  1:1122201/1‑70 (MQ=255)
                                    cAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTACTGTTTAAAGAAATGc                                    <  1:3527874/70‑1 (MQ=255)
                                      attaattCGATAAACAGACCGTGACACATCACAGc                                                                      >  1:953804/1‑35 (MQ=255)
                                      attaattCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:2987508/1‑66 (MQ=255)
                                      attaattCGAAAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:127601/1‑66 (MQ=255)
                                        taattCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGcac                                  >  1:372331/1‑69 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGt                                                                  >  1:3596767/1‑36 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                                             >  1:1893129/1‑41 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                                             >  1:2514878/1‑41 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                                             >  1:4029323/1‑41 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                                              >  1:1682024/1‑40 (MQ=255)
                                         aattCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAg                               >  1:1100724/1‑71 (MQ=255)
                                          attCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAGTGCACAAGt                              <  1:1129950/71‑1 (MQ=255)
                                          attCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAg                               <  1:4095664/70‑1 (MQ=255)
                                           ttCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTa                             >  1:3249589/1‑71 (MQ=255)
                                           ttCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTa                             >  1:558509/1‑71 (MQ=255)
                                           ttCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTa                             >  1:3207228/1‑71 (MQ=255)
                                             cGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGt                                                >  1:1077456/1‑50 (MQ=255)
                                             cGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGt                                                >  1:135433/1‑50 (MQ=255)
                                               aTAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                                                      >  1:3909296/1‑42 (MQ=255)
                                               aTAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                                                      >  1:888810/1‑42 (MQ=255)
                                                 aaaCAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGtga                       >  1:887728/1‑71 (MQ=255)
                                                  aaCAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGtgat                      >  1:568125/1‑71 (MQ=255)
                                                    cAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGATATGcac                                  >  1:1976010/1‑57 (MQ=255)
                                                    cAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGtgatt                     >  1:3261147/1‑70 (MQ=255)
                                                       aCCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTgattgatt                 >  1:729718/1‑71 (MQ=255)
                                                          gTGACACATCACAGCCTGTTTATTTTCTGTTAACt                                                  <  1:627010/35‑1 (MQ=255)
                                                           tGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGaa                                        >  1:1059539/1‑44 (MQ=255)
                                                           tGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGaa                                        >  1:1835804/1‑44 (MQ=255)
                                                           tGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGAtttttt             <  1:1947184/71‑1 (MQ=255)
                                                           tGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGAtttttt             <  1:3353416/71‑1 (MQ=255)
                                                             acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:3951594/71‑1 (MQ=255)
                                                             acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:91278/71‑1 (MQ=255)
                                                             acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:4052812/71‑1 (MQ=255)
                                                              cacaTCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGTAATGCACCAGTATTGTGATTGATTTTTTAgt          <  1:3104955/71‑1 (MQ=255)
                                                               acaTCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgt          >  1:1747749/1‑70 (MQ=255)
                                                                caTCACAGCCTGTTTATTTTCTGTTAACTGTTTaat                                           <  1:2393512/36‑2 (MQ=255)
                                                                 gtCACGGTCTGTTTATTTTCTGTTAACTGTTTaaa                                           >  1:974015/2‑35 (MQ=25)
                                                                 gtCACGGTCTGTTTATTTTCTGTTAACTGTTTaaa                                           >  1:2050147/2‑35 (MQ=25)
                                                                 aTCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttgt       <  1:1986963/71‑1 (MQ=255)
                                                                 aTCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           >  1:1560221/1‑67 (MQ=255)
                                                                  tCACGGTCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttgt       <  1:729088/70‑1 (MQ=255)
                                                                  tCACGGTCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttgt       <  1:829509/70‑1 (MQ=255)
                                                                  tCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATg                                     <  1:3322025/40‑1 (MQ=255)
                                                                   cacaGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACaa                                >  1:3165920/1‑44 (MQ=255)
                                                                    acaGGCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGtttt    >  1:2809644/1‑71 (MQ=255)
                                                                    acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:99650/1‑36 (MQ=255)
                                                                    acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:3494038/1‑36 (MQ=255)
                                                                    acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:3001392/1‑36 (MQ=255)
                                                                    acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGtttt    >  1:1388098/1‑71 (MQ=255)
                                                                     caGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGttttt   >  1:1416834/1‑71 (MQ=255)
                                                                     caGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGttttt   >  1:1106181/1‑71 (MQ=255)
                                                                      aGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTc  >  1:3229344/1‑71 (MQ=255)
                                                                      aGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTc  <  1:3334963/71‑1 (MQ=255)
                                                                      |                                                                      
AAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTC  >  minE/1417732‑1417872

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: