Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,062 G→C 9.2% intergenic (+39/+64) 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,0620GC 9.2% 134.8 / 2.5 50intergenic (+39/+64)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base C (3/2);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCC  >  NC_002947/5536915‑5537207
                                                                                                                                                   |                                                                                                                                                 
gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                  >  7:251841/1‑149 (MQ=255)
  tCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                <  4:121827/149‑1 (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                                                                                                                                           <  1:232707/122‑1 (MQ=255)
                                  cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                           >  2:232707/1‑122 (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                                                           <  7:167795/108‑1 (MQ=255)
                                                                                                                                ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                           >  8:167795/1‑108 (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)
                                                                                                                                                cGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc  <  2:342184/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCC  >  NC_002947/5536915‑5537207

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: ‑