Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,258 A→C 9.1% E176A (GAA→GCA) ‡ 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,2580AC 9.1% 130.6 / 8.3 55E176A (GAA→GCA) ‡yejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base C (3/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTG  >  minE/1408195‑1408328
                                                               |                                                                      
aTCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAattatt                                                                  <  1:717644/70‑1 (MQ=255)
aTCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAattatt                                                                  <  1:399649/70‑1 (MQ=255)
aTCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAAttatta                                                                 <  1:1971364/71‑1 (MQ=255)
   aGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTaa                                                                >  1:1110707/1‑69 (MQ=255)
   aGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAAtt                                                              <  1:3020082/71‑1 (MQ=255)
   aGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAAtt                                                              >  1:2278228/1‑71 (MQ=255)
   aGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAAtt                                                              >  1:588467/1‑71 (MQ=255)
        tcCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGa                                                         >  1:1337274/1‑71 (MQ=255)
            gcggcgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                     <  1:39836/71‑1 (MQ=255)
                cgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGa                                                 >  1:335207/1‑71 (MQ=255)
                 gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAatt                                                                     >  1:568819/1‑50 (MQ=255)
                 gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAatt                                                                     >  1:2551363/1‑50 (MQ=255)
                 gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                     <  1:1887324/66‑1 (MQ=255)
                 gAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGaa                                                     <  1:3605202/66‑1 (MQ=255)
                   aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTaa                                                                >  1:1124960/1‑53 (MQ=255)
                   aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGacca                                              >  1:2566559/1‑71 (MQ=255)
                   aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGacc                                               >  1:2130392/1‑70 (MQ=255)
                   aCGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGac                                                >  1:1864905/1‑69 (MQ=255)
                   aCGGCAGCGGGTGATGATAGCGATGGCGCTGTTAACGCGACCGGAATTATTaa                                                                >  1:3733882/1‑53 (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)
                         gCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGa                                                                        <  1:3835204/39‑1 (MQ=255)
                         gCGAGTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGcgc                                        >  1:1874009/1‑69 (MQ=255)
                          cgagTGATGATTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGCGCt                                       >  1:506781/4‑71 (MQ=255)
                          cGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGa                                                                        <  1:3427926/38‑1 (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:1230163/70‑1 (MQ=255)
                                       cGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGt                           <  1:1986457/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:2448138/71‑1 (MQ=255)
                                                               aaTTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTg  <  1:2766835/71‑1 (MQ=255)
                                                               |                                                                      
ATCAGCTCTCCGGCGGCGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTG  >  minE/1408195‑1408328

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: