Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,280 G→A 10.6% P183P (CCG→CCA 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,2800GA10.6% 127.4 / 6.3 47P183P (CCG→CCAyejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (3/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCG  >  minE/1408211‑1408344
                                                                     |                                                                
cgAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGa                                                                 >  1:335207/1‑71 (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:3474781/1‑71 (MQ=255)
               gatgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACg                                                  >  1:397596/1‑71 (MQ=255)
                atgatTGCGATGGCGCTGTCATGCCGCCCGGCGGTGCTGATTGCCGATGAGCCAACCACCGCGCTGGatg                                                  <  1:2923797/70‑3 (MQ=255)
                 tgatTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtc                                                >  1:479258/1‑71 (MQ=255)
                    ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGa                                                    >  1:2963857/1‑64 (MQ=255)
                    ttGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGtc                                                >  1:3608624/1‑68 (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)
                            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)
                                                                     |                                                                
CGAACGGCAGCGGGTGATGATTGCGATGGCGCTGTTAACGCGACCGGAATTATTAATTGCCGATGAACCGACCACCGCACTGGACGTCTCTGTCCAGGCGCAGATTTTACAGCTGTTGCGCGAACTGCAAGGCG  >  minE/1408211‑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: