Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,183 G→A 14.1% D190N (GAT→AAT)  yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,449,1830GA14.1% 111.1 / 6.6 63D190N (GAT→AAT) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/32);  new base A (4/5);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAG  >  minE/2449116‑2449249
                                                                   |                                                                  
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTg                                                                    >  1:17899/1‑68 (MQ=255)
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                                                   >  1:1980186/1‑69 (MQ=255)
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                                                   >  1:674143/1‑69 (MQ=255)
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                                                   >  1:2251165/1‑69 (MQ=255)
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                                                   >  1:2586502/1‑69 (MQ=255)
cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                                                   >  1:2918755/1‑69 (MQ=255)
     atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                                              <  1:2537010/69‑1 (MQ=255)
     atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                                              <  1:1299031/69‑1 (MQ=255)
     atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                                              <  1:2928382/69‑1 (MQ=255)
     atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                                              <  1:3057307/69‑1 (MQ=255)
     atcTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTATCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:3109224/69‑1 (MQ=255)
      tcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                                              <  1:1299442/68‑1 (MQ=255)
          cAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATcc                                                         >  1:2246788/1‑69 (MQ=255)
                   tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc                                                  <  1:2165061/67‑1 (MQ=255)
                   tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc                                                  <  1:1473579/67‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCTAACGCCCGAAACCAACTGATGCAg                                                              <  1:2078615/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAATCCAACTGATGCAg                                                              <  1:2858125/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:1609168/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              >  1:1631100/1‑51 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              >  1:3230825/1‑51 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:412390/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              >  1:590504/1‑51 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:2663206/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:612463/51‑1 (MQ=255)
                       cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                                              <  1:2497355/51‑1 (MQ=255)
                           tCTGGAAGTTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                                        >  1:3361057/1‑69 (MQ=255)
                           tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                                        >  1:3177122/1‑69 (MQ=255)
                           tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                                        >  1:1494341/1‑69 (MQ=255)
                           tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                                        >  1:1228663/1‑69 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGATACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:665338/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:1906360/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:3194067/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:11116/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:1160450/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATACCGGGCGCTACAACATGATCTCAGt                                 <  1:2326511/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGACGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:1222405/69‑1 (MQ=255)
                                  ggTAAGCACGCTCTCGAACGCCCGAAAACAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:2556804/69‑1 (MQ=255)
                                   gTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt                                 <  1:1788888/68‑1 (MQ=255)
                                    tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa                               >  1:96538/1‑69 (MQ=255)
                                    tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa                               >  1:879272/1‑69 (MQ=255)
                                    tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa                               >  1:818564/1‑69 (MQ=255)
                                    tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa                               >  1:1483174/1‑69 (MQ=255)
                                    tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa                               >  1:1840727/1‑69 (MQ=255)
                                     aaGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGAt                              <  1:2368894/69‑1 (MQ=255)
                                        cACGCTCTCGAACGCCCGAACCCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt                           >  1:2658903/1‑69 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc                                                  >  1:121899/1‑46 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt                           >  1:2570388/1‑69 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt                           >  1:2937750/1‑69 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt                           >  1:3059592/1‑69 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt                           >  1:1396788/1‑69 (MQ=255)
                                        cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCAGGGCGCTAAAACATGATCTCAGTGATGtt                           >  1:1782639/1‑69 (MQ=255)
                                         aCGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTa                          <  1:1982651/69‑1 (MQ=255)
                                         aCGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTa                          <  1:1566981/69‑1 (MQ=255)
                                                           aaCCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGaa        <  1:1489544/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:2472533/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:183824/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:2357125/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  >  1:1281545/1‑69 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:1378014/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:1634217/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:1693596/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:2779690/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:1832340/69‑1 (MQ=255)
                                                                 cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg  <  1:2913066/69‑1 (MQ=255)
                                                                   |                                                                  
CGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAG  >  minE/2449116‑2449249

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: