Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 92,124 G→T 39.4% L276L (CTG→CTT lpxC → UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE92,1240GT39.4% 46.9 / 54.4 66L276L (CTG→CTTlpxCUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base T (19/7);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 7.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

TACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGC  >  minE/92070‑92183
                                                      |                                                           
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:1019343/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:875195/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:819066/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:378198/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:1133431/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:3214288/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:2302250/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:2153587/1‑69 (MQ=255)
tACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGc                                               >  1:204294/1‑69 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:992448/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:2977642/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:2457139/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:2450446/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:2448997/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:1140258/69‑1 (MQ=255)
       tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGaa                                        <  1:291614/69‑1 (MQ=255)
        ataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGaa                                        <  1:220986/68‑1 (MQ=255)
          acaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:3085873/66‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:504932/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:3091871/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:2349877/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:2989938/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:50765/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:574192/68‑1 (MQ=255)
          aaaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAata                                      <  1:1027837/68‑1 (MQ=255)
           aaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAatat                                     >  1:1168011/1‑68 (MQ=255)
           aaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAatat                                     >  1:3425152/1‑68 (MQ=255)
           aaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAatat                                     >  1:2930797/1‑68 (MQ=255)
           aaTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAAAAGGAAGACTGGGAatat                                     >  1:84334/1‑68 (MQ=255)
              ccGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGa                                 <  1:1776961/69‑1 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAATACAGGAAGCCTGCGAATATGTGAc                                >  1:3249161/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGGATATGTGAc                                >  1:121838/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:3302302/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:2955809/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:3066856/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1003648/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:638319/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:743378/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:931287/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1600639/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:2872488/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:2797483/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1907686/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1173207/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1380715/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1983851/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:1423754/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:2923839/1‑69 (MQ=255)
               cGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTTGCAAAACAGGAAGCCTGGGAATATGTGAc                                >  1:2224614/1‑69 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:399485/69‑1 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:1862603/69‑1 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:1036309/69‑1 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:1112917/69‑1 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:2441562/69‑1 (MQ=255)
                                aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAg               <  1:2584343/69‑1 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:1580584/1‑69 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:539221/1‑69 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:199871/1‑69 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:68437/1‑69 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:1949240/1‑69 (MQ=255)
                                     ctgctgCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTg          >  1:1768786/1‑69 (MQ=255)
                                             ggCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTg    <  1:3255222/67‑1 (MQ=255)
                                             ggCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTg    <  1:2456466/67‑1 (MQ=255)
                                             ggCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTg    <  1:189796/67‑1 (MQ=255)
                                             ggCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTg    <  1:2682133/67‑1 (MQ=255)
                                             ggCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGc  >  1:1404052/1‑69 (MQ=255)
                                                      |                                                           
TACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGC  >  minE/92070‑92183

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