Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 162,223 G→A 45.2% A312A (GCG→GCA cdaR → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE162,2230GA45.2% 5.4 / 30.5 42A312A (GCG→GCAcdaRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (20/3);  new base A (6/13);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 3.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGGCGGCGTGGTTTCGCCACAA  >  minE/162160‑162288
                                                               |                                                                 
aCCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGc                                                               >  1:1052497/1‑68 (MQ=255)
aCCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGc                                                               >  1:1862797/1‑68 (MQ=255)
aCCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGc                                                               >  1:3381898/1‑68 (MQ=255)
aCCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGc                                                               >  1:1487679/1‑68 (MQ=255)
       ttACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACg                                                       >  1:482609/1‑69 (MQ=255)
       ttACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACg                                                       >  1:1000042/1‑69 (MQ=255)
         aCTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGa                                                      >  1:291111/1‑68 (MQ=255)
         aCTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGATGAAAACGa                                                      >  1:2063285/1‑68 (MQ=255)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:699573/1‑69 (MQ=39)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:1525687/1‑69 (MQ=39)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:307219/1‑69 (MQ=39)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2138601/1‑69 (MQ=39)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2840280/1‑69 (MQ=39)
                                 gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2750812/1‑69 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTTAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:704628/67‑1 (MQ=255)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:67224/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:3390959/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2832687/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2960168/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2857349/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2290912/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:197122/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:1908447/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:1289811/67‑1 (MQ=39)
                                   aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2771442/67‑1 (MQ=39)
                                    ggCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2448465/66‑1 (MQ=39)
                                      ccAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGcg                                <  1:234220/61‑1 (MQ=255)
                                       cAACGAACTGGCGCGACCGCTGGCGCGGCTTAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCt                       <  1:923367/69‑1 (MQ=255)
                                       cAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCt                       >  1:1976369/1‑69 (MQ=255)
                                       cAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCt                       <  1:3055377/69‑1 (MQ=255)
                                        aaCGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2880726/62‑1 (MQ=38)
                                           gAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTggcg                   >  1:631234/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGcgacg                             >  1:1736110/1‑56 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:1611935/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:522211/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:605677/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:1702766/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:1324543/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:3138386/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:780534/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:306036/1‑69 (MQ=255)
                                              cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:95542/1‑69 (MQ=255)
                                                               gcgGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGGCGGCGTGGTTTCGCCACaa  >  1:1418348/1‑66 (MQ=255)
                                                               |                                                                 
ACCTGTGTTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGGCGGCGTGGTTTCGCCACAA  >  minE/162160‑162288

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