Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,635,541 T→C 64.4% N226N (AAT→AAC ugpB → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,635,5410TC64.4% 26.4 / 21.2 42N226N (AAT→AACugpBglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (5/10);  new base C (13/14);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTCTGCCGTTTGCCAGCAAAAACAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAGCTA  >  minE/2635477‑2635606
                                                                |                                                                 
aCGGTCTGCCGTTTGCCAGCAAAAACAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGc                                                               <  1:1462678/69‑1 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGa                                                   <  1:2486446/69‑1 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGa                                                   <  1:1933096/69‑1 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:980859/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:1016145/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:979940/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:337109/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:312943/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:2996154/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:2676853/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:26628/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:1961976/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:1020079/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:115283/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:1467981/1‑69 (MQ=255)
            ttGCCAGCAAAAACAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGa                                                   >  1:1530980/1‑69 (MQ=255)
                     aaaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:1773103/69‑1 (MQ=255)
                     aaaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:1055047/69‑1 (MQ=255)
                     aaaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:2772623/69‑1 (MQ=255)
                     aaaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:1562176/69‑1 (MQ=255)
                     aaaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:1676960/69‑1 (MQ=255)
                      aaaCAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGc                                         <  1:2907028/69‑1 (MQ=255)
                      aaaCAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGc                                         <  1:2999672/69‑1 (MQ=255)
                      aaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCTGAGCAGGTGAAACACATCg                                          <  1:3321058/68‑1 (MQ=255)
                      aaaCAACGGCTTTGACGGCACAGACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCg                                          <  1:2480595/68‑1 (MQ=255)
                       aacaacGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGcc                                        <  1:1466742/69‑1 (MQ=255)
                       aacaacGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGcc                                        <  1:2775028/69‑1 (MQ=255)
                                    aCGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCgagga                             <  1:2272833/67‑1 (MQ=255)
                                           agACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:1734601/68‑1 (MQ=255)
                                           agACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:401411/68‑1 (MQ=255)
                                           agACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:1964569/68‑1 (MQ=255)
                                           agACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:2768299/68‑1 (MQ=255)
                                           agACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:1633244/68‑1 (MQ=255)
                                            gACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:2260697/68‑1 (MQ=255)
                                            gACGCGGTGCTTGAGTTCAACAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAAGAGATGAACaa                    <  1:1734034/68‑1 (MQ=255)
                                                cgGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACaagaag                >  1:322717/1‑68 (MQ=255)
                                                cgGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAgg               >  1:733947/1‑69 (MQ=255)
                                                          gTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAg     <  1:1970751/69‑1 (MQ=255)
                                                           ttCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAg     <  1:3193690/68‑1 (MQ=255)
                                                             cAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACaaga                  >  1:2926947/1‑53 (MQ=255)
                                                             cAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAGCTa  >  1:3019830/1‑69 (MQ=255)
                                                             cAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAGCTa  >  1:2761491/1‑69 (MQ=255)
                                                                |                                                                 
ACGGTCTGCCGTTTGCCAGCAAAAACAACGGCTTTGACGGCACGGACGCGGTGCTGGAGTTCAATAAGCCGGAGCAGGTGAAACACATCGCCATGCTCGAGGAGATGAACAAGAAGGGCGACTTCAGCTA  >  minE/2635477‑2635606

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: