Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,960,617 G→A 12.7% intergenic (‑233/‑96) 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_0029473,960,6170GA12.7% 114.7 / 8.4 47intergenic (‑233/‑96)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (3/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGA  >  NC_002947/3960470‑3960762
                                                                                                                                                   |                                                                                                                                                 
gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                   <  3:154957/148‑1 (MQ=255)
gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTCCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGa                                                                                                                                                  >  7:64611/1‑149 (MQ=255)
       cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTTCTCCACCAGCCCGCGCTGGCCTCAAACAAGGTAAAAACGCCCGGTACCCATTCACAATTGTCGATTGTGc                                                                                                                                            >  1:247998/1‑148 (MQ=255)
           cGCAACGGATTGGGGAGGATCCCGCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGc                                                                                                                                       <  5:247729/149‑1 (MQ=255)
            gCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGc                                                                                                                                       >  3:262929/1‑148 (MQ=255)
                   aTTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCg                                                                                                                                >  2:279428/1‑148 (MQ=255)
                      gcggaggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGGGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCAt                                                                                                                             <  2:243437/146‑1 (MQ=255)
                         gaggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGc                                                                                                                           >  7:26015/1‑147 (MQ=255)
                            gaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGACGATTGTGCAAGGCAAATTCGCATGCCAGGa                                                                                                                      >  7:212901/1‑149 (MQ=255)
                                     gcagGTTCCCCCGATTGTGGGGACGTAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAt                                                                                                             <  2:269468/149‑1 (MQ=255)
                                        ccgtCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcct                                                                                                          >  8:248972/4‑149 (MQ=255)
                                              cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCg                                                                                                    <  1:279428/149‑1 (MQ=255)
                                                  aaTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACt                                                                                                <  4:262929/149‑1 (MQ=255)
                                                       ggggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                           <  5:104609/149‑1 (MQ=255)
                                                       ggggACGGAATCCAGTGGGGGCAAAGGCTGCTCCTCAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                           <  8:64611/149‑1 (MQ=255)
                                                         ggACGGAATCCCGTGGGGGCAAAGGCTACTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACt                                                                                         >  1:182790/1‑149 (MQ=255)
                                                             ggAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCAc                                                                                     >  6:111615/1‑149 (MQ=255)
                                                                  cccGTGGGGGCAAAGGCTGCTCGACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCGGCGATCCTCCAGCGAACTCCAGACTTAACTTgg                                                                                 >  7:205648/1‑148 (MQ=255)
                                                                   ccGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGa                                                                                <  7:248972/148‑1 (MQ=255)
                                                                          ggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                        <  2:182790/149‑1 (MQ=255)
                                                                             aaaGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGCCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAAc                                                                     <  8:212901/149‑1 (MQ=255)
                                                                                  ctgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAg                                                                 >  3:319082/1‑148 (MQ=255)
                                                                                     ctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAgg                                                             <  8:26015/149‑1 (MQ=255)
                                                                                      tCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAgg                                                             >  3:181228/1‑148 (MQ=255)
                                                                                         acaAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGgag                                                          >  3:29027/1‑148 (MQ=255)
                                                                                          caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGaa                                                        <  8:51716/149‑1 (MQ=255)
                                                                                                       gCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGca                                           >  7:229541/1‑149 (MQ=255)
                                                                                                        ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGGTTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGc                                                     <  7:219767/138‑1 (MQ=255)
                                                                                                        ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGc                                                     >  8:219767/1‑138 (MQ=255)
                                                                                                        ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGca                                           >  3:285658/1‑148 (MQ=255)
                                                                                                         cTCACACAAGGTAAATACGCACGGTAACCTTTCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCATCCAGAGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGcaca                                         <  4:181228/149‑1 (MQ=255)
                                                                                                                  ggTATAATCGGACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAc                                <  6:184268/149‑1 (MQ=255)
                                                                                                                    tAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                              <  2:247998/149‑1 (MQ=255)
                                                                                                                     aaaaaCGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAcccc                             >  2:207804/1‑149 (MQ=255)
                                                                                                                         acggacggTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGgca                         <  4:167993/149‑1 (MQ=255)
                                                                                                                              cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                                       <  4:227615/66‑1 (MQ=255)
                                                                                                                              cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                                       >  3:227615/1‑66 (MQ=255)
                                                                                                                                    ccTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCa              >  6:176210/1‑149 (MQ=255)
                                                                                                                                     cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCgg                           <  8:184184/135‑1 (MQ=255)
                                                                                                                                     cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCAt             <  5:176210/149‑1 (MQ=255)
                                                                                                                                     cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCAgg                           >  7:184184/1‑135 (MQ=255)
                                                                                                                                      ttGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGc                                                                >  6:312208/1‑97 (MQ=255)
                                                                                                                                      ttGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGc                                                                <  5:312208/97‑1 (MQ=255)
                                                                                                                                        gCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg            <  8:205648/147‑1 (MQ=255)
                                                                                                                                          acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                        <  6:3271/85‑1 (MQ=255)
                                                                                                                                          acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                        >  5:3271/1‑85 (MQ=255)
                                                                                                                                                gTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGcc    >  2:3933/1‑147 (MQ=255)
                                                                                                                                                 tCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGa  <  1:3933/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGA  >  NC_002947/3960470‑3960762

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