Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,563,207 C→T 6.6% Q1072* (CAG→TAG)  evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,563,2070CT 6.6% 326.5 / 8.7 106Q1072* (CAG→TAG) evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base C (52/47);  new base T (3/4);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 108 total aligned reads displayed.

TTAAGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAAATAATA  >  minE/1563138‑1563276
                                                                     |                                                                     
ttAAGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCa                                                                      <  1:1809092/71‑1 (MQ=255)
ttAAGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCa                                                                      <  1:529166/71‑1 (MQ=255)
  aaGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAg                                                                     <  1:1788358/70‑1 (MQ=255)
  aaGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAg                                                                     <  1:3345294/70‑1 (MQ=255)
   aGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAg                                                                     <  1:3942425/69‑1 (MQ=255)
   aGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGtt                                                                   <  1:2127843/71‑1 (MQ=255)
   aGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGt                                                                    <  1:1316478/70‑1 (MQ=255)
   aGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGt                                                                    <  1:3949620/70‑1 (MQ=255)
   aGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGt                                                                    <  1:2616303/70‑1 (MQ=255)
    gTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTc                                                                       <  1:2597227/66‑1 (MQ=255)
    gTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTTAGt                                                                    <  1:1080082/69‑1 (MQ=255)
      tGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTaca                                                                >  1:2163671/1‑71 (MQ=255)
        cGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACAcc                                                              <  1:3332718/71‑1 (MQ=255)
         ggCATGAACTTATGTTAGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACAcc                                                              >  1:910190/1‑70 (MQ=255)
          gCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTaca                                                                >  1:1442249/1‑67 (MQ=255)
          gCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACAcc                                                              >  1:2182926/1‑69 (MQ=255)
           cATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAg                                                           >  1:2726695/1‑71 (MQ=255)
              gAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTg                                                        >  1:2135060/1‑71 (MQ=255)
               aaCTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTgc                                                       <  1:1795587/71‑1 (MQ=255)
                 cTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTgcg                                                      >  1:2432976/1‑70 (MQ=255)
                  ttATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTa                                                                  <  1:2438755/57‑1 (MQ=255)
                  ttATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTa                                                                  <  1:492727/57‑1 (MQ=255)
                  ttATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAg                                                           <  1:143206/64‑1 (MQ=255)
                  ttATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCa                                                    >  1:190024/1‑71 (MQ=255)
                      gtttgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCat                                                   >  1:3748772/1‑68 (MQ=255)
                       tttgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTc                                                                       <  1:2580254/47‑1 (MQ=255)
                       tttgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTg                                               <  1:1264873/71‑1 (MQ=255)
                        ttgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCata                                                  >  1:2464493/1‑67 (MQ=255)
                        ttgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTg                                               >  1:2779436/1‑70 (MQ=255)
                         tgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCa                                             >  1:605899/1‑71 (MQ=255)
                         tgttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCa                                             >  1:2081223/1‑71 (MQ=255)
                          gttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCata                                                  >  1:3587154/1‑65 (MQ=255)
                          gttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCAc                                            >  1:1145142/1‑71 (MQ=255)
                           ttCAAAGCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCa                                                    <  1:477692/62‑1 (MQ=255)
                           ttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCa                                                    <  1:1989057/62‑1 (MQ=255)
                           ttCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCAc                                            <  1:1137790/70‑1 (MQ=255)
                            tCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTacac                                                               >  1:66281/1‑50 (MQ=255)
                             cAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCata                                                  >  1:1116418/1‑62 (MQ=255)
                             cAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCata                                                  >  1:3290294/1‑62 (MQ=255)
                             cAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCt                                          <  1:882958/70‑1 (MQ=255)
                              aaaCCGTTGAACCTGGATGTACTGAAAACACATTTAAGTTAGTTaca                                                                >  1:3078444/1‑47 (MQ=255)
                               aaCCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATAtt                                                <  1:1060534/62‑1 (MQ=255)
                                aCCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTgc                                                       <  1:162750/54‑1 (MQ=255)
                                aCCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAg                                       >  1:2127914/1‑70 (MQ=255)
                                 ccGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTa                                                                  >  1:614713/1‑42 (MQ=255)
                                 ccGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTa                                     <  1:1178041/71‑1 (MQ=255)
                                 ccGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTa                                     <  1:1312725/71‑1 (MQ=255)
                                 ccGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTa                                     <  1:346991/71‑1 (MQ=255)
                                    ttGACCCTGGATGTACTGAAAACACATTTAAGTTAGTTACACCAAGTTGCGCATATTGCACCTCAGt                                      <  1:742446/67‑1 (MQ=255)
                                    ttGACCCTGGATGTACTGAAAACACATTTAAGTTAGTTACACCAAGTTGCGCATATTGCACCTCAGt                                      <  1:3946439/67‑1 (MQ=255)
                                      gACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCa                                                             <  1:2640827/42‑1 (MQ=255)
                                      gACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGt                                                          >  1:3542506/1‑45 (MQ=255)
                                      gACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGc                                 >  1:2459251/1‑70 (MQ=255)
                                        cccTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAg                                                           >  1:249361/1‑42 (MQ=255)
                                        cccTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCAc                              <  1:1257040/71‑1 (MQ=255)
                                        cccTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCAc                              >  1:1387746/1‑71 (MQ=255)
                                         ccTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTa                                     <  1:217738/63‑1 (MQ=255)
                                          cTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGtt                                                         >  1:1002473/1‑42 (MQ=255)
                                          cTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCt                            >  1:2774143/1‑71 (MQ=255)
                                             gATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCtt                           <  1:1076163/69‑1 (MQ=255)
                                             gATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTg                          >  1:3578512/1‑70 (MQ=255)
                                              aTGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTgc                                                       >  1:3377341/1‑40 (MQ=255)
                                              aTGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTgc                                                       >  1:3137743/1‑40 (MQ=255)
                                              aTGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGa                         >  1:1939147/1‑70 (MQ=255)
                                                 tACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTAt                                    >  1:313961/1‑56 (MQ=255)
                                                 tACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTAt                                    >  1:2618291/1‑56 (MQ=255)
                                                 tACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTAt                                    >  1:3903025/1‑56 (MQ=255)
                                                   cTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCata                                                  >  1:2695048/1‑40 (MQ=255)
                                                   cTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCTCCACCTTGATATCGa                   <  1:1287696/71‑1 (MQ=255)
                                                   cTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCg                    <  1:889170/70‑1 (MQ=255)
                                                   cTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCg                    <  1:2176/70‑1 (MQ=255)
                                                   cTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGa                   <  1:2413759/71‑1 (MQ=255)
                                                    tGAAAACACATTTAAGTTAGTTACACCAAGTTGCGCATATTGCACCTCAGTATc                                   >  1:770539/1‑54 (MQ=255)
                                                    tGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCAcc                             <  1:653324/60‑1 (MQ=255)
                                                    tGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCAcc                             <  1:2349382/60‑1 (MQ=255)
                                                    tGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCg                    >  1:2751699/1‑69 (MQ=255)
                                                     gAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCAcc                             <  1:1671229/59‑1 (MQ=255)
                                                     gAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAg                 >  1:2904130/1‑71 (MQ=255)
                                                     gAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAg                 >  1:3653875/1‑71 (MQ=255)
                                                     gAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAg                 >  1:2582975/1‑71 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCAc                                            >  1:252217/1‑43 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCAc                                            >  1:3093390/1‑43 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGaa                  <  1:1580829/69‑1 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGc                >  1:2300132/1‑71 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGc                >  1:1369256/1‑71 (MQ=255)
                                                      aaaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGc                >  1:1123387/1‑71 (MQ=255)
                                                       aaaCACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCt                                          >  1:348365/1‑44 (MQ=255)
                                                          cacaTTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTg            <  1:3742309/71‑1 (MQ=255)
                                                             aTTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCa                               >  1:951910/1‑49 (MQ=255)
                                                                tAAGTTAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGaaaaa       >  1:103182/1‑70 (MQ=255)
                                                                tAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGaaaaa       <  1:3298692/70‑1 (MQ=255)
                                                                tAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGaaaaa       <  1:2961779/70‑1 (MQ=255)
                                                                tAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGaaa         >  1:3152966/1‑68 (MQ=255)
                                                                 aaGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAaat      <  1:85406/70‑1 (MQ=255)
                                                                 aaGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAaat      <  1:266936/70‑1 (MQ=255)
                                                                  aGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAaataa    <  1:812816/71‑1 (MQ=255)
                                                                  aGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAaataa    <  1:2706180/71‑1 (MQ=255)
                                                                   gTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAaataa    >  1:551907/1‑70 (MQ=255)
                                                                    tCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAAataata  >  1:3437052/1‑71 (MQ=255)
                                                                    tCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAAataata  >  1:271826/1‑71 (MQ=255)
                                                                     |                                                                     
TTAAGTTGCGGCATGAACTTATGTTTGTTCAAACCGTTGACCCTGGATGTACTGAAAACACATTTAAGTCAGTTACACCAAGTTGCGCATATTGCACCTCAGTATCGCCACCTTGATATCGAAGCCCTGAAAAATAATA  >  minE/1563138‑1563276

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: