Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,448,005 G→C 41.8% V137V (GTG→GTC atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,448,0050GC41.8% 36.5 / 81.3 79V137V (GTG→GTCatoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base G (8/38);  new base C (23/10);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 4.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTC  >  minE/1447940‑1448067
                                                                 |                                                              
ccTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATg                                                             >  1:3262995/1‑69 (MQ=255)
ccTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATg                                                             >  1:578882/1‑69 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:446577/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:1896519/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:1865800/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:2380118/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:1949621/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:342782/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:1412770/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:589131/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:733369/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:2276625/69‑1 (MQ=255)
        gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:3078092/69‑1 (MQ=255)
         cTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:1199131/68‑1 (MQ=255)
          tataTCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACaaa                                                     <  1:2306267/67‑1 (MQ=255)
          tataTCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATa                                                   <  1:510787/69‑1 (MQ=255)
          tataTCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATa                                                   <  1:2867908/69‑1 (MQ=255)
                     aatGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGtt                                        <  1:1380926/69‑1 (MQ=255)
                     aatGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGtt                                        <  1:3342928/69‑1 (MQ=255)
                     aatGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGtt                                        <  1:3019717/69‑1 (MQ=255)
                      atGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAAt                                                    >  1:1749955/1‑56 (MQ=255)
                      atGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGttt                                       >  1:2922666/1‑69 (MQ=255)
                      atGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGttt                                       >  1:2793410/1‑69 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCTAAGTGATGCGTACa                                                       <  1:2743899/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:2377536/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:2964222/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:271060/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:306429/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:1886682/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:2108205/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:210550/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:3242119/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:1842944/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:1778715/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:3332757/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:445920/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:499054/51‑1 (MQ=255)
                        gTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACa                                                       <  1:1331356/51‑1 (MQ=255)
                                 aCCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggca                              >  1:3228305/1‑67 (MQ=255)
                                 aCCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggca                              >  1:2189962/1‑67 (MQ=255)
                                 aCCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggca                              >  1:1980572/1‑67 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:1187411/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:795925/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:2751140/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:888163/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:247958/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:2423988/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:2012068/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:1912230/69‑1 (MQ=255)
                                           cAGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGAAGGTGcgcg                  <  1:738754/69‑1 (MQ=255)
                                            aGATCACCCTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGcgcg                  <  1:3196113/68‑1 (MQ=255)
                                               tCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGa              <  1:2045655/69‑1 (MQ=255)
                                               tCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGa              <  1:1431051/69‑1 (MQ=255)
                                               tCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGa              <  1:1553962/69‑1 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTcag                                 >  1:2146092/1‑47 (MQ=39)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:458754/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:425732/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1482120/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1437915/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1388857/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:6941/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:706705/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1279222/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:123288/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:917085/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:447125/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:149379/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:2724978/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1516086/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:337266/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:1680629/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:228005/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:2409446/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGatat           >  1:2948528/1‑69 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGata            >  1:1508011/1‑68 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGata            >  1:649477/1‑68 (MQ=255)
                                                  cccTGGTCGCGAAGTCATGCGTACAAATACCCCTTTGGTTTATTCAGGCAGGCAGGTGCGCGGCGata            >  1:1076057/1‑68 (MQ=255)
                                                            gAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTc  <  1:2709577/68‑1 (MQ=255)
                                                             aaGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGGGCGCGGCGATATTTTGAATTc  <  1:1674981/67‑1 (MQ=255)
                                                                 |                                                              
CCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTC  >  minE/1447940‑1448067

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: