Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,277 A→G 10.2% E182E (GAA→GAG yejF → fused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,408,2770AG10.2% 125.1 / 6.0 49E182E (GAA→GAGyejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (3/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCG  >  minE/1408207‑1408344
                                                                      |                                                                   
gcggcgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                                     <  1:39836/71‑1 (MQ=255)
    cgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGa                                                                 >  1:335207/1‑71 (MQ=255)
     gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                                     <  1:3605202/66‑1 (MQ=255)
     gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                                     <  1:1887324/66‑1 (MQ=255)
       aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGacca                                                              >  1:2566559/1‑71 (MQ=255)
       aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGacc                                                               >  1:2130392/1‑70 (MQ=255)
       aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGac                                                                >  1:1864905/1‑69 (MQ=255)
        cgccAGCGAGTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAacca                                                              <  1:3784002/67‑1 (MQ=255)
           cAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCg                                                           >  1:1159318/1‑70 (MQ=255)
            aGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGaccaccact                                                         <  1:2966377/71‑4 (MQ=255)
            aGCGAGTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGc                                                          >  1:1209337/1‑70 (MQ=255)
             gCGAGTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGcgc                                                        >  1:1874009/1‑69 (MQ=255)
              cgagTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGCGCt                                                       >  1:506781/4‑71 (MQ=255)
                ggTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTgg                                                     <  1:1813691/71‑1 (MQ=255)
                ggTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTgg                                                     <  1:3110895/71‑1 (MQ=255)
                   gatgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACg                                                  >  1:397596/1‑71 (MQ=255)
                   gatgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACg                                                  >  1:3474781/1‑71 (MQ=255)
                    atgatTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGCGCTGGatg                                                  <  1:2923797/70‑3 (MQ=255)
                     tgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtc                                                >  1:479258/1‑71 (MQ=255)
                        ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAAcc                                                                   <  1:573229/49‑1 (MQ=255)
                        ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGa                                                    >  1:2963857/1‑64 (MQ=255)
                        ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtc                                                >  1:3608624/1‑68 (MQ=255)
                         tGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAAcc                                                                   <  1:1703096/48‑1 (MQ=255)
                          gCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGt                                           >  1:1692609/1‑71 (MQ=255)
                           cGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGt                                           <  1:1986457/70‑1 (MQ=255)
                           cGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGt                                           <  1:1230163/70‑1 (MQ=255)
                           cGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTc                                          <  1:465261/71‑1 (MQ=255)
                            gATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTc                                          >  1:3569246/1‑70 (MQ=255)
                             aTGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGac                                                                <  1:74249/47‑1 (MQ=255)
                              tGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAg                                       >  1:2500184/1‑71 (MQ=255)
                               ggCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                                     >  1:2045517/1‑40 (MQ=255)
                                gcgcTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGgc                                     <  1:2669850/71‑1 (MQ=255)
                                  gcTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtc                                                >  1:1251168/1‑58 (MQ=255)
                                     gTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTc                                          >  1:773352/1‑61 (MQ=255)
                                      ttAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTg                                                      <  1:2404572/48‑1 (MQ=255)
                                      ttAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGa                                <  1:3698854/70‑1 (MQ=255)
                                       tAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGAtt                              >  1:1948707/1‑71 (MQ=255)
                                         aCGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGAttt                             <  1:913819/70‑1 (MQ=255)
                                           gcgACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTAc                          <  1:929024/71‑1 (MQ=255)
                                           gcgACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTAc                          <  1:1334743/71‑1 (MQ=255)
                                               ccGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtct                                               >  1:2922257/1‑46 (MQ=255)
                                                 ggAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGt                    <  1:3439397/71‑1 (MQ=255)
                                                  gAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGtt                   <  1:3076783/71‑1 (MQ=255)
                                                   aaTTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTg                  <  1:2766835/71‑1 (MQ=255)
                                                   aaTTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTg                  <  1:2448138/71‑1 (MQ=255)
                                                     ttattaATTGCCAATGAACCGACCACCGCACTGGACGTCTCTGTCCAgg                                      <  1:438210/49‑1 (MQ=255)
                                                       attaATTGCCGATGAACCGACCACCGCACTGGACGTctct                                             >  1:451486/1‑40 (MQ=255)
                                                        ttaATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGcgcg              <  1:3823627/70‑1 (MQ=255)
                                                         taATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGaa            >  1:1096674/1‑71 (MQ=255)
                                                           aTTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACt          <  1:3293405/71‑1 (MQ=255)
                                                                  aTGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTa                           >  1:1044962/1‑47 (MQ=255)
                                                                  aTGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGc   <  1:339931/71‑1 (MQ=255)
                                                                   tGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCg  <  1:2272539/71‑1 (MQ=255)
                                                                      |                                                                   
GCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCG  >  minE/1408207‑1408344

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: