Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,905:4 +C 25.9% 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,9054.C25.9% 158.4 / 69.5 54intergenic (+64/+64)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base . (19/21);  new base C (8/6);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAC  >  minE/1467835‑1467972
                                                                           |                                                                       
tAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAt                                                                             <  1:3314348/71‑1 (MQ=255)
    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:2865534/71‑1 (MQ=0)
                                   aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:3152114/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:3428344/71‑1 (MQ=255)
                                                                    aaaT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc  >  1:437211/1‑71 (MQ=255)
                                                                           |                                                                       
TAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAC  >  minE/1467835‑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: