Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 162,233 A→G 48.7% T316A (ACG→GCG)  cdaR → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE162,2330AG48.7% ‑5.8 / 33.2 39T316A (ACG→GCG) cdaRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/3);  new base G (6/13);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 1.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

TTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGGCGGCGTGGTTTCGCCACAA  >  minE/162167‑162288
                                                                  |                                                       
ttACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACg                                                       >  1:1000042/1‑69 (MQ=255)
ttACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACg                                                       >  1:482609/1‑69 (MQ=255)
  aCTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGa                                                      >  1:291111/1‑68 (MQ=255)
  aCTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGATGAAAACGa                                                      >  1:2063285/1‑68 (MQ=255)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2750812/1‑69 (MQ=39)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2840280/1‑69 (MQ=39)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:2138601/1‑69 (MQ=39)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:307219/1‑69 (MQ=39)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:1525687/1‑69 (MQ=39)
                          gCAGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             >  1:699573/1‑69 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTTAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:704628/67‑1 (MQ=255)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:3390959/67‑1 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2960168/67‑1 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:67224/67‑1 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2857349/67‑1 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:2832687/67‑1 (MQ=39)
                            aGGCCAACGAACTGGCGCGACCGCTGGCACGGCTGAAAGCGATGGACAATAACGGCTTGCTgcggcg                             <  1:1289811/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: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:1702766/1‑69 (MQ=255)
                                       cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:522211/1‑69 (MQ=255)
                                       cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:605677/1‑69 (MQ=255)
                                       cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:3138386/1‑69 (MQ=255)
                                       cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:1611935/1‑69 (MQ=255)
                                       cTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGgcggcg                >  1:1324543/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)
                                                                  |                                                       
TTACTCGACAGTTTGCGTGGCGACTGGCAGGCCAACGAACTGGCGCGACCGCTGGCGCGGCTGAAAACGATGGACAATAACGGCTTGCTGCGACGAACGCTGGCGGCGTGGTTTCGCCACAA  >  minE/162167‑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: