Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,467,910:2 +T 27.6% 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.T27.6% 176.0 / 76.4 58intergenic (+69/+59)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base . (26/16);  new base T (8/8);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCAC  >  minE/1467843‑1467980
                                                                         |                                                                      
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:5374231/71‑1 (MQ=0)
                           aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                               >  1:4963371/1‑71 (MQ=0)
                           aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                               >  1:4890722/1‑71 (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:4485858/1‑71 (MQ=0)
                             gCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                             >  1:29656/1‑71 (MQ=0)
                              cGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                             <  1:4567717/70‑1 (MQ=0)
                                aCGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                     >  1:1858644/1‑71 (MQ=255)
                                aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                          >  1:3608400/1‑71 (MQ=0)
                                aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                          >  1:1204526/1‑71 (MQ=0)
                                aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                          <  1:3533806/71‑1 (MQ=0)
                                 cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                         <  1:4720923/71‑1 (MQ=0)
                                 cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                         <  1:5289091/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:4990484/1‑71 (MQ=255)
                                         gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                            >  1:5225554/1‑71 (MQ=255)
                                         gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                            >  1:4659247/1‑71 (MQ=255)
                                         gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                            >  1:338519/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:2112006/1‑70 (MQ=255)
                                              ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                        >  1:5660861/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:1430356/1‑71 (MQ=255)
                                                          cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa           >  1:5600342/1‑71 (MQ=255)
                                                          cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa           >  1:1377418/1‑71 (MQ=255)
                                                           aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                          <  1:3969461/39‑1 (MQ=255)
                                                           aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                          <  1:108964/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)
                                                                    cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt     <  1:2905643/71‑1 (MQ=255)
                                                                    cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt     <  1:3436856/71‑1 (MQ=255)
                                                                    cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt     <  1:3545106/71‑1 (MQ=255)
                                                                    cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt     <  1:3545542/71‑1 (MQ=255)
                                                                       c‑‑ccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCAc  >  1:5001259/1‑71 (MQ=255)
                                                                         |                                                                      
GTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTTCAC  >  minE/1467843‑1467980

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: