Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,325 G→T 12.8% intergenic (+49/+105) ubiF → / ← glnX 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913696,3250GT12.8% 156.0 / 12.7 47intergenic (+49/+105)ubiF/glnX2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base T (3/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAG  >  NC_000913/696201‑696414
                                                                                                                            |                                                                                         
tAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                               >  3:116869/1‑133 (MQ=255)
tAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                               <  4:116869/137‑5 (MQ=255)
          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTctta                                                                                <  2:679809/126‑4 (MQ=255)
          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTctta                                                                                >  1:679809/1‑123 (MQ=255)
          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                    >  4:1068609/1‑138 (MQ=255)
          ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                    <  3:1068609/138‑1 (MQ=255)
                 cggcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                    >  4:829015/1‑115 (MQ=255)
                 cggcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                    <  3:829015/115‑1 (MQ=255)
                    cggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                        >  3:83738/1‑108 (MQ=255)
                    cggAGCGTACTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                        <  4:83738/108‑1 (MQ=255)
                                      tGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCgg                                                                                          <  1:409652/88‑1 (MQ=255)
                                      tGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCgg                                                                                          >  2:409652/1‑88 (MQ=255)
                                          aCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                           >  4:211267/1‑131 (MQ=255)
                                          aCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                           <  3:211267/131‑1 (MQ=255)
                                             tCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                    >  4:1091282/1‑87 (MQ=255)
                                             tCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                    <  3:1091282/87‑1 (MQ=255)
                                               aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     <  1:316162/148‑1 (MQ=255)
                                                      gAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     >  4:377896/1‑141 (MQ=255)
                                                      gAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     <  3:377896/141‑1 (MQ=255)
                                                      gAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCaaaa              >  2:682849/1‑148 (MQ=255)
                                                         atatGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCACCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaatga                              >  4:935467/1‑129 (MQ=255)
                                                         atatGCGTTAGGGTTATAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaatga                              <  3:935467/129‑1 (MQ=255)
                                                           atGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     >  2:997459/1‑136 (MQ=255)
                                                           atGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     <  1:997459/136‑1 (MQ=255)
                                                            tGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                           >  4:379574/1‑113 (MQ=255)
                                                            tGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAtt                                           <  3:379574/113‑1 (MQ=255)
                                                                  aGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAg  <  1:682849/148‑1 (MQ=255)
                                                                              ttACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                              >  1:1025675/1‑60 (MQ=255)
                                                                              ttACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                              <  2:1025675/60‑1 (MQ=255)
                                                                                 caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACg                                                          <  4:831891/77‑1 (MQ=255)
                                                                                 caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACg                                                          >  3:831891/1‑77 (MQ=255)
                                                                                   acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAt                                            >  1:296933/1‑89 (MQ=255)
                                                                                   acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAt                                            <  2:296933/89‑1 (MQ=255)
                                                                                    caTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGa                                        <  2:342642/92‑1 (MQ=255)
                                                                                    caTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGa                                        >  1:342642/1‑92 (MQ=255)
                                                                                      ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     >  2:949284/1‑109 (MQ=255)
                                                                                      ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaa                     <  1:949284/109‑1 (MQ=255)
                                                                                                  tGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                 >  4:1063328/1‑85 (MQ=255)
                                                                                                  tGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                 <  3:1063328/85‑1 (MQ=255)
                                                                                                               tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGa                                                                     <  4:415335/36‑1 (MQ=255)
                                                                                                               tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGa                                                                     >  3:415335/1‑36 (MQ=255)
                                                                                                               tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAt                                     >  2:865586/1‑68 (MQ=255)
                                                                                                               tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAt                                     <  1:865586/68‑1 (MQ=255)
                                                                                                                 ggATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                  <  1:469427/69‑1 (MQ=255)
                                                                                                                 ggATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                  >  2:469427/1‑69 (MQ=255)
                                                                                                                   aTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACAt                                                 >  1:969356/1‑52 (MQ=255)
                                                                                                                   aTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACAt                                                 <  2:969356/52‑1 (MQ=255)
                                                                                                                    tgcgACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                  >  1:829794/4‑66 (MQ=25)
                                                                                                                    tgcgACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtga                                  <  2:829794/63‑1 (MQ=25)
                                                                                                                            |                                                                                         
TAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAG  >  NC_000913/696201‑696414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: