Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,905:1 +T 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,9051.T25.9% 159.2 / 67.6 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 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

TAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAC  >  minE/1467835‑1467972
                                                                        |                                                                          
tAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAt                                                                             <  1:3314348/71‑1 (MQ=255)
    tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTc                                                                   <  1:1857861/71‑1 (MQ=255)
          tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCg                                                             >  1:1061835/1‑69 (MQ=255)
           gATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          >  1:1395669/1‑71 (MQ=255)
            aTGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          <  1:3460347/70‑1 (MQ=255)
                ttGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCg                                                      >  1:21774/1‑70 (MQ=255)
                   gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa                                                  <  1:959386/71‑1 (MQ=255)
                   gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa                                                  <  1:2080018/71‑1 (MQ=255)
                       cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGc                                              <  1:2404202/71‑1 (MQ=255)
                          ggttGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGc                                              >  1:3763323/4‑71 (MQ=2)
                          cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCTGATAAGGCGtt                                           <  1:105602/71‑1 (MQ=255)
                           aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGttt                                          <  1:3945874/71‑1 (MQ=255)
                           aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGtt                                           >  1:785002/1‑70 (MQ=255)
                            aTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa                                         >  1:552824/1‑71 (MQ=255)
                             tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGCGt                                            >  1:155206/1‑70 (MQ=11)
                             tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa                                         <  1:3159968/70‑1 (MQ=255)
                              gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       >  1:1053630/1‑71 (MQ=255)
                               ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:58221/70‑1 (MQ=255)
                                  gATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt                                                   >  1:1678345/1‑55 (MQ=255)
                                  gATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt                                                   >  1:2093085/1‑55 (MQ=255)
                                  gATGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAc                                        >  1:3139993/1‑71 (MQ=0)
                                   aTGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:3890194/1‑71 (MQ=255)
                                   aTGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:2865534/71‑1 (MQ=0)
                                   aTGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:3152114/71‑1 (MQ=0)
                                      cGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                <  1:3065057/70‑1 (MQ=255)
                                      cGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg                                    <  1:1147494/71‑1 (MQ=0)
                                       gACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg                                    >  1:110259/1‑70 (MQ=0)
                                        aCGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  <  1:651183/71‑1 (MQ=0)
                                        aCGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:1521629/1‑71 (MQ=0)
                                         cGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                  >  1:725712/1‑70 (MQ=0)
                                          gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:2562409/1‑71 (MQ=0)
                                          gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:2531117/1‑71 (MQ=0)
                                          gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:3740559/1‑71 (MQ=0)
                                           cTTGCCGCGTCTTATCTGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCg                                                             <  1:3155041/36‑1 (MQ=255)
                                              gCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg                                                          >  1:3487816/1‑36 (MQ=255)
                                               acgcgTCTTATCAGGCCTAC‑AAATTCCCG‑‑CGCCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                            <  1:2864009/70‑1 (MQ=0)
                                                cgcgTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                            >  1:1977681/1‑70 (MQ=2)
                                                cgcgTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                           <  1:2418182/71‑1 (MQ=2)
                                                 gcgTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                    <  1:400610/71‑1 (MQ=255)
                                                       tATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                >  1:3449640/1‑69 (MQ=255)
                                                        aTCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTa              <  1:486145/70‑1 (MQ=255)
                                                          cAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGcc                   >  1:916471/1‑63 (MQ=255)
                                                            ggCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAg         <  1:3686530/71‑1 (MQ=255)
                                                             gCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt        >  1:3476971/1‑71 (MQ=255)
                                                             gCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGt        >  1:3221245/1‑71 (MQ=255)
                                                              ccTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTg       <  1:3586403/71‑1 (MQ=255)
                                                               cTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2792831/71‑1 (MQ=255)
                                                                  c‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg                                       >  1:495872/1‑35 (MQ=255)
                                                                  c‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                             >  1:610608/1‑45 (MQ=255)
                                                                  c‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2453992/68‑1 (MQ=255)
                                                                  c‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGa      <  1:2996247/68‑1 (MQ=255)
                                                                  c‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt     <  1:2304244/69‑1 (MQ=255)
                                                                  c‑AAAT‑‑‑‑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)
                                                                        |                                                                          
TAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑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: ‑