Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,417,804 C→T 5.3% intergenic (+85/‑248) 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,8040CT 5.3% 243.4 / 3.2 76intergenic (+85/‑248)proL/narPtRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base C (48/23);  new base T (2/2);  total (51/25)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTC  >  minE/1417736‑1417872
                                                                    |                                                                    
accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                    >  1:3546563/1‑71 (MQ=255)
accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                    >  1:1799754/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:2298954/70‑1 (MQ=255)
          aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                                           <  1:630439/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:558509/1‑71 (MQ=255)
                                       ttCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTa                             >  1:3249589/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:3353416/71‑1 (MQ=255)
                                                       tGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGAtttttt             <  1:1947184/71‑1 (MQ=255)
                                                         acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:4052812/71‑1 (MQ=255)
                                                         acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:3951594/71‑1 (MQ=255)
                                                         acacATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg           <  1:91278/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:829509/70‑1 (MQ=255)
                                                              tCACGGTCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttgt       <  1:729088/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:3494038/1‑36 (MQ=255)
                                                                acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:3001392/1‑36 (MQ=255)
                                                                acaGCCTGTTTATTTTCTGTTAACTGTTTAAAGaaa                                       >  1:99650/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:3334963/71‑1 (MQ=255)
                                                                  aGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTc  >  1:3229344/1‑71 (MQ=255)
                                                                    |                                                                    
ACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTC  >  minE/1417736‑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: