Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,328 A→C 12.9% intergenic (+52/+102) ubiF → / ← glnX 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol hydroxylase/tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913696,3280AC12.9% 105.4 / 9.7 51intergenic (+52/+102)ubiF/glnX2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol hydroxylase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base C (5/2);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

CGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGA  >  NC_000913/696190‑696467
                                                                                                                                          |                                                                                                                                           
cGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGACTTACAACATTGCCGGATCCGTGCCCACCGTAGGCCGGATCAGACGCGGCAGCGt                                                                                                                                         >  2:501329/1‑143 (MQ=255)
    ttATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGGCGCa                                                                                                                                     >  1:85699/1‑143 (MQ=255)
      aTGCGTACCCTCGGGTCAGTGGCGGCGGTGTGTGCTGGCGTGTTGAACCGTCAGGCGCTGACGTATGCGTTAGCGTTGTCGCCTTACACCATTGCCGGATGCGTGCCACCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                                                                   <  2:345126/143‑1 (MQ=255)
               ctcGGGTTAATGGCGGCGGAGCGTGCTGGCGCGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTCTTATCCGAcctacg                                                                                                                          <  1:649979/143‑6 (MQ=37)
                        aTGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGa                                                                                                                 >  2:46094/1‑143 (MQ=255)
                           gcggcgGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACgt                                                                                                              >  1:605747/1‑143 (MQ=255)
                           gcggcgGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACgt                                                                                                              >  1:672362/1‑143 (MQ=255)
                            cgggggAGCGTGCTGGGGTGGTGCAACGTCAGGCGGTGAAATGTGGGTTAGGGTTGTGGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACgtg                                                                                                             <  2:585436/143‑1 (MQ=38)
                             ggcggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGAAgtgt                                                                                                            >  2:119874/1‑143 (MQ=255)
                               cggAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGGTGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCCTCAGACATTGAAGGATAAGACGGGTCa                                                                                                          >  2:833751/1‑143 (MQ=255)
                                 gAGCGGGGTGGGGTGGTGAAACGGCAGGCGCTGAAGTATGCGTGAGGGTTGTAGCCGTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAAc                                                                                                        <  2:200305/143‑1 (MQ=255)
                                      tGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACAGTGAAGGATAAGACGTGTCAACATCGc                                                                                                   >  1:94069/1‑143 (MQ=255)
                                             gtgtGGAAACGTCAGGCGCGGAAATATGGGTTAGGGGTGGAGCCTTACAACAGTGCCGGATTCGTGCCAACCGTAGGGCGGAGAAGACGCGTCAGCGTCGCAGCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGAc                                                                                            <  1:646061/143‑1 (MQ=255)
                                              tgtTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAAGATTCCTGTATACTTGGCAACCTTAGGTTTGAGAAGACGCGGCAGCGTCGCATGCGACATTGAAGGATAAGACGTGTCAACATCGCACTCGACa                                                                                           >  2:920010/1‑143 (MQ=255)
                                               gtTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCATTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGAAAt                                                                                          >  2:312677/1‑143 (MQ=255)
                                                ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAtt                                                                                         >  2:670342/1‑143 (MQ=255)
                                                ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAtt                                                                                         >  2:287244/1‑143 (MQ=255)
                                                ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCAGACATTGAAGGAAAAGACGCGTCAACATCGCATTCGACAtt                                                                                         >  1:888572/1‑143 (MQ=255)
                                                  gAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGGTAGGACGTGTCAACATCGCATTCGACATtga                                                                                       >  1:860003/1‑143 (MQ=255)
                                                   aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGCCGGATAAGACGCGGCAGCGTCGCCCCCGACATCGAAGGCTAAGACGCGTCAACCTCGCATTCGACCTcgaa                                                                                      >  2:14949/1‑143 (MQ=255)
                                                   aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCGTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACGTtgaa                                                                                      >  1:876504/1‑143 (MQ=255)
                                                                tGAAGTATGCGTTAGGGTTGTCGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                                                         <  1:345704/143‑1 (MQ=255)
                                                                tGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                                                         <  2:651647/143‑1 (MQ=255)
                                                                                  tGTAGTCTTACAACAGTGCCGGATGCGTGCCACCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCTGACATTGCAGGATAAGACGTGTGAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAg                                                       <  2:254180/143‑1 (MQ=39)
                                                                                     aGCCTCACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGACAAGACGCGGCAGCGCCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGCAAAGCAAAAAGCTCGCCGAAGCGa                                                    <  1:554573/143‑1 (MQ=255)
                                                                                      gcgtgcaaacgTTGCCGGATGGGGGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAAGGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAg                                                   <  2:729709/136‑1 (MQ=38)
                                                                                           acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAAACGCGGCAACGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAATCAATGAACACAGAAAAGCAAAAACATCGCCCAAACGAGCtttt                                              >  2:90816/1‑143 (MQ=255)
                                                                                             aacaTTGCCGGATGCGTGCCCACCGTAGGTCGTATAAGACTCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCtttttt                                            <  1:731648/143‑1 (MQ=255)
                                                                                             aacaTTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGGCAACATCGCATTCGACACTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCgttttt                                            >  1:827697/1‑143 (MQ=255)
                                                                                                  tGCCGGAGGCGGCCCAACCGTAGGGCGGATCAGACGCGGCAGCGTCGCATCCGACGTTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtgt                                       <  2:784604/143‑1 (MQ=255)
                                                                                                         tGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTggg                                <  1:820678/143‑1 (MQ=255)
                                                                                                            gTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTa                             >  2:531548/1‑143 (MQ=255)
                                                                                                               ccAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTGCGa                          >  1:884158/1‑143 (MQ=255)
                                                                                                                   ccGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGAt                      >  2:314507/1‑143 (MQ=255)
                                                                                                                    cGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGAtt                     >  1:814524/1‑143 (MQ=255)
                                                                                                                    cGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTCTTTTAATGTGGCTGGGGTACGAGGAtt                     >  2:436238/1‑143 (MQ=255)
                                                                                                                     gTAGGGCGGATGAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGGCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTc                    <  2:888572/143‑1 (MQ=255)
                                                                                                                       aGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGa                  <  1:363573/143‑1 (MQ=255)
                                                                                                                        ggTCGGATAAGACGCGGCAGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGaa                 <  1:944563/143‑1 (MQ=255)
                                                                                                                         gTCGTATAAGCCGCGGCAGCGTCGCATCCGACATCGCACGTTAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAAc                <  1:32191/143‑1 (MQ=255)
                                                                                                                           cGGATGCGACGCTGCCGCGTCTTATCCGACATGGAAGGATAAGACGTGTCAACATCGCATTCGACATTGCATGCACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGagact              >  2:490497/1‑139 (MQ=255)
                                                                                                                             ggtgaGACGCGGCAGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg            <  1:585402/139‑1 (MQ=255)
                                                                                                                             gatgcgACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg            >  2:889825/6‑143 (MQ=255)
                                                                                                                             gatgcgACGCTGCCGCGTCTTATCCGACATGGAAGGATAAGACGTGTAAACATAGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCCTTATTAATGTGGCTGGGGTACGATGCAGCGAACCTCg            >  2:71171/6‑143 (MQ=255)
                                                                                                                             gatgcgACGCTGCCGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg            >  1:537814/6‑143 (MQ=255)
                                                                                                                                  gACGCTGCCGCGTCTTATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGAAAGAGGATACGAACCTCGGaaaa       >  2:921375/1‑141 (MQ=255)
                                                                                                                                   aCGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc      >  1:321276/1‑143 (MQ=255)
                                                                                                                                   aCGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGACGATTCGAACCTCGGAATGc      >  1:945175/1‑143 (MQ=255)
                                                                                                                                    cgcgGCCGCGTCGCATCCGACGTTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGcc     <  2:672362/143‑1 (MQ=255)
                                                                                                                                    cgcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGcc     <  2:591306/143‑1 (MQ=255)
                                                                                                                                       ggCAGCCTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGa  <  2:250073/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
CGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGA  >  NC_000913/696190‑696467

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: