Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,986,542 A→G 19.9% L19L (TTG→CTG)  arcA ← DNA‑binding response regulator in two‑component regulatory system with ArcB or CpxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,986,5420AG19.9% 113.1 / 23.7 54L19L (TTG→CTG) arcADNA‑binding response regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base A (20/23);  new base G (3/7);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTA  >  minE/2986473‑2986608
                                                                     |                                                                  
tCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAAc                                                                   >  1:1225905/1‑71 (MQ=255)
     tGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGtt                                                              <  1:1186229/71‑1 (MQ=255)
       cATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTg                                                             <  1:1862607/70‑1 (MQ=255)
       cATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGc                                                            >  1:1682736/1‑71 (MQ=255)
         tttCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCgt                                                          <  1:1810332/71‑1 (MQ=255)
           tCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCgtg                                                         <  1:1560071/70‑1 (MQ=255)
              gcgcCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTAc                                                     <  1:768294/71‑1 (MQ=255)
              gcgcCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTAc                                                     <  1:209367/71‑1 (MQ=255)
              gcgcCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTAc                                                     <  1:2093047/71‑1 (MQ=255)
               cgcCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTAcc                                                    <  1:1286350/71‑1 (MQ=255)
                    tCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACt                                                >  1:1913679/1‑70 (MQ=255)
                    tCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTc                                               >  1:425315/1‑71 (MQ=255)
                      tGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGt                                             >  1:780400/1‑71 (MQ=255)
                      tGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGt                                             >  1:546437/1‑71 (MQ=255)
                        tCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCGTGTTACCaa                                                  >  1:1731080/1‑64 (MQ=255)
                        tCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACgt                                                                 <  1:99055/49‑1 (MQ=255)
                         cGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCg                                              <  1:886853/67‑1 (MQ=255)
                          gCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACgtgt                                                               >  1:920760/1‑49 (MQ=255)
                          gCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTc                                         <  1:1782035/71‑1 (MQ=255)
                            ttCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTc                                         <  1:225259/69‑1 (MQ=255)
                              cGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAAc                                      >  1:1459172/1‑70 (MQ=255)
                              cGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAAc                                      <  1:501526/70‑1 (MQ=255)
                               gAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGa                                    <  1:963523/71‑1 (MQ=255)
                                 aaaCATCATAGCCTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCa                                        >  1:1143601/1‑65 (MQ=255)
                                    catcatAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGa                               <  1:480807/71‑1 (MQ=255)
                                      tcatAGCCTTCCGCTTCGAAAATACTTTTCAACgtgt                                                               >  1:2028752/1‑37 (MQ=255)
                                      tcatAGCCTTCCGCTTCGAAAATACTTTTCAACgt                                                                 <  1:2124758/35‑1 (MQ=255)
                                       catAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCg                                              <  1:1696955/53‑1 (MQ=255)
                                        atAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCg                                              <  1:1882693/52‑1 (MQ=255)
                                         tAGCCTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGaa                              <  1:984659/67‑1 (MQ=255)
                                           gCCTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCg                        <  1:1832618/71‑1 (MQ=255)
                                           gCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCTTCTTCAACGATAAGAATGTGCg                        >  1:1439767/1‑71 (MQ=255)
                                             cTTCCGCTTCGAAAATACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCggg                      <  1:1010438/71‑1 (MQ=255)
                                              ttCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAAtgtg                          <  1:2131307/66‑1 (MQ=255)
                                               tCCGCTTCGAAAATACATTTCAACGTGTTGAGTCTTACCAACTCGTCTTCAACGATa                                  >  1:1714275/1‑57 (MQ=255)
                                                ccGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGt                    >  1:1385889/1‑70 (MQ=255)
                                                ccGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGt                    >  1:2104864/1‑70 (MQ=255)
                                                 cGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCt                  <  1:2137459/71‑1 (MQ=255)
                                                     tCGAAAATACTTTTCATCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCAt              >  1:2206205/1‑71 (MQ=255)
                                                        aaaaTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGtt           <  1:1072759/71‑1 (MQ=255)
                                                          aaTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTg         >  1:877303/1‑71 (MQ=255)
                                                          aaTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTg         >  1:752308/1‑71 (MQ=255)
                                                           aTACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGc        <  1:965826/71‑1 (MQ=255)
                                                           aTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCtgg                      >  1:1438436/1‑57 (MQ=255)
                                                           aTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCggg                      >  1:711192/1‑57 (MQ=255)
                                                           aTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCggg                      >  1:397169/1‑57 (MQ=255)
                                                           aTACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGttt          >  1:1297955/1‑69 (MQ=255)
                                                            tACTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCt       <  1:93667/71‑1 (MQ=255)
                                                            tACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCt       <  1:1668564/71‑1 (MQ=255)
                                                            tACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCt       <  1:1673173/71‑1 (MQ=255)
                                                             aCTTTTCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGAt                                   >  1:129631/1‑42 (MQ=255)
                                                               ttttCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGt            <  1:983573/63‑1 (MQ=255)
                                                                 ttCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGc        >  1:2114122/1‑65 (MQ=255)
                                                                  tCAGCGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTa  <  1:19054/70‑1 (MQ=255)
                                                                     |                                                                  
TCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTA  >  minE/2986473‑2986608

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: