Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 246,374 C→T 17.2% S6L (TCG→TTG)  phoR → sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE246,3740CT17.2% 132.9 / 19.2 58S6L (TCG→TTG) phoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base T (2/8);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTGCATTTTT  >  minE/246304‑246440
                                                                      |                                                                  
cTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTc                                                                    <  1:2885222/71‑1 (MQ=255)
cTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTc                                                                    <  1:866032/71‑1 (MQ=255)
cTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTc                                                                    <  1:542488/71‑1 (MQ=255)
cTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTc                                                                    <  1:1634387/71‑1 (MQ=255)
 ttGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCg                                                                   >  1:57444/1‑71 (MQ=255)
  tGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGt                                                                  >  1:3950549/1‑71 (MQ=255)
   gCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTg                                                                 >  1:33705/1‑71 (MQ=255)
     tCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTg                                                                 <  1:1984397/69‑1 (MQ=255)
       aTCGGACGCAGAGCAGGGCTTATGATTTATTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGaa                                                              >  1:3789726/1‑70 (MQ=255)
          ggACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAgg                                                          >  1:3888839/1‑71 (MQ=255)
             cGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTg                                                       <  1:3002050/71‑1 (MQ=255)
             cGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTg                                                       <  1:3155606/71‑1 (MQ=255)
              gCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTgg                                                      <  1:2091357/71‑1 (MQ=255)
               cAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGt                                                     <  1:1375235/71‑1 (MQ=255)
                agagCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTc                                                                    >  1:3358678/1‑55 (MQ=255)
                agagCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTg                                                                 >  1:2808771/1‑58 (MQ=255)
                agagCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGaaaa                                                            >  1:3301495/1‑63 (MQ=255)
                  agCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGc                                                   <  1:694464/70‑1 (MQ=255)
                      gggCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAg                                              <  1:1709634/71‑1 (MQ=255)
                      gggCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAg                                              <  1:1574182/71‑1 (MQ=255)
                         cTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAg                                                           <  1:536742/55‑1 (MQ=255)
                          ttATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAgctgc                                          >  1:1602932/1‑71 (MQ=255)
                            aTGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAgctgc                                          <  1:1595265/69‑1 (MQ=255)
                              gATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAg                                              >  1:463861/1‑63 (MQ=255)
                              gATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACt                                      >  1:2310713/1‑71 (MQ=255)
                                tttCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTTGTGGaaaa                                                            >  1:1612679/1‑47 (MQ=255)
                                tttCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTTGTGGaaaa                                                            >  1:1036979/1‑47 (MQ=255)
                                tttCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTgg                                                >  1:1139668/1‑59 (MQ=255)
                                  tCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACttt                                    <  1:259942/69‑1 (MQ=255)
                                   cTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTtgct                                 >  1:1176390/1‑71 (MQ=255)
                                   cTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTtgc                                  <  1:234671/70‑1 (MQ=255)
                                    ttAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGaa                                                              >  1:992520/1‑41 (MQ=255)
                                       aCTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGcc                              <  1:1393504/70‑1 (MQ=255)
                                        cTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTc                            <  1:1836602/71‑1 (MQ=255)
                                          ggAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTccc                          >  1:3866456/1‑71 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTgg                                                <  1:981442/45‑1 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTgg                                                <  1:3323964/45‑1 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTGGAGCTGCTACtt                                     <  1:3281677/56‑1 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACttt                                    <  1:2149313/57‑1 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACtt                                     <  1:1470141/56‑1 (MQ=255)
                                              tATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTtgctgc                               >  1:3456202/1‑62 (MQ=255)
                                                tCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCttt                    >  1:2155168/1‑71 (MQ=255)
                                                tCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCttt                    >  1:3069421/1‑71 (MQ=255)
                                                tCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCttt                    >  1:1833958/1‑71 (MQ=255)
                                                  ttACGTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCt                                                  <  1:1604272/39‑1 (MQ=255)
                                                  ttACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAgc                                             <  1:761692/44‑1 (MQ=255)
                                                    aCGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCttt                    <  1:116762/67‑1 (MQ=255)
                                                    aCGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCAt                 >  1:2566389/1‑70 (MQ=255)
                                                     cGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGTTGGAGTTGCTACTTTGCTGCCTCCCGGCTTTCATc                <  1:3754317/70‑1 (MQ=255)
                                                     cGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAgc                                             <  1:622306/41‑1 (MQ=255)
                                                     cGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATc                <  1:352945/70‑1 (MQ=255)
                                                      gTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCt              <  1:3741745/71‑1 (MQ=255)
                                                      gTGCTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCt              <  1:2254351/71‑1 (MQ=255)
                                                      gTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCt              <  1:2078662/71‑1 (MQ=255)
                                                       tGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCgg                        <  1:3490766/60‑1 (MQ=255)
                                                         cTGGAACGGCTGTTGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCa                  <  1:1568894/64‑1 (MQ=255)
                                                         cTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTgg            <  1:1573568/70‑1 (MQ=255)
                                                          tGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTgg            <  1:64517/69‑1 (MQ=255)
                                                            gAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTg         >  1:648318/1‑70 (MQ=255)
                                                              aCGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTGCAt      >  1:2454284/1‑71 (MQ=255)
                                                                ggCTGTTGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTGCAttt    <  1:2686792/71‑1 (MQ=255)
                                                                  cTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTtgc                                  >  1:981290/1‑39 (MQ=255)
                                                                  cTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTGCAttttt  <  1:3748502/71‑1 (MQ=255)
                                                                      |                                                                  
CTTGCTCATCGGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCCTCCCGGCTTTCATCCTGGGTGCATTTTT  >  minE/246304‑246440

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: