Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 804,017 G→A 51.6% R246R (CGG→CGA dhaR → predicted DNA‑binding transcriptional regulator, dihydroxyacetone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE804,0170GA51.6% ‑3.7 / 53.2 60R246R (CGG→CGAdhaRpredicted DNA‑binding transcriptional regulator, dihydroxyacetone
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (17/14);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCG  >  minE/803952‑804080
                                                                 |                                                               
ttGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGTGACGGGCa                                                              >  1:2919863/1‑69 (MQ=255)
ttGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGTGACGGGCa                                                              >  1:2923751/1‑69 (MQ=255)
ttGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCa                                                              >  1:237378/1‑69 (MQ=255)
ttGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCa                                                              >  1:3408486/1‑69 (MQ=255)
       ttattaATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCa                                                          <  1:232387/66‑1 (MQ=255)
               gCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCtt                                                <  1:77267/68‑1 (MQ=255)
               gCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCtt                                                <  1:541111/68‑1 (MQ=255)
               gCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCtt                                                <  1:45875/68‑1 (MQ=255)
               gCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCtt                                                <  1:276127/68‑1 (MQ=255)
               gCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCtt                                                <  1:1625923/68‑1 (MQ=255)
                       ggcgCGGGACTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCa                                                          <  1:681833/50‑1 (MQ=255)
                          gcgGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACccgc                                    >  1:244846/1‑69 (MQ=255)
                                cTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTAtt                              <  1:2497559/69‑1 (MQ=255)
                                 ttGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTAtt                              <  1:1135275/68‑1 (MQ=255)
                                   gcgcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:3427097/68‑1 (MQ=255)
                                   gcgcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1334529/68‑1 (MQ=255)
                                     gcgttGACGCGACGGCAAGTCAGGGTCGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1752054/63‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCATGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:2307987/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:301860/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1655969/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1685450/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:269073/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:2677717/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:3427311/66‑1 (MQ=255)
                                     gcCTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1222761/66‑1 (MQ=255)
                                      ccTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:1918640/65‑1 (MQ=255)
                                      ccTTGACGCGACGGCAAGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:2222546/65‑1 (MQ=255)
                                           aCGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTAtt                              <  1:105841/58‑1 (MQ=255)
                                                    cAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgg          <  1:1271817/69‑1 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTa                                >  1:479988/1‑46 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCa                    >  1:1842801/1‑58 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:589786/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:605446/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:413355/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:1330681/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:319176/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:684003/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:3160736/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:3102226/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:725514/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:1384008/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:1419155/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:2707744/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:3047804/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:2876965/1‑67 (MQ=255)
                                                     aaGTCAGGGACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAgt          >  1:1775074/1‑67 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGcaca      >  1:3086086/1‑68 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGcaca      >  1:312842/1‑68 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGcaca      >  1:335680/1‑68 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:1371617/1‑69 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:4993/1‑69 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:2757087/1‑69 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:2796872/1‑69 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:1201772/1‑69 (MQ=255)
                                                         cAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACAt     >  1:2217507/1‑69 (MQ=255)
                                                          aGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATc    <  1:1025030/69‑1 (MQ=255)
                                                          aGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATc    <  1:1259902/69‑1 (MQ=255)
                                                            ggACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAAGAAGCAATAAAACAGGCACATCCg  <  1:1359683/69‑1 (MQ=255)
                                                            ggACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCg  <  1:1774821/69‑1 (MQ=255)
                                                             gACGAGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                            <  1:3172390/42‑1 (MQ=255)
                                                                 ggCAATCACTGAACTCTTAACGTTACCCGCCGTAtt                              <  1:3039700/36‑1 (MQ=255)
                                                                 |                                                               
TTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCG  >  minE/803952‑804080

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: