Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,908 T→A 24.9% intergenic (+67/+61) ubiG → / ← yfaL bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,467,9080TA24.9% 113.1 / 28.6 54intergenic (+67/+61)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base A (8/5);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAC  >  minE/1467839‑1467972
                                                                            |                                                                  
tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTc                                                                   <  1:1857861/71‑1 (MQ=255)
      tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCg                                                             >  1:1061835/1‑69 (MQ=255)
       gATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          >  1:1395669/1‑71 (MQ=255)
        aTGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          <  1:3460347/70‑1 (MQ=255)
            ttGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCg                                                      >  1:21774/1‑70 (MQ=255)
               gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa                                                  <  1:959386/71‑1 (MQ=255)
               gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa                                                  <  1:2080018/71‑1 (MQ=255)
                   cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGc                                              <  1:2404202/71‑1 (MQ=255)
                      ggttGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGc                                              >  1:3763323/4‑71 (MQ=2)
                      cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCTGATAAGGCGtt                                           <  1:105602/71‑1 (MQ=255)
                       aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGttt                                          <  1:3945874/71‑1 (MQ=255)
                       aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGtt                                           >  1:785002/1‑70 (MQ=255)
                        aTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa                                         >  1:552824/1‑71 (MQ=255)
                         tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGCGt                                            >  1:155206/1‑70 (MQ=11)
                         tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa                                         <  1:3159968/70‑1 (MQ=255)
                          gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       >  1:1053630/1‑71 (MQ=255)
                           ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:58221/70‑1 (MQ=255)
                              gATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt                                                   >  1:1678345/1‑55 (MQ=255)
                              gATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt                                                   >  1:2093085/1‑55 (MQ=255)
                              gATGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAc                                        >  1:3139993/1‑71 (MQ=0)
                               aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:3890194/1‑71 (MQ=255)
                               aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:3152114/71‑1 (MQ=0)
                               aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:2865534/71‑1 (MQ=0)
                                  cGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                <  1:3065057/70‑1 (MQ=255)
                                  cGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg                                    <  1:1147494/71‑1 (MQ=0)
                                   gACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg                                    >  1:110259/1‑70 (MQ=0)
                                    aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  <  1:651183/71‑1 (MQ=0)
                                    aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:1521629/1‑71 (MQ=0)
                                     cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:725712/1‑70 (MQ=0)
                                      gCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:2562409/1‑71 (MQ=0)
                                      gCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:2531117/1‑71 (MQ=0)
                                      gCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:3740559/1‑71 (MQ=0)
                                       cTTGCCGCGTCTTATCTGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCg                                                             <  1:3155041/36‑1 (MQ=255)
                                          gCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          >  1:3487816/1‑36 (MQ=255)
                                           acgcgTCTTATCAGGCCTACAAATTCCCG‑‑CGCCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                            <  1:2864009/70‑1 (MQ=0)
                                            cgcgTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                            >  1:1977681/1‑70 (MQ=2)
                                            cgcgTCTTATCAGGCCTACAAATTCCCG‑‑CACCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                           <  1:2418182/71‑1 (MQ=2)
                                             gcgTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                    <  1:400610/71‑1 (MQ=255)
                                                   tATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                >  1:3449640/1‑69 (MQ=255)
                                                    aTCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTa              <  1:486145/70‑1 (MQ=255)
                                                      cAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGcc                   >  1:916471/1‑63 (MQ=255)
                                                        ggCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg         <  1:3686530/71‑1 (MQ=255)
                                                         gCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt        >  1:3476971/1‑71 (MQ=255)
                                                         gCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt        >  1:3221245/1‑71 (MQ=255)
                                                          ccTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTg       <  1:3586403/71‑1 (MQ=255)
                                                           cTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2792831/71‑1 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       >  1:495872/1‑35 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                             >  1:610608/1‑45 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2453992/68‑1 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2996247/68‑1 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt     <  1:2304244/69‑1 (MQ=255)
                                                              cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt     <  1:1145492/69‑1 (MQ=255)
                                                                aaaT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc  >  1:437211/1‑71 (MQ=255)
                                                                aaaT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc  <  1:3428344/71‑1 (MQ=255)
                                                                            tCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt     <  1:1516659/62‑1 (MQ=255)
                                                                            |                                                                  
TGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAC  >  minE/1467839‑1467972

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