Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,332,245 C→G 11.1% intergenic (+28/‑21) dksA → / → gluQ RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,332,2450CG11.1% 119.7 / 7.8 45intergenic (+28/‑21)dksA/gluQRNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base G (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCG  >  NC_002947/5332106‑5332386
                                                                                                                                           |                                                                                                                                              
gTGCGGCATCGAGATCGTCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATg                                                                                                                                             >  4:136805/1‑143 (MQ=255)
gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATg                                                                                                                                             <  3:136805/143‑1 (MQ=255)
gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCAGAACGGGGCGCATCATGCGcccc                                                                                                                                       <  6:39477/149‑1 (MQ=255)
gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCAGAACGGGGCGCATCATGCGcccc                                                                                                                                       <  2:171453/149‑1 (MQ=255)
    ggCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCTCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGCCCCGttc                                                                                                                                   <  3:105430/149‑2 (MQ=255)
     gCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGttt                                                                                                                                   >  1:71602/1‑148 (MQ=255)
       aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                <  8:19106/149‑1 (MQ=255)
       aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                                <  8:174864/149‑1 (MQ=255)
                             tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                           >  4:148134/1‑148 (MQ=255)
                             tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                           >  4:144673/1‑148 (MQ=255)
                                    ccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                                                   >  3:101198/1‑149 (MQ=255)
                                    ccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                                                   >  6:261393/1‑149 (MQ=255)
                                            ccgccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCt                                                                                            <  7:44087/147‑1 (MQ=255)
                                             cgccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTAc                                                                                          >  6:215866/1‑148 (MQ=255)
                                               ccgACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATc                                                                                       <  8:109479/149‑1 (MQ=255)
                                                   ccTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGgc                                                                                   <  4:101198/149‑1 (MQ=255)
                                                      gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑c                                                                                                 >  5:111984/1‑132 (MQ=255)
                                                      gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑c                                                                                                 <  6:111984/132‑1 (MQ=255)
                                                       tgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACTGTCGGCAAAGGCTGATCTTTCTTCCAGCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGgcgc                                                                                 <  1:266294/147‑1 (MQ=255)
                                                                  cAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGcccccagc                                                                     <  3:144673/148‑3 (MQ=255)
                                                                  cAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGcccccacc                                                                     <  3:148134/148‑1 (MQ=255)
                                                                   aaGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGcccccacccc                                                                   >  2:74211/1‑149 (MQ=255)
                                                                       aaaTTGCCGAGAAAAAGGAAAACACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGc                                                               <  8:116012/149‑1 (MQ=255)
                                                                            gCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCt                                                           <  5:261393/148‑1 (MQ=255)
                                                                               gagaAAAAGGATAAAACGGTCGGCAAAGGCTGATCTTTCATCCACCTGAACGCGGCGCATCATGCGCGCCGGTTCATTTATATGCCCAGCCTGCCCATGACCGACT‑CCAGCTACATCGGGCGCTTCTCCCCCACCCCCAGCGGGATCCt                                                       >  4:200436/1‑149 (MQ=255)
                                                                               gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAg                                                                                              <  3:212906/110‑1 (MQ=255)
                                                                               gagaAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                       <  5:215866/149‑1 (MQ=255)
                                                                               gagaAAAAGGAAAAAACCCTCGGCAAACGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGAGCACCGTTTCATTTATATGACCAGCCTGACCATGACCGAAT‑CCAg                                                                                              >  4:212906/1‑110 (MQ=255)
                                                                                       ggAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCg                                               <  2:4008/149‑1 (MQ=255)
                                                                                             aaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                              >  6:107359/1‑81 (MQ=255)
                                                                                             aaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAc                                                                                                              <  5:107359/81‑1 (MQ=255)
                                                                                             aaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGc                                         <  2:115631/149‑1 (MQ=255)
                                                                                              aCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGc                                         <  1:74211/149‑1 (MQ=255)
                                                                                                ggTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTgg                                      >  2:176808/1‑149 (MQ=255)
                                                                                                     gCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgcc                                 <  2:71602/149‑1 (MQ=255)
                                                                                                                     cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCt                            >  5:215783/1‑138 (MQ=255)
                                                                                                                     cttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCt                            <  6:215783/138‑1 (MQ=255)
                                                                                                                        ccACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                <  1:176808/147‑1 (MQ=255)
                                                                                                                          aatgaaaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                      >  4:84294/6‑139 (MQ=255)
                                                                                                                          aatgaaaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                      <  3:84294/134‑1 (MQ=255)
                                                                                                                           ccTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATg           >  5:214603/1‑149 (MQ=255)
                                                                                                                            tgaaaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGc          >  4:128524/4‑149 (MQ=255)
                                                                                                                               aaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcg       >  6:41300/1‑149 (MQ=255)
                                                                                                                               aaCGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcg       >  7:137938/1‑149 (MQ=255)
                                                                                                                                    ggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCg  >  4:152193/1‑149 (MQ=255)
                                                                                                                                        catcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc               <  1:161787/132‑1 (MQ=255)
                                                                                                                                        catcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc               >  2:161787/1‑132 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACT‑CCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCG  >  NC_002947/5332106‑5332386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑