Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,011,495 T→C 16.5% V308A (GTG→GCG)  rstB → sensory histidine kinase in two‑component regulatory system with RstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,011,4950TC16.5% 124.9 / 13.6 78V308A (GTG→GCG) rstBsensory histidine kinase in two‑component regulatory system with RstA
Reads supporting (aligned to +/- strand):  ref base T (41/24);  new base C (3/10);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAAT  >  minE/1011429‑1011559
                                                                  |                                                                
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACTCCCGATAAAACGGTACGGATAAAAACGCTCGt                                                                  >  1:2296773/1‑67 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGt                                                                  >  1:3151402/1‑67 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGt                                                                  >  1:1617407/1‑67 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGt                                                                  >  1:1517188/1‑67 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGc                                                                <  1:203768/69‑1 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGc                                                                <  1:986991/69‑1 (MQ=255)
tGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGc                                                                <  1:2375596/69‑1 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:2017954/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:1641745/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:1632056/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:2331829/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:646632/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:1532656/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:627842/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:2780216/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:491746/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:985740/1‑68 (MQ=255)
    aaCGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAg                                                             >  1:993591/1‑68 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCa                                                               >  1:2042750/1‑62 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:282089/1‑69 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:379720/1‑69 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:3209153/1‑69 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:2220570/1‑69 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:130253/1‑69 (MQ=255)
        cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCAt                                                        >  1:1580998/1‑69 (MQ=255)
           cTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTAt                                                     >  1:1292955/1‑69 (MQ=255)
           cTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTAt                                                     >  1:2673428/1‑69 (MQ=255)
           cTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTAt                                                     >  1:866123/1‑69 (MQ=255)
                       cAGGCAGTAACGCCTGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:296539/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:2967706/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:2980703/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:1381189/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:744842/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:1575672/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:751651/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGat                                         >  1:2393627/1‑69 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1264419/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1287067/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1274053/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:849800/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:225986/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1751078/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2384452/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:520287/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:3042967/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2788712/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2126103/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2886068/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2893155/68‑1 (MQ=255)
                       cAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1697704/68‑1 (MQ=255)
                        aGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:1775959/67‑1 (MQ=255)
                        aGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGCGCAAGGCCATTATGCGGCGTTGGa                                          <  1:2429379/67‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCTCtt                                 <  1:2713410/69‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:3096250/69‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:1991432/69‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:3317450/69‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:2073417/69‑1 (MQ=255)
                               aaCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:1918264/69‑1 (MQ=255)
                                aCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGa                                          <  1:3199951/59‑1 (MQ=255)
                                aCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 >  1:443775/1‑68 (MQ=255)
                                aCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 >  1:1137315/1‑68 (MQ=255)
                                aCGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 >  1:2896039/1‑68 (MQ=255)
                                  gCCCGATAAAACGGTACGGATAAAAACGCTCGTGc                                                                <  1:859381/35‑1 (MQ=255)
                                          aaaCGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGata                                        <  1:1570394/51‑1 (MQ=255)
                                               gTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTg                 >  1:2959359/1‑69 (MQ=255)
                                               gTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTg                 >  1:1382128/1‑69 (MQ=255)
                                               gTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTg                 >  1:2880910/1‑69 (MQ=255)
                                               gTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTg                 >  1:2365258/1‑69 (MQ=255)
                                                    gATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATaa            >  1:539657/1‑69 (MQ=255)
                                                    gATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATaa            >  1:1214685/1‑69 (MQ=255)
                                                    gATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATaa            >  1:1515031/1‑69 (MQ=255)
                                                    gATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATaa            >  1:2307433/1‑69 (MQ=255)
                                                    gATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATaa            >  1:3432533/1‑69 (MQ=255)
                                                        aaaaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCtt                                 <  1:505354/44‑1 (MQ=255)
                                                        aaaaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAAtt          >  1:2450985/1‑67 (MQ=255)
                                                        aaaaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTg        >  1:2785049/1‑69 (MQ=255)
                                                        aaaaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTg        >  1:3400920/1‑69 (MQ=255)
                                                        aaaaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTg        >  1:1154401/1‑69 (MQ=255)
                                                          aaCGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGc       <  1:3020609/68‑1 (MQ=255)
                                                              cTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAAt  <  1:2615857/69‑1 (MQ=255)
                                                              cTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAAt  <  1:3404979/69‑1 (MQ=255)
                                                              cTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAAt  <  1:109783/69‑1 (MQ=255)
                                                                  |                                                                
TGTCAACGCATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAAT  >  minE/1011429‑1011559

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: