Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,905:2 +C 29.1% intergenic (+64/+64) 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,9052.C29.1% 162.6 / 81.0 55intergenic (+64/+64)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base . (26/13);  new base C (8/8);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAT  >  minE/1467838‑1467975
                                                                     |                                                                          
gTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCt                                                                       <  1:1940781/71‑1 (MQ=255)
     gTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTcc‑‑cc                                                                 <  1:843942/70‑1 (MQ=255)
       tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt                                                               >  1:2447871/1‑70 (MQ=255)
          tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc                                                           >  1:3144410/1‑71 (MQ=255)
          tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc                                                           <  1:2023096/71‑1 (MQ=255)
              tGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGa                                                       >  1:3473915/1‑71 (MQ=255)
                  cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg                                                   <  1:3733454/71‑1 (MQ=255)
                  cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg                                                   <  1:468940/71‑1 (MQ=255)
                    cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGc                                                 <  1:150269/71‑1 (MQ=255)
                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt                                               >  1:1534616/1‑71 (MQ=255)
                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt                                               >  1:5322616/1‑71 (MQ=255)
                       cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                              <  1:1838129/71‑1 (MQ=255)
                       cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                              <  1:647468/71‑1 (MQ=255)
                        aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                              >  1:125423/1‑70 (MQ=255)
                            ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt                                                               >  1:5569488/1‑49 (MQ=255)
                              tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                <  1:2097987/46‑1 (MQ=255)
                                aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     >  1:2453231/1‑71 (MQ=255)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                          >  1:4890722/1‑71 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                          >  1:4963371/1‑71 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                          <  1:5374231/71‑1 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAc                                           <  1:4719260/70‑1 (MQ=0)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                >  1:1700515/1‑43 (MQ=255)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                >  1:37603/1‑43 (MQ=255)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     <  1:4427823/70‑1 (MQ=255)
                                 tGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcgc                                         >  1:1345902/1‑71 (MQ=0)
                                  gCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                        >  1:29656/1‑71 (MQ=0)
                                  gCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                        >  1:4485858/1‑71 (MQ=0)
                                   cGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                        <  1:4567717/70‑1 (MQ=0)
                                     aCGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                >  1:1858644/1‑71 (MQ=255)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     >  1:1204526/1‑71 (MQ=0)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     >  1:3608400/1‑71 (MQ=0)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     <  1:3533806/71‑1 (MQ=0)
                                      cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                    <  1:5289091/71‑1 (MQ=0)
                                      cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                    <  1:4720923/71‑1 (MQ=0)
                                       gCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                               >  1:3375370/1‑70 (MQ=255)
                                       gCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                   <  1:445947/71‑1 (MQ=0)
                                           tgcgcgTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                >  1:3564206/3‑71 (MQ=2)
                                            acgcgTCTTATCAGGCCTACAAATTCCCGCGCCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                 <  1:2023343/68‑1 (MQ=0)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                       >  1:4659247/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                       >  1:338519/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                       >  1:5225554/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                       >  1:4990484/1‑71 (MQ=255)
                                               cgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGcc                      <  1:3636583/71‑1 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGa                    >  1:4731841/1‑69 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                   >  1:5660861/1‑70 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                   >  1:2112006/1‑70 (MQ=255)
                                                             tACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt        >  1:545890/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAc                            >  1:2670849/1‑49 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa      >  1:1377418/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa      >  1:1430356/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa      >  1:5600342/1‑71 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     <  1:108964/39‑1 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                     <  1:3969461/39‑1 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc     >  1:695619/1‑71 (MQ=255)
                                                                  aT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAACGTGCCGACTAGACAGTGATTACCa   >  1:5178548/1‑71 (MQ=255)
                                                                   t‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt  <  1:2224276/71‑1 (MQ=255)
                                                                     |                                                                          
GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAT  >  minE/1467838‑1467975

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: