Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,169 G→A 29.3% R185H (CGC→CAC)  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,1690GA29.3% 79.4 / 32.6 78R185H (CGC→CAC) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base A (18/5);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTA  >  minE/2449105‑2449225
                                                                |                                                        
tGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGccc                                                       >  1:2997321/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACAccc                                                       >  1:1152753/1‑65 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:1337024/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:2908961/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:2698910/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:3219560/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:3027641/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:201694/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:2035343/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:1279500/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:1231485/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:2187993/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:3096494/1‑68 (MQ=255)
   cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:914199/1‑68 (MQ=255)
   cTGATAACCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                                                    >  1:1991/1‑68 (MQ=255)
     gATAGCCGATCATCTGCATGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:2002059/69‑1 (MQ=255)
     gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:263278/69‑1 (MQ=255)
     gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:481784/69‑1 (MQ=255)
     gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:3292829/69‑1 (MQ=255)
     gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:2816202/69‑1 (MQ=255)
     gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:3062869/69‑1 (MQ=255)
         gCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACt                                             >  1:3370659/1‑69 (MQ=255)
         gCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACt                                             >  1:2115206/1‑69 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTg                                            >  1:17899/1‑68 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                           >  1:2918755/1‑69 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                           >  1:1980186/1‑69 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                           >  1:674143/1‑69 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                           >  1:2586502/1‑69 (MQ=255)
           cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                                           >  1:2251165/1‑69 (MQ=255)
                atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                      <  1:2537010/69‑1 (MQ=255)
                atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                      <  1:2928382/69‑1 (MQ=255)
                atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                      <  1:3057307/69‑1 (MQ=255)
                atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg                                      <  1:1299031/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)
                              tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                                                 <  1:2669019/44‑1 (MQ=255)
                              tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc                          <  1:2165061/67‑1 (MQ=255)
                              tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc                          <  1:1473579/67‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCTAACGCCCGAAACCAACTGATGCAg                                      <  1:2078615/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGaaa                                                   >  1:2569504/1‑38 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAATCCAACTGATGCAg                                      <  1:2858125/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      <  1:2497355/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      >  1:3230825/1‑51 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      <  1:2663206/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      <  1:412390/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      >  1:590504/1‑51 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      >  1:1631100/1‑51 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      <  1:1609168/51‑1 (MQ=255)
                                  cGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAg                                      <  1:612463/51‑1 (MQ=255)
                                      tCTGGAAGTTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                >  1:3361057/1‑69 (MQ=255)
                                      tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                >  1:3177122/1‑69 (MQ=255)
                                      tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                >  1:1228663/1‑69 (MQ=255)
                                      tCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAGTATCCCGGGCGCTACAACATGa                >  1:1494341/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:1160450/69‑1 (MQ=255)
                                             ggTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGt         <  1:11116/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:1840727/1‑69 (MQ=255)
                                               tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa       >  1:1483174/1‑69 (MQ=255)
                                               tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa       >  1:818564/1‑69 (MQ=255)
                                               tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa       >  1:879272/1‑69 (MQ=255)
                                               tAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGa       >  1:96538/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:3059592/1‑69 (MQ=255)
                                                   cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt   >  1:2937750/1‑69 (MQ=255)
                                                   cACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGtt   >  1:2570388/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)
                                                                |                                                        
TGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTA  >  minE/2449105‑2449225

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: