Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,960,521	A→T	7.7%	intergenic (‑137/‑192)	mnxG ← / → PP_3491	manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,960,521	0	A	T	7.7%	139.1 / 3.3	52	intergenic (‑137/‑192)	mnxG/PP_3491	manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base T (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTC  >  NC_002947/3960375‑3960669                                                                                                                                                   |                                                                                                                                                     gcgTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATc                                                                                                                                                    >  4:88937/1‑149 (MQ=255)  cgTCGACATGCTGCACCTCCGTGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGCATACCCCAGCAGGATTCCCCGATTCg                                                                                                                                                   >  8:159808/1‑149 (MQ=255)   gTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCgg                                                                                                                                                  <  3:16003/149‑1 (MQ=255)        cATGCTGCAACTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTTCCCCGATTCGGGGGaa                                                                                                                                             <  7:159808/149‑2 (MQ=255)         aTGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGATTCGGGGGaac                                                                                                                                            <  3:341515/149‑3 (MQ=255)                     cGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCTCGAATCGGGGGACGGAATCCCGTg                                                                                                                                  <  8:52617/147‑1 (MQ=255)                          cacaGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCa                                                                                                                            >  8:82395/1‑148 (MQ=255)                           acaGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCa                                                                                                                            >  2:186371/1‑147 (MQ=255)                             aGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCGAGAAGATAAGTTCTTTATTTATCTTTTACTTGTAGGATGAACCGCAACGGATTGGGGAGGATCCCCCAGCAGCTTCCCGCGAATCGGGGGACGGAATCCCGTGGGGGCATAgg                                                                                                                        >  8:197807/1‑149 (MQ=255)                                 tctcTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGgctgc                                                                                                                    >  2:138738/1‑149 (MQ=255)                                        ggggAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACaa                                                                                                             <  4:135723/149‑1 (MQ=255)                                          ggAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAgc                                                                                                           >  2:319015/1‑149 (MQ=255)                                          ggAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAgc                                                                                                           <  7:82395/149‑1 (MQ=255)                                             gCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAgcgcg                                                                                                        >  2:126460/1‑149 (MQ=255)                                                 gAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCt                                                                                                    >  1:163910/1‑149 (MQ=255)                                                                cAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCggg                                                                                                                                                 <  1:301865/89‑1 (MQ=255)                                                                cAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCggg                                                                                                                                                 >  2:301865/1‑89 (MQ=255)                                                                   aagaaAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCACGAATCGGGGGACGGAATACCGTGGGGGCAAAGGCTGCTCCACAACCGCGCGCTGGCCTCACACAAGGTaa                                                                                   >  5:199041/1‑148 (MQ=255)                                                                    agaaAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTaaa                                                                                  >  2:314905/1‑148 (MQ=255)                                                                           ttCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACg                                                                          >  6:196853/1‑149 (MQ=255)                                                                             cAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGTGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCTTGGGGGCAAAGGATGCTCCACAAGCGCGCTCTGGCGTCACACAAGGTAATAACGCACgg                                                                         >  1:251057/1‑148 (MQ=255)                                                                                tatTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTAcc                                                                     >  3:199580/1‑149 (MQ=255)                                                                                    aaTCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCtt                                                                  <  1:319015/148‑1 (MQ=255)                                                                                    aaTCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTg                                                                 >  6:320998/1‑149 (MQ=255)                                                                                       cAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGca                                                               >  2:175217/1‑148 (MQ=255)                                                                                         aTAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACaa                                                            >  1:283540/1‑149 (MQ=255)                                                                                            aCTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAAtt                                                          <  1:175217/148‑1 (MQ=255)                                                                                                gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGa                                                     >  3:341971/1‑149 (MQ=255)                                                                                                 tACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                      <  8:58915/147‑1 (MQ=255)                                                                                                 tACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                      >  7:58915/1‑147 (MQ=255)                                                                                                   cGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtg                                                  <  1:314905/149‑1 (MQ=255)                                                                                                    gagacaCCGCAACGGATTGGGGTGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtgt                                                 <  6:199041/144‑1 (MQ=255)                                                                                                       cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTc                                                       <  4:110568/140‑1 (MQ=255)                                                                                                       cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTc                                                       >  3:110568/1‑140 (MQ=255)                                                                                                             cAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCaa                                        <  3:88937/149‑1 (MQ=255)                                                                                                                   aTTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGc                                  <  5:320998/149‑1 (MQ=255)                                                                                                                    ttGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGTGCAAAGGCTTCACCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCa                                 <  7:197807/149‑1 (MQ=255)                                                                                                                         gaggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCa                            <  5:196853/149‑1 (MQ=255)                                                                                                                                ccccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCt                     >  2:264016/1‑149 (MQ=255)                                                                                                                                 cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTg                    <  1:264016/149‑1 (MQ=255)                                                                                                                                 cccAGCAGGTTCCCCCGAACCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTg                    <  2:283780/149‑1 (MQ=255)                                                                                                                                  ccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc                   >  5:228551/1‑149 (MQ=255)                                                                                                                                        ccgtCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTgg                                                                                                  >  8:135177/4‑64 (MQ=255)                                                                                                                                         gTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACACGCGCGCGCTGGCATCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcctc            >  1:10827/1‑149 (MQ=255)                                                                                                                                         gTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcctc            >  5:151724/1‑149 (MQ=255)                                                                                                                                         gTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcctc            >  2:152542/1‑149 (MQ=255)                                                                                                                                             ccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGc        <  3:227785/149‑1 (MQ=255)                                                                                                                                              cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcct             <  3:221941/143‑1 (MQ=255)                                                                                                                                              cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcct             >  4:221941/1‑143 (MQ=255)                                                                                                                                               ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                            <  6:296986/111‑1 (MQ=255)                                                                                                                                               ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                            >  5:296986/1‑111 (MQ=255)                                                                                                                                                cGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa     <  1:152542/149‑1 (MQ=255)                                                                                                                                                   aTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTc  >  2:185110/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGTCGACATGCTGCACCTCCGGGACACAGGATCTCTGCGGGGAGCTGGAGCAAGGTGTGTGCCAGAAGAAAAGTTCAATATTAATCAATAACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTC  >  NC_002947/3960375‑3960669

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