Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,058 A→G 8.1% intergenic (+35/+68) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,537,0580AG 8.1% 135.4 / 2.9 52intergenic (+35/+68)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACC  >  NC_002947/5536915‑5537200
                                                                                                                                               |                                                                                                                                              
gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                           >  7:251841/1‑149 (MQ=255)
  tCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                         <  4:121827/149‑1 (MQ=255)
    gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATa                                                                                                                                                >  7:288032/1‑140 (MQ=255)
    gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATa                                                                                                                                                <  8:288032/140‑1 (MQ=255)
     aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAc                                                                                                                                               <  4:161495/140‑1 (MQ=255)
     aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAc                                                                                                                                               >  3:161495/1‑140 (MQ=255)
     aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                      <  6:25180/149‑1 (MQ=255)
        ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTAACCCGCGCAAACGAGGGTGAAacc                                                                                                                                   >  6:185957/1‑149 (MQ=255)
         cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcaca                                                                                                                                  >  2:42184/1‑149 (MQ=255)
            ccTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                <  1:154220/148‑1 (MQ=255)
                 aTGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCTCGTATTCGCGGGTGAAcccgc                                                                                                                                 <  1:201565/142‑3 (MQ=255)
                 aTGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCTCGTATTCGCGGGTGAAcccgc                                                                                                                                 >  2:201565/1‑140 (MQ=255)
                  tGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                         >  3:197814/1‑149 (MQ=255)
                       aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                    <  4:197814/149‑1 (MQ=255)
                        gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                   <  5:203775/149‑1 (MQ=25)
                        gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                   >  2:11411/1‑149 (MQ=255)
                        gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                   <  7:299405/149‑1 (MQ=255)
                            gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCCGGAGCGGGTTGTCCCGCGAAGACGAGGGGGGCACCACCGACGCATTCGCggggg                                                                                                                >  8:235645/1‑146 (MQ=255)
                             ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                              <  1:11411/149‑1 (MQ=255)
                                  cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                    >  2:232707/1‑122 (MQ=255)
                                  cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                    <  1:232707/122‑1 (MQ=255)
                                  cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACCCAGTAGCGGGTGAACCCGc                                                                                                         >  5:127977/1‑149 (MQ=255)
                                    ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCCGCACAAGCGCGACCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGGGAACCCACCGACGCATTCGCGGGTGAACCCGCTc                                                                                                       <  3:328988/149‑1 (MQ=255)
                                    ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTc                                                                                                       <  8:251841/149‑1 (MQ=255)
                                     gCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCCCGCATTCGCGGGTGGACCCGCtgc                                                                                                      >  7:78211/1‑147 (MQ=255)
                                      cAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCt                                                                                                     >  1:5723/1‑149 (MQ=255)
                                          gTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAg                                                                                                 >  1:32916/1‑149 (MQ=255)
                                             cGACATCATCGGCGCCACCAACCCCAAGCCCCAGCACAAGCCCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGCGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAggg                                                                                               <  7:235645/148‑1 (MQ=38)
                                                catcatGGTGGACAACTACCTCATCACCCAGAACTAGCTCCAACTGGCGTTGGATCTTTGATATGGGTGCTGCTTGTTGGGGTTCCACGGCTATTACGAGGTTTAACCCACCGTGTCATTCGCGGTTTAACCCGCTCCTGCAGGGTTCg                                                                                           <  8:73342/149‑1 (MQ=255)
                                                attcatcGTCGCCACCTAGCGCAAGACCCTGCACAAGGGCGAACTGCCTTGCGGTCCTTTATAAGTGGCCTGCAGGTGGGGGTTCACCGGCTAATAGGTGGGTGAACCCACCGACGCTTTGGCGGGTGAACCCGCTCCTGCAGGGTTc                                                                                            <  8:78211/146‑1 (MQ=255)
                                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                   >  1:121689/1‑122 (MQ=18)
                                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                   <  2:121689/122‑1 (MQ=18)
                                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAACCCCCCGACGCATTCGCGGGTGGACCCGCGCCTGGAGGGTTAGCCg                                                                                        >  3:94262/1‑149 (MQ=255)
                                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCGAATACGAGGGTTCAACCCCCCACGCCTTAGCGGGTGAACCCCCAACAGCACGGTTGGCCg                                                                                        >  8:199883/1‑149 (MQ=255)
                                                     tcgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTACGGCAGGGGTCGCCGt                                                                                       >  2:230132/1‑148 (MQ=255)
                                                     tcgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCCAAAACCACGGTGAACCCACCCACCCATTCGCGGGTGTACACCCGCCAGCACGGTTAGCCGt                                                                                       >  7:197465/1‑148 (MQ=255)
                                                        tcgCCACATACCGCAAGACCCAGCACACGCGCTAACTGCCGTTCGCTCTTTTATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGc                                                                                   <  4:94262/149‑1 (MQ=38)
                                                            caacTACCGCAAGACCCTCCAAATGGGCTAAGTGCCGTTCGCTCGTTTATATGTGGCCTGCTGTTGGGGGTTCACCCGCTAAAACTAGGTTTAACCCACCGACGCATTGGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCttt                                                                               <  6:127977/146‑1 (MQ=255)
                                                             aCCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTCCAGGCGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTg                                                                              <  5:185957/149‑1 (MQ=255)
                                                                     aaaGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTg                                                                      <  7:199883/148‑1 (MQ=255)
                                                                           ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                     <  5:261559/144‑1 (MQ=255)
                                                                                 aaggagcgAACGGCAGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAACCCACCCCCCCATTCGCGGGTGAACCCCCTCCTGGAGGGTTCGCCGTCAGCCGTTGATTACTTTACGAAGACCATg                                                          >  8:91713/6‑149 (MQ=255)
                                                                                    gcgcGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCt                                                       <  2:32916/149‑1 (MQ=255)
                                                                                      gcgAACGGCAGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGACCCCACCGACGCCTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGTTTACTTGACGATGACCTTGCCttt                                                     >  5:245634/1‑149 (MQ=38)
                                                                                                        cTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTgg                                     <  1:230132/147‑1 (MQ=255)
                                                                                                                 tatccTCCAGGAGCGGGTTCACCCGCGAATACGCGGTTTTACCCACGTTCCCATGCGCGGGCTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTCCCTTTCATCATCGAGATGTGGCCCGGGAAg                            <  1:29409/144‑1 (MQ=255)
                                                                                                                                ggTTCACCCGCTAATACGAGGGGGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa           <  2:5723/149‑1 (MQ=255)
                                                                                                                                ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                    >  8:167795/1‑108 (MQ=255)
                                                                                                                                ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                    <  7:167795/108‑1 (MQ=255)
                                                                                                                                  ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg          >  4:83971/1‑148 (MQ=255)
                                                                                                                                  ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt         <  3:83971/149‑1 (MQ=255)
                                                                                                                                    cACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca       >  8:231817/1‑149 (MQ=25)
                                                                                                                                      cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac      >  1:342184/1‑148 (MQ=255)
                                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCAGCTCCTGCAGGGTGCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc  >  8:132593/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACC  >  NC_002947/5536915‑5537200

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