Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,910:2 +T 25.9% intergenic (+69/+59) 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,9102.T25.9% 158.6 / 67.4 54intergenic (+69/+59)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 T (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 = 1.00e+00

TGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTT  >  minE/1467845‑1467977
                                                                          |                                                                   
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:3221245/1‑71 (MQ=255)
                                                   gCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt             >  1:3476971/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:2996247/68‑1 (MQ=255)
                                                        cAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa           <  1:2453992/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)
                                                                      tCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt          <  1:1516659/62‑1 (MQ=255)
                                                                       cc‑‑ccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt  <  1:2240166/69‑1 (MQ=255)
                                                                          |                                                                   
TGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTT  >  minE/1467845‑1467977

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: