Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,885,303 G→A 54.0% Q255Q (CAG→CAA barA → hybrid sensory histidine kinase, in two‑component regulatory system with UvrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,885,3030GA54.0% 12.6 / 96.9 113Q255Q (CAG→CAAbarAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base G (37/15);  new base A (37/24);  total (74/39)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 128 total aligned reads displayed.

ACGGTATCAACTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTG  >  minE/1885251‑1885370
                                                    |                                                                   
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:620218/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:31404/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:328271/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:3083988/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:2920590/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:2424815/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:2345160/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:229851/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:581536/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:1914622/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:1892410/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:1642801/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:856500/1‑68 (MQ=255)
aCGGTATCAACTCGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:1218025/1‑68 (MQ=255)
aCGGTATCAACACGATGGCAATGTCGCTGGCGGCTTATCACGAAGAGATGCAACACAATATCGACCAg                                                      >  1:381592/1‑68 (MQ=255)
  ggTATCAACTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCg                                                   <  1:2913718/69‑1 (MQ=255)
  ggTATCAACTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCg                                                   <  1:2753309/69‑1 (MQ=255)
  ggTATCAACTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCg                                                   <  1:2553955/69‑1 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGATGAGATGCAGCACAATATCGACCAGGc                                                    >  1:2540350/1‑62 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAat                                                               >  1:1209209/1‑51 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGc                                                    >  1:2966504/1‑62 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGc                                                    >  1:343718/1‑62 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGc                                                    >  1:2572742/1‑62 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:1681723/1‑69 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:2200291/1‑69 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:1781816/1‑69 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:2974572/1‑69 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:265565/1‑69 (MQ=255)
        aaCTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTcc                                             >  1:1012231/1‑69 (MQ=255)
                   aaTGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGa                                  <  1:3352639/69‑1 (MQ=255)
                       tCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACg                              >  1:707461/1‑69 (MQ=255)
                       tCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACg                              >  1:368895/1‑69 (MQ=255)
                            ggctgCTTATCACGAAGTGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTgg                          >  1:2168026/1‑68 (MQ=255)
                            ggctgCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTgg                          >  1:2643329/1‑68 (MQ=255)
                                      cACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATgg                 <  1:1089017/67‑1 (MQ=255)
                                      cACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaa               >  1:1190963/1‑69 (MQ=255)
                                      cACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaa               >  1:858734/1‑69 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:671252/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:267489/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:1040420/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:1126437/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:2097116/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:114001/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:2942001/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:3042186/1‑36 (MQ=255)
                                      cACGAAGAGATGCAACACAATATCGACCAGGCGACg                                                >  1:1542542/1‑36 (MQ=255)
                                       aCGAAGAGATGTAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:1575240/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:2913878/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:1034751/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:1251612/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:101312/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:1003889/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:2585290/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:2434486/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:899387/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:836559/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:1657254/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaaa              >  1:943706/1‑69 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaa               >  1:1054027/1‑68 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGaa               >  1:1905180/1‑68 (MQ=255)
                                       aCGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGATGGAGCAGATGGaa               >  1:1549187/1‑68 (MQ=255)
                                            gagaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:2484796/69‑1 (MQ=255)
                                            gagaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:1273536/69‑1 (MQ=255)
                                            gagaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:2255608/69‑1 (MQ=255)
                                            gagaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:2781407/69‑1 (MQ=255)
                                            gagaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:467104/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:638381/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:34720/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:1011562/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:797858/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:3302086/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:919022/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:1489687/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:2072023/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:1533362/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:2291549/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:2395070/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:2766328/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:3033006/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:3042889/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:1281781/69‑1 (MQ=255)
                                            gagaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:3103116/69‑1 (MQ=255)
                                             agaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAg         <  1:2896964/68‑1 (MQ=255)
                                             agaTGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGa        >  1:3042839/1‑69 (MQ=255)
                                             agaTGCAACACAATATCGACCAGGCGACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:487750/68‑1 (MQ=255)
                                             agaTGCAACACAATATCGACCAGGCCACGTCCGATCTGCGTGAGACGCTGGAGCAAATGGAGATTCAg         <  1:635745/68‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:1483329/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:943820/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:1213445/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:3220265/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:1995382/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:357507/56‑1 (MQ=255)
                                                  cagcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGa                <  1:3130186/56‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTTAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:3319507/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAAGGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:2121460/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:1900904/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:1518173/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:3411190/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:3324607/68‑1 (MQ=255)
                                                    gcaCAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTg  <  1:1149039/68‑1 (MQ=255)
                                                    |                                                                   
ACGGTATCAACTCGATGGCAATGTCGCTGGCTGCTTATCACGAAGAGATGCAGCACAATATCGACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTG  >  minE/1885251‑1885370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: